The following notes have been formated nicely in www.lstc.com/sdb/358 ===================================================================== For some slides on this topic, see http://ftp.lstc.com/anonymous/outgoing/support/FAQ_docs/hourglass.pdf . Guidelines provided herein for hourglass control are not to be considered as hard and fast rules, but instead may be considered as a starting point. The user is responsible for observing the deformed shapes of the elements (to see if hourglassing modes are evident) and for checking hourglass energy in a manner described herein and in http://ftp.lstc.com/anonymous/outgoing/support/FAQ/energy_balance . Trial and error may be necessary to determine the most effective hourglass control for a particular problem. Starting with 971 R3.1, the default hourglass control for underintegrated elements is: - Shells: viscous type (1=2=3) for explicit, stiffness type (4=5) for implicit - Solids: type 2 for explicit, type 6 for implicit (only 6 is implemented for implicit) See the User's Manual (*HOURGLASS) and sections 3.2 and 7.4 of the 2006 LS-DYNA Theory Manual. Hourglass (HG) modes are nonphysical, zero-energy modes of deformation that produce zero strain and no stress. Hourglass modes occur only in under-intetgrated (single integration point) solid, shell, and thick shell elements. LS-DYNA has various algorithms for inhibiting hourglass modes. The first algorithm (type 1), while the cheapest, is generally not the most effective algorithm. A way to entirely eliminate hourglass concerns is to switch to element formulations with fully-integrated or selectively reduced (S/R) integration. There can be a downside to this approach. For example, Type 2 solids are much more expensive than the single point default solid. Secondly, they are much more unstable in large deformation applications (negative volumes much more likely). Third, type 2 solids have some tendency to 'shear-lock' and thus behave too stiffly in applications where the element shape is poor. To overcome the shear locking problem, selectively-reduced ELFORMs -1 and -2 were developed and are generally recommended over ELFORM 2. Triangular shells and tetrahedral solid elements do not have hourglassing modes but have drawbacks with regard to overly stiff behavior in many applications. A good way to reduce hourglassing is to refine your mesh. The method of loading can affect the degree of hourglassing. A pressure loading is preferred over loading individual nodes as the latter approach is more likely to excite hourglassing modes. To evaluate hourglass energy, set HGEN to 2 in *control_energy and use *database_glstat and *database_matsum to report the HG energy for the system and for each part, respectively. The point is to confirm that the nonphysical HG energy is small relative to peak internal energy for each part (<10% as a rule-of-thumb). For shells only, you can fringe hourglass energy density by first setting SHGE=2 in the LS-DYNA input deck (*database_extent_binary). Then, in LS-Prepost, choose Fcomp > Misc > hourglass energy. For solids only, see HYDRO=4 in *database_extent_binary. (r99645, 7/23/15). For fluid parts, the default HG coefficient is generally inappropriate (too high). Thus for fluids, the hourglass coefficient should generally be scaled back several orders of magnitude. Use only viscosity-based HG control for fluids. The default HG formulation (type 1) is generally ok for fluids. Please note that in 971 R3, the default hourglass coefficient for ALE parts (ELFORM 11) is 1.e-6. To override that default, as might be appropriate for non-fluid materials, use *HOURGLASS and HGID in *PART. Check the hourglass energy via matsum. Stiffness-based HG control (types 4,5) is generally more effective than viscous HG control for structural parts. Usually, when stiffness-based HG control is invoked, I like to reduce the HG coefficient, usually in the range of 0.03 to 0.05, soas to minimize nonphysical stiffening of the response and at the same time effectively inhibiting hourglass modes. For high velocity impacts, viscosity-based HG control (types 1,2,3) is recommended even for solid/structural parts. NOTE: Quasi-static simulations may be considered as low speed for which stiffness based hourglass control is recommended. Problems involving high explosives may be considered as high speed for which viscous-based hourglass control is recommended. Impact velocities typical of automotive crash would tend toward the low speed end of the spectrum and so stiffness-based hourglass control is commonly used, but see the note above about "Guidelines provided herein". Type 8 HG control applies only to shell formulation 16. This HG type activates warping stiffness in type 16 shells so that warping of the element does not degrade the solution. Type 16 shells will solve the so-called Twisted Beam problem correctly if HG type 8 is invoked. Type 6 HG control invokes an assumed-strain co-rotational formulation for type 1 solid elements and under-integrated 2D solids (shell types 13 and 15). With the HG type set to 6 and the hourglass coefficient set to 1.0, an elastic part need only be modeled with a single type 1 solid through its thickness to achieve the exact bending stiffness. Type 6 HG control is automatically invoked for type 1 solids in implicit simulations. Also applicable to type 1 hex elements are: - Hourglass type 9 (Puso; similar to but more advanced than type 6) - Hourglass type 10 (Cosserat Point Element; see Remark under *CONTROL_HOURGLASS) The hourglass coefficient for type 6 HG control will typically range from 0.1 (default) to 1.0. For elastic material, use 1.0. For other materials, the choice of HG coefficient is not obvious. Even looking at results, it may be difficult to quantify the 'goodness' of the hourglass coefficient used. Too low a value may result in visible hourglass modes of deformation (unlikely). Too high a value may result in overly stiff behavior. It may be necessary to run the model twice to see if the results exhibit any sensitivity to the hourglass coefficient. Checking the hourglass energy is a good idea. The default hourglass coefficient of 0.1 is superseded by any nonzero value given for QH in *control_hourglass. I see nothing in the manual to contradict this interpretation. The manual does say that the default hourglass type in *hourglass is 1 regardless of what's given in *control_hourglass. Unless I missed something, no such note appears with regard to hourglass coefficient. The lesson here is that users should specify a nonzero hourglass coefficient wherever *hourglass is used. Otherwise, the user may, as you did, inadvertently change the intended coefficient by use of *control_hourglass. _______________________________________________________________________ RE: Viscous damping coefficient (VDC) in *HOURGLASS Note that HG types 6 and 7 have an optional damping coefficient VDC in *hourglass. For 20% damping, set VDC=0.2. In cases documented in bugs 10438 and 13980, tshell type 5 benefited from VDC=0.2 and VDC=1.0, resp. Lee said, "The hourglass control for thick shell 5 is a modified version of hourglass form 6." and "For bricks and also thick shell 5, VDC is a viscous damping coefficient.".