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

LS-DYNA R10.1.0 (R10.123355; 2018/02) released

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

Herein are summarized new features and enhancements in version 10.1.0 since the release of version R10.0.0.

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

Airbag
  • Fix an airbag bug that can cause MPP to crash when an airbag reference geoemtry is used in a huge model.
  • Improvements to *AIRBAG_PARTICLE:
    • Extend the exterior air drag force capability (CD_EXT) to UP airbag after switching CPM (corpuscular particle method) to UP (uniform pressure) via TSW variable.
    • FRIC, PFRIC parameter for surface roughness:
LT.-1                 Absolute value is ID of curve defining FRIC vs. part pressure.
GT.-1 .AND. .LT.1     Specifies the factor.
GT.1                  ID of *DEFINE_FUNCTION defining FRIC.  Bag pressure, bag volume, and current time are passed as arguments into the function.
    • NP_RELX:
      • This value will be automatically set to -abs(NP_RELX) for stability while using IAIR=4.
  • *MAT_ADD_AIRBAG_POROSITY_LEAKAGE now works with *AIRBAG_PARTICLE but only for VOPT=7/8. Data on this card has higher priority than data specified in *MAT_FABRIC.
    • FAC:
LT.0                   ID of *DEFINE_FUNCTION. Absolute part pressure and current time are passed as arguments into this function that determines the leakage factor.
              -or-
                       ID of load curve for leakage velocity vs. absolute pressure.
GT.0                   Leakage factor (a constant).
  • *DEFINE_CPM_BAG_INTERACTION treats airbags with chambers and allows interaction with different gas components in each bag.

.

ALE
  • Fixed bug in Structured ALE for cases with ESORT=1 in *CONTROL_SOLID. Element ID mismatch would cause ALE element volume fraction to be assigned wrong value.

.

BOUNDARY
  • *BOUNDARY_MCOL: Correct printing to files "mcol" and "mcolkine" in both SMP and MPP. Solution is unaffected by this correction.
  • *BOUNDARY_SPC_SYMMETRIC_SET: Fix a bug related to option _SET, so the option works with solid elements.

.

Blast
  • Fixed *PARTICLE_BLAST bug if detonation point is defined using a node (MPP only).

.

CESE (Compressible Fluid Solver)
  • Corrected a prescribed boundary input reading error. This fix will only affect calculations if the user accidentally used the 'SEGMENT' option in the '*CESE_BOUNDARY_PRESCRIBED_OPTION' card.
  • Corrected geometry deduction errors in the 2D rotation case of CESE cyclic boundary conditions.
  • Internal sorted element indexing is now taken into account for the structural elements involved in the newer CESE conjugate heat transfer solvers.
  • CESE immersed-boundary method (IBM) FSI solvers: For wedge (prism) elements, fixed some internal geometry calculations for the 3D immersed boundary conjugate heat transfer solver.

.

