Negative Volumes in Foams (or other soft materials) In materials that undergo extremely large deformations, such as soft foams, an element may become so distorted that the volume of the element is calculated as negative. This may occur without the material reaching a failure criterion. There is an inherent limit to how much deformation a Lagrangian mesh can accommodate without some sort of mesh smoothing or remeshing taking place. A negative volume calculation in LS-DYNA will cause the calculation to terminate unless ERODE in *control_timestep is set to 1 and DTMIN in *control_termination is set to any nonzero value in which case the offending element is deleted and the calculation continues. See http://ftp.lstc.com/anonymous/outgoing/support/FAQ_kw/negvol.k. See also PSFAIL in *control_solid. If PSFAIL is nonzero, ERODE is not used, rather, solid elements of part set PSFAIL will erode based on negative volume. Furthermore, all solids elements will be subject to time-step-based erosion based on DTMIN. See http://ftp.lstc.com/anonymous/outgoing/support/FAQ_kw/erode_solids_by_negvol.k and http://ftp.lstc.com/anonymous/outgoing/support/FAQ_kw/erode_solids_by_tsmin.k . These two input decks require the word "case" on the execution line and run different combinations of ERODE and PSFAIL. Some approaches that can help to overcome negative volumes include the following. - Stiffen up the material stress-strain curve at large strains. This approach can be quite effective. - Sometimes tailoring the initial mesh to accomodate a particular deformation field will prevent formation of negative volumes. Again, negative volumes are generally only an issue for very severe deformation problems and typically occur only in soft materials like foam. - Reduce the timestep scale factor. The default of 0.9 may not be sufficient to prevent numerical instabilities. - Avoid fully-integrated solids (e.g., formulations 2,3,-1,-2) which tend to be less stable in situations involving large deformation or distortion. (The fully integrated element is less robust than a 1-point element when deformation is large because a negative Jacobian can occur at one of the integration points while the element as a whole maintains a positive volume. The calculation with fully integrated element will therefore terminate with a negative Jacobian much sooner than will a 1-point element. (lpb)) - Use the default element formulation (1 point solid) with type 4 or 5 hourglass control (will stiffen response). A better choice of hourglass control for foams modeled with type 1 solids may be: type 6 HG control with coef. = 0.5 if low velocity impact, types 2 or 3 HG control if high velcocity impact. The hourglass type and hourglass coefficient may warrant modification based on observed hourglass modes and on reported hourglass energy in matsum. - Model the foam with tetrahedral elements using solid element formulation 10 (see http://ftp.lstc.com/anonymous/outgoing/support/FAQ_docs/dubois-foam-tets.pdf). - Increase the DAMP parameter (foam model 57) to the maximum recommended value of 0.5. - Use optional card B of *contact to turn shooting node logic off for contacts involving foam. - Use *contact_interior. A part set defines the parts to be treated by contact_interior. Attribute 4 (DA4 = 5th field of Card 1) of the part set defines the TYPE of contact_interior used. The default TYPE is 1 which is recommended for uniform compression. In version 970, solid formulation 1 elements can be assigned TYPE=2 which treats combined modes of shear and compression. - If mat_126 is used, try ELFORM = 0. - Try EFG formulation (*section_solid_EFG). Use only where deformations are severe as this formulation is very expensive. Use only with hex elements. - An ALE approach is often a preferred approach for modeling fluids or solids undergoing large deformations. - An SPH approach may be a viable alternative, although SPH is generally not recommended for situations where tensile behavior is important. *********************************************************** See also: instability.tips ale.basic morefaq's more_on_negative_volume_in_foams <<< maybe should fold this in directly to FAQ jpd