LS-DYNA V971 R5.0 (R5.59419) released

Version 971 R5 includes the following main developments and updates from the previous release R4.2.1

Airbag

• Fix for bug related to 8 digit load curve ID for LCDT in *AIRBAG_WANG_NEFSKE option.
• Extended the number of gases permitted for hybrid inflators from 9 gases to 17.
• For molar-ratio airbag inflator input: correct a bug when calculating original gas number. It could be zero when it is firstly called. That will lead to a “defined before used” bug.
• Fix for bug related to using *AIRBAG_HYBRID in a deck with more than 9999 load curves.
• Fix for large load curve ID’s related to *AIRBAG_SIMPLE_AIRBAG_MODEL.
• Fixed bug – *INCLUDE_TRANSFORM and *AIRBAG_HYBRID_JETTING_CM: Node number NREACT was not offset.
• Add part pressure of alone/away shell normal to abstat_cpm in order to see the pressure on both side from p2f collision.
• Improved air particles thermal equilibrium for smp code.
• Additional development for *AIRBAG_PARTICLE chamber interaction.
• Airbag Corpuscular Particle Method: enhancement in consistency of venting option between SMP and MPP.

• Add *MAT_ALE_INCOMPRESSIBLE.
• Disable node sharing check between ALE/nonALE parts for Steady State Rolling.
• ALE *BOUNDARY_NON_REFLECTING: Calculate the bulk modulus every time step instead of using the initial value. For HE material the bulk modulus becomes hundreds times smaller after detonated so using the initial value would make it like hitting a still wall and hence kill the timestep size.
• Bug fix in ALE tracer particles.
• Bug fix for ALE 2D reset the shared number of history variables for general ALE 3D application.
• Support of ALE IVFG to set initial velocity in MPP.
• Fix ALE Leakage coupling force bug.
• Fix ALE EBC=2 bug.
• Fix for ALE automatic BC. In the first cycle, x was updated before ALE BC was enforced.
• ALE2d: use old search code of the neighbor elements in mapping for windows.
• Bug Fix for ALE 2d: uninitialized variables and mismatched declarations

Boundary Conditions

• *BOUNDARY_PRESCRIBED_FINAL_GEOMETRY:
  Nodes are displaced from the initial geometry to the final geometry along a straight line trajectory.
• *LOAD_SEISMIC_SSI:
  *BOUNDARY_FREE_FIELD_GROUND_MOTION has been renamed to *LOAD_SEISMIC_SSI

• Fixed bug: MCOL restart failed for NCPU > 1. Actually, MCOL never worked for NCPU > 1. This bug is fixed and NCPU >1 works now too for MCOL.
• Suppress SPC printing with NPOPT=1 like nodes as suggested by Len Schwer.
• Add damping energy calculation to 2-D non-reflecting boundary conditions.
• ALE *BOUNDARY_NON_REFLECTING: Calculate the bulk modulus every time step instead of using the initial value. For HE material the bulk modulus becomes hundreds times smaller after detonated so using the initial value would make it like hitting a still wall and hence kill the timestep size.
• Additional changes to new keyword *BOUNDARY_PRESCRIBED_FINAL_GEOMETRY.
• To fix the MCOL restart for both SMP and MPP.
• Added ID and death time to *BOUNDARY_PRESCRIBED_FINAL_GEOMETRY.
• Add support for *BOUNDARY_SPC_BIRTH_DEATH to implicit.
• Fix for bug related to *BOUNDARY_THERMAL_BULKFLOW_SET.
• Default ATYPE on *BOUNDARY_PORE_FLUID should be the user-input ATYPE on *CONTROL_PORE_FLUID. Ensure that this is always so.
• Fix bug that prevented more than one *BOUNDARY_PORE_FLUID_SET cards.
• Added MPP support for *BOUNDARY_PRESCRIBED_FINAL_GEOMETRY.
• Implement *SECTION_BEAM_AISC to incorporate AISC standard steel sections. Rename *BOUNDARY_FREE_FIELD_GROUND_MOTION as *LOAD_SEISMIC_SSI.
• Implementation of *INCLUDE_TRANSFORM to *BOUNDARY_PWP_option.
• Implemented *BOUNDARY_SPH_FLOW.
• A fix for *BOUNDARY_PRESCRIBED_ACCELEROMETER.
• Allow *MAT_160 to constrain nodes in local coordinate systems defined by *BOUNDARY_SPC_NODE.
• Changed keyword *BOUNDARY_ELEMENT_METHOD_ACOUSTIC to *FREQUENCY_DOMAIN_ACOUSTIC_BEM (the old keyword still works).
• Fixed bug for which the time step was being used when imposing boundary conditions was incorrect on the first cycle, so some position and velocity boundary conditions could be wrong on the first cycle, resulting in displacements that were off by about 1 cycle worth of motion.

Constrained Options

• *CONSTRAINED_JOINT_STIFFNESS_…_TITLE
  Allows for a description of the joint which is then written to the D3HSP and JNTFORC files.
• *CONSTRAINT_JOINT_COOR_…
  Defines a joint using coordinates instead of nodes.
• *CONSTRAINED_LAGRANGE_IN_SOLID, CTYPE=4
  Implemented in MPP ALE 2D
• *CONSTRAINED_LAGRANGE_IN_SOLID, CTYPE=11
  Output of flow rate in dbfsi for *CONSTRAINED_LAGRANGE_IN_SOLID with ctype=11
• *CONSTRAINED_SPR3
  New material model for spotwelds in *CONSTRAINED_SPR3 environment.

• Fix bug (MPP restart with *CONSTRAINED_TIED_NODES_FAILURE).
• A bug fix in mesh coarsening with *CONSTRAINT_EXTRA_NODES.
• Fix for 2d adaptivity problem related to *CONSTRAINED_LOCAL.
• Fix masses of shared nodes in case of lumped masses being added via RBE3 constraints.
• Fix the interaction of joints and implicit springback.
• Preventing mass scaling of nodal rigid bodies if selective mass scaling is turned on.
• New material model in *CONSTRAINED_SPR environment.
• Add a new option to *CONSTRAINED_COORDINATES 4,5,6 and 7 for xy,yz,zx,and xyz.
• Added no failure option to *CONSTRAINED_SPR3 model.
• Added second damage value to *CONSTRAINED_SPR3 model.
• Update SPR for option 2 to include spotweld thinning.
• Add SPOTHIN option and symmetric behavior to *CONSTRAINED_SPR3 model.
• For reasons not fully understood, *CONSTRAINED_LOCAL was only being allowed in 2d adaptivity, and 3d solid adaptivity — not adaptivity in general. This restriction has been removed.
• Fix *CONSTRAINED_INTERPOLATION_LOCAL: when the local option is invoked with the local coordinate systems set to 0.
• Added _AUTO_JOINT option to simplify joint input.
• Added translational and rotational mass to the *CONSTRAINED_JOINT_COOR option.
• Added TITLE option to *CONSTRAINED_JOINT_STIFFNESS.

Contact Options

• *CONTACT (Optional Card D)
  Various options for SOFT=2: split quads in trias, time interval btw penetration reports, scale factor of neighbor segment contact Option for tied contacts: incremental displacement.
• *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE
  Enhanced treatment of contact between beams and surfaces.
• *CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR
  To use a single surface mortar contact, implicit (preferred) or explicit analysis, including edge treatment.
• *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR
  To use a surface to surface mortar contact, implicit (preferred) or explicit analysis, including edge treatment.
• *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_MORTAR_TIED
  To use a tied surface to surface mortar contact, implicit (preferred) or explicit analysis, including edge treatment.
• *CONTACT_TIED_NODES_TO_SURFACE_OFFSET, _TIED_SURFACE_TO_SURFACE_OFFSET, _AUTOMATIC_TIEBREAK (options 1,2,3,5)
  The option is turned on by placing a 1 in field 7 of optional contact card D. It results in contact displacements being computed incrementally, which should be more stable. (TIEBREAK options > 5 are already incremental). MPP only.
• *CONTACT_…, SOFT=2, SFNBR>0
  For segment based (SOFT=2) contact, a neighbor segment contact check was added. When solving compression folded airbags, sharp folds can result in neighbor segments coming into contact. The new neighbor segment checking option can reduce penetration between neighbors.  It is turned on by a new parameter, SFNBR in the 3rd field of optional card D of *CONTACT. SFNBR is a scale factor for the stiffness with 0.5 or 1.0 as a reasonable value. Setting SFNBR=0. makes the check inactive. For most applications, folding does not occur and this option is not helpful.
• *CONTACT…_…MPP, GRPABLE
  This flag turns on a new set of MPP contact routines that are more highly parallelized than the traditional routines. If used, it should be used for as many contacts as possible, as the main benefit is in consolidating communication and parallelizing across contacts, and there is no point in having just one groupable contact. MPP only.
• *CONTACT_AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE_TIEBREAK, *DATABASE_INTFOR
  Visualization of delamination via “contact gap” in intfor database: instead of only showing a tied/released flag, now a damage value ranging from 0 to 1 is shown. This is now working for SMP and tiebreak options 6/7/9 (between solids) and 8/10/11 (between shells).