CONTACT
  • Fix reported forces and energies for *CONTACT_TIED_SHELL_EDGE_TO_SOLID in MPP.
  • Add command line option "soft=" for switching value of SOFT in contacts:
    • "soft=2to1" converts all SOFT=2 contacts to SOFT=1;
    • "soft=1to2" converts SURFACE_TO_SURFACE and SINGLE_SURFACE contacts with SOFT=1 to SOFT=2. This option does not set other variables related to SOFT=2, e.g., SBOPT and DEPTH.
  • Fix *DEFORMABLE_TO_RIGID switching using contact forces, CODE=2 & 4, when the contact type is *CONTACT_FORCE_TRANSDUCER.
  • Fix SMP consistency mode (ncpu<0) for *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_COMPOSITE.
  • Fixed ISLIDE option of *CONTACT_2D_AUTOMATIC in MPP version. Some data was incorrectly accessed leading to possible termination or bad results.
  • Fixed an error in segment based (SOFT=2) eroding contact that caused the contact to not create new segments on interior elements when there were no active segments at the start of the run. This error was most likely to occur in MPP runs when decomposition can result in a partition that has no active segments at the start of the run.
  • Changed the behavior of segment based contact when the contact keyword indicates a surface to surface contact, but the master surface has no segments. In older versions, the contact would use the slave segments and do a single surface contact with them. Going forward, the contact will not do anything because the slave surface has nothing to contact.
  • Fixed some contact options when using linear 3D solids or linear shell elements in implicit analysis. These elements were causing a zero explicit time step, and some contact algorithms use this explicit step in the denominator of the stiffness matrix. These include SOFT=1 contact, SOFT=2 contact, and tied contact by the penalty method.
  • Fixed SMP segment based contact with CNTCO.ne.0 on *CONTROL_SHELL. Nodal coordinates were being changed.
  • Fixed a flaw in the support of SSF on *PART_CONTACT when using segment based (SOFT=2) contact with DEPTH=5. The consequence was to sometimes choose the wrong part when looking up the scale factor so the stiffness was not predictable.
  • Fix MPP problem where initial coordinates of shell nodes in contact may be moved inappropriately when NLOC and CNTCO are nonzero.
  • Fix output of FORCE_TRANSDUCERs that collect groupable beam-to-beam contact forces.
  • Fix loss of continuity in contact forces for several time steps after a bucket sort occurs in *CONTACT_AUTOMATIC_SINGLE_SURFACE.
  • Fix *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE to detect beam contact from both sides of master segments.

.

CONSTRAINED
  • For the case of *CONSTRAINED_INTERPOLATION where two independent nodes had the same coordinates, a divide by zero could occur during initialization. This was corrected. A further correction was made to the colinearity check of *CONSTRAINED_INTERPOLATION in MPP.
  • Add option to define a part ID for visualization beams used with *CONSTRAINED_SPR2.
  • *CONSTRAINED_SPR2/SPR3: Remove visualization beams from automatic contact. Undesired contact based deformations could occur.
  • Fix for *CONSTRAINED_SPR2 with more than two sheets (XPIDi>0 on Card 4). Errors such as "Beam # ... has zero length" were possible.

.

CONTROL
  • Fixed bugs in Staged Construction (*CONTROL_STAGED_CONSTRUCTION).
  1. Tshells were not writing element history data to dynain file.
  2. When starting the analysis at a stage with non-zero start time (*CONTROL_STAGED_CONSTRUCTION, STGS>1), output data was being sent to ASCII and binout files every time step, no matter what intervals were set on *DATABASE_.... This could have slowed down the code, or made the output files very large. The bug did not affect d3plot or d3thdt files.
  3. When LS-DYNA automatically creates beam elements duplicating *ELEMENT_DISCRETE (for plotting purposes), these beams should not be written to the dynain files.
  4. *LOAD_GRAVITY_PART and *DEFINE_STAGED_CONSTRUCTION_PART were unnecessarily issuing an error termination if parts containing beam elements with ELFORM=3 were referenced.
  5. *DEFINE_STAGED_CONSTRUCTION_PART automatically generates gravity loading by creating *LOAD_GRAVITY_PART, but it should not do this if the user sets ACCEL=0 on *CONTROL_STAGED_CONSTRUCTION. This would have no effect on results (because the applied gravity acceleration is zero) but computational effort was wasted.
  6. Fixed bug affecting large format dynain files (as used in Staged Construction, or if FTYPE=10 on *INTERFACE_SPRINGBACK). The *INITIAL_STRAIN_SHELL data was written incorrectly under these conditions:
    • "Large" dynain file format is used (as is the case in Staged Construction),
    • MPP on multiple processors (the bug does not affect SMP),
    • The model contains shell elements with a mix of numbers of in-plane integration points, for example ELFORM=2 (1 in-plane point) and ELFORM=16 (4 in-plane points).
    If the dynain file was then used as input to a second analysis, LS-DYNA could not read it. The work-around was to delete the *INITIAL_STRAIN_SHELL data from the dynain file: this strain data is for post processing only and does not affect other results.
  • Added option to evaluate *DEFINE_CURVE_FUNCTION every Nth cycle by setting position 6 on *CONTROL_SOLUTION to an integer value N. This setting applies to all *DEFINE_CURVE_FUNCTION keywords in the model.
  • Allow ADPOPT=7 of *CONTROL_ADAPTIVE to be used with *MAT_103. This is 3D r-adaptive remeshing with tetrahedral solid elements type 13.
  • *CONTROL_ADAPTIVITY, *CONTROL_IMPLICIT: A bug that resulted in apparent deformation of rigid bodies in implicit solutions involving adaptivity has now been fixed
  • The IACC flag of *CONTROL_ACCURACY is now echoed to d3hsp.

