Our website uses cookies. By using the website you agree ot its use. More information can be found in our privacy policy.

LS-DYNA R10.2.0 (R10.135267) released

New version of LS-DYNA is released for all common platforms.

Release notes for LS-DYNA R10.2.0

Herein are summarized new features and enhancements in version 10.2.0.  Some bug fixes are also described.

The items are arranged by category.  Understand that in many cases, a particular item could fall under more than one category, but in the interest of brevity, each item is listed only once, under a single category.

Excluding the "Miscellaneous" category, the categories are arranged alphabetically.


    • Improve p2f and p2p algorithms. In general CPM is speeded up by 3x.

    • OpenMP enabled.

    • Support autoliv mat_34 VOPT=-1/-2 option.

    • Impulse decay scale factor can be controlled by each bag.

    • If the internal parts are solids, the volume from those parts is removed from bag volume.

    • External drag CD_ext can be applied as time dependent curve.

    • New option allows energy exchange between gas and surrounding structure parts. To turn on the feature, provide Cp of target parts. The initial temperature of the parts is same as the resident air (not ambient temperature). The time history of part temperature is reported to abstat_cpm under part_temp

    • If chamber definition is used, the code will error terminate when inflator nodes are not not associated with any chamber ID.


  • Control Volume (Uniform Pressure) *AIRBAG types:

    • Fix bug in *AIRBAG_HYBRID_JETTING with vent coefficient defined with a load curve.

    • *AIRBAG_INTERACTION with LCID will error terminate if used with *AIRBAG_LOAD_CURVE or *AIRBAG_LINEAR_FLUID.



  • *DATABASE_BINARY_FSILNK, *LOAD_SEGMENT_FSILNK: Fix a bug in the routine reading the binary file that created empty files with random numbers appended to their names.

  • *ALE_REFERENCE_SYSTEM_GROUP: If PRTRYPE=9 (the ALE mesh can expand and move with a structure), correct truncation errors in the moving mesh logic

  • *INITIAL_VOLUME_FRACTION_GEOMETRY: If VX or VY or VZ or VR > 0 (initialize the velocities for an ALE group), fix an issue related to the order in which the velocities of the ALE groups are initialized.

  • *CONTROL_ALE: If DCT>-1, add warning messag (KEY+1199) to advise the user to set DCT=-1



  • Fix old bug in *INTERFACE_LINKING that would cause incorrect behavior for displacements scaled via a *DEFINE_FUNCTION if (and only if) the function was defined in terms of 3 variables (spatial displacement dependence only).

  • Fix SMP *INTERFACE_LINKING which was broken in r109945.

  • Corrected the local rigid body rotation in implicit prescribed by *BOUNDARY_PRESCRIBED_MOTION_RIGID_LOCAL which could be wrong when combined with translational motion.

  • Fix processing of prescribed motion constraints for explict computation with superelements (*ELEMENT_DIRECT_MATRIX_INPUT).

  • Add a more robust solution for using birth and death times for prescribed motion constraints when using explicit LaGrange Multiplier Formulation for joints and rigid bodies.

  • Fix seg fault when using *INTERFACE_LINKING_NODE_SET.

  • Fix *BOUNDARY_PRESCRIBED_ORIENTATION_RIGID_DIRCOS which does not orientate the rigid body when angles are close to PI when using double precision.


    • Fix segmentation fault.

    • Fix adaptive nodes which were not constrained.

    • Allow for more than 2 symmetry planes to be defined.