• Added frictional work calculation to 2D automatic contact.
• Fixed problem in *SET_SEGMENT_GENERAL where first part ID is 0 causing all parts to be included in the contact.
• Fix for MPP soft=4 contact: a parameter typo, plus slight difference in thickness treatment in some cases between the bucketsort and normal cycles, together causing some contacts to be missed.
• Bug fix for SOFT=4 forming contact (smp only).
• Fixed segment based contact option to use the time averaged current time step in the stiffness calculation.
• Added an option to geometric contact entities to account for the thickness of the slave nodes. Only for entity types 1-3 so far (plane, sphere, cylinder).
• Fix bug: Thickness updates for the single surface contact option on SMP. The MPP version was working as expected.
• Fix bug: Untied nodes are now printed in *CONTACT_TIEBREAK_SURFACE_TO_SURFACE in SMP needs to be added in MPP version.
• Fix for bug related to contact ID output.
• Add stiffness for implicit in *CONTACT_AUTOMATIC_GENERAL.
• Fix for lack of consistency related to automatic tiebreak option type 8.
• Fix for units transformation (*INCLUDE_TRANSFORM) for *MAT_138 and *CONTACT_AUTOMATIC_…_TIEBREAK options 9 and 11.
• Improvements in _AUTOMATIC_TIEBREAK option 4 with sliding.
• Obsoleting the use of *LOAD_SURFACE_STRESS in conjunction with *CONTACT_…_THERMAL_FRICTION.
• Fix for MPP smooth contact: in some cases, triangles were being improperly extended, resulting in strange contact behavior.
• Fix for MPP single surface contact in case PENOPT=4 is used. This is NOT normal for single surface contact, and appears to have never been used, because it was resulting in 0 forces due to the wrong value being passed in as the # of segments to the stiffness rescaling code.
• Fixed possible memory error when using SMP parallel consistency option with 2D automatic contact.
• Added new option *CONTACT_REFERENCE_GEOMETRY.
• Fixes related to constrained offset contact.
• Fixes for MPP implementation of *CONTACT_CONSTRAINT.
• Neglect birth and death for tied contacts when arclength us used.
• alias *CONTACT_SPOTWELD_*_OFFSET to *CONTACT_TIED_SHELL_EDGE_TO_SURFACE_*_OFFSET
• MPP changes to support _AUTOMATIC_TIEBREAK in implicit.
• Fix for bug related to “remaining thickness” reports from contact. The fix ensures that SMP and MPP give similar results.
• Fix for auto tiebreak contact on SMP. Consistency option failed to work correctly leading to erroneous answers.
• Fix for proper initialization of *CONTACT_AUTOMATIC_TIEBREAK in MPP implicit: nodes were not being tied due to confusion about the definition of “the first cycle”.
• Fix decomposition options for distributing/isolating contact interfaces so they will work in case the slave side has no segments (beam or _NODE_TO_SURFACE type contacts) — previously they only looked at the slave segment list. If that list is empty, it now considers the slave node list.
• Add interface surface energy output for type 4, 13 and 26 contacts in MPP.
• Fixed a neighbor segment/node data problem in MPP segment based contact when the quad to triangle switch flag is set. This fix should improve the behavior of the sliding option and edge-to-edge checking when the quad to triangle option is active.
• Added a neighbor segment contact check to segment based contact. When an airbag folds, neighboring segments may pass through each other if this check is not active. It is turned on by a new parameter, SFNBR in the 3rd field of optional card D of *CONTACT. SFNBR is a scale factor for the stiffness with 0.5 or 1.0 as a reasonable value. Setting SFNBR makes the check inactive.
• New feature for thermal contact with FRICTION: conductance a function of the gap defined by a load curve.
• Change to the way MPP treats “initial penetration tracking” in contact. Previously it was generating energy (it was showing up as negative sliding energy).
• In segment based contact, support the 2 surface force transducer option to add forces to all matching Force-Transducers instead of only first match. This is only useful if 2-surface force transducers overlap on both the master and slave side.
• Fix bug: incorrect behavior of tied-offset contact under full-deck restart. After the restart, the contact was being (re-)initialized using the coordinates at the end of the initial run, which caused the nodes to be attached to different segments than they were to in the initial run, thus  changing the contact stiffnesses.
• Add soft=1 option to old node to surface eroding contact options.
• Add MPP support for non-linear stiffness scaling to type 3/5/10 contacts, as well as the groupable versions of 3/5/10 and 13 contacts.
• Fix for remapping of master force histories in MPP eroding contact: due to a variable name typo, some might have been missed, or a segfault could have occurred.
• Add damping in mortar contact.
• Changing mpp implementation of implicit tied offset contact to behave similar to smp.
• Fix for implicit tied beam offset.
• MPP edge to surface contact: forces were being applied to the beam nodes with the wrong weights (they were swapped), causing incorrect behavior and energy issues.
• Added new contact option for SURFACE_TO_SURFACE contacts: On *CONTACT optional card D, in field 2, the flag TCSO. If set, thicknesses for contact nodes are computed based only on the attached segments that are actually part of the contact, not ALL attached segments. It only considers shell elements, and only works in MPP so far (but SMP implementation should be simple) TCSO = Thickness from Contact Shells Only.
• Fix bug in MPP adaptivity: in some cases, the code was hanging if soft=4 contact was used (or groupable contact).
• Change normal calculation for type 7 offset penalty contact to use the simpler “flat segment” normal, instead of the actual normal of the bilinear element. This change was made the other day for the type 2/6 tied contacts since it seems to help stability.
• Steady Static Rolling: modify contact subroutine, improve the calculation of contact node velocity such that SSR can get stable result quickly.
• A bug fix for smooth contact.
• Fixed an MPP bug in neighbour segment checking option of segment based contact that caused jobs to hang if any process did not participate in an interface that uses the option.
• Fixed *CONTACT_2D_AUTOMATIC_SURFACE_IN_CONTINUUM.
• Added new MPP contact option for _TIED_ and _AUTOMATIC_TIEBREAK contacts to compute displacements incrementally — this should take care of lots of issues with noise in non-moving parts.
• Visualization of delamination via “contact gap” in INTFOR database (*CONTACT_AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE_TIEBREAK): Instead of only showing a tied/released flag, now a damage value ranging from 0 to 1 is shown. This is now working for SMP and tiebreak options.
• Added error termination if null set is defined in contact.
• Fix implicit for case where the ignore option is on in contact. Do not add stiffness as the bodies separate.
• Improvements to handling of initial penetration tracking in MPP contacts.
• Add support for viscous damping and soft=1 for MPP penalty based tied contacts of type 2 and 6.
• MPP contact: added support for soft=1 to the penalty based tied contacts (types 2 and 6 only). Also added support for viscous damping to these contacts, as well as the automatic tiebreak contact types.
• Change to 3d thermal contact. The use of the parameter characteristic length (chl, chlmul, char_len) is superseded by the 3d bucket sort algorithm.
• Don’t do *INTERFACE_LINKING / tied contact conflict check for penalty based tied contacts.
• Changed the output of interface force file, used contact id’s for part number instead of the sequential numbers. Assigned positive contact id for the slave side and negative contact id for the master side.
• Fix MPP soft=4 contact: nodal constraints were not being honoured.
• Fixed “input order” dependency from bucketsort results, and also return a better collection of possible contacts, which should reduce spurious force spikes right after bucketsorts.
• Fix in PARMAX for SOFT=4 MPP.
• Fix problem of MPP groupable contact hanging in cases of mixed penalty and constraint groupable contacts
• Switch MPP tied contacts all to type GROUPABLE if the problem (which was not previously the case). Also, make sure some important tied contact related variables get reinitialized each adaptive step.
• Fix for MPP groupable contact (including soft=4 contact). It was not properly supporting frictional energy output to the interface force file (s=intfor).

Control Options

• *CONTROL_IMPLICIT_EXPLICIT_HYBRID
  To run an explicit simulation and specifying a (small) set of elements that are treated by implicit time integration in order to increase the critical timestep.
• *CONTROL_IMPLICIT_SOLUTION
  To invoke modified Crisfield arclength method, supports both load and prescribed motion, and can be activated at specified time.
• *CONTROL_TIMESTEP
  Selective mass scaling supports rigid walls, and selective mass scaling treatment can now include rotational motion for improved accuracy.

• Fix a bug when a box for activation of SPH elements defined on the *CONTROL_SPH is used and when the consistency flag is activated.
• Suppress SPC printing with NPOPT=1 like nodes as suggested by Len Schwer.
• Fixed tolerance bug associated with cntc0 in *CONTROL_SHELL.
• Added laminated shell theory option for thick shell type 5 for all available materials. It is turned on by setting LAMSHT=4 or LAMSHT=5 on *CONTROL_SHELL.
• Add a new keyword *CONTROL_FORMING_PARAMETER_READ. The purpose is to read some parameter values from other files, which might be the output of some other previous runs.
• Added laminated shell theory for most shell materials by assumed strain. It is made active by setting LAMSHT=3 or LAMSHT=5 on *CONTROL_SHELL. The elements that support it are:

  type 1 (Hughes-Liu), type 2 (Belytschko-Tsay),   type 4 (C0 triangle), type 6 (S/R Hughes-Liu), type 7 (S/R co-rotational Hughes-Liu),   type 8 (Belytschko-Leviathan), type 10 (Belytschko-Wong-Chiang),   type 11 (fast co-rotational Hughes-Liu), type 16 (fully integrated),   type 25 (Belyschko-Tsay shell with thickness stretch),   type 26 (fully integrated shell with thickness stretch), and   type 27 (triangular shell with thickness stretch).

• Fixed the C0 triangular shell when used for implicit analysis when the invariant node numbering option is switched off (INN=0 or 1 on *CONTROL_ACCURACY).
• Extend *CONTROL_IMPLICIT_AUTO to allow specification of explicit time interval for implicit-explicit switching to be in terms of the number of explicit cycles.
• Fix for mass scaling bug related to MS1ST=1. Occurs when dynamic relaxation is active.
• Fixes related to selective mass scaling (SMS):
  • introducing flag for accounting rotations in SMS scheme
  • elements that are in a nodal rigid body are not mass scaled
  • accounting for ESORT flag when mass scaling pentahedra in SMS

• Fix value of RBC in explicit to implicit switching handoff. Bug: Discrepancies in output motion associated with implicit/explicit switching.
• Preventing mass scaling of nodal rigid bodies if selective mass scaling is turned on.
• Fix hand off of dynamic relaxation implicit to implicit. D3plot were not written when IDRFLG=5 and DRTERM=ENDTIM.
• Fix of bug related to shell 26 and ESORT.
• Add SPOTHIN option and symmetric behavior to *CONSTRAINED_SPR3 model.
• Default ATYPE on *BOUNDARY_PORE_FLUID should be the user-input ATYPE on *CONTROL_PORE_FLUID. Ensure that this is always so.
• Support in SMP the option to output contact forces to RCFORC for all 2 surface force transducers when slave and master surfaces of the force transducers overlap. The flag is FTALL on the 5th field of the 6th line of *CONTROL_CONTACT.
• Fix for *CONTROL_CHECK_SHELL: remove unnecessary constraints after the node has been deleted.
• Fix for bug in *CONTROL_SOLUTION LCINT=1000 in restart.
• Fix writing of inertia matrix in binary form for  *CONTROL_IMPLICIT_MODES.
• Implicit now honors the ncycle limit from either the command line or *CONTROL_TERMINATION.
• Add *INCLUDE_TRANSFORM compatibility in keyword reader for *CONTROL_PORE_FLUID.
• Add new keyword: *CONTROL_FORMING_NORMALS.
• A bug fix for *CONTROL_FORMING_NORMALS.
• Add new option in *CONTROL_SPH to allow LS-PREPOST to treat deactivate SPH as deleted element
• Support dynamic implicit for *CONTROL_IMPLICIT_FORMING.
• Added *CONTROL_FORMING_BLANKMESH: when cid <0, the center coordinate means local coordinate.