.

Discrete Element Method (DEM)
  • Fixed *DEFINE_DE_MASSFLOW_PLANE output error if multiple planes are defined.
  • Fixed *DEFINE_DE_INJECTION bug if multiple injection planes are defined.
  • Fixed wear depth calculation error for DEM interface force file (MPP only).
  • Fixed bug in calculation of DEM-DEM capillary force (MPP only).
  • Fixed bug in calculation of *DEFINE_DE_TO_SURFACE_COUPLING capillary force.
  • Fixed DES injection bug such that DEM particles with different radii can be injected.

.

ELEMENT
  • Fixed bug that could occasionally cause LS-DYNA to crash during the input stage when automatically-sorted solid elements are present (wedge or tetrahedron elements with ESORT=1 on *CONTROL_SOLID).
  • Fixed bug affecting *PART_COMPOSITE. If the first part in the model had very low density, or was a part for Discrete Elements that have no density, large rotational inertias could be added to the nodes of the *PART_COMPOSITE. The inertia could artificially resist the correct rotational motion of the elements. In this context, "first part in the model" means the first part encountered by the keyword reader, not the part with the lowest PID. For this reason the bug might be observed or not according to the order in which Include files were defined. For the effect of the bug to be noticeable, the "very low density" would have to be at least one or two orders of magnitude less than the density of the *PART_COMPOSITE itself.
  • Enhance SPC constraint handling for explicit using superelements to properly handle the impicit/explicit switching case.  Also modify explicit handling of superelements to recognize *BOUNDARY_SPC definitions instead of just those on the *NODE cards.
  • Fix storage allocation for *ELEMENT_DIRECT_MATRIX_INPUT when the matrices have different orders.
  • Correct error message of KEY+109 which mis-printed the section number and material number for a 2D seatbelt.
  • Corrected the thinning of shell formulation 24. It was thinning about 33% too much under large tensile strain.
  • Enabled adaptive tet meshing (ADPOPT=7 on *CONTROL_ADAPTIVE) to work with any reasonable solid formulation at the start of the calculation. Previously, the part could only be tet meshed with form 13 elements, or else EFG solid formulations 41 or 42, or the job would crash at the first remap step.
  • Fixed a segmentation fault that occurred in models with beam elements when some elements in a part use *ELEMENT_BEAM_ORIENTATION, and some do not. The error affected beam formulations are 1, 2, 11, 12 and 14.
  • Improved the accuracy of tet formulation 13.  During rigid body translation, some pressure could develop.  Now it can translate pressure-free, and therefore stress-free.
  • Fixed parallel consistency for beam formulation 14, the elbow beam.
  • Fixed dynain writing and reading of *INITIAL_STRESS_SHELL for the fully integrated C0 shell (shell form 20).
  • Fixed adaptivity (and restart) for tet formulation 13 by porting nodal averaged data between adaptive steps. The bug caused spurious peaks in cross section forces, which could be confused with contact pressure peaks.
  • Enable solid element type 23 (20-node hexa) to be used with *MAT_ADD_EROSION.
  • Fix 2D axisymmetric adaptivity failure.
  • Fix bug in adaptivity triggered by presence of *ELEMENT_DISCRETE.
  • Fix calculation of eroded hourglass energy in MPP.

.

Forming Analysis
  • *CONTROL_FORMING_AUTOPOSITION_PARAMETER: Fix for possible divide by zero in case THICK=0.0.
  • *CONTROL_FORMING_OUTPUT: Fix bug preventing d3plot data from being output between TBEG and TEND if LCID and Y1~Y4 are undefined.
  • Fix bug in MPP in which nonnzero stresses were written to *INITIAL_STRESS_SOLID in dynain for rigid solids.
  • *DEFINE_CURVE_TRIM: Bug fix that prevented trimming of solids in MPP.
  • Fix failure of *CONTACT_FORMING_NODES_TO_SURFACE with SOFT=6.
  • *INTERFACE_SPRINGBACK_COMPENSATION: Bug fix related to the stack size, which caused a segmentation fault to occur.