CESE (Compressible Fluid Solver)

  • Expanded a scratch memory array in the CESE 3D FSI-ibm solver in order to avoid issues in very rare cases.

  • For the *CESE_DATABASE cards, corrected these issues:

    • Any used segments that are not in the CESE mesh are reported.

    • Element sets look up is now guaranteed to work.

    • Drag calculation output that is attempted for the incorrect FSI method, is now reported. Here are the illegal cases:

      • *CESE_DATABASE_FSIDRAG is only available with the CESE immersed boundary FSI solvers.

      • *CESE_DATABASE_SSETDRAG is not available with the CESE immersed boundary FSI solvers. In other words, it only works with the non-FSI solvers and moving mesh FSI solvers.

    • Failure to find an internal node number in *CESE_DATABASE commands is now reported.

    • Output was not being performed for the initial conditions. Instead, it is now done after the first time step.



  • Fix possible hanging of SMP single precision when interforce file output (*DATABASE_BINARY_INTFOR) is requested.

  • Undo r108248, which was not a good idea. "Interface Pressure" in the intfor file is now returned to always being positive.

  • Fix initialization error related to *CONTACT_AUTOMATIC_GENERAL_INTERIOR that resulted in every subsequent *CONTACT_AUTOMATIC_GENERAL also being treated as INTERIOR.

  • Fix for spotwelds improperly deleted due to rigid body conflict when an IPBACK contact interface is in effect. MPP only.

  • Fix to recompute rotational node mass increment for BEAM_OFFSET contact after full deck restart.

  • Corrected the collection of resultant forces in implicit for tied contact using the OFFSET option. Other tied contact options such as CONSTRAINED_OFFSET were correct.

  • Made correction to bndout output, including that from prescribed motion, for Implicit Mechanics.

  • Account for inertia relief modes in MPP for the computation of resulting forces from constraints for Implicit Mechanics.

  • Ignore spotwelds with inactive parts during implicit dynamic relaxation computations.

  • Implement consistency, i.e. ncpu<0, when using *CONTACT_SURFACE_TO_SURFACE_THERMAL_FRICTION where static and dynamic coefficients of friction are function of temperature.

  • Fix incorrect rigidwall forces & stonewall energy which will result in incorrect rigidwall frictional behavior when using *RIGIDWALL_PLANAR_ORTHO.

  • Fix *CONTACT_AUTOMATIC_GENERAL for spotweld beams when using SSID=0, i.e. all parts included in the contact, and CPARM8=2.

  • Implement unloading curve, UNLCID, for options FCM=2/3 in *CONTACT_RIGID_(OPTION).

  • Fix zero frictional energy output to glstat and sleout when using *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_ORTHO_FRICTION.

  • Fix input error for duplicate part if multiple *RIGIDWALL_GEOMETRIC_..._DISPLAY are used with some with specified part id, PID, while some are not.

  • Enabled entity contact (*CONTACT_ENTITY) to work with thick shell elements. It previously error terminated in the input phase.

  • Fixed *CONTACT_2D_AUTOMATIC_TIED when used with selective mass scaling, which is made active by setting IMSL=1 and DT2MS<0 on *CONTROL_CONTACT. The problem occcured when selective mass scaleing was applied to the nodes that are tied. Prior to the change, kinematic constraints were failing.

  • Enabled *CONTACT_2D_AUTOMATIC to work with element erosion when using cohesive element forms 29 and -29.

  • Enabled *CONTACT_2D_AUTOMATIC to work with element erosion when using material types 81 or 82.

  • Fixed MPP *CONTACT_2D_AUTOMATIC when beams were in contact, or when solids were used with nonzero values of SLDSOS or SLDSOM. Both cases could lead to using invalid memory space and a segmentation fault.

  • Added support for death time in *CONTACT_2D_TIED. When the death time is reached, the ties are released. To prevent penetration of the released nodes, add another contact definition such as *CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE with a birth time equal to the death time of the tied contact.

  • Enabled multiple instances of *CONTACT_2D_AUTOMATIC_TIED to work correctly. Runs were error terminating.

  • Fixed an MPP segment based (SOFT=2) contact error that could occur if the contact uses SBOXID to eliminate slave segments. This could cause a segmentation fault during initialization.

  • Fixed MPP implicit segment based contact which was failing to properly reset the segment and node data leading to possible "floating invalid" error or possible convergence failure.

  • Fixed an error in MPP segment based eroding contact that could cause force spikes and inappropriate contact detection after solid elements erode.

  • Improve spotweld beams in MORTAR contact, by making sure the end is not subject to contact when already tied or merged to another element, thus avoiding large spurious penetrations and bad convergence in implicit.

  • Fixed interface force file (*DATABASE_BINARY_INTFOR) for SMP version if the file size is larger than 7M byte.



  • Correct the MPP communication during the testing of co-linear nodes in *CONSTRAINED_INTERPOLATION. Promote linear algebra processing for explicit use of *CONSTRAINED_INTERPOLATION to always use double precision even in the single precision executables in order to remove round-off error accumulation.
  • Add additional error checking on coincident nodes for REVOLUTE and SCREW joints.

  • Implement handling of failure of *CONSTRAINE_TIE_BREAK for implicit. This required saving of the failure flag to use when constraint matrix structure has to be the same as the last analyze phase.

  • Fix kinetic energy dependence on cpu count when using *CONSTRAINED_JOINT_(OPTION) with LMF=1 in *CONTROL_RIGID for MPP.

  • Several updates for *CONSTRAINED_INTERPOLATION_SPOTWELD (new card 5):

    • Allow definition of separate stiffnesses STIFF, STIFF2, and STIFF3 for tension, shear, and bending.

    • Optional exponential damage law via LCDEXP.

    • Alternative shear kinematics treatment, variable SROPT.

    • Output of history variables with NEIPB=7 on *DATABASE_EXTENT_BINARY.


  • *CONSTRAINED_COORDINATE: Fix a machine dependent error when the defined coordinates are far away from the actual part.



    • If NLVL>1 (refine more than one level) and CRITRM>0 (a coarsening criterion exists), fix a condition that prevents the coarsening of a given element if one of its child elements is itself still refined.

    • If CRITRM>0 (a coarsening criterion exists), deactivate contact segments of child elements that become inactive after coarsening.

  • *CONTROL_REFINE_SOLID, *CONTROL_REFINE_ALE: If NLVL>1 (refine more than one level), fix the logic merging elements refined at different levels that only share sides.

  • *CONTROL_ADAPTIVE: If IADPOPT=8 (r-adaptive remeshing) and ELFORM=12 (plane stress elements), map the shell thicknesses from old to new mesh.

  • Improve the cleanup of temporary files during adaptive (*CONTROL_ADAPTIVE) steps.

  • Fix excessive memory growth during adaptive problems, and reduce memory requirements overall.

  • Expand the meaning of *CONTROL_ADAPTIVE's variable CNLA: Avoid simplified boundary description in 2D r-adaptive remeshing (ADPOPT=8) by setting CNLA less than zero, recommended value is CNLA = -110.

  • Fix tolerance in selective mass scaling (*CONTROL_TIMESTEP) to avoid non-convergence in single precision.

  • Fixed bug in reading multiple entries for *CONTROL_REQUIRE_REVISION.

  • Fixed bug in reading *CONTROL_REQUIRE_REVISION if free format is used.


    • Fix bug in models featuring shell adaptivity but which also had beam elements (such as used in drawbeads). The bug caused an error termination.

    • Allow a parameter file name in *INCLUDE.