Damping Options

• Fixed energy calculation bug in *DAMPING_PART_STIFFNESS for 2-D.
• *DAMPING_PART_STIFFNESS for the Belyschko-Schwer beam no longer limits the value of the damping coefficient. Also, the energy calculation due to damping is now fixed.

Database/Output

• *DATABASE_ALE_MAT
  This keyword allows for each material group to ouput the volume, mass, internal energy and kinetic energy.
• *DATABASE_TRACER
  Implemented in MPP ALE 2D

• Fix related to bug related to database history with *DATABASE_HISTORY_DISCRETE_SET option.
• Fix sign error in BNDOUT forces for implicit.
• Bug fix for CMPFLG=1 in combination with IHYPER=-1 for user materials.
• Bug fix for CMPFLG=1 and IHYPER=-1.
• Bug fix MPP SPH D3PLOT database SPH element smooth length output.
• Fixed stress output for *MAT_022 when used with 2D continuum elements (shell form 13,14, or 15) and CMPFLG=1 on *DATABASE_EXTENT_BINARY.
• Add contact moments to RCFORC file.
• Add _ID option to *DATABASE_NODAL_FORCE_GROUP. The user defined ID number and title appear in the NODFOR file in the usual “LEGEND” at the top of the file, as well as in a “setid =” field on the ends of the data line for each node. I’m not sure why the “nodal output group number” isn’t a better place to put this ID, but that is how the BNDOUT file was done and I wanted to be compatible.
• Fix output of accelerations of rigid body slave nodes to file NODOUT.
• Shell ELFORM=14 and cross section forces output: the volume weighted forces for ELFORM=14 to be consistent with the resultants from ELFORM=15
• Bug fix for INTFOR file for MPP.
• Fixed a memory error in cross section forces of form 2,3,4,12,13, and 14 fabric elements on systems that use nlq less than 128.
• Fix for bug (set ids wrong in SPCFORC output to LSDA file).
• Fix reporting of moments to SPCFORC for implicit.
• Some fix so that random vibration can go through even if *DATABASE_POWER_SPECTRAL_DENSITY is not present.
• Added LSDA/BINOUT support for SPRFORC data. This is currently only output to the binout file (no ascii file support yet). It is output if and only if there ARE self-piercing rivets in the model AND the SWFORC file is output (this is output at the same frequency).
• Support in SMP the option to output contact forces to RCFORC for all 2 surface force transducers when slave and master surfaces of the force transducers overlap. The flag is FTALL on the 5th field of the 6th line of *CONTROL_CONTACT.
• Fix for bug related to the D3PLOT file not being written.
• Fixed bug for missing binary files (D3PLOT, D3PART,,…) if the file index is greater than 2G.
• Bug fix related to a heading meant for the DEFORC file ending up in the ELOUT file instead.
• Fix for bug related to the output of added mass into the D3PLOT file.
• Standardized SMP/MPP screen output for dynamic relaxation info and also added ‘the dynamic relaxation info header for the D3HSP file.
• Fixed transformation of fabric material during output to D3HSP and ELOUT. It was slightly wrong for large deformation.
• Fix implicit’s control output to SPCFORC.
• Structured write and read for *DATABASE_CURVOUT.
• Fix database corruption for Dynamic Relaxation run for SMP and MPP.
• Some output format changes in file TPRINT.
• Add a flag to the header of BLSTFOR to allow LS-PREPOST to distinguish it from INTFOR. Switch from 21 to 16 variables.
• Visualization of delamination via “contact gap” in INTFOR database (*CONTACT_AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE_TIEBREAK): Instead of only showing a tied/released flag, now a damage value ranging from 0 to 1 is shown. This is now working for SMP and tiebreak options.
• Fix material numbers in output of solid elements to ELOUT file: the internal rather than the user material numbers were being output in ALL versions of this file (LSDA and ascii, SMP and MPP).
• Fix for D3HSP output of LCID from *INITIAL_AXIAL_FORCE_BEAM.
• Minor fix related to Nodout and Node_Merge.
• Corrected thick shell output to the ELOUT file when Gauss integration is not used and there are 3 or 5 integration points through the thickness. The output of stress was in the wrong order when looping over points.
• Add part pressure of alone/away shell normal to abstat_cpm in order to see the pressure on both side from p2f collision.
• Fix for SPHOUT.
• Added new keyword *DATABASE_EXTEND_INTFOR to control the output of INTFOR.
• Changed the output of interface force file, used contact id’s for part number instead of the sequential numbers. Assigned positive contact id for the slave side and negative contact id for the master side.
• Added output to D3HSP for *MAT_088.
• Fixed SWFORC output for spotwelds and also the resultant failure criteria for hex spotwelds.
• *DATABASE_FSI now properly account for wedges and tetras
• Fixed bug in reading *DATABASE_EXTENT_INTFOR.
• Fixed bug: corrupted D3PLOT data from full restart, due to memory issues

Define Options

• *DEFINE_TABLE_2D *DEFINE_TABLE_3D
  Define bi- and tri-dimensional tables.

• Fix for *DEFINE_SEGMENT_GENERAL related to shells and solids sharing identical nodes.
• Fix for *MAT_187. Fix conflicts with curve generation and *DEFINE_CURVE_FUNCTION.
• Fixed bug in shell belt element capability related to curve generation when the *DEFINE_CURVE_FUNCTION options are used.
• MPP fix for parameters used in *DEFINE_CURVE_FUNCTION.
• Changes to *ELEMENT_SEATBELT_RETRACTOR: code added to allow *DEFINE_FUNCION to be used to specify the force vs pull-out curve. This new capability is as yet untested, but many old problems were run to insure these changes did not break existing capability nor change results.
• Fix for bug related to *LOAD_SEGMENT_SET used with *DEFINE_FUNCTION.
• Fix for bug related to load curve offsets.
• Fix for bug related to *DEFINE_TRANSFORMATION combined with *PART_MOVE.
• Fix for bug related to offsets and output into dyna.inc for *DEFINE_CURVE_FUNCTION.
• MPP implementation for *DEFINE_CURVE_FUNCTION , and .
• Implemented *DEFINE_TRACER_PARTICLES_2D.
• Add _SET option for staged construction keywords:

   *LOAD_STIFFEN_PART_SET    *LOAD_REMOVE_PART_SET    *LOAD_GRAVITY_PART_SET    *DEFINE_STAGED_CONSTRUCTION_PART_SET    *INITIAL_STRESS_DEPTH_SET    *INITIAL_PWP_DEPTH_SET

• Bug fix in read/write for mpp, related to *DEFINE_TRACER_PARTICLES_2D.
• Fix for bug: shell element thinning strain from *DEFINE_CURVE_FUNCTION not working correctly in MPP.
• Add missing type declaration related to *DEFINE_FUNCTION_TABULATED.
• A fix for *DEFINE_CURVE_FUNCTION temperature extraction.
• Developed *DEFINE_ADAPTIVE_SOLID_TO_SPH as coupling interface.
• Implementation of *DEFINE_CURVE_FUNCTION TEMP().
• Add support for offsets and scale factors in *DEFINE_CURVE_FUNCTION.

Deformable

• Fix bug about rigid to deformable switching causes I/O problem.
• Fix for *DEFORMABLE_TO_RIGID, code=3 or 4 only compare dist with offset when offset>0; otherwise switch will happen right away since dist=2.*offset and switch condition is automatically met.

Delete Options

• Fix for small deck restart related to bug. Restart issue involving *DELETE_PART.

• MPP support for EFG shell type 41 restart.
• MPP support for EFG solid element restart.
• Bug fix in EFG local adaptivity.
• Fix for a restart in 3D EFG adaptivity with formulation 7.
• Bug fixed in implicit EFG for *MAT_027.
• Bug fixed in EFG *INITIAL_SPRINGBACK_SOLID read/write.
• Added maximum entropy approximation to EFG shell type 42. Fixed options IEBT and IDIM according to the manual.
• Added one option ISPM=2 for EFG *MAT_027 in pressure smoothing.
• Fixed surface meshing problem in adaptive EFG with pressure smoothing formulation.
• Activated *MAT_ADD_EROSION capability to EFG shell type 42.
• Bug fixed in adaptive EFG when there is a formulation switch.
• Added interactive adaptivity into 3D EFG.
• Fixed a bug in strain output for EFG tetrahedra.
• A minor change in EFG crack propagation criterion.
• Bug fixed in adaptive EFG using stabilized method.
• Added one more control parameter in adaptive EFG.
• Bug fix related to the switch between 8-noded and 4-noded cells in adaptive EFG.

Elements

• Fixed an error in the MCID option of *ELEMENT_SHELL that causes incorrect fiber directions. Also, changed behavior such that elements that use the MCID option ignore the material system that is determined by AOPT.
• Added mass scaling to seat belt elements. This scaling was somehow missed in the mass scaling implementation, but is not known to have caused any problems since belt elements generally have a large time step size.
• Changes to *ELEMENT_SEATBELT_RETRACTOR: code added to allow *DEFINE_FUNCION to be used to specify the force vs pull-out curve. This new capability is as yet untested, but many old problems were run to insure these changes did not break existing capability nor change results.
• Fix for bug related to *MAT_103 and *ELEMENT_SHELL, ELFORM=13, NIP=4.
• Added *MAT_145 to 2-D solid elements.
• Fix for bug related to floating point exception in *MAT_015 and shell elements.
• Fix for *ELEMENT_INERTIA to allow a local coordinate system to be defined by using nodal points.
• Bug fix related to *ELEMENT_MASS_PART_SET which caused a memory clobber.
• Fix problem with MPP decomposition related to seatbelt sensors of type 2 and 4, that was causing a memory clobber that led to long initialization times for some models.
• Fix material numbers in output of solid elements to ELOUT file: the internal rather than the user material numbers were being output in ALL versions of this file (LSDA and ascii, SMP and MPP).
• Further accuracy considerations in axisymmetric elements, and also some bugfixes for the implicit implementation of the same element.
• Add solid and SPH coupling.
• Fix for seatbelt mass scaling: belt mass should be scaled based on the ratio of its own dt (dtt) and dt2mas, not global dt (dt2) and dt2mas.

