This file constitutes revision 0 of the release notes for LS-DYNA version R12.1.0.
The string ‘REVISION 12’ must appear in the LS-DYNA license file in order to run version R12.1.0. Please contact your LS-DYNA distributor or your Ansys sales representative if you have to have your license updated.
Documentation of R12.1.0 is provided in the R12.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 R12.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 R12.1.0.
The remainder of this file provides a brief description of new features, enhancements, and bug fixes in version R12.1.0. Some of the bug fixes may have also been included in other recent releases. Many of the new features and enhancements are also in R13.0.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_SPR2 or by *CONSTRAINED_INTERPOLATION_SPOTWELD (SPR3).*EM_RANDLES_EXOTHERMIC_REACTION to batmac (*EM_RANDLES_BATMAC).*DEFINE_MATERIAL_HISTORIES.*DATABASE_FREQUENCY_ASCII_NODOUT_SSD) for direct SSD (*FREQUENCY_DOMAIN_SSD_DIRECT).*NODE or *ELEMENT data.*CONTROL_ADAPTIVE) is not tetrahedrons.*CONSTRAINED_BEAM_IN_SOLID definitions.*CONSTRAINED_BEAM_IN_SOLID_PENALTY.*BOUNDARY_AMBIENT AMBTYP = 5 (*LOAD_BLAST_ENHANCED), we’ve made the following enhancements:
*LOAD_BLAST_ENHANCE for better accuracy*INITIAL_VOLUME_FRACTION_GEOMETRY to fill ALE elements with the material of a solid part in 3D (or shell part with 2D). The boundaries of the part specify the geometry of the container to be filled. The material in the part can optionally fill this container. This feature is useful for setting up an ALE model with a complex dummy structure that outlines a volume of ALE elements to fill.*CHEMISTRY with *CESE.*CONSTRAINED_INTERPOLATION, like for instance if independent nodes are co-aligned and thus unable to represent a rotation with respect to the corresponding axis.*CONSTRAINED_INTERPOLATION_SPOTWELD can now refer to part sets, not only single parts. Set to negative part ID number to do so, that is, PID1 < 0 and/or PID2 < 0. This could be used for in-plane composed sheets, such as Taylor welded blanks.*CONSTRAINED_INTERPOLATION_SPOTWELD (SPR3):
*CONSTRAINED_TIED_NODES_FAILURE now supports damage values from *MAT_107 and *MAT_110 to be used as failure variable for EPPF.*CONSTRAINED_BEAM / SHELL_IN_SOLID are the PENALTY version in the d3hsp file.*CONSTRAINED_BEAM_IN_SOLID and *CONSTRAINED_SHELL_IN_SOLID in the same way as for *CONTACT_OPTION.*CONTACT_AUTOMATIC_GENERAL*DEFINE_FRICTION.*PART_CONTACT is used. If a segment pair in contact has nonzero friction parameters from only one of the segments, the nonzero friction parameters will be used for the pair. If friction parameters are defined for both parts, then the average is used. Previously, the average was used without checking values, so zero values could be used in the average causing the Coulomb coefficient to by 1/2 of the expected value.*CONTACT_OPTION to work with MPP contact when SOFT = 0 or 1.*ELEMENT_SHELL_SOURCE_SINK segment-based to work with contact.*CONTACT_TIEBREAK_NODES_TO_SURFACE. TIEDIED = 1 enables incremental displacement update of the contact normal. Support was added for this contact to fix spurious stresses even with no loading.*CONTACT_OPTION to scale the computed volume of tet elements for the mass calculation in SOFT=2 contact. Better agreement is achieved with TETFAC between 3 and 5 when contact forces are compared between hex and tet meshes. This input field is supported for 4-node and 10-node tets.*DEFINE_STAGED_CONSTRUCTION_PART and *CONTROL_STAGED_CONSTRUCTION) did not work for thick shells (*ELEMENT_TSHELL) with ELFORM = 3 or 5. Previously, the staged construction keywords had no effect on these elements, meaning they were always active.*CONTROL_TIMESTEP.*CONTROL_CONSTRAINED for setting constraint related properties. Currently it only sets the behavior of the initialization check for SPR2 / SPR3.*CONTROL_ENERGY. Currently this feature is supported for material models *MAT_003, *MAT_004, *MAT_015, *MAT_019, *MAT_024, *MAT_063, *MAT_081, *MAT_082, *MAT_098, *MAT_104, *MAT_105, *MAT_106, *MAT_107, *MAT_123, *MAT_124, *MAT_188, *MAT_224, *MAT_225, *MAT_240, and *MAT_251 for shell and solid elements. Setting MATEN = 2 causes the detailed material energies to be computed. For this calculation, internal energy is split into the elastic, plastic, and damage portions. They are computed and reported as ‘mat_energy_elastic ‘, ‘mat_energy_plastic ‘, and ‘mat_energy_damage ‘ in glstat and matsum*CONTROL_MAT to potentially issue an error if only unmodified user material subroutines are called.*CONTROL_SHELL to define a part set, for which the stresses from *INITIAL_STRESS_SHELL are with respect to the local coordinate system.*CONTROL_TERMINATION: A curve with ID = |ENDMAS| provides ENDMAS as a function of total mass.*CONTROL_SOLID for *MAT_169. This is automatic deletion if neighboring elements fail. It is already available for cohesive elements but now also available for this special material.*MESH capabilities to permit use of *MESH_BL, etc.*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.*ELEMENT_BEAM to work with beam form 9, spot weld beams.*ELEMENT_DIRECT_MATRIX_INPUT:
*PART_COMPOSITE. This can be defined via an irregular optional Card 2 (OPTCARD).*SECTION_SHELL with first new parameter THKSCL. This allows the scaling of shell element thickness for all elements in the corresponding section at once.*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_SOLID, ELFORM = 19 / 20) to detect inverted elements in the first cycle.*SECTION_SOLID_EFG is now the default.*INCLUDE_STAMPED_PART.*DEFINE_CURVE_STRESS allows you to change the number of output stress-strain pairs based on the value of LCINT in *CONTROL_SOLUTION.*ICFD_BOUNDARY_FSWAVE.*CONTROL_IMPLICIT_EIGENVALUE for MPP execution. These controls were initially only available for SMP.*PART_MODES to work in single precision.*CONTROL_IMPLICIT_INERTIA_RELIEF. These contacts are now enabled.*CONTROL_IMPLICIT_EIGENVALUE, rewrote how data is output to the d3plot , d3eigv , and similar files, to vastly improve performance of the output. This saves 15*CONTROL_IMPLICIT_MODAL_DYNAMICS) when reading the modes from d3eigv .*INITIAL_STRESS_SOLID to work when the data has more points than the element form has. For example, the data has stress data for 8 stress points, but the element has only 1 or 2 integration points.*INITIAL_VELOCITY_GENERATION to work with solid element forms 16, 17, 45, and 61. Only the first 8 nodes were being initialized.*INITIAL_HISTORY_NODE(_SET). So far ELFORMS = -2, -1, 1, 2, 10, and 13 are supported.*CONTROL_PARALLEL, available for *IGA keywords.*LOAD_THERMAL_VARIABLE_NODE in case of dynamic relaxation. Now we accept that the keyword is defined twice, where in one definition the load curve is active during DR and in the other it is active during the transient phase.*MAT_ADD_COHESIVE now supports hyperelastic foam materials *MAT_057 and *MAT_083.*MAT_ADD_EROSION so that it works correctly with failure criteria in *MAT_ADD_DAMAGE_DIEM/GISSMO and *MAT_ADD_EROSION itself. Before it only worked with failure in materials, such as FAIL in *MAT_024.*DEFINE_HAZ_PROPERTIES. It is necessary to use *MAT_ADD_DAMAGE_GISSMO_STOCHASTIC to make it work.*MAT_ADD_GENERALIZED_DAMAGE.*MAT_ADD_GENERALIZED_DAMAGE.*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_024 option VP = 3 in implicit analysis for solid and shell elements.*MAT_024 with VP = 3 now supports *DEFINE_TABLE_{X }D / COMPACT up to a level of 9, meaning yield stress can be a function of plastic strain, strain rate, and up to seven history variables. If *MAT_ADD_DAMAGE_DIEM is used together with *MAT_024, then P1 of *MAT_ADD_DAMAGE_DIEM also supports a table up to level 9.*MAT_034 (FORM = 14, -14, and 24).*MAT_036 as potential weld partner material for PROPRUL = 2/3 of *DEFINE_CONNECTION_PROPERTIES.*MAT_068:
*MAT_SIMPLIFIED_JOHNSON_COOK (*MAT_098). It should improve robustness.*MAT_100 and *MAT_100_DA. In addition to setting the transverse shear stresses to zero (E < 0), the transverse shear strain increments get zeroed as well. This can reduce occasional instabilities.*MAT_100_DA / *MAT_SPOTWELD_DAIMLERCHRYSLER:
*MAT_JOHNSON_COOK / *MAT_110 for shells.*MAT_123 (VP = 1) for shell and solid elements. Viscoelastic approach adopted from *MAT_187 with LCEMOD > 0.*MAT_124. This enables strain rate filtering using exponential moving average. It works for LCSRC or LCSRT being nonzero and SRFLAG = 0 or 1.*MAT_136 to improve response for kinematic and mixed hardening.*MAT_138 as element size depende4tn curves.*MAT_SAMP_LIGHT / *MAT_187L.*MAT_SAMP_LIGHT (*MAT_187L) to choose type of effective strain rate calculation.*MAT_SAMP_LIGHT (*MAT_187L) to better deal with highly nonlinear inputs, meaning sharp curves and logarithmic strain rate interpolation.*MAT_BOLT_BEAM / *MAT_208, added feature (see input parameter HOLSHR) to enable modeling of a shear deformation mode in which the bolt shank tears through the plates. After the clearance gap has been closed, further shear deformation would enlarge the hole in the plates rather than deform the bolt itself. The forcedeformation relation of this mechanism is still governed by LCSHR, but the deformation (i.e. enlargement of the hole) is tracked separately in each of four principal directions. Thus, enlargement of the hole in the positive Y -direction has no effect on the position of the edge of the hole in the negative Y -direction.*MAT_224_GYS with the one in *MAT_224, e.g., NUMINT=-200 should prevent element erosion at all.*MAT_240 to include connection partner properties. With FUNCTIONS, the following parameters can/must be defined as *DEFINE_FUNCTION IDs: EMOD, GMOD, G1C_0, T0, FG1, G2C_0, S0, and FG2 (and GMOD3, G3C_0, R0, and FG3 if keyword option 3MODES is also used). The arguments of these functions include thicknesses and maximum yield stresses of both weld partners (from tied contact) as well as the strain rate and element area. A similar option is already available for *MAT_169.*MAT_240. Parameter INICRT can now be given as negative value, invoking a flexible exponent in the yield and damage initiation criteria.*MAT_240: exponential moving average of the effective plastic strain rate (RFILTP > 0) or of the effective total strain rate (RFILTP < 0).*MAT_240, namely number 16 is FG1 and number 17 is FG2.*MAT_251 in implicit analysis for solid and shell elements.*MAT_251 as a potential weld partner material for PROPRUL = 2 / 3 of *DEFINE_CONNECTION_PROPERTIES. Therefore use values from constant inputs PHASE1 and PHASE2 to find the associated curve in a *DEFINE_TABLE_3D or *DEFINE_TABLE_4D.*MAT_254, moved extended Koistinen-Marburger model to a new transformation law 9. Transformation 1 is the standard Koistinen Marburger again to ensure backward compatibility.*MAT_261 (*MAT_LAMINATED_FRACTURE_DAIMLER_PINHO):
*MAT_262 (*MAT_LAMINATED_FRACTURE_DAIMLER_CAMANHO), add a warning message in case SL > YC, as this may lead to unexpected behavior especially in fiber compression. IF SL is big enough, the approximation for ‘etal’ (longitudinal friction coefficient) will result in values >1. That will then in turn lead to phifc always being 0 (eq.(21) in CamI-paper).*MAT_280: FMOD = 10, 11, and 12 use a slightly modified condition for compressive failure.*MAT_GLASS (*MAT_280). It is invoked by FT < 0, with curve ID |FT| specifying tensile strength as a function of strain rate.*MAT_280 to allow FT < 0 to be used with FTSCL > 0. In that case, |FT| defines a curve for tensile strength as a function of strain rate and FTSCL scales the strength values from that curve as long as the material is intact. If cracked, neighbors get non-scaled values from that curve.*MAT_280 with ENGCRT and RADCRT. If used together with *MAT_ADD_EROSION (where these criteria are also available), *MAT_280 only ‘wins’ if those parameters are nonzero.*CONTROL_MPP_REBALANCE to control features associated with dynamic load balancing.*CONTROL_STRUCTURED.*DATABASE_BINARY_OPTION works for keyword options INTFOR, D3DUMP, and DEMFOR.*DATABASE_RECOVER_NODE, added support for solid element formulations – 1 and -2 for the elemental extrapolation method (METHOD = 1).*LOAD_ERODING_PART_SET to nodfor output.*DATABASE_CROSS_SECTION_PLANE. When RADIUS < 0.0, the radius is the absolute value of RADIUS, and XCT and XCH are node IDs. The center of the circular cut plane is given by the location of XCT. The normal vector of the plane is the vector pointing from XCT to XCH.*DATABASE_EXTENT_BINARY.*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 ).*CONSTRAINED_INTERPOLATION_SPOTWELD (SPR3) to binout / swforc output.*DEFINE_CURVE_TITLE and *DEFINE_COORDINATE_…_TITLE to d3hsp .*VENDOR_LICENSE scheme, both with the FlexLM server and with the lstc_server licensing system. Documentation for this will be released soon.*DEFINE_HAZ_TAILOR_WELDED_BLANK to allow load curves, ISW and IFW, from Card 2 of *DEFINE_HAZ_PROPERTIES to either increase monotonically (IMONFLAG = 0) or arbitrarily (IMONFLAG = 1).*DEFINE_TABLE_COMPACT. It allows a more compact notation of 2D-, 3D-, 4D-, … tables. Everything takes place in the keyword reader, where a *DEFINE_TABLE_{X }D is created internally.*DEFINE_MATERIAL_HISTORIES. This will allow you to specify names for history variables.*DEFINE_HEX_SPOTWELD_ASSEMBLY). Up to now, this error or warning (depending on SPOTSTP from *CONTROL_CONTACT) was only issued for single spot weld solids.*USER_LOADING interface to support thick shell elements.*FREQUENCY_DOOMAIN_FRF, *FREQUENCY_DOMAIN_SSD, and *FREQUENCY_DOMAIN_RANDOM_VIBRATION). The results before this fix were incorrect when pressure loading a curved surface.*FREQUENCY_DOMAIN_RESPONSE_SPECTRUM Fixed a bug in combining ground motion in response spectrum analysis (*FREQUENCY_DOMAIN_RESPONSE_SPECTRUM) when the load is defined by time history (LCTYP = 10, 11, or 12).*FREQUENCY_DOMAIN_FRF Fixed bug for frf (*FREQUENCY_DOMAIN_FRF) with response in all x, y and z directions. Before this fix, the output file names are incorrect.*FREQUENCY_DOMAIN_SSD) in MPP which only shows up on IBM and Windows.*RIGIDWALL_GEOMETRIC_FLAT_MOTION_ID_DISPLAY with adaptivity.*USER_INTERFACE_FRICTION, support h-adaptivity (MPP), meaning if *CONTROL_ADAPTIVE with ADPTYP = 1 or 2 is used, then history variables uhnew of subroutine usrfrc are transferred correctly from the old to new mesh in each adaptive refinement step.*CONTROL_ADAPTIVE). The issue could lead to a segmentation fault.*AIRBAG_SHELL_REFERENCE_GEOMETRY. This bug could cause difficulties in the element normal calculations.*ALE_STRUCTURED_MESH_TRIM) for the case of trimming both sides of a structure by specifying a part set. Previously, only one side of the structure was being trimmed, even if both sides were requested.*ALE_STRUCTURED_MESH_TRIM) that may lead to NaN.*MAT_HIGH_EXPLOSIVE_BURN / *MAT_008. The bug caused an incorrect detonation time.*ALE_STRUCTURED_FSI with Lagrange structure containing triangles. It might cause erroneous results, mostly leakage. It sometimes caused crashes.*ALE_STRUCTURED_FSI. The bugs were in certain cases causing the time step to drop or crashes.*BOUNDARY_NON_REFLECTING_2D. The implementation was not accounting for all the 2D segments.*ALE_ESSENTIAL_BOUNDARY did not work properly in R12.0. It is fixed in this version.*BOUNDARY_PRESCRIBED_MOTION_SET_BOX. The time was not being offset correctly when TOFFSET = 1.*BOUNDARY_PRESCRIBED_MOTION_SET_BOX. It occurred when the node set is part of a rigid body.*BOUNDARY_PRESCRIBED_MOTION.*BOUNDARY_PRESCRIBED_MOTION not work with *DEFINE_BOX_LOCAL.*BOUNDARY_PRESCRIBED_MOTION‘s that supersede *BOUNDARY_SPC‘s.*BOUNDARY_MCOL:
*BOUNDARY_NON_REFLECTING now uses the correct velocity in implicit so that explicit results are matched.*CONSTRAINED_NODAL_RIGID_BODY that occurs when using *INITIAL_VELOCITY_GENERATION with NX = -999. The rotational axis was incorrect.*CONSTRAINED_JOINT_GEARS) for implicit. We were not accounting for the helix angle.*CONSTRAINED_INTERPOLATION. It has been broken for some time before the release of R12.0.0.*CONSTRAINED_INTERPOLATION_SPOTWELD (SPR3) with huge number of associated nodes. A memory issue was likely to occur if more than 200 nodes were found in the domain of influence of upper or lower shell element part.*CONSTRAINED_SPR2 and *CONSTRAINED_INTERPOLATION_SPOTWELD (SPR3): If nodes in the domain of influence get completely free due to parent element failure, these nodes should not get any forces or moments distributed anymore.*CONSTRAINED_BEAM_IN_SOLID:
*CONTROL_SHELL).*CONTACT). The values were too low.*CONTACT_ATUOMATIC_SINGLE_SURFACE and *CONTACT_AUTOMATIC_*_TO_SURFACE. It covers both groupable and non groupable.*CONTACT_RIGID_SURFACE initialization which has in some cases been using the wrong coordinates in the single precision version of the code since 08/10/2016.*CONTACT_DRAWBEAD) behavior when the contact is off the edge of a thin triangle. It was producing NaNs during the friction calculation.*CONTACT_DRAWBEAD). It had an overly restrictive bucketsort, so it could miss some of the ramp up force a drawbead should create.*PART_CONTACT which could result in some nodes not being tied that should be.*CONTACT) when used with PSTIFM > 0 and selective mass scaling. Segment masses were being scaled incorrectly causing excessive stiffness.*ELEMENT_SHELL_COMPOSITE is in the input deck. Not all parts were being included.*TERMINATION_CONTACT to make it work when used with THRES > 0 or DOF = 2 or 3. The threshold was being compared with the maximum nodal contact force instead of the resultant. For DOF = 2 and DOF = 3, the z and y force, respectively, were being checked instead of y and z. Now resultant force is used and all options are working.*CONTACT_AUTOMATIC_GENERAL, SURFACE_TO_SURFACE, NODES_TO_SURFACE, and ONE_WAY_SURFACE_TO_SURFACE. Forces were sometimes being assigned to the wrong interfaces or to nonexistent interfaces, in which case a memory error could have cause as segmentation fault.*CONTACT_2D_FORCE_TRANDUDCER) that have both a slave and master surface defined.*CONTACT_AUTOMATIC_GENERAL in SMP missed the contact between beams generated with *ELEMENT_BEAM_SOURCE.*CONTACT_TIED_SURFACE_TO_SURFACE after dynamic relaxation when death time is set to -9999.*DEFINE_BOX_LOCAL abd the master box (MBOXID) is defined with *DEFINE_BOX as specified on Card 1 of *CONTACT_OPTION.*CONTACT_SLIDING_ONLY_PENALTY and running with SMP and NCPU =|| = 1.*CONTACT_SURFACE_TO_SURFACE. This contact is two-way, but one-way contact was being performed instead.*CONTACT_TIED_SURFACE_TO_SURFACE when using implicit in MPP. Before the resultant forces were being reported as zeroes for this contact type.*CONTACT_TIEBREAK_NODES_TO_SURFACE_ID is not supported for implicit computations.*CONTACT_AUTOMATIC_BEAMS_TO_SURFACE. Prior to this fix, unreasonable penetration may have been computed that would have led to wrong contact forces.*CONTACT_ADD_WEAR for MPP.*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, now resolved.THERMAL_FRICTION option of *CONTACT). Bug could result in neglecting heat transfer and throw a warning even if curve was correctly defined.*CONTROL_OUTPUT and only wrote the corner nodes of the elements to d3plot for the connectivity. Now it accepts the setting of TET10S8, so the full connectivity can be written to the d3eigv database.*CONTROL_STAGED_CONSTRUCTION). The scaling down had been erroneously omitted for the hourglass force pertaining to the element types and hourglass setting mentioned above.*CONTROL_PORE_FLUID):
*CONTROL_PORE_FLUID) together with suction limits (see SUCLIM on *BOUNDARY_PORE_FLUID). If a node on the tied contact reached the suction limit, the suction limit pressure then spread along the tied contact in a non-physical way.*CONTROL_PORE_FLUID) together with time-dependent consolidation (ATYPE = 3). Pore pressures could become unstable at nodes shared across MPP processors where the node belongs to more than one tied contact.*BOUNDARY_PORE_FLUID to set the analysis type to time-dependent consolidation (ATYPE = 3) for at least 100 time steps before switching to ATYPE = 4.*ELEMENT_BEAM_THICKNESS. The local coordinate ID specified by PARM3 was ignored by the discrete beam material models *MAT_069, *MAT_070, *MAT_071, *MAT_074, *MAT_094, *MAT_197, *MAT_205 and *MAT_208. The global coordinate system was used instead.*ELEMENT_SEATBELT:
*SECTION_SHELL).*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.*ELEMENT_SHELL. A redundant conversion from degrees to radians was changing the material direction.*ELEMENT_SHELL_COMPOSITE. This affects MPP only.*ELEMENT_SHELL_COMPOSITE and when *MAT_ADD_THERMAL_EXPANSION is present in the input.*ELEMENT_DIRECT_MATRIX_INPUT).*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.*EM_CONTROL).*EM_CONTROL) slower.*FATIGUE_MULTIAXIAL). The fatigue damage ratio was incorrect (multiaxial correction is not enabled) before this fix.*ICFD_BOUNDARY_FSI_EXCLUDE. Solid parts were not being excluded from FSI, only shells.*MESH_BL). The surface mesh was deforming when boundary layer nodes were clamped.*MESH_EMBEDSHELL).*ICFD_MODEL_POROUS). It was using the permeability instead of the porosity.*MAT_ADD_THERMAL_EXPANSION for implicit.*CONTROL_IMPLICIT_BUCKLING. It was broken in R11.*CONTROL_IMPLICIT_SOLUTION). Prior to the bug fix, the statistics were additive, meaning the statistics for the new time step were added to the ones from the previous time step. Now the statistics are reset at the beginning of the new time step, so they only reflect the current time step.*CONTROL_IMPLICIT_EIGENVALUE,
*CONTROL_IMPLICIT_MODES in MPP when nodes involved were shared across processors.*CONTROL_IMPLICIT_DYNAMICS and *CONTROL_IMPLICIT_AUTO with the options DYN or SPR enabled. Before if the curves IDs did not match the input order (such as an ID of 97125 for the curve input), then LS-DYNA would given an error. This is not how curves are normally processed. Now they are being processed like all other curves in LS-DYNA where the curve ID and order of input are independent.*CONTROL_IMPLICIT_MODES) in MPP for the case where a constraint mode is also a shared mode.*CONTROL_IMPLICIT_GENERAL, adjust logic of when the last intermittent eigenvalue computation is performed during an explicit simulation if the last one is performed near the termination time. For instance, if the termination time is 10.00 and the explicit time step is 0.01, the the explicit calculation terminates between 9.99 and 10.00. Therefore, if you wanted an eigenvalue computation at 10.00, it would not happen because the calculation would have already terminated. This change causes the intermittent eigenvalue computation to occur before the termination time.*INTERFACE) may be shared by multiple processes. We fixed a bug in implicit where we were double counting the interface forces for the shared nodes.*CONTROL_IMPLICIT_MODAL_DYNAMICS) was still trying to read the nonexistent data for rotational dofs and became lost. This revision corrects that oversight.*CONTROL_ACCURACY. This applies to implicit only. There was numerical noise for no deformation.*CONTROL_RIGID. The outputs were wrong.*LOAD_SEGMENT_FILE and *LOAD_THERMAL_D3PLOT in implicit when time step is cut.*INITIAL_STRESS_SECTION, fixed segmentation fault for ISTIFF < 0 when running SMP on many cores.*INITIAL_VELOCITY_GENERATION with PHASE = 1 and NX = -999. This affects MPP only.*INITIAL_VELOCITY_GENERATION yields incorrect velocity for *CONSTRAINED_NODAL_RIGID_BODY_INERTIA when ICID is not zero.*INITIAL_HISTORY_NODE(_SET), single precision executables for Windows or IBM linux may have terminated unexpectedly.*INITIAL_STRESS_SHELL buffers. It is needed especially when high order shells are used.*INITIAL_AXIAL_FORCE_BEAM with beam element 1 (Hughes-Liu beam). Available since R12.0, this feature could have lead to various side effects, such as severe beam deformations.*CONTROL_PARALLEL) and contact is involved for IGA elements.*LOAD_MOVING_PRESSURE.*MAT_ADD_DAMAGE_DIEM. It could have resulted in a memory clobber for certain inputs.*MAT_ADD_EROSION and *MAT_ADD_DAMAGE_DIEM/GISSMO. Values from *MAT_ADD_EROSION‘s optional cards 4 and 5 could have been ignored.*MAT_ADD_DAMAGE_GISSMO: LCSRS being a table and using logarithmic strain rates (first value negative) did not work correctly before.*MAT_ADD_DAMAGE_GISSMO with HISVN = -(ND+19). It did not work correctly for shell elements.*MAT_ADD_DAMAGE_GISSMO with LCSDG < 0 to avoid incorrect element deletion (only SMP/Hybrid with ncpu > 1).*MAT_ADD_DAMAGE_GISSMO with LCSDG being a 3D table (*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 *MAT_ADD_GENERALIZED_DAMAGE in full deck restart. Element size for LCREGD and FADEXP < 0 was re-computed after restart, but it should not.*MAT_ADD_GENERALIZED_DAMAGE with long=s. The long input format was not working with this keyword before.*MAT_ADD_THERMAL_EXPANSION with solid element formulation 13. Incorrect elastic strain was used resulting in incorrect stress.*MAT_NONLOCAL when used with eroding elements. The smoothing calculation was using uninitiialized data due to a failure to set initial values for the eroded elements prior to smoothing. This error could be subtle or cause an error termination.*MAT_ELASTIC_FLUID output to d3hsp. The bulk modulus, k, was described incorrectly as ‘e’ in d3hsp .*MAT_ORTHOTROPIC_ELASTIC. The response was not linear.*MAT_PIECEWISE_LINEAR_PLASTICITY. Sometimes (rarely) the interpolation between load curves failed.*MAT_024 (shells) with TDEL > 0 and GISSMO. History variables conflict could lead to wrong results.*MAT_026, *MAT_126, or *MAT_201 with thick shell elements when AOPT = 2 or 3. The material directions were being miscalculated causing wrong stress.*MAT_030:
*DEFINE_MATERIAL_HISTORIES is used with *MAT_FABRIC. History data was being clobbered.*MAT_USER_DEFINED_MATERIAL_MODELS / *MAT_041-050. Before this fix user some defined features ( usrshl, usrsld, usrfrc, usrtie, usrtbrk, usr_nunonl, mortar_*, rebar_* ) using *MODULE did not work correctly.*DEFINE_MATERIAL_HISTORIES is used with *MAT_USER_DEFINED_MATERIAL_MODELS. The bug led to an inconsistent count of history variables.*MAT_053, extrapolate load curve for von Mises stress as a function of negative volumetric strain using the first 2 points of the load curve if the volumetric strain is less than the first point of the curve. This fix prevents an error termination.*MAT_054 (*MAT_ENHANCED_COMPOSITE_DAMAGE):
*AIRBAG definitions were present in the input.*MAT_058 (*MAT_LAMINATED_COMPOSITE_FABRIC):
*MAT_076 with temperature shift. Initial temperature in an uncoupled analysis is not always zero (such as that coming from *LOAD_THERMAL).*MAT_077 / *MAT_181 with VFLAG = 1 and *INCLUDE_TRANSFORM. Normalized shear relaxation moduli Gi (unitless) should not be converted.*MAT_077, fixed bug for shell elements that could cause instability.*MAT_FU_CHANG_FOAM / *MAT_083 on initial time step size.*MAT_WINFRITH_CONCRETE_REINFORCEMENT / *MAT_084_REINF did not work in MPP when the reinforced elements were defined with EID1 and EID2 (Card 1a). It worked correctly when defined by the PID method (Card 1b).*MAT_100 (*MAT_SPOTWELD) for OPT = 0/-1 when using *DEFINE_CURVEs for defining force / moment resultants as a function of the effective strain rate, by setting NRR, NRS, NRT, MRR, MSS, MTT as negative values. Added a more meaningful error message in case load curve is not defined.*MAT_SPOTWELD_DAMAGE-FAILURE with beams, fixed bug that could have led to earlier damage evolution than anticipated.*MAT_106 for shell elements. In a structural-only analysis, the ‘old’ temperature should be set to the initial temperature, such as from *LOAD_THERMAL_LOAD_CURVE.*MAT_110 / *MAT_JOHNSON_HOLMQUIST_CERAMICS when you iput phel = 0 in the data cards.*MAT_110 / *MAT_JOHNSON_HOLMQUIST_CERAMICS at t = 0.0.*MAT_122_3D, updated plasticity algorithm due to implicit problems.*MAT_122_3D and also add missing d3hsp echo of material parameters.*MAT_MODIFIED_HONEYCOMB / *MAT_126 with solid element types 0 and 9 to work with *CONSTRAINED_TIED_NODES_FAILURE.*MAT_126 and *MAT_ADD_THERMAL_EXPANSION. It only worked for ELFORM = 0, but now also works for others.*MAT_136. Bug resulted in wrong strain output for post-processing.*MAT_142 from ‘anisotropic crushable foam’ to ‘transversly isotropic crushable foam’ in d3hsp to be consistent with the Users’ Material Manual.*MAT_155 will now call subroutine matusr_24 for FS < 0.*MAT_157 (*MAT_ANISOTROPIC_ELASTIC_PLASTIC):
*MAT_SIMPLIFIED_RUBBER/FOAM / *MAT_181 for solid formulation 13, tetrahedron with nodal pressure.*MAT_186: variable TSLC2 was not working in free format or with *PARAMETER.*MAT_187 to increase accuracy of results if LCEMOD > 0.*MAT_196 to avoid error termination.*MAT_208 (*MAT_BOLT_BEAM) bug affected R12.0 only. When AXSHFL = 0 and LCSHR was a *DEFINE_TABLE, the shear force did not correctly follow the relation with displacement defined in the *DEFINE_TABLE. It worked correctly when LCSHR was a *DEFINE_CURVE, and it worked correctly in R11 with both *DEFINE_CURVE and *DEFINE_TABLE.*MAT_224_GYS work correctly, which was not the case before.*MAT_PIECEWISE_LINEAR_PLASTIC_THERMAL (*MAT_255):
*MAT_258: repair visco-plasticity, use correct thickness of shell type 16, and eliminate potential dependence on ncpu.*MAT_264 / *MAT_TABULATED_JOHNSON_COOK_ORTHO_PLASTICITY which gives incorrect results when running with adaptivity.*MAT_BERGSTROM_BOYCE_RUBBER.*MAT_BERGSTROM_BOYCE_RUBBER according to the rules of IRATE on *CONTROL_IMPLICIT_DYNAMICS.*MAT_274, fixed bug where shell thickness strain was not updated properly.*MAT_278 / *MAT_CF_MICROMECHANICS.*PART_COMPOSITE or *INTEGRATION_SHELL was not used to define an integration rule.*TERMINATION_BODY. It was never properly implemented in MPP.*DEFINE_CURVE_FUNCTION. Specifically, if the built-in function PIDCTL was used, along with MPP predecomposition, the behavior was unpredictable due to a memory clobber.*DATABASE_CURVOUT, fix thermal-only solutions to write the curvout file in LSDA format. The ASCII version was working, but if the LSDA format was requested, then this data was not being written at all.*SECTION_SHELL, ELFORM=18), the moments are now correctly written to the cross section output generated by *DATABASE_SECFORC. Before, incorrect values were reported.*DATABASE_RECOVER_NODE:
*DATABASE_CROSS_SECTION, fixed a bug for cross sections that contain 2D seat belts.*DATABASE_EXTENT_BINARY is negative, reset it to positive and issue warning message, KEY+1396.*DATABASE_EXTENT_D3PART and STRFLG = 0 in *DATABASE_EXTENT_BINARY.*DAMPING_RELATIVE.*DATABASE_SSSTAT_MASS_PROPERTIES).*INITIAL_VELOCITY_GENERATION*DATABASE_NCFORC_FILTER in MPP.*DATABASE_RCFORC_MOMENT). This bug only affects post-processing, not the analysis itself.*DATABASE_EXTENT_BINARY), the d3thdt output file was corrupted. Results could still be obtained via the ascii and binout files, only the d3thdt file was affected.*DATABASE_CROSS_SECTION and one or more of the elements eroded.*DATABASE_EXTENT_BINARY and *DEFINE_MATERIAL_HISTORIES. Prior to this fix, the output stresses might not have been properly transformed into the material frame. This only affects the post-processing and not the analysis itself.*DEFINE_MATERIAL_HISTORIES. An incorrect number of variables and associated labels were reported.*DEFINE_MATERIAL_HISTORIES. The values were incorrect.*CONTROL_OUTPUT.*DATABASE_DCFAIL.*DATABASE_DISBOUT output files.*DATABASE_BINARY_INTFOR and full deck restart:
*INTERFACE_LINKING file on processor 0 in MPP. The other processors were trying to open the file with a nonsense name, causing a non-obvious ‘file open’ error message.*DEFINE_MATERIAL_HISTORIES and *DEFINE_FUNCTION_TABULATED for newformat=long option*CONTACT_TIEBREAK_SURFACE_TO_SURFACE. This affects SMP only.*SENSOR_DEFINE_FORCE with TYPE = JOINTSTIF. It didn’t trace the force correctly when more than one core is used.*SENSOR_DEFINE_ELEMENT with TYPE = STRAIN. It occurred when STRFLG of *DATABASE_EXTENT_BINARY is 0.*SENSOR_DEFINE_FORCE with CRD \u0338 = 0 and long=s. Error in structured read happened (STR+211) with MPP and np > 1.*DATABASE_EXTENT_BINARY now honors IEVERP in thermal only calculations.*CONSTRAINED_NODAL_RIGID_BODY cards in a thermal-only analysis, the d3plot files were corrupted and could not be read by LS-PrePost. The solution found by the solver was correct.*CONSTRAINED_ADAPTIVE defined on solid/beam nodes when thermal thick shells are also present in the model.*CONSTRAINED_ADAPTIVITY and *CONSTRAINED_SPOTWELD are excluded from the *BOUNDARY_TEMPERATURE constraint.*BOUNDARY_FLUX condition.*BOUNDARY_THERMAL_WELD.*SECTION_SHELL_XFEM) now works for plastic materials.*BOUNDARY_PRECRACK. Gave warning message when precrack points are changed to avoid the precrack passing through nodal points.*SET_SEGMENT_GENERAL, which could have resulted in incorrect values being applied for the DA1-DA4 parameters, or possibly a segmentation fault during input processing.*INTERFACE_LINKING that could lead to segmentation faults.*DEFINE_FUNCTION in structured input. Also made sure they work properly with predecomposition.*PARAMETER$ PRMR1 VAL1 PRMR2 VAL2RuPLscXXa1 &R1uPllSARuPLscXXt1 0.0*PARAMETER if long format is used.*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.*PARAMETER that occurs if there is no space between a symbol and a value like in Listing 1.*DEFINE_CURVE_FUNCTION:
*MODULE_PATH did not work correctly for user defined features.*DEFINE_PRESSURE_TUBE:
*DEFINE_MATERIAL_HISTORIES to be used with Rayleigh damping. The history data and damping variables were conflicting, causing early failure of elements.*NODE_TRANSFORM when used with spot weld assembly generation. Spot weld elements were getting distorted.*DEFINE_TABLE_2D/3D are contained in a different include files after the include file that contains the *DEFINE_TABLE_2D/3D statement.*INTEGRATION_BEAM with *INCLUDE_TRANSFORM.*DEFORMABLE_TO_RIGID_AUTOMATIC when CODE = 4 and OFFSET > 0.0. Offset was not working correctly.*RIGIDWALL_PLANAR.