lsdyna R10.2 リリースノート « ACS NEWS

lsdyna R10.2 リリースノート

2019 年 4 月 1 日 by yo yamgatta


This file constitutes revision 0 of the release notes for LS-DYNA version R10.2.0.


The string "REVISION 10" must appear in the LS-DYNA license file in order to run
version R10.2.0. Please contact your LS-DYNA distributor ( if
you're a direct customer of LSTC) if you need to have your license updated.



Documentation of R10.2.0 is provided in the R10.0 User's Manuals which may
be downloaded from
For features mentioned in these release notes but which appear to be missing
from the R10.0 User's Manuals, please refer to the DRAFT User's Manuals at Please note though that not all the features
in the DRAFT User's Manuals are available in version R10.2.0.



The remainder of this file describes what's been added or fixed in version R10.2.0
since the release of version R10.1.0. The changes are primarily just bug fixes.

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.

The categories are:

________________________ *AIRBAG ___________________
________________________ ALE ___________________
________________________ *BOUNDARY ___________________
________________________ CESE (Compressible Fluid Solver) ________________
________________________ CONTACT ___________________
________________________ *CONSTRAINED ___________________
________________________ *CONTROL ___________________
________________________ Discrete Element Method ___________________
________________________ EFG (Element Free Galerkin) _____________
________________________ *ELEMENT ___________________
________________________ Forming Analysis __________
________________________ *FREQUENCY_DOMAIN ___________________
________________________ ICFD (Incompressible Fluid Solver) ______________
________________________ Implicit (Mechanical) Solver _______
________________________ *INITIAL ___________________
________________________ Isogeometric Elements ___________________
________________________ *LOAD ___________________
________________________ *MAT and *EOS ______________
________________________ MPP ___________________
________________________ Output ___________________
________________________ Restarts ___________________
________________________ *SENSOR ___________________
________________________ SPG (Smooth Particle Galerkin) __________________
________________________ SPH (Smooth Particle Hydrodynamics) ______________
________________________ Thermal Solver ____________
________________________ Miscellaneous ___________________


________________________ *AIRBAG ___________________

Exclude *AIRBAG_LOAD_CURVE and *AIRBAG_LINEAR_FLUID airbags from the
instability stablilization added to *AIRBAG_INTERACTION in
r117308/R9 and 117309/R10.

- 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

________________________ ALE ___________________

Fix a bug in the routine reading the binary file that created empty files with
random numbers appended to their names.

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

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.

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

________________________ *BOUNDARY ___________________

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
combined with translational motion.

Fix processing of prescribed motion constraints for explict computation with superelements

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.

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:
a) Any used segments that are not in the CESE mesh are reported.
b) Element sets look up is now guaranteed to work.
c) Drag calculation output that is attempted for the incorrect FSI method,
is now reported. Here are the illegal cases:
i) *CESE_DATABASE_FSIDRAG is only available with the CESE immersed
boundary FSI solvers.
ii) *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.
d) Failure to find an internal node number in *CESE_DATABASE commands is now
e) Output was not being performed for the initial conditions. Instead, it
is now done after the first time step.

________________________ CONTACT ___________________

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

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

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.

________________________ *CONSTRAINED ___________________

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.


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

________________________ *CONTROL ___________________

- 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.

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

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

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.

________________________ *ELEMENT ___________________

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.

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.

Fix bug of not filling the corner with elements properly.

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

- 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.

Fix duplicate node set problem, which caused the simulation to crash.

________________________ *FREQUENCY_DOMAIN ___________________

Fixed a bug in computing RMS (Root Mean Square) values of Von Mises stress for
solid elements.

Fixed nodal displacement output to d3spcm. Only the nodal displacement (u) itself
(not updated nodal coordinates, which is u+x0) is written to d3spcm.

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

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

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.

________________________ *INITIAL ___________________

Fix seg fault or incorrect stresses when initializing stresses using *INITIAL_STRESS_SOLID

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 ___________________

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

________________________ *LOAD ___________________

Fix output to nodfor when using *LOAD_SSA.

Fix ineffective *LOAD_THERMAL_CONSTANT when used on *MAT_106_THERMAL/


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

i.e. using user-defined failure subroutine, matusr_24 in dyn21.F, to determine

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

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/

Fix contact penetration after element erosion when using *MAT_110/

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/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

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

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

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

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.

- 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.

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

________________________ MPP ___________________

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.

________________________ Output ___________________

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
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.


Fixed bug for d3plot if DCOMP=5 or 6 in *DATABASE_EXTENT_BINARY and

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.

________________________ Restarts ___________________

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.

Fix bug in accessing ground motion ID (*DEFINE_GROUND_MOTION) from

________________________ *SENSOR ___________________

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

________________________ 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

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

Improved reporting of the thermal time step.

________________________ Miscellaneous ___________________

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.

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:
Fix bug which resulted in a corrupt node list.

====================== END of Release Notes for version R10.2.0 =============================