• Allow *MAT_224 to use EOS as an option.

Frequency options

• Implemented *FREQUENCY_DOMAIN_ACOUSTIC_FEM to run frequency domain acoustic computation.
• Enabled free edge constraints for *FREQUENCY_DOMAIN_ACOUSTIC_BEM
• Extended *FREQUENCY_DOMAIN_ACOUSTIC_BEM to study acoustic panel contribution
• Changed keyword *BOUNDARY_ELEMENT_METHOD_ACOUSTIC to *FREQUENCY_DOMAIN_ACOUSTIC_BEM (the old keyword still works).

Hourglass Options

• Add Hourglass option for *MAT_125 with solid element 1.

Implicit

• Implementation of *MAT_091 for tet element 13 with implicit analysis. Slightly modification of implicit *MAT_091 for hex element 2.
• Properly handle dead nodes and tied constraints for implicit.
• Bug fixed in implicit EFG for *MAT_027.
• Fix sign error in BNDOUT forces for implicit.
• Implementation of implicit *MAT_129 for hex and tet13 element Link explicit *MAT_128 to tet13 element.
• Add stiffness for implicit in *CONTACT_AUTOMATIC_GENERAL.
• Add *IMPLICIT_LINEAR_PARTS.
• Fixed the C0 triangular shell when used for implicit analysis when the invariant node numbering option is switched off (INN=0 or 1 on *CONTROL_ACCURACY).
• Extend *CONTROL_IMPLICIT_AUTO to allow specification of explicit time interval for implicit-explicit switching to be in terms of the number of explicit cycles.
• Turn on implicit treatment of joints when *INTERFACE_SPRINGBACK_SEAMLESS is active.
• Fix the interaction of joints and implicit springback.
• Add support for *BOUNDARY_SPC_BIRTH_DEATH to implicit.
• MPP changes to support _AUTOMATIC_TIEBREAK in implicit.
• Fix reporting of moments to SPCFORC for implicit.
• Fix for proper initialization of *CONTACT_AUTOMATIC_TIEBREAK in MPP implicit: nodes were not being tied due to confusion about the definition of “the first cycle”.
• Fix implicit for switching from rigid back to deformable.
• Fix value of RBC in explicit to implicit switching handoff. Bug: Discrepancies in output motion associated with implicit/explicit switching.
• Fix hand off of dynamic relaxation implicit to implicit. D3plot were not written when IDRFLG=5 and DRTERM=ENDTIM.
• A bug fix in combining gravity loading and binder wrapping simulation with implicit method.
• Fixed bug for *MAT_181 tangent for implicit.
• Fix for bug related to erroneous spotweld stresses. Implicit calculations with *MAT_100 solid spotwelds exhibited the error.
• Fix for Implicit MPP cyclic symmetry.
• Changing mpp implementation of implicit tied offset contact to behave similar to smp.
• Fix for implicit spotweld beam material.
• Fix for implicit tied beam offset.
• Fix writing of inertia matrix in binary form for *CONTROL_IMPLICIT_MODES.
• Fixed the material constitutive matrix used in the implicit element stiffness calculation for *MAT_022 and *MAT_059. It is now orthotropic.
• Enable implicit capability for unloading model of *MAT_181 (solids).
• Fix implicit termination when dynamics is active and endtim = 0.0.
• Fix implicit’s control output to SPCFORC.
• Updated checking of thick shell forms for implicit solutions. Now, forms 1 and 4 will switch to 2, and type 5 will switch to 3.
• Support dynamic implicit for *CONTROL_IMPLICIT_FORMING.
• Further accuracy considerations in axisymmetric elements, and also some bugfixes for the implicit implementation of the same element.
• Added material tangents for solids and shells to make *MAT_225 (*MAT_VISCOPLASTIC_MIXED_HARDENING) ready for implicit computations.
• Fix bug in *MAT_003 affecting convergence in implicit computations for ‘shell’-elements using BETA.ne.1.0.
• Fix a bug associated with implicit linearized parts.

Include Options

• Update to include box transformation when *INCLUDE_TRANSFORM is used.
• Fix for units transformation (*INCLUDE_TRANSFORM) for *MAT_138 and *CONTACT_AUTOMATIC_…_TIEBREAK options 9 and 11.
• Fixed error in readin *INCLUDE_NASTRAN.
• A bug fix in *INCLUDE_TRIM.
• Fixed error in *INCLUDE_STAMPED_PART for filename that is longer than 80 characters or file is in *INCLUDE_PATH directory
• Fix for bug when *MAT_ADD_EROSION and *INCLUDE… with INCOUT=1 was used.
• Add *INCLUDE_TRANSFORM compatibility in keyword reader for *CONTROL_PORE_FLUID.
• Fixed: *INCLUDE_TRANSFORM and discrete beam offset.
• Implementation of *INCLUDE_TRANSFORM to *BOUNDARY_PWP_option.
• Fix a bug for *INCLUDE_TRANSFORM for DEFINE_BOX: transformation is needed even when no rotation is involved.
• Added new option to *INCLUDE_STAMPED_PART. New value PERCELE on optional Card 4, Column 3 defines percentage of elements which should be mapped, otherwise error termination occurs.
• Fixed bug – *INCLUDE_TRANSFORM and *AIRBAG_HYBRID_JETTING_CM: Node number NREACT was not offset.

Initial Conditions

• *INITIAL_STRESS_SOLID
  Allow for initialization of specific history variables such as thermal history variables.
• *INITIAL_STRESS_SHELL, *INITIAL_STRESS_SOLID, *INTERFACE_SPRINGBACK_LSDYNA
  Thermal user material history variables are processed in dynain files.
• *INITIAL_ALE_MAPPING
  Restart for the mapping: a d3dump file can be created to restart a job creating a mapping file. The ALE domain of a previous run can be rotated before being mapped to the current mesh. The rotation angle is the 5th parameter on the 2nd line of.  MPP only.

• Bug fixed in EFG *INITIAL_SPRINGBACK_SOLID read/write.
• Bug fix for *INITIAL_ALE_EOS cards.
• Added scale factor to *INITIAL_AXIAL_FORCE_BEAM input.
• Compile a list of nodes in each segment set for *INITIAL_IMPULSE_BLAST.
• Thin shell support for *INITIAL_IMPULSE_BLAST. Westine’s buried mine model.
• New optional field LCK0 on *INITIAL_STRESS_DEPTH allowing input of K0 vs Z-coordinate.
• Add _SET option for staged construction keywords:

   *LOAD_STIFFEN_PART_SET    *LOAD_REMOVE_PART_SET    *LOAD_GRAVITY_PART_SET    *DEFINE_STAGED_CONSTRUCTION_PART_SET    *INITIAL_STRESS_DEPTH_SET    *INITIAL_PWP_DEPTH_SET 

• Fix for D3HSP output of LCID from *INITIAL_AXIAL_FORCE_BEAM.
• Modified *INITIAL_TEMPERATURE behavior: NODES are written before SETS. This enables a methodology to model a sequence of hot stampings where we want to change the blank temperature in each successive stamping by using the keyword *SET_SEGMENT_GENERAL.
• Added internal energy density to *INITIAL_STRESS_SOLIDS.

Interface

• *INITIAL_STRESS_SHELL, *INITIAL_STRESS_SOLID, *INTERFACE_SPRINGBACK_LSDYNA
  Thermal user material history variables are processed in dynain files.
• *INTERFACE_LINKING_NODE, *INTERFACE_LINKING_NODE_SET
  (Just like the _DISCRETE_ variants, but without that term). Does node-to-node matching for linked nodes, matching each to the nearest node in the linking file, independent of ordering. Many customers are misusing *INTERFACE_LINKING_EDGE to get this effect, which generally works but is technically incorrect and might give strange behavior. These users should be able to replace _LINKING_EDGE with LINKING_NODE_SET with no other changes.
• *INTERFACE_SSI_AUX
  To store ground motion data at a free surface in LSDA format, with the intent of using it as free-field ground motion data in a subsequent soil-structure interaction analysis via *INTERFACE_SSI.
• *INTERFACE_SSI_…_OFFSET; _CONSTRAINED_OFFSET

  Allow _OFFSET and _CONSTRAINED_OFFSET options in *INTERFACE_SSI, which defines the soil-structure interface as TIED_SURFACE_TO_SURFACE_OFFSET and TIED_SURFACE_TO_SURFACE_CONSTRAINED_OFFSET contacts respectively.

• Bugfix to *INTERFACE_SSI: resolve confusion about initializing coordinates of soil-side nodes.
• Allow *INTERFACE_SSI to read more than one binary ground motion file for the same interface.
• Turn on implicit treatment of joints when *INTERFACE_SPRINGBACK_SEAMLESS is active.
• Allow _OFFSET and _CONSTRAINED_OFFSET options in *INTERFACE_SSI.
• Fixed a possible segmentation fault in MPP phase 4 when *INTERFACE_LINKING_SEGMENT is used.
• Added new *INTERFACE_LINKING_NODE option that searches for the nearest node to tie to — order doesn’t matter. There was existing checks to make sure master nodes of a tied (constrained) interface were not ALSO slaves due to an *INTERFACE_LINKING_EDGE keyword. This checking is now extended to include nodes slave due to an *INTERFACE_LINKING_NODE keyword.
• Don’t do *INTERFACE_LINKING / tied contact conflict check for penalty based tied contacts.

Keyword

• Fixed a bug in SPF.
• Added *KEYWORD_JOBID This option allows the user to enter a complete jobid after *KEYWORD_JOBDID. The original *KEYWORD_JOD splits the input into 4 sections and combines them with “_”.

Loads

• *LOAD_SEGMENT_GEOMETRY
  Allow the pressure to be defined as a function of the coordinates of the integration points and reference node.
• *LOAD_SEISMIC_SSI
  *BOUNDARY_FREE_FIELD_GROUND_MOTION has been renamed to *LOAD_SEISMIC_SSI.