.

FREQUENCY_DOMAIN
  • *FREQUENCY_DOMAIN_FRF: Fix a bug in writing frf output for rotational dof. This corrects the FRF results when loading is defined as base angular velocity, base angular acceleration or base angular displacement, and absolute response (RELATV=0) is requested.
  • Fixed a bug when the loading and response are both defined in z-direction in frequency response function (FRF) computation.
  • Frequency domain (BEM Acoustics) Fixed a bug in reading *FREQUENCY_DOMAIN_ACOUSTIC_BEM keywords in MPP.

.

ICFD (Incompressible Fluid Solver)
  • Fixed the FSI computation with moving walls which was broken in R10.0.
  • Fixed the recognition of inflows in level set that allows for volumes to be empty when inflows of free-surface are present.
  • Fix the viscosity as a function of temperature when using Non-Newtonian fluids with NNID=6,7,8 in Free-Surface problems.
  • Fixed the assignment of PID for the boundary layer nodes that are inside the volume.
  • Fixed conjugate heat transfer solver with adaptivity.
  • Added support of separate part id drag to steady state solver.
  • Add non inertial reference frame feature to steady state solver.
  • Fixed heatflux computation.

.

Implicit (Mechanical) Solver
  • Redo implicit's treatment of nodal rotational inertias to match other changes in LS-DYNA.
  • Backport the Implicit Linear Multistep feature (NSOLVR=-1) from Dev to R10. Correct Implicit Linear Multi-step to get onlinear elements to work correctly.
  • Include pertinent damping-related data in dynain for *INTERFACE_SPRINGBACK in implicit static analysis. Such data is needed if dynain is subsequently used in a dynamic analysis.
  • Update Implicit's collection of damping terms for discrete elements to account for user specified coordinate system and the case of no second node.
  • Correct typing of disp array in implicit modal dynamics. It was not typed so it was real*4 instead of real*8 in single precision. This led to incorrect output to nodout file.
  • Make correction to get proper output to the binary format of superlement file generated by *IMPLICIT_STATIC_CONDENSATION.
  • The separable component report in implicit was extended to include rigid bodies. Also reverted to format that produces a prettier table instead of having the spaces squeezed out by lsmsg.
  • Logic to handle *CONSTRAINED_INTERPOLATION with large number of independent nodes was extended for *CONTROL_IMPLICIT_MODES.
  • *INITIAL_STRESS_SECTION, *CONTROL_IMPLICIT_GENERAL Fix bug in implicit treatment of *INITIAL_STRESS_SECTION. The prestress might have been zeroed out if the implicit solver required a cut-back of the time step.
  • In MPP, the last file of the d3eigv family of files was missing an end-of-file marker but is now correctly marked.
  • Fix bug in writing of d3eigv data for model with *CONTROL_IMPLICIT_INERTIA_RELIEF.

.

*INITIAL
  • Enabled multiple *INITIAL_VELOCITY_GENERATION keywords to be used with *ELEMENT_SHELL_COMPOSITE or *ELEMENT_TSHELL_COMPOSITE. Only one velocity generation keyword was supported previously. Fixed a similar problem for *INITIAL_VELOCITY_GENERATION when using it with both *ELEMENT_SHELL_COMPOSITE and *ELEMENT_TSHELL_COMPOSITE in the same model.
  • Fix for *INITIAL_AXIAL_FORCE_BEAM with KBEND=1. Wrong internal index could have led to unintended memory overwrite, thus arbitrary problems.
  • Format of *INITIAL_CONTACT_WEAR, when written to dynain, was broken in early versions of R10, but is now fixed.
  • *INITIAL_STRAIN_SHELL data written to dynain for rigid shell(s) is now just ignored rather than triggering a fatal error.

.