Discrete Element Method

  • Fix bug when both *DEFINE_DE_MASSFLOW_PLANE and *DATABASE_BINARY_DEMFOR are defined.
  • Fix output bug if multiple *PARTICLE_BLAST are defined ---MPP only.


EFG (Element Free Galerkin)

  • Fixed bug in reading mutiple *SECTION_SHELL_EFG.

  • 3D adaptivity EFG now supports output associated with variables STRFLG and INTOUT in *DATABASE_EXTENT_BINARY.

  • EFG shell (ELFORM=41) parts are distributed to all CPUs by default in MPP.

  • Fixed incorrect stresses in EFG plane strain formulation 43 for both implicit and explicit analyses and EFG axisymmetric solid formulation 44 for explicit analysis only.



  • Fix a MPP bug for type 9 *ELEMENT_SEATBELT_PRETENSIONER that could be triggered by having more than one type 9 pretensioner.

  • Fixed mass and kinetic energy output in d3hsp and matsum files for type 8 (axisymmetric) beams.

  • Enabled thick shell form 7 to output data for 4 points per layer to the ELOUTDET file.

  • Added support of the NINTSLD variable on *DATABASE_EXTENT_BINARY for solid pentehedral element forms 15 and 115. Also, when NINTSLD=1, form 15 now averages the data from the two points for output rather than printing the data of the first point.

  • Fixed bulk viscosity of thick shell forms 1, 2, and 6. Unstable behavior was likely.

  • Fixed the behavior of CDL and TDL on *SECTION_DISCRETE when one of the nodes is constrained by an SPC. Also, corrected the DEFORC output when CDL and TDL limits are reached.

  • Fixed thick shell forms 5 and 7 when used with total Lagrangian material models such as *MAT_002. The thick shells use a corotational coordinate system which rotates when the element rotates. However, the calculation of the deformation gradient, F for the strain measure also includes this rotation, so we were accounting for the rotation twice. To fix this, the thick shells now calculate F using the symmetric rate of deformation, D, rather than the velocity gradient, L. With this change, large rotation of elements is possible and the strain and stress remain accurate.

  • *SECTION_SOLID, ELFORM=19,20,21,22: Suppress unnecessary writing of warning messages concerning element type switching for cohesive pentas. This change does not affect the analysis!

  • Triangular shell element type 4 that combines thermal expansion and thick thermal shells, see *MAT_ADD_THERMAL_EXPANSION and TSHELL on *CONTROL_SHELL, now correctly treat temperature gradient through the thickness to create bending moments. This is already available for shell element types 2 and 16 since R10.0.0.

  • Fix *ELEMENT_BEAM_PID when offsets from *INCLUDE_TRANSFORM are used for PID1 and PID2.

  • Fix seg fault in Cosserat element, solid formulation 1 or 16 with hourglass type 10, when hourglass energy is not asked for, and fix hourglass energy output for when it is asked for.

  • Fix spurious read error for *SECTION_BEAM_AISC.

  • Fix error in reading Card 5 of *SECTION_SHELL for user-defined shell elements.