• *LOAD_BLAST_ENHANCED
  Implemented in ALE 2D. It can be used in 2D with 4-node segments

• New thermal keyword *LOAD_HEAT_GENERATION_SHELL implemented.
• Obsoleting the use of *LOAD_SURFACE_STRESS in conjunction with *CONTACT_…_THERMAL_FRICTION.
• Fix for *LOAD_HEAT_CONTROLLER translation.
• Fix for bug related to *LOAD_SEGMENT_SET used with *DEFINE_FUNCTION.
• Implement *SECTION_BEAM_AISC to incorporate AISC standard steel sections. Rename *BOUNDARY_FREE_FIELD_GROUND_MOTION as *LOAD_SEISMIC_SSI.
• Add _SET option for staged construction keywords:

   *LOAD_STIFFEN_PART_SET    *LOAD_REMOVE_PART_SET    *LOAD_GRAVITY_PART_SET    *DEFINE_STAGED_CONSTRUCTION_PART_SET    *INITIAL_STRESS_DEPTH_SET    *INITIAL_PWP_DEPTH_SET 

• Add an option to disable the negative phase in *LOAD_BLAST_ENHANCED.
• Write non-reflecting boundary scale factor for *LOAD_BLAST_ENHANCED to D3HSP.
• Fix an undocumented bug with *LOAD_BLAST_ENHANCED BLAST=3. The gradient vector used in computing the blast surface normal was not properly transformed to the global system.
• Extend *LOAD_VOLUME_LOSS to solid tet and wedge (ELFORM=10 and 15).
• Add a new option to *LOAD_BODY_VECTOR to allow the load in any director

Material Models

• *MAT_ADD_EROSION
  DIEM, a Damage Initiation and Evolution Model, to arbitrarily combine a damage initiation criterion with an evolution criterion.
• *MAT_ADD_THERMAL_EXPANSION
  To define orthotropic thermal expansion properties for anisotropic materials.
• *MAT_RIGID_DISCRETE
  Rigid material which is discretize into multiple disjoint pieces.
• *MAT_006
  Added temperature dependent relaxation parameters.
• *MAT_034
  To include yarn interaction stiffening through biaxial stress-strain curves, strain rate dependency and hysteresis in cyclic loading.T_041-050 – erosion can be based on number of failed integration points
• *MAT_036
  Hardening tables for strain rate dependent hardening, Chaboche-Roussilier kinematic hardening supported, temperature dependend hardening, can use different hardening and r-values in different directions (0, 45 and 90)
• *MAT_041-050
  Erosion can be based on number of failed integration points.
• *MAT_077
  Frictional (frequency independent) damping can be based on multiple spring-slider systems in parallel.
• *MAT_133
  Hardening tables for strain rate dependent hardening, Chaboche-Roussilier kinematic hardening supported
• *MAT_103
  Can be used with a user defined damage model and hardening can depend on heterogenous multiphase material decomposition
• *MAT_106
  Hardening can depend on plastic strain, plastic strain rate and temperature through *DEFINE_TABLE_3D
• *MAT_135
  Option to use constant strain at necking for transition to power-law hardening, supported for shells 25-27
• *MAT_276
  General viscoelastic Maxwell model, with up to 6 terms in the prony series expansion.
• *MAT_SPOTWELD
  Added a new user-defined weld failure subroutine call uweldfail12. With this new subroutine, it is possible to initiate damage rather than failing the weld immediately. Damage can progress using built in code to do linear or exponential damage, or the damage variable can be updated directly by the user. Added new feature to *MAT_100 beams. If parameter SIGY is given a negative number, then a yield curve/table is assigned by abs(SIGY).
• *MAT_224
  New criterion for Lode angle dependent failure in 3D. Added NUMINT option to material 224 (Card 1, Column 8)
• *MAT_181
  Added unloading options to *MAT_181. Therefore, 3 new parameters are on a new optional card 3/4, columns 1-3. Reloading option for *MAT_034, FORM=14. New parameter RLMAX on Card 6, Column 7, governs reloading behavior. RLMAX=0.0: Reloading on Unloading Curve (default as before) RLMAX>0.0: Reloading on Linear Slope betw. Loading and Unloading Curve
• *MAT_054
  New option for *MAT_054. Number of in-plane integration points for layer failure: 1 (PFL>0) or 4 (PFL<0).
• *MAT_019
  New option RDEF=3 for material *MAT_019. Now we have: RDEF=2: max. principal stress and abs. value of min. principal stress RDEF=3: max. principal stress
• *MAT_183
  Enable *MAT_183 to work with tet type 13.
• *MAT_ADD_EROSION
  Enable *MAT_ADD_EROSION to work with solid type 13. Added new failure criterion to *MAT_ADD_EROSION. If EFFEPS (Crd 1, Clmn 5) has a negative value, then the absolute value is the plastic strain to failure. Added new *MAT_ADD_EROSION damage model called “GISSMO”.
• *MAT_225
  New material model *MAT_VISCOPLASTIC_MIXED_HARDENING or *MAT_225. This is like material *MAT_024, VP=1 but with additional mixed hardening like in *MAT_003.