*LOAD
  • Fixed *LOAD_THERMAL_VARIABLE when used with shell formulation 2 that has default warping stiffness (BWC=2 on *CONTROL_SHELL) and the objective stress updated turned on (OSU=1 on *CONTROL_ACCURACY).
  • Fix bug in *LOAD_HEAT_GENERATION in MPP. Some of the parameters had been ignored.
  • *LOAD_HEAT_CONTROLLER now works in MPP.
  • Fix broken follower force option of *LOAD_NODE in MPP.
  • Fix MPP bug for *LOAD_NODE_POINT follower force with M3<=0.

.

*MAT and *EOS
  • Fixed bug that could potentially affect *MAT_021, *MAT_076, *MAT_172, *MAT_215, *MAT_219, *MAT_249, *MAT_277, *MAT_278, *MAT_293. Where two *MAT cards are almost identical, LS-DYNA could sometimes ignore the differences and treat the elements as if they all were made of one of the two *MAT cards. Whether the bug occurred or not depended on which input fields were different. For example with *MAT_172, the bug occurred if the only differences were on Cards 7, 8 or 9.
  • *MAT_208/*MAT_BOLT_BEAM:  Add erosion and failure due to axial compression.
  • Fix bug in *MAT_172 that could cause a crash on certain operating systems if the input data is pure reinforcement in the y-direction (FRACRX=0 and FRACRY=1), and concrete compressive strength FC=0. Because the fraction of concrete is zero, it should not matter that FC=0. A work-around is to define non-zero FC, which will have no effect on results.
  • *MAT_172: Thick shells (*ELEMENT_TSHELL) were not working with *MAT_172 reinforcement. Even if the input field FRACRX or FRACRY > 0, it behaved as if no reinforcement was present. Now fixed.
  • Fixed bug that could cause a crash when *MAT_173 was used with ALE Multi-Material.
  • *MAT_169/*MAT_ARUP_ADHESIVE was causing LS-DYNA to crash when elements were deleted due to reaching the failure criteria, and the bond thickness card was defined in the input file. This bug affected R10.0 but not R9.
  • Fix *MAT_ADD_THERMAL_EXPANSION which did not work for hyperelastic materials in 2D analysis for fully-integrated 2D elements, e.g., shell type 15 with NIP=4.
  • Fix immediate deletion of elements at time 0 when using *MAT_PIECEWISE_LINEAR_PLASTICITY with FAIL>0.0 and IACC=1 in *CONTROL_ACCURACY.
  • Added *MAT_221 options NERODE=10 and NERODE=11 to thick shell formulations 3, 5, and 7. They are available in bricks too.
  • Fixed laminated shell theory with shell *MAT_172.
  • Fixed *MAT_024 plastic strain failure for beam element forms 4 and 5. The elements were not failing.
  • Fixed thermal strains in *MAT_021 when used by thick shells. The material directions were not processed correctly causing incorrect thermal strains.
  • Fixed thick shell form 3 when used with hyperelastic materials. The L matrix was in the global coordinate system instead of the local system, resulting in an incorrect strain field.
  • Fixed shell formulation 18 when used with *MAT_058 . A divide by zero occurred.
  • Fixed ESORT>0 when used with shell formulation 18 and *MAT_077. The triangular element is now switched to form 17 because the other fully integrated triangular formulation is incompatible.
  • Fixed shell form 17 when used with *MAT_077. There was no call to the stress update routine.
  • Fixed the *EOS_GASKET model when used with thick shell form 2. It was giving the wrong stress.
  • Fixed the compressive failure strain (CFS) for *EOS_GASKET with bricks or thick shell formulations 5 and 7.  It was setting CFS=TFS.
  • Fixed an error in brick spot weld assemblies when used with strain rate sensitive failure functions. For assemblies, the strain rate is a volume weighted average strain rate of the elements in the assembly. However, the volume was likely wrong resulting in a strain rate that is likely too small. This error affected failure options OPT=1 and 7 on *MAT_100, and also weld failure for *MAT_100_DA.
  • For *MAT_261/*MAT_262, set default strengths to large values if they are not defined.
  • Enable GISSMO (*MAT_ADD_EROSION) to be used with beam element type 1.
  • Fix for GISSMO with non-zero values for SHRF and/or BIAXF. Rarely occuring problems with element failure were possible.
  • Fix for *MAT_034 (FABRIC) with FORM = 14, 24, or -14. Number of history variables could have been too low for the case that tables are used for parameters LCA, LCB, LCAA, or LCBB. And that could have led to noisy stress response.
  • Add plane stress option for *MAT_024 and shell elements. This "_2D" option invokes actual plane stress treatment, i.e, transverse shear stresses are not part of the yield condition but updated elastically.
  • Fix for rarely occurring negative plastic strain increments/rates in *MAT_224.
  • Add optional interpretation of EPSTHIN in *MAT_123. The sign of the input value governs the type of strain that is used:
    • EPSTHIN>0: thinning total strain (as in ISTUPD=1)
    • EPSTHIN<0: thinning plastic strain (as in ISTUPD=4)
  • Fix rarely occurring NaN issue with LOG_INTERPOLATION option of *MAT_024, *MAT_083, and *MAT_181.
  • Correct contact stiffness for *MAT_030.
  • Add option for *MAT_280: New parameter FTSCL on Card 2, field 8, is an optional scale factor for the tensile strength: FT_mod = FTSCL * FT. As soon as first crack happens in a part, tensile strength drops to its original value, FT.
  • Activate *MAT_280 for thick shell formulations 1 and 2.
  • Reduce the default number of history variables of *MAT_036 from 29 (as in R10.0) to 8 (as in versions before).
  • Add mid-plane failure option for GISSMO (MIDFAIL on *MAT_ADD_EROSION card 4).
  • *DEFINE_MATERIAL_HISTORIES was broken in early versions of R10 MPP, but is now fixed.
  • *DEFINE_FABRIC_ASSEMBLIES was broken in early versions of R10 MPP, but is now fixed.
  • Fix problem in which encrypted *MAT_075 input is echoed in d3hsp.
  • If an encrypted material is using load curve(s), the name of the material model is not echoed in d3hsp.
  • *MAT_293/*MAT_COMPRF is now available. This material models the behavior of pre-impregnated (prepreg) composite fibers during the high temperature preforming process. In addition to providing stress and strain, it also provides warp and weft yarn directions and stretch ratios after the forming process. The major applications of the model are for materials used in light weight automobile parts.
  • Fix 'floating overflow' error when TDEL=10.0E+20 (default) for *MAT_CONCRETE_BEAM/*MAT_195.
  • Fixed shell *MAT_091 when used with shell form 16. Also, *MAT_091 and *MAT_092 now support the CMPFLG=1 on *DATABASE_EXTENT_BINARY. The material directions were being calculated incorrectly causing the stress to be wrong. The scheme for calculating the material direction has been swapped out for one that is used with other orthotropic materials. This causes a very slight change in results with shell forms other than 16.

.

MPP
  • Store keyword partset info in the dyna.lsda file for later use in MPP decomposition, so partsets can be used in the pfile even if keyword->structured conversion is used before decomposition.
  • Update shell and solid decomposition cost parameters for the decomposition option "decomp { newcost }".
  • Decrease the memory footprint of MPP during keyword processing and decomposition.
  • Implement the checking for rigid body node replication already in SMP for MPP. Also added proper error termination for analyzing the constraint matrix in MPP.
  • Nodes defined in the following three keywords share the same database in LS-DYNA:
    1. *AIRBAG_REFERENCE_GEOMERY,
    2. *AIRBAG_SHELL_REFERENCE_GEOMERY,
    3. *INITIAL_FOAM_REFERENCE_GEOMETRY.

When *CONTROL_MPP_DECOMPOSITION_BAGREF is used, dummy with (3) is also treated in the special decomposition algorithm which gives very poor MPP performance. A fix was made to do the special decomposition to those shells with (1) or (2).

  • Allow use of parameters (defined via *PARAMETER) in *CONTROL_MPP_PFILE.
  • Fix a couple of issues in the MPP show decomposition option.

.

Output
  • Fix bug in output of first state to the d3part database when running double precision with the environment variable LSTC_BINARY set to 32ieee.
  • Fix for improper MPP input processing of set ID in nodfor file in some cases, which could have resulted in strange LEGEND sections in the nodfor output file.
  • Fix the marking of the end of the d3mode and d3eigv files when writing modal stresses.
  • Enabled cross section force output (secforc) for higher order shell forms 23 and 24.
  • Fixed stress output for shell forms 13, 14 and 15 when NIP=4 on *SECTION_SHELL and MAXINT<0 on *DATABASE_EXTENT_BINARY.  The stress output was incorrect.
  • Fixed an error in the output of fully integrated triangles with 7 or 9 through-thickness integration points. The data at the first integration point and strain at the last point were wrong.  This error had no effect on the solution.
  • Implement optional output of node sets and element sets used in *DATABASE_CROSS_SECTION. See variable ICRFILE in *CONTROL_OUTPUT. This allows for convenient display of those sets in LS-PrePost.
  • Fix bug in writing *SECTION_*_TITLE to d3prop file.
  • Fix missing frictional energy in glstat's sliding interface energy and total energy.

.

Restarts
  • Properly start up Explicit LaGrange Multiplier treatment of joints for restart.

.

SPH (Smooth Particle Hydrodynamics)
  • *ELEMENT_SPH_VOLUME allows user to input SPH particle volume instead of particle mass. This has the same effect as using plain *ELEMENT_SPH with MASS set to a negative value.
  • Fixed bug in reading of input for *DEFINE_SPH_TO_SPH_COUPLING_TITLE with optional second card.
  • Fixed incorrect output to binary (binout) rcforc in the case of *CONTACT_2D_NODE_TO_SOLID. Similarly, fixed incorrect output to ASCII rcforc (SMP only) in the case of multiple definitions of *CONTACT_2D_NODE_TO_SOLID.
  • Fix seg fault when using *CONTACT_AUTOMATIC_NODES_TO_SURFACE for SPH elements impacting on shell plate.

.

Thermal Solver
  • *DEFINE_CURVE_FUNCTION should now be working for thermal-solver-related commands that use a curve.

.

Miscellaneous
  • Fix incorrect handling of symmetric load curves when checking curve rediscretization errors, which resulted in bogus and misleading error messages.
  • Fix bug in *INTERFACE_LINKING that would only be hit if the user was using time-dependent functions for the y or z displacment scale factors.
  • Correct line count for error message IMP+31.
  • Fix a bug for PID controller function PIDCTL of *DEFINE_CURVE_FUNCTION. The bug occurred when using "0" as the sampling rate.
  • LS-DYNA now creates $JOBNAME.group_file, which can be loaded into LS-PrePost using Model>Groups so that groups of parts corresponding to each *SET_PART can be easily displayed.
  • Fix *DAMPING_RELATIVE when used with *DAMPING_PART_MASS on parts with shared nodes.
  • Fix long format errors for long=s option if the following keywords are used in the input deck:
*CONTACT_*_TITLE,
 *EOS_RATIO_OF_POLYNOMIALS,
 *CONSTRAINED_LAGRANGE_IN_SOLID_TITLE,
 *CONSTRAINED_BEAM_IN_SOLID,
 *CONTACT_ENTITY,
 *DEFINE_FILTER,
 *ELEMENT_NURB_SOLID_PATH,
 *ELEMENT_SHELL_OFFSET_COMPOSITE,
 *INITIAL(_STRESS,_STRAIN)_TSHELL,
 *MAT_002,
 *MAT_SPRING_MUSCLE,
 *MAT_INV_HYPERBOLIC_SIN,
 *SECTION_BEM, or
 if multiple include files are used.
  • Fix bug for *INCLUDE_TRANSFORM if material ID is input in alpha (non-numeric) format.
  • Fix bug in *INCLUDE_TRANSFORM if section offset is specified.
  • Fix bugs in *INTERFACE_SSI in reading of the keyword card and of intermediate files.
  • Fix option to dump all nodes into dynain; see NDFLAG in *INTERFACE_SPRINGBACK.

.