Forming Analysis

    • Fix bug for trimming with a vector (option _2D), which was not recognized.

    • Fix bug in missing strain tensors after trimming.

    • Fix bug for "Normal Termination" not written into mes0000.

    • Fix bug for trimming of tshell with a defined vector, which was not recognized.

    • Fix bug for trimming of sandwiched parts with a defined vector, which was not recognized.

    • Fix bug for some lost history information during trimming of shells.

    • Fix bug in 3D trimming, where specifying a vector which is not needed causes an error termination.

    • Fix bug in trimming of multi-layer sandwiched sheets, where it would generate unnecessary shells at the top and bottom of every solid.

    • Fix stress output issue.


    • Fix bug when the smallest distance in the curve referenced by *CONTROL_FORMING_OUTPUT is equal to the total travel distance of the tool, LS-DYNA outputs d3plot files every time step.

    • Fix bug where only the 2nd NOUT (forming) is recognized and 1st NOUT (closing) is missed.

  • *ELEMENT_BLANKING: Fix bug of not filling the corner with elements properly.

  • *CONTROL_FORMING_REMOVE_ADAPTIVE_CONSTRAINTS: Fix bug for trianglar elements whose adaptive constraints were not removed.


    • Type=2: fix bug for implicit calcuation convergence issue; the problem was caused by truncation error causing solution to oscillate leading to non-convergence.

    • Type=1: fix truncation error which caused solution to diverge.

    • Type=1: fix energy explosion and division by zero error, and fix truncation error when blank and binder are far apart.

  • Fix bug where *PART_MOVE with IFSET=1 does not move parts if set is defined using *SET_PART_ADD, and/or *SET_PART_LIST_COLLECT.


    • Output an error message: "Vector: xx was not defined" in case a vector is prescribed to a fixed rigid tool by mistake.

    • Check constraints for the rigid clamper and free the constraints from the moving clampers.

  • *DEFINE_MULTI_DRAWBEADS_IGES: Fix duplicate node set problem, which caused the simulation to crash.



  •  *FREQUENCY_DOMAIN_RANDOM_VIBRATION: Fixed a bug in computing RMS (Root Mean Square) values of Von Mises stress for solid elements.
  • *FREQUENCY_DOMAIN_RESPONSE_SPECTRUM: Fixed nodal displacement output to d3spcm. Only the nodal displacement (u) itself (not updated nodal coordinates, which is u+x0) is written to d3spcm.

  • *FREQUENCY_DOMAIN_SSD: Fixed dumping nodout_ssd for the load cases of base velocity (vad=2) and base acceleration (vad=3), when local damping is used (dmpflg=1)


ICFD (Incompressible Fluid Solver)

  • Fixed HCSFLCID and TCSFLCID issue in *ICFD_MAT. Those two optional load curve IDs were ignored in MPP.

  • Fixed erroneous stabilization of temperature field. Results should now be close again to R9 and R11 versions.

  • Brought in R11 fix for multiphase cases (eg air+water problem) which should now be more accurate and stable.

  • Changed the format of ICFD_BOUNDARY_FSWAVE keyword to reflect the subsequent changes in R11 and Dev.