• Implementation of *MAT_091 for tet element 13 with implicit analysis. Slightly modification of implicit *MAT_091 for hex element 2.
• New material model *MAT_VISCOPLASTIC_MIXED_HARDENING or *MAT_225. This is like material *MAT_024, VP=1 but with additional mixed hardening like in *MAT_003.
• Fix for bug related to *MAT_ARUP and restarts.
• More changes for new *MAT_225.
• Fixed 2D beam element failure so that beams of *MAT_003 and *MAT_015 will fail.
• Fix for bug related to material model *MAT_086 for shell elements.
• Bug fixed in implicit EFG for *MAT_027.
• Fixed thick shell initialization for material types *MAT_54, *MAT_, and *MAT_. Also, added a correction to thick shell formulation 5 to prevent instability due to the assumed z-direction strain.
• Fixed an error in the MCID option of *ELEMENT_SHELL that causes incorrect fiber directions. Also, changed behavior such that elements that use the MCID option ignore the material system that is determined by AOPT.
• Added geometric stiffness for Belytschko-Schwer beams.
• Fixed a type 5 thick shell bug that would cause incorrect stress if the element used mixed materials and some of those materials are isotropic.
• Modified thick shell form 7 for mixed materials. The assumed strain field in the thickness direction now uses the inverse of the 3-3 term of the stiffness matrix to scale the thickness strain such that the strain is correct in each layer.
• Enable *MAT_ADD_EROSION to with thick shell types 3 and 5.
• Fixed form 3 and 5 thick shell elements. For type 3, anisotropic material types *MAT_, *MAT_, *MAT_0, *MAT_0, *MAT_0, *MAT_, *MAT_, *MAT_, *MAT_, *MAT_, and *MAT_ now use the proper material coordinate system for the stress update. For type 5, material types *MAT_, *MAT_, *MAT_0, *MAT_0, *MAT_, *MAT_0, *MAT_, *MAT_, and *MAT_ had unnecessary Jaumann updates on the strains which caused error in the stress to grow in elements that underwent significant rotation.
• Bug fix related to *MAT_153.
• Bug fix for CMPFLG=1 in combination with IHYPER=-1 for user materials.
• Bug fix for CMPFLG=1 and IHYPER=-1.
• Add *MAT_ALE_INCOMPRESSIBLE.
• Implementation of implicit *MAT_129 for hex and tet13 element Link explicit *MAT_128 to tet13 element.
• Updates and fixes for *MAT_ADD_EROSION damage model (GISSMO).
• Further fix for the type 3 solid and *MAT_015.
• Fixed bug where material parameter BETA is read from the wrong card.
• Correction to *MAT_022 and *MAT_059 for SPH elements in order to take care of rotations of SPH elements.
• Implement *MAT_SOIL_CONCRETE (*MAT_078) for SPH elements.
• Fix for bug for 2d solid elements and *MAT_159.
• Fix for *MAT_187. Fix conflicts with curve generation and *DEFINE_CURVE_FUNCTION.
• Fix bug for SPH and *MAT_059.
• Fix problem in output of history variables for SPH when *USER_DEFINED material subroutine is used. Implement *USER_EOS for SPH elements.
• Fix for *MAT_123 internal energy calculation.
• Added NUMINT option to *MAT_224 (Card 1, Column 8).
• Enable solid *MAT_132 to output stress and strain in the material system.
• Fixed failure of materials *MAT_022, *MAT_054, and *MAT_055 with laminated shell theory active. *MAT_022 had an incorrect calculation of the out-of-plane failure function. Also, prior to this fix, all 3 materials had out-of-plane stress that remained non-zero from the time of layer failure until element deletion.
• Enable 2D analysis (shells with ELFORM=13,14,15)  with staged construction and material types *MAT_173, *MAT_192, *MAT_193 and *MAT_198. 
• Fixed shell elements 25 to 27 when used with anisotropic materials when IDOF=1 on *SECTION_SHELL. Materials *MAT_022, *MAT_054, *MAT_055, *MAT_059, *MAT_103, *MAT_104, *MAT_220, and *MAT_221, now transform the strain and stress correctly.
• Fixed stress output for *MAT_022 when used with 2D continuum elements (shell form 13,14, or 15) and CMPFLG=1 on *DATABASE_EXTENT_BINARY.
• Fixed thick shell type 5 and laminated shell theory by assumed strain in thick shells in some cases with mixed materials.
• ALE *BOUNDARY_NON_REFLECTING: Calculate the bulk modulus every time step instead of using the initial value. For HE material the bulk modulus becomes hundreds times smaller after detonated so using the initial value would make it like hitting a still wall and hence kill the timestep size.
• Enabled *MAT_059 with thick shell forms 3 and 5.
• Fixed bug in looking up shear moduli for solid *MAT_040.
• Add stiffness matrix for *MAT_040.
• Fix for units transformation (*INCLUDE_TRANSFORM) for *MAT_138 and *CONTACT_AUTOMATIC_…_TIEBREAK options 9 and 11.
• Fix for failure of truss material *MAT_098.
• Fix for solid material *MAT_019 with VP=1.
• New option RDEF=3 for material *MAT_019.
• New option for *MAT_054: number of in-plane integration points for layer failure: 1 (PFL>0) or 4 (PFL<0).
• Supporting orthotropic expansion properties in *MAT_ADD_THERMAL_EXPANSION.
• Modifications to *MAT_190.
• A bug fix for *MAT_190 when local coordinate system is defined.
• *MAT_135 supported for shells 25-27 (STM and WTM option).
• Reloading option for *MAT_034, FORM=14.
• Add *MAT_226 (Barlat Yield surface + Yoshida kinematic hardening).
• Fixed a memory error in cross section forces of form 2,3,4,12,13, and 14 fabric elements on systems that use nlq less than 128.
• New *MAT_RIGID_DISCRETE.
• Number of minor fixes to *MAT_189.
• Enabled the use of thick shell type 3 and 5 with *PART_COMPOSITE and *MAT_059. It was terminating with an improper error message.
• Fix for *MAT_098 related to Bug report and SIGMAX and SIGSAT.
• Bug fix for joint planes in *MAT_173 and *MAT_198 when tetrahedra are used. The joint plane angles were wrongly calculated for this element type.
• Fix of bug related to *MAT_034 (*MAT_ADD_THERMAL_EXPANSION).
• Further corrections for thermal expansion with *MAT_034 forms 2,3,4 and 12,13,14.
• Option to prescribe plastic strain at necking in *MAT_13 for the s100 option.
• For *MAT_244 calculate and print the martensite, ferrite, pearlite, and bainite start temperatures.
• Fix bug in beam to hex spotweld conversion in *MAT_100.
• Fix *MAT_181 initialization bug for 2-D.
• Added one option ISPM=2 for EFG *MAT_027 in pressure smoothing.
• Added unloading model to *MAT_181 shell materials.
• Fix bug due to roundoff errors with *MAT_100 on some systems.
• Fixed bug for solid *MAT_064.
• Fixed bug on the material erosion of shell elements.
• Enable *MAT_GENERAL_SPRING_DISCRETE_BEAM to read load curves id with more than 7 significant digits.
• Fixed bug in *MAT_119.
• Added a new TSRFAC option for *MAT_034 (*MAT_FABRIC).
• Added damage evolution for *MAT_100 when using OPT=1, stress based failure criteria. FVAL is used as the amount plastic strain the spotweld can take under damage before it fails. If FVAL=0 no damage evolution and the spotweld behaves as before. It fails immediately when the stress exceeds the trigger stress.
• Add option for stress initialization in *MAT_192 (*MAT_SOIL_BRICK) to achieve user-specified horizontal stress (as well as vertical stress), by unloading past the required vertical stress, then reloading. This is enabled by setting THEORY=104.
• Minor fix in curve output format and write the original and fitted curve to file (curveplot_#) for all materials that need curve fitting, such as *MAT_178, *MAT_176, *MAT_077, *MAT_175, *MAT_027 # is the part number.
• Activated *MAT_ADD_EROSION capability to EFG shell type 42.
• Added new feature to *MAT_100 beams. If parameter SIGY is given a negative number, then a yield curve/table is assigned by abs(SIGY).
• Fix for *MAT_ADD_EROSION damage model GISSMO with initial stresses/history.
• Added new film-model *MAT_267.
• Fix for bug related to *MAT_103 and *ELEMENT_SHELL, ELFORM=13, NIP=4.
• Fixed bug for hex spotwelds.
• Bug fix for *MAT_224.
• Fixed bug for *MAT_181 tangent for implicit.
• *PERTURBATION_MAT: change *MAT_238 name from P_PIECEWISE_LINEAR_PLASTICITY to PERT_PIECEWISE_LINEAR_PLASTICITY.
• Added *MAT_145 to 2-D solid elements.
• Fix for *MAT_220.
• Add Hourglass option for *MAT_125 with solid element 1.
• Fix for bug related to floating point exception in *MAT_015 and shell elements.
• Fix for bug related to inconsistency with users manual for *MAT_039. The 1st history variable is now the strain ratio.
• Fix for bug related to erroneous spotweld stresses. Implicit calculations with *MAT_100 solid spotwelds exhibited the error.
• Fix for *MAT_OGDEN, in which the energy calculation was inconsistent with the nodal pressure formulation.
• Fix for *MAT_196 when more than 10000 load curves are present.
• Fixed a bug that prevented the use of user defined weld failure option 22 with hex weld assemblies.
• Fix bug in *MAT_197 (seismic isolator). The response could be noisy when velocity-dependent friction is used in combination with large DISPY. This applies to the friction-type isolators only.
• Fix for bug related to *MAT_015 and the damage calculation.
• Fix for implicit spotweld beam material.
• Bug fix for *MAT_123 with load curve LCTSRF.
• Bug fix for *MAT_224: avoid infinity of triaxiality.
• A bug fix for the shear term of *MAT_226
• Added new failure criterion to *MAT_ADD_EROSION. If EFFEPS (Crd 1, Clmn 5) has a negative value, then the absolute value is the plastic strain to failure.
• Bug fix for *MAT_065. Avoid division by zero.
• Bug fix for *MAT_168. Output real plastic strain.

• Improved convergence of 3D brick elements other than element form 1 by fixing the constitutive matrix used in the element stiffness. Those materials affected are *MAT_011, *MAT_014, *MAT_017, *MAT_019, *MAT_022, *MAT_033, *MAT_035, *MAT_036, *MAT_038, *MAT_051, *MAT_052, *MAT_053, *MAT_059, *MAT_061, *MAT_062, *MAT_064, *MAT_065, *MAT_070, *MAT_072, *MAT_073, *MAT_075, *MAT_079, *MAT_080, *MAT_083, *MAT_084, *MAT_087, *MAT_088, *MAT_089, *MAT_098, *MAT_102, *MAT_105, *MAT_109, *MAT_110, *MAT_111, *MAT_112, *MAT_124, *MAT_126, *MAT_127, *MAT_141, *MAT_142, *MAT_143, *MAT_144, *MAT_145, *MAT_164, *MAT_192, *MAT_193, *MAT_221, and *MAT_223.

• Fixed the material constitutive matrix used in the implicit element stiffness calculation for *MAT_022 and *MAT_059. It is now orthotropic.
• Enable *MAT_ADD_EROSION to work with solid type 13.
• Minor fix for *MAT_ADD_EROSION damage model GISSMO.
• Fix for bug when *MAT_ADD_EROSION and *INCLUDE_ with INCOUT=1 was used.
• Support of extended frictional damping in *MAT_077.
• Improvement for *MAT_155 to solve the result discrepancy between single and double precision.
• Fixed stress output for fabric material forms 12, 13, and 14. Only half the correct stress was being output.
• Fixed bugs in *MAT_121 for UNLOAD OPT = 2.
• New feature for *MAT_054 (thin shells): optional minimum stress limits similar to *MAT_058.
• Supporting temperature dependent hardening in *MAT_036.
• Fixed transformation of fabric material during output to D3HSP and ELOUT. It was slightly wrong for large deformation.
• Enable implicit capability for unloading model of *MAT_181 (solids).
• Added failure to *MAT_067.
• Add a new material: *MAT_242: Barlat 2000 + Yoshida nonlinear hardening.
• Add a new option to *MAT_226 to allow the predicted stress be corrected by the test curve.
• Implemented *MAT_ALE_INCOMPRESSIBLE.
• Enable *MAT_183 to work with tet type 13.
• A bug fix for *MAT_242.
• Fix numerical instability of *MAT_072 for SPH.
• Fix for *MAT_067 caused when failure was added.
• Minor fix for shell materials *MAT_123 (VP=1.0) and *MAT_224. Tightened tolerance.
• Bug fix for *MAT_003 (VP=1.0) with BETA.ne.1.
• Fix for *MAT_123, VP=1, shells (shl123svp). Option NUMINT<0 was broken.
• Minor fix for D3HSP output of *MAT_123 NUMINT value.
• Added NUMINT option to *MAT_187.
• Major update of *MAT_267.
• Bug fix for MPP decomposition with SPH elements and support symmetric BC for SPH.
• Fix a bug affected initial variable in viscoplastic materials.
• Updates related to 3D table lookup and the inclusion of 3D tables in *MAT_106.
• Added effective plastic strain rate as hisvar #6 to *MAT_036 and *MAT_133 for postprocessing.
• Fixed bug. Initialization of *MAT_181 and *MAT_ damage curve.
• Make convergence criteria for *MAT_106 the same for shell and solid implementations.
• Update for *MAT_267: Load curve can now be used for hardening.
• Added missing rotation of back stress for *MAT_VISCOPLASTIC_MIXED_HARDENING (*MAT_225).
• Added physical hourglass control to *MAT_267.
• Add a new option to *MAT_037: to predict material failure with non-linear strain path.
• Fix material numbers in output of solid elements to ELOUT file: the internal rather than the user material numbers were being output in ALL versions of this file (LSDA and ascii, SMP and MPP).
• Fixed possible bug in *MAT_267.
• Fixed the weld assembly ID number that is passed in the user defined weld failure routines.
• Fixed SPH nodal stress update for *MAT_ELASTIC_PLASTIC_HYDRO when failed.
• Adding support for shells 25,26,27 in user materials.
• Updates to *MAT_267. added functionality: viscosity+viscoplasticity more stable
• Bug fix for combination of truss elements with *MAT_003 and BETA=1.0.
• For thick shell form 2, enabled also the use of materials *MAT_036, *MAT_065, *MAT_190, and *MAT_233.
• For thick shell forms 1 and 2, enabled type use of materials *MAT_020, *MAT_030, *MAT_034, *MAT_060, *MAT_076, *MAT_089, *MAT_101, *MAT_106, *MAT_107, *MAT_108, *MAT_109, *MAT_113, *MAT_114, *MAT_115, *MAT_120, *MAT_125, *MAT_134, *MAT_136, *MAT_137, *MAT_153, *MAT_157, *MAT_165, *MAT_170, *MAT_172, *MAT_175, *MAT_176, *MAT_182, *MAT_187, *MAT_188, *MAT_194, *MAT_203, *MAT_204, *MAT_224, *MAT_225, *MAT_226, *MAT_232, *MAT_234, *MAT_235, *MAT_236, *MAT_253 and *MAT_276.
• Now a much more stable *MAT_267.
• Added material tangents for solids and shells to make *MAT_225 (*MAT_VISCOPLASTIC_MIXED_HARDENING) ready for implicit computations.
• Fix bug in *MAT_003 affecting convergence in implicit computations for ‘shell’-elements using BETA.ne.1.0.
• Enable *MAT_123 (VP=1.0/’shell’) to deal with logarithmic strain rates.
• Corrected an error in spot weld material input that cause an erroneous missing curve error. Also, some minor changes for new user failure option 12.
• Added user weld failure option 12. It allows the user to pass back values of eps_to_fail and beta such that they can define different damage growth for different elements. It supports linear damage, exponential decay, and sudden failure. If requested, I will added support for directly updating the damage variable. It can be used by beam welds and hex weld assemblies. Currently none of the user failures are supported for individual hex welds.
• Also, correct the output of stress for fabric material forms 12, 13, and 14 when no reference geometry was used.
• Allow *MAT_224 to use EOS as an option.
• Fixed uninitialized variable in *MAT_224.
• Fix for bug, axisymmetric beam (form=8) with *MAT_003 (BETA.eq.0).
• Added a new feature to *MAT_088 improving accuracy of modeling Tantalum.
• Added output to D3HSP for *MAT_088.
• Fixed a problem in beam routine for spotweld failure option 12. The stress calculation was incorrect.
• Fixed history variable in *MAT_088 – Tantalum option.
• Make *MAT_ADD_EROSION damage model GISSMO available for more shell element types and correction for sorted DKT triangles in combination with NUMINT<0 (*MAT_120, *MAT_123, …).
• Allow *MAT_160 to constrain nodes in local coordinate systems defined by *BOUNDARY_SPC_NODE.
• Fix for *MAT_ADD_EROSION damage model GISSMO. Ensure that the initial plastic strain increment is zero.
• A fix for *MAT_016 using ALE to remove NaNs due to roundoff in advection.
• Fix for *MAT_224 with CP=0.0 to avoid division by zero.
• Fix for *MAT_224. It is now possible to give only one of those load curves defining plastic failure strain – as it should be. Before, curves defining plastic failure strain – as it should be. Before, input of LCI was mandatory, otherwise failure never happened.
• Fixed – Error message for *MAT_168 now types correctly.
• Fix for bug related to mass calculation for material model *MAT_220.

• *NODE_MERGE, *NODE_MERGE_SET
  Nodes with the same coordinates (_TOLERANCE) are replaced during the input phase by the first node encountered that share the coordinate.

• Add MPP support for new *NODE_MERGE option.
• Minor fix related to NODOUT and *NODE_MERGE.

Parameter

• *PARAMETER_LOCAL, *PARAMETER_EXPRESSION_LOCAL, *PARAMETER_DUPLICATION
  Historically, all parameters, once defined, exist throughout the duration of the input processing, and duplicating an existing parameter name was not allowed. The _LOCAL modification is intended to be used in *INCLUDE files. A parameter defined at the top of an *INCLUDE file, and using the _LOCAL modifier, will be deleted when that *INCLUDE file is done being processed. By default, it is still not allowed to have multiple *PARAMETERs with the same name. However, the *PARAMETER_DUPLICATION keyword can be used to modify this behavior, allowing a _LOCAL parameter to temporarily mask a previously defined *PARAMETER or *PARAMETER_LOCAL — the previous value will be restored when the *INCLUDE file is finished being read.

• Fixed bug. Integer values defined with *PARAMETER cards having digits greater than 7 are not being read correctly into single precision versions.
• Add new *PARAMETER_LOCAL and *PARAMETER_EXPRESSION_LOCAL, which define parameters effective only in the current include file (and those included by it). The old *PARAMETER statements define values available globally, so there should be no compatibility problems with old inputs.
• Improve handling of duplicate parameter definitions with new parameter handling code, and better honor the setting of *PARAMETER_DUPLICATION.
• Fix for use of & along with *PARAMETER_EXPRESSION_TITLE.

Part Options

• Added *PART_COMPOSITE_TSHELL option.
• Enabled the use of thick shell type 3 and 5 with *PART_COMPOSITE and *MAT_059. It was terminating with an improper error message.
• Fix for *SET_NODE_GENERAL OPTION=PART where the parts are defined by *PART_COMPOSITE.
• Added *PART_ANNEAL.
• Added *PART_DELETE option for SPH restart.
• Fix for bug related to *DEFINE_TRANSFORMATION combined with *PART_MOVE.
• Bug fix for bug related to *PART_COMPOSITE and *SET_PART_LIST_GENERATE.
• Fix report of masses for materials involved in *PART_MASS.

Pertubation Options

• *PERTURBATION
  Spherical coordinate systems are supported. Some components of *MAT_PIECEWISE_LINEAR_PLASTICITY can be perturbed. Coordinate systems can be mixed in a definition: One can create perturbation in one coordinate system and apply it in another.

• *PERTURBATION: perturbations can now be created Cartesian but applied in spherical.
• *PERTURBATION: perturbations can now be created Cartesian but applied in cylindrical.
• *PERTURBATION: fix bug allowing only one function definition for spectral fields.
• *PERTURBATION fix bug for noisy spectral perturbations for large value of CF1.
• *PERTURBATION fix bug for spectral type = 5 (2D isotropic).
• *PERTURBATION_MAT: change *MAT_238 name from P_PIECEWISE_LINEAR_PLASTICITY to PERT_PIECEWISE_LINEAR_PLASTICITY.
• Fix for *PERTUBATION_MAT: change material id to be part id.

Pore Options

• Change in behavior for *PORE_AIR. Calculation terminates when the min termination time in *CONTROL_TERMINATION or in *PORE_AIR is reached. *PORE_WATER behavior is unchanged.

Rigidwall

• *RIGIDWALL_GEOMETRIC_DISPLAY
  Define a rigid body which represents the shape and physical position of the wall. The motion is updated if the _MOTION option is specified.

• Added DISPLAY option to *RIGIDWALL_GEOMETRIC.

Section Options

• *SECTION_SHELL
  Element types 2 and 16 can have nonzero stress through thickness from contact pressure, user defined elements can take arbitrary number of history variables and extra degrees of freedom.
• *SECTION_SOLID
  User defined elements can take arbitrary number of history variables and extra degrees of freedom.
• *SECTION_BEAM_AISC
  Incorporates the AISC standard steel sections database (Shapes Database v13.1.1) into LS-DYNA. The user can specify a section label, and the appropriate *SECTION_BEAM and *INTEGRATION_BEAM cards will be defined automatically.
• *SECTION_TSHELL
  For thick shell forms 1 and 2, enabled tye use of materials *MAT_020, *MAT_030, *MAT_034, *MAT_060, *MAT_076, *MAT_089, *MAT_101, *MAT_106, *MAT_107, *MAT_108, *MAT_109, *MAT_113, *MAT_114, *MAT_115, *MAT_120, *MAT_125, *MAT_134, *MAT_136, *MAT_137, *MAT_153, *MAT_157, *MAT_165, *MAT_170, *MAT_172, *MAT_175, *MAT_176, *MAT_182, *MAT_187, *MAT_188, *MAT_194, *MAT_203, *MAT_204, *MAT_224, *MAT_225, *MAT_226, *MAT_232, *MAT_234, *MAT_235, *MAT_236, *MAT_253 and *MAT_276. For thick shell form 2, enabled also the use of materials *MAT_036, *MAT_065, *MAT_0190 and *MAT_233. For thick shell form 4, enable all regular shell materials.  Previously, only a subset available to thick shell forms 1 and 2 were available for 4.
• *SECTION_SOLID, ELFORM=17
  ELFORM 17 solid element is deleted under the conditions that ERODE and DTMIN are nonzero, and any of the element’s subelements develops a negative volume.

• Implementation of *MAT_091 for tet element 13 with implicit analysis. Slightly modification of implicit *MAT_091 for hex element 2.
• Fixed 2D beam element failure so that beams of *MAT_003 and *MAT_015 will fail.
• Fix for bug related to material model *MAT_086 for shell elements.
• Fixed thick shell initialization for material types *MAT_054, *MAT_055, and *MAT_059. Also, added a correction to thick shell formulation 5 to prevent instability due to the assumed z-direction strain.
• Enable *MAT_ADD_EROSION to with thick shell types 3 and 5.
• Fixed form 3 and 5 thick shell elements. For type 3, anisotropic material types *MAT_022, *MAT_033, *MAT_054, *MAT_055, *MAT_059, *MAT_103, *MAT_104, *MAT_142, *MAT_143, *MAT_220, and *MAT_221 now use the proper material coordinate system for the stress update. For type 5, material types *MAT_006, *MAT_051, *MAT_052, *MAT_061, *MAT_062, *MAT_079, *MAT_141, *MAT_151, and *MAT_225 had unnecessary Jaumann updates on the strains which caused error in the stress to grow in elements that underwent significant rotation.
• Added maximum entropy approximation to EFG shell type 42. Fixed options IEBT and IDIM according to the manual.
• Implementation of implicit *MAT_129 for hex and tet13 element Link explicit *MAT_128 to tet13 element.
• Fix for bug related to integration point temperatures in the fully integrated solid element.
• Fix for internal energy calculation for type 3 solid element.
• Further fix for the type 3 solid and *MAT_015.
• Fix for solid element type 3 related to energy calculation and equations of state.
• Implement the LEAF option for SPH elements in 2D.
• Implement LEAF for SPH 3D elements.
• Correction to *MAT_022 and *MAT_059 for SPH elements in order to take care of rotations of SPH elements.
• Implement *MAT_SOIL_CONCRETE (*MAT_078) for SPH elements.
• Fix for bug for 2d solid elements and *MAT_159.
• Add initial stress capability for cohesive elements.
• Added laminated shell theory option for thick shell type 5 for all available materials. It is turned on by setting LAMSHT=4 or LAMSHT=5 on *CONTROL_SHELL.
• Fix for bug with Belytschko-Schwer beam. The third node used for orientation had a bad coordinate in the last state of D3PLOT file.
• Enable 2D analysis (shells with ELFORM=13,14,15)  with staged construction and material types *MAT_173, *MAT_192, *MAT_193 and *MAT_198. 
• Fixed shell elements 25 to 27 when used with anisotropic materials when IDOF=1 on *SECTION_SHELL. Materials *MAT_022, *MAT_054, *MAT_055, *MAT_059, *MAT_103, *MAT_104, *MAT_220, and *MAT_221, now transform the strain and stress correctly.
• Fixed stress output for *MAT_022 when used with 2D continuum elements (shell form 13,14, or 15) and CMPFLG=1 on *DATABASE_EXTENT_BINARY.
• Fixed thick shell type 5 and laminated shell theory by assumed strain in thick shells in some cases with mixed materials.
• Enabled *MAT_059 with thick shell forms 3 and 5.
• Added laminated shell theory for most shell materials by assumed strain. It is made active by setting LAMSHT=3 or LAMSHT=5 on *CONTROL_SHELL. The elements that support it are:

  type 1 (Hughes-Liu), type 2 (Belytschko-Tsay),   type 4 (C0 triangle), type 6 (S/R Hughes-Liu), type 7 (S/R co-rotational Hughes-Liu),   type 8 (Belytschko-Leviathan), type 10 (Belytschko-Wong-Chiang),   type 11 (fast co-rotational Hughes-Liu), type 16 (fully integrated),   type 25 (Belyschko-Tsay shell with thickness stretch),   type 26 (fully integrated shell with thickness stretch), and   type 27 (triangular shell with thickness stretch). 

• Fixed the calculation of number of thru-thickness integration points triangle elements when either type 9 membrane elements are sorted to C0 triangles, or when type 18, 20, or 21 elements are sorted to DKT triangles.
• Corrected the material directions for anisotropic materials when used with shell types 12 and 51.
• Shell ELFORM=14 and cross section forces output: the volume weighted forces for ELFORM=14 to be consistent with the resultants from ELFORM=15
• Enabled the use of thick shell type 3 and 5 with *PART_COMPOSITE and *MAT_059. It was terminating with an improper error message.
• Bug fix for joint planes in *MAT_173 and *MAT_198 when tetrahedra are used. The joint plane angles were wrongly calculated for this element type.
• Fix bug in stress initialization for type 16 and 17 solids.
• Fix for solid element elform3 & elform4 restart related to bug.
• Fix for tet type 4 strain tensor is wrong when NCPU=-2,-3, etc.
• Fixed bug for solid *MAT_064.
• Fix of bug related to shell 26 and ESORT.
• Activated *MAT_ADD_EROSION capability to EFG shell type 42.
• Fix for bug related to *MAT_103 and *ELEMENT_SHELL, ELFORM=13, NIP=4.
• Add Hourglass option for *MAT_125 with solid element 1.
• Bug fix for bug related to tets type 17 and eigenvalue analysis.
• Improved convergence of 3D brick elements other than element form 1 by fixing the constitutive matrix used in the element stiffness. Those materials affected are *MAT_011, *MAT_014, *MAT_017, *MAT_019, *MAT_022, *MAT_033, *MAT_035, *MAT_036, *MAT_038, *MAT_051, *MAT_052, *MAT_053, *MAT_059, *MAT_061, *MAT_062, *MAT_064, *MAT_065, *MAT_070, *MAT_072, *MAT_073, *MAT_075, *MAT_079, *MAT_080, *MAT_083, *MAT_084, *MAT_087, *MAT_088, *MAT_089, *MAT_098, *MAT_102, *MAT_105, *MAT_109, *MAT_110, *MAT_111, *MAT_112, *MAT_124, *MAT_126, *MAT_127, *MAT_141, *MAT_142, *MAT_143, *MAT_144, *MAT_145, *MAT_164, *MAT_192, *MAT_193, *MAT_221 and *MAT_223.
• Enable *MAT_ADD_EROSION to work with solid type 13.
• Implement *SECTION_BEAM_AISC to incorporate AISC standard steel sections. Rename *BOUNDARY_FREE_FIELD_GROUND_MOTION as *LOAD_SEISMIC_SSI.
• *DAMPING_PART_STIFFNESS for the Belyschko-Schwer beam no longer limits the value of the damping coefficient. Also, the energy calculation due to damping is now fixed.
• Enable *MAT_183 to work with tet type 13.
• Updated checking of thick shell forms for implicit solutions. Now, forms 1 and 4 will switch to 2, and type 5 will switch to 3.
• Now ELFORM 17 solid element is deleted under the conditions that ERODE and DTMIN are nonzero, and any of the element’s subelements develops a negative volume.
• Adding support for shells 25,26,27 in user materials.
• Bug fix for combination of truss elements with *MAT_003 and BETA=1.0.
• For thick shell form 4, enable all regular shell materials.  Previously, only a subset available to thick shell forms 1 and 2 were available for 4.
• For thick shell form 2, enabled also the use of materials *MAT_036, *MAT_065, *MAT_190 and *MAT_233.
• For thick shell forms 1 and 2, enabled type use of materials *MAT_020, *MAT_030, *MAT_034, *MAT_060, *MAT_076, *MAT_089, *MAT_101, *MAT_106, *MAT_107, *MAT_108, *MAT_109, *MAT_113, *MAT_114, *MAT_115, *MAT_120, *MAT_125, *MAT_134, *MAT_136, *MAT_137, *MAT_153, *MAT_157, *MAT_165, *MAT_170, *MAT_172, *MAT_175, *MAT_176, *MAT_182, *MAT_187, *MAT_188, *MAT_194, *MAT_203, *MAT_204, *MAT_224, *MAT_225, *MAT_226, *MAT_232, *MAT_234, *MAT_235, *MAT_236, *MAT_253 and *MAT_276.
• Minor bugfixes in *SECTION_BEAM_AISC.
• Added parabolic shear to the thick shell formulation 5 when laminated shell theory is not used. Previously, it was constant through the thickness so it is now more accurate.
• Extend *LOAD_VOLUME_LOSS to solid tet and wedge (ELFORM=10 and 15).
• Fix when unreferenced *SECTION_SOLID with elform 11 appears in the keyword file for a non-ale run.
• Fix for bug, axisymmetric beam (form=8) with *MAT_003 (BETA.eq.0).
• Fixed brick spotwelds and assemblies with failure type other than option 8 and 12 when swforc output is not requested. A failure to allocate sufficient memory caused spurious stress at the start of the calculation.
• Improve accuracy of tetrahedron type 13 by making the averaged nodal volumes (current and initial) always double precision.

Set Options

• *SET_…_COLLECT
  All set components with same SID will be collected into one big SET with the same setid. Applies to

    *SET_PART   *SET_SEGMENT   *SET_NODE   *SET_SOLID   *SET_SHELL   *SET_TSHELL   *SET_BEAM

• Fixed problem in *SET_SEGMENT_GENERAL where first part ID is 0 causing all parts to be included in the contact.
• Fix related to bug where the behavior of option BOX_… in *SET_SEGMENT_GENERAL is completely contrary to the users manual. 
• Fix for *SET_SEGMENT generation with the _GENERATE option.
• Problem with *SET_PART_LIST_GENERATE is now fixed.
• Fixed bug for *SET_PART_LIST_GENERATE if total number of sets exceeds 500.
• Fix for *SET_NODE_GENERAL OPTION=PART where the parts are defined by *PART_COMPOSITE.
• Bug fix for bug related to *PART_COMPOSITE and *SET_PART_LIST_GENERATE.
• Fix a bug for *SET_SEGMENT_GENERAL.
• Modified *INITIAL_TEMPERATURE behavior: NODES are written before SETS. This enables a methodology to model a sequence of hot stampings where we want to change the blank temperature in each successive stamping by using the keyword *SET_SEGMENT_GENERAL.

• Bug fix MPP SPH D3PLOT database SPH element smooth length output.
• Fix a bug when a box for activation of SPH elements defined on the *CONTROL_SPH is used and when the consistency flag is activated.
• Fix bug for SPH and *MAT_059.
• Fix problem in output of history variables for SPH when *USER_DEFINED material subroutine is used. Implement *USER_EOS for SPH elements.
• Fix for bug related to SPH mass and inertia missing from the body output.
• Fix numerical instability of *MAT_072 for SPH.
• Add new option in *CONTROL_SPH to allow LS-PREPOST to treat deactivate SPH as deleted element
• Add support of DTMIN and erode options for SPH
• Added *PART_DELETE option for SPH restart.
• Support SPH delete parts for MPP.
• Fix element order dependency in SPH: elements are sorted in user ID and then part ID to ensure the order same to get consistent answer.
• Fixed SPH nodal stress update for *MAT_ELASTIC_PLASTIC_HYDRO when failed.
• Bug fix for SPH readin memory allocation.
• Add solid and SPH coupling.
• Improved speed and accuracy of axisymmetric SPH (SMP only).

Termination

• Fix for *TERMINATION_CONTACT which did not seem to be implemented completely as described in the manual.

• Fix problem in output of history variables for SPH when *USER_DEFINED material subroutine is used. Implement *USER_EOS for SPH elements.

Restart_Input_Data

• Fix bug (MPP restart with *CONSTRAINED_TIED_NODES_FAILURE).
• MPP support for EFG shell type 41 restart.
• MPP support for EFG solid element restart.
• Fix for a restart in 3D EFG adaptivity with formulation 7.
• Fix for bug related to *MAT_ARUP and restarts.
• Fix bug (MPP restart with *CONSTRAINED_TIED_NODES_FAILURE).
• To fix the MCOL restart for both SMP and MPP.
• Fix bug: incorrect behavior of tied-offset contact under full-deck restart. After the restart, the contact was being (re-)initialized using the coordinates at the end of the initial run, which caused the nodes to be attached to different segments than they were to in the initial run, thus  changing the contact stiffnesses.
• Fix for bug in *CONTROL_SOLUTION LCINT=1000 in restart.
• Fixed bug: corrupted D3PLOT data from full restart, due to memory issues
• Fix for small deck restart related to bug. Restart issue involving *DELETE_PART.
• Added *PART_DELETE option for SPH restart.
• Fix for solid element elform3 & elform4 restart related to bug.
• Added a restart option to add contribution from higher frequencies in FRF computation.

General

• Fixed bug related to length of jobid greater than 72.
• Fix for seatbelt mass scaling: belt mass should be scaled based on the ratio of its own dt (dtt) and dt2mas, not global dt (dt2) and dt2mas.
• Fix some logic errors that were causing MPP element costs to be incorrect for some elements, which could affect the decomposition and possibly the load balance.
• Fix some cases where MPP files were not properly copied back from the working directories if jobid= was used.
• For indirect coupling with MADYMO 7.2
  • Fixed exit problem
  • Sync d3plot
  • Added restart

• Bug fix: the code was not working if “i=” was NOT on the command line