This file constitutes revision 0 of the release notes for LS-DYNA version R13.1.0.
.
The string “REVISION 13” must appear in the LS-DYNA license file in order to run version R13. Please contact your LS-DYNA distributor or your Ansys sales representative if you need to have your license updated.
.
Documentation for R13.1.0 is provided in the R13.0 User’s Manuals which may be downloaded from www.lstc.com/download/manuals. For features mentioned in these release notes but which appear to be missing from the R13.0 User’s Manuals, please refer to the DRAFT User’s Manuals at www.lstc.com/download/manuals. Please note that not all features in the DRAFT User’s Manuals are available in version R13.1.0.
The remainder of this file provides a brief description of new features, enhancements, and bug fixes in version R13.1.0. Some of the bug fixes may have also been included in R12.1.0.
The items are arranged by category. Understand that in many cases, a particular item may pertain to more than one category. In the interest of brevity, each item is listed only one, under a single category.
*CONSTRAINED_INTERPOLATION_SPOTWELD
(SPR3).*CONTACT_EXCLUDE_INTERACTION
which allows a set of segments to be excluded from contact with another set of segments or a set of segments to be excluded from all contact. The types of sets supported are segment sets, shell element sets, and part sets.*ICFD_MODEL_SPECIES_TRANSPORT
*ICFD_BOUNDARY_PRESCRIBED_SPTRANSP_CONC
*ICFD_INITIAL_SPTRANSP
*ICFD_DEFINE_SPTRANSPSOURCE
.*MAT_213/*MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE
is available. It is an orthotropic, combined visco-elastic-plastic and damage material model. You input tabulated stress-strain and damage curves with the option of providing quasi-static curves only or stress-strain and damage curves that are temperature and strain rate dependent. Element erosion criteria can be spcified with the Puck Failure Criterion, Tsai-Wu Failure Criterion, or a Generalized Tabulated Failure Criteria. Both solid and thin shell elements are supported. Probabilistic modeling of some material parameters is available. References are in Volume II of the User’s Manual.*VENDOR_LICENSE
scheme both with the FlexLM server and with the lstc_server licensing system.*AIRBAG_PARTICLE
:
*SENSOR_SWITCH_SHELL_TO_VENT
with IAIR = 1.*DATABASE_BINARY_CPMFOR
).*SENSOR_CPM_AIRBAG
.*DEFINE_CPM_VENT
to allow you to change blockage treatment for each vent individually.*DEFINE_PARTICLE_BLAST
.*CONSTRAINED_NODAL_RIGID_BODY_OVERRIDE
, added optional synchronization card to automatically turn off related contact, airbag, and prescribed motion when the overriding nodal rigid body is on and turn them back on when the related overriding rigid body is off.*CONSTRAINED_INTERPOLATION_SPOTWELD
can now refer to part sets, not just single parts. A value of less than zero indicates that it is a set, meaning PID1 < 0 references a part set with ID |PID1|. This feature could be used for in-plane composed sheets, such as Taylor welded blanks.*CONSTRAINED_INTERPOLATION_SPOTWELD
(SPR3) and *MAT_265_SPR3
for MODEL = 1. The peel ratio is the ratio of the bending moment to the resultant axial force. This value is kept constant when plastic yield is reached. Three fields can depend on the peel ratio:
*CONTACT_TIED_SURFACE_TO_SURFACE
for multistage analysis using FTYPE = 3 on *INTERFACE_SPRINGBACK_LSDYNA
.*CONTROL_ACCURACY
for compatibility between explicit and implicit when performing multistage analysis through FTYPE = 3 on *INTERFACE_SPRINGBACK_LSDYNA
. In particular, it invokes strongly objective tied contact algorithms, even for explicit analysis, in addition to selected “implicit” features.*CONTROL_TERMINATION
. A curve with ID = |ENDMAS| provides ENDMAS as a function of total mass.*DEFINE_DE_INJECT_BONDED
) for better consistency with changing core counts.*SECTION_SHELL
) or OFFSET (*ELEMENT_SHELL
), the time step size of those shell elements is reduced to fix instabilities. The reduction of the time step size is based on numerical tests which show a dependence on the offset distance and the ratio of shell thickness to edge length (T /L). Instabilities were most likely observed for aspect ratios of T /L > 0.5.*CONTROL_SHELL
to decide if the time step size is reduced when the shell reference surface is offset. NLOCDT = 0 is the default and reduces the time step size to fix instabilities. NLOCDT = 1 restores the prior behaviour and does not reduce the time step size at the risk of instabilities.*ELEMENT_BEAM
to work with beam form 9 (spot welds).*SECTION_SOLID
) to suppress the influence of non-physical resistance.*SECTION_SHELL
) or OFFSET (*ELEMENT_SHELL
), the time step size of those shell elements is reduced to fix instabilities. The reduction of the time step size is based on numerical tests which show a dependence on the offset distance and the ratio of shell thickness to edge length (T /L). Instabilities were most likely observed for aspect ratios of T /L > 0.5.*SECTION_SOLID
. The first new parameter is cohesive thickness, COHTHK. This allows section / partwise definition of cohesive thickness. This value supersedes the one from *MAT_240
or *MAT_ADD_COHESIVE
.*SECTION_SHELL_XFEM
).*FREQUENCY_DOMAIN_SSD
:
*CONTROL_IMPLICIT_BUCKLE
or *CONTROL_IMPLICIT_EIGENVALUE
.*CONTROL_IMPLICIT_SOLUTION
. This option ensures that prescribed motions are automatically satisfied, regardless of step size. This should avoid the infamous message “Convergence prevented due to unfulfilled bc…”*IGA_SHELL
, check the trimming curves and limit the generated sampling points to the parametric space of the underlying surface description. This reduces bad parameterizations given by “dirty” geometry definitions.*IGA_SHELL
and *IGA_SOLID
) to be treated as rigid bodies.*IGA_SHELL
and *IGA_SOLID
) erosion: Add element-type information to the message file specifying IGA elements and related interpolation elements.*IGA_POINT_UVW
is within a valid knot-span. If the location of an *IGA_POINT_UVW
(connected to a *NODE
) is not within a valid knot-span of the underlying face, a warning message will be written and this point/node pair will be removed.*LOAD_NODE
to evaulate the load function only once. The value is stored and applied for the rest of the run.*LOAD_SEISMIC_SSI_DECONV
for unstructured meshes or for meshes with small element height.*MAT_ADD_DAMAGE_GISSMO
. |LCREGD| can now refer to a *DEFINE_TABLE_3D
. This makes the regularization factor a function of Lode parameter, triaxiality, and element size.*MAT_ADD_DAMAGE_GISSMO
. |FADEXP| can refer to *DEFINE_TABLE
(fading exponent is a function of element size and triaxiality) or *DEFINE_TABLE_3D
(fading exponent is a function of element size, triaxiality, and Lode parameter).*MAT_ADD_DAMAGE_GISSMO
. This flag governs the behavior of instability measure, F , and fading exponent, FADEXP:
*MAT_ADD_DAMAGE_GISSMO
by adding new function arguments: strain rate, temperature, history variable, and element size. So now we have epf = epf(triaxiality, Lode parameter, plastic strain rate, temperature, history, element size).*MAT_ADD_DAMAGE_GISSMO
. Therefore, parameters LCSDG > 0 and ECRIT < 0 can now refer to *DEFINE_TABLE_4D
to make failure / critical strain a function of plastic strain rate (TABLE_4D), temperature (TABLE_3D), Lode parameter (TABLE), and triaxiality (CURVE).*MAT_ADD_GENERALIZED_DAMAGE
with IFLG2 = 1. It now supports the following additional material models for solid elements: *MAT_133
, *MAT_199
, and *MAT_233
.*MAT_ADD_INELASTICITY
. The improved algorithm is more efficient and robust.*DEFINE_CONNECTION_PROPERTIES
: types *MAT_003
, *MAT_081
, *MAT_133
, *MAT_187
, *MAT_224
, and *MAT_243
.*MAT_HYPERELASTIC/OGDEN_RUBBER
(*MAT_077
) compatible when used with multistage analysis using FTYPE = 3 on *INTERFACE_SPRINGBACK_LSDYNA
.*MAT_FU_CHANG_FOAM
(*MAT_083
) for implicit/explicit multistage analysis. Explicit and implicit were using different algorithms for damage causing an incorrect transition between stages. This is remedied now by setting IACC = 2 on *CONTROL_ACCURACY
.*MAT_083
. This is a new path dependent, incremental formulation for foams. It supports *DEFINE_TABLE_3D
for TBID with additional dependence on the volumetric strain (1 – J).*MAT_091/*MAT_092
modified initialization of material directions for solid elements when FANG is used. There was an invisible option BETA for many years in field 8 of card 3. Before this update solids used the invisible BETA, but shells used FANG to input a material rotation angle. With this update, FANG will be used by solids too if it is nonzero. However, if FANG = 0.0, then the hidden BETA will be used for backward compatibility.*MAT_SIMPLIFIED_JOHNSON_COOK_ORTHOTROPIC_DAMAGE/*MAT_099
to avoid numerical issues such as division by zero and floating point overflows.*MAT_187
(SAMP-1) and *MAT_187L
(LIGHT) with LCEMOD > 0. VERM can activate / deactivate (0 / 1) the equilibrium yield condition.*MAT_THERMO_ELASTO_VISCOPLASTIC_CREEP
(*MAT_188
).*MAT_GENERAL_SPRING_DISCRETE_BEAM/*MAT_196
which activates spot weld thinning when discrete beams that use the material are tied to shell elements. This behavior was turned on by default, but now this parameter must be used to activate it.*MAT_224
. If keyword option LOG_INTERPOLATION is used, and |BETA| refers to a *DEFINE_TABLE_4D
, then logarithmic interpolation is used in the strain rate direction.*MAT_240
: exponential moving average of the effective plastic strain rate (RFILTP>0) or of the effective total strain rate (RFILTP<0).*MAT_240
to govern yield under compression (COMPY) and enable smooth treatment of asymptotic limits (SMOLIM). These fields are currently only available if none of the keyword options THERMAL, 3MODES, or FUNCTIONS is used.*MAT_UHS_STEEL
(*MAT_244
) for implicit. The transition between transformation and plasticity caused a discontinuity in the response.*MAT_252
to apply an exponential moving average.*MAT_258
: If input parameters A and GAMMA (exponents) are integers, then use them like this internally as well (real exponents are quite expensive).*MAT_BERGSTROM_BOYCE_RUBBER
(*MAT_269
) that should reduce numerical roundoff sensitivity.*MAT_280
to allow FT < 0 to be used with FTSCL > 0. In that case, |FT| defines a curve for tensile strength vs. strain rate, and FTSCL scales the strength values from that curve as long as material is intact. If cracked, neighbors get non-scaled values from that curve.*CONTROL_MPP_REBALANCE
).*INCLUDE_COSIM
) now uses a simplified MPI command as follows:*INCLUDE_COSIM
) is supported.*DATABASE_BINARY_D3DUMP
, *DATABASE_BINARY_RUNRSF
, *DATABASE_BINARY_DEMFOR
, and *DATABASE_BINARY_INTFOR
.*DATABASE_BINARY_D3DUMP
.*DATABASE_EXTENT_BINARY
.*INITIAL_STRAIN_TSHELL
with *INTERFACE_SPRINGBACK_EXCLUDE
.*CONTROL_OUTPUT
). We now replace & , ” , < , > , and ‘ with & , " , < , > , and ' to avoid problems with XML file readers.*DATABASE_EXTENT_BINARY
to write last available results to d3plot and d3part for deleted elements.*CONSTRAINED_INTERPOLATION_SPOTWELD
) to binout / swforc database: “resultant_axial” (manual: F n ) and “resultant_shear” (F s ).*DEFINE_CURVE_TITLE
to d3hsp .*DEFINE_COORDINATE_..._TITLE
to d3hsp .*SENSOR_DEFINE_FORCE
.*PARAMETER
and *PARAMETER_EXPRESSION
. If “_NOECHO” is appended, the defined parameters will not be echoed to the d3hsp file.*VENDOR_BEGIN/*VENDOR_END
licensing. The string CUST_DATE is in the first field, the MM/DD/YYYY expiration date is in columns 11-20, and a string customer ID is in columns 21-80. With the existing LST license types (lstc_server and nodelocked), the customer ID will be checked against a program OPTION named “CUSTOMER_NUM”. If they do not match (or CUSTOMER_NUM is not set), the VENDOR block will be skipped. We plan to eventually add this for the FlexLM licenses.*NODE_TRANSFORM
. This field allows you to control the processing sequence of *NODE_TRANSFORM
.*DEFINE_TABLE_COMPACT
.*DEFINE_CURVE_FUNCTION
. They now work with EXACC = 1 on *CONTROL_ACCURACY
.*FREQUENCY_DOMAIN_FRF
, fixed a bug in reading the d3eigv database when modal stress is present. The bug caused an error termination.*FREQUENCY_DOMAIN_ACOUSTIC_BEM
in MPP with predecomposition. The bug resulted in an error termination.*CONTROL_ADAPTIVE
:
*CONTROL_ADAPTIVE
). The bug caused the boundary condition to be applied incorrectly.*AIRBAG_WANG_NEFSKE
when used with *AIRBAG_INTERFACTION
. The bug might result in a segmentation fault.*SENSOR_CPM_AIRBAG
birth time shift to make it work for all shell element formulations that work with *MAT_034
. It was only working for element formulation 12. It should now work for shell element formulations 12, 13, 14, 22, 23, and 24.*ALE_STRUCTURED_FSI
:
*ALE_STRUCTURED_MESH_MESH_MOTION
). The bug caused the detonation time to be kept the same after advection.*DEFORMABLE_TO_RIGID_AUTOMATIC
with *BOUNDARY_PRESCRIBED_MOTION
. This bug intermittently caused segmentation faults. The SMP version does not have this issue.*BOUNDARY_PZEPOT
are in the new or old format (NSID, LCID, SF, ID). The check was initially just based on the existence of ‘.’ in the SF field of the old format, but this check missed certain cases, particularly when the input deck came from LS-PrePost.*BOUNDARY_PRESCRIBED_MOTION_RIGID(_LOCAL)
). It was not implemented but included in the documentation.*CONSTRAINED_INTERPOLATION
in single precision that were leading to an error termination.*CONSTRAINED_GENERALIZED_WELD_BUTT
which was created when *CONSTRAINED_NODAL_RIGID_BODY_OVERRIDE
.*CONSTRAINED_NODAL_RIGID_BODY_OVERRIDE
, fixed an MPP bug that happens when trying to resume a contact that was initially turned off.*CONSTRAINED_INTERPOLATION
in single precision that were leading to an error termination.*CONSTRAINED_SOIL_PILE
for elements on the boundaries between MPP domains.*CONSTRAINED_NODAL_RIGID_BODY
and *CONSTRAINED_NODAL_RIGID_BODY_INERTIA
. Incorrect logic caused the inertia properties to be put on the wrong nodal rigid bodies.*CONSTRAINED_JOINT_GEARS
in implicit. It was not being accounted for.*CONSTRAINED_JOINT_RACK_AND_PINION
in implicit.*CONSTRAINED_INTERPOLATION
: failure flag in swforc was not showing failure correctly.*CONSTRAINED_INTERPOLATION_SPOTWELD
). Failed beam elements in a deformed model were flying away. Now they stay with the disconnected sheets.*CONSTRAINED_INTERPOLATION_SPOTWELD
with MODEL > 10, where selected material parameters are defined as functions (*DEFINE_FUNCTION
). This did not work correctly for function IDs > 9999.*CONTACT_TIEBREAK_SURFACE_TO_SURFACE_ONLY
in MPP.*CONTACT_AUTOMATIC_TIEBREAK
with OPTION = 4 when the tied nodes slide off the edge of the constraining surface.*CONTACT_TIED_SHELL_EDGE_TO_SOLID
which was not computing the master side energy correctly.*CONTACT_..._INTERFERENCE
which was basically bypassing the interference part completely.*CONTACT_RIGID_SURFACE
initialization which has in some cases been using the wrong coordinates in the single precision version of the LS-DYNA starting with R10.0.*CONTACT_OPTION
.*CONTACT
when the following conditions are met:
*CONTACT
). The values were too low.*DATABASE_EXTENT_BINARY
). Tributary areas for the nodal forces were not being accounted for, so the calculated peak was too high.*CONTACT
option to use SSID = 0 with single surface contact to include all parts when *ELEMENT_SHELL_COMPOSITE
is used in the model.*CONTACT_AUTOMATIC_NODES_TO_SURFACE
, the slave side contact thickness offset was considered 0 unless SHLTHK was set to 2 on *CONTROL_CONTACT
. Since it is an automatic contact, the contact thickness offset should always be considered. The contact thickness offset is now considered regardless of the setting of SHLTHK.*CONTACT_SURFACE_TO_SURFACE
. This contact is two-way, but one-way contact was being performed.*CONTACT_AUTOMATIC_BEAMS_TO_SURFACE
) for triangular segments. The initial penetrations were not computed correctly which led to incorrect penetration reports and wrong contact forces. For larger models mainly comprised of quadrilateral segments, the error was hardly noticeable.*CONTROL_OUTPUT
is used. Too many (incorrect) penetrations were reported in explicit analysis.*CONTACT_OPTION
(MPP only).*CONTACT_AUTOMATIC_BEAMS_TO_SURFACE
. Beams with *MAT_NULL
and RO < 1.e-11 get zero nodal mass, but that could lead to later problems in contact. It is now resolved.*CONTROL_ACCURACY
, making sure stresses are exactly zero at time zero of a simulation.*CONTROL_TIMESTEP
).*CONTROL_ENERGY
: detailed material energies were not transferred from one step to the next.*CONTROL_ENERGY
: for solid elements, energy results could have been wrong.*CONTROL_TIMESTEP
and *MAT_187
leading to problems with curve generation.*CONTROL_TIMESTEP
.*DEFINE_DE_TO_BEAM_COUPLING
to get the correct beam thickness.*PART_AVERAGED
, and prevented execution of related code on processors that have none of these beams (MPP only).*SECTION_SHELL
). The issue could cause a seg fault.*DEFINE_ELEMENT_EROSION_TSHELL
to work with thick shell forms 3, 5, and 7 that use material type 54 with DFAILT and DFAILC.*CONTROL_SHELL
is set to sort triangular elements into triangle forms, and the sorting changes the number of in-plane integration points or the ordering of the data. This feature was not working.*PART_STACKED_ELEMENTS
: memory overwrite. Models could fail with “Error 10183 (KEY+183), part … not defined,” even if that part did not exist at all.*FATIGUE_MODAL_DYNAMIC
). The control data for the number of variables for each element was not correct.*FATIGUE_MULTIAXIAL
) with solid elements was not being performed when requested.*ICFD_MODEL_POROUS
). It was using the permeability instead of the porosity.*ICFD_MODEL_POROUS
). A factor of 1/2 was introduced while testing around mid-2020.*CONTROL_IMPLICIT_EIGENVALUE
).*CONTROL_IMPLICIT_MODAL_DYNAMICS
to correctly use *LOAD_BODY_OPTION
.*CONSTRAINED_RIGID_BODY_STOPPER
.*CONTROL_IMPLICIT_GENERAL
. It was also not accounting for the time step limitations on *CONTROL_IMPLICIT_AUTO
.*INITIAL_AXIAL_FORCE_BEAM
used with KBEND = 1 in implicit. If the beam did not align with the X, Y , or Z-directions in the global coordinate system, the results were not consistent.*INITIAL_AXIAL_FORCE_BEAM
when KBEND = 1. The transformation between the local and global system was not being accounted for correctly.*INITIAL_VELOCITY_GENERATION
with ID = 0 and STYP = 0 with another *INITIAL_VELOCITY_GENERATION
with nonzero ID and STYP.*INITIAL_VELOCITY_GENERATION
yields the wrong velocity for *CONSTRAINED_NODAL_RIGID_BODY_INERTIA
when CID is nonzero.*IGA_SOLID
was not passed to the subroutine which calculates the internal energy reported in glstat .*DATABASE_CROSS_SECTION
) occurred when the cross section had more than one IGA patch. It worked correctly for one patch.*IGA_TIED_EDGE_TO_EDGE
, where more than two faces are tied along the same edge. The bug showed up only in MPP. Prior to this fix, the tying behavior might be wrong, depending on the number of cores used.*IGA_SHELL
and *IGA_SOLID
) that may have led to error termination for very large models.*IGA_SHELL
and *IGA_SOLID
). Prior to this fix, models with surface-to-surface or node-to-surface contacts may have terminated unexpectedly.*IGA_SHELL/SOLID
and *ELEMENT_SHELL/SOLID_NURBS_PATCH
. This bug only showed up for a multi-patch analysis with different polynomial degrees in the patches.*IGA_SHELL
:
*RIGIDWALL
in MPP together with IGA elements. Prior to the fix, if some processors did not contain IGA elements, deadlock could occur.*LOAD_VOLUME_LOSS
did not work correctly when COORD is nonzero, and the coordinate system had certain characteristics. The symptom was warning message SOL+1127 in the message file.*LOAD_SEISMIC_SSI_DECONV
. The errors indicated that the load curves had unequal time steps (equal time steps are required) when they did not.*LOAD_SEISMIC_SSI
in SMP or hybrid. It was giving different results for different values of ncpu < 0.*MAT_ADD_DAMAGE_GISSMO
with LCSDG being a *DEFINE_TABLE_3D
. The first value in the sub-table was interpreted as logarithmic strain rate, and logarithmic interpolation was used in that case. But it is the Lode parameter and interpolation should be linear.*MAT_ADD_DAMAGE_GISSMO
and *PART_COMPOSITE
. If the first layer in the composite is a non GISSMO material, failure treatment could have been wrong.*MAT_ADD_EROSION
. LS-DYNA issued a missing curve error.*MAT_ADD_INELASTICITY
. If used in conjunction with a plastic material in the first place, it was incorrectly overwritten.*INCLUDE_TRANSFORM
) for PERMi / LCPGDi on *MAT_ADD_PORE_AIR
.*MAT_ADD_PZELECTRIC
that occurred when the eigenvalue calculation was requested. The bug caused a seg fault.*MAT_005
to reduce noisy stresses.*MAT_022
and *MAT_054
which is used in implicit analysis. Previously, if elements deviated from a perfect square shape the stiffness and results like eigenvalues changed.*MAT_PIECEWISE_LINEAR_PLASTICITY/*MAT_024
with VP = 1 when used with tetrahedra elements (ie tet13) after dynamic relaxation. The problem showed up when the part reached plastic deformation during dynamic relaxation.*MAT_PIECEWISE_LINEAR_PLASTICITY
(*MAT_024
). For certain tables the incorrect value of yield stress was extracted.*MAT_024
(shells) with TDEL > 0 and GISSMO. History variables conflict could lead to wrong results.*MAT_024
with VP = 3 for H-L beam elements now supports LCSS being a *DEFINE_TABLE_{X}D
with dimension X > 2. Prior to this fix, 3D tables and larger only worked for shells and solids.*MAT_030
where plasticity routine failed to converge for shell elements.*MAT_033
with N = 0.*MAT_037/*MAT_TRANSVERSELY_ANISOTROPIC_ELASTIC_PLASTIC
, when used for shells with thickness stretch, there was a bug in the numerator for plastic strain update, leading to incorrect results.*MAT_LAMINATED_COMPOSITE_FABRIC
(*MAT_058
), if EA and EB were specified as *DEFINE_CURVES
, the scale and offset factors (SFA,OFFA) were not taken into account in an error check.*MAT_077
that occurred for shell elements only. The bug caused a seg fault.*MAT_HYPERELASTIC/OGDEN_RUBBER
(*MAT_077
).*MAT_079
(*MAT_HYSTERETIC_SOIL
) exhibited spurious stress changes for elements that rotate significantly. Errors are small for typical seismic or civil/structural applications but could potentially become problematic, such as in highly deforming regions of the soil.*MAT_122_3D
was not working for this material. Also added missing d3hsp echo of material parameters for this material.*MAT_126
and *MAT_ADD_THERMAL_EXPANSION
. It only worked for solid ELFORM = 0, but now it also works for other element formulations.*MAT_127
: variable N (number of statistical links) should be > 0.*MAT_135
: the manual mentions files Contour_1# etc. in Remark 2, but they have been output as empty files for many years. The files now contain data again. Also, added jobid to the file names.*MAT_WOOD/*MAT_143
.*MAT_157_IHIS
.*MAT_161
and *MAT_162
licenses that was broken some time ago. Also, changed the required environment variable which for some strange reason contained the license version number embedded in it, meaning it was different for every version of the code. Now these are set to ANSYS_LSMAT161 and ANSYS_LSMAT162.*MAT_172
(*MAT_CONCRETE_EC2
) with TYPEC= 1, 2, 4, 5, 7, or 8. For these settings of TYPEC, the mechanical properties of the concrete vary with temperature, in particular the compressive strength reduces with increasing temperature. When subjected to increasing temperature under confined/restrained conditions, there is a complex situation whereby the compressive stresses build up due to thermal expansion, but the compressive strength is reducing due to the effect of the increasing temperature on the mechanical properties. In this situation, the bug caused the compressive stress to build up more slowly than expected.*MAT_173
(*MAT_MOHR_COULOMB
) calculates a tensile pressure cut-off which is intended to occur at the stress state where the yield surface reaches zero shear stress. This was done correctly for friction angle PHI > 0, but when PHI = 0 (meaning purely cohesive behavior), the tensile cut-off was incorrectly set to zero. Now, an arbitrary very large tensile cut-off is applied when PHI = 0.*MAT_173
and *MAT_198
, curves used for these material were not being passed correctly for solid element form 62. If the curves were not in use, then these materials worked with this element formulation.*MAT_186
: variable TSLC2 was not working in free format or with *PARAMETER
.*MAT_4a_MICROMEC/*MAT_215
was not correct.*MAT_224
combined with *LOAD_THERMAL
. Temperatures in hisvar #14 (solids) and #10 (shells) can be set with *INITIAL_STRESS_SOLID/SHELL
. Simultaneous use of *LOAD_THERMAL
on other parts in the model should not alter these values.*MAT_244
to use the same Cowper Symonds expression as *MAT_106
.*MAT_248
to be consistent with *MAT_106
:
*MAT_254
:
*MAT_255
added missing integration point coordinates needed for temperature.*MAT_280
and *PART_COMPOSITE
. Some history variables (1,2,3) are copied across thickness and that should not be the case for integration points with non-glass material.*CONTROL_MPP_REBALANCE
) that were resulting in seg faults.FORCE_TRANSDUCER
type contacts, which were not previously supported.*TERMINATION_BODY
which was never properly implemented in MPP.*INCLUDE_COSIM
).*INCLUDE_COSIM
.*CONTROL_ADAPTIVITY
):
*DATABASE_HISTORY_NODE_SET
if i10 format is used ( i10=yes on the command line).*DATABASE_RECOVER_NODE
that occurs for solid elements with more than one integration point. The results were incorrect.*DATABASE_HISTORY_SEATBELT
) for the specified elements.*DATABASE_SECFORC
that caused SECFORC force to be too large. This issue occurred for MPP with more than one cpu.*DEFINE_MATERIAL_HISTORIES
to be used with Rayleigh damping. The history data and damping variables were conflicting, causing early failure of elements.*DATABASE_BINARY_D3THDT
) was corrupt if NINTSLD = 8 on *DATABASE_EXTENT_BINARY
. Element stress output to elout and binout files was unaffected by the bug.*MAT_FABRIC
should not be written to the *INITIAL_STRESS_SHELL
cards in the dynain file. This was done correctly for dynain files triggered by the *INTERFACE_SPRINGBACK
card, but not for dynain files triggerd by the staged construction keywords (see *CONTROL_STAGED_CONSTRUCTION
and *DEFINE_CONSTRUCTION_STAGES
).*INTERFACE_SPRINGBACK
, for some anisotropic material models, the requested *ELEMENT_SOLID_ORTHO
was not written properly.*DATABASE_SSSTAT
) and could cause errors.*DEFINE_MATERIAL_HISTORIES
:
*DEFINE_MATERIAL_HISTORIES
is used with *MAT_FABRIC
(*MAT_034
) to make sure that transformation quantities stored among history variables are not affected or affecting results.*DATABASE_EXTENT_BINARY
.*DATABASE_D3MAX
, fixed bugs in the single precision versions which caused crashes in MPP runs and 0 von Mises stress in SMP runs.*CONTACT_AUTOMATIC_GENERAL
. This has been broken since January of 2019.*DATABASE_BINARY_INTFOR
and full deck restart:
*CONTROL_PORE_FLUID
) did not work with full-deck restarts.*SENSOR_DEFINE_FORCE
. It occurred when the related *DATABASE_CROSS_SECTION
was not defined. The bug was leading to a seg fault.*SENSOR_CONTROL
. It occurred when a core did not own any part of the rigid body.*SENSOR_DEFINE_FORCE
. The results were incorrect when more than one core was used.*SENSOR_CONTROL
that occurred when TYPE = ELESET and ESTYP = DISC. The bug cuased discrete elements to not be deleted as designed.*INCLUDE_TRANSFORM
with *SENSOR_DEFINE_FORCE
when TYPEID = X-SECTION, CONTACT, CONTACT2D, JOINT, JOINTSTIF, RWALL or SPC and ID offsets are used for the transformed parts.*CONTROL_THERMAL_SOLVER
to d3hsp.*CONSTRAINED_ADAPTIVE
for constrained solid / beam nodes in case thermal thick shells are also present in the model. This bug could lead to LS-DYNA crashing.*INTERFACE_LINKING
interface cannot be found in the interface file from the previous run.*SET_SEGMENT_GENERAL
. It could have resulted in incorrect values being applied for the DA1-DA4 parameters, or possibly a seg fault during input processing.*TERMINATION_DELETED_SOLIDS
. The simulation was sometimes continuing to run even though the termination criteria was met.*DEFINE_FUNCTION
definitions.*PARAMETER
) were not conforming to the usual convention for comma-delimited input.*PARAMETER
if long format is used.*PART
type and *ELEMENT
type against the formulation defined in *SECTION
.*PARAMETER
) in a larger string as part of an include file name. The bug caused an extra space in the file name which led to an error termination.*DEFINE_CURVE_FUNCTION
,
*INCLUDE_TRANSFORM
,
*PART_INERTIA
or *CONSTRAINED_NODAL_RIGID_BODY_INERTIA
. The bug was leading to an incorrect warning message.*DEFINE_TRANSFORM
. The model was not correctly transformed.*NODE_TRANSFORM
and *SET_NODE_GENERAL
. It occured when BOX or VOL wass used for *SET_NODE_GENERAL
in a model that also included *NODE_TRANSFORM
. The nodes were not being transformed correctly.*MODULE_USE
did not work correctly for user material numbers in the extended range 1001 to 2000.*MODULE_PATH
or LD_LOAD_LIBRARY_PATH
did not work correctly.*MODULE_LOAD
or the command line option ” module= ” could not be used with the default usermat shared object filename.*DEFINE_PRESSURE_TUBE
where automatic generation of type 25/26 shells did not work.*RIGIDWALL
code to prevent memory errors in very large models. Small models were working.*RIGIDWALL_PLANAR_MOVING
from structured file in case of long = s. Mass and velocity were not being read back correctly.*PART_COMPOSITE
as a label.*DEFINE_DE_ACTIVE_REGION
with TYPE = 1 (box) and long=s: coordinates were read in the wrong format from the structured input.