Implicit (Mechanical) Solver

  • Corrected Multistep linear (NSOLVR=-1) for the case of just one step.

  • Properly size an array used during the computation of participation factors associated with eigenmodes for the very small model case.

  • Assure that Implicit Mechanics is turned off for thermal only solutions.

  • Enhanced the MMD ordering for Implicit Mechanics ISOLVR=2 option.

  • Corrected an inefficiency in Implicit Mechanics when there are numerous AutoSPC constraints in MPP. A linear search through the AutoSPC constraints was replaced with a binary search.

  • Fix divide by zero in power iteration for buckling problems with inertia relief.

  • Fix to properly handle an exactly singular constraint matrix for *CONSTRAINED_INTERPOLATION running in fast/vector mode.

  • Extend logic for implicit constraint handling to tied contact with CONSTRAINED_OFFSET to skip any SPC constraints on the rotational dofs of the slave nodes. We already have such logic to skip the SPC constraints on the translational dofs for the slave nodes for all tied contact.

  • Restrict a block of code executed during *PART_MODES initialization to only run on Processor 0 in MPP.

  • Enhance the logic for dumping matrices in the presence of rotational dynamics.

  • Correct the issue where use of *LOAD_BODY is applied to a model with rigid bodies. If the vector of nodes having the load applied belong to a rigid body then the vector of elemental stiffness matrices is null. Added logic to skip the call to the implicit matrix assembly subroutine for that case.

  • Solution control was unexpectedly returning to implicit after a "sw1" sense switch. This termination is now cleanly handled.

  • Correct the dumping of matrices from implicit when MTXDMP > 1.

  • Added an error test for improper DRTERM when using implicit for dynamic relaxation. Enhanced implicit key point logic during dynamic relaxation phase.

  • Enable Implicit to collect resultant forces for SPC constraints in local coordinates.

  • Fixed implicit element stiffness of shell element when used with laminated shell theory, and when also using anisotropic material models 22, 34, 40, 54, 55, or 58. It was likely that the shear stiffness was too small.

  • Fixed implicit treatment of beam spot weld elements (*SECTION_BEAM ELFORM=9) when used with NF>0 on *MAT_SPOTWELD. The averaging of the stress was causing the model to converge to the wrong solution with beam forces too low.

  • Fix stiffness matrix for Belytschko-Tsay shell to yield proper rigid body modes in an eigenvalue analysis.

  • Fix bug that made rigid bodies deform in implicit-explicit switching problems.



  • Fix seg fault or incorrect stresses when initializing stresses using *INITIAL_STRESS_SOLID for *MAT_107/*MAT_MODIFIED_JOHNSON_COOK.
  • Fixed *INITIAL_STRESS_SHELL and *INITIAL_STRESS_TSHELL when used with *INTEGRATION_SHELL. The integration rule was getting lost leading to unecessary interpolation of data.

  • Fix beam ID in message file when *INITIAL_AXIAL_FORCE_BEAM preload is completed.


Isogeometric Elements

  •  *ELEMENT_SHELL_NURBS_PATCH: Fixed bug in reduced integration rule INT=2 that appeared if continuity within a quadratic patch may have been reduced to C0 locally due to duplicated knot-vector-values.



  • Fix output to nodfor when using *LOAD_SSA.



  • Fixed bug for *LOAD_SEGMENT_FILE when jobid is used.


*MAT and *EOS

  • Fixed bug in *MAT_203 / *MAT_HYSTERETIC_REINFORCEMENT affecting beam elements with EPDAM1, EPDAM2 non-zero, and DRESID equal to zero. When the beam reached its erosion limit, error termination could occur or there could be unexpected failure of the whole model.
  • Fixed bug in *MAT_211 / *MAT_SPR_JLR. Load curve IDs 8 or 9 digits long were not read correctly from *MAT card in single precision executables.

  • Bug fix to resolve a divergence issue when using *MAT_242 / *MAT_KINEMATIC_HARDENING_BARLAT2000.

  • Fix ineffective load curves LCSRA LCSRB LCSRC LCSRAB LCSRBC LCSCA in card 7 of *MAT_MODIFIED_HONEYCOMB / *MAT_126.

  • Fix ineffective FAIL<0 flag for *MAT_MODIFIED_PIECEWISE_LINEAR_PLASTICITY / *MAT_123, i.e. using user-defined failure subroutine, matusr_24 in dyn21.F, to determine failure.

  • Fix output to elout for solids using *MAT_PAPER / *MAT_274 and CMPFLG>0 in *DATABASE_EXTENT_BINARY.

  • Fix seg fault when using *MAT_ADD_EROSION for *PART_COMPOSITE when the first layer of the *PART_COMPOSITE is not using *MAT_ADD_EROSION.

  • Fix incorrect direction of thermal expansion when using *MAT_NONLINEAR_ORTHOTROPIC / *MAT_040 with *MAT_ADD_THERMAL_EXPANSION.

  • Fix contact penetration after element erosion when using *MAT_110 / *MAT_JOHNSON_HOLMQUIST_CERAMICS with eroding contact.

  • Fix incorrect stress/strain output to d3plot when using *MAT_TEMPERATURE_DEPENDENT_ORTHOTROPIC / *MAT_023 and running 2D analysis and MAXINT>3 in *DATABASE_EXTENT_BINARY.

  • Fix incorrect stress initialization when using *MAT_005 / *MAT_SOIL_AND_FOAM with LCID (pressure vs. volumetric strain), combined with *LOAD_DENSITY_DEPTH.

  • Fix strain rate effects on *MAT_157 / *MAT_ANISOTROPIC_ELASTIC_PLASTIC for implicit static analysis.

  • Fixed *MAT_128 / *MAT_HEART_TISSUE for thick shell forms 3, 5, and 7. Initializaion of the F tensor was missing causing termination in the first cycle.

  • Fixed use of *MAT_FABRIC form -14 with Rayleigh damping. Runs were failing in the first cycle.

  • Improved the behavior of hourglass form 6 with used with *MAT_089. The material parameter used to scale the hourglass stiffness was wrong leading to bad results.

  • Fixed *MAT_126 when used with type 3, 5, or 7 thick shells. The angle initialization was incorrect leading to bad stress values.


    • Fixed bug when using *MAT_058 with load curves defined for EA, EB, GAB together with thick shell elements type=1 and 2. Prior to this fix, the analysis with this combination may have gone unstable.


    • Add flag "DSF" to *MAT_262 (Card 3, Column 7) that controls the failure of an integration point based on in-plane shear failure. This feature was added to reinstate the behavior of this material model prior to a fix made in r117876.

    • Note: The behavior prior to r117876 was not really correct due to the theory of the model, but customers adjusted their model with this setting. Using the new DSF flag allows to enforce the same behavior as before.


    • Fix bug in *MAT_261 / *MAT_262 (shells) that may have appeared if this material was used in *PART_COMPOSITE with a different material model located at first integration point (layer). Prior to this fix, the deletion of the layered element was not triggered correctly.


    • Fix bug when using table (2D) for defining strain rate dependent fracture toughness values with respect to characteristic element length. Prior to this fix, a wrong value for the fracture toughness could have been used.


    • In case EPSR and EPSF are defined, correct computation of transverse shear strains for solids and tshells.


    • Added a criterion to avoid possible snapback behavior in *MAT_262 by only allowing certain "softening" modulus with respect to the elastic stiffness. This might have lead to problems when using rate-dependent strength limits together with fracture toughnesses values that are not properly set.

  • Bug fixes for *MAT_264 that resolve issues with AWG Test Case 16 (https://awg.lstc.com/tiki-index.php?page=AWG+ERIF+Test+Example+16).

  • Fix NaN problem with shells using *MAT_024_LOG_INTERPOLATION with VP=1 and dynamic relaxation.

  • Improve algorithm for the LOG_INTERPOLATION option of *MAT_024 when used with VP=1 for shell elements. CPU time is reduced and results better match the alternative approach for log interpolation (negative first value in table).

  • Fix for combination of *MAT_ADD_EROSION and thick shell type 2. Strain-based criteria (e.g. MXEPS) did not work correctly before.

  • Fix for combination of *MAT_ADD_EROSION and *PART_COMPOSITE by clearly detecting if erosion/damage is used for a layer.

  • Add failure function terms (normal, bending, shear) as arguments of functions in *DEFINE_CONNECTION_PROPERTIES with PRUL.ge.2.

  • Fix for GISSMO to be used in adaptive remeshing. Evaluation of damage coupling flag could go wrong due to averaging of history variables in rezone step.



  • Fix for failure treatment in *MAT_224_GYS: Update of failure function F could have been wrong, because the plastic strain increment was not always zero for elastic integration points.

  • Fix bug in *MAT_ADD_EROSION, DIEM option, that would otherwise result in spurious damage when more than damage model is used (IDAM less than -1)

  • Fix *MAT_ADD_EROSION, DIEM option, in combination with *INTEGRATION_SHELL that would otherwise not work.

  • Support logarithmic strain rates in tables for *MAT_ADD_EROSION, DIEM option, activated by having the first strain rate negative.

  • Fix incorrect internal energy calculations in *MAT_ADD_THERMAL_EXPANSION.

  • Change the stress update in *MAT_PML_ELASTIC_FLUID to make it stable.

  • A first-order low-pass filter is implemented to smooth the effective strain rate used for MXEPS<0 in *MAT_ADD_EROSION. All high frequency response above the cutoff frequency is filtered out. The variable for invoking this filter is DTEFLT.

  • *MAT226_NLP:

    • Fix truncation error problems which caused the simulation to crash.

    • Add a lower limit for strain ratio after smoothing, which will be the the strain ratio of the first point on the FLD curve input.

  • *MAT_036_NLP: Bug fix for non-linear strain path forming limit; reset beta value to zero.



  • Fix problem with "%proc" option in the region based decomposition in the pfile, which was always putting the indicated region onto a range of processors starting at processor 0, even if that would result in poor load balance.

  • Improve MPP decomposition processor assignment for regions assigned to a subset of the processors.

  • Fix for possible incorrect synchronization of shared nodal coordinates after initial penetration checking in MPP contact.

  • Fix missing initialization that broke d3plot files when MPP predecomposition was used with models having *DEFINE_CURVE_FUNCTIONS that used the PIDCTL function.

  • Fix for MULTISCALE_SPOTWELD in double precision

  • Support Intel MPI 2018 on Linux Cent OS 6.5 and later.

  • AVX2 MPP binaries should now get identical results from AMD and Intel chip.

  • A listing of part sets defined in the keyword input is output by LS_DYNA to an ASCII file called group_file. To display the parts in any part set, read the keyword input into LS-PrePost and then go to "Model>Groups". group_file is automatically loaded and a list of part sets appears. You can highlight any of those part sets and click "Select" to display the parts in those selected part set(s).

  • Output user defined *PART_SET to group_file with title.

  • Collect and print massless nodes from all processors to mes0000 and d3hsp files.

  • Report time history of compute node memory, max/min memory to binary glstat for performance analysis.

  • MPP option VSPEED will test and report network latency for performance debugging.

  • *PART_AVERAGED, which pertains to truss elements, is enhanced for the better MPP scalability without special MPP decomposition.

  • A new option "together" for MPP decompostion is implemented to keep a region together on one processor, as a super block as "lumped". It is also kept for the decomposition with other parts to reduce the MPP communication.



  • Corrected internal energy output for *DAMPING_FREQUENCY_RANGE_DEFORM. The bug affected solid elements only, and related only to the output values of internal energy. The solution itself was not affected, but users may have been concerned that the energy did not seem to balance.

  • Corrected internal energy and external work outputs for staged construction (*CONTROL_STAGED_CONSTRUCTION, *DEFINE_STAGED_CONSTRUCTION_PART). The bug affected only the output values of internal energy and external work. The solution itself was not affected by the bug, but users may have been concerned that the energy did not seem to balance.

  • Corrected internal energy output for pore fluid (*CONTROL_PORE_FLUID) with multi-integration-point solid elements such as ELFORM=2. The bug affected only the output values of internal energy. The solution itself was not affected, but users may have been concerned that the energy did not seem to balance.

  • Fix bug, the effect of which is not obvious, but could possibly have led to incorrect output of beam data in elout.

  • Fix a MPP bug for nodal stress/strain output to eloutdet (see variable NODOUT in *DATABASE_EXTENT_BINARY), which could occur when multiple parts share the same node where stress/strain output is requested.

  • Fix a bug for 2d seatbelts whereby incorrect ID of 2d seatbelt retractor was output to binary sbtout data.

  • Improved the reporting of hourglass energy to the glstat and matsum files. Hourglass energy of eroded bricks was added to the eroded hourglass energy without removing it from the regular hourglass energy. Also, thick shells were not reporting eroded hourglass energy. Fixing these issues does not affect the solution, but improves the energy balance reported to glstat.

  • Fixed a memory error that occurred when there were spot weld assemblies that failed and the input file did not include the keyword *DATABASE_SWFORC.

  • Fixed bug for d3plot output if both high order and ALE elements are present in a model.

  • Echo *DATABASE_EXTENT_BINARY_COMP flags in d3hsp.

  • Fixed bug for d3plot if DCOMP=5 or 6 in *DATABASE_EXTENT_BINARY and PSETID in *DATABASE_BINARY_D3PLOT is used.

  • Fixed bug for *DATABASE_HISTORY_NODE_ID in free format if the heading is more than 10 characters.

  • Fixed part id output for discrete element in d3hsp.



  • Fix MPP full deck restart issue that caused hanging if the new input file contained encrypted input.
  • Corrected old memory pointers to work correctly in implicit on restart. Also correct dump and restart length for implicit common blocks to enable new nonlinear solution features.

  • Fix seg fault when using *DELETE_CONTACT for restart when running with SMP.

  • Fix error termination for full deck restart including *DEFINE_ELEMENT_DEATH.

  • Fixed corruption of d3part caused by a simple restart.

  • Fixed bug for corrupted d3part on small restart.

  • Fix bug in full-deck restart of tied contact so that contact forces are continuous across the restart.

  • *LOAD: Fix bug in accessing ground motion ID (*DEFINE_GROUND_MOTION) from *LOAD_SEISMIC_SSI.



  • Fix problem that the presence of *DATABASE_BNDOUT could cause *SENSORs in the model to misbehave, with unexpected results.

  • Fix a bug for the ELESET option of *SENSOR_CONTROL.

  • Fix a bug for *SENSOR_DEFINE_FORCE, TYPE=JOINT, triggered by having a joint ID > 9999999.

  • Fix a bug for *SENSOR_DEFINE_..._SET that could occur if the *SENSOR command appears before the corresponding *SET command in the input file.


SPG (Smooth Particle Galerkin)

  • Fixed problem of zero stress being output to d3plot for SPG with isotropic materials when CMPFLG=1 in *DATABASE_EXTENT_BINARY.


SPH (Smooth Particle Hydrodynamics)

  •  SPH performance improvement for OpenMP and HYBRID.
  • Support full deck restart for SPH.

  • Improve performance of *DEFINE_ADAPTIVE_SOLID_TO_SPH.

  • Support up to 6 SPH particles in each direction of each solid element in *DEFINE_ADAPTIVE_SOLID_TO_SPH.

  • Identify element type as SPH when echoing part information in d3hsp:
       element type ..................... = 4
                eq.0: 4, 6, 8, 10-node solid element
                eq.1: 2-node beam or truss or 2D shell element
                eq.2: 3, 4-node membrane/shell or 2D continuum element
                eq.3: 8-node thick shell element
                eq.4: SPH element


Thermal Solver

  • Correct logic on acceptable thermal solver options for SMP.

  • Fixed bug where dynain file was missing some thermal history variables when using *INTERFACE_SPRINGBACK_LSDYNA_THICKNESS with NTTHSV>0.

  • Improved reporting of the thermal time step.



  • Fix functions PIDCTL and DELAY in *DEFINE_CURVE_FUNCTION when used in a full restart.

  • Fix a bug for MIRROR option of *DEFINE_TRANSFORMATION that was triggered when used together with other transformation options.

  • Fix wrong pointer used for section id/properties when generating part for visualization of rigid wall during adaptivity. This caused error message, KEY+137, during adaptivity.

  • *DEFINE_FUNCTION: Corrected output of function definition to d3hsp for double precision executables.

  • Fixed segmentation fault when using *DEFINE_BOX_ADAPTIVE in a large model.

  • Fixed bug affecting node sets created by *SET_NODE_ADD.

  • Fixed bug where use of *SET_PART_COLLECT caused parts in PSETID of of *DATABASE_BINARY_D3PLOT not to be excluded from d3plot. This was related to the duplicated ID in *SET_PART_COLLECT. Now, an error message will be issued.

  • Fixed bug in newformat option for *NODE_SCALAR.

  • For newformat=i10 option, "%" sign is only added to the modified keyword.

  • Fixed command line option: plabel=(yes,no). This command (plabel=no) will disable the checking for alpha IDs file to speed up the input phase or disable any bug in the checking of labels.

  • Fixed legend output in nodout file if long format is used.

  • Fixed bug in reading long format if "*KEYWORD long=yes" is used in include file.

  • Fixed bug for *DATABASE_HISTORY_NODE_LOCAL_ID when long=s is used.

  • Fixed bug for reading long format input in the following keywords:

  • *SET_NODE_LIST_GENERATE: Fix bug which resulted in a corrupt node list.