Home

LS-DYNA KEYWORD USER`S MANUAL

1. Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID Part identifier for the surface element This is the surface with sym metry LS DYNA R7 0 6 3 MESH MESH MESH_EMBEDSHELL MESH_EMBEDSHELL Purpose Define surfaces that the mesher will embed inside the volume mesh These surfaces will have no thickness and will conform to the rest of the volume mesh having matching nodes on the interface Card 1 1 2 3 4 5 6 7 8 Variable VOLID Type I Default none Cards 2 3 4 the next card terminates the input 1 2 3 4 5 6 7 8 Variable PIDI PID2 PID3 PID4 PIDS PID6 PID7 PID8 Type I I I I I I I I Default none none none none none none none None VARIABLE DESCRIPTION VOLID ID assigned to the new volume in the keyword MESH VOLUME The surface mesh size will be applied to this volume PIDn Part IDs for the surface elements that will be embedded in the volume mesh 6 4 MESH LS DYNA R7 0 MESH_INTERF MESH MESH_INTERF Purpose Define the surfaces that will be used by the mesher to specify fluid interfaces in multi fluid simulations Card 1 1 2 3 4 5 6 7 8 Variable VOLID Type I Default none Cards 2 3 4 the next card terminates the input 1 2 3 4 5 6 7 8
2. Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface where a free slip boundary condition is applied 5 4 ICFD LS DYNA R7 0 ICFD BOUNDARY FSI ICFD ICFD BOUNDARY FSI Purpose This keyword defines which fluid surfaces will be considered in contact with the solid surfaces for fluid structure interaction FSI analysis This keyword should not be defined if ICFD_CONTROL FSI is not defined Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface in contact with the solid domain LS DYNA R7 0 5 5 ICFD ICFD ICFD BOUNDARY NONSLIP ICFD BOUNDARY NONSLIP Purpose Specify the fluid boundary with a non slip boundary condition Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface where a non slip boundary condition is applied 5 6 ICFD LS DYNA R7 0 ICFD BOUNDARY PRESCRIBED MOVEMESH CFD ICFD BOUNDARY PRESCRIBED MOVEMESH Purpose Allows the node of a fluid surface to translate in certain directions using an ALE ap proach This is useful in piston type applications or can also be used in certain cases to avoid big mesh deformation Card 1 1 2 3 4 5 6 7 8 Variable PID dofx dofy dofz Type
3. Card 1 1 2 3 4 5 6 7 8 Variable PARTID CHEMID COMPID Type I I Default none none Card 2 1 2 3 4 5 6 7 8 Variable UIC VIC WIC RHOIC PIC TIC HIC Type F F F F F F F Default none none none none VARIABLE DESCRIPTION PARTID Identifier of the CESE part on which to initialize CHEMID Identifier of chemistry control card to use COMPID Identifier of chemical composition to use UIC X component of the fluid velocity VIC Y component of the fluid velocity WIC Z component of the fluid velocity RHOIC Initial fluid density PIC Initial fluid pressure LS DYNA R7 0 2 29 CESE CESE VARIABLE TIC HIC CESE_INITIAL_CHEMISTRY_PART DESCRIPTION Initial fluid temperature Initial fluid enthalpy However when CHEMID refers to a ZND 1 step reaction card this is the progressive variable degree of combustion 2 30 CESE LS DYNA R7 0 CESE INITIAL CHEMISTRY SET CESE CESE INITIAL CHEMISTRY SET Purpose Initializes the chemistry and fluid state in every element of the specified element set in the CESE mesh that has not already been initialized by CESE INITIAL CHEMISTRY _ ELEMENT cards This is only used when chemistry is being solved with the CESE solver Card 1 1 2 3 4 5 6 7 8 Variable CHEMID COMPID Type I I Default none none Card 2 1 2 3 4 5 6 7 8 Variable Type Default VARIABLE DESCRIPTION SETID Identifier of the CESE element set to initialize
4. Card 1 1 2 3 4 5 6 7 8 GAMMA THETA C3 none none none Card 2 1 2 3 4 5 6 7 8 Variable C4 K EXPON LGTUNIT TIMUNIT ADJUST Type F F I F F I I Default none none none none none none none In the following UUS stands for User Units System and BUS for Burgess Units VARIABLE DESCRIPTION EOSID ID of the EM_EOS specified by an EM_MAT card Reference specific volume Vo UUS GAMMAO Reference Gruneisen value Yo no units THETA Reference melting temperature in eV BUS LF Latent heat of fusion in kJoule mol BUS CI C1 constant BUS 4 16 EM LS DYNA R7 0 EOS BURGESS DESCRIPTION C2 C2 constant no units C3 C3 constant no units C4 C4 constant no units K Parameter k no units EXPON Exponent in equations 2 see remarks LGTUNIT Length units for UUS relative to meter i e 1 3 if UUS in mm TIMUNIT Time units for UUS relative to seconds TEMUNIT Temperature units EQ 1 temperature in Celsius EQ 2 temperature in Kelvins ADJUST Conductivity modification EQ 0 default the conductivity is given by the Burgess formula EQ 1 The conductivity is adjusted so that it is equal to the conduc tivity defined in EM_MAT card ma at room temperature 0 Opurgess 8 Omat OBurgess room Remarks The Burgess model gives the electrical resistivity vs temperature and density for the solid phase li
5. Card 5 1 Variable FILE4 Type A VARIABLE DESCRIPTION ID Identifier for this one dimensional detonation solution XYZD Position of the detonation front in the DETDIR direction DETDIR Detonation propagation direction 1 gt X 2 gt Y 3 gt Z FILE1 Name of the LSDA file containing the one dimensional solution FILE2 Name of the file containing the Chemkin compatible input FILE3 Name of the file containing the chemistry thermodynamics database FILE4 Name of the file containing the chemistry transport properties database 3 6 CHEMISTRY LS DYNA R7 0 CHEMISTRY_CONTROL_FULL CHEMISTRY_CONTROL_FULL CHEMISTRY Purpose Computes the full chemistry specified by a chemkin chemistry model This card can be used for general purpose chemical reaction calculations Card 1 1 2 4 5 6 7 8 Variable ID ERRLIM Type I F Default none none Card 2 1 Variable FILEI Type A Card 3 1 Variable FILE2 Type A Card 4 1 Variable FILE3 Type A VARIABLE DESCRIPTION ID Identifier for this full chemistry calculation LS DYNA R7 0 3 7 CHEMISTRY CHEMISTRY CHEMISTRY CONTROL FULL VARIABLE DESCRIPTION ERRLIM Error tolerance for the full chemistry calculation FILE1 Name of the file containing the Chemkin compatible input FILE2 Name of the file containing the chemistry thermodynamics database FILE3 Name of the file containing the chem
6. CHEMID Identifier of chemistry control card to use COMPID Identifier of chemical composition to use UIC X component of the fluid velocity VIC Y component of the fluid velocity WIC Z component of the fluid velocity RHOIC Initial fluid density PIC Initial fluid pressure LS DYNA R7 0 2 31 CESE CESE VARIABLE TIC HIC CESE INITIAL CHEMISTRY SET DESCRIPTION Initial fluid temperature Initial fluid enthalpy However when CHEMID refers to a ZND 1 step reaction card this is the progressive variable degree of combustion 2 32 CESE LS DYNA R7 0 CESE MAT GAS CESE CESE MAT GAS Purpose Define the fluid gas properties in a viscous flow for the CESE solver Card 1 1 2 3 4 5 6 7 8 i Default 0 72 1 VARIABLE DESCRIPTION MID Material identifier C1 C2 Two coefficients in the Sutherland s formula for viscosity i e L C T T C5 where C and are constants for a given gas For example for air at moderate temperatures C 71 458x10 kg m s 110 4 Prnd The Prandtl Number used to determine the coefficient of thermal con ductivity It is approximately constant for most gases For air at stand ard conditions Prnd 0 72 Remarks 1 These material coefficients are only used to calculate the viscosity in viscous flows so for invisid flows this material card is not needed 2 As with other solvers in LS DYNA the user is responsible for unit consis
7. ttt ttt ttt ttt 5 26 DEITA 5 27 TC BID AA saat tha BE 5 28 ICED PART tee re ce uot unicum M M LAE 5 30 VOL OPTION ee st dide e tado li tena dp 5 31 WICH SECTION ta he E e EAE coe das 5 33 PIC HSE NOI ec tef he rcs 5 34 EMESE ues tie st Mtt ID tM eI 6 1 IMESH BA eas ete DEL eS 6 2 IMESH BES Moe is ih Dict ala No tan 6 3 oh oes bale oh de Baal 6 4 AMESA See eye rcd i Sep 6 5 IMESID SIZE 6 6 MESH SIZE SHAPE sssesssesssssscsscsssssssssssssssssvevsssssssscssessssssensnsssssnsvensesssseseesessssesesensusnsnsnensentsee 6 7 MESH_SURFACE_ELEMENT ste ato Seca cies 6 9 MESH SURFACE NODE ioris chet i ial epe Er I MEE 6 10 VES EV OUI ies ten aed ub puc x AR EE ala ala DEL ans 6 11 MESH VOLUME EBEMENT 6 12 MESH VOLUME NODE ettet ttt ttti ttd 6 13 MESH_VOLUME_PART ode tenebat detected e s cuu AM DEL 6 14 STOCHASTIC A PM uL LT ML KDE KL MEME NM ME MEN DL 7 1 STOCHASTIC SPRAY PARTICLES ettet titt 72 SSTOCHASTICOIBX sic nd donc De ome ise M set as i 7 5 eae ec xd MI A nU uM E A C 8 1 ESO DOMAINS 8 2 ESO POINT SET anrea aaa
8. fe Aan tet saith tide te tie studeas 3 3 CEMISIRY CONTROL enone dale enh tame Ohi Ald a 3 5 CHBEMISTR Y CONTROL FUEL hacen etu NAR gu ROI 3 7 OHEMISTRY CONIRULCSZND G Aue e onto 3 9 CHEMISTRY_DET_INITIATION S ttd ba dete fa b RUE 3 10 CHEMISTRY_PATH cesta e uU ls e cu E 3 12 N MN VETTER NER TES 4 1 ioi rp o 4 3 SEM BOUNDARY cero uaque a AM mda M M D EE 4 4 nc 4 5 LS DYNA R7 0 0 3 TABLE OF CONTENTS FEM CIRC iatis tense cole e eere ee EM Rc MU cesta 4 8 EM CONDUCT nann uM EU nn e c pe 4 9 VENE CONTACT RESISTANCE 2 itte td s ub d obs t ch o DAE 4 10 4 12 4 13 CONTROL TIMESTEB o oe tbt b a olio oot 4 14 EOS 8 ttt aisnean natt 4 16 SEMCEOSCNIBADON suene ira D EORR dE Med dota ME awed 4 19 EM_EOS_PERMEABILITY baia P c bt eu nod ois 4 22 EOS TABULATED Ia S ER DOR 4 23 KEMCEXTERNAL FIELD 26 litis Malls vse 4 24 asd tfo etras anatomia uates 4 25 sues cochon tpe dashes toc tae ane 4 26 Gone teet bas A eeu mu Mt DDR la E 4 27 EM ROTATION AXIS sette ttt 4 29 oa dde eub bantE S toe Mc d SUE 4 30 FEM SODVER BEMM AS iie dea Steele Me 4 32 FEMUSOBVE
9. impose flow variable s velocity density pressure temperature etc at the centroid of the element connected with this boundary segment OPTION SEGMENT or at the centroid of the elements connected with each segment in a set of boundary segments OPTION SET or at the centroid of the elements created by MESH VOLUME cards that are connected with each face in a surface part PID defined with MESH SURFACE ELEMENT cards OPTION PART For the PART option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SURFPRT IDCOMP Type I I Default none none For the SET option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SSID IDCOMP Type I I Default none none 2 8 CESE LS DYNA R7 0 CESE_BOUNDARY_ PRESCRIBED For the SEGMENT option define the following card CESE Card 1 1 2 3 4 5 6 7 8 Variable N3 N4 IDCOMP Type I I I I I Default none none none none none Card 2 1 2 3 4 5 6 7 8 Variable LC U LC V LCW LCRHO LCP LC T Type I I I I I I Remarks 12 1 2 1 2 1 2 1 2 1 2 Card 3 1 2 3 4 5 6 7 8 Variable SF U SF V SF W SFRHO SFP SF Type F F F F F F Default 1 0 1 0 1 0 1 0 1 0 1 0 Remarks 1 1 1 1 1 1 VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards SS
10. 0 5 21 ICFD ICFD ICFD_CONTROL_TIME ICFD_CONTROL_TURBULENCE Purpose This keyword enables the user to modify the default values for the turbulence model Card 1 1 2 3 4 5 6 7 8 Variable TMOD Type I Default 0 This card is optional and may be used if TMOD 1 Card 2 1 2 3 4 5 Variable Cel Ce2 oe ok Cu Type F F F F F Default 1 44 1 92 13 1 0 0 09 This card is optional and may be used if TMOD 2 Card 2 1 2 3 4 5 Variable Cs Type F Default 0 2 VARIABLE DESCRIPTION TMOD Indicates what turbulence model will be used EQ 0 Turbulence model based on a variational multiscale approach is used by default EQ 1 RANS k e approach 5 22 ICFD LS DYNA R7 0 ICFD CONTROL TURBULENCE CFD VARIABLE DESCRIPTION EQ 2 LES Smagorinsky sub grid scale model Cel Ce2 ce ck k e model constants Cu Cs Smagorinsky constant LS DYNA R7 0 5 23 ICFD ICFD ICFD DATABASE DRAG ICFD DATABASE DRAG Purpose This keyword enables the computation of drag forces over given parts of the model If multiple keywords are given the forces over the PID surfaces are given in separate files and are also added and output in a separate file Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default None VARIABLE DESCRIPTION PID Part ID of the surface where the drag force will be
11. 6 7 8 Variable PID LCID SF DEATH BIRTH Type I I F F F Default none none 1 1 E 28 0 0 VARIABLE DESCRIPTION PID PID for a fluid surface LCID Load curve ID to describe the temperature value versus time see DEFINE CURVE SF Load curve scale factor default 1 0 DEATH Time at which the imposed temperature is removed EQ 0 0 default set to 1028 BIRTH Time at which the imposed temperature is activated starting from the initial abscissa value of the curve LS DYNA R7 0 5 11 ICFD ICFD ICFD CONTROL ADAPT ICFD CONTROL ADAPT Purpose This keyword will activate the adaptive mesh refinement feature The solver will use an a posteriori error estimator to compute a new mesh size bounded by the user to satisfy a maxi mum perceptual global error Card 1 1 2 3 4 5 6 7 8 Variable MINH MAXH ERR MTH NIT Type F F F I I Default none none none 0 0 VARIABLE DESCRIPTION MINH Minimum mesh size allowed to the mesh generator The resulting mesh will not have an element smaller than MINH even if the minimum size does not satisfy the maximum error MAXH Maximum mesh size ERR Maximum perceptual error allowed in the whole domain MTH Specify if the mesh size is computed based on function error or gradient error EQ 0 Function error EQ 1 Gradient error NIT Number of iterations before a re meshing is forced Default forces a re meshing at every timestep 5 12 ICFD LS DYNA R7 0 ICFD CONTROL ADAP
12. CONJ HEAT ICFD BOUNDARY FREESLIP ICFD BOUNDARY FSI ICFD BOUNDARY NONSLIP ICFD BOUNDARY PRESCRIBED VEL ICFD BOUNDARY PRESCRIBED ICFD BOUNDARY PRESCRIBED TEMP ICFD CONTROL ADAPT ICFD CONTROL ADAPT SZIE ICFD CONTROL FSI ICFD CONTROL MESH ICFD CONTROL MESH MOV ICFD CONTROL OUTPUT ICFD CONTROL PARTITION ICFD CONTROL SPLIT ICFD CONTROL SURFMESH ICFD CONTROL TIME ICFD CONTROL TURBULENCE ICFD DATABASE AVERAGE ICFD DATABASE DRAG ICFD DEFINE POINT ICFD INITIAL ICFD MAT ICFD PART ICFD PART VOL ICFD SECTION ICFD SET NODE LIST LS DYNA R7 0 5 ICFD ICFD ICFD An additional option _TITLE may be appended to all ICFD keywords If this option is used then an 80 character string is read as a title from the first card of that keyword s input At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks 5 2 ICFD LS DYNA R7 0 ICFD BOUNDARY CONJ HEAT ICFD ICFD BOUNDARY CONJ HEAT Purpose Specify which boundary of the fluid domain will exchange heat with the solid Card 1 1 2 3 4 5 6 7 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface in contact with the solid LS DYNA R7 0 5 3 ICFD ICFD ICFD BOUNDARY FREESLIP ICFD BOUNDARY FREESLIP Purpose Specify the fluid boundary with free slip boundary condition
13. Define the type of EOS EQ 1 Permeability defined by a function of B curve H B p EQ 2 Permeability defined by a B function of H curve B uH LCID Load curve ID 4 22 EM LS DYNA R7 0 EOS TABULATEDI EM EM EOS TABULATEDI Purpose Define the electrical conductivity as a function of temperature by using a load curve Card 1 1 2 3 4 5 6 7 8 EOSID LCID Type I I Default VARIABLE DESCRIPTION EOSID ID of the EM_EOS LCID Load curve ID Remarks 1 The load curve describes the electrical conductivity ordinate vs the temperature abscissa The user needs to make sure the temperature and the electrical conductivity given by the load curve are in the correct units Also it is advised to give some bounds to the load curve con ductivities at very low and very high temperatures to avoid bad extrapolations of the con ductivity if the temperature gets out of the load curve bounds LS DYNA R7 0 4 23 EM EXTERNAL FIELD EXTERNAL FIELD Purpose Define the components of a time dependent exterior field uniform in space applied on the conducting parts Card 1 1 2 3 4 5 6 7 8 FIELDID FTYPE FDEF Default 0 VARIABLE DESCRIPTION FIELDID External Field ID FTYPE Field type EQ 1 Magnetic field EQ 2 Electric field not available yet FDEF Field defined by EQ 1 Load Curves LCID X Y Z Load curve ID defining the X Y Z component of the field function of
14. Nz Remarks DESCRIPTION Segment set ID Node ID s defining a segment Load curve ID to define this solid wall boundary movement lt BLANK gt velocity vector of the solid wall If LCID 0 it is defined by Vx Vy Vz itself Otherwise it will be defined by both of the load curve and Vx Vy Vz Nx Ny Nz are not used in this case ROTAT x y amp z coordinates of a point on the rotating axis Unit vector of the rotating axis for the 2D case this is not used The rotating frequency Hz is given by the load curve 1 In this solid wall condition SBC the boundary movement can only be in the tangential di rection of the wall and should not affect the fluid domain size and mesh during the calcula tion otherwise an FSI or moving mesh solver should be used Also this moving SBC only affects viscous flows no slip BC 2 If LCID 0 and Vx Vy Vz 0 0 default this will be a regular solid wall BC 3 For rotating SBC LCID gt 0 must be used to define the rotating speed frequency Hz Also in the 2D case Nx Ny Nz does not need to be defined because it is not needed LS DYNA R7 0 2 15 CESE CESE CESE_CONTROL_LIMITER CESE CONTROL LIMITER Purpose Sets some stability parameters used in the CESE scheme for this CESE compressible flow solver Card 1 1 2 3 4 5 6 7 8 Variable Type Default Remarks VARIABLE DESCRIPTION IDLMT Set the stability limiter option See CESE
15. PRE Purpose Impose a fluid pressure on the boundary Card 1 1 2 3 4 5 6 7 8 Variable PID LCID SF DEATH BIRTH Type I I F F F Default none none 1 1 E 28 0 0 VARIABLE DESCRIPTION PID PID for a fluid surface LCID Load curve ID to describe the pressure value versus time see DEFINE CURVE SF Load curve scale factor default 1 0 DEATH Time at which the imposed motion constraint is removed EQ 0 0 default set to 1028 BIRTH Time at which the imposed pressure is activated starting from the initial abscissa value of the curve LS DYNA R7 0 5 9 ICFD ICFD ICFD BOUNDARY FLUX TEMP ICFD BOUNDARY FLUX TEMP Purpose Impose a temperature flux on the boundary Card 1 1 2 3 4 5 6 7 8 Variable PID LCID SF DEATH BIRTH Type I I F F F Default none none 1 1 E 28 0 0 VARIABLE DESCRIPTION PID PID for a fluid surface LCID Load curve ID to describe the temperature flux value versus time see DEFINE CURVE SF Load curve scale factor default 1 0 DEATH Time at which the imposed motion constraint is removed EQ 0 0 default set to 1028 BIRTH Time at which the imposed pressure is activated starting from the initial abscissa value of the curve 5 10 ICFD LS DYNA R7 0 ICFD BOUNDARY PRESCRIBED TEMP ICFD ICFD_ BOUNDARY PRESCRIBED TEMP Purpose Impose a fluid temperature on the boundary Card 1 1 2 3 4 5
16. computed Remarks 1 The file name for this database is icfdragi for instantaneous drag and icfdraga for the drag computed using average values of pressure and velocities 5 24 ICFD LS DYNA R7 0 ICFD DATABASE AVERAGE CFD ICFD DATABASE AVERAGE Purpose This keyword enables the computation of average variable values at given time inter vals Card 1 1 2 3 4 5 6 7 8 Variable DT Type F Default None VARIABLE DESCRIPTION DT Over each DT time interval an average of the different fluid variables will be calculated and then reset when moving to the next DT interval Remarks 1 The file name for this database is icfdavg dat with the different averaged variable values copied in a ASCII format LS DYNA R7 0 5 25 ICFD ICFD ICFD DEFINE POINT ICFD DEFINE POINT Purpose This keyword defines a point in space that could be used for multiple purposes Card 1 1 2 3 4 5 6 7 8 Variable POID x Y Z Type I F F F Default none none none None VARIABLE DESCRIPTION POID Point ID X x coordinate for the point Y y coordinate for the point Z z coordinate for the point 5 26 ICFD LS DYNA R7 0 ICFD INITIAL ICFD ICFD_INITIAL Purpose Simple initialization of velocity and temperature within a volume Card 1 1 2 3 4 5 6 7 8 Variable PID Vx Vy Vz T Type I F F F F Default none none none No
17. generation of the volume mesh This zone does not need to be entirely defined in the volume mesh Card 1 1 2 3 4 5 6 f 8 Variable Sname Type C Default none if Sname box Cards 2 1 2 3 4 5 6 7 8 Variable MSIZE PminX PminY PminZ PmaxX PmaxY PmaxZ Type F F F F F F F Default none none none none none none none if Sname sphere Cards 2 1 2 3 4 6 7 8 Variable MSIZE Radius CenterX CenterY CenterZ Type F F F F F Default none none none none none LS DYNA R7 0 6 7 MESH MESH if Sname cylinder MESH_SIZE_SHAPE Cards 2 1 2 3 4 5 6 8 Variable MSIZE Radius PminX PminY PminZ PmaxX PmaxY PmaxZ Type F F F F F F F Default none none none none none none none VARIABLE DESCRIPTION Sname Shape name Possibilities include box cylinder and sphere MSIZE Mesh size that needs to be applied in the zone of the shape defined by Sname Pmin X Y Z X Y AZ for the point of minimum coordinates Pmax X Y Z X Y AZ for the point of maximum coordinates Center X Y Z Coordinates of the sphere center in cases where Sname is Sphere Radius Radius of the sphere if Sname is Sphere or of the cross section disk if Sname is Cylinder 6 8 MESH LS DYNA R7 0 MESH_SURFACE_ELEMENT MESH MESH SURFACE ELEMENT Purpose Specify a set of surface elements quadrilateral or tr
18. proval of LSTC requests to reproduce the contents hereof should be sent to sales g lstc com TABLE OF CONTENTS TABLE OF CONTENTS INTRODUCTION A dod oe rR au npe MALE a ra MU ME MD ADAM DM 1 1 M H 2 1 CESE BOUNDARY AXISYMMETRIC 2 2 CESE_BOUNDARY_FSI_OPTION e nene anbtt is mti I ih dtd ndi 2 4 CESE BOUNDARY NON REFLECTIVE OPTION 2 6 CESE BOUNDARY PRESCRIBED OPTION ene 2 8 CESE BOUNDARY REFLECTIVE OPTION eene 2 11 CESE BOUNDARY SOLID WALL OPTIONI OPTIOND2 2 13 CESE CONTROL unum 2 16 SCESE CONTROL MESH e t tct cu d 2 17 CONTROL SOLVER oce ecco a ta EM EI am ae E 2 18 2 20 CESE EOS HOMOG EQUILIB 2 21 GESE EOS DEAL GAS 2 22 2 23 CESE INIVIAE CHEMISTRY consectetue eo atta te pa Tad tede 2 25 CESE INITIAL CHEMISTRY ELEMENT ccccssssssssssseeccsssssesecvessesesseecensccesscenrecsessseseces 2 27 CESE INITIAL CHEMISTRY PART ei ERE E ees 2 29 CHEMISTRY a hid tetto 2 31 2 33 PCESE PAR ee ates cele 2 34 FOHEMISDRY redis ced ic al E ea Qe MER Ea 3 1 SCHEMISTRY COMPOSI LOI tt saa tach te aed E LE 3 2 SOBMISTRYSCONTROE UD oi
19. saturated vapor Remarks 1 Once a cavitation EOS is used the cavitation flow solver will be triggered 2 In this homogeneous equilibrium cavitation model a barotropic equation of state is used This model can be used in small scale amp high speed cavitation flows and it is not good for large scale low speed cavitation calculations LS DYNA R7 0 2 21 CESE CESE CESE EOS IDEAL GAS CESE EOS IDEAL GAS Purpose Define the coefficients Cv and Cp in the equation of state for an ideal gas in the CESE fluid solver Card 1 1 2 3 4 5 6 7 8 Variable EOSID Cv Cp Type I F F Default none 717 5 1004 5 VARIABLE DESCRIPTION EOSID Equation of state identifier Specific heat at constant volume Cp Specific heat at constant pressure Remarks 1 As with other solvers in LS DYNA the user is responsible for unit consistency For example if a user wants to use dimensionless variables Cv amp Cp above also should be replaced by the corresponding dimensionless ones 2 22 CESE LS DYNA R7 0 CESE_INITIAL CESE CESE_INITIAL Purpose Specify constant initial conditions ICs for flow variables at the centroid of each fluid element Card 1 1 2 3 4 5 6 7 8 i i Default 1 225 0 0 VARIABLE DESCRIPTION U V W x y z velocity components respectively RHO density p P pressure P T temperature T Remarks 1 Usually only two of p P amp T are needed to be specified be
20. theory manual EQ 0 limiter format 1 Re weighting EQ 1 limiter format 2 Relaxing ALFA Re weighting coefficient See CESE theory manual BETA Numerical viscosity control coefficient See CESE theory manual EPSR Stability control coefficient See CESE theory manual Remarks 1 020 0 larger values give more stability but less accuracy Usually 0 2 0 or 4 0 will be enough for normal shock problems 2 0 0 lt B lt 1 0 larger values give more stability For problems with shock waves B 1 0 is rec ommended 3 20 0 larger values give more stability but less accuracy 2 16 CESE LS DYNA R7 0 CESE CONTROL MESH MOV CESE CESE CONTROL MESH MOV Purpose For moving mesh CESE this keyword is used to choose the type of algorithm to be used for calculating mesh movement Card 1 1 2 3 4 5 6 7 8 Variable MMSH LIM_ITER RELTOL ABSTOL Type I I F F Default 1 100 1 0e 3 1 0e 3 VARIABLE DESCRIPTION MMSH Mesh motion selector EQ 1 mesh moves using an implicit ball vertex spring method EQ 9 the IDW scheme is used to move the mesh LIM ITER Maximum number of linear solver iterations for the ball vertex linear system RELTOL Relative tolerance to use as a stopping criterion for the iterative linear solver conjugate gradient solver with diagonal scaling preconditioner ABSTOL Absolute tolerance measure for the size of mesh displacement changes to use as a stopping criterion fo
21. time 4 24 EM LS DYNA R7 0 MAT 001 EM MAT 001 Purpose Define the electromagnetic material type and properties for a material whose permea bility equals the free space permeability Card 1 1 2 3 4 5 6 7 8 MID MTYPE SIGMA EOSID i Default none none VARIABLE DESCRIPTION MID Material ID refers to MID in the PART card MTYPE Defines the electromagnetism type of the material EQ 0 Air or vacuum EQ 1 Insulator material these materials have the same electromag netism behavior as 0 EQ 2 Conductor carrying a source In these conductors the eddy cur rent problem is solved which gives the actual current density Typi cally this would correspond to the coil EQ 4 Conductor not connected to any current or voltage source where the Eddy current problem is solved Typically this would cor respond to the workpiece SIGMA Initial electrical conductivity of the material EOSID ID of the EOS to be used for the electrical conductivity see EM EOS card LS DYNA R7 0 4 25 EM 002 002 Purpose Define an electromagnetic material type and properties whose permeability is different than the free space permeability Card 1 1 2 3 4 2 6 7 8 MID MTYPE SIGMA i Default none VARIABLE DESCRIPTION MID Material ID refers to MID in the PART card MTYPE Defines the electromagnetism type of the material EQ 0 Air or vacuum EQ 1 Insulator
22. using an implicit coupling when the FSI coupling is strong In addition to being able to handle free surface flows there is also a bi phasic flow capability that involves modeling using a conservative Lagrangian interface tracking technique Basic turbulence models are also supported This solver is the first in LS DYNA to make use of a new volume mesher that takes surface meshes bounding the fluid domain as input keywords In addition LS DYNA R7 0 1 1 INTRODUCTION INTRODUCTION INTRODUCTION during the time advancement of the incompressible flow the solution is adaptively re meshed as an automatic feature of the solver Another important feature of the mesher is the ability to create boundary layer meshes These anisotropic meshes become a crucial part of the model when shear stresses are to be calculated near fluid walls The ICFD solver is also coupled to the solid thermal solver using a monolithic approach for conjugate heat transfer problems The third solver is an electromagnetics EM solver This module solves the Maxwell equations in the Eddy current induction diffusion approximation This is suitable for cases where the propagation of electromagnetic waves in air or vacuum can be considered as instantaneous Therefore the wave propagation is not solved The main applications are Magnetic Metal Form ing bending or welding induced heating ring expansions and so forth The EM module allows the introduction of a sou
23. 0 170 0 465 0 233 0 330 C 1 145 1 191 1 178 1 226 1 210 0 4133 EXPON 1 1 1 1 1 0 0 700 0 672 0 673 0 670 0 638 0 089 k 0 964 0 910 1 08 1 0 878 2 Table 4 1 4 18 EM LS DYNA R7 0 EOS MEADON EM EM EOS MEADON Purpose Define the parameters for a Meadon model giving the electrical conductivity as a function of the temperature and the density see T J Burgess Electrical resistivity model of metals 4 International Conference on Megagauss Magnetic Field Generation and Related Topics Santa Fe NM USA 1986 Card 1 1 2 3 4 5 6 7 8 i i Card 2 1 2 3 4 5 6 7 8 LGTUNIT TIMUNIT ADJUST In the following UUS stands for User Units System and BUS for Burgess Units none VARIABLE DESCRIPTION EOSID ID of the EM_EOS Cl C1 constant BUS C2 C2 constant no units C3 C3 constant no units TEMUNIT Temperature units temperature in Celsius 2 temperature in Kelvins LS DYNA R7 0 4 19 EM EOS MEADON VARIABLE DESCRIPTION Reference specific volume 005 Reference Gruneisen value Yo no units EXPON Exponent in equations 2 LGTUNIT Length units for UUS relative to meter i e 1 e 3 if UUS in mm TIMUNIT Time units for UUS relative to seconds ADJUST EQ 0 default the conductivity is given by the Burgess formula EQ 1 The conductivity is adjusted so that it is equal to the conductivity defined in the EM MAT card at roo
24. 1 Basic assembly steps EQ 2 Basic assembly steps percentage completed final statistics EQ 3 Basic assembly steps percentage completed statistics at each percentage of completion SOLS Level of solver output on the screen 0 No output EQ 1 Global information at each FEM iteration EQ 2 Detailed information at each FEM iteration SOLF Level of solver output to the messag file 0 No output EQ 1 Global information at each FEM iteration EQ 2 Detailed information at each FEM iteration MESH Controls the output of the mesh data to the d3hsp file 0 No mesh output EQ 1 Mesh info is written to the d3hsp file MEMORY Controls the output of information about the memory used by the EM solve to the messag file LS DYNA R7 0 4 27 EM EM OUTPUT VARIABLE DESCRIPTION EQ 0 no memory information written EQ 1 memory information written TIMING Controls the output of information about the time spent in the different parts of the EM solver to the messag file EQ 0 no timing information written EQ 1 timing information written 4 28 EM LS DYNA R7 0 ROTATION AXIS EM EM ROTATION AXIS Purpose Define a rotation axis for the EM solver This is used with the 2D axisymmetric solver The axis is defined by a point and a direction Card 1 1 2 3 4 5 6 7 8 Variable ZP XD Type F F VARIABLE DESCRIPTION NP X Y Z coordinate of the point ND X Y Z coordinate of direction of the axi
25. 14 LS DYNA R7 0 ICFD_CONTROL_MESH ICFD ICFD_CONTROL_MESH Purpose This keyword modifies default values for the automatic volume mesh generation Only used in 3D cases Card 1 1 2 3 4 5 6 7 8 Variable MGSF Type F Default 1 41 VARIABLE DESCRIPTION MGSF Mesh Growth Scale Factor Specifies the maximum mesh size that the volume mesher is allowed to use when generating the volume mesh based on the mesh surface element sizes defined in MESH SURFACE ELEMENT Values between 1 and 2 are allowed Values closer to 1 will result in a finer volume mesh 1 means the volume mesh is not allowed to be coarser than the element size from the closest surface meshes and val ues closer to 2 will result in a coarser volume mesh 2 means the volume can use elements as much as twice as coarse as those from the closest surface mesh LS DYNA R7 0 5 15 ICFD ICFD ICFD CONTROL MESH ICFD CONTROL MESH MOV Purpose With this keyword the user can choose the type of algorithm for mesh movement Card 1 1 2 3 4 5 6 7 8 Variable MMSH LIM ITER RELTOL ABSTOL Type I I F F Default 2 100 1 0e 3 5 0e 4 VARIABLE DESCRIPTION MMSH Mesh motion selector EQ 1 mesh moves based on the distance to moving walls EQ 2 mesh moves by solving a linear elasticity problem using the element sizes as stiffness EQ 3 mesh uses a Laplacian smoothing with stiffne
26. BLE DESCRIPTION TSTYPE Time Step type EQ 1 constant time step given in 5 EQ 2 time step vs time given by a load curve specified in LCID EQ 3 automatic time step computation depending on the solver type This time step is then multiplied by FACTOR DTCONST Constant value for the time step for TSTYPE 1 LCID Load curve ID giving the time step vs time for TSTYPE 2 FACTOR Multiplicative factor applied to the time step for TSTYPE 3 Remarks 1 For an eddy current solver the time step is based on the diffusion equation for the magnetic field elg bodas x Vx It is computed as the minimal elemental diffusion time step over the elements For a given element the elemental diffusion time step is given as _ le dt 2D where D the diffusion coefficient D is element electrical conductivity 4 14 EM LS DYNA R7 0 EM CONTROL TIMESTEP EM Lig is the permeability of free space le is the minimal edge length of the element minimal size of the element LS DYNA R7 0 4 15 EM EOS BURGESS EOS BURGESS Purpose Define parameters a Burgess model giving electrical conductivity as func tion of the temperature and the density see T J Burgess Electrical resistivity model of metals 4 International Conference on Mega gauss Magnetic Field Generation and Related Topics Santa Fe NM USA 1986
27. HOoxi FaceActe FaceAfilm CIRCID Type I I F F F F I Default none none none none none none none none VARIABLE DESCRIPTION CRID Resistive contact ID CONTID EM contact ID defined in EM CONTACT RHOprobe Probe resistivity Pprob RHOsub Substrate resistivity Psup RHOoxi Film resistivity FaceActe Scale factor on the constriction area when calculating the constriction resistance If negative the factor is time dependent defined by the load curve absolute value FaceActe FaceAfilm Scale factor on the constriction area when calculating the film resistance If negative the factor is time dependent defined by the load curve abso lute value FaceAfilm CIRCID When defined the contact resistance will be added to the corresponding circuit ID total resistance and taken into account in the circuit equations Remarks 1 The contact resistance calculation is based on the book by Ragmar Holm s Electric Con tacts Pprob Psub Reonstriction e 4 faceActe x ContactArea 4 10 EM LS DYNA R7 0 EM_CONTACT_RESISTANCE EM Poxy R Ium faceAfilm x ContactArea Reontact Reonstriction LS DYNA R7 0 4 11 EM EM_CONTROL EM CONTROL Purpose Enable the EM solver and set its options Card 1 1 2 3 4 5 6 7 8 Variable EMSOL NUMLS DTINIT DTMAX T INIT T END NCYLFE NCYLBE Type I I F F F F I I Default 0 25 none none 0 0 E
28. I I I I Default none 1 1 1 VARIABLE DESCRIPTION PID PID for a fluid surface dofx dofy dofz Degrees of freedom in the X Y and Z directions EQ 0 degree of freedom left free Surface nodes can translate in the chosen direction EQ 1 prescribed degree of freedom Surface nodes are blocked LS DYNA R7 0 5 7 ICFD ICFD ICFD BOUNDARY PRESCRIBED VEL ICFD BOUNDARY PRESCRIBED VEL Purpose Impose the fluid velocity on the boundary Card 1 1 2 3 4 5 6 7 8 Variable PID DOF VAD LCID SF VID DEATH BIRTH Type I I I I F I F F Default none None 1 none 1 0 1 E 28 0 0 VARIABLE DESCRIPTION PID PID for a fluid surface DOF Applicable degrees of freedom EQ 1 x degree of freedom EQ 2 y degree of freedom EQ 3 z degree of freedom VAD Velocity flag EQ 1 linear velocity EQ 2 angular velocity EQ 3 parabolic velocity profile LCID Load curve ID used to describe motion value versus time see DEFINE CURVE DEFINE CURVE FUNCTION or DEFINE FUNCTION See BIRTH below SF Load curve scale factor default 1 0 VID Point ID for angular velocity application point see ICFD DEFINE POINT DEATH Time at which the imposed motion constraint is removed EQ 0 0 default set to 1028 BIRTH Time at which the imposed motion constraint is activated starting from the initial abscissa value of the curve 5 8 LS DYNA R7 0 ICFD BOUNDARY PRESCRIBED PRE ICFD ICFD BOUNDARY PRESCRIBED
29. ID Segment set ID NI Node ID s defining a segment IDCOMP For inflow boundaries in problems involving chemical reacting flows LS DYNA R7 0 2 9 CESE CESE VARIABLE LC_U LC_V LC_W LC RHO LC P LC T SF_U SF_V SF_W SF_RHO SF P SF T Remarks CESE BOUNDARY PRESCRIBED DESCRIPTION the chemical mixture of the fluid entering the domain as defined with a CHEMISTRY COMPOSITION card Load curve ID to describe the x component of the velocity versus time see DEFINE CURVE Load curve ID to describe the y component of the velocity versus time Load curve ID to describe the z component of the velocity versus time Load curve ID to describe the density versus time Load curve ID to describe the pressure versus time Load curve ID to describe the temperature versus time Scale factor for LC U default 1 0 Scale factor for LC V default 1 0 Scale factor for LC default 1 0 Scale factor for LC RHO default 1 0 Scale factor for LC P default 1 0 Scale factor for LC T default 1 0 1 On each centroid or set of centroids the variables x vel y vel z vel p P T that are given values must be consistent and make the model well posed 1 e be such that the solution of the model exists is unique and physical 2 If any of the load curves are 0 the corresponding variable will take the constant value of the corresponding scale factor For instance if LC 0 then the
30. IPTION DIR Directory path to add to the search set 3 12 CHEMISTRY LS DYNA R7 0 The keyword cards provide input for a new electromagnetism module for solving 3D eddy current inductive heating or resistive heating problems coupled with mechanical and thermal solvers Typical applications include magnetic metal forming and welding A boundary element method in the air is coupled to finite elements in the conductor in order to avoid meshing the air EM 2DAXI EM BOUNDARY EM CIRCUIT EM CIRCUIT ROGO EM CONTACT EM CONTACT RESISTANCE EM CONTROL EM CONTROL CONTACT EM CONTROL TIMESTEP EM EOS BURGESS EM EOS MEADON EM EOS PERMEABILITY EM EOS TABULATEDI EM EXTERNAL FIELD EM MAT 001 EM MAT 002 EM OUTPUT EM ROTATION AXIS EM SOLVER BEM EM SOLVER BEMMAT EM SOLVER FEM EM SOLVER FEMBEM LS DYNA R7 0 4 1 EM An additional option TITLE may appended to all keywords If this option is used then an 80 character string is read as a title from the first card of that keyword s input At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks 4 2 EM LS DYNA R7 0 20 2DAXI Purpose Sets up electromagnetism solver as 2D axisymmetric instead of 3D on given part in order to save computational time as well as memory The electromagnetism is solved in 2D
31. LL 6 12 MESH LS DYNA R7 0 MESH VOLUME NODE MESH MESH VOLUME NODE Purpose Define a node and its coordinates This keyword is only used in cases where the fluid Volume Mesh is provided by the user and is not automatically generated It serves the same pur pose as the NODE keyword for solid mechanics Only Tets are supported Card 1 1 2 3 4 5 6 7 8 Variable NID X Y Z H Type I F F F F Default none 0 0 0 0 VARIABLE DESCRIPTION NID Node ID A unique number with respect to the other volume nodes X x coordinate Y y coordinate Z z coordinate LS DYNA R7 0 6 13 MESH MESH MESH VOLUME PART MESH VOLUME PART Purpose Associate a volume part number created by a MESH VOLUME card with the part number of a part card from a selected solver designated by the SOLVER field 1 Variable Type Default VOLPRT VARIABLE VOLPRT SOLPRT SOLVER SOLPRT SOLVER DESCRIPTION Part ID of a volume part created by a MESH_VOLUME card Part ID of a part created using SOLVER s part card Name of a solver using a mesh created with MESH cards 6 14 MESH LS DYNA R7 0 STOCHASTIC STOCHASTIC STOCHASTIC The keyword STOCHASTIC is used to describe the particles and numerical details for solving a set of stochastic PDEs Currently there are two types of stochastic PDE models in the code a model of embedded part
32. LS DYNA KEYWORD USER S MANUAL VOLUME III Multi Physics Solvers February 2013 Version R7 0 LIVERMORE SOFTWARE TECHNOLOGY CORPORATION LSTC Corporate Address Livermore Software Technology Corporation P O Box 712 Livermore California 94551 0712 Support Addresses Livermore Software Technology Corporation 7374 Las Positas Road Livermore California 94551 Tel 925 449 2500 Fax 925 449 2507 Email sales g lstc com Website www lstc com Disclaimer Livermore Software Technology Corporation 1740 West Big Beaver Road Suite 100 Troy Michigan 48084 Tel 248 649 4728 Fax 248 649 6328 Copyright O 1992 2013 Livermore Software Technology Corporation All Rights Reserved LS DYNA LS OPT and LS PrePost are registered trademarks of Livermore Software Tech nology Corporation in the United States other trademarks product names and brand names belong to their respective owners LSTC reserves the right to modify the material contained within this manual without prior no tice The information and examples included herein are for illustrative purposes only and are not in tended to be exhaustive or all inclusive LSTC assumes no liability or responsibility whatsoever for any direct of indirect damages or inaccuracies of any type or nature that could be deemed to have resulted from the use of this manual Any reproduction in whole or in part of this manual is prohibited without the prior written ap
33. N SIDVOUT PARTID Type I I I I Default none none none none VARIABLE DESCRIPTION CIRCID Circuit ID CIRCTYP Circuit type EQ 1 Imposed current vs time defined by a load curve EQ 2 Imposed voltage vs time defined by a load curve EQ 3 R L C V0 circuit EQ 4 Induced circuit Obsolete LS DYNA R7 0 4 5 EM EM VARIABLE LCID R F L A VO SIDCURR SIDVIN SIDVOUT PARTID Remarks EM CIRCUIT DESCRIPTION EQ 11 Imposed current defined by an amplitude A frequency F and initial time ty 1 Asin 2zF t 60 EQ 12 Imposed voltage defined by an amplitude A frequency F and initial time tg V Asin 27F t t9 Load curve ID for CIRCTYP 1 or 2 Value of the circuit resistance for CIRCTYP EQ 3 Value of the Frequency for CIRCTYP EQ 11 or 12 Value of the circuit inductance for CIRCTYP EQ 3 Value of the Amplitude for CIRCTYP EQ 11 or 12 Value of the circuit capacity for CIRCTYP EQ 3 Value of the initial time t0 for CIRCTYP EQ 11 or 12 Value of the circuit initial voltage for CIRCTYP EQ 3 Segment set ID for the current It uses the orientation given by the nor mal of the segments To use the opposite orientation use a minus sign in front of the segment set 14 CIRCTYP EQ 1 11 The current is imposed through this segment set CIRCTYP EQ 3 The current needed by the circuit equations is measured through this segment set Segment se
34. NDTIM 1 1 VARIABLE DESCRIPTION EMSOL Electromagnetism solver selector EQ 1 eddy current solver EQ 2 induced heating solver EQ 3 resistive heating solver NUMLS Number of local EM steps in gt period for EMSOL 2 Not used for EMSOL 1 DTINIT Initial electromagnetism time step DTMAX Maximum electromagnetism time step T INIT Time when electromagnetism solver is turned on default is at the begin ning of the run T END Time when electromagnetism solver is turned off default is at the end of the run NCYCLFEM Number of electromagnetism cycles between the recomputation of EM FEM matrices NCYCLBEM Number of electromagnetism cycles between the recomputation of EM BEM matrices 4 12 EM LS DYNA R7 0 CONTROL CONTACT EM EM CONTROL CONTACT Purpose Turns on the electromagnetism contact algorithms to check for contact between con ductors and allow the electromagnetic fields to flow from one conductor to another when detect ed as in contact Card 1 1 2 3 4 5 6 7 8 Variable EMCT Type I Default 0 VARIABLE DESCRIPTION EMCT EQ 0 no contact detection EQ 1 contact detection LS DYNA R7 0 4 13 EM CONTROL TIMESTEP EM CONTROL TIMESTEP Purpose Allows to control the EM time step and its evolution Card 1 1 2 3 4 5 6 7 8 Variable TSTYPE DTCONS LCID FACTOR Type I F I F Default none none none 1 0 VARIA
35. NT Purpose Initializes the chemistry and fluid state in every element of the list of CESE elements This is only used when chemistry is being solved with the CESE solver Card 1 1 2 3 4 5 6 7 8 Variable CHEMID COMPID Type Default none Card 2 1 2 3 4 5 6 7 8 Variable Type Default Card 3 Variable Type Default VARIABLE DESCRIPTION CHEMID Identifier of chemistry control card to use COMPID Identifier of chemical composition to use UIC X component of the fluid velocity LS DYNA R7 0 2 27 CESE CESE VARIABLE VIC WIC RHOIC PIC TIC HIC ELEI CESE INITIAL CHEMISTRY ELEMENT DESCRIPTION Y component of the fluid velocity Z component of the fluid velocity Initial fluid density Initial fluid pressure Initial fluid temperature Initial fluid enthalpy However when CHEMID refers to a ZND 1 step reaction card this is the progressive variable degree of combustion User element numbers to initialize 2 28 CESE LS DYNA R7 0 CESE INITIAL CHEMISTRY PART CESE CESE INITIAL CHEMISTRY PART Purpose Initializes the chemistry and fluid state in every element of the specified CESE part that has not already been initialized CESE INITIAL CHEMISTRY ELEMENT CESE INITIAL CHEMISTRY SET cards This is only used when chemistry is being solved with the CESE solver
36. OpenDX format A directory named output dx has to be created one level above the executable EQ 3 output a file with mesh statistics and the fluid results in OpenDX GMV and the mesh visualization tool medit format Three new directories have to be created output dx output gmv and out put medit each one level above the executable DTOUT Time interval to print the output when OUTL is different than 0 LS DYNA R7 0 5 17 ICFD ICFD ICFD CONTROL PARTITION ICFD CONTROL PARTITION Purpose This keyword changes the default option for the partition in MPP thus it is only valid in MPP Card 1 1 2 3 4 5 6 7 8 Variable PTECH Type I Default 1 VARIABLE DESCRIPTION PTECH Indicates the type of partition EQ 1 the library Metis is used EQ 2 partition along the axis with higher aspect ratio EQ 3 partition along X axis EQ 4 partition along Y axis EQ 5 partition along Z axis 5 18 ICFD LS DYNA R7 0 ICFD CONTROL SURFMESH CFD ICFD CONTROL SURFMESH Purpose This keyword enables automatic surface re meshing The objective of the re meshing is to improve the mesh quality on the boundaries It should not be used on a regular basis Card 1 1 2 3 4 5 6 7 8 Variable RSRF Type I Default 0 VARIABLE DESCRIPTION RSRF Indicates whether or not to perform a surface re meshing EQ 0 no re meshing is applied EQ 1 allows the sur
37. PTION X component of the initial particle velocity the TBX explosive Y component of the initial particle velocity the TBX explosive Z component of the initial particle velocity the TBX explosive 1 Ifradiation heat transfer is being modeled then EMISS and BOLTZ are required LS DYNA R7 0 7 7 STOCHASTIC LSO LSO LSO These cards provide a general data output mechanism causing the creation of a sequence of LSDA files This facility is intended to allow several different time sequences of data to be out put in the same simulation In addition any number of domains and any number of variables on those domains may be specified within each time sequence The keyword cards in this section are defined in alphabetical order LSO DOMAIN LSO POINT SET LSO TIME SEQUENCE An additional option TITLE may be appended to all LSO keywords If this option is used then an 80 character string is read as a title from the first card of that keyword s input At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks LS DYNA 7 0 8 1 LSO LSO LSO DOMAIN LSO DOMAIN Purpose This command provides a way to output data for specific variables at a frequency de fined by LSO TIME SEQUENCE and for a specific set of data structure Points nodes ele ments parts etc See LSO POINT SET in a specific binary LSO file Card 1 DOMAIN TYPE Card 2 1 Variable SOLVE
38. R NAME Type C When DOMAIN TYPE is THIST POINT the following cards apply Card 3 1 2 3 4 5 6 7 8 Variable OUTID SETID Type I I Default none none Card 4 VARIABLE NAME 8 2 LSO LS DYNA R7 0 LSO_DOMAIN VARIABLE DOMAIN TYPE SOLVER NAME OUTID SET ID VARIA BLE NAME Remarks LSO DESCRIPTION The type of domain for which LSO output may be generated Accepted entries so far are thist point Selects the solver from which data is output on this domain Accepted entries so far are em cese and icfd Output frequency ID associated to the domain See LSO TIME SEQUENCE Support set ID See LSO POINT SET Name of a single output variable See Remarks 1 When DOMAIN TYPE is THIST POINT the following variable names are accepted but note that some of these variable names may change in future releases more variables will al so be added Solver Name Variable type EM ICFD CESE Vectors currentDensity2 point velocity thist velocity thist ScalarPotential pt Pressure thist pressure thist Scalars temperature thist temperature thist density thist density thist thist Table 8 1 LS DYNA R7 0 8 3 LSO LSO LSO POINT SET LSO POINT SET Purpose Define a list of points used to sample variables in time Of the different sampling meth ods the most common one is to specify point
39. RCEBEM uode qutd PD UE d t RR t 4 33 SEMSOLVER EEMBENGs ui te Ao beth b m eta rod re eR 4 34 VICE eat a Rd x eA a LR AM M ERU ut un dard 5 1 ICFD_BOUNDARY_CONJ_ 5 3 WICH BOUNDARY EREESLIPA sd sedet 5 4 ICFD_BOUNDARY ES lia ia 5 5 ICRF BOUNDARY NONSUIP ute past te av uiros a 5 6 BOUNDARY PRESCRIBED MOVEMESH cnn 5 7 PIC PD BOUNDARY PRESCRIBED harten 5 8 ICFD BOUNDARY PRESCRIBED PRE sette 5 9 CFD BOUNDARY FLUX TEMP ese ttt ttn 5 10 BOUNDARY PRESCRIBED TEMP 5 11 FIGED CONTROL ADAPT dotes beta eo Send elo i calo 5 12 ICFD CONTROL ADAPT SIZE Luciae p c o p ER P Uu 5 13 CONTROL scheinen aedi ode uidet de Midi 5 14 SICED CONTBOE AMESEA qe dint sian teta dab bebe arena fida iens 5 15 ICFD_CONTROL_MESH_MOV leer tetas tars 5 16 ICFD CONTROL OUTPUT tite ear toti 5 17 CFD CONTROL 5 18 0 4 LS DYNA R7 0 TABLE OF CONTENTS MICED CONTROL SURE MESH ot e a leit ae des 5 19 CONTROL SPLIT ena tasca ese tal alto 5 20 ICFD_CONTROL_TIME sd e causa rie alien cus eh pu a hl Moe LEES 5 01 ICFD CONTROL TURBULENCE ttti rtl ebd nud 5 22 Lo eee t 5 24 PICEA TAB ASH AVERAGE eit ot set na lea ac tL 5 25 CFD DEFINE POINT
40. T SIZE CFD ICFD CONTROL ADAPT SIZE Purpose This keyword controls the re meshing of elements taking into account the element qual ity and distortion overwriting the default which only checks for inverted elements Card 1 1 2 3 4 5 6 7 8 Variable ASIZE Type I Default 0 VARIABLE DESCRIPTION ASIZE EQ 0 only re mesh in cases where elements invert EQ 1 re mesh if elements invert or if element quality deteriorates LS DYNA R7 0 5 13 ICFD ICFD ICFD CONTROL FSI ICFD CONTROL FSI Purpose This keyword modifies default values for the fluid structure interaction coupling Card 1 1 2 3 4 5 6 7 8 Variable OWC BT DT Type I F F Default 0 0 1 20 VARIABLE DESCRIPTION OWC Indicates the coupling direction to the solver EQ 0 two way coupling Loads and displacements are trans ferred across the FSI interface and the full non linear problem is solved EQ 1 one way coupling The solid solver transfers displace ments to the fluid solver EQ 2 one way coupling The fluid solver transfers stresses to the solid solver BT Birth time for the FSI coupling Before BT the fluid solver will not pass any loads to the structure but it will receive displacements from the solid solver DT Death time for the FSI coupling After DT the fluid solver will not trans fer any loads to the solid solver but it will continue to deform with the solid 5
41. VAP IPULSE LIMPRT IDFUEL STOCHASTIC TILTXZ CONE DCONE ANOZ AMPO DESCRIPTION Spray particle size distribution EQ 1 uniform EQ 2 Rosin Rammler default EQ 3 Chi squared degree of 2 EQ 4 Chi squared degree of 6 Type of particle breakup model EQ 0 off no breakup EQ 1 TAP EQ 2 KHRT Turn collision modeling on or off Turn evaporation on or off Type of injection EQ 0 continuous injection EQ 1 sine wave EQ 2 square wave Upper limit on the number of parent particles modeled in this spray This is not used with the continuous injection case IPULSE 0 Selected spray liquid fuels EQ 1 Default H20 EQ 2 Benzene C6H6 EQ 3 Diesel 2 C12H26 EQ 4 Diesel 2 C13H13 EQ 5 Ethanol 2 50 EQ 6 Gasoline C8H18 EQ 7 Jet A C12H23 EQ 8 Kerosene C12H23 EQ 9 Methanol CH3OH EQ 10 N dodecane C12H26 LS DYNA R7 0 7 3 STOCHASTIC STOCHASTIC STOCHASTIC_SPRAY PARTICLES VARIABLE RHOP TIP PMASS PRTRTE STRINJ DURINJ XORIG YORIG ZORIG SMR VELINJ DRNOZ DTHNOZ TILTXY TILTXZ CONE DCONE ANOZ AMPO DESCRIPTION Particle density Initial particle temperature Total particle mass Number of particles injected per second for continuous injection Start of injection s Duration of injection s X coordinate of center of a nozzle s exit plane Y coordinate of center of a nozzle s exit plane Z coordinate of c
42. VEL YVEL ZVEL Type F F F F F F Default none none none 0 0 0 0 0 0 LS DYNA R7 0 7 5 STOCHASTIC STOCHASTIC STOCHASTIC_TBX_PARTICLES VARIABLE PCOMB NPRTCL MXCNT PMASS SMR RHOP TICP T_IGNIT INITDST AZIMTH ALTITD CPS CVS HVAP EMISS BOLTZ XORIG YORIG ZORIG DESCRIPTION Particle combustion mode EQ 0 no burning EQ 1 K model EQ 2 Hybrid Initial total number of parent particles discrete particles for calculation Maximum allowed number of parent particles during the simulation Total particle mass Sort mean particle radius Particle density Initial particle temperature Particle ignition temperature Initial particle distribution EQ 1 spatially uniform EQ 2 Rosin Rammler EQ 3 Chi squared Angle in degrees from X axis in X Y plane of reference frame of TBX explosive 0 lt AZMITH lt 360 Angle in degrees from Z axis of reference frame of TBX explosive 0 lt ALTITD lt 180 Heat coefficient Latent heat of vaporization Particle emissivity Boltzmann coefficient X coordinate of the origin of the initial reference frame of the TBX ex plosive Y coordinate of the origin of the initial reference frame of the TBX ex plosive Z coordinate of the origin of the initial reference frame of the TBX ex plosive 7 6 STOCHASTIC LS DYNA R7 0 STOCHASTIC_TBX_PARTICLES VARIABLE XVEL YVEL ZVEL Remarks STOCHASTIC DESCRI
43. Variable PIDI PID2 PID3 PID4 5 PID6 PID7 PID8 Type I I I I I I I I Default none none none none none none none None VARIABLE DESCRIPTION VOLID ID assigned to the new volume in the keyword MESH_VOLUME The interface meshes will be applied to this volume PIDn Part IDs for the surface elements LS DYNA R7 0 6 5 MESH MESH MESH SIZE MESH SIZE Purpose Define the surfaces that will be used by the mesher to specify a local mesh size inside the volume If no internal mesh is used to specify the size the mesher will use a linear interpola tion of the surface sizes that define the volume enclosure Card 1 1 2 3 4 5 6 7 8 Variable VOLID Type I Default none Cards 2 3 4 the next card terminates the input 1 2 3 4 5 6 7 8 Variable PIDI PID2 PID3 PID4 PIDS PID6 PID7 PID8 Type I I I I I I I I Default none none none none none none none None VARIABLE DESCRIPTION VOLID ID assigned to the new volume in the keyword MESH VOLUME The mesh sizing will be applied to this volume PIDn Part IDs for the surface elements that are used to define the mesh size next to the surface mesh 6 6 MESH LS DYNA R7 0 MESH_SIZE_SHAPE MESH MESH_SIZE_SHAPE Purpose Defines a local mesh size in specific zones corresponding to given geometrical shapes box sphere cylinder The solver will automatically apply the conditions specified during the
44. a 8 4 TIME SEOUENCE2 2 2 e eae p lat nn tela eee d chus 8 5 LS DYNA R7 0 0 5 INTRODUCTION INTRODUCTION LS DYNA MULTIPHYSICS USER S MANUAL INTRODUCTION In this manual there are three main solvers a compressible flow solver an incompressible flow solver and an electromagnetism solver Each of them implements coupling with the structural solver in LS DYNA The keywords covered in this manual fit into one of three categories In the first category are the keyword cards that provide input to each of the multiphysics solvers that in turn couple with the structural solver In the second category are keyword cards involving extensions to the basic solvers Presently the chemistry and stochastic particle solvers are the two solvers in this catego ry and they are used in conjunction with the compressible flow solver discussed below In the third category are keyword cards for support facilities A volume mesher that creates volume tet rahedral element meshes from bounding surface meshes is one of these tools Another is a new data output mechanism for a limited set of variables from the solvers in this manual This mech anism is accessed through LSO keyword cards The CESE solver is a compressible flow solver based upon the Conservation Element Solution Element CE SE method originally proposed by Chang of the NASA Glenn Research Center This method is a novel numerical framework for conservation laws It has m
45. any non traditional features including a unified treatment of space and time the introduction of separate conserva tion elements CE and solution elements SE and a novel shock capturing strategy without us ing a Riemann solver This method has been used to solve many types of flow problems such as detonation waves shock acoustic wave interaction cavitating flows supersonic liquid jets and chemically reacting flows In LS DYNA it has been extended to also solve fluid structure inter action FSI problems It does this with two approaches The first approach solves the compress ible flow equations on an Eulerian mesh while the structural mechanics is solved on a moving mesh that moves through the fixed CE SE mesh In the second approach new with this version the CE SE mesh moves in a fashion such that its FSI boundary surface matches the correspond ing FSI boundary surface of the moving structural mechanics mesh This second approach is more accurate for FSI problems especially with boundary layers flows Another new feature with the CESE moving mesh solver is conjugate heat transfer coupling with the solid thermal solver The chemistry and stochastic particle solvers are two addon solvers that further extend the CESE solver The second solver 15 incompressible flow solver that is fully coupled with the solid mechanics solver This coupling permits robust FSI analysis via either an explicit technique when the FSI is weak or
46. constant value of the densi ty for this boundary condition will be SF RHO 2 10 CESE LS DYNA R7 0 CESE_BOUNDARY_REFLECTIVE CESE CESE BOUNDARY REFLECTIVE OPTION Available options are PART SET SEGMENT Purpose Define a reflective boundary condition RBC for the CESE compressible flow solver This boundary condition can be applied on a symmetrical surface or a solid wall for inviscid flow only of the computational domain For the PART option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SURFPRT Type I Default none For the SET option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SSID Type I Default none LS DYNA R7 0 2 11 CESE CESE CESE BOUNDARY REFLECTIVE For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 7 8 Variable NI N2 N3 N4 Type I I I I Default none none none none VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards SSID Segment set ID N1 N2 Node ID s defining a segment Remarks 1 This boundary condition has the same effect as a solid wall boundary condition for inviscid flows 2 12 CESE LS DYNA R7 0 CESE BOUNDARY SOLID WALL CESE CESE BOUNDARY SOLID WALL OPTION OPTION2 For OPTIONI the choices are PART SET SEGMENT For OPTION the choices are lt BLANK
47. enter of a nozzle s exit plane Sauter mean radius Injection velocity Nozzle radius Azimuthal angle in degrees measured counterclockwise of the injector nozzle from the j 1 plane Rotation angle in degrees of the injector in the x y plane where 0 0 points towards the 3 o clock position j 1 line and the angle increases counterclockwise from there Inclination angle in degrees of the injection in the x z plane where 0 0 points straight down x gt 0 0 points in the positive x direction and x lt 0 0 points in the negative x direction Spray mean cone angle in degrees for hollow cone spray spray cone angle in degrees for solid cone spray Injection liquid jet thickness in degrees Area of injector Initial amplitude of droplet oscillation at injector 7 4 STOCHASTIC LS DYNA R7 0 STOCHASTIC_TBX_PARTICLES STOCHASTIC STOCHASTIC_TBX_PARTICLES Purpose Specify particle and other model details for stochastic PDEs that model embedded par ticles in TBX explosives Card 1 1 2 3 4 5 6 7 8 Variable PCOMB NPRTCL MXCNT PMASS SMR RHOP TICP T_IGNIT Type I I I F F F F F Default 0 none none none none none none none Card 2 1 2 3 4 5 6 7 8 Variable INITDST AZIMTH ALTITD CPS CVS HVAP EMISS BOLTZ Type I F F F F F F Default 1 none none none none none none Remarks 1 1 Card 3 1 2 3 4 5 6 7 8 Variable XORIG YORIG ZORIG X
48. etric CESE fluid solver LS DYNA R7 0 2 3 CESE CESE CESE BOUNDARY FSI OPTION Available options are PART SET SEGMENT CESE BOUNDARY FSI Purpose Define an FSI boundary condition for the CESE compressible flow solver This boundary condition must be applied on a surface of the CESE computational domain that is shared with surfaces of the outside boundary of the structural mechanics solver The nodes of the two meshes will generally not be shared For the PART option define the following card Card 1 1 2 3 8 Variable SURFPRT Type I Default none For the SET option define the following card Card 1 1 2 3 8 Variable SSID Type I Default none 2 4 CESE LS DYNA R7 0 CESE_BOUNDARY FSI CESE For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 7 8 Variable N2 N3 N4 Type I I I I Default none none none none VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards SSID Segment set ID NI N2 Node ID s defining a segment Remarks 1 This boundary condition card is also needed for conjugate heat transfer problems with the CESE solver LS DYNA R7 0 2 5 CESE CESE CESE BOUNDARY NON REFLECTIVE CESE BOUNDARY NON REFLECTIVE OPTION Available options are PART SET SEGMENT Purpose Define a passive boundary conditi
49. face mesh to be re meshed LS DYNA R7 0 5 19 ICFD ICFD ICFD CONTROL SPLIT ICFD CONTROL SPLIT Purpose This keyword provides an option to trigger an iterative procedure on the fluid system This procedure aims to bring more precision to the final pressure and velocity values but is often very time consuming It must therefore be used with caution in specific cases For stability pur poses this method is automatically used for the first ICFD time step Card 1 1 2 3 4 5 6 7 8 Variable NIT TOL Type I F Default none none VARIABLE DESCRIPTION NIT Maximum Number of iterations of the system for each fluid time step If TOL criteria is not reached after NIT iterations the run will proceed TOL Tolerance Criteria for the pressure residual during the fluid system solve 5 20 ICFD LS DYNA R7 0 ICFD CONTROL TIME CFD CONTROL TIME Purpose This keyword is used to change the defaults related to time parameters in the fluid prob lem Card 1 1 2 3 4 5 6 7 8 Variable TTM DT CFL Type F F F Default 1 28 0 1 VARIABLE DESCRIPTION TTM Total time of simulation for the fluid problem DT Time step for the fluid problem If different from zero the time step will be set constant and equal to this value If DT 0 then the time step is automatically computed based on the CFL condition CFL Scale factor that multiplies DT LS DYNA R7
50. gment set are eliminated from the BEM calculations used for example for the rear or side faces of a workpiece EQ 10 EQ 11 used by the 2D axisymmetric solver to make the connection between two corresponding boundaries on each side of a slice when the model is a slice of the full 360 circle See EM ROTATION AXIS card A boundary of type 10 and LTYPE n is connected to the boundary of type 11 and LTYPE n LTYPE Loadtype used to identify two boundaries of BTYPE 10 and 11 The same value of loadtype identifies two boundaries connected together 4 4 EM LS DYNA R7 0 EM_CIRCUIT EM EM_CIRCUIT_ OPTION Available options include SOURCE Purpose Define an electrical circuit For the SOURCE option the current will be considered uniform in the circuit This can be use ful in order to save computational time in cases with a low frequency current and where the dif fusion of the EM fields is a very fast process This option is in contrast with the general case where the current density in a circuit is completed in accordance with the solver type defined in EMSOL of EM_ CONTROL For example if an eddy current solver is selected the diffusion of the current in the circuit is taken into account Card 1 1 2 3 4 5 6 7 8 Variable CIRCID CIRCTYP LCID R F L A C to 1 1 1 Default none none none none none none none Card 2 1 2 3 4 5 6 7 8 Variable SIDCURR SIDVI
51. gt ROTAT Purpose Define a solid wall boundary condition SBC for this CESE compressible flow solver This boundary condition can be applied at a solid boundary that is the physical boundary for the flow field For inviscid flow this will be a slip boundary condition while for viscous flows it is no slip boundary condition For the PART option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SURFPRT LCID Vx Vy Vz Nx Ny Nz Type I I F F F F F F Default none 0 0 0 0 0 0 0 0 0 0 0 0 0 Remarks 2 3 2 2 2 3 3 3 LS DYNA R7 0 2 13 CESE CESE For the SET option define the following card CESE BOUNDARY SOLID WALL Card 1 1 2 3 4 5 6 7 8 Variable SSID LCID Vx Vy Vz Nx Ny Nz Type I I F F F F F F Default none 0 0 0 0 0 0 0 0 0 0 0 0 0 Remarks 2 3 2 2 2 3 3 3 For the SEGMENT option define the following two cards Card 1 1 2 3 4 5 6 7 8 Variable N2 N3 N4 LCID Vx Vy Vz Type I I I I I F F F Default none none none none 0 0 0 0 0 0 0 Remarks 2 3 2 2 2 Card 2 1 2 3 4 5 6 7 8 Variable Ny Nz Type F F Default 0 0 0 0 Remarks 3 3 VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards 2 14 CESE LS DYNA R7 0 CESE BOUNDARY SOLID WALL CESE VARIABLE SSID N1 N2 LCID If OPTION2 Vx Vu Vz If OPTION2 Vx Vy Vz Nx Ny
52. hen an 80 character string is read as a title from the first card of that keyword s in put At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks LS DYNA R7 0 3 1 CHEMISTRY CHEMISTRY CHEMISTRY_COMPOSITION CHEMISTRY_COMPOSITION Purpose Provides a general way to specify a chemical composition via a list of species mole numbers in the context of a chemkin database model Card 1 1 2 3 4 5 6 7 8 Variable ID Type I Default none Card 2 1 2 Variable MOLFR SPECIES Type F A Default none none VARIABLE DESCRIPTION ID A unique identifier among all chemistry compositions MOLFR The mole number corresponding to the species named in the SPECIES field SPECIES The chemkin compatible name of a chemical species 3 2 CHEMISTRY LS DYNA R7 0 CHEMISTRY CONTROL 0D CEMISTRY CONTROL 0D CHEMISTRY Purpose Performs an isotropic chemistry calculation that operates standalone does not call the CESE solver This is for ISOBARIC or ISOCHORIC cases Card 1 1 2 3 4 5 6 7 8 Variable ID COMPID SOLTYP ERRLIM Type I I I F Default none none none 1 0 3 Card 2 1 2 3 4 5 6 7 8 Variable DT TLIMIT TIC PIC RIC EIC Type F F F F F F Default none none none none none none Card 3 Variable FILEI T
53. iangular in 3 d and linear segments in 2 d that will be used by the mesher to construct a volume mesh These surface elements may be used to define the enclosed volume to be meshed or alternatively they could be used to apply different mesh sizes inside the volume see card MESH SIZE Card 1 1 2 3 4 5 6 7 8 T VARIABLE DESCRIPTION EID Element ID A unique number with respect to all MESH SURFACE ELEMENTS cards PID Part ID A unique part identification number NI Nodal point 1 N2 Nodal point 2 N3 Nodal point 3 N4 Nodal point 4 Remarks 1 The convention is the same used by the keyword ELEMENT_SHELL In the case of a tri angular face N3 N4 In 2 d N2 N3 N4 LS DYNA R7 0 6 9 MESH MESH MESH_ SURFACE NODE MESH SURFACE NODE Purpose Define a node and its coordinates These nodes will be used in the mesh generation process by the MESH SURFACE ELEMENT keyword Card 1 1 2 3 4 5 6 7 8 Variable NID x Y Z Type I F F F Default none 0 0 0 VARIABLE DESCRIPTION NID Node ID A unique number with respect to the other surface nodes X coordinate Y y coordinate Z z coordinate 6 10 MESH LS DYNA R7 0 MESH_VOLUME MESH MESH_ VOLUME Purpose This keyword defines the volume space that will be meshed The boundaries of the vol ume are the surfaces defined by MESH SURFACE ELEMENT The list of surfaces have to be non overlapping and should not
54. icles in TBX explosives and a spray model The cards for using these models are STOCHASTIC SPRAY PARTICLES STOCHASTIC TBX PARTICLES An additional option TITLE may be appended to all STOCHASTIC keywords If this op tion is used then an 80 character string is read as a title from the first card of that keyword s in put At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks LS DYNA R7 0 7 1 STOCHASTIC STOCHASTIC STOCHASTIC SPRAY PARTICLES STOCHASTIC SPRAY PARTICLES Purpose Specify particle and other model details for spray modeling using stochastic PDEs that approximate such processes A pair of cards is required to specify the characteristics of each nozzle cards 3 and 4 describe the first nozzle Card 1 1 2 3 4 5 6 7 Variable INJDIST IBRKUP ICOLLDE IEVAP IPULSE LIMPR IDFUEL Type I I I I I I I Default 1 none none none none none 1 Card 2 1 2 3 4 5 6 7 Variable RHOP TIP PMASS PRTRTE STRINJ DURINJ Type F F F F F F Default none none none none none none Card 3 5 1 2 3 4 5 6 7 Variable XORIG YORIG ZORIG SMR VELINJ DRNOZ DTHNOZ Type F F F F F F F 7 2 STOCHASTIC LS DYNA R7 0 STOCHASTIC SPRAY PARTICLES Card 4 6 Card N Variable Type Default TILTXY VARIABLE INJDIST IBRKUP ICOLLDE IE
55. inear system EQ 1 starts from previous solution normalized by the rhs change 4 30 EM LS DYNA R7 0 SOLVER BEM Remarks 1 Using USELAST 1 can save many iterations in the further solves if the vector solution of the present solve is assumed to be nearly parallel to the vector solution of the previous solve as usually happens in time domain eddy current problems LS DYNA R7 0 4 31 EM EM SOLVER BEMMAT EM SOLVER BEMMAT Purpose Define the type of BEM matrices as well as the way they are assembled Card 1 1 2 3 4 5 6 7 8 ee i VARIABLE DESCRIPTION MATID Defines which BEM matrix the card refers to EQ 1 P matrix EQ 2 Q matrix RELTOL Relative tolerance on the sub blocks of the matrix when doing low rank approximations The default values are 1 e 3 for the P matrix and 1 2 for the Q matrix The user should try to decrease these tolerances if the results are not accurate enough More memory will then be needed 4 32 EM LS DYNA R7 0 EM_SOLVER_FEM EM EM SOLVER FEM Purpose Define some parameters for the EM FEM solver Card 1 1 2 3 4 5 6 7 8 RELTOL MAXITE STYPE i VARIABLE DESCRIPTION RELTOL Relative tolerance for the iterative solvers PCG or GMRES The user should try to decrease this tolerance if the results are not accurate enough More iterations will then be needed MAXITER Maximal number of iterations STYPE Solver type EQ 1 Direct so
56. ion of the CESE_INITIAL_CHEMISTRY cards Card 1 1 2 3 4 5 6 7 8 Variable ID COMPID NMESH DLEN CFL TLIMIT XYZD DETDIR Type I I I F F F F I Default none none none none none none none none Card 2 1 Variable FILEI Type A Card 3 1 Variable FILE2 Type A Card 4 1 Variable FILE3 Type A 3 10 CHEMISTRY LS DYNA R7 0 CHEMISTRY_DET_INITIATION CHEMISTRY Card 5 1 Variable FILE4 Type A VARIABLE DESCRIPTION ID Identifier for this one dimensional detonation computation COMPID Chemical composition identifier of composition to use NMESH Number of equal width elements in the one dimensional domain DLEN Length of the one dimensional domain CFL Time step limiting factor TLIMIT Time limit for the simulation XYZD Position of the detonation front in the DETDIR direction DETDIR Detonation propagation direction 1 gt X 2 gt Y 3 gt Z FILEI Name of the LSDA file in which to write the one dimensional solution FILE2 Name of the file containing the Chemkin compatible input FILE3 Name of the file containing the chemistry thermodynamics database FILE4 Name of the file containing the chemistry transport properties database LS DYNA R7 0 3 11 CHEMISTRY CHEMISTRY CHEMISTRY_PATH CHEMISTRY PATH Purpose To specify one or more search paths to look for chemistry database files Card 1 1 VARIABLE DESCR
57. istry transport properties database 3 8 CHEMISTRY LS DYNA R7 0 CHEMISTRY_CONTROL_ZND CHEMISTRY CHEMISTRY_CONTROL_ZND Purpose Computes the one dimensional reduced chemistry of a ZND model It is then used in the initialization of the chemistry part of the CESE solver When this card is used the CESE_INITIAL_CHEMISTRY cards must specify the progressive variable degree of com bustion in the HIC field Card 1 1 2 3 4 5 6 7 8 Variable ID Type I Default none Card 2 1 2 3 4 5 6 7 8 Variable F EPLUS 00 XYZD DETDIR Type F F F F F I Default none none none none none none VARIABLE DESCRIPTION ID Identifier for this full chemistry calculation F Overdriven factor EPLUS EPLUS parameter of the ZND model Q0 QO parameter of the ZND model GAM GAM parameter of the ZND model XYZD Position of the detonation front in the DETDIR direction DETDIR Detonation propagation direction 1 gt X 2 gt Y 3 gt Z LS DYNA R7 0 3 9 CHEMISTRY CHEMISTRY CHEMISTRY_DET_INITIATION CHEMISTRY_DET_INITIATION Purpose Performs a one dimensional detonation calculation based upon a chemical composition and initial conditions It is then available for use immediately in the CESE solver for initializing a computation or it can be subsequently used by CHEMISTRY_CONTROL_1D card in a later run In the product regions this card overrides the initializat
58. l directions respectively 4 Currently the FSI solver only runs in 3D cases LS DYNA R7 0 2 19 CESE CESE CESE_CONTROL_TIMESTEP CESE CONTROL TIMESTEP Purpose Sets the time step control parameters for the CESE compressible flow solver Card 1 1 2 3 4 5 6 7 8 IDDT CFL DTINT Default 1 0E 3 VARIABLE DESCRIPTION IDDT Sets the time step option EQ 0 Fixed time step size DTINT 1 given initial time step size Otherwise the time step size will be calculated based on the given CFL number and the flow solution at the previous time step CFL CFL number Courant Friedrichs Lewy condition 0 0 lt CFL lt 1 0 DTINT Initial time step size 2 20 CESE LS DYNA R7 0 CESE EOS EQUILIB CESE CESE EOS HOMOG EQUILIB Purpose Define the coefficients in the equation of state EOS for the homogeneous equilibrium cavitation model Card 1 1 2 3 4 5 6 7 8 Variable EOSID p liq a vap a liq u liq P sat vap Type I F F F F F F F Default none 0 8 880 0 334 0 1386 0 1 435e 5 1 586e 4 1 2e 4 VARIABLE DESCRIPTION EOSID Equation of state identifier p_vap density of the saturated vapor density of saturated liquid sound speed of saturated vapor a liq sound speed of the saturated liquid dynamic viscosity of vapor u liq dynamic viscosity of the liquid P_sat_vap pressure of the
59. leave any gaps or open spaces between the surface bounda ries The nodes on the boundary of two neighbor surfaces have to be uniquely defined by the keyword MESH SURFACE NODE and should match exactly on the interface This card will be ignored if the volume mesh is specified by the user and not generated automatically Card 1 1 2 3 4 5 6 7 8 Variable VOLID 4 Default Cards 2 3 4 the next card terminates the input PID1 PID2 PID3 PID4 PID6 PID7 PID8 Default none none None VARIABLE DESCRIPTION VOLID ID assigned to the new volume PIDn Part IDs for the surface elements that are used to define the volume LS DYNA R7 0 6 11 MESH MESH MESH VOLUME ELEMENT MESH VOLUME ELEMENT Purpose Specify a set of volume elements for the fluid volume mesh in cases where the volume mesh is specified by the user and not generated automatically The Nodal point are specified in the MESH VOLUME NODE keyword Only Tets are supported triangles in 2D Card 1 1 2 3 4 5 6 7 8 Variable EID PID Nl N2 N3 N4 Type I I I I I I Default none none none none none None VARIABLE DESCRIPTION EID Element ID A unique number with respect to all MESH_VOLUME_ELEMENTS cards PID Part ID A unique part identification number NI Nodal point 1 N2 Nodal point 2 N3 Nodal point 3 N4 Nodal point 4 Remarks 1 The convention is the same used by the keyword ELEMENT_SHE
60. luid volumes LS DYNA R7 0 5 31 ICFD ICFD ICFD PART VOL PID 3 CN H Jp 5 m PID 3 Figure 5 1 VARIABLE DESCRIPTION SECID Section identifier defined by the SECTION card MID Material identifier SPIDn Part IDs for the surface elements that define the volume mesh 9999999999999999999999999999999999999999999999999999599999999999999999999999999 5585 NICFD PART VOL 5359995999399999955999999999999999995999999999999999999999999999999999959999599 5 5 PART ID 5 is defined by the surfaces that enclose it 5 ICFD PART VOL 5 5 secid mid 5 1 1 Seid fled i eee ei ee Mel Bul M MM PE 5 pidl pid2 pia3 194 195 196 7 pia8 1 2 3 4 5 32 ICFD LS DYNA R7 0 ICFD SECTION ICFD SECTION Purpose Define a section for the incompressible flow solver ICFD Card 1 1 2 3 4 5 8 Variable SID Type I Default none VARIABLE DESCRIPTION SID Section identifier LS DYNA R7 0 5 33 ICFD ICFD ICFD SET NODE ICFD SET NODE Purpose Only used in cases where the mesh is specified by the user See MESH VOLUME ELEMENT Defines a set of nodes associated to a part ID on which boundary conditions can be applied Card 1 1 2 3 4 5 6 gi 8 Card 2 1 3 4 8 NID3 NID4 NID8 2 T VARIABLE DESCRIPTION SID Set ID PID Associated Part ID NIDI Node ID n Remarks 1 The convention is the similar to the one used by the ke
61. lve EQ 2 Conditioned Gradient Method PCG PRECON Preconditioner type for PCG EQ 0 preconditioner EQ 1 Diagonal line USELAST This is used only for iterative solvers PCG EQ 1 starts from 0 as initial solution of the linear system starts from previous solution normalized by rhs change Remarks 1 Using USELAST 1 can save many iterations in the further solves if the vector solution of the present solve is assumed to be nearly parallel to the vector solution of the previous solve as usually happens in time domain eddy current problems LS DYNA R7 0 4 33 EM SOLVER FEMBEM EM SOLVER FEMBEM Purpose Define some parameters for the coupling between the EM FEM and EM BEM solv ers Card 1 1 2 3 4 5 6 7 8 Variable RELTOL MAXITE FORCON Type I Default 0 VARIABLE DESCRIPTION RELTOL Relative tolerance for the solver The user should try to decrease this tol erance if the results are not accurate enough More iterations will then be needed MAXITER Maximal number of iterations FORCON 0 the code stops with an error if no convergence the code continues to the next time step even if the RELTOL convergence criteria has not been reached 4 34 EM LS DYNA R7 0 ICFD ICFD ICFD The keyword ICFD covers all the different options available in the incompressible fluid solver The keyword cards in this section are defined in alphabetical order ICFD_ BOUNDARY
62. m temperature Omat 9 OBurgess room Remarks The Meadon model is a simplified Burgess model with the solid phase equations only The electrical resistivity is given by V ns C20 fe 1 where is the temperature V is the specific volume and Vp is the reference specific volume ze ro pressure solid phase In 1 the volume dependence is given by fe ca if EXPON 1 most materials 2a f C ore if EXPON 1 tungsten 2b EG 2y if EXPON 0 stainless steel 2c fe 1 if Vo is not defined or equal to 0 2d In this last case only EOSID C2 C3 TEMUNIT TIMUNIT and LGTUNIT need to be defined with Y Yo 20 7 2 3 4 20 EM LS DYNA R7 0 EM EOS MEADON The following table reports some sets of parameters given by Burgess in his paper EM Parameter Cu Ag Au W Al 2024 SS 304 Vo cm gm 0 112 0 0953 0 0518 0 0518 0 370 0 1265 Yo 2 00 2 55 3 29 1 55 2 13 2 00 C BUS 4 12e 5 3 37e 5 4 95e 5 9 73e 5 5 35e 5 0 C 0 113 0 131 0 170 0 465 0 233 0 330 C 1 145 1 191 1 178 1 226 1 210 0 4133 EXPON 1 1 1 1 1 0 Table 4 2 4 21 EM LS DYNA R7 0 EM EOS PERMEABILITY EM EOS PERMEABILITY Purpose Define the parameters for the behavior of a materials permeability Card 1 1 2 3 4 5 6 7 8 EOSID EOSTYPE LCID VARIABLE DESCRIPTION EOSID ID of the EOS EOSTYPE
63. material these materials have the same electromag netism behavior as 0 EQ 2 Conductor carrying a source In these conductors the eddy cur rent problem is solved which gives the actual current density Typi cally this would correspond to the coil EQ 4 Conductor not connected to any current or voltage source where the Eddy current problem is solved Typically this would cor respond to the workpiece SIGMA Initial electrical conductivity of the material EOSID ID of the EOS to be used for the electrical conductivity see EM EOS card MUREL Relative permeability Is the ratio of the permeability of a specific me dium to the permeability of free space EOSMU ID of the EOS to be used to define the behavior of u by an equation of state Note 1f EOSMU 15 defined MUREL will be used for the initial value only 4 26 EM LS DYNA R7 0 EM OUTPUT EM EM OUTPUT Purpose Define the level of EM related output on the screen and in the messag file Card 1 1 2 3 4 5 6 7 8 Variable MATS MATF SOLS SOLF MEM TIMING Default 0 0 i i VARIABLE DESCRIPTION MATS Level of matrix assembly output to the screen 0 No output EQ 1 Basic assembly steps EQ 2 Basic assembly steps percentage completed final statistics EQ 3 Basic assembly steps percentage completed statistics at each percentage of completion MATF Level of matrix assembly output to the messag file 0 No output EQ
64. ne None VARIABLE DESCRIPTION PID Part ID for the volume elements where the values are initialized see ICFD PART VOL Vx x coordinate for the velocity Vy y coordinate for the velocity Vz z coordinate for the velocity T Initial temperature LS DYNA R7 0 5 27 ICFD ICFD ICFD MAT ICFD MAT Purpose Specify physical properties for the fluid material Card 1 1 2 3 4 5 6 7 8 Variable MID FLG RO VIS ST THD Type I I F F F F Default None 1 none 0 0 0 Card 2 1 2 3 4 5 6 7 8 Variable HC TC Beta Type F F F Default 0 0 0 Card 3 1 2 3 4 5 6 7 8 Variable NNID k n nmin nmax Type F F F F F Default 0 0 0 0 1E30 VARIABLE DESCRIPTION MID Material ID FLG Flag to choose between fully incompressible slightly compressible or barotropic flows So far only fully incompressible flow is supported RO Flow density 5 28 ICFD LS DYNA R7 0 ICFD VARIABLE DESCRIPTION VIS Dynamic viscosity ST Surface tension coefficient THD Thermal diffusion used in the solution of the thermal problem Does not need to be defined if HC and TC exist HC Heat capacity TC Thermal conductivity Beta Thermal expansion coefficient used in the Boussinesq approximation for buoyancy NNID Non Newtonian flows model ID EQ 1 Power Law k Power Law input parameter Measure of the average velocity of the fluid con
65. on a given cross section of the part defined by a segment set with a symmetry axis defined by its direction at this time it can be the x y or z axis The EM forces and Joule heating are then computed over the full 3D part by rotations The part needs to be compatible with the symmetry i e each node in the part needs to be the child of a parent node on the segment set by a rotation around the axis Only the conductor parts with a MAT of type 2 or 4 should be defined as 2D axisymmetric Card 1 1 2 3 4 5 6 7 8 Variable PID SSID DIR 2DOPT Type I I I I Default none none none none VARIABLE DESCRIPTION PID Part ID of the part to be solved using 2D axisymmetry SSID Segment Set ID Segment that will define the 2D cross section of the part where the EM field is solved DIR Defines the symmetry axis At this time only the x y or z axis can be used EQ 1 X axis EQ 2 Y axis EQ 3 Z axis Remarks At this time either all or none of the conductor parts should be 2D axisymmetric In the future a mix between 2D axisymmetric and 3D parts will be allowed LS DYNA R7 0 4 3 EM BOUNDARY BOUNDARY Purpose Define some boundary conditions for electromagnetism problems Card 1 1 2 3 4 5 6 7 8 Variable SSID BTYPE LTYPE Type I I I Default none none none e VARIABLE DESCRIPTION SSID Segment Set Id BTYPE EQ 9 The faces of this se
66. on for CESE compressible flows This non reflective boundary condition NBC provides an artificial computational boundary for an open boundary that is passive For the PART option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SURFPRT Type I Default none For the SET option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SSID Type I Default none 2 6 CESE LS DYNA R7 0 CESE_BOUNDARY_NON_REFLECTIVE CESE For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 7 8 Variable NI N2 N3 N4 Type I I I I Default none none none none VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards SSID Segment set ID N2 Node ID s defining a segment Remarks 1 This boundary condition is usually imposed on an open surface that is far from the main dis turbed flow the farther away the better 1 the flow on that boundary surface should be almost uniform 2 If any boundary segment has not been assigned a boundary condition by any of the CESE BOUNDARY cards then it will automatically be assigned this non reflective boundary condition LS DYNA 7 0 2 7 CESE CESE CESE_BOUNDARY_PRESCRIBED CESE BOUNDARY PRESCRIBED OPTION Available options include PART SET SEGMENT Purpose For the CESE compressible flow solver
67. pose Define Rogowsky coils to measure a global current vs time through a segment set or a node set Card 1 1 2 3 4 5 6 7 8 Variable ROGOID SETID SETTYPE CURTYP Type I I I I Default 0 0 0 0 VARIABLE DESCRIPTION ROGOID Rogowsky coil ID SETID Segment or node set ID SETTYPE Type of set EQ 1 Segment set EQ 2 Node set not available yet CURTYP Type of current measured EQ 1 Volume current EQ 2 Surface current not available yet Remarks 1 A node set allows computing only a surface current An ASCII file xxx with xxx representing the rogold is generated for each EM_ROGO card giving the value of the cur rent vs time 4 8 EM LS DYNA R7 0 CONTACT CONTACT Purpose EM Define an electromagnetic contact between two sets of parts Used with the EM CONTACT RESISTANCE card Card 1 1 2 3 4 8 Variable CONTID COTYPE PSIDM PSIDS Type I I I I Default none none none none VARIABLE DESCRIPTION CONTID Electromagnetic contact ID COTYPE Type of EM contact EQ 1 Face to face PSIDM Master part set ID PSIDS Slave part set ID LS DYNA R7 0 4 9 EM EM CONTACT RESISTANCE EM CONTACT RESISTANCE Purpose Calculate the contact resistance of a previously defined EM contact in EM CONTACT Card 1 1 2 3 4 5 6 7 8 Variable CRID CONTID RHOprobe RHOsub R
68. quid phase and vapor phase At this time only the solid and liquid phases are implemented To check which elements are in the solid and in the liquid phase a melting temperature is first computed by 1 V 3 ee 0 amo 7 Dao Vo f T lt solid phase model The solid phase electrical resistivity corresponds to the Meadon model V ns fe 1 where is the temperature V is the specific volume and Vo is the reference specific volume ze ro pressure solid phase In 1 the volume dependence is given by X 1 if EXPON 1 most materials 2a f m if EXPON 1 tungsten 2b LS DYNA R7 0 4 17 EM EM EM EOS BURGESS f GO Gr if EXPON 0 stainless steel 2c with 1 V Y o gt Ug ob 3 If T liquid phase model n G 4 with Mem An ns e where An ke999LrF 8mif k gt 0 An 1 0 0772 2 Om if k 1 tungsten Sb An 1 0 106 0 846 Om if k 2 stainless steel SS 304 5c The following table reports some sets of parameters given by Burgess in his paper Parameter Cu Ag Au W Al 2024 SS 304 0 112 0 0953 0 0518 0 0518 0 370 0 1265 Yo 2 00 2 55 3 29 1 55 2 13 2 00 BUS 0 117 0 106 0 115 0 315 0 0804 0 156 Lr BUS 0 130 0 113 0 127 0 337 0 107 0 153 C BUS 4 12e 5 3 37 5 4 95e 5 9 73 5 5 35e 5 0 C 0 113 0 131
69. r the iterative linear solver LS DYNA R7 0 2 17 CESE CESE CESE CONTROL SOLVER CESE CONTROL SOLVER Purpose Set general purpose control variables for the CESE compressible flow solver Card 1 1 2 3 4 5 6 7 8 Variable IGEOM Type I Default 0 Remarks 1 2 3 4 VARIABLE DESCRIPTION IFRAME Sets the framework of the CESE solver EQ 0 Fixed Eulerian EQ 100 Moving Mesh FSI EQ 200 Immersed boundary FSI IFLOW Sets the compressible flow types EQ 0 Viscous flows laminar EQ 1 Invisid flows IGEOM Sets the geometric dimension EQ 0 2D or 3D will be selected based upon the input mesh and the given boundary conditions EQ 2 two dimensional 2D problem EQ 3 three dimensional 3D problem EQ 101 2D axisymmetric Remarks 1 If the user wants to use the 2D IGEOM 2 or 2D axisymmetric IGEOM 101 solver the mesh should only be distributed in the x y plane one layer with the boundary conditions given only at the x y domain boundaries Otherwise a warning message will be given and the 3D solver will be triggered instead 2 If IGEOM 0 default the code will automatically check the mesh and the given boundary conditions to decide the problem geometry type 2D or 3D 2 18 CESE LS DYNA R7 0 CESE CONTROL SOLVER CESE 3 The 2D axisymmetric case will work only if the 2D mesh and corresponding boundary condi tions are properly defined with the x and y coordinates corresponding to the axial and radia
70. rce of electrical current into solid conductors and the computation of the associated magnetic field electric field as well as induced currents The EM solver is coupled with the structural mechanics solver the Lorentz forces are added to the mechanics equations of motion and with the structural thermal solver the ohmic heating is added to the thermal solver as an extra source of heat The EM fields are solved using a Finite Element Method FEM for the conductors and a Boundary Element Method BEM for the surrounding air insulators Thus no air mesh is necessary As stated above the CHEMISTRY and STOCHASTIC cards are only used in the CESE solv er at this time 1 2 INTRODUCTION LS DYNA R7 0 CESE The keyword CESE provides input data for the Conservation Element Solution Element CESE compressible fluid solver CESE BOUNDARY AXISYMMETRIC OPTION CESE BOUNDARY FSI OPTION CESE BOUNDARY NON REFLECTIVE OPTION CESE BOUNDARY PRESCRIBED OPTION CESE BOUNDARY REFLECTIVE OPTION CESE BOUNDARY SOLID WALL OPTION2 CESE CONTROL LIMITER CESE CONTROL MESH MOV CESE CONTROL SOLVER CESE CONTROL TIMESTEP CESE EOS CAV HOMOG EQUILIB CESE EOS IDEAL GAS CESE INITIAL CESE INITIAL OPTION CESE INITIAL CHEMISTRY CESE INITIAL CHEMISTRY ELEMENT CESE INITIAL CHEMISTRY PART CESE INITIAL CHEMISTRY SET CESE MAT GAS CESE PART Note that when performing a chemistry calculation with the CESE solver ini
71. s NUMSEC Number of sectors ratio of the full circle to the angular extension of the mesh This has to be a power of 2 For example NUMSEC 4 means that the mesh represents one fourth of the full 360 degrees circle LS DYNA R7 0 4 29 EM EM SOLVER BEM EM SOLVER BEM Purpose Define the type of linear solver and pre conditioner as well as tolerance for the EM BEM solve Card 1 1 2 3 4 5 6 7 8 RELTOL STYPE VARIABLE DESCRIPTION RELTOL Relative tolerance for the iterative solvers PCG or GMRES The user should try to decrease this tolerance if the results are not accurate enough More iterations will then be needed MAXITER Maximal number of iterations STYPE Solver type EQ 1 Direct solve the matrices will then be considered as dense EQ 2 Pre Conditioned Gradient method PCG this allows to have block matrices with low rank blocks and thus reduce memory used EQ 3 GMRES method this allows to have block matrices with low rank blocks and thus reduce memory used The GMRES option only works in Serial for now PRECON Preconditioner type for PCG or GMRES iterative solves EQ 0 no preconditioner EQ 1 Diagonal line EQ 2 diagonal block EQ 3 broad diagonal including all neighbor faces EQ 4 LLT factorization The LLT factorization option only works in serial for now USELAST This is used only for iterative solvers PCG or GMRES EQ 1 starts from 0 as initial solution of the l
72. s for time history output Card 1 1 2 3 4 3 6 7 8 Variable SETID Type I Default none Card 2 3 4 1 2 3 4 5 6 7 8 Variable x Y Z Type F F F Default none none none VARIABLE DESCRIPTION SETID Identifier for this point set Called by LSO DOMAIN X Y Z Coordinates of the point As many points as desired can be specified Remarks 1 For the ICFD and CESE solvers the points have to remain inside the fluid mesh For the EM solver the points can be defined inside the conductors or in the air In the latter case the fields will be computed using a Biot Savart type integration 8 4 LSO LS DYNA R7 0 LSO TIME SEQUENCE LSO LSO TIME SEQUENCE Purpose This command provides users with maximum flexibility in specifying exactly the fre quency of output of the LSO data points Each instance of the LSO TIME SEQUENCE command creates a new time sequence with an independent output frequency Card 1 1 Card 2 1 2 3 4 5 6 7 8 Variable DT LCDT LCOPT NPLTC TBEG TEND Type F I I I F F Default 0 0 0 1 0 0 0 0 0 Card 3 1 2 3 4 5 6 7 8 Variable SSID Type I Default none VARIABLE DESCRIPTION SOLVER NAME Selects the solver from which data is output in this time sequence Accepted entries so far cese and icfd DT Time interval between outputs LCDT Op
73. sides the velocity If all three are given only p and P will be used 2 These initial condition will be applied in those elements that have not been assigned a value by CESE_INITIAL_OPTION cards for individual elements or sets of elements LS DYNA R7 0 2 23 CESE CESE CESE_INITIAL 2 24 CESE LS DYNA R7 0 CESE_INITIAL_CHEMISTRY CHEMISTRY CESE_INITIAL_CHEMISTRY Purpose Initializes the chemistry and fluid state in every element of the CESE mesh that has not already been initialized by one of the other CESE_INITIAL_CHEMISTRY cards This is on ly used when chemistry is being solved with the CESE solver Card 1 1 2 3 4 5 6 7 8 Variable CHEMID COMPID Type I Default none Card 2 1 2 3 4 5 6 7 8 Variable Type Default VARIABLE DESCRIPTION CHEMID Identifier of chemistry control card to use COMPID Identifier of chemical composition to use UIC X component of the fluid velocity VIC Y component of the fluid velocity WIC Z component of the fluid velocity RHOIC Initial fluid density PIC Initial fluid pressure TIC Initial fluid temperature HIC Initial fluid enthalpy However when CHEMID refers to a ZND 1 step LS DYNA R7 0 2 25 CESE CHEMISTRY CESE_INITIAL_CHEMISTRY VARIABLE DESCRIPTION reaction card this is the progressive variable degree of combustion 2 26 CESE LS DYNA R7 0 CESE INITIAL CHEMISTRY ELEMENT CESE CESE INITIAL CHEMISTRY ELEME
74. sistency index n Measure of the deviation of the fluid from Newtonian nmin Lower n value limit of the Power Law nmax Upper n value limit of the Power Law LS DYNA R7 0 5 29 ICFD ICFD ICFD_PART ICFD PART OPTION Available options include TITLE Purpose Define parts for this incompressible flow solver The TITLE option allows the user to define an additional line with a HEADING in order to as sociate a name to the part Optional Variable HEADING Type C Default none Remark 1 Card 1 1 2 3 4 5 6 7 8 i VARIABLE DESCRIPTION PID Part identifier for fluid surfaces SECID Section identifier defined with ICFD_SECTION card MID Material identifier defined with the ICFD_MAT card 5 30 ICFD LS DYNA R7 0 ICFD PART VOL CFD ICFD PART VOL OPTION Purpose Define parts for a volume enclosed by surfaces This keyword is used when no volume mesh is provided by the user The solver will build a volume mesh and assign to it the PID de fined by this card Optional Variable HEADING Type C Default none Remark 1 Card 1 1 2 3 4 5 6 7 8 Variable PID SECID MID Type I I I Default none none none Card 2 3 1 2 3 4 5 6 7 8 Variable SPIDI SPID2 SPID3 SPID4 SPIDS SPID6 SPID7 SPID8 Type I I I I I I I I Default none none none none none none none None VARIABLE DESCRIPTION PID Part identifier for f
75. ss on edges and from node to opposite faces Very robust but costly EQ 4 full Lagrangian The mesh moves with the velocity of the flow EQ 11 mesh moves using an implicit ball vertex spring method LIM_ITER Maximum number of linear solver iterations for the ball vertex linear system RELTOL Relative tolerance to use as a stopping criterion for the ball vertex meth od iterative linear solver conjugate gradient solver with diagonal scaling preconditioner ABSTOL Absolute tolerance measure for the size of mesh displacement changes to use as a stopping criterion for the ball vertex iterative linear solver 5 16 ICFD LS DYNA R7 0 ICFD CONTROL OUTPUT ICFD ICFD_CONTROL_OUTPUT Purpose This keyword modifies default values for screen and file outputs related to this fluid solver only Card 1 1 2 3 4 5 6 7 8 Variable MSGL OUTL DTOUT Type I I F Default 0 0 0 VARIABLE DESCRIPTION MSGL Message level EQ 0 only time step information is output EQ 1 first level solver information EQ 2 full output information with details about linear algebra and convergence steps OUTL Output the fluid results in other file formats apart from d3plot EQ 0 only d3plot output EQ 1 output a file with mesh statistics and fluid results in GMV format A directory named output gmv has to be created one level above the executable EQ 2 output a file with mesh statistics and the fluid results in
76. t ID for input voltage when CIRCTYP EQ 2 12 or CIRCTY P EQ 3 It is considered to be oriented as going into the struc tural mesh irrespective of the orientation of the segment Segment set ID for output voltage when CIRCTYP EQ 2 12 or CIRCTY P EQ 3 It is considered to be oriented as going out of the struc tural mesh irrespective of the orientation of the segment Part ID associated to the Source Circuit Only needed if Op tion SOURCE 1 When defining a circuit with an imposed current Type 1 or 11 and in cases of a closed loop geometry torus SIDVIN and SIDVOUT cannot be defined and thus only SIDCURR is necessary 2 When defining a circuit with an imposed tension Type 2 or 12 it is possible to also define SIDCURR This can be useful in circuits where various flow paths are possible for the cur rent in order to force the entire current to go through SIDCURR 4 6 EM LS DYNA R7 0 CIRCUIT Correspondence between circuit type and card entries EM CIRCTYP 1 2 3 12 imp curr imp volt F A t0 F A t0 LCID M M 5 5 R L C VO M 5 R M M SIDCURR M M M SIDVIN M M M M SIDVOUT M M M M M PARTID M mandatory M mandatory with certain exceptions See Remark 1 O optional See Remark 2 not taken into account LS DYNA R7 0 4 7 EM EM CIRCUIT ROGO EM CIRCUIT ROGO Pur
77. tency For example if dimensionless variables are used and should be replaced by the corresponding di mensionless ones LS DYNA R7 0 2 33 CESE CESE CESE_PART CESE_PART Purpose Define CESE solver parts 1 connect CESE material and EOS information Card 1 1 2 3 4 5 6 7 8 VARIABLE DESCRIPTION PID Part identifier MID Material identifier defined by a CESE MAT card EOSID Equation of state identifier defined by a CESE EOS card Remarks 1 Since material coefficients are only used in viscous flows the MID can be left blank for in visid flows 2 34 CESE LS DYNA R7 0 CHEMISTRY CHEMISTRY CHEMISTRY The keyword CHEMISTRY is used to access chemistry databases that include Chemkin based descriptions of a chemical model as well as to select a method of solving the model The key word cards in this section are defined in alphabetical order CHEMISTRY COMPOSITION CHEMISTRY CONTROL 0D CHEMISTRY CONTROL 1D CHEMISTRY CONTROL FULL CHEMISTRY CONTROL TBX CHEMISTRY CONTROL ZND CHEMISTRY DET INITIATION CHEMISTRY PATH Currently the following cards may be used only once in a given problem CHEMISTRY CONTROL 1D CHEMISTRY CONTROL ZND and CHEMISTRY DET INITIATION Also CHEMISTRY CONTROL 00 is only used in a standalone fashion That is it does not involve any other solvers An additional option TITLE may be appended to all CHEMISTRY keywords If this op tion is used t
78. tialization should only be done with the CESE INITIAL CHEMISTRY cards not the CESE INITIAL cards An additional option TITLE may be appended to all CESE keywords If this option is used then an 80 character string is read as a title from the first card of that keyword s input At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks LS DYNA R7 0 2 1 CESE CESE CESE_BOUNDARY_AXISYMMETRIC CESE BOUNDARY AXISYMMETRIC OPTION Available options are PART SET SEGMENT Purpose Define an axisymmetric boundary condition on the axisymmetric axis for the special 2D axisymmetric CESE compressible flow solver For the PART option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SURFPRT Type I Default none For the SET option define the following card Card 1 1 2 3 4 5 6 7 8 Variable SSID Type I Default none 2 2 CESE LS DYNA R7 0 CESE_BOUNDARY_AXISYMMETRIC CESE For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 7 8 Variable NI N2 N3 N4 Type I I I I Default none none none none VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH SURFACE ELEMENT cards SSID Segment set ID N1 N2 Node ID s defining a segment Remarks 1 This boundary condition can only be used on the axisymmetric axis for the 2D axisymm
79. tional load curve ID specifying the time interval between dumps LCOPT Flag to govern behavior of plot frequency load curve LS DYNA R7 0 8 5 LSO LSO LSO TIME SEQUENCE VARIABLE DESCRIPTION EQ 1 At the time each plot is generated the load curve value is added to the current time to determine the next plot time this is the default behavior EQ 2 At the time each plot is generated the next plot time T is computed so that T the current time plus the load curve value at the time T EQ 3 A plot is generated for each ordinate point in the load curve definition The actual value of the load curve is ignored NPLTC DT ENDTIM NPLTC overrides DT specified in the first field TBEG The problem time at which to begin writing output to this time se quence TEND The problem time at which to terminate writing output to this time sequence SETID Output set ID defining the domain over which variable output is to be performed in this time sequence Each SETID refers to the output set defined by a LSO DOMAIN keyword command Remarks 1 If LCDT is nonzero then it is used and DT and NPLTC are ignored If LCDT is zero and NPLTC is non zero then NPLTC determines the snapshot time increment If LCDT and NPLTC are both zero then the minimum non zero time increment specified by DT is used to determine the snapshot times 8 6 LSO LS DYNA R7 0 LSO TIME SEQUENCE LSO LS DYNA 7 0 8 7 LSO
80. ype Card 4 Variable FILE2 Type LS DYNA R7 0 3 3 CHEMISTRY CHEMISTRY CHEMISTRY CONTROL 0D Card 5 1 Variable FILE3 Type A VARIABLE DESCRIPTION ID Identifier for this 0D computation COMPID Chemical composition identifier of composition to use SOLTYP Type of 0D calculation EQ 1 Isochoric EQ 2 Isobaric ERRLIM Error tolerance for the calculation DT Initial time step TLIMIT Time limit for the simulation TIC Initial temperature PIC Initial pressure RIC Initial density EIC Initial internal energy FILE1 Name of the file containing the Chemkin compatible input FILE2 Name of the file containing the chemistry thermodynamics database FILE3 Name of the file containing the chemistry transport properties database 3 4 CHEMISTRY LS DYNA R7 0 CHEMISTRY CONTROL 1D CEMISTRY CONTROL 1D CHEMISTRY Purpose Loads a previously computed one dimensional detonation It is then available for use in the CESE solver for initializing a computation In the product regions this card overrides the ini tialization of CESE INITIAL CHEMISTRY cards Card 1 1 Variable ID Type I Default none Card 2 XYZD none DETDIR none Variable FILEI Type Card 3 Variable FILE2 Type Card 4 Variable FILE3 Type LS DYNA R7 0 3 5 CHEMISTRY CHEMISTRY CHEMISTRY CONTROL 1D
81. yword SET NODE LIST and serves a similar purpose 5 34 ICFD LS DYNA R7 0 MESH MESH MESH The keyword MESH is used to create a mesh that will be used in the analysis So far only tetra hedral or triangular in 2 d elements can be generated The keyword cards in this section are de fined in alphabetical order MESH BL MESH BL SYM MESH EMBEDSHELL MESH INTERF MESH SIZE MESH SIZE SHAPE MESH SURFACE ELEMENT MESH SURFACE NODE MESH VOLUME MESH VOLUME ELEMENT MESH VOLUME NODE An additional option TITLE may be appended to all MESH keywords If this option is used then an 80 character string 15 read as a title from the first card of that keyword s input At present LS DYNA does not make use of the title Inclusion of titles gives greater clarity to in put decks LS DYNA R7 0 6 1 MESH MESH MESH BL MESH BL Purpose Specify the part ID meshed by a boundary layer mesh and how many elements to place in the thickness Card 1 1 2 3 4 5 6 7 8 Variable PID NELTH Type I I Default none None VARIABLE DESCRIPTION PID Part identifier for the surface element NELTH Number of elements normal to the surface in the boundary layer 6 2 MESH LS DYNA R7 0 MESH BL SYM MESH MESH BL SYM Purpose Specify the part ID that will have symmetry conditions for the boundary layer On these surfaces the boundary layer mesh follows the surface tangent

LS-DYNA KEYWORD USER`S MANUAL

Contents

Download Pdf Manuals

Related Search

Related Contents