Home

LS-DYNA R7.1.1 Release notes

1. recommended for ALE multi material use e Upgraded MAT_172 MAT_CONCRETE_EC2 e Corrections to DEGRAD option e Fix bug affecting shells not beams When an element has two cracks and one crack is under compression while the other is still widening the element could erroneously shrink in the compression direction Another symptom of the same bug would be that if two cracks are opening at the same time in the same element one or both might not follow the expected stress crack opening curve Ove Arup and Partners Limited June 16 June 2014 Page 13 LS DYNA Software Release Notes Version 971 R7 1 1 e Concrete and reinforcement types 7 and 8 have been added to reflect changes to Eurocode 2 e Extra history variables for reinforcement stress and strain are now output as zero for zero fraction reinforcement directions e Added RCDC model for solid MAT_082 e Added Feng s failure model to solid MAT_021 e Added MAT 027 for beams e Added DEFINE_HAZ PROPERTIES and DEFINE_TAILOR_WELDED_BLANK for modifying material behaviour near a spot weld e Added fourth rate form to viscoplastic Johnson Cook model MAT_015 e Added option to MAT_224 to not delete the element if NUMINT 200 e New damage initiation option 3 in multi fold damage criteria in MAT ADD EROSION Very similar to option 2 but insensitive to pressure e Rotational resistance in MAT_034 MAT_FABRIC Optionally the user may specify the stiffness yield a
2. CONSTRAINED_NODAL_RIGID_BODY nodes are shared by several processors with mass scaling on the added mass is not summed up across processors This results in an instability of the NRB MPP only 3 9 CONTROL e ALE REFINE has been replaced and expanded upon by the CONTROL_REFINE family of commands These commands invoke local mesh refinement of shells solids and ALE elements based on various criteria Shells or solids in a region selected for refinement parent element are replaced by 4 shells or 8 solids respectively CONTROL_REFINE_SHELL applies to shells CONTROL_REFINE_SOLID applies to solids CONTROL_REFINE_ALE and CONTROL_REFINE_ALE2D applies to ALE elements Ove Arup and Partners Limited June 16 June 2014 Page 7 LS DYNA Software Release Notes Version 971 R7 1 1 Each keyword has up to 3 lines of input If only the Ist card is defined the refinement occurs during the initialization The 2nd card defines a criterion CRITRE to automatically refine the elements during the run If the 3rd card is defined the refinement can be reversed based on a criterion CRITM All commands are implemented for MPP CONTROL_REFINE_MPP_DISTRIBUTION distributes the elements required by the refinement across the MPP processes Fixed bug in CONTROL_STAGED_CONSTRUCTION whereby staged construction was not working with the combination of Hughes Liu beams and CONTROL_ACCURACY 3 10 DATABASE Eliminate automatic writing of a d3plot
3. MAT_076 related to Laminated shell theory Extra prony series used in implicit shells Ove Arup and Partners Limited June 16 June 2014 Page 44 LS DYNA Software Release Notes Version 971 R7 1 1 Temperature scaling e Fix bad tangent matrix in MAT_027 for 2D elements e Fix incorrect stresses strains when using MAT_PIECEWISE_LINEAR_PLASTIC_THERMAL MAT_255 with LOAD_THERMAL_ OPTION e Fix output to elout for shell elements with MAT_FABRIC when EOCS 2 for CONTROL_OUTPUT Stresses strains were still output in local or lamina coordinate system instead of global coordinate system e Fix convergence error when using MAT_ORTHOTROPIC_SIMPLIFIED_DAMAGE MAT_221 for implicit e Fix divide by zero for MAT_HILL_FOAM MAT_177 during the stress tensor eigenvector iterations which resulted in elements disappearing in d3plot e Fix MAT_SIMPLIFIED_JOHNSON_COOK_ORTHOTROPIC_DAMAGE MA T_99 shell elements which were not deleted despite reaching the rupture strain and number of intg points failed e Fix MAT_PLASTIC_NONLINEAR_KINEMATIC MAT_165 for implicit analysis e Fix incorrect large time step when using MAT_MODIFIED_PIECEWISE_LINEAR_PLASTICITY without any failure criteria FAIL TDEL EPSMAJ for solid elements e Fix incorrect viscous force to elout in binout for FMAT_MUSCLE MAT_ 156 e Fix nonzero z strain when using MAT_USER_DEFINED_MATERIAL_MODELS and shell element type 13 with nip 4 e Fix zero stress strains when
4. PARAMETER_TYPE according to how that ID is used e Allow use of load curve to specify damping as a function of time in DAMPING_RELATIVE e Adda segment based SOFT 2 contact option to include the overlap area in the contact stiffness calculation This is good for improving the friction calculation and possibly for implicit convergenece The option is turned on by setting FNLSCL gt 0 and DNLSCL 0 As DNLSCL 0 the contact stiffness is not nonlinear This new option is also useful when used with another improvement that was made to the FS 2 friction coefficient by table lookup option in segment based contact When the above mentioned FNLSCL gt 0 option is used the FS 2 option is now very accurate e Add a new RCDC damage option MAT_PLASTICITY_WITH_DAMAGE_ORTHO_RCDC1980 which is consistent with the WILKINS paper It uses the principal values of stress deviators and a different expression for the A_d term e Adda TIETYP option to CONTACT 2D AUTOMATIC By default the tied contact automatically uses constraint equations when possible for 2D tied contact If a conflict is detected with other constraints or to avoid 2 way constraints penalty type ties are used when constraints are not possible The TIETYP option when set to 1 causes all ties to use the penalty method This Ove Arup and Partners Limited June 16 June 2014 Page 23 LS DYNA Software Release Notes Version 971 R7 1 1 is useful if in spite of the code s best effo
5. 1 1 6 R7 1 1 EXECUTABLES MPP Executables SMP Executables Ove Arup and Partners Limited June 16 June 2014 Release Notes Version 971 R7 1 1 34 35 38 41 42 43 43 44 51 52 52 53 53 55 56 58 59 60 61 62 64 65 66 67 LS DYNA Software Release Notes Version 971 R7 1 1 1 Introduction This document lists the update and release notes for version R7 1 1 of LS DYNA As many users will be coming to version R7 of LS DYNA directly from LS971 R6 1 1 this document includes the following three sets of updates e Updates in version R7 0 0 from Version LS971 R6 1 1 Updates in version R7 1 0 from Version R7 0 0 Updates in version R7 1 1 from Version R7 1 0 Please note the LS971 naming has been removed as from version R7 0 0 onwards This is reflected in the licenses file needed to run an R7 version of LS DYNA The license file must contain the following lines LS DYNA 12312013 OPTION REVISION 7 MPPDYNA 12312013 OPTION REVISION 7 Ove Arup and Partners Limited June 16 June 2014 Page 1 LS DYNA Software Release Notes Version 971 R7 1 1 2 Release notes for LS DYNA version R7 New solvers and mesher in the double precision version of LS DYNA R7 e Electromagnetic EM solver Compressible CFD CESE solver Incompressible CFD ICFD e Volume mesher supporting the ICFD and CESE solvers Ove Arup and Partners Limited June 16 June 2014 Page 2 LS DYNA Software Release Notes Version 971
6. 16 June 2014 Page 54 LS DYNA Software Release Notes Version 971 R7 1 1 Correct the computation of superelements using CONTROL_IMPLICIT_MODES Enhance implementation of PART_MODES for parts made up of beam elements Rectify the node count due to differences with beam reference nodes Correct PART_MODES implementation for reading mode shapes from d3mode database Fix convergence problem when using idrflg 5 for CONTROL_DYNAMIC_RELAXATION and MAT_MODIFIED_PIECEWISE_LINEAR_PLASTICITY MAT123 Add support for IGAP gt 2 in implicit segment based SOFT 2 contact Fix a serious error in MPP implicit 2D automatic contact The node numbers passed to the stiffness routine were completely wrong and problems would error terminate Fix the section force calculation DATABASE_SECFORC to work with the implicit formulation 3 solid element Allow birth and death time in CONTROL IMPLICIT DYNAMIC to be effective when used together with CONTROL_IMPLICIT_FORMING Correct the computation of resultant forces for rigid bodies during implicit mechanics computation 4DATABASE BNDOUT BOUNDARY_PRESCRIBED_MOTION_RIGID For nonlinear implicit solver 12 convergence is prevented if boundary conditions BOUNDARY_PRESCRIBED _ are not properly fulfilled For this the message Convergence prevented due to unfulfilled bc is issued and the solver keeps iterating See CONTROL_IMPLICIT_SOLUTION Allow CONTROL IMPLIICT FORMING to work
7. 39 LS DYNA Software Release Notes Version 971 R7 1 1 outputting the wrong strain to d3plot and elout files This fix otherwise has no effect on results e Fix parallel consistency for triangular shell formulation 24 Both implicit and explicit jobs were failing due to incorrect forces calculated for mid side nodes e Fix the secforc data DATABASE_SECFORC for tet formulation 4 The sign was wrong on the forces and moments e Fix the mass calculation in the initialization phase to determine the mass of thick shell rigid bodies so that the total mass reported is correct This likely had no effect on results e Fix parallel consistency for warped beam formulation 11 It was giving incorrect results for all but small problems and could cause convergence failure of implicit problems e Fix the option to use a negative connection ID on DEFINE CONNECTION PROPERTIES with solid weld assemblies The failure function was evaluated but the weld did not fail as it should e Fix possibility of a segmentation fault in thick shell formulation 3 when used with isotropic materials in an implicit run e Improve cohesive solid element formulation 20 e Correct orientation check for cohesive elements Solids elform 19 20 MPP only e Improve implicit capability of pulley elements ELEMENT_BEAM_PULLEY e Improve characteristic length Le calculation for 10 noded tet elements ELFORM 16 which is used in time step calculation For flat
8. Arup and Partners Limited June 16 June 2014 Page 4 LS DYNA Software Release Notes Version 971 R7 1 1 Loose or strong coupling for FSI Fluid structure interaction Exact boundary condition imposition for FSI problems Solver Applications Non exhaustive External aerodynamics for incompressible flows Internal aerodynamics for incompressible flows Sloshing Slamming and Wave impacts FSI problems Conjugate heat transfer problems 3 4 New Volume Mesh generator Keyword family MESH_ The keywords starting with MESH refer to and control the tools for the automatic volume mesh generator Mesh Generator Characteristics Automatic Robust Generic Tetrahedral elements for 3D Triangles in 2D Closed body fitted mesh surface mesh needs to be provided for volume generation Mesh Generator Main Features e Automatic remeshing to keep acceptable mesh quality for FSI problems ICFD only Adaptive meshing tools ICFD only Anisotropic boundary layer mesh Mesh element size control tools Remeshing tools for surface meshes to ensure mesh quality Mesh Generator Applications e Used by the Incompressible CFD solver ICFD e Used by the Compressible CFD solver CESE 3 5 AIRBAG e Add new parameter VNTOPT to AIRBAG_HYBRID that allows user more control on bag venting area calculation e Allow heat convection between environment and CPM bag FAIRBAG_PARTICLE bag Apply proper probability density function to part
9. IGNORE 0 it now reports the distance based on the initial configuration which is the only meaningful value It was reporting either 0 or the distance before the last iteration neither of which is meaningful Ove Arup and Partners Limited June 16 June 2014 Page 30 LS DYNA Software Release Notes Version 971 R7 1 1 Decrease memory used by keyword reader while processing drawbeads Change spatial extension to the segment used in MPP tied contact The MAXPAR parameter used by SMP has long been ignored with MPP using PARMAX from the optional _MPP card with a default value of 1 0005 For non tied contacts this PARMAX value is used as given but for tied contacts a value of max PARMAX 1 035 was used instead After due consideration about not wanting to change behaviour for old user decks the following input is now accepted if 0 0 lt PARMAX lt 0 5 use 1 PARMAX if 1 0 lt PARMAX lt 1 0004 use PARMAX else use max PARMAX 1 035 This allows any specific value gt 1 to be specified without risking breaking user inputs that had values of 1 0005 to 1 035 Also change groupable tied contact to honour this new scheme Implement splitting pinball contact Belytschko amp Yeh 1992 1993 This new contact option is invoked by setting SOFT 2 SBOPT 3 and DEPTH 45 A penetration check method based on LS PrePost version 4 0 is implemented for the new bilinear patch based contact SOFT 2 DEPTH 45 amp Q2TRI 0 The new me
10. INITIAL STRESS BEAM e Correct contact stiffness calculation for MAT 224 with E lt 0 curve and MAT_244 with E lt 0 table e Fixed MAT_VISCOELASTIC MAT_006 for type 13 tet elements e Switch off rate effects in MAT_083 for static implicit analyses or dynamic analyses with IRATE 1 or 2 This was already working for the viscous Ove Arup and Partners Limited June 16 June 2014 Page 49 LS DYNA Software Release Notes Version 971 R7 1 1 damping part DAMP but now is also working for the table TBID i e only the first load curve is used e Add shell element IDs as additional output to messag file for MAT_036 s warning plasticity algorithm did not converge e Enable varying temperatures over shell thickness for MAT_106 and MAT 188 e Fix MAT_224 when used with an equation of state EOS Internal energy and relative volume were not correctly initialized e Modification for MAT_244 TRIP effect should only be active for martensite phase 5 In addition TRIP effect can totally be disabled by setting input parameter CRSH 2 e Add plastic strain rate as history variable 9 to MAT_106 e Fix failure behaviour of MAT_PLASTICITY_POLYMER for multi integration point solid elements They are now eroded as soon as the first integration point fails e For MAT_ADD_EROSION damage model GISSMO parameter FADEXP allow values gt 0 Before only values gt 1 were possible e MAT ADD THERMAL EXPANSION is n
11. INITIAL_ VELOCITY GENERATION The local coordinate system ICID was not used to transform the velocity vectors and this has been fixed Fix incorrect initial kinetic energy in d3hsp and screen output when INITIAL VELOCITY_RIGID_BODY is used Ove Arup and Partners Limited June 16 June 2014 Page 42 LS DYNA Software Release Notes Version 971 R7 1 1 Initial velocities were not initialized during transient phase when idrflg 5 for CONTROL_DYNAMIC_RELAXATION despite presence of INITIAL_VELOCITY_GENERATION keyword MPP Only Fix INITIAL_VELOCITY_GENERATION not generating correct initial velocities when ID 0 Fix incorrect error when using scale factor SCALE in INITIAL_AXIAL_FORCE_ BEAM Improve output of INITIAL_STRESS_BEAM data to dynein via INTERFACE_SPRINGBACK_LSDYNA Now large format can be chosen history variables are written and local axes vectors are included Correct dynain output of INITIAL_STRAIN_TSHELL with np gt 1 MPP Only Fix spurious error message pertaining to INITIAL_VELOCITY_GENERATION when NSID 0 and STYP 1 Improve excessive input amp initialization time when using multiple INITIAL STRESS_SHELL with INITIAL_STRESS_SHELL_SET The same improvement applies to INITIAL_STRAIN_SHELL and INITIAL STRAIN_SHELL SET 4 12 LOAD Fix zero pressure load problem when using LOAD_SEGMENT with nodes 3 and 4 of the segment set to 0 for shell elements and running 2D analysis MPP Only Add war
12. Limited June 16 June 2014 Page 28 LS DYNA Software Release Notes Version 971 R7 1 1 e Trap input error if part id is 0 for BOUNDARY_PRESCRIBED_MOTIION_RIGID with error message KEY 927 as this could cause a crash later in the run e Change BOUNDARY_ACOUSTIC_COUPLING_MISMATCH to rank order opposing acoustic faces and structural segments by proximity thereby accelerating the preprocessing stage enhancing reliability and allowing some liberalization of the search parameters e Fix a bug where BOUNDARY_PRESCRIBED_ACCELEROMETER_RIGID failed when three or more rigid bodies were driven e Re active the error stop when DEFINE_CURVE_FUNCTION with simulation response is used for displacement boundary condition e Fix for boundary conditions for mid side nodal points generated by ELEMENT_SHELL_SHL4_TO_SHL8 e Fix incorrect motion of slave rigid body when using BOUNDARY_PRESCRIBED_MOTION_RIGID with vad 4 and the nodel node2 vector is orthogonal to the master rigid body motion e Fix bug wherein BOUNDARY_PRESCRIBED_ORIENTATION_VECTOR viewed small angle excursions as states of rest and produced unexpected motions 4 4 BLAST e Fix time of arrival of the Mach wave for BLAST 4 in LOAD BLAST ENHANCED e Correct the way forces are handled for triangular segments subjected to LOAD_BLAST_ENHANCED e Implement hemispherical geometry for particle blast DEFINE_PBLAST_GEOMETRY e Add explosive type for PARTICLE_BLAST e For
13. This combination is not recommended and generates a warning message but such models could previously generate unpredictable error messages e Added DEFINE_VECTOR_NODES A vector is defined using two node IDs e Add sense switch prof to output current timing profile to messag SMP file or mes MPP files Also for MPP only collect timing information from processor and output to prof out when sense switch prof is detected Ove Arup and Partners Limited June 16 June 2014 Page 21 LS DYNA Software Release Notes Version 971 R7 1 1 4 Release notes for LS DYNA version R7 1 0 4 1 NEW FEATURES GENERAL e Add MUTABLE option for PARAMETER so that parameter values can be redefined later in the input deck e Add new command line option Idir for setting a local working directory In MPP this has the same effect as setting the directory local pfile option and it overrides that option For SMP it indicates a directory where local working files should be placed e Add support for SMOOTH option in MPP groupable contact e Add new keyword card CONTROL_REQUIRE_REVISION to prevent the model from being run in old versions of LS DYNA e Add part set specification for dynamic relaxation with implicit using CONTROL_DYNAMIC_RELAXATION This feature specified with idrflg 6 allows implicit to be used for the dynamic relaxation phase for models involving parts being modelled with SPH and or ALE while excluding those
14. based contact for a tied constraint based contact in the input IPBACK on Card E of CONTACT e New option for CONTACT_ENTITY If variable SO is set to 2 then a constraint like option is used to compute the forces in the normal direction Friction is treated in the usual way e CONTACT_ENTITY allow friction coefficient to be given by a coefficient vs time load curve input lt 0 gt absolute value is the load curve ID Also if the friction coefficient gt 1 0 the node sticks with no sliding at all e Fix for MPP implementation of CONTACT_ORTHO_FRICTION e Minor tweak to the way both MPP and SMP handle nodes sliding off the ends of beams in CONTACT_GUIDED_CABLE e Fix MPP groupable tied CONTACT energy calculations The contact was working correctly but the energy calculation was incorrect since the changes of 171470 e Frictional energy output in sleout DATABASE _SLEOUT supported for CONTACT MORTAR e Tiebreak damage parameter output as contact gap in intfor file for CONTACT AUTOMATIC SURFACE TO SURFACE TIEBREAK MO RTAR OPTION 9 e Added MPP support for CONTACT_2D_AUTOMATIC_SINGLE_SURFACE and CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE 3 8 CONSTRAINED e Added keyword CONSTRAINED_MULTIPLE_GLOBAL for defining multi node constraints for imposing periodic boundary conditions e Enhancement for CONSTRAINED_INTERPOLATION SPOTWELD SPR3 calculation of bending moment is more accurate now e If
15. by CONSTRAINED_INTERPOLATION e Add a length unit to the tolerance used for the checking of non coincident nodes in CONSTRAINED_JOINTs excluding spherical joints The old tolerance was 1 e 3 e The new tolerance is 1 e 4 times the distance between nodes and 3 The error messages were changed to warnings since this change might otherwise cause existing models to stop running Ove Arup and Partners Limited June 16 June 2014 Page 33 LS DYNA Software Release Notes Version 971 R7 1 1 e Update single precision rigid bodies to improve the initialization of rigid body velocity so that imposed initial translational velocity does not create spurious rotational velocity e Add d3hsp output for CONSTRAINED_INTERPOLATION_SPOTWELD SPR3 and e CONSTRAINED_SPR2 Can be deactivated by setting NPOPT 1 on CONTROL_OUTPUT e Add warning message if dependent node of CONSTRAINED_INTERPOLATION is on a rigid body e Fix SMP butt weld tied constraints CONSTRAINED_BUTT_WELD The constraints were imposed on the accelerations before the accelerations had been updated from forces to accelerations As a result the welds could open up The welds are now properly constrained e CONSTRAINED_TIEBREAK and CONSTRAINED_BUTT_WELD supported for selective mass scaling See CONTROL_TIMESTEP 4 7 CONTROL e Fix a bug that could happen when a part set for thickness update is used PSSTUPD in CONTROL_SHELL e For initialization by prescribed
16. e Add new failure model OPT 11 to MAT SPOTWELD MAT_ 100 for beam elements e Add three new failure criteria for shell elements to MAT_ ADD EROSION on optional card 4 columns 6 8 LCEPS 12 load curve in plane shear strain limit vs element size LCEPS 13 load curve cross thickness shear strain limit vs element size LCEPSMX load curve in plane major strain limit vs element size e Add new capability to MAT_ADD_EROSION damage model GISSMO Strain rate scaling curve LCSRS can now contain natural logarithm values of Ove Arup and Partners Limited June 16 June 2014 Page 25 LS DYNA Software Release Notes Version 971 R7 1 1 strain rates as abscissa values This is automatically assumed when the first value is negative e Add new parameter NHMOD to MAT_266 The constitutive model for the isotropic part can now be chosen NHMOD 0 original implementation modified Neo Hooke NHMOD 1 standard Neo Hookeon as in umat45 New keyword DEFINE_TABLE_MATRIX is an alternative way of defining a table and the curves that the table references from a single unformatted text file e g as saved from an Excel spreadsheet e Change long format so that all data fields are 20 columns and each line of input can hold up to 200 columns In this way the number of input lines is the same for long format as for standard format 8 variables per line in long format 160 columns 10 variables per line in long format 200 columns e Add a new
17. element shapes midnodes exactly on straight lines between corner nodes Le is the same as before But for curved shapes Le gets reduced e Reduce noise in ELEMENT_BEAM_PULLEY when used together with stiff truss cable elements by making relevant quantities double precision e Add MAT_103 and MAT_187 for tet type 13 e Fix for combination of DEFINE_ ELEMENT DEATH and INCLUDE_TRANSFORM e Offset for EID SID ideoff idsoff was not correct e Improvements and fixes for ELEMENT_BEAM_PULLEY tighten tolerances for force equilibrium at the pulley slip correct history update after element passes through pulley swap correct velocity of swapped node correct swap calculations for cable elements with MAT_071 and E lt 0 e Fix for tetrahedral solid element type 4 SECTION SOLID ELFORM 4 The combination with DAMPING_PART_STIFFNESS was not working correctly leading to wrong stresses or even error termination Ove Arup and Partners Limited June 16 June 2014 Page 40 LS DYNA Software Release Notes Version 971 R7 1 1 e Shell formulations 25 26 and 27 are now fully supported in writing to dynain file INTERFACE _SPRINGBACK_LSDYNA e Shell formulations 23 quad and 24 triangle can now be mixed in a single part When ESORT 1 in CONTROL_SHELL triangular shells assigned by SECTION_SHELL to be type 23 will automatically be changed to type 24 e New element formulation ELFORM 45 in SECTION SOLID Tied Meshf
18. in SECTION_BEAM ELEMENT_BEAM_ ELBOW Ove Arup and Partners Limited June 16 June 2014 Page 26 LS DYNA Software Release Notes Version 971 R7 1 1 NEIPB in DATABASE_EXTENT_BINARY to direct output of elbow loop stresses to d3plot Otherwise output goes to ASCII file elbwls k Supported by a subset of material models including mats 3 4 6 153 195 e Add discrete element option DE to DATABASE_TRACER Includes variable RADIUS average result of all e RADIUS gt 0 Reports the average result of all DE particles in a spherical volume having radius RADIUS and centred at the tracer e RADIUS lt 0 Reports result of the closest particle to the tracer If a tracer node NID is given then the tracer moves with this node The node must belong to a DES e Add option LONG to PART_COMPOSITE and ELEMENT SHELL COMPOSITE In contrast to COMPOSITE one integration point is defined per card This is done to allow for more information e g new variable ply id e Add support of MAT_ADD_EROSION option NUMFIP lt 0 for standard non GISSMO failure criteria Only for shells e Add viscoplastic behaviour to MAT_157 i e parameter LCSS can now refer to a table with strain rate dependent yield curves e Add singular finite element with midside nodes for 2D plane strain fracture analysis ELFORM 55 in SECTION_SHELL This is an 8 noded element and can induce a singular displacement field by moving mid side nodes to quarter locations 4 2 A
19. only Fixed bug in d3part output triggered by MAXINT lt 0 in DATABASE_EXTENT_D3PART Fixed bad report of beam integration rules in d3hsp Ove Arup and Partners Limited June 16 June 2014 Page 8 LS DYNA Software Release Notes Version 971 R7 1 1 Fixed bug for encrypted data in adativity run if KEYWORD line exists in the include files Fixed bug in intfor database DATABASE_BINARY_INTFOR if geometric CONTACT_ENTITY exists Fixed bug in small restart whereby intfor database is overwritten Fixed bug in status out file if no sw2 or glstat DATABASE_GLSTAT is requested Fixed bug in Material Model Driver if number of load curves is less than 9 Fixed bug whereby correct IRID of INTEGRATION_BEAM was not output to d3hsp Increased default binary file size scale factor x from 7 to 1024 That means the default binary file size will be 1 Gb for single version and 2 Gb for double version Add echo of new max frequency of element failure summaries flag FRFREQ in CONTROL_OUTPUT to d3hsp file Support in LSDA binout output for new pllyout file DATABASE PLLYOUT ELEMENT_BEAM_PULLEY 3 11 ELEMENT Allow degenerated hexahedrons pentas for cohesive solid elements ELFORM 19 20 that evolve from an extrusion of triangular shells The input of nodes on the element cards for such a pentahedron is given by N1 N2 N3 N3 N4 N5 N6 N6 Fixed error in reading ELEMENT_SHELL COMPOSITE Fixed bug in SECTION_SHE
20. only but both for explicit and implicit simulations Using an implicit solution when damage is activated may trigger a slow convergense IMFLAG 4 or 5 can be useful e Added an option in MAT_266 MAT_TISSUE_DISPERSED so that the user can tailor the active contribution with a time dependent load curve instead of using the internal hardcoded option See ACT 10 in the User s Manual e Fixed a problem related to the shell thickness update ISTUPD in CONTROL SHELL and MAT_255 e MAT_173 MAT_MOHR_ COULOMB is available in 2D e Fix a bug affecting evolution type failure IDAM of MAT_ADD_EROSION lt X0 is applied to shell element of MAT 015 e Fix a bug in the general Darcy s law pore air flow calculation triggered when the pore air pressure gradient fell outside the range of the input curve FMAT_ADD_PORE_AIR e Fix a bug in the user defined failure subroutine matusr_24 optional for MAT 024 and others whereby the failure flag fail i was stored incorrectly e Enable MAT_103 and MAT_104 to discretize the material load curves according to the number of points specified by LCINT in CONTROL SOLUTION e Implement Prony series up to 18 terms for shells using MAT 076 MAT GENERAL VISCOELASTIC e Added DEFINE_STOCHASTIC_VARIATION and the STOCHASTIC option for MATSs 10 15 24 81 98 for shells solids and type 13 tets This feature defines a stochastic variation in the yield stress and damage failure of the afore
21. particle based blast PARTICLE_BLAST Include random distribution of initial air molecules Modify algorithm to account for the non thermally equilibrated state of high velocity gas e Improve particle contact method for particle based blast loading PARTICLE_BLAST e Fixed bug which affected the way forces were handled for triangular segments subjected to LOAD_BLAST_ENHANCED 4 5 CONTACT e Increase MPP search distance for tied contacts to include slave and master thicknesses Ove Arup and Partners Limited June 16 June 2014 Page 29 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix problem of MPP CONTACT_ _TO_SURFACE missing some nodes during bucketsort which led to unexpected penetration e Fix summation of master node forces for output to ncforc file This affects only MPP constraint based contacts e Fix CONTACT_2D_AUTOMATIC_TIED and TIED_ONE_WAY to tie only to exterior surfaces of solid elements with materials that may erode This is mainly an issue for adaptive remeshing when after remeshing nodes may lie behind the surface to be tied and nodes to be tied lie on or near the line between surface nodes and those behind the surface e Fix a memory error in SMP segment based contact that could happen occasionally when reading the history data back from the adaptivity file when h adaptivity is active e Fix anerror in MPP CONTACT 2D AUTOMATIC SURFACE TO SURFACE and SINGLE SURFACE that could cause forces that
22. plot state after each 3D tetrahedral remeshing operation CONTROL_REMESHING to reduce volume of output Generate disbout output DATABASE_DISBOUT for MPP and SMP binout files Extend DATABASE MASSOUT to include option to output mass information on rigid body nodes Bug fixes for D3PLOT e Fixed bug in title of d3plot for double precision and LSTC_BINARY ieee32 format e Fixed bug in d3plot if a solid part is removed e Fixed bug in d3plot if any part is removed for 10 noded tetrahedral element e Fixed bug where hiding a part from d3plot removed its effect from the model A removed part was also removed from contact definition e Fixed bug in binary file output d3plot d3drlf if the size of first geometry file is bigger than the defined size e Fixed improper warning message if any part is removed from d3plot Added new keyword CHANGE_OUTPUT for full deck restart to override default behaviour of overwriting existing ASCII files For small restart this option has no effect since all ASCII output is appended to the result of previous run already Added new option NEWLENGD to 2nd field of 3rd card of CONTROL_OUTPUT to write more detailed legend in ASCII output files At present only rcfore and jntforc are implemented Fixed bug in writing dyna inc INCLUDE_TRANSFORM for ELEMENT_DISCRETE ELEMENT_SOLID and ELEMENT_TSHELL Fixed bug in d3part files DATABASE_BINARY_D3PART triggered by MAT_ADD_EROSION SMP
23. running double precision executable for frequency domain features when 32 IEEE version databases are used e FREQUENCY_DOMAIN_FRF e Remove redundant scaling of the modal stress and nodal force e Add output of modal contribution ratio in d3hsp for FRF e Fix in running FRF in MPP when response direction is given in VID e FREQUENCY_DOMAIN_ACOUSTIC_BEM e Fix Bug involving boundary elements defined with shell elements e Fix IUNITS for impedance boundary condition e Add database D3ACS for collocation acoustic BEM to show the surface pressure and normal velocities Ove Arup and Partners Limited June 16 June 2014 Page 41 LS DYNA Software Release Notes Version 971 R7 1 1 Fix running ATV Acoustic Transfer Vector Fix running FFT for acoustic BEM in MPP Add frequency domain nodal or element velocity output for acoustic BEM Fix locating multi connection edges for triangular elements for acoustic BEM e Update acoustic BEM to allow using DEFINE_CURVE to define the output frequencies e Implement a scheme to map velocity boundary condition from dense BEM mesh to coarse mesh to speed up the computation e Add user node ID for acoustic field points in D3ATV Now D3ATV is given for multiple field points and multiple frequencies e Add database D3ATV for acoustic transfer vector binary plot e Fix bug in random vibration with cross PSD in MPP e Implement an option to read in user specified nodal velocity history data for
24. s temperature created by the particle impact e AIRBAG PARTICLE SENSOR_SWITCH_SHELL_TO_VENT allows user to input load curve to control the venting using choking flow equation to get proper probability function for vents Therefore this vent will have the same vent rate as real vent hole e Add new option NP2P in CONTROL_CPM to control the repartition frequency of CPM particles among processors MPP only Ove Arup and Partners Limited June 16 June 2014 Page 5 LS DYNA Software Release Notes Version 971 R7 1 1 Enhance AIRBAG_PARTICLE to support a negative friction factor FRIC or PFRIC in particle to fabric contact Particles are thus able to rebound at a trajectory closer to the fabric surface after contact Use heat convection coefficient HCONV and fabric thermal conductivity KP to get correct effective heat transfer coefficient for heat loss calculation in CPM method h_effective 1 1 H shell_thickness kP If KP is not given H will be used as effective heat transfer coefficient Extend CPM inflator orifice limit from 100 to unlimited AIRBAG_ PARTICLE Support dm_in_dt and dm_out_dt output to CPM chamber database DATABASE_ABSTAT to allow user to study mass flow rate between multiple chambers 3 6 BOUNDARY Previously the number of ships rigid bodies in BOUNDARY_MCOL as specified by NMCOL was limited to 2 Apparently this was because the code had not been validated for more than 2 rigid bodies bu
25. the combined macro element instead of the individual lengths of each element e Shell formulation 14 area weighted axisymmetric gave bad stress distributions for LOAD_BODY_Y for implicit e Enable writing of midside nodes to d3plot or 6 and 8 node quadratic shell elements e Automatically renumber shells in a consistent manner for ELEMENT_SHELL_SOURCE_SINK Ove Arup and Partners Limited June 16 June 2014 Page 38 LS DYNA Software Release Notes Version 971 R7 1 1 e In CONTROL_ SHELL add options for deletion of shells based on diagonal stretch ratio STRETCH w mode amplitude in degrees W MODE e Fix bugs with null shells for quadratic element formulations e Fix bug for spotweld beam type 9 when attached to rigid body e Fix incorrect rotational displacements for accelerometers in nodout when intopt 1 in ELEMENT_SEATBELT_ACCELEROMETER e Fix incorrect stresses output to eloutdet for thick shell formulation 5 e Fix incorrect moments in secforc when using beam forumulation 5 e Fix part not found error during keyword input phase when using ELEMENT_SHELL_COMPOSITE e Fix input error when using standard cross sections for Nastran beams input via INCLUDE_NASTRAN e Fix dynain output for solid element type 3 and 4 e Fix incorrect y stress to dynain file when using iop 3 for thick shells e Fix no ELEMENT_SHELL output to dynain file if the NOTHICKNESS option is present in INTERFACE SPRINGBACK LSDYNA_NOTH
26. well as shell formulations 2 16 IDOF 3 e Use SEGANG in CONTROL_REMESHING to define positive critical angle unit is radian to preserve feature lines in 3D tetrahedral remeshing ADPOPT 2 in PART e Fix remesher failure problem in 3D axisymmetric orbital remeshing ADPOPT 3 in PART when jobid is specified by user e For 3D solid adaptive remeshing including ADPOPT 2 and ADPOPT 3 PART the old mesh will be used automatically if the remesher fails generating a new mesh e Fix bug for pore pressure analysis FCONTROL_PORE_FLUID The bug affected time dependent or steady state analysis types when elements had reached the suction limit A symptom was occasional sudden non physical water flows at a particular node or nodes e Add option INTPERR on CONTROL SHELL Optional Card 3 Column 8 By default warning messages are written in case of non matching number of integration points between INITIAL_STRESS_SHELL and SECTION_SHELL Now with INTPERR 1 LS DYNA can terminate with an error e Add variable DITRACE on CONTROL ADAPTIVE The user can now force a plot state to d3plot just before and just after an adaptive step This option is necessary for tracing particles across adaptive steps using LS PrePost e The flag RBSMS on CONTROL_RIGID is now active for regular and selective mass scaling to consistently treat interfaces between rigid and deformable bodies e By putting MINFO 1 on CONTROL_OUTPUT penetration info is
27. were not normal to the surface This would mostly be a problem with sliding parts meshed with 2D solid elements Set CT2CN of CONTACT_AUTOMATIC_SURFACE_TIEBREAK to default of 1 0 as given in the User s Manual e CONTACT now works for parts refined using CONTROL_REFINE_SOLID or CONTROL_REFINE_SHELL e MPP fixes for contact option SRNDE CONTACT optional card E field 4 e Improve calculation of shell element contact segment thicknesses particularly at material boundaries e Fix CONTACT_ORTHO_FRICTION in SMP which was behaving incorrectly in some cases e Fix ordering inconsistencies between the writing of contact data in the keyword and reading of that data in the structured reader In particular ORTHO_FRICITON could not be used in the same input deck as either MORTAR contact or INTERFERENCE contact e MPP Add output to rcforc file for CONTACT_AUTOMATIC_TIEBREAK to record the of nodes tied and the total tied area e MPP Add calculation of contact gap for master side of FORMING contact e MPP groupable contact Remove deleted slave nodes from contact e MPP Add support for table based friction FS 2 0 to groupable contact e Minor fixes to MPP single surface contact initial penetration checking Skip nodes with parametric coordinates more than 50 outside the segment These are ignored during contact anyway so no point in reporting initial penetrations for them Also fix remaining distance to midsurface report if
28. work well with INITIAL_FOAM_REFERENCE_GEOMETRY so that initial hourglass energy is properly calculated and foam will spring back to the initial geometry 3 14 LOAD Accommodate erosion of thin shells in FLOAD _BLAST_ENHANCHED Ove Arup and Partners Limited June 16 June 2014 Page 10 LS DYNA Software Release Notes Version 971 R7 1 1 LOAD_VOLUME_LOSS has been changed such that after the analysis time exceeds the last point on the curve of volume change fraction versus time the volume change is no longer enforced LOAD_BODY_POROUS new option AOPT added to assign porosity values in material coordinate system Added LOAD_SEGMENT_FILE 3 15 SENSOR Add new sensor definition SENSOR DEFINE ANGLE This card traces the angle formed between two lines SENSOR_DEFINE_NODE can be used to trace the magnitude of nodal values coordinate velocity or accleration when VID is 0 or undefined Add two new parameters to SENSOR_DEFINE_ELEMENT scale factor and power so that user can adjust the element based sensor values strain stress force 3 16 MAT Change history variables 10 12 in MAT_054 MAT_ENHANCED_COMPOSITE_DAMAGE thin shells only to represent strains in material coordinate system rather than in local element coordinate system This is a lot more helpful for postprocessing issues This change should not lead to different results other than due to different round off errors New features and enhancements t
29. ATABASE_BINARY_D3PLOT Ove Arup and Partners Limited June 16 June 2014 Page 36 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix seg fault error when using DATABASE_DISBOUT for large number of discrete beam elements MPP only e Fix zero external work in glstat when using only one BOUNDARY_PRESCRIBED_MOTION_RIGID keyword card e Fix incorrect forces in bndout when using BOUNDARY_PRESCRIBED_MOTION_RIGID e Fix missing eroded hourglass energy from glstat in binout file e Fix incorrect beam resulants in d3part MPP only e Fix incorrect output for local resultant normal t to disbout DATABASE_DISBOUT for discrete beams MPP only e Fix incorrect forces in bndout file for rigid bodies when using LOAD_RIGID_BODY e Enable CONTROL_OUTPUT flag EOCS which wasn t having any effect on the shells output to elout file e Fix incorrect label in elout file when doing implicit linear analysis MPP Only e Fix incorrect stonewall energy in gldata due to nodes belonging to deleted elements impacting on the RIGIDWALL e Fix incorrect acceleration output to d3plot amp nodout file when using e CONSTRAINED_NODAL_RIGID_BODY with DRFLG set to nonzero e Fix difference between nodout and nodouthf results when using iaccop 1 or 2 in CONTROL_OUTPUT and intopt 0 in ELEMENT_SEATBELT_ACCELEROMETER e Fix d3drlf for MPP double precision Data could not be postprocessed e Fix bug in option for part removal from d3p
30. ATIC MORTAR the contact thickness is now consistent with other automatic contacts i e thickness on PART CONTACT overrides the thickness on CONTACT _ which in turn overrides the physical thickness e For CONTACT MORTAR the bucket sort is made much more efficient resulting in significant speedup in large scale contact simulations e For CONTACT AUTOMATIC MORTAR the mortar contact now supports contact with the lateral surface of beam elements e For CONTACT MORTAR IGAP GT 1 stiffens the mortar contact for large penetrations The mortar contact has a maximum penetration depth DMAX that depends on geometry and input parameters if penetration is larger than this value the contact is released To prevent this release which is unwanted the user may put IGAP GT 1 which stiffens the behaviour for penetrations larger than 0 5 DMAX without changing the behaviour for small penetrations This should hopefully not be as detrimental to convergence as increasing the overall contact stiffness e Turn off normal check for CONTACT_DRAWBEAD when ORIEN 4 CONTROL_CONTACT 4 6 CONSTRAINED e Fix for MPP selective mass scaling when used together with CONSTRAINED_TIE BREAK e Fix incorrect mass amp inertia for rigid body part when IFLAG 1 in CONSTRAINED_EXTRA_NODES e Allow CONSTRAINED_INTERPOLATION to use node set to define the independent nodes e Fix incorrect inertia output to d3hsp for discrete parts defined
31. Boyce viscoelastic rubber model is now available in explicit and implicit analysis as MAT_269 MAT_BERGSTROM_BOYCE_RUBBER The Arruda Boyce elastic stress is augmented with a Bergstrom Boyce viscoelastic stress Ove Arup and Partners Limited June 16 June 2014 Page 14 LS DYNA Software Release Notes Version 971 R7 1 1 corresponding to the response of a single entangled chain in a polymer gel matrix Added a new parameter IEVTS to MAT USER DEFINED MATERIAL MODELS MAT_041 050 IEVTS is optional and is used only by thick shell formulation 5 It points to the position of E a in the material constants array Following E a the next 5 material constants must be E b E c v ba v ca and v cb This data enables thick shell formulation 5 to calculate an accurate thickness strain otherwise the thickness strain will be based on the elastic constants pointed to by IBULK and IG Implemented enhancements to fabric material MAT_034 FORM 14 as follows Stress strain curves may include a portion for fibres in compression When unload reload curves with negative curve ID are input curve stretch options the code that finds the intersection point now extrapolates the curves at their end rather than simply printing an error message if an intersection point cannot be found before the last point in either curve Fix bug in output of crack data for TMAT 084 MPP only 3 17 ALE Map 1D to 3D by beam volume averaging the 1D data o
32. C This is an orthotropic material with temperature dependent properties defined by load curves Change in structured file format for control card 27 first thermal control card Several input variables used 15 format limiting their value to 99 999 A recent large model exceeded this limit The format was changed to 110 This change is not backward compatible Old structured input files will no longer run unless control card 27 is changed to the new 110 format This change does not affect the KEYWORD file Fix coupled thermal solution for problems using thick thermal shells Fix computation of assembled stiffness matrix for the unsymmetrical case Should only affect coupled thermal implicit mechanical problems in MPP Thermal material MAT_T07 MAT_THERMAL_CWM for welding simulations to be used in conjunction with MAT_270 MAT_CWM mechanical counterpart 3 22 MPP Modify decomposition costs of MAT_181 and MAT_183 Fix transfer of mes files at code termination Fix problem in CONTROL_MPP_DECOMPOSTION_ARRANGE_ PARTS Ove Arup and Partners Limited June 16 June 2014 Page 17 LS DYNA Software Release Notes Version 971 R7 1 1 Fix for MPP CONTACT_GUIDED_CABLE in case some processors had no slave nodes Fixes to MPP reporting of element controlling the time step Introduction of new timing routines and summary at termination Echo MPP contact is groupable flag to d3hsp For MPP groupable tied contact only output of part
33. DAIMLER_CAMANHO e Add pentahedra cohesive solid element types TYPE 21 amp 22 e Type 21 is the pentahedra version of Type 19 e Type 22 is the pentahedra versionof Type 20 e Using ESORT gt 0 in CONTROL_SOLID will automatically sort out the pentahedra elements 19 to 21 and 20 to 22 e Add DEFINE_DE_BY_PART to define control parameters for DES by part ID including damping coefficient friction coefficient spring constant etc If defined it will overwrite the parameters in CONTROL_DISCRETE_ELEMENT e Add new feature for MAT_030 MAT_SHAPE MEMORY as optional 3rd card Curves or tables strain rate dependency can be defined to describe plastic loading and unloading behaviour e New feature for ELEMENT_BEAM_PULLEY Beam elements BIDI and BID2 can now both be defined as 0 zero In that case adjacent beam elements are automatically detected Therefore the first two beam elements with nodal distance lt 1 0e 6 to the pulley node PNID will be chosen Ove Arup and Partners Limited June 16 June 2014 Page 24 LS DYNA Software Release Notes Version 971 R7 1 1 e Add new feature to MAT_ADD_EROSION s damage model GISSMO By default damage is driven by equivalent plastic strain Now users can optionally define another history variable as driving quantity by setting DMGTYP according to the manual e Add volumetric plastic strain to MAT_187 as history variable 6 e Add internal energy calculation for ELEMENT_BEAM_
34. EMBEDSHELL can coexist Add the option to give a title to ICFD parts using ICFD_PART_TITLE ICFD_INITIAL Add a reference pressure pressurization pressure for when no pressure is imposed on the boundaries Add the initialization of all nodes at once by setting PID 0 Add the non inertial reference frame implementation defined by the keyword ICFD_DEFINE_NONINERTIAL Add several new state variables to LSO Please refer to the LSO manual to see how to print out the list of supported variables Add title option to the material keyword using ICFD_MAT_TITLE Fix the automatic time step calculation in the presence of boundary layer meshes Add support for FSI with thick shells 2D shells are now supported for FSI in MPP In the past only beams could be used in MPP and beams and shells could be used in SMP The keyword ICFD_CONTROL_FSI has a new field to control the sensitivity of the algorithm to find the solid boundaries used in FSI calculations The 2D mesh now generates semi structured meshes near the boundaries Fix the generation of boundary layer meshes on multiple embedded shells that are connected with each other or to the walls Add heat flux boundary condition using ICFD_BOUNDARY_FLUX_TEMP Add divergence free and Space Correlated Synthetic Turbulence Inlet Boundary Condition for LES Smirnov et al using ICFD_BOUNDARY_PRESCRIBED_VEL ICFD_BOUNDARY_PRESCRIBED_ VEL Add inflow velocities using the wall normal and
35. ESE BOUNDARY_CYCLIC_SET CESE BOUNDARY_CYCLIC_PART Add code for 2D CESE sliding boundary conditions For 2D correct CESE point output using LSO Correct CESE moving mesh restart Correct CESE chemistry communications involving species data MPP Only Add support in CESE FSI for 2D shells MPP Only Add support for CESE FSI with thick shells Add 2D amp 2D axisymmetric cases in the CESE FSI solver including both immersed boundary method amp moving mesh method 4 25 CHEMISTRY Add the CSP reduced chemistry model with 0D 2D and 3D combustion The 2D and 3D combustion cases couple with the CESE compressible flow solver Add the G scheme reduced chemistry model only for 0D combustion Add two different reduced chemistry models Ove Arup and Partners Limited June 16 June 2014 Page 61 LS DYNA Software Release Notes Version 971 R7 1 1 e The Computational Singular Perturbation CSP reduced model is implemented with existing compressible CESE solver The CSP is now working on 0 dimensional constant volume and pressure combustion 2 D and 3 D combustion problems e The new reduced chemistry model G scheme is implemented but currently works only 0 dimensional problems such as constant combustors 4 26 MISCELLANEOUS Fix spurious error KEY 141 when using RIGIDWALL_GEOMETRIC _DISPLAY and the optional card 6 for the _DISPLAY option is not defined Fix incorrect nodes generated by SET NODE GENERAL with O
36. ICKNESS e Write complete history variables to dynain file for 2D solids using MAT_NULL and equation of state e Fix incorrect transformation when using DEFINE_TRANSFORMATION with option ROTATE and FELEMENT_BEAM_ORIENTATION e Fix thick shell full deck restart for thick shell formulation 5 or for any thick shell form that uses laminated shell theory e Fix contact and mass and inertia reporting to d3hsp of composite parts that are created by ELEMENT_SHELL_COMPOSITE or ELEMENT_TSHELL_COMPOSITE The contact was failing when composite parts or composite parts in a part set were used to define the master or slave surface The composite parts were not included in the surface leading to pass through This affects all contact types e Fix a problem with using ELEMENT_SHELL_COMPOSITE with ELEMENT_BEAM_PID The PID option allows the user to specify the parts to which welds should be tied When using ELEMENT_SHELL_COMPOSITE we create new part numbers for the elements so we fail to find a match with the PID and the ties failed to form e Fix strain output of shell formulations 13 15 when 1 point integration is used with objective stress update CONTROL _ACCURACY The output to d3plot and elout was not correct for elements that rotate e Fix strain output of shell formulation 13 2D solid with 2x2 in plane integration when objective strain update is active Rotating elements were Ove Arup and Partners Limited June 16 June 2014 Page
37. IDs for MAT_100_DA spotwelds DEFINE_CONNECTION_PROPERTIES was missing Bodies using MAT_RIGID_DISCRETE were never expected to share nodes with non rigid bodies but this now works in MPP Fix variable initialization that resulted in MPP writing bad intfor data DATABASE_BINARY_INTFOR on some systems There is no longer any built in limitation on the number of processors that may be used in MPP Echo contents of the MPP pfile including keyword additions to the d3hsp and mes0000 files Add new keyword CONTROL_MPP_PFILE which allows for insertion of text following this command to be inserted into the MPP pfile p pfile Change in MPP treatment of CONSTRAINED_TIE BREAK They now share a single MPI communicator and a single round of communication This should improve performance for problems with large numbers of these without affecting the results Fix MPP decomposition error for LOAD_BODY_GENERALIZED that could result in forces being applied incorrectly Fix summation order for CONSTRAINED LAGRANGE IN SOLID CTYPE 2 for rebar This will avoid the round off error causing slightly different shared nodal velocity on echo processor 3 23 FORMING Added two input variables for CONTROL_FORMING_ONESTEP simulation e TSCLMIN is a scale factor limiting the thickness reduction e EPSMAX defines the maximum plastic strain allowed e CONTROL_FORMING_ONESTEP added output of strain and stress tensors for onestep solver to
38. LL for seatbelt where SECID was input as a label Add new option to activate drilling constraint force for shells in explicit calculations This can be defined by parameters DRCPSID part set and DRCPRM scaling factor on CONTROL_SHELL 3 12 FREQUENCY_DOMAIN FREQUENCY_DOMAIN_ACOUSTIC_BEM e Addedanoptionto output real part of acoustic pressure in time domain e Enabled BEM acoustic computation following implicit transient analysis e Implemented coupling between steady state dynamics and collocation acoustic BEM e Implemented Acoustic Transfer Vector ATV to variational indirect BEM acoustics Enabled boundary acoustic mapping in BEM acoustics FREQUENCY_DOMAIN_ACOUSTIC_FEM Added boundary nodal velocity to binary plot file d3acs Implemented pentahedron elements in FEM acoustics Enabled using boundary acoustic mapping in FEM acoustics Ove Arup and Partners Limited June 16 June 2014 Page 9 LS DYNA Software Release Notes Version 971 R7 1 1 FREQUENCY_DOMAIN_FRF Updated FRF to include output in all directions VAD2 4 Added treatment for FRF with base acceleration node id can be 0 FREQUENCY_DOMAIN_RANDOM_ VIBRATION Updated calculating PSD and RMS von Mises stress in random vibration environment based on the Report of Sandia National Laboratories 1998 FREQUENCY_DOMAIN_RANDOM_VIBRATION_FATIGUE e Implemented an option to incorporate initial damage ratio in random vibration fatigue e FREQUENC
39. LS DYNA Software Version 971 R7 1 1 Release Notes LS DYNA Software Contents 1 Introduction 2 Release notes for LS DYNA version R7 3 Release notes for LS DYNA version R7 0 0 3 1 New Electromagnetic EM solver 3 2 New Compressible CFD CESE solver 3 3 New Incompressible CFD ICFD 3 4 New Volume Mesh generator 3 5 AIRBAG 3 6 BOUNDARY 3 7 CONTACT 3 8 CONSTRAINED 3 9 CONTROL 3 10 DATABASE 3 11 ELEMENT 3 12 FREQUENCY_DOMAIN 3 13 INITIAL 3 14 LOAD 3 15 SENSOR 3 16 MAT 3 17 ALE 3 18 SPH 3 19 EFG 3 20 IMPLICIT 3 21 THERMAL 3 22 MPP 3 23 FORMING 3 24 DISCRETE ELEMENT METHOD 3 25 ISOGEOMETRIC ELEMENTS 3 26 MISCELLANEOUS 4 Release notes for LS DYNA version R7 1 0 4 1 NEW FEATURES GENERAL 4 2 AIRBAG 4 3 BOUNDARY 4 4 BLAST 4 5 CONTACT 4 6 CONSTRAINED Ove Arup and Partners Limited June 16 June 2014 Release Notes Version 971 R7 1 1 Page OO 0 NY CS A Vi Li HWW W N N e e o e o cr cnr O OO O uN kx i Uli ai KF DO CO 22 22 27 28 29 29 33 LS DYNA Software 4 7 CONTROL 4 8 DATABASE 4 9 ELEMENT 4 10 FREQUENCY_DOMAIN 4 11 INITIAL 4 12 LOAD 4 13 SENSOR 4 14 MAT 4 15 ALE 4 16 SPH 4 17 DISCRETE ELEMENT METHOD 4 18 EFG 4 19 IMPLICIT 4 20 THERMAL 4 21 MPP 4 22 FORMING 4 23 ISOGEOMETRIC ELEMENTS 4 24 ICFD INCOMPRESSIBLE FLOW SOLVER 4 25 CHEMISTRY 4 26 MISCELLANEOUS 5 Release notes for LS DYNA version R7
40. MPP when using staged construction e MPP spotweld deletion better handling of bad user input A spotweld node should not be defined as slave to more than one tied interface If it is it gets deleted from all but one of them Previous code would delete the spotweld if the slave node of any interface reported the slave node not tied and SPOTDEL was active Now only delete the spotweld if ALL interfaces in which it appears report it as not tied e Fix tied node checking in MPP for the IPBACK flag that affected non groupable tied contacts e MPP synchronize rigid body flags for shared nodes during DEFORMABLE_TO_RIGID switching so that these nodes are handled consistently across processors e Fix SMP bug that causes too many opened files in problems with many many adaptive steps a file was being opened every cycle but not closed Ove Arup and Partners Limited June 16 June 2014 Page 64 LS DYNA Software Release Notes Version 971 R7 1 1 6 R7 1 1 EXECUTABLES Naming convention for LS DYNA Ls dyna_Type_Precision_Revision_Machine_OS_Compiler_Misc_Options Where e Type smp mpp hyb e Precision Single s Double d e Revision r7_1_1_build_id e Machine x64 1a64 e OS redhat54 suse102 winx64 redhat59 e Compiler Intel Fortran 13 1 ifort131 PGI 605 pgi605 Intel Fortran 10 1 ifort101 e Misc Options Instruction set e sse2 e avx e avx2 DYNALIB e sharelib MPI Libraries e platformmpi pmpi intelmpi im
41. Mises stress SIGVM can now be a function of strain rate by specifying a negative load curve ID Rotational treatment in MAT_034 MAT_FABRIC supported with airbag initial reference geometry 4 15 ALE Implement DATABASE_ALE to write time history data volume fractions stresses for a set of ALE elements Not to be confused with DATABASE ALE MAT For INITIAL_ALE_ MAPPING write the mapping file if the sense switches sw1 sw2 stop or quit are invokes or if the criteria set by TERMINATE stop the run SECTION ALEID and SECTION_ALE2D now work on multiple processors Ove Arup and Partners Limited June 16 June 2014 Page 51 LS DYNA Software Release Notes Version 971 R7 1 1 For CONSTRAINED_ LAGRANGE _IN_SOLID ctype 5 support tshell elements CONSTRAINED_LAGRANGE_IN_SOLD ctype 4 5 now converts friction energy to heat Note it only works for ALE elform 12 Clean up MPP memory allocation for CONSTRAINED_LAGRANGE_IN_SOLID CTYPES 4 and 5 Implement DELETE_PART in small restarts for ALE2D parts 4 16 SPH Fix output suppression PSETID in DATABASE BINARY_D3PLOT and DATABASE BINARY_D3PART for SPH parts Correct nodal contact stiffness for SPH nodes in MPP NODES_TO_SURFACE type contacts both SOFT 0 and 1 Results are now consistent with SMP Enable SPH tracer particle DATABASE _TRHIST output for SMP in Isda database Add conversion of frictional contact energy into heat when doing a coupled therm
42. NTROL_FORMING_TIPPING which tipped all parts with the same amount regardless of the PID SID specified Fix a bug in DEFINE_CURVE_FLC Ove Arup and Partners Limited June 16 June 2014 Page 58 LS DYNA Software Release Notes Version 971 R7 1 1 Allow handling of nano sized structure and assorted bugs in CONTROL_FORMING_AUTOPOSITION Positioning distances were off Fix various bugs in linear unflanging to improve convergence for CONTROL_FORMING_UNFLANGING Fix a bug in multiple bead definition related to error message duplicate node set found etc This is related to NBEAD parameter in CONTACT_DRAWBEAD Fix segmentation fault when using DEFINE_FORMING_BLANKMESH Fix a bug for trimming when the trimming curves are not self closed in DEFINE_CURVE_TRIM_NEW Fix a bug in CONTROL_FORMING_ONESTEP with encrypted material input 4 23 ISOGEOMETRIC ELEMENTS Add error message if polynomial order of NURBS Patch exceeds 9 ELEMENT_SHELL_NURBS_PATCH Set INT 1 full integration as default if linear NURBS elements are used ELEMENT_SHELL_NURBS_PATCH Otherwise it may lead to instability Fix transverse shear bug in isogeometric shell form 1 Isogeometric analysis contact is available in MPP Normalize tangent vectors for local coordinate system for the rotation free isogeometric shells Add support for dumping shell internal energy density for isogeometric shells ELEMENT_SHELL_NURBS_PATCH via interpolation shell
43. PTION PART when used together with PART_MOVE and IFSET 1 i e part set id is used Fix bug affecting TERMINATION_DELETED_SOLIDS_SET and TERMINATION_DELETED_SHELLS_SET Only the first part in the part set was counted and used to determine the termination Fix a problem of random termination when PART_COMPOSITE is used with adaptive remeshing Fixes to NODE_TRANSFORM which was not working correctly with SET_NODE_GENERAL Several cases of input processing speedups due to the use of better sorting algorithms Fix integer overflow in reporting of memory usage if the values are gt 1GW report in MW instead of in words Fix some bugs for INCLUDE_NASTRAN The translated LS DYNA deck is now written to dyna inc Fix input error during keyword input phase when using alphanumeric label for SECTION and SET NODE GENERAL e Fix input error during keyword input phase when PLOTEL 1 in CONTROL_RIGID and number of CONSTRAINED_NODAL_RIGID_BODIES exceeds 500 Fix input error while reading in structured file during restart and using long s Fix skipping of the line following an inadvertent blank line when reading ELEMENT_SOLID Fix warning message for empty part set when using SET_PART_LIST_GENERATE Fix extra DEFINE_CURVE card in dyna inc file when using DEFINE_CURVE_SMOOTH with INCLUDE_TRANSFORM Fix input error when using SET PART ADD to define part set for DEFORMABLE_TO_RIGID Ove Arup and Partners L
44. PULLEY e Add viscoplastic option to MAT_157 new parameter VP on Card 5 Column 6 e Add new keyword MAT_ADD_COHESIVE which is intended to make 3D material models available for cohesive elements e Add new parameters to MAT_CABLE_DISCRETE MAT_071 MXEPS Card 2 Column 4 maximum strain at failure and MXFRC Card 2 Column 5 maximum force at failure e Add MAT_124 as potential weld partner material for PROPRUL 2 3 of DEFINE_CONNECTION_PROPERTIES e Add new material MAT_TOUGHENED_ADHESIVE_POLYMER TAPO or MAT_252 for epoxy based toughened ductile adhesives e Add new option to MAT_002_ANIS parameter IHIS on Card 4 Column 8 e IHIS 0 terms C11 C12 from Cards 1 2 and 3 are used e JHIS 1 terms C11 C12 initialized by INITIAL_STRESS_SOLID s extra history variables e Add new option to MAT_102 Instead of constant activation energy Q one can define a load curve LCQ on Card 2 Column 7 LCQ can be defined as a function of plastic strain or of temperature e Add new option to MAT_071 MAT_CABLE_DISCRETE_BEAM New parameter FRACLO Card 2 Column 3 is fraction of initial length that should be reached over time period of TRAMP That means the cable element length gets modified from LO to FRACLO LO between t 0 and t TRAMP e Add internal energy calculation for SPR models CONSTRAINED_INTERPOLATION_SPOTWELD SPR3 and CONSTRAINED_SPR2 Their contribution was missing in energy reports like glstat
45. R7 1 1 3 Release notes for LS DYNA version R7 0 0 3 1 New Electromagnetic EM solver Keyword family EM_ The keywords starting with EM refer to and control the Electromagnetic solver problem set up Solver Characteristics Implicit Double precision Dynamic memory handling SMP and MPP 2D axisymmetric solver 3D solver Automatic coupling with structural and thermal LS DYNA solvers FEM for conducting pieces only no air mesh needed FEM BEM system Solid elements for conductors shells can be insulators Solver Main Features Eddy Current a k a Induction Diffusion solver Induced heating solver Resistive heating solver Imposed tension or current circuits Exterior field Magnetic materials beta version Electromagnetic contact EM Equation of states Conductivity as a function of temperature Solver Applications Non exhaustive Electromagnetic forming Electromagnetic welding Electromagnetic bending Inductive heating Resistive heating Rail gun Ring expansions 3 2 New Compressible CFD CESE solver Keyword family CESE_ The keywords starting with CESE refer to and control the Compressible CFD solver problem set up Solver Characteristics e Explicit e Double precision Ove Arup and Partners Limited June 16 June 2014 Page 3 LS DYNA Software Release Notes Version 971 R7 1 1 Dynamic memory handling SMP and MPP 3D solver special case 2D solver and 2D axisymmetric solver Automatic coupli
46. Release Notes Version 971 R7 1 1 Change MAT_ACOUSTIC elform 8 time step calculation to more conservatively handle highly distorted tets Fix internal energy calculation for MAT_ACOUSTIC elform 8 tetrahedrons and pentahedrons wedges Add an error termination if MAT_SPOTWELD is used with OPT 6 or OPT 9 and the necessary DEFINE_SPOTWELD_RUPTURE_STRESS or DEFINE_SPOTWELD_RUPTURE_PARAMETER do not exist This prevents the use of uninitialized memory Add a warning if a user tries to use INTEGRATION_SHELL or PART_COMPOSITE with MAT_FABRIC It becomes an error if there is not thickness data provided by a SECTION card Add error message for MAT USER DEFINE MATERIAL MODELS if IBULK and IG gt LMC Extend recently developed features for MAT_054 shells to MAT_054 thick shells TYPE 3 and 5 Including SLIM input variables Orthogonal crashfront Strain rate dependent strength via DEFINE_CURVE Separate subroutine for MAT_055 Enable MAT_036 and MAT_243 to be used with user defined failure Input parameter USRFAIL 0 off 1 on is on new Optional Card 6 Column 1 Add MAT_007 for type 13 tetrahedral element SECTION_SOLID ELFORM 13 Enable regularization curve LCREGD of MAT_ADD_EROSION to be used with standard non GISSMO failure criteria Users can now define a failure criterion plus IDAM 0 plus LCREGD scaling factor vs element size to get a regularized failure criterion MAT_ADD_EROSION equivalent von
47. SURFACE_TIED The implementation includes bending and torsion Implement keyword DEFINE_DE_HBOND to define heterogeneous bond for discrete element spheres DES DES ELEMENT_DISCRETE_SPHERE with different material models can be bonded Implement keyword INTERFACE DE BOND to define multiple failure models for various bonds within one part or between different parts through the keyword DEFINE_DE_HBOND Implement DEFINE_DE_TO_BEAM_COUPLING for coupling of discrete element spheres to beam elements Add support of triangular segments for DEFINE_DE_TO_SURFACE_COUPLING Add variable MAXGAP in DEFINE_DE_BOND to give user control of distance used in judging whether to bond two DES together or not based on their initial separation 4 18 EFG Add IAT 3 in CONTROL_REMESHING_EFG which uses FEM remapping scheme in EFG adaptivity Compared to IAT 2 1 1 2 IAT 3 is faster and more robust but less accurate In SECTION SHELL EFG set the default values of DX DY to 1 1 Add control flag MM in CONTROL_REMESHING_EFG to turn on off monotonic mesh resizing for EFG 3D general remeshing ADPOPT 2 in PART 4 19 IMPLICIT Fix a bug in CONTACT_TIED_SHELL_EDGE_TO_SURFACE_CONSTRAINED_OFF SET in implicit The bug caused unreasonable contact forces when dealing with slave nodes offset from warped master segments Upgrade the implicit implementation of rack and pinion and screw joints CONSTRAINED_JOINT so the joint is driven
48. TABASE BINARY_FSILNK This feature stores coupling pressure from CONSTRAINED_LAGRANGE_IN_SOLID in a binary time history file for use in a separate model that does not include ALE e Add LOAD SEGMENT_FSILNK Use pressure loads stored in aforementioned binary time history file to load model that does not have ALE elements e Add new keyword DEFINE_SPH_DE COUPLING to allow SPH particles to contact discrete element spheres DES e Add MOISTURE option to MAT_076 solids Allows moisture content to be input as a function of time Material parameters are then scaled according to the moisture and a moisture strain is also introduced e Add RIGIDWALL_FORCE_TRANSDUCER to output forces from rigidwalls acting on node sets e Add LOG_INTERPOLATION option to MAT_024 This offers an alternate means of invoking logarithmic interpolation for strain rate effects The other way is to input the natural log of strain rate in the table LCSS e Add capability in MAT_ADD_EROSION NUMFIP lt 100 to set stress to zero in each shell integration point as it reaches the failure criterion When INUMFIPI 100 integration points have failed the shell is eroded In contrast when NUMFIP gt 0 failed integration points continue to carry full load as though they were unfailed until element erosion occurs e Add new keyword PARAMETER_TYPE for use by LS PrePost when combining keyword input files The appropriate offset is applied to each ID value defined using e
49. TABASE EXTENT BINARY e Fix incorrect strains for thick shells in eloutdet and also fix INTOUT amp NODOUT options in DATABASE_EXTENT_BINARY for thick shells e Fix NaN in secforc file when using DATABASE_CROSS_SECTION with NID specified in DEFINE_COORDINATE_VECTOR MPP only e Include eroded hourglass energy in hourglass energy in glstat file to be consistent with KE amp IE calculations so that e total energy kinetic energy internal energy hourglass energy rigidwall energy e Fix zero output to upper strains to d3plot and elout when nip 1 for shell elements e Fix d3drlf output with duplicate states at t 0 0 e Fix corrupted d3plot when using DCOMP 2 4 in DATABASE EXTENT_ BINARY in combination with DEFORMABLE_TO_RIGID SMP only e Fix corrupted d3thdt when DCOMP 3 4 in DATABASE_EXTENT_BINARY e Fix incorrect mass properties for solids in SSSTAT file when using DATABASE_SSSTAT_MASS_PROPERTY e Fix incorrect zero values output to curvout file for DEFINE_CURVE_FUNCTION using functions VM and WM e Fix output of zeros NaNs to curvout for DEFINE_CURVE_FUNCTIONs AXIAYIAZIWXIW YIWZ e Fix incorrect output to curvout file when using RCFORC function for DEFINE_CURVE_FUNCTION for MPP with ncpu gt 1 e Fix zero strain values output to curvout for DEFINE_CURVE_FUNCTION using function ELHIST for shell strains e Fix input error undefined beam when reading discrete element connectivity cards if BEAM 1 in D
50. TRBAG e Enhance AIRBAG PARTICLE Fix blockage calculation to be less sensitive to number of processors Ath field of CONTROL CPM set to 1 will check bag chamber connectivity s and error terminate if problem is detected If HOONV lt 0 IHCONVI is a curve of heat convection coefficient vs time Error terminate if chamber ID CHM in AIRBAG_PARTICLE is referenced but not defined e Add new option DECOMPOSITION for AIRBAG_PARTICLE This automatically invokes the recommended decomposition commands CONTROL_MPP_DECOMPOSITION_BAGREF if applicable and CONTROL_MPP_DECOMPOSITION_ARRANGE_ PARTS for the bag e Add new blockage option for vents in AIRBAG_PARTICLE Blockage considered eq 0 no eq l yes eq 2 yes exclude external vents eq 3 yes exclude internal vents eq 4 yes exclude all vents Ove Arup and Partners Limited June 16 June 2014 Page 27 LS DYNA Software Release Notes Version 971 R7 1 1 e Fixed a bug in AIRBAG_PARTICLE where the porosity and vent blockage option was not treated correctly in MPP and the bag behaviour was sensitive to the number of cores e Add option in CONTROL CPM to consider CPM in the time step size calculation e When using AIRBAG_PARTICLE with IAIR 2 user should keep mole particle similar between inflator gas and initial air particles to ensure the correct elastic collision If different by more than 10 code will issue warning message and provide the suggested initial
51. UTOMATIC e Fix SMP and MPP segment based contact when there are thick shells in contact and also solid element form 1 or 2 in the data The contact was assuming contiguous connectivity data for the solids and thick shells but this assumption does not work when there are solid form 1 or 2 elements e Improve the segment based single surface contact search for thick segment pairs that are too close together The code was not working well with triangluar segments This change affects models with shell segments that have thickness greater than about 2 3 of the segment length e Fix CONTACT_2D_AUTOMATIC_TIED and TIED_ONE_WAY in two ways First the positioning of slave nodes that are offset from the master surface was not done with sufficient accuracy and the offset distance could change a little Second it was possible for master nodes constrained by SPC s to move e Fix an MPP segment based contact error that occurred for 2 surface contact in both the initial penetration check and the periodic intersection check e Enable segment based quad splitting options to work when shell sets or segment sets are used to define the surface that will be split This is really a bug fix because there was no check to prevent this and the result was writing past the allocated memory for segment connectivity e Fix the detailed searching in segment based contact when edge extensions are turned off The search was failing to detect some penetration soon enou
52. Y_DOMAIN_RESPONSE_SPECTRUM e Implemented double sum methods based on Gupta Cordero coefficient modified Gupta Cordero coefficient and Rosenblueth Elorduy coefficient Updated calculating von Mises stress in response spectrum analysis Implemented treatment for multi simultaneous input spectra Improved double sum methods by reducing number of loops FREQUENCY_DOMAIN_SSD Added the option to output real and imaginary parts of frequency response to d3ssd e Added the option to output relative displacement velocity and acceleration in SSD computation in the case of base acceleration Previously only absolute values were provided Implemented keyword FREQUENCY_DOMAIN_MODE_ OPTION so that user can select the vibration modes to be used for frequency response analysis Implemented keyword SET MODE OPTION so that user can define a set of vibration modes to be used for frequency response analysis Implemented keyword FREQUENCY_DOMAIN_PATH to define the path of binary databases containing mode information used in restarting frequency domain analysis e g frf ssd random vibration 3 13 INITIAL Compute normal component of impulse for oblique plates in INITIAL MINE_IMPULSE The feature is no longer limited to horizontal plates Disable license security for INITIAL_IMPULSE_MINE The feature is no longer restricted Fix a bug of INITIAL_STRESS_SOLID for 4 noded tetrahedral elements Enabled hourglass type 7 to
53. _88920_x64_redhat54_ifort131_avx2_platformmpi tar gz 1_88920_x64_redhat54_ifort131_avx2_platformmpi_sharelib tar gz 1_88920_x64_redhat54_ifort131_avx2_platformmpi tar gz 1 1 1 88920_x64_redhat54_ifort131_avx2_intelmpi tar gz 88920_x64_redhat54_ifort131_avx2_openmpi165 tar gz 88920_x64_redhat54_ifort131_avx2_platformmpi_sharelib tar gz m RR a i n ek Ove Arup and Partners Limited June 16 June 2014 Page 66 LS DYNA Software Release Notes Version 971 R7 1 1 ls dyna_hyb_d_r7_1_1_88920_x64 redhat54_ifort131_avx2_platformmpi tar gz ls dyna_hyb_s_r7_1_1_88920_x64_redhat54_ifort131_avx2_intelmpi tar gz ls dyna_hyb_s_r7_1_1_88920_x64_redhat54_ifort131_avx2_openmpi165 tar gz ls dyna_hyb_s_r7_1_1_88920_x64_redhat54_ifort131_avx2_platformmpi_sharelib tar gz ls dyna_hyb_s_r7_1_1_88920_x64_redhat54_ifort131_avx2_platformmpi tar gz ls dyna_mpp_d_r7_1_1_88920_x64_redhat54_ifort131_avx2_intelmpi tar gz ls dyna_mpp_d_r7_1_1_88920_x64_redhat54_ifort131_avx2_mpich126 tar gz ls dyna_mpp_d_r7_1_1_88920_x64_redhat54_ifort131_avx2_mpich2 121 tar gz ls dyna_mpp_d_r7_1_1_88920_x64_redhat54_ifort131_avx2_openmpi165 tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_avx2_intelmpi tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_avx2_mpich126 tar gz ls dyna_mpp_s_r7_1_1_88920_x64_redhat54_ifort131_avx2_mpich2 121 tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_avx2_openmpi165 tar gz ope
54. a velocity magnitude using the 3rd field VAD Add the activation of synthetic turbulence using the 3rd field VAD ICFD CONTROL LOAD use to reset the body load provided by the LOAD_BODY keyword Ove Arup and Partners Limited June 16 June 2014 Page 60 LS DYNA Software Release Notes Version 971 R7 1 1 Add the option to control the re meshing frequency in both keywords see ICFD CONTROL ADAPT SIZE and ICFD_CONTROL_ADAPT ICFD CONTROL IMPOSED MOVE Impose a rigid body motion to the mesh bypassing the solid solver to avoid the FSI variable transfer in applications that need simple motions ICFD_CONTROL_TURB_SYNTHESIS control parameters for the synthetic turbulence inflow ICFD_CONTROL_OUTPUT Added the VTK option to OUTL using the number 6 By default vtk files will be written to a directory named vtk in the working directory ICFD_DATABASE_AVERAGE Allows the computation of time average variable values at given time intervals ICFD_BOUNDARY_PRESCRIBED_MOVEMESH Allows the mesh to slide on the boundaries following the Cartesian axis CESE COMPRESSIBLE FLOW SOLVER Add a PART SET option for CESE_BOUNDARY_ _PART cards Bring in more 2D mesh support both from the PFEM mesher and a user input 2D mesh via ELEMENT_SOLID with 0 for the last 4 of 8 nodes Enable the 2D ball vertex mesh motion solver for the 2D CESE solver Fix d3plot output for the 2D Eulerian and moving mesh cases Add new input cards C
55. air particle number e Fix aconflict in MAT_FABRIC when FORM 14 is used with the negative ID for unload curves and the Anagonye and Wang porosity equation for airbag leakage area on MAT_FABRIC is being used with a ATRBAG WANG _NEEFSKE control volume The conflict caused a bad venting calculation leading to NaN pressure e Enable DEFINE_CURVE_FUNCTION for SECTION_POINT_SOURCE_MIXTURE and SECTION_POINT_SOURCE 4 3 BOUNDARY e Make BOUNDARY_PRESCRIBED_MOTION_SET compatible with CONTROL_REFINE e Fix SET option for BOUNDARY_PWP_TABLE e Fix a bug for BOUNDARY_PAP that could happen in MPP when CVMAS gt 0 e Fix a bug for BOUNDARY_PRESCRIBED_ORIENTATION_RIGID_ANGLES that could be triggered when a high sampling rate is used for LCIDQ e Fix incorrect ID number in d3hsp for BOUNDARY_CYCLIC e Fix output to bndout for BOUNDARY_PRESCRIBED_ORIENTATION_RIGID e Fix incorrect velocity boundary conditions if BOUNDARY_PRESCRIBED_MOTION with dof 141 or I 81 is used with DAMPING GLOBAL e Fix input error when using more than one BOUNDARY_ACOUSTIC_COUPLING keyword card e Fix input error when using BOUNDARY_CYCLIC with the _ID option when e there s only one BOUNDARY_CYCLIC card e Fix input error when using BOUNDARY_THERMAL BULKNODE with number of segments NBNSEG gt 1 e Fix long computation time when there s a large number of velocity boundary conditions BOUNDARY_PRESCRIBED_MOTION Ove Arup and Partners
56. al mechanical problem for SPH variable FRCENG in CONTROL_CONTACT This option applies to all 3D contact types supported by SPH particles Add support of MAT_027 MAT_MOONEY_RIVLIN_RUBBER for SPH For keyword DEFINE_ADAPTIVE_SOLID_TO_SPH add support of explicit SPH thermal solver for the newly generated SPH particles which were converted from solid elements The temperatures of those newly generated SPH particles are mapped from corresponding solid elements Instead of writing all extra history variables for SPH particles to d3plot write only NEIPH extra history variables NEIPH is given in DATABASE EXTENT_BINARY Add support of MAT_THERMAL_ISOTROPIC_TD_LC for SPH with thermal option Add support of MAT_188 MAT_THERMO_ELASTO_VISCOPLASTIC_CREEP for SPH Add support of BOUNDARY_SPH_SYMMETRY for SPH formulation with Lagrangian kernel IFORM 7 or 8 in CONTROL SPH Add support of LOAD_REMOVE_PART for SPH particles which allows users to deactivate the SPH particles of a part 4 17 DISCRETE ELEMENT METHOD Disable particle to particle contact between particles that are bonded via DEFINE DE BOND or DEFINE_DE_HBOND Ove Arup and Partners Limited June 16 June 2014 Page 52 LS DYNA Software Release Notes Version 971 R7 1 1 Add time step calculation for DES bond When DEFINE_DE_BOND is defined skip the contact between two DES nodes if they are bonded Implement DE to surface tied contact DEFINE_DE_TO_
57. allow better evaluation of formability e Improved CONTACT AUTO MOVE e Before Changes the termination time causing problems when several tools needed to be moved e Now Does not change the termination time but changes the current time In this way several tools can be moved without the need to worry about the other tool s move This is especially useful in multi flanging and hemming simulations INTERFACE_BLANKSIZE improvements including adding the options _INITIAL_TRIM and _INITIAL ADAPTIVE This keyword was developed Ove Arup and Partners Limited June 16 June 2014 Page 18 LS DYNA Software Release Notes Version 971 R7 1 1 for blank size development in sheet metal forming Generally for a single forming process only the option DEVELOPMENT is needed and inputs are an initial estimated blank shape a formed blank shape and a target blank shape in either mesh or boundary coordinates Output will be the calculated corrected initial blank shape Initial blank mesh and formed blank mesh can be different e g adaptive For a multi stamping process involving draw trimming and flanging all three options are needed Related commands for blank size estimation are CONTROL_ FORMING_ONESTEP and for trim line development CONTROL_FORMING_UNFLANGING improvements This keyword unfolds flanges of a deformable blank e g flanged or hemmed portions of a sheet metal part onto a rigid tooling mesh using the implicit static sol
58. and changed the wedge element case to use 2 pt integration e Add MAT_ACOUSTIC elform 8 support for pyramid element case using 5 pt integration e Add support to MAT_219 FMAT_CODAM for negative AOPT values which point to coordinate system ID s e Add and reorder extra history variables in MAT_155 per user request Also corrected an apparent problem with history data for optional viscous stress and reduced the memory allocated at the same time e Add support of MAT_PLASTICITY_WITH_DAMAGE_ORTHO_RCDC for the element formulation 13 This element formulation was apparently overlooked when RCDC work was done e Fix the orientation of aniso orthotropic materials that have AOPT 0 when used with 2D plane strain plane stress or axisymmetric elements and the element normal are defined in the z direction We automatically correct the orientation by changing the order of the element connectivity but now we preserve the 1 2 edge so that AOPT 0 works properly Ove Arup and Partners Limited June 16 June 2014 Page 47 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix damage calculation of MAT_221 MAT_ORTHOTROPIC_SIMPLIFIED DAMAGE e The strain calculation used to determine the damage and failure was incorrect causing failure to occur at too small a strain value e Fix SMP and hybrid parallel consistency when using MAT_239 for shells e Fix the use of MAT ADD _ EROSION with thick shell forms 1 and 2 The messages report
59. b tar gz ls dyna_mpp_s_r7_1_1_88920_x64_suse102_pgil05_sse2_platformmpi tar gz SMP Executables Linux ia64 ls dyna_ smp d r7 1 1 x64_redhat59_ifort131 gz ls dyna_smp_s_r7_1_1_x64_redhat59_ifort131 gz ls dyna_ smp_d_r7_1_1_x64_susell_pgil34 gz ls dyna_ smp_s_r7_1_1_x64_susell_pgil34 gz Windowx X64 ls dyna_smp_s_R711_winx64_ifort31 zip ls dyna_ smp_d_R711_winx64_ifort31 zip Ove Arup and Partners Limited June 16 June 2014 Page 67
60. by relative motion of the assembly instead of absolute motion CONTROL_IMPLICIT_BUCKLING Extend Implicit Buckling Feature to allow for Implicit problems using Inertia Relief This involves adding the Power Method as a solution technology for buckling eigenvalue problems Using the power method as an option for buckling problems that are not using inertia relief has been added as well Ove Arup and Partners Limited June 16 June 2014 Page 53 LS DYNA Software Release Notes Version 971 R7 1 1 e Extend Implicit Buckling to allow for intermittent extraction by using negative values of NMODE on CONTROL_IMPLICIT_BUCKLING similar to using negative values of NEIG on CONTROL_IMPLICIT_EIGENVALUE e Extend implicit explicit switching specified on CONTROL_DYNAMIC_RELAXATION to allow explicit simulation for the dynamic relaxation phase and implicit for the transient phase e New implementation for extracting resultant forces due to joints for implicit mechanics e New implementation of extracting resultant forces due to prescribed motion for implicit mechanics e Add support for DEFINE_FUNCTION_TABULATED to implicit mechanics for prescribed motion constraints e Add the treatment of CONTACT_1D to implicit mechanics constraint processing e Correct handling of tied groupable contact with implicit MPP Only e Correct the generation of the header for the d3eigv d3iter and d3mode files for the case of MPP with predecomposition
61. cy is controlled by the new keyword Add new feature DEFINE_DE_INJECTION to allow DEM particle dropping from user defined plane Add new option VOLUME to ELEMENT_DISCRETE_SPHERE This will allow DEM input based on per unit density and use MAT card to get consistent material properties Ove Arup and Partners Limited June 16 June 2014 Page 19 LS DYNA Software Release Notes Version 971 R7 1 1 3 25 ISOGEOMETRIC ELEMENTS Added FORM 4 for ELEMENT_SHELL_NURBS_PATCH Rotational dofs are automatically set at control points at the patch boundaries whereas in the interior of the patch only translational dofs are present This helps for joining multiple nurbs patches at their CO boundaries Disabled FORM 2 and 3 for ELEMENT_SHELL_NURBS_PATCH These formulations are experimental and not fully validated yet Added energy computation for isogeometric shells ELEMENT_SHELL_NURBS_PATCH to matsum Allow isogeometric shells ELEMENT _SHELL_NURBS_PATCH to behave as rigid body MAT_RIGID Fix for degenerated isogeometric shells ELEMENT_SHELL_NURBS_PATCH that share multiple control points at the same physical location which results in a not well defined element normal This caused stability problems and error terminations 3 26 MISCELLANEOUS Fix a bug affecting INCLUDE_TRANSFORM that refers to a transformation DEFINE_TRANSFORMATION involving SCALE Added g as abbreviation for gigawords in specification of me
62. e Convert the computation of modal stresses to use the orthonormal form of the eigenvectors instead of the one scaled for LS PrePost visualization e Fix gathering of forces for MPP Superelement implementation Forces on shared nodes were counted multiple times e Fix gathering of forces for MPP PART_MODES implementation Forces on shared nodes were counted multiple times e Fix implementation of LaGrange Multiplier joint formulation CONSTRAINED_JOINT CONTROL RIGID to properly handle prescribed motion constraints MPP Only e Fix handling of multiple cyclic boundary conditions by implicit mechanics Some dynamic memory was being freed after the first set causing a memory error for subsequent sets SMP Only e Correct the computation of external work for implicit mechanic problems where the external work is coming from prescribed motion on rigid bodies e Correct the time of the output to spcforc for implicit mechanics e Fix the computation of resultant forces for implicit MPP Only e Add synchronization for implicit contact penetration checking used by implicit forming MPP Only e Fix computation of inertia computation for superelements created by CONTROL_IMPLICIT_MODES e Fix guided contact for implicit applications SMP Only e Correct setting of drilling rotation stiffness for quadratic shell and triangle elements formulations 23 amp 24 for implicit eigenvalue problems Ove Arup and Partners Limited June
63. e precision INTERFACE_LINKING when 32ieee yes is used Stability improvements for tied offset contact when TIEDID 1 CONTACT or when groupable CONTROL_MPP_CONTACT_GROUPABLE Ove Arup and Partners Limited June 16 June 2014 Page 56 LS DYNA Software Release Notes Version 971 R7 1 1 e Reduce MPP memory usage on clusters e Add support for ELEMENT_SOURCE_SINK e Fix a bug that could cause a run with a large number of cores to hang when the model includes CONTACT_ENTITY and solid elements e Fix memory error messages after Normal Termination for when using DAMPING_RELATIVE e Add support for MAT_ACOUSTIC in combination with BOUNDARY_USA_SURFACE e Support hybrid with parallel consistency option for shell formulations 23 and e Fix hex assemblies by DEFINE_HEX_SPOTWELD_ASSEMBLY when more than 8 elements are in the assembly There was an error in the decomposition code that caused an error termination e Fix a bug in segment based contact A problem was found when contact occurs between segments on different partitions It was likely that the scale factor used to scale the stiffness for the passed segments was the SFM parameter even if the segment was on the slave side or for single surface contact If SFM SFS in the data then there was no problem e Fix CONTACT_2D_AUTOMATIC when used with mass scaling The problem resulted in a failed run due to segmentation fault e Fix an error that could cause a hang in sleout
64. ed include element ID s Either the element ID s were wrong or a segmentation fault occurred when writing the messages e Make major update of spot weld failure option 9 in MAT_SPOTWELD The nugget fracture calculation at the middle of a 3 sheet weld was incorrectly calculated using the difference of forces and moments e Fix type 6 7 9 and 10 hourglass controls for MAT_076 MAT_GENERAL_VISCOELASTIC These hourglass types all use material constants and for brick elements of MAT_076 a failure to define the shear modulus caused a divide by zero and NaN s e Fix bugs related to material direction definitions for aniso orthotropic materials When BETA is set to zero when using ELEMENT_SHELL BETA or ELEMENT_SOLID_ORTHO we now correctly override the BETA value on the MAT card Also when the MCID option was used on ELEMENT_SHELL_MCID we failed to convert degrees to radians resulting in the wrong material direction e Fix MAT ADD THERMAL EXPANSION for thick shell formulation 5 The strain was totally wrong Also for thick shell formulations 3 and 5 add support for thermal expansion of total Lagrangian materials that use the deformation gradient e Fix stress and strain output of solid elements when CMPFLG 1 and the material type is 2 21 23 40 91 92 128 132 189 266 or 267 and the elements rotate The outputted values were not being rotated so output was in the initial material direction not the current direction e Changes
65. ehaviour when using MAT_FINITE_ELASTIC_STRAIN_PLASTICITY with LSCR defined for the strain rate effects Ove Arup and Partners Limited June 16 June 2014 Page 46 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix incorrect strains when using MAT USER DEFINED MATERIAL MODELS in Material Model Driver e Fix spurious element deletion when using MAT_081 MAT_PLASTICITY_WITH_DAMAGE with MAT ADD EROSION e Fix incorrect failure strains and yield stress when using DEFINE_HAZ PROPERTIES for MAT_PIECEWISE_LINEAR_PLASTICITY_STOCHASTIC shells e Fix incorrect results when using MAT_SOFT_TISSUE MAT_091 with AOPT lt 0 e Fix incorrect behaviour for MAT MODIFIED HONEYCOMB when LCSR 1 e Incorrect error message was displayed for MAT_HYDRAULIC_GAS_DAMPER_DISCRETE_BEAM when the length of the beam is less than then the length of the gas column CO e Fix incorrect behaviour for materials 22 54 55 and 213 when LAMSHT is set to 1 2 or 3 in CONTROL_SHELL and thick shell type 1 or 2 are being used e Check that MAT_ELASTIC_PLASTIC_THERMAL MAT_ 004 input temperature is either increasing or decreasing continuously e Add warning message MSG_SOL 1231 if load curve has negative slope for MAT_024 or MAT_251 e Fix internal energy calculation by correcting uninitialized data with MAT ACOUSTIC elform 14 e Fix MAT_ACOUSTIC elform 8 time step calculation for wedge elements to reliably handle distorted configurations
66. fer to a load curve stress vs stretch ratio or a table stress vs stretch ratio vs normalized strain rate e Expand list of variables for MAT_USER_DEFINED_MATERIAL_MODELS by characteristic element size and element id e Enable MAT_USER_DEFINED_MATERIAL_MODELS to be used with tetrahedron element type 13 New sample routines umat41_t13 and umat41v_t13 show corresponding pressure calculation in the elastic case e Fix for MAT_119 with initial displacements IUR IUS IUT IWR IWS IWT which were erroneously applied twice e Fix for MAT_083 with damage hysteresis TBID lt 0 or HU gt 0 Combination with FMATRX 2 on CONTROL_SOLID default in implicit could lead to uninitialized energies and consequential incorrect unloading behaviour e Fix for MAT_104 For FLAG 1 anisotropic damage it was possible that the load curve number of LCSS was changed internally e Fix for combination of INCLUDE_TRANSFORM with FCTTIM ne 1 0 and material tables with logarithmic strain rates The transformation of these logarithmic strain rates is correct now e Add anew feature to MAT_125 allowing C1 and C2 to be used in calculation of back stress When plastic strain lt 0 5 C1 is used otherwise C2 is used as described in Yoshida s paper e Extend non linear strain path _ NLP FAILURE in MAT_037 to implicit e MAT_173 MAT_MOHR_COULOMB now works in ALE A new option has been added to suppress the tensile limit on hydrostatic stress
67. geometry in dynamic relaxation IDRFLG 2 CONTROL_DYNAMIC_RELAXATION add an option where displacements are not imposed linearly but rather according to a polar coordinate system This option was added to accommodate large rotations e Disable improper restart of dynamic relaxation after adaptive remeshing e Fix input error when using CONTROL_CHECK_SHELL and there are solid elements in the input e Fix a typo that causes seg fault when using CONTROL_ADAPTIVE_CURVE e Fix incorrect response of shells triggered by the following three concurrent conditions DAMPING_PART_STIFENESS is applied to the shells and Damping energy calculation is not turned on in CONTROL_ENERGY and Shell bulk viscosity is turned on in CONTROL_BULK_VISCOSITY e Fix CONTROL_REMESHING settings VF_LOSS MFRAC DT_MIN which were not triggering the remesh at all e Display warning message if more than one CONTROL_DYNAMIC_RELAXATION or CONTROL_DAMPING keywords are used e Suppress two more items during d3hsp printing if NPOPT 1 CONTROL OUTPUT Initial strains and adaptive constraints Ove Arup and Partners Limited June 16 June 2014 Page 34 LS DYNA Software Release Notes Version 971 R7 1 1 e Support NFAILI and NFAIL4 of CONTROL_SHELL in coupled thermal mechanical analysis i e erode distorted elements instead of error termination e PTSCL on CONTROL CONTACT can be used to scale contact force exerted on shell formulations 25 26 27 as
68. gh to account properly for thickness offsets when one segment is thick and the other thin This will change results when SHLEDG 1 or when DEPTH gt 23 for which edge extensions are always turned off e Fix two bugs related to 2D automatic tied contact When the stiffness scaling was determined by part stiffness and dimensions ISTIFF 1 default for implicit the tied contact created NaN s which destroyed the solution Fix for the IGAPCL option to close constraint gaps when R adaptive remeshing is active Gaps were allowed to grow during remeshing e Fix an apparent memory error in the segment based contact surface to surface bucket sort that could result in a memory overrun e Fix DELETE_CONTACT_2DAUTO when the contact interfaces did not have sequential ID numbers that begin with 1 e For CONTACT_2D_TIED and CONTACT_2D_TIED_ONE_WAY enforce penalty constraints if either side of the constraint belongs to a rigid body Ove Arup and Partners Limited June 16 June 2014 Page 32 LS DYNA Software Release Notes Version 971 R7 1 1 e Improve convergence in implicit for CONTACT_TIED_SHELL EDGE_TO_SURFACE_BEAM_ OFFSET e Fix for OPTION 7 of CONTACT AUTOMATIC TIEBREAK The failure equation adopted from the paper of Lemmen and Meijer was wrong This is now corrected in the code and in the manual e Correct out of range forces problem with CONTACT AUTOMATIC TIEBREAK OPTIONS 7 and 10 Only MPP e For CONTACT AUTOM
69. ime is the current solution time The other arguments are a numerical average of the values at the node points defining the element 4 21 MPP Change treatment of two sided CONTACT_FORCE_TRANSDUCER so that proper mass and moment values can be output to the rcforc file Support for non zero birthtime for CONTACT_SINGLE_EDGE Fix d3hsp echo of bucketsort frequency Fix some CONTROL_MPP_DECOMPOSITION keywords Fix adaptivity error termination in the case where processors do not all share the same local directory Add new pfile decomposition region option partsets Takes a list of part sets SET PART from the keyword input and uses them to define a region e g region partsets 102 215 sy 1000 This example would take partsets scale y by 1000 and decompose them and distribute them to all processors Changed the behaviour when a decomposition transformation is used with CONTROL_MPP_DECOMPOSITION_PARTS_DISTRIBUTE CONTROL_MPP_DECOMPOSITION_PARTSET_DISTRIBUTE and CONTROL_MPP_DECOMPOSITION_ARRANGE_PARTS Previously any such regions were distributed without the transformation being applied This has been fixed so that any given transformation applies to these regions also Improvements to velocity initialization in MPP full deck restart CHANGE_VELOCITY_RIGID_BODY now works CHANGE_VELOCITY can be applied to rigid body nodes Rigid body velocity is calculated taking into account those nodal velocities Fix for doubl
70. imited June 16 June 2014 Page 62 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix error when using DEFINE_CURVE_FUNCTION with RCFORC contact function and local coordinate system local is defined e Fix SET_NODE_LIST_GENERATE_INCREMENT which is generating incorrect node list e Fix error message when using SET PART LIST GENERATE COLLECT e Add warning message if MID1 MID2 0 for PART_COMPOSITE e Fix incorrect node set generated when using NODE_TRANSFORM with SET_NODE_GENERAL with BOX option e Add warning message if some of the nodes in SET_NODE referenced by e EDGSET of SECTION SHELL do not belong to seatbelt elements e Fix bug whereby INCLUDE_TRANSFORM created an ELEMENT_DISCRETE and an ELEMENT_BEAM for each ELEMENT_DISCRETE that it transformed No ELEMENT BEAM should be created e Jobid can now be changed in a restart by including jobid on the restart execution line Previously the jobid stored in d3dump could not be overwritten e All dynain output INTERFACE_SPRINGBACK_LSDYNA is disabled for any model in which any nodes or elements are encrypted in the input and suppressed from output via PSETID in DATABASE_BINARY_D3PLOT e Part labels PID can be up to 8 characters in standard format 20 characters in long format Labels for sections SID materials MID equations of state EOSID hourglass IDs HGID and thermal materials TMID can be up to 10 characters in standard f
71. lot PSETID in DATABASE_BINARY_D3PLOT when model contained an SPH part and or a DEM part e Fix bug in reading 64 bit d3eigv d3mode data using a single precision executable Single precision LS DYNA can now read single precision and double precision d3eigv d3mode data e Fix corruption of d3part database by a small restart e DATABASE FSI SENSOR Create sensors at solid faces in 3D and at shell sides in 2D e DATABASE_TRACER If the input deck has EOS_IDEAL_GAS and SECTION ALEID or SECTION_ALE2D write the temperature instead of z coordinate to trhist e DATABASE PROFILE Implement the option DIRA to plot data with curvilinear distributions and the flag UPDLOC to update the profile positions e Fix incorrect stress output to d3plot for the in plane integration points 2 3 4 The stresses were not in the local coordinate system when CMPFLG 1 and Ove Arup and Partners Limited June 16 June 2014 Page 37 LS DYNA Software Release Notes Version 971 R7 1 1 MAXINT lt 0 in FDATABASE EXTENT_BINARY for fully integrated shell type 16 Further to this fix duplicate the stresses amp history variables from in plane intg point 1 to in plane intg points 2 3 4 for shell element formulations with only in plane intg point This is to handle models that use both 1 and 4 in plane integration point shell formulations and which have MAXINT lt 0 e Fix cross section forces output for solid elements when DATABASE_CROSS_SECTION is
72. mentioned material models e Minor modification for DEFINE CONNECTION PROPERTIES PROPRUL 2 thinner weld partner is first partner PROPRUL 3 bottom nodes 1 2 3 4 weld partner is first partner e Add spotweld area to debug output of DEFINE_CONNECTION_PROPERTIES which is activated by CONTROL_DEBUG e Add support of MAT_ADD_EROSION option NUMFIP lt 0 for standard non GISSMO failure criteria Only for shells e Improve implicit convergence of MAT ADD EROSION damage model GISSMO by adding damage scaling 1 D to the tangent stiffness matrix e Provide plastic strain rates tension compression shear biaxial as history variables no 16 17 and 18 for MAT_187 Ove Arup and Partners Limited June 16 June 2014 Page 12 LS DYNA Software Release Notes Version 971 R7 1 1 e Add new variables to user failure routine matusr_24 activated by FAIL lt 0 on MAT 024 and other materials integration point numbers and element id e Add new energy based nonlocal failure criterion for MAT_ADD_EROSION parameters ENGCRT critical energy and RADCRT critical radius after EPSTHIN Total internal energy of elements within a radius RADCRT must exceed ENGCRT for erosion to occur Intended for windshield impact e Add new option to MAT_054 for thin shells Load curves for rate dependent strengths and a rate averaging flag can be defined on new optional card 9 e Add new option for MAT_MUSCLE Input parameter SSP lt 0 can now re
73. mory on execution line e g memory 16g is 16 billion words Suppress non printing characters in COMMENT output Add command line option pgpkey to output the current public PGP key used by LS DYNA The output goes to the screen as well as a file named Istc_pgpkey asc suitable for directly importing into GPG When reading the NAMES file allow a anywhere on a line to indicate there will be a following line not just at the end This was never intended but worked before r73972 and some customers use it that way Check for integer overflow when processing command line arguments and the memory value on the KEYWORD card Added new capability for INTERFACE_LINKING_NODE to scale the displacements of the moving interface Several fixes for internal CASE driver Fix problem of LSDA file names getting mangled in the restart file which was causing file reopening problems during restarts Support for KEYWORD_JOBID with internal CASE driver DAMPING_FREQUENCY_RANGE now works for implicit dynamic solutions An error check has been added to ensure that the timestep is small enough for the damping card to work correctly Added new option DAMPING_FREQUENCY_RANGE_DEFORM to damp only the deformation instead of the global motion Ove Arup and Partners Limited June 16 June 2014 Page 20 LS DYNA Software Release Notes Version 971 R7 1 1 e Fixed bug affecting models containing both DAMPING_GLOBAL and DAMPING_FREQUENCY_RANGE
74. n64_455 ls dyna_hyb_d_r7_1_1_88920_x64_suselll_open64455_avx_platformmpi_sharelib tar gz ls dyna_hyb_d_r7_1_1_88920_x64_suselll_open64455_avx_platformmpi tar gz ls dyna_hyb_s_r7_1_1_88920_x64_susel11_open64455_avx_platformmpi_sharelib tar gz ls dyna_hyb_s_r7_1_1_88920_x64_suse1l11_open64455_avx_platformmpi tar gz ls dyna_mpp_d_r7_1_1 88920_x64 suselll_open64455_avx_platformmpi_sharelib tar gz ls dyna_mpp_d_r7_1_1 88920_x64_susell1_open64455_avx_platformmpi tar gz ls dyna_mpp_s_r7_1_1 88920_x64_susel11_open64455_avx_platformmpi_sharelib tar gz ls dyna_mpp_s_r7_1_1_88920_x64_suse111_open64455_avx_platformmpi tar gz pgi_105 ls dyna_hyb_d_r7_1_1_88920_x64_suse102_pgil05_sse2_platformmpi_sharelib tar gz Is dyna_hyb_d_r7_1_1_88920_x64_suse102_pgil05_sse2_platformmpi tar gz ls dyna_hyb_s_r7_1_1_88920_x64_suse102_pgi105_sse2_platformmpi_sharelib tar gz ls dyna_hyb_s_r7_1_1_88920_x64_suse102_pgi105_sse2_platformmpi tar gz ls dyna_mpp_d_r7_1_1_88920_x64_suse102_pgil05_sse2_mpich126 tar gz ls dyna_mpp_d_r7_1_1_88920_x64_suse102_pgil05_sse2_openmpi165 tar gz ls dyna_mpp_d_r7_1_1_88920_x64_suse102_pgil05_sse2_platformmpi_sharelib tar gz ls dyna_mpp_d_r7_1_1_88920_x64_suse102_pgil05_sse2_platformmpi tar gz ls dyna_mpp_s_r7_1_1_88920_x64_suse102_pgil05_sse2_mpich126 tar ls dyna_mpp_s_r7_1_1_ 88920_x64_suse102_pgil05_sse2_openmpi165 tar gz ls dyna_mpp_s_r7_1_1_88920_x64_suse102_pgi105_sse2_platformmpi_shareli
75. nd thickness of and elastic perfectly plastic coated layer of a fabric that results in a rotational resistance during the simulation e FLDNIPF lt 0 in MAT_190 MAT_FLD_3 PARAMETER_ BARLAT for shell elements means that failure occurs when all integration points within a relative distance of FLDNIPF from the mid surface has reached the fld criterion e A computational welding mechanics MAT_270 MAT_CWM material is available that allows for element birth based on a birth temperature as well as annealing based on an annealing temperature The material is in addition a thermo elasto plastic material with kinematic hardening and temperature dependent properties e MAT_271 MAT_POWDER is a material for manufacturing 1 e compaction and sintering of cemented carbides It is divided into an elastic plastic compaction model that is supposed to be run in a first phase and a viscoelastic sintering model that should be run in a second phase This model is for solid elements e For IHYPER 3 on a MAT_USER_ DEFINED shell material the deformation gradient is calculated from the geometry instead of incremented by the velocity gradient The deformation gradient is also passed to the user defined subroutines in the global system together with a transformation matrix between the global and material frames This allows for freedom in how to deal with the deformation gradient and its transformations in orthotropic layered materials e The Bergstrom
76. ng with structural and thermal LS DYNA solvers Eulerian fixed mesh or moving mesh Either type input with ELEMENT_SOLID cards or using MESH cards Solver Main Features e CESE method Conservation Element Solution Element used for CFD solver Highly accurate shock wave capturing Cavitation model Embedded immersed boundary approach or moving fitting approach for FSI problems Coupled stochastic fuel spray solver See STOCHASTIC keywords Coupling with chemistry See CHEMISTRY keywords solver Solver Applications Non exhaustive Shock wave capturing Shock acoustic wave interaction Cavitating flows Conjugate heat transfer problems Many different kinds of stochastic particle flows e g dust water fuel Chemically reacting flows e g detonating flow supersonic combustion 3 3 New Incompressible CFD ICFD Keyword family ICFD_ The keywords starting with ICFD refer to and control the incompressible CFD solver problem set up Solver Characteristics Implicit Double precision Dynamic memory handling SMP and MPP 2D solver 3D solver Makes use of an automatic volume mesh generator for fluid domain See MESH keywords e Coupling with structural and thermal LS DYNA solvers Solver Main Features Incompressible fluid solver Thermal solver for fluids Free Surface flows Two phase flows Turbulence models Transient or steady state problems Non Newtonian fluids Boussinesq model for convection Ove
77. ning message when more than one LOAD_DENSITY_DEPTH is used Fix the use of LOAD_ BODY with ELEMENT_ MASS MATRIX The applied load was scaled using the nodal mass without the contribution of the added mass matrix so the load was too small Fix error in LOAD_THERMAL_ VARIABLE BEAM SET Fix bug in staged construction keywords LOAD_STIFFEN_PART and LOAD REMOVE PART These did not work with PART COMPOSITE 4 13 SENSOR Add BAGVENTPOP for SENSOR_CONTROL This allows user more flexibility controlling the pop up of the venting hole of AIRBAG HYBRID and AIRBAG_WANG_NEFSKE Add command SENSOR_DEFINE_FUNCTION Up to 15 DEFINE_SENSORs can be referenced in defining a mathematical operation LAYER of SENSOR_DEFINE_ELEMENT can now be an integer T representing the Ith integration point at which the stress strain of the shell or tshell element will be monitored Ove Arup and Partners Limited June 16 June 2014 Page 43 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix a bug when SENSOR_DEFINE_FORCE and SENSOR_CONTROL involves a CONSTRAINED_JOINT e Add control of LOAD MOVING PRESSURE by using SENSOR_CONTROL e Add thick shells to the ETYPE option list of SENSOR_DEFINE_ELEMENT e Fix a bug that could happen when VID of SENSOR_DEFINE_FORCE is larger than 100000 e Fix problem where sensor switch did not turn off when using SENSOR_DEFINE_FORCE with ftype PRESC MOT when running implicit e Fix input err
78. o MAT_244 MAT UHS STEEL e Added implicit support for MAT_244 e Changed the influence of the austenite grain size in Mat244 according to Li et al e Changed the start temperatures to fully follow WATT et al and Li et al e Hardness calculation is now improved when noncontinuous cooling is applied i e tempering e Added temperature dependent Poisson ratio and advanced reaction kinetics e Added new advanced option to describe the thermal expansion coefficients for each phase e Added option to use Curve ID or a Table ID for describing the latent heat generation during phase transormations e Added support for table definition for Youngs modulus Now you can have one temperature dependent curve for each of the 5 phases Bugfixes to MAT_244 Fixed a bug related to thick shell formulation 1 Fixed an error in load curves when using temperature dependent variables Youngs modulus and poissons ratio Fixed a problem when using dynain files Added support for implicit to MAT_188 Ove Arup and Partners Limited June 16 June 2014 Page 11 LS DYNA Software Release Notes Version 971 R7 1 1 e Added material model MAT_273 MAT_CDPM MAT_CONCRETE_DAMAGE_PLASTIC_MO DEL This model is aimed at simulations where failure of concrete structures subjected to dynamic loadings is sought The model is based on effective stress plasticity and has a damage model based on both plastic and elastic strain measures Implemented for solids
79. option SOFT 6 in CONTACT_FORMING_NODES_TO_SURFACE for blank edge and guide pin contact e Add user defined criteria for mesh refinement or coarsening in CONTROL_REFINE_ e Add new contact option that currently only works for MPP SINGLE_SURFACE contact with SOFT 0 or 1 if SRNDE field 4 of optional card E is a 1 then free edges of the contact definition will be rounded WITHOUT extending the segments Rather than having cylindrical caps on the ends of the segments the corners of the squared off thickness are rounded over e Add geometric contact entity type 3 finite cylinder e Add irate 2 to CONTROL_IMPLICIT DYNAMICS to turn off rate effects for both implicit and explicit e Add quadratic 8 node and 6 node shells shell formulations 23 and 24 e Add LOG_LOG_INTERPOLATION option for table defining strain rate effects in MAT_083 MAT_181 and MAT_ 183 e Add automatic generation of null shells for quadratic shell contact PART_DUPLICATE_NULL_OVERLAY e Add beam contact forces to rcforc output DATABASE RCFORC e Add SHIA TO SHL8 option to ELEMENT_SHELL to automatically convert 4 node shells to 8 node quadratic shells e Add 3 node beam element with quadratic interpolation that is tailored for the piping industry It includes 12 degrees of freedom including 6 ovalization degrees of freedom per node for a total of 36 DOF An internal pressure can be given that can stiffen and elongate the pipe ELFORM 14
80. or when using SENSOR_DEFINE_FORCE with INCLUDE_TRANSFORM due to the idtype not being updated as specified with IDDOFF in INCLUDE _TRANSFORM 4 14 MAT e Fix bug in MAT_015 where EROD flag was inactive for VP 1 e Fix spurious deletion of beam elements that could happen when using MAT ADD EROSION with exclusion number EXCL e Fix initialization error for MAT_108 MAT_ORTHO_ELASTIC_ PLASTIC with thick shell formulation 2 e Update MAT_214 MAT_DRY_FABRIC to allow fibres to rotate independently e Enable regularization curve LCREGD of MAT_ADD_EROSION to be used with FLD criterion i e load curve LCFLD Ordinate values major strain will be scaled with the regularization factor e Correct initialization of MAT_CABLE in Staged Construction when starting from a nonzero time Fix a bug for MAT_255 that could happen when Poisson s ratio is input as a function of temperature e Add CONTROL_MPP_MATERIAL_MODEL_DRIVER in order to enable the Material Model Driver for MPP 1 core e Fix incorrect nodal pressure scaling w r t the damage parameter in MAT_081 for type 13 tet elements This caused error termination at start of damage e Add table input of thermal expansion coefficient for MAT_270 Supports temperature dependent curves arranged according to maximum temperature e Add table input of heat capacity for MAT T07 Supports temperature dependent curves arranged according to maximum temperature e Fix various bugs in
81. ormat 20 characters in long format e Create bg switch and kill_by_pid for SMP Both files will be removed at the termination of the run e Increase the overall length of command line to 1000 characters and length of each command line option to 50 characters e Fix input error if SET SHELL GENERAL created shell set of more than 250 shells e Fix incorrect train forces when using RAIL_TRAIN RAIL_ TRACK and the tracks are not in the global XY plane Ove Arup and Partners Limited June 16 June 2014 Page 63 LS DYNA Software Release Notes Version 971 R7 1 1 5 Release notes for LS DYNA version R7 1 1 e Fix softening of MAT_BRITTLE_DAMAGE MAT_096 behaviour caused by typo in r86847 e Fix error in handling of CONTROL_REQUIRE_REVISION e Fix problem with parameter SCLMRR of DEFINE_CONNECTION_PROPERTIES e New encryption keys now in effect There are two a 1024 bit key and a new 2048 bit key The 2048 bit key has a bit of a performance hit during keyword input so should be generally avoided e Fix problem generating RIGIDWALL visualization nodes if beam orientation nodes are also generated e Fix MPP memory overwrite during DATABASE_BINARY_INTFOR initialization that could have caused a segmentation fault e Fix bug and clean up d3thdt output FDATABASE BINARY _D3THDT for SPH e Fix memory problem while using BOUNDARY_SPH_SYMMETRY_PLANE option MPP only e Fix output of INITIAL_STRESS_SHELL to dynain for
82. output initialization when CONTACT_2D_AUTOMATIC was used e Fix rcforc output for CONTACT_2D_AUTOMATIC The contact ID was incorrect e Fix a memory error in MPP segment based eroding contact that occurred with thick shells in the contact definition This error caused a segmentation fault in phase 3 of initialization e Fix bug affecting reading of ELEMENT_SOLID_ORTHO data when there is a mix of 8 and 10 noded elements e Add new pfile options decomp d2r_as_rigid decomp d2ra_as_rigid which cause materials appearing in DEFORMABLE_TO_RIGID and DEFORMABLE_TO_RIGID_AUTOMATIC to have their computational costs set as if they were rigid materials during the decomposition e Modify decomposition related to massless nodes previously they were all dumped on processor 0 under the assumption that there weren t very many of them For topological optimization there can be a LOT of them They are now distributed across all processors on a round robin basis Ove Arup and Partners Limited June 16 June 2014 Page 57 LS DYNA Software Release Notes Version 971 R7 1 1 4 22 FORMING Fix no output to trimcurves ibo when using INTERFACE BLANKESIZE DEVELOPMENT and defining JOBID Add option ISRCOUT to INCLUDE_STAMPED_PART to dump out the transformed source stamp mesh Fix bugs and make improvements for trimming DEFINE_CURVE_TRIM Fix a bug in trimming using single precision solver Fix problem in 2 D trimming
83. ow supported for solid element formulation 13 e For MAT USER DEFINED MATERIAL MODELS the subroutines crvval and tabval can be called with negative curve table id which will extract values from the user input version of the curve or table instead of the internally converted 100 point curve table e Inthe damage initiation and evolution criteria of MAT_ADD_EROSION invoked by IDAM lt 0 add the option Q1 lt 0 for DETYP 0 Here IQ1I is the table ID defining the ufp plastic displacement at failure as a function of triaxiality and damage value i e ufp ufp eta D as opposed to being constant which is the default e In MAT_RHT ONEMPA 6 generates parameters in g cm and microS and ONEMPA 7 generates parameters in mm ms and mg e In MAT_SIMPLIFIED_RUBBER FOAM STOL gt 0 invokes a stability analysis and warning messages are issued if an unstable stretch point is found within a logarithmic strain level of 100 e Fix bug for MAT_190 when FLD level is low element won t fail e Fix a bug in MAT_037 when normal stress is applied e Fix spontaneous failure of elements for MAT_123 MAT_MODIFIED_PIECEWISE_LINEAR_PLASTICITY if EPSTHIN is specified e Fix incorrect von Mises stress for MAT_123 MAT_MODIFIED_PIECEWISE_LINEAR_PLASTICITY after restart when using stress strain input via the 8 data points specified in cards 3 amp 4 of the material cards Ove Arup and Partners Limited June 16 June 2014 Page 50 LS DYNA Software
84. parts from the dynamic relaxation phase e Add new feature for implicit automatic time step control to cooperate with thermal time step control On CONTROL_IMPLICIT_AUTO IAUTO 2 is the same as IAUTO 1 with the extension that the implicit mechanical time step is limited by the active thermal time step e On CONTROL_IMPLICIT_SOLUTION add negative value of MAXREF for implicit mechanics Nonlinear iteration will terminate after IMA XREFI iterations With MAXREF lt 0 convergence is declared with warning Simulation will continue Positive values of MAXREF still cause failure of convergence to be declared leading to either a time step reduction or an error termination e Add CONTROL_IMPLICIT_MODAL_DYNAMIC keywords and features This elevates the modal dynamic features of IMASS 2 on CONTROL_IMPLICIT_DYNAMICS It also adds additional features of damping and mode selection and stress computations e New material model MAT DRY FABRIC MAT 214 which can be used in modelling high strength woven fabrics with transverse orthotropic behaviour e Add ALE COUPLING_NODAL PENALTY penalty based nodal coupling with ALE e Add type 8 ELEMENT_SEATBELT_PRETENSIONER which takes energy time curve instead of pull in or force curve e Add type 9 ELEMENT_SEATBELT_PRETENSIONER for energy based buckle anchor pretensioner Ove Arup and Partners Limited June 16 June 2014 Page 22 LS DYNA Software Release Notes Version 971 R7 1 1 e Add DA
85. pi openmpi mpich msmpi The names of the LS DYNA R7 1 1 executables are listed below The most common executables are highlighted in red and can be downloaded from www oasys software com dyna All other executables can be obtained by contacting the Oasys support team or downloaded directly from the LSTC ftp site Ove Arup and Partners Limited June 16 June 2014 Page 65 Release Notes LS DYNA Software Version 971 R7 1 1 MPP Executables IA64 ls dyna_mpp_d_r7_1_1_88920_ia64_suse102_ifort101_platformmpi tar gz ls dyna_mpp_s_r7_1_1_88920_ia64_suse102_ifort101_platformmpi tar gz Windowx x64 ls dyna_mpp_s_R711_winx64_ifort31_pmpi zip ls dyna_mpp_d_R711_winx64_ifort31_pmpi zip ls dyna_mpp_s_R711_winx64_ifort31_anlmpi zip ls dyna_mpp_d_R711_winx64_ifort31_anlmpi zip ls dyna_mpp_s_R711_winx64_ifort31_msmpi zip ls dyna_mpp_d_R711_winx64_ifort31_msmpi zip ls dyna_mpp_s_R711_winx64_ifort31_impi zip ls dyna_mpp_d_R711_winx64_ifort31_impi zip Linux x86 64 ifort_131 ls dyna_mpp_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_platformmpi_sharelib tar gz ls dyna_mpp_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_platformmpi tar gz ls dyna_mpp_s_r7_1_1_ 88920_x64 _redhat54_ifort131_sse2_platformmpi_sharelib tar gz ls dyna_mpp_s_r7_1_1 88920_x64 _redhat54_ifort131_sse2_platformmpi tar gz ls dyna_hyb_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_in
86. r mass scaling when using MAT_COHESIVE_MIXED_MODE e Output contour plots for MAT_CAZACU_BARLAT MAT_233 when FIT 2 0 e Add materials MAT CONCRETE DAMAGE REL3 MAT_072R3 and MAT CSCM_ CONCRETE MAT_ 159 to Interactive Material Model Driver e Enable MAT_JOHNSON_COOK MAT_015 for shell elements to work with coupled structural thermal analysis e Allow shells used in 2D analysis plane strain amp axisymmetric to be deleted when meeting failure criteria in the constitutive routine for materials 9 26 72 110 111 145 210 221 241 and 272 e Fix input error when using parameter for card 3 fields 4 to 8 of MAT_FLD_3 PARAMETER_BARLAT e Allow MAT_SOIL_AND_FOAM MAT_005 to use positive or negative abscissa values for load curve input of volumetric strains e Fix spurious element deletions when using LAMSHT 3 laminated shell theory on CONTROL_SHELL with shell elements using MAT_NULL MAT_009 e Fix spurious element failures when using MAT_072R3 and DAMPING_PART_STIFENESS in explicit analysis e Fix incorrect DATABASE SSSTAT MASS PROPERTIES and DATABASE_GLSTAT_MASS_ PROPERTIES output to binout e Add warning message to warn user that if more than one user material of the same type i e more than one MAT_041 only parameters from the cards 1 amp 2 first one will be used e Allow load curves for MAT_GURSON MAT_GURSON_JC MAT_GURSON_RCDC in card 6 to be read as integers or real numbers e Fix incorrect b
87. ree enriched FEM MEFEM This element is based on the 4 noded MEFEM element ELFORM 43 SECTION_SOLID Combined with CONSTRAINED_TIED_NODES_FAILURE SET_NODE_LIST and cohesive model this element can be used to model dynamic multiple crack propagation along the element boundaries e Fix the volume calculation for thick shell form 3 Prior to the fix the calculated volume was 2x the correct value This caused a wrong pressure calculation for equation of states 2 3 and 14 which use the element volume in the denominator of the pressure calculation It also caused the mass scaling calculation to use 2x the added mass that was needed for effective mass scaling e Fix mass scaling of thick shell 5 which added 8x the mass that was needed e Enable hyperelastic materials those that use Green s strain to be used with thick shell form 5 Previously use of these materials 2 7 21 23 27 30 31 38 40 112 128 168 and 189 with thick shell 5 has been an input error e New high order tetrahedron CPE3D10 based on Cosserat Point theory can be invoked by specifying element formulation ELFORM 16 and combining this with hourglass formulation IHQ 10 See SECTION SOLID and HOURGLASS e Improve ELEMENT_LANCING Fix mismatch between local and user curve IDs in lancing when interval is defined with a curve Allow multiple curves during lancing Fix mismatch of global and local element IDs in lancing 4 10 FREQUENCY_DOMAIN e Enable
88. rts to avoid problems there is still a conflict in the model e Adda scale factor for scaling the frictional stiffness for contact The parameter is e FRICSF on optional card E and it s only supported for segment based SOFT 2 contact This was motivated by a rubber vs road skidding problem where the friction coefficient had static dynamic and decay parameters defined The growth of the frictional force was too slow so the static coulomb value could not be achieved By scaling the frictional stiffness higher the coulomb value could approach the static value e Add keyword CONTACT_2D_AUTOMATIC_FORCE_TRANSDUCER Like the 3D force transducers it does no contact calculation but only measures the contact forces from other contact definitions When only a slave side is defined the contact force on those segments is measured Currently two surface force transducers are not available e Add options to MAT 058 e Load curves for rate dependent strain values E11C ELIT can be defined on new optional card 9 e Load curves for rate dependent strengths XC XT etc and a rate averaging flag can be defined on new optional card 8 e Abscissa values in above curves are taken to be natural log of strain rate when the first value is negative e Add optional transverse shear damage to MAT_058 e Add material models MAT 261 and MAT_262 e MAT 261 is MAT_LAMINATED_ FRACTURE DAIMLER _PINHO e MAT 262 is MAT_LAMINATED_FRACTURE_
89. running BEM acoustics Extend Kirchhoff acoustic method to MPP e Implement acoustic panel contribution analysis to collocation BEM and dual collocation BEM e Implement panel contribution analysis capability to Rayleigh method asse _PANEL CONTRIBUTION e FREQUENCY _DOMAIN_RESPONSE_SPECTRUM e Fix bug in response spectrum analysis when damping is defined by DEFINE_CURVE e Fix bug in response spectrum analysis when different type input spectra are present e Enable FREQUENCY_DOMAIN_MODE in response spectrum analysis e Extend response spectrum analysis to multiple load spectra cases FREQUENCY_DOMAIN_RANDOM_VIBRATION Fix running random vibration with time history load Fix running random fatigue with sets of elements instead of parts Modify the definition of snlimit Now it is defined on each SN curve Implement biased spacing for output frequencies for random vibration Fix dumping NODOUT_PSD and ELOUT_PSD to binout in MPP Fix running random vibration with base acceleration Fix support of dumping D3FTG database BOUNDARY_ACOUSTIC_MAPPING e Implement boundary acoustic mapping to acoustic BEM in MPP This is enabled only for segment sets at present e Fix running boundary acoustic mapping with multiple field points in MPP 4 11 INITIAL Fix a bug that could be triggered in an adaptive simulation involving INITIAL STRAIN_SHELL Fix incorrect initial velocities if PART INERTIA is used with IRIGID 1 and ICID is defined
90. s Add support for dumping of strain tensor STRFLG eq 1 for isogeometric shells 4FELEMENT_SHELL_NURBS_PATCH via interpolation shells Fix thickness change output istupd ne 0 for ELEMENT_SHELL_NURBS_PATCH amp generalized shells Fix post processing using interpolation shells ELEMENT_SHELL_NURBS_PATCH amp ELEMENT_INTERPOLATION_SHELL Fix writing mass properties to d3hsp file in case of isogeometric elements ELEMENT_SHELL_NURBS_PATCH Correct stress computation in interpolation elements for isogeometric shells ELEMENT_SHELL_NURBS_PATCH Correct writing of control net to igaplot for ELEMENT_SHELL_NURBS_PATCH formulation FORM 4 Improve the critical time step computation for NURBS shells ELEMENT_SHELL_NURBS_PATCH Ove Arup and Partners Limited June 16 June 2014 Page 59 LS DYNA Software Release Notes Version 971 R7 1 1 Correct energy computation reported to glstat for isogeometric shell elements ELEMENT_SHELL_NURBS_PATCH EM ELECTROMAGNETIC SOLVER Change the EM_2DAXI card in order to incorporate the ssidIn and ssidOut which were previously in the EM_BOUNDARY card The use of EM_BOUNDARY is thus no longer necessary for 2D axi simulations Add H field magnetization and relative permeability to d3plot output 4 24 ICFD INCOMPRESSIBLE FLOW SOLVER Fix a bug that was causing the level set interface to duplicate in the presence of embedded shells Now MESH_INTERF and MESH_
91. s now also correct in elastic regions e Fix for failure criterion LCFLD of MAT_ADD_EROSION Reincarnation of an already failed integration point was possible Modify MAT 224 so it uses the temperatures from the thermal solution for a coupled thermal mechanical problem e Switch off rate effects in MAT_123 shells for static implicit analyses or dynamic analyses with IRATE 1 or 2 e Switch off rate effects in MAT_181 and MAT_183 for static implicit analyses or dynamic analyses with IRATE 1 or 2 This was not yet working for the table option e Correct contact stiffness calculation for MAT_106 with LCE gt 0 and LCPR gt 0 e Enable MAT_133 to be used with user defined failure matusr_24 Input parameter USRFAIL 0 off 1 on is on Card 6 7 Column 7 e Modify MAT_ADD_EROSION parameter EPSTHIN EPSTHIN gt 0 individual thinning for each IP from z strain as before EPSTHIN lt 0 averaged thinning strain from element thickness new e Add alternative solution method Brent for MAT_015 and MAT_157 in case standard iteration fails to converge e Fix principal strain calculation for MAT_ADD_EROSION option LCFLD Tensorial shear strain was used but engineering shear strain is correct e Fix for MAT_ADD EROSION parameter MXEPS lt 0 The strain rates in the corresponding load curve were not correctly transformed when INCLUDE_TRANSFORM with FCTTIM ne 1 0 was used e Enable MAT_ 071 to use initial axial force F11 from
92. t it is believed that it should not be a problem to remove this restriction Consequently this limit has been raised to 10 with the caveat that the user should verify the results for NMCOL gt 2 Implemented a structural acoustic mapping scheme BOUNDARY_ACOUSTIC_MAPPING for mapping transient structural nodal velocity to acoustic volume surface nodes This is useful if the structure finite element mesh and the acoustic boundary finite element mesh are mismatched 3 7 CONTACT CONTACT FORMING ONE WAY SURFACE TO SURACE ORTHO_ FRICTION can now be defined by part set IDs when supplemented by DEFINE FRICTION ORIENTATION Segment sets with orientation per DEFINE FRICTION ORIENTATION are generated automatically Contact force of CONTACT ENTITY is now available in intfor DATABASE BINARY INTFOR Fix for CONTACT AUTOMATIC TIEBREAK options 9 and 11 element size dependent failure stresses NFLS SFLS lt 0 were not working correctly in MPP Fix for CONTACT_AUTOMATIC_ _TIEBREAK options 7 and 10 in MPP no more unexpected high forces and subsequent instability CONTACT_FORCE_TRANSDUCER_PENALTY will now accept node sets for both the slave and master sides which should allow them to work correctly for eroding materials BOTH sides should use node sets or neither Ove Arup and Partners Limited June 16 June 2014 Page 6 LS DYNA Software Release Notes Version 971 R7 1 1 e Added option to create a backup penalty
93. telmpi 322 tar gz ls dyna_hyb_d_r7_1_1 88920_x64 _redhat54_ifort131_sse2_openmpi165 tar gz ls dyna_hyb_d_r7_1_1_88920_x64 redhat54_ifort131_sse2_platformmpi_sharelib tar gz ls dyna_hyb_d_r7_1_1 88920_x64 _redhat54_ifort131_sse2_platformmpi tar gz ls dyna_hyb_s_r7_1_1_88920_x64 redhat54_ifort131_sse2_intelmpi 322 tar gz ls dyna_hyb_s_r7_1_1_88920_x64 redhat54_ifort131_sse2_openmpi165 tar gz ls dyna_hyb_s_r7_1_1_88920_x64 _redhat54_ifort131_sse2_platformmpi_sharelib tar gz ls dyna_hyb_s_r7_1_1_88920_x64_redhat54_ifort131_sse2_platformmpi tar gz ls dyna_mpp_d_r7_1_1_88920_x64 redhat54_ifort131_sse2_fjmpi_ompi tar gz ls dyna_mpp_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_intelmpi 322 tar gz ls dyna_mpp_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_mpich126 tar gz ls dyna_mpp_d_r7_1_1 88920_x64_redhat54_ifort131_sse2_mpich2 121 tar gz ls dyna_mpp_d_r7_1_1_88920_x64_redhat54_ifort131_sse2_openmpi165 tar gz ls dyna_mpp_s_r7_1_1_88920_x64_redhat54_ifort131_sse2_fjmpi_ompi tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_sse2_intelmpi 322 tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_sse2_mpich126 tar gz ls dyna_mpp_s_r7_1_1_ 88920_x64_redhat54_ifort131_sse2_mpich2 121 tar gz ls dyna_mpp_s_r7_1_1 88920_x64_redhat54_ifort131_sse2_openmpi165 tar gz ifort_131_avx2 ls dyna_mpp_d_r7 ls dyna_mpp_d_r7 ls dyna_mpp_s_r7 ls dyna_mpp_s_r7 ls dyna_hyb_d_r7 ls dyna_hyb_d_r7 ls dyna_hyb_d_r7 1_88920_x64_redhat54_ifort131_avx2_platformmpi_sharelib tar gz 1
94. thod provides more accurate intersection information when Q2TRI 0 Fix incorrect CONTACT_AUTOMATIC TIEBREAK interface forces if OPTION 4 and PARAM and viscous damping coefficient VDC gt 0 0 Fix full deck restart for CONTACT AUTOMATIC ONE WAY SURFACE TO SURFACE TIE BREAK for SMP only Fix full deck restart for CONTACT AUTOMATIC ONE WAY SURFACE TO SURFACE TIE BREAK Fix intermittent seg fault when using CONTACT_FORMING_ SMOOTH in SMP with ncpu gt 1 Fix incorrect numbering of last intfor file DATABASE_BINARY_INTFOR Fix input error when SSTYP 6 in CONTACT and there are gt 9999 parts in the part set defined in SSID Fix seg fault problem on some platforms when using CONTACT_NODES_TO_SURFACE with FS 2 0 and THKOPT ne 0 in optional card B or SALTHK ne 0 in CONTROL_CONTACT Fix seg fault problem when using CONTACT_AUTOMATIC_SURFACE_TO_SURFACE_ORTHO_FRICTI ON Fix contact instability when contact thickness OPTT is specified in PART_CONTACT Contact instability may occur if the merged nodes of shells have large differences in shell thicknesses Ove Arup and Partners Limited June 16 June 2014 Page 31 LS DYNA Software Release Notes Version 971 R7 1 1 e Add warning message INI 364 that erosion has been enabled when CONTACT ERODING is used e Fix contact thickness for beam types 9 11 and 12 e Add support for birth time for CONTACT_2D_AUTOMATIC_TIED e Fix MPP rcforc output for CONTACT_2D_A
95. tion of EFG in an implicit analysis no longer requires a special license 3 20 IMPLICIT Fix sign for resultant moments for joints using implicit jntforc Fix numerous problems with stress recovery for mode shapes especially for intermittent eigenvalue analysis Fix output of forces pertaining to prescribed motion bndout for implicit Ove Arup and Partners Limited June 16 June 2014 Page 16 LS DYNA Software Release Notes Version 971 R7 1 1 Fix INTERFACE LINKING problem for SMP implicit Fix implicit output of spcforc for CONSTRAINED_COORDINATE Fix implicit BOUNDARY_SPC with birth time and death time Add SENSOR_CONTROL for prescribed motion constraints in implicit Update INTERFACE LINKING NODE in implicit to catch up with explicit including adding scaling factors Add support for DATABASE_RCFORC_MOMENT for implicit Enhance Iterative solvers for Implicit Mechanics Add after the first implicit time step the output of projected cpu and wall clock times This was already in place for explicit Also echo the termination time Fix an issue with implicit storing the resultant forces for CONSTRAINED_SPOTWELD 3 21 THERMAL Add variable MXDMP in CONTROL_ THERMAL SOLVER to write thermal conductance matrix and right hand side every MXDMP time steps Add keyword CONTROL_THERMAL_EIGENVALUE to calculate eigenvalue s of each thermal conductance matrix New thermal material MAT THERMAL ORTHOTROPIC_TD_L
96. to MAT_244 MAT_UHS_STEEL The phase change kinetics in MAT_244 for solids and shells has been rewritten This is a major change and there will be a difference when comparing new phase fractions with old results The differences will be small at fast cooling rates The code changes affect results at slow cooling rates the most Fix bug in calculating cooling phase change latent heat for solids and shells Add latent heat calculation to heating phase change kinetics Fix bug in HEAT 2 option for shells to allow both heating and cooling Fix bug when using a dynain file for re start e Fix stability problems in MAT_058 using higher Poisson s ratio together with CONTROL_BULK_VISCOSITY for shells e Fix for damage evolution in fibre compression for MAT_261 MAT_LAMINATED_FRACTURE_DAIMLER_PINHO Ove Arup and Partners Limited June 16 June 2014 Page 48 LS DYNA Software Release Notes Version 971 R7 1 1 e Improve handling of nonphysical material data input for MAT_262 by setting reasonable default fracture toughnesses A warning is printed if default parameters are used e Put classical equivalent plastic strain von Mises type to plastic strain field of MAT_187 e Simplify usage of solid element tet type 13 SECTION SOLID ELFORM 13 with user materials MAT_USER_ DEFINED MATERIAL MODELS Subroutines umat41_t13 and umat41v_t13 of version R7 0 are obsolete e Triaxiality on hisvar 12 of MAT_224 shell elements i
97. upling between two SPH parts and thereby reduce the relative velocity between the parts Added BOUNDARY_CONVECTION and BOUNDARY_RADIATION for explicit SPH thermal solver Fix bug for SPH symmetric boundary conditions CONTROL_SPH 3 19 EFG CONTROL_REMESHING_EFG e Add eroding failed surface elements and reconstructing surface in EFG adaptivity e Adda control parameter for monotonic mesh resizing in EFG adaptivity e Add searching and correcting self penetration for adaptive parts in 3D tetrahedron remeshing e Enhance 3D axisymmetric remeshing CONTROL_REMESHING with 6 and 8 node elements e Use RMIN RMAX along with SEGANG to determine element size e Remove the restriction that the reference point of computational model has to be at original point 0 0 0 e Rewrite the searching algorithm for identifying the feature lines of cross sections in order to provide more stable remeshing results e Improve rigid body motion in EFG shell type 41 Support EFG pressure smoothing in solid type 42 for MAT_ELASTIC_VISCOPLASTIC_THERMAL Add visco effect for implicit EFG solid type 42 Add new EFG solid type 43 called Meshfree Enriched FEM MEFEM for both implicit and explicit This element formulation is able to relieve the volumetric locking for nearly incompressible material eg rubber and performs strain smoothing across elements with common faces EFG shell adaptivity no longer requires a special license Applica
98. used with the PLANE option The normal direction calculation was broken This caused wrong forces as well e Remove DATABASE BINARY_XTFILE since it is obsolete 4 9 ELEMENT e Fix a bug that caused an immediate crash when shell seatbelt elements ELEMENT_SEATBELT were used in the same model with the BIRTH option of AIRBAG REFERENCE GEOMETRY e Fix and enhance ELEMENT_ BEAM PULLEY Correct swapping when node topology of beam1 and beam2 does not form a sequential chain Increase accuracy of some internal variables Account for damping stress correctly when used with MAT_ELASTIC truss elements e Fix a bug for ELEMENT_SEATBELT_SLIPRING that occurred when node topology of sbl and sb2 does not form a sequential chain For example sb1 has node 100 200 as Ist and 2nd connectivity nodes and sb2 has nodes 101 200 e Fix a bug for one way slipring that prevents slipping along the designed direction e Issue error message when seatbelt is initially inside an undefined retractor e Allow SBRNID of 2D slipring and retractor to have up to 8 digits it was previously limited to 7 digits e Correct deforc output for discrete elements whose initial offset is negative that could sometimes be caused by unintentional nonzero discrete element length e Fix bug in linear beam formulation 13 for frequency analysis e When using PART_AVERAGED for truss elements beam formulation 3 calculate the time step based on the total length of
99. using IORTHO 1 for MAT_USER_DEFINED_MATERIAL_MODELS and shell element type 13 e Fix incorrect strains output for shell type 13 and using MAT USER DEFINED MATERIAL MODELS e Fix keyword input error if using MAT CONCRETE DAMAGE PLASTIC MODEL instead of equivalent commands MAT CDPM or MAT 273 e Output warning message KEY 4954 if all data are zero for MAT_WINFRITH_CONCRETE_REINFORCEMENT e Add two more kinematic hardening terms for MAT_DAMAGE _3 MAT_153 c2 amp gamma2 e Include MAT_091 092 MAT SOFT TISSUE in tet solid formulation type 10 e Fix input error for MAT_CONCRETE_DAMAGE_REL3 MAT_72R3 when using large integer value for load curve LCRATE which was read in as real number Change input format for LCRATE from real to integer Ove Arup and Partners Limited June 16 June 2014 Page 45 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix incorrect results when using MAT_GENERAL_VISCOELASTIC with the WLF shift functions for solid elements e Fix incorrect stresses for MAT_COMPOSITE_DAMAGE MAT_022 when INTOUT ALL in DATABASE EXTENT_BINARY for solid elements e Fix time step for MAT_72R3 for 2D analysis e Increase the number of significant digits of the cofficients mu and alpha to 20 in d3hsp for MAT_OGDEN_RUBBER Higher precision is needed if the coefficients are used directly for the material input to achieve the same results as the experimental data e Fix incorrect added mass fo
100. ver using the implicit static solver It is typically used in trim line mapping during a draw die development process The roots of the flanges or hemmed edges are automatically processed based on a user input of a distance tolerance between the flanges hemmed edges and rigid tool It includes the ability to handle a vertical flange wall Other keywords related to blank size development are CONTROL_FORMING_ONESTEP and INTERFACE_BLANKSIZE_DEVELOPMENT Added keyword CONTROL_FORMING_OUTPUT which allows control of d3plot output by specifying distances to tooling home It works with automatic position of stamping tools using CONTROL_FORMING_AUTOPOSITION_PARAMETER Added the LOCAL_SMOOTH option to INTERFACE_COMPENSATION_NEW which features smoothing of a tool s local area mesh which could otherwise become distorted due to e bad coarse mesh of the original tool surface e tooling pairs for example flanging post and flanging steel do not maintain a constant gap e several compensation iterations This new option also allows for multiple regions to be smoothed Local areas are defined by SET_LIST_NODE_SMOOTH 3 24 DISCRETE ELEMENT METHOD Added output to rcfore for DEFINE_DE_TO_SURFACE_COUPLING Implement traction surface for DEFINE_DE_TO_SURFACE_COUPLING Add keyword DATABASE_BINARY_DEMFOR with command line option dem dem_int_force This will turn on the DEM interface force file for DEM coupling option The output frequen
101. ver the 3D elements INITIAL_ALE_MAPPING In a 3D to 3D mapping INITIAL_ALE_ MAPPING map the relative displacements for the penalty coupling in CONSTRAINED_LAGRANGE_IN_SOLID The xy files associated with DATABASE_ ALE MAT are now created when sense switches sw1 sw2 quit or stop are issued ALE ESSENTIAL BOUNDARY is available in 2D DATABASE_FSI is available for 2D MPP only ALE ESSENTIAL BOUNDARY implemented to apply slip only velocity BC along ALE mesh surface CONTROL_ ALE flag INIJWL 2 option added to balance initial pressure state between ALE Soil and HE 3 18 SPH Include SPH element ELEMENT_SPH in time step report Time step and internal energy of 2D axisymmetric SPH elements are calculated in a new way more consistent with the viscosity force calculation Only apply viscosity force to x and y components of 2D axisymmetric SPH element not on hoop component MAXV in CONTROL_SPH can be defined as a negative number to turn off velocity checking Ove Arup and Partners Limited June 16 June 2014 Page 15 LS DYNA Software Release Notes Version 971 R7 1 1 Improve calculation of 2D axisymmetric SPH contact force in DEFINE_SPH_TO_SPH_COUPLING Added the following material models for SPH particles MAT_004 MAT_ELASTIC_PLASTIC_THERMAL 3D only MAT_106 MAT_ELASTIC_VISCOPLASTIC_THERMAL Added a new parameter DFACT for DEFINE_SPH_TO_SPH_COUPLING DFACT invokes a viscous term to damp the co
102. when the trimming curve point spacing is too dense Allow nodal temperature to be considered in trimming Allow nodal temperature and CONTROL CHECK SHELL to be used together Fix a missed trimming problem for an outer trim loop CONTROL_FORMING_OUTPUT Allow NTIMES to be zero support birth and death time support scale factor in curve definition Add a new option INTFOR to CONTROL_FORMING_OUTPUT to control the output frequency of the INTFOR database Add new features instant and progressive lancing in ELEMENT_LANCING for sheet metal lancing simulation Fix adaptive issues in lancing Allow multiple curve intersections during lancing Allow multi lancing locations Set adaptive level to be 1 around the lancing curve Allow lancing boundary to be a closed loop Improve adaptive constraint treatment during lancing Fix problems for symmetry boundary conditions in springback compensation Add a new keyword CONTROL_FORMING_INITIAL_THICKNESS Add a new option for springback compensation INCLUDE_COMPENSATION_ORIGINAL_TOOLS Add a new keyword INTERFACE_COMPENSATION_NEW_PART_CHANGE Fix a bug in springback compensation with local smoothing in INTERFACE COMPENSATION LOCAL SMOOTH Add a new keyword DEFINE_CURVE_BOX_ADAPTIVITY to provide better control of mesh refinement along two sides of the curve Fix a bug in CONTROL FORMING TIPPING where one stress tensor output was incorrect Fix a bug in CO
103. with CONTROL_IMPLICIT_DYNAMICS by activating the DEATH parameter 4 20 THERMAL Add constraint based thermal nodal coupling for CONSTRAINED_LAGRANGE_IN_SOLID HMIN lt 0 turns it on Fix incorrect thermal time step when using iauto lt 0 load curve to determine time step size for CONTROL_IMPLICIT_AUTO together with CASE Add FRCENG 2 on CONTROL_CONTACT keyword FRCENG 1 convert contact frictional energy to heat FRCENG 2 do not convert contact frictional energy to heat Fix bug that prevented thermal solver from correctly handling CONTACT_ENTITY input Add effect of thermal time scaling TSF in CONTROL_THERMAL_SOLVER to 2D contact Fix problem with INITIAL_TEMPERATURE keyword Now bottom and top fictitious nodes belonging to thick thermal shell elements are given the Ove Arup and Partners Limited June 16 June 2014 Page 55 LS DYNA Software Release Notes Version 971 R7 1 1 temperature initial condition of the mid plane node when nsid 0 NSID 0 means initialize all nodes Fix bug that prevented thermal code from correctly handling RIGIDWALL_GEOMETRIC input Fix bug that prevented Ctrl C from interrupting the thermal solver Now Ctrl C will allow the user to enter a sense switch The parameter TGRLC thermal generation rate curve number in MAT_THERMAL option can now reference DEFINE_FUNCTION The function arguments are heat_generation function x y z vx vy vz temperature time The argument t
104. written to message files for mortar contact see also CONTACT_ Good for debugging implicit models not available for explicit e Fix spurious error termination when there is encrypted input and shell adaptivity is invoked e Fix skinny long triangle in CONTROL_CHECK_SHELL 4 8 DATABASE e Change the default scale factor for binary file sizes back to 70 This value can be changed using x on the execution line In version R7 0 0 the default value of x is 4096 and that sometimes leads to difficulty in postprocessing owing to the large size of the d3plot file s e Remove static linking for 12a as many systems do not have the required static libraries Ove Arup and Partners Limited June 16 June 2014 Page 35 LS DYNA Software Release Notes Version 971 R7 1 1 e Fix incorrect strains output to d3plot for thick shells if STRFLG 1 EPSFLG 2 and ENGFLG 2 in DATABASE_EXTENT_BINARY e Fix jntforc forces DATABASE_JNTFORC which were output in wrong coordinate system for rigid bodies when LMF 1 in CONTROL_RIGID e Add IELOUT in CONTROL_REFINE to handle how child element data is handled in elout 4FDATABASE_HISTORY_SOLID and DATABASE HISTORY SHELL Child element data are stored if IELOUT 1 or if refinement is set to occur only during initialization e Prevent d3drlf file from being overwritten when dynamic relaxation is enabled and adaptivity is turned on e Fix incorrect stresses to eloutdet when MAXINT lt 0 in DA

LS-DYNA R7.1.1 Release notes

Contents

Download Pdf Manuals

Related Search

Related Contents