Home

LS-DYNA®

Contents

1. Card 1 1 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 DEV February 19 2013 revision 2754 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 1e20 VARIABLE DESCRIPTION OWC Indicates the coupling direction to the solver EQ 0 two way coupling Loads and displacements are transferred across the FSI interface and the full non linear problem is solved EQ 1 one way coupling The solid solver transfers displacements 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 transfer any loads to the solid solver but it will continue to deform with the solid 5 14 CFD LS DYNA DEV February 19 2013 revision 2754 ICFD_CONTROL_ LOAD ITCFD ICFD_CONTROL_LOAD Purpose This keyword can reset the body load in the fluid problem provide
2. Type I I I Default none none none e VARIABLE DESCRIPTION SSID Segment Set Id BTYPE EQ 9 The faces of this segment set are eliminated from the BEM cal culations used for example for the rear or side faces of a work piece EQ 10 used by the 2D axisymmetric solver to make the connection be tween 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 EQ 11 See EQ 10 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 DEV February 19 2013 revision 2754 EM_CIRCUIT KM 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 useful in order to save computational time in cases with a low frequency current and where the diffusion 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
3. In the following UUS stands for User Units System and BUS for Burgess Units VARIABLE DESCRIPTION EOSID Cl C2 C3 TEMUNIT VO ID of the EM_EOS C1 constant BUS C2 constant no units C3 constant no units Temperature units EQ 1 temperature in Celsius EQ 2 temperature in Kelvins Reference specific volume VO UUS 4 18 EM LS DYNA DEV February 19 2013 revision 2754 EM_EOS_MEADON KM VARIABLE DESCRIPTION GAMMAO 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 Omatat room temperature o 0 Opurgess 0 Omat Opurgess 8room Remarks 1 The Meadon model is a simplified Burgess model with the solid phase equations only The electrical resistivity is given by ns C1 0 fC 1 where 0 is the temperature V is the specific volume and Vo is the reference specific volume zero pressure solid phase In 1 the volume dependence is given by fe G A if EXPON 1 most materials 2a fe G ANE if EXPON 1 tungsten 2b kG G 2y if EXPON 0 stainless steel 2c fe G 1 if Vo is not defined or equal
4. Card 1 1 2 3 4 5 6 7 8 Variable RELTOL MAXITE FORCON Type F I I Default 1E 2 50 0 VARIABLE DESCRIPTION RELTOL Relative tolerance for the solver The user should try to decrease this toler ance if the results are not accurate enough More iterations will then be needed MAXITER Maximal number of iterations FORCON EQ 0 the code stops with an error if no convergence EQ 1 the code continues to the next time step even if the RELTOL convergence criteria has not been reached LS DYNA DEV February 19 2013 revision 2754 4 33 EM 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_CONJ_HEAT ICFD_BOUNDARY_FREESLIP ICFD_BOUNDARY_FSI ICFD_BOUNDARY_NONSLIP ICFD_BOUNDARY_PRESCRIBED_VEL ICFD_BOUNDARY_PRESCRIBED_PRE 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 DEV February 19 2013 revision 2754 5 1 ICFD ICFD ICFD An additional option
5. Card 1 1 2 3 4 5 6 7 8 Variable XP YP ZP XD YD ZD NUMSEC Type F F F F F F I Default none none none none none none none VARIABLE DESCRIPTION NP X Y Z coordinate of the point ND X Y Z coordinate of direction of the axis 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 4 28 EM LS DYNA DEV February 19 2013 revision 2754 EM_SOLVER_BEM KM 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 Variable RELTOL MAXITE STYPE PRECON USELAS Type I I I I I Default 1E 3 1000 2 2 1 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 PREC
6. 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 DEV February 19 2013 revision 2754 4 13 EM EM 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 VARIABLE DESCRIPTION TSTYPE Time Step type EQ 1 constant time step given in DTCONST 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 L E ee 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 12 dt where 2D e Dis the diffusion coefficient D Hoe e Gis the element electrical conductivity e uois the permeability of free space e 1 is the minimal edge length of the elemen
7. LS DYNA DEV February 19 2013 revision 2754 8 5 LSO LSO LSO_TIME_SEQUENCE VARIABLE LCOPT NPLTC TBEG TEND DOMID Remarks DESCRIPTION Flag to govern behavior of plot frequency load curve EQ 1 At the time each plot is generated the load curve value is add ed 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 com puted 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 DT ENDTIM NPLTC overrides DT specified in the first field The problem time at which to begin writing output to this time sequence The problem time at which to terminate writing output to this time se quence Output set ID defining the domain over which variable output is to be performed in this time sequence Each DOMID refers to the domain identified in an LSO_DOMAIN keyword card 1 IfLCDT 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 DEV February 19 2013 revision 2754 LSO_TIME_SEQUENCE LSO LS DYNA
8. 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 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 7 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 5 6 7 8 Variable MSIZE Radius CenterX CenterY CenterZ Type F F F F F Default none none none none none LS DYNA DEV February 19 2013 revision 2754 6 7 MESH MESH MESH_SIZE if Sname cylinder Cards 2 1 2 3 4 5 6 7 8 Variable MSIZE Radius PminX PminY PminZ PmaxX PmaxY PmaxZ Type F F F F F F F F Default none 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 Z for the point of minimum coordinates Pmax X Y Z X Y Z for the point of maximum coordinates Center X Y Z Coordinates of the sphere center in cases where Snam
9. 2 18 CESE LS DYNA DEV February 19 2013 revision 2754 CESE_CONTROL SOLVER CESE 2 IfIGEOM 0 default the code will automatically check the mesh and the given boundary con ditions to decide the problem geometry type 2D or 3D 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 radial di rections respectively 4 Currently the FSI solver only runs in 3D cases LS DYNA DEV February 19 2013 revision 2754 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 5 6 7 8 Variable IDDT CFL DTINT Type I F F Default 0 0 9 1 0E 3 VARIABLE DESCRIPTION IDDT Sets the time step option EQ 0 Fixed time step size DTINT i e given initial time step size NE O 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 DEV February 19 2013 revision 2754 CESE_EOS_CAV_HOMOG_EQUILIB CESE CESE_EOS_CAV_HOMOG_EQUILIB Purpose Define the coefficients in the equation of state EOS for the homogeneous equilibrium cavitation m
10. LS DYNA DEV February 19 2013 revision 2754 5 19 ICFD ICFD ICFD_CONTROL_SURFMESH 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 surface mesh to be re meshed 5 20 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_CONTROL_ADAPT ICFD 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 purposes this method is automatically used for the first ICFD time step Card 1 1 2 3 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 LS DYNA DEV February 19 2013 revision
11. Card 1 1 Variable MID MTYPE SIGMA EOSID Type I I F I Default none none 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 EQ 0 EQ 2 Conductor carrying a source In these conductors the eddy cur rent problem is solved which gives the actual current density Typically 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 corre spond to the workpiece SIGMA Initial electrical conductivity of the material EOSID ID of the EOS to be used for the electrical conductivity see FEM_EOS card 4 24 EM LS DYNA DEV February 19 2013 revision 2754 EM_MAT_002 EM EM_MAT_002 Purpose Define an electromagnetic material type and properties whose permeability is different than the free space permeability Card 1 1 2 3 4 5 6 8 Variable MID MTYPE SIGMA EOSID MUREL EOSMU Type I I F I F I Default none none none none 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 el
12. DLEN none CFL none TLIMIT none XYZD none DETDIR none Card 2 Variable Type FILE1 Card 3 Variable Type FILE2 Card 4 Variable Type FILE3 3 10 CHEMISTRY LS DYNA DEV February 19 2013 revision 2754 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 FILE 1 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 DEV February 19 2013 revision 2754 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 DIR Type A VARIABLE DES
13. 2754 5 21 ICFD ICFD ICFD_CONTROL_TIME ICFD_CONTROL_TIME Purpose This keyword is used to change the defaults related to time parameters in the fluid problem Card 1 ll 2 3 4 5 6 7 8 Variable TTM DT CFL Type F F F Default 1 E28 0 1 VARIABLE DESCRIPTION TIM 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 auto matically computed based on the CFL condition CFL Scale factor that multiplies DT 5 22 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_ CONTROL _TURBULENCE ITCFD 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 6 7 8 Variable Cel Ce2 Xe Xk Cu Type F F F F F Default 1 44 192 1 3 1 0 0 09 This card is optional and may be used if TMOD 2 Card 2 1 2 3 4 5 6 7 8 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 amp approach EQ 2 L
14. FSI problems It does this with two approaches The first approach solves the compressible flow equations on an Eu lerian 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 corresponding 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 extend the CESE solver The second solver is the incompressible flow solver ICFD that is fully coupled with the solid me chanics solver This coupling permits robust FSI analysis via either an explicit technique when the FST is weak or 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 MESH keywords In addition during the time advancement of the incompressible flow the solution
15. 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 CONE Spray mean cone angle in degrees for hollow cone spray spray cone an gle in degrees for solid cone spray DCONE Injection liquid jet thickness in degrees ANOZ Area of injector AMPO Initial amplitude of droplet oscillation at injector LS DYNA DEV February 19 2013 revision 2754 7 5 STOCHASTIC STOCHASTIC STOCHASTIC_TBX_PARTICLES STOCHASTIC_TBX_PARTICLES Purpose Specify particle and other model details for stochastic PDEs that model embedded particles 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 l Card 3 1 2 3 4 5 6 7 8 Variable XORIG YORIG ZORIG XVEL YVEL ZVEL Type F E F F F F Default none none none 0 0 0 0 0 0 VARIABLE DESCRIPTION 7 6 STOCHASTIC LS DYNA DEV February 19 2013 revision 2754 STOCHASTIC_TBX_PARTICLES STOCHASTIC VARIABLE PCOMB NPRTCL MXCNT PMASS SMR RHOP TICP T_IGNIT INITDST AZIMTH ALTITD CPS CVS HVAP EMISS BOLTZ XO
16. 2 3 4 5 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 DEV February 19 2013 revision 2754 LSO_TIME_SEQUENCE LSO LSO_TIME_SEQUENCE Purpose This command provides users with maximum flexibility in specifying exactly the frequen cy of output of the LSO data points Each instance of the LSO_ TIME_SEQUENCE command cre ates a new time sequence with an independent output frequency Card 1 1 Variable SOLVER_NAME Type C 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 DOMID Type I Default none VARIABLE DESCRIPTION SOLVER Selects the solver from which data is output in this time sequence Ac _NAME cepted entries so far are em cese and icfd DT Time interval between outputs LCDT Optional load curve ID specifying the time interval between dumps
17. February 19 2013 revision 2754 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 DEV February 19 2013 revision 2754 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 axisymmet ric CESE fluid solver LS DYNA DEV February 19 2013 revision 2754 2 3 CESE CESE CESE_ BOUNDARY FSI CESE_BOUNDARY_FSI_OPTION Available options are PART SET SEGMENT Purpose Define an FSI boundary condition for the CESE compressible flow solver This boundary condition must be
18. LSTC All requests to reproduce the contents hereof should be sent to sales Istc com TABLE OF CONTENTS TABLE OF CONTENTS TABCE OF CONTENTS cas da iaa 0 3 INTRODUCTION 0 a iia 1 1 PE ESE 0 A A A A 2 1 CESE BOUNDARY AUS MINE TRIO anti Mellen Aen 2 2 CESE BOUNDARY PS Lc 2 4 CESE BOUNDARY NON REFLECTIVE cerais an da 2 6 C ESE BOUNDARY PRESCRIBED nas lied 2 8 CESE BOUNDARY REFLECTIVE quisas 2 11 CESE_BOUNDAR Y SOLID WALL 2 13 CESE CONTROL EMITE sui a as 2 16 CESE CONTROL MESH MOV ai ia 2 17 IC ESE CONTROL SOLVER utili ptoits 2 18 CESE CONTROL TIMESTEP cia E E A T OE 2 20 CESE BEOS CAV HOMOG EQUI IB omiso 2 21 CESE EOS IDEAL GAS in ia 2 22 AC ESE INITIAL on ato 2 23 CESE INITIAL cr A A N R E 2 24 CESE INITIAL CHEMISTRY as aaa Ea a N dues 2 25 CESE INITIAL CHEMISTRY ELEMENT coup iia 2 26 CESE INITIAL CHEMISTRY PART ontario 2 28 CESE MATGAS io 2 32 TC ESE PAR asta dao 2 33 A A aid 3 1 CHEMISTRY COMPOSITION uri ias 3 2 CHEMES TR CONTROL DD sica ia 3 3 CHEMISTRY CONTROL ID ui 3 5 CHENESTRY CONTROL FUEL dodne oa 3 7 CHEMESTEY CONTROL IND lies 3 9 CHEMISTRY DET INCOSTI N qa aos 3 10 CHEMISTRY PATH iD 3 12 PEM orir r i i E ENE OTOT OOO NR N EE O N O OT OO E E COO eee Pe 4 1 A A EN E E E enn eee E 4 3 EM BOUNDARY spp alain ica E 4 4 LS DYNA DEV February 19 2013 revision 2754 0 3 TABLE OF CONTENTS TABLE OF CONTENTS SI CRU sio 4 5 PIV CIRCUIT ROGO ncini 4 8
19. Me G IA 4 with Mem ANMs em where An ke0 22Lr 8m if k gt 0 Sa An 1 0 0772 2 Om if k 1 tungsten 5b 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 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 Omo 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 37e 5 4 95e 5 9 73e 5 5 35e 5 0 C2 0 113 0 131 0 170 0 465 0 233 0 330 C3 1 145 1 191 1 178 1 226 1 210 0 4133 EXPON 1 1 1 1 1 0 Cs 0 700 0 672 0 673 0 670 0 638 0 089 k 0 964 0 910 1 08 1 0 878 2 LS DYNA DEV February 19 2013 revision 2754 4 17 EM EM EM_EOS_MEADON EM_EOS_MEADON Purpose Define the parameters for a Meadon model giving the electrical conductivity as as a function of the temperature and the density see T J Burgess Electrical resistivity model of metals 4th 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 Variable EOSID Cl C C3 TEMUNI VO GAMMA EXPON Type I F F F I F F I Default none none none none none none none none Card 2 1 2 3 4 5 6 7 8 Variable LGTUNIT TIMUNIT ADJUST Type F F I Default none none none
20. On these surfaces the boundary layer mesh follows the surface tangent Card 1 1 2 3 4 5 6 T 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID Part identifier for the surface element This is the surface with symmetry LS DYNA DEV February 19 2013 revision 2754 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 PID1 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 DEV February 19 2013 revision 2754 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
21. ZND 1 step re action card this is the progressive variable degree of combustion LS DYNA DEV February 19 2013 revision 2754 2 25 CESE CESE CESE INITIAL CHEMISTRY ELEMENT CESE_INITIAL_CHEMISTRY_ELEMENT 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 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 none none none Card 3 1 2 3 4 5 6 7 8 Variable ELE1 ELE2 ELE3 ELE4 ELES ELE6 ELE ELE8 Type I I I I I I I I Default none none none none none none none none VARIABLE DESCRIPTION CHEMID Identifier of chemistry control card to use COMPID Identifier of chemical composition to use 2 26 CESE LS DYNA DEV February 19 2013 revision 2754 CESE INITIAL CHEMISTRY ELEMENT CESE VARIABLE DESCRIPTION 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 re action card this is the progressive v
22. _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 DEV February 19 2013 revision 2754 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 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface in contact with the solid LS DYNA DEV February 19 2013 revision 2754 5 3 ICFD ICFD ICFD BOUNDARY FREESLIP ICFD_BOUNDARY_FREESLIP Purpose Specify the fluid boundary with free 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 free slip boundary condition is applied 5 4 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_BOUNDARY_FSI ICFD ICFD_BOUNDARY_FSI Purpose This keyword defines which fluid surfaces will be considered in contact with the solid sur faces for fluid structure interaction FSI analysis This keyword should not be defined if ICFD _CONTROL_ESI is not defined Card 1
23. applied on a surface of the CESE computational domain that is shared with sur faces 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 4 5 6 7 8 Variable SURFPRT Type I Default none For the SET option define the following card Card 1 il 2 3 4 5 6 7 8 Variable SSID Type I Default none 2 4 CESE LS DYNA DEV February 19 2013 revision 2754 CESE_ BOUNDARY FSI CESE For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 7 8 Variable N1 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 card is also needed for conjugate heat transfer problems with the CESE solver LS DYNA DEV February 19 2013 revision 2754 2 5 CESE CESE CESE BOUNDARY NON REFLECTIVE CESE_BOUNDARY_NON_REFLECTIVE_OPTION Available options are PART SET SEGMENT Purpose Define a passive boundary condition for CESE compressible flows This non reflective boundary condition NBC provides an artificial computational boundary for an open boundary that is passive
24. 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 combustion 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 E EPLUS Q0 GAM 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 QO 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 DEV February 19 2013 revision 2754 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 com putation or it can be subsequently used by the CHEMISTRY _CONTROL_ 1D card in a later run In the product regions this card overrides the initialization of the CESE_INITIAL_CHEMISTRY cards Card 1 6 Variable Type Default ID none COMPID none NMESH none
25. 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 NI 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_SHELL 6 12 MESH LS DYNA DEV February 19 2013 revision 2754 MESH_VOLUME_NODE MESH MESH_VOLUME_NODE Purpose Define a node and its coordinates This keyword is only used in cases where the fluid Vol ume Mesh is provided by the user and is not automatically generated It serves the same purpose as the NODE keyword for solid mechanics Only Tets are supported Card 1 1 2 3 8 Variable NID X BG Type I F F Default none 0 0 VARIABLE DESCRIPTION NID Node ID A unique number with respect to the other volume nodes x coordinate y coordinate z coordinate LS DYNA DEV February 19 2013 revision 2754 6 13 MESH MESH MESH_VOLUME_PART MESH_VOLUME_PART Purpose Associate a volume part number created by a MESH
26. none Card A 1 Variable VARIABLE_NAME Type C VARIABLE DESCRIPTION DOMAIN The type of domain for which LSO output may be generated Accepted en _TYPE tries so far are thist_point 8 2 LSO LS DYNA DEV February 19 2013 revision 2754 LSO_DOMAIN LSO VARIABLE SOLVER _NAME OUTID SET_ID VARIABLE _NAME Remarks DESCRIPTION Selects the solver from which data is output on this domain Accepted en tries 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 Either the name of a single output variable or a variable group See remarks 1 When DOMAIN_TYPE is THIST_POINT the following variable names are accepted Solver Name Variable type EM ICFD CESE magneticField_point Vectors electricField_point velocity_point velocity_point vecpotField_point currentDensity2_point ScalarPotential_point pressure_point pressure_point Scalars temperature_point temperature _point density_ point density_ point Iset_point LS DYNA DEV February 19 2013 revision 2754 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 methods the most common one is to specify points for time history output Card 1 1
27. 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 DEV February 19 2013 revision 2754 CHEMISTRY CONTROL_FULL CHEMISTRY CHEMISTRY_CONTROL_FULL 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 3 4 5 6 7 8 Variable ID ERRLIM Type I F Default none none Card 2 1 Variable FILE1 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 ERRLIM Error tolerance for the full chemistry calculation FILE1 Name of the file containing the Chemkin compatible input LS DYNA DEV February 19 2013 revision 2754 3 7 CHEMISTRY CHEMISTRY CHEMISTRY CONTROL_FULL VARIABLE DESCRIPTION FILE2 Name of the file containing the chemistry thermodynamics database FILE3 Name of the file containing the chemistry transport properties database 3 8 CHEMISTRY LS DYNA DEV February 19 2013 revision 2754 CHEMISTRY_CONTROL_ZND CHEMISTRY CHEMISTRY_CONTROL_ZND Purpose Computes the one dimensional reduced
28. this CESE compressible flow solver Card 1 1 2 3 4 5 6 7 8 Variable IDLMT ALFA BETA EPSR Type I F F F Default 0 0 0 0 0 0 0 Remarks 1 2 3 VARIABLE DESCRIPTION IDLMT Set the stability limiter option See CESE 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 20 0 larger values give more stability but less accuracy Usually a 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 recom mended 3 20 0 larger values give more stability but less accuracy 2 16 CESE LS DYNA DEV February 19 2013 revision 2754 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 sol
29. wants to use dimensionless variables Cv amp Cp above also should be replaced by the cor responding dimensionless ones 2 22 CESE LS DYNA DEV February 19 2013 revision 2754 CESE_ INITIAL CESE CESE_INITIAL Purpose Specify constant initial conditions ICs for flow variables at the centroid of each fluid el ement Card 1 1 2 3 4 5 6 7 8 Variable U V W RH P T Type F F F F F F Default 0 0 0 0 0 1 225 0 0 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 besides 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 DEV February 19 2013 revision 2754 2 23 CESE CESE CESE INITIAL CESE_INITIAL_OPTION Available options include SET ELEMENT Purpose Specify initial conditions for the flow variables at the centroid of each element in a set of elements or at the centroid of a single element Card 1 1 2 3 4 5 6 7 8 Variable EID ESID U V W RHO P T Type I F F F F F F Default none 0 0 0 0 0 0 25 0 0 0 0 Remarks 1 1 1 VARIABLE DESCRIPTION EID ESID Solid element ID EID or solid element
30. 0 Card 2 1 2 3 4 5 6 8 Variable HC TC Beta Type F F F Default 0 0 0 Card 3 1 2 3 4 5 6 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 ba rotropic flows So far only fully incompressible flow is supported RO Flow density VIS Dynamic viscosity ST Surface tension coefficient LS DYNA DEV February 19 2013 revision 2754 5 29 ICFD ICFD ICFD_MAT VARIABLE DESCRIPTION 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 consistency 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 5 30 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_PART ICFD 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 associ ate a name to the part Optional 1 Variable HEADING Type C Default none Rem
31. 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 PID1 PID2 PID3 PID4 PID5 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 DEV February 19 2013 revision 2754 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 interpolation 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 PID1 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 DEV February 19 2013 revision 2754 MESH_SIZE MESH
32. 1 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID PID of the fluid surface in contact with the solid domain LS DYNA DEV February 19 2013 revision 2754 5 5 ICFD ICFD ICFD_BOUNDARY_NONSLIP ICFD_BOUNDARY_NONSLIP Purpose Specify the fluid boundary with a non slip boundary condition Card 1 ll 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 DEV February 19 2013 revision 2754 ICFD_BOUNDARY_PRESCRIBED_MOVEMESH ICFD ICFD_BOUNDARY_PRESCRIBED_MOVEMESH Purpose Allows the node of a fluid surface to translate in certain directions using an ALE approach This is useful in piston type applications or can also be used in certain cases to avoid big mesh de formation Card 1 1 2 3 4 8 Variable PID dofx dofy dofz Type 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 DEV February 19 2013 revision 2754 5 7 ICED ICFD ICFD_BOUNDARY_PRESCRIBED_VEL ICFD_BOUNDARY_PRESCRIBED_VEL Purpose Impose the fluid velocity on the bou
33. 4 CESE LS DYNA DEV February 19 2013 revision 2754 CESE BOUNDARY SOLID WALL CESE VARIABLE DESCRIPTION N1 N2 Node ID s defining a segment LCID Load curve ID to define this solid wall boundary movement If OPTION2 lt BLANK gt Vx Vy Vz velocity vector of the solid wall LCID EQ 0 it is defined by Vx Vy Vz itself LCID NE 0 it will be defined by both of the load curve and Vx Vy Vz Nx Ny Nz are not used in this case If OPTION2 ROTAT Vx Vy Vz x y amp z coordinates of a point on the rotating axis Nx Ny Nz Unit vector of the rotating axis for the 2D case this is not used The rotating frequency Hz is given by the load curve Remarks 1 In this solid wall condition SBC the boundary movement can only be in the tangential direc tion of the wall and should not affect the fluid domain size and mesh during the calculation oth erwise an FSI or moving mesh solver should be used Also this moving SBC only affects vis cous 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 DEV February 19 2013 revision 2754 2 15 CESE CESE CESE_CONTROL_LIMITER CESE_CONTROL_LIMITER Purpose Sets some stability parameters used in the CESE scheme for
34. 5 6 7 8 Variable CIRCID CIRCTYP LCID R F L A C to vo Type I I I F F F F Default none none none none none none none Card 2 1 2 3 4 5 6 7 8 Variable SIDCURR SIDVIN SIDVOUT PARTID Type I I I I Default none none none none VARIABLE DESCRIPTION CIRCID Circuit ID LS DYNA DEV February 19 2013 revision 2754 4 5 EM EM EM_CIRCUIT VARIABLE CIRCTYP LCID R F L A C t0 VO SIDCURR SIDVIN SIDVOUT PARTID DESCRIPTION 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 VO circuit EQ 4 Induced circuit Obsolete EQ 11 Imposed current defined by an amplitude A frequency F and in itial time t I Asin 2mF t to EQ 12 Imposed voltage defined by an amplitude A frequency F and in itial time t V Asin 27F t to 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 normal of the segments To use the opposite orientation use a mi
35. CRIPTION DIR Directory path to add to the search set 3 12 CHEMISTRY LS DYNA DEV February 19 2013 revision 2754 EM EM EM The EM 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 solv ers 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_TABULATED1 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 DEV February 19 2013 revision 2754 4 1 EM EM EM An additional option _TITLE may be appended to all EM 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 DEV February 19 2013 revision 2754 EM_2DAXI KM EM_2DAXI Purpose Sets up the electromagnetism solver as 2D axisymmetric instead of 3D on a given par
36. DEV February 19 2013 revision 2754 EM_ CONTACT RESISTANCE KM Reontact Reonstriction Rfim LS DYNA DEV February 19 2013 revision 2754 4 11 EM EM EM_CONTROL EM_CONTROL Purpose Enable the EM solver and set its options Card 1 1 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 ENDTIM 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 2 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 DEV February 19 2013 revision 2754 EM CONTROL CONTACT KM EM_CONTROL_CONTACT Purpose Turns on the electromagnetism contact algorithms to check for contact between conductors and allow the electromagnetic fields to flow from one conductor to another when detected as in con tact
37. DEV February 19 2013 revision 2754 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 output in the same simulation In addition any number of domains and any number of variables on those do mains 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 DEV February 19 2013 revision 2754 8 1 LSO LSO LSO_DOMAIN LSO_DOMAIN Purpose This command provides a way to output data for specific variables at a frequency defined by LSO_TIME_SEQUENCE and for a specific set of data structure Points nodes elements parts etc See LSO_POINT_SET in a specific binary LSO file Card 1 1 Variable DOMAIN_TYPE Type C Card 2 1 Variable SOLVER_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
38. DEV February 19 2013 revision 2754 8 7 LSO
39. ES Smagorinsky sub grid scale model Cel Ce2 Me k X model constants Xk Cu LS DYNA DEV February 19 2013 revision 2754 5 23 ICFD ITCFD ICFD_CONTROL_ TURBULENCE VARIABLE DESCRIPTION Cs Smagorinsky constant 5 24 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_ DATABASE _DRAG ITCFD ICFD_DATABASE_DRAG Purpose This keyword enables the computation of drag forces over given parts of the model If mul tiple 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 T 8 Variable PID Type I Default none VARIABLE DESCRIPTION PID Part ID of the surface where the drag force will be computed Remarks 1 The file name for this database is icfdragi for instantaneous drag and icfdraga for the drag com puted using average values of pressure and velocities LS DYNA DEV February 19 2013 revision 2754 5 25 ICFD ICFD ICFD DATABASE _ DRAG ICFD_DATABASE_AVERAGE Purpose This keyword enables the computation of average variable values at given time intervals 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 T
40. 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 coordinate z coordinate 6 10 MESH LS DYNA DEV February 19 2013 revision 2754 MESH_VOLUME MESH MESH_VOLUME Purpose This keyword defines the volume space that will be meshed The boundaries of the volume are the surfaces defined by MESH_SURFACE_ELEMENT The list of surfaces have to be non overlapping and should not leave any gaps or open spaces between the surface boundaries 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 Type I Default none Cards 2 3 4 the next card terminates the input 1 2 3 4 5 6 7 8 Variable PID1 PID2 PID3 PID4 PID5 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 PIDn Part IDs for the surface elements that are used to define the volume LS DYNA DEV February 19 2013 revision 2754 6 11 MESH MESH MESH_VOLUME_ELEMENT MESH_VOLUME_ELEMENT Purpose Specify a set of volume elements for the
41. 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 DEV February 19 2013 revision 2754 CESE_BOUNDARY_NON_REFLECTIVE CESE For the SEGMENT option define the following card Card 1 1 2 3 4 8 Variable N1 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 is usually imposed on an open surface that is far from the main dis turbed flow the farther away the better i e the flow on that boundary surface should be al most 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 DEV February 19 2013 revision 2754 2 7 CESE CESE CESE BOUNDARY PRESCRIBED CESE_BOUNDARY PRESCRIBED_OPTION Available options include PART SET SEGMENT Purpose For the CESE compressible flow solver impose flow variable s velocity density pres sure temperature et
42. LS DYNA KEYWORD USER S MANUAL VOLUME III Multi Physics Solvers February 19 2013 revision 2754 R7 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 Livermore Software Technology Corporation 7374 Las Positas Road 1740 West Big Beaver Road Livermore California 94551 Suite 100 Tel 925 449 2500 M Fax 925 449 2507 Troy Michigan 48084 Email sales Istc com Tel 248 649 4728 M Fax 248 649 6328 Website www lstc com Disclaimer Copyright 1992 2010 Livermore Software Technology Corporation All Rights Reserved LS DYNA LS OPT and LS PrePost are registered trademarks of Livermore Software Technolo gy Corporation in the United States All 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 notice The information and examples included herein are for illustrative purposes only and are not intended 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 result ed from the use of this manual Any reproduction in whole or in part of this manual is prohibited without the prior written approval of
43. LS DYNA DEV February 19 2013 revision 2754 CESE CESE_BOUNDARY PRESCRIBED VARIABLE IDCOMP LC_U LC_V LC_W LC_RHO LCP LC_T SF_U SF_V SF_W SF_RHO SF_P SF_T Remarks DESCRIPTION For inflow boundaries in problems involving chemical reacting flows 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_W 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 Oneach centroid or set of centroids the variables x vel y vel z vel p P T that are given val ues must be consistent and make the model well posed i e be such that the solution of the model exists is unique and physical 2 Ifany of the load curves are 0 the corresponding variable will take the constant value of the cor responding scale factor For instance if LC_RHO 0 then the constant val
44. NA DEV February 19 2013 revision 2754 2 29 CESE CESE CESE INITIAL CHEMISTRY SET 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 SETID CHEMID COMPID Type I I I Default none none none Card 2 1 2 3 4 5 6 7 8 Variable UIC VIC WIC RHOIC PIC TIC HIC Type F E F F F F F Default none none none none none none none VARIABLE DESCRIPTION SETID Identifier of the CESE element set 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 TIC Initial fluid temperature 2 30 CESE LS DYNA DEV February 19 2013 revision 2754 CESE INITIAL CHEMISTRY SET CESE VARIABLE DESCRIPTION HIC Initial fluid enthalpy However when CHEMID refers to a ZND 1 step re action card this is the progressive variable degree of combustion LS DYNA DEV February 19 2013 revision 2754 2 31 CESE CESE CESE_MAT_GAS CESE_MAT_GAS P
45. ON 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 linear system EQ 1 starts from previous solution normalized by the rhs change LS DYNA DEV February 19 2013 revision 2754 4 29 EM EM EM_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 4 30 EM LS DYNA DEV February 19 2013 revision 2754 EM_SOLVER_BEMMAT KM 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 Variable MATID RELTOL Type I F Default none pe ons 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 e 2 for the Q matrix The user should try to decrease these tolerances if the results are not ac
46. Q 1 Face to face PSIDM Master part set ID PSIDS Slave part set ID LS DYNA DEV February 19 2013 revision 2754 4 9 EM EM EM_CONTACT_RESISTANCE EM_CONTACT_RESISTANCE Purpose Calculate the contact resistance of a previously defined EM contact in EM_CONTACT Card 1 1 Variable CRID CONTID Pprobe Psub Poxi 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 Pprobe Probe resistivity Pprob Psub Substrate resistivity Psp Poxi Film resistivity Poxy FACEACTE Scale factor on the constriction area when calculating the constriction re sistance 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 absolute 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 Contacts Pprob Psub Reonstriction CO y faceActe x ContactArea Poxy R firm gt rim faceAfilm x ContactArea 4 10 EM LS DYNA
47. RIG YORIG ZORIG XVEL YVEL ZVEL 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 ex plosive 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 explo sive Y coordinate of the origin of the initial reference frame of the TBX explo sive Z coordinate of the origin of the initial reference frame of the TBX explo sive 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 LS DYNA DEV February 19 2013 revision 2754 7 7 STOCHASTIC STOCHASTIC STOCHASTIC_TBX_PARTICLES Remarks 1 Ifradiation heat transfer is being modeled then EMISS and BOLTZ are required 7 8 STOCHASTIC LS DYNA
48. TION2 For OPTION the choices are PART SET SEGMENT For OPTIONZ2 the choices are lt BLANK 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 a 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 DEV February 19 2013 revision 2754 2 13 CESE CESE CESE_BOUNDARY_SOLID_WALL For the SET option define the following card 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 Nl 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 Nx Ny Nz Type F F F Default 0 0 0 0 0 0 Remarks 3 3 3 VARIABLE DESCRIPTION SURFPRT Surface part ID referenced in MESH_SURFACE_ELEMENT cards SSID Segment set ID 2 1
49. Variable EOSID LCID Type I I Default none none 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 conductivities at very low and very high temperatures to avoid bad extrapolations of the conductivity 1f the temperature gets out of the load curve bounds 4 22 EM LS DYNA DEV February 19 2013 revision 2754 EM_EXTERNAL_FIELD EM EM_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 8 Variable FIELDID FTYPE FDEF LCIDX LCIDY LCIDZ Type I I F I I I Default 0 0 0 0 0 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 time LS DYNA DEV February 19 2013 revision 2754 4 23 EM EM EM_MAT_001 EM_MAT_001 Purpose Define the electromagnetic material type and properties for a material whose permeability equals the free space permeability
50. _VOLUME card with the part num ber of a part card from a selected solver designated by the SOLVER field REQ 1 2 3 5 6 7 8 Variable VOLPRT SOLPRT SOLVER Type I A Default VARIABLE DESCRIPTION VOLPRT Part ID of a volume part created by a MESH_VOLUME card SOLPRT Part ID of a part created using SOLVER s part card SOLVER Name of a solver using a mesh created with MESH cards 6 14 MESH LS DYNA DEV February 19 2013 revision 2754 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 particles 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 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 DEV February 19 2013 revision 2754 7 1 STOCHASTIC STOCHASTIC STOCHASTIC_SPRAY_PARTICLES STOCHASTIC_SPRAY_ PARTICLES Purpose Specify particle and other model details for spray modeling using stochastic PDEs that ap proximate such proces
51. 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 solver at this time 1 2 INTRODUCTION LS DYNA DEV February 19 2013 revision 2754 CESE CESE 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_ OPTION _ 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 initialization 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 DEV
52. ariable degree of combustion ELEI User element numbers to initialize LS DYNA DEV February 19 2013 revision 2754 2 27 CESE CESE CESE INITIAL CHEMISTRY PART 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 by CESE_INITIAL CHEMISTRY _ ELEMENT or CESE _INITIAL_CHEMISTRY_SET 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 PARTID CHEMID COMPID Type I I I Default none 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 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 TIC Initial fluid temperature 2 28 CESE LS DYNA DEV February 19 2013 revision 2754 CESE INITIAL CHEMISTRY PART CESE VARIABLE DESCRIPTION HIC Initial fluid enthalpy However when CHEMID refers to a ZND 1 step re action card this is the progressive variable degree of combustion LS DY
53. ark 1 Card 1 1 2 3 4 5 6 a 8 Variable PID SECID MID Type I I I Default none none none VARIABLE DESCRIPTION PID Part identifier for fluid surfaces SECID Section identifier defined with the ICFD_ SECTION card MID Material identifier defined with the ICFD_MAT card LS DYNA DEV February 19 2013 revision 2754 5 31 ICFD ICFD ICFD_PART_VOL ICFD_PART_VOL_ OPTION Available options include TITLE 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 defined by this card The TITLE option allows the user to define an additional line with a HEADING in order to associ ate a name to the part Optional 1 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 5 32 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_PART_VOL ICFD Card Dig doe 1 2 3 4 5 6 7 8 Variable SPID1 SPID2 SPID3 SPID4 SPID5 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 fluid volumes SECID Section identifier defined by the ICFD_SECTION card MID Material identifier SPIDn Part IDs for
54. c at the centroid of the element connected with this boundary segment OP TION 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 DEV February 19 2013 revision 2754 CESE_BOUNDARY_PRESCRIBED For the SEGMENT option define the following card Card 1 1 2 3 4 5 6 Variable NI N1 N3 N4 IDCOMP Type I I I I I Default none none none none none Card 2 1 2 3 4 5 6 Variable LC_U LC_V LC W LC RHO LC P LC IT Type I I I I I I Remarks 162 1 2 172 142 172 1 2 Card 3 1 2 3 4 5 6 Variable SESU SF_V SF_W SF_RHO SF_P SES 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 A surface part ID referenced in MESH_SURFACE_ELEMENT cards SSID Segment set ID N1 N2 Node ID s defining a segment
55. curate enough More memory will then be needed LS DYNA DEV February 19 2013 revision 2754 4 31 EM EM EM_SOLVER_FEM EM_SOLVER_FEM Purpose Define some parameters for the EM_FEM solver Card 1 1 2 3 4 5 6 7 8 Variable RELTOL MAXITE STYPE PRECON USELAS Type I I I I I Default 1E 3 1000 2 2 1 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 EQ 2 Conditioned Gradient Method PCG PRECON Preconditioner type for PCG EQ 0 no 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 EQ 1 starts from previous solution normalized by the 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 4 32 EM LS DYNA DEV February 19 2013 revision 2754 EM_SOLVER_FEMBEM EM EM_SOLVER_FEMBEM Purpose Define some parameters for the coupling between the EM_FEM and EM_BEM solvers
56. d by keyword LOAD _BODY It is useful in problems where the gravity acceleration may be neglected for the fluid prob lem but not for the solid Card 1 1 2 3 4 5 6 T 8 Variable ABL Type I Default 1 VARIABLE DESCRIPTION ABL EQ 0 the body load provided in LOAD_BODY is reset to zero only for the fluid analysis LS DYNA DEV February 19 2013 revision 2754 5 15 ICFD ICFD ICFD_CONTROL_MESH 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 E Default 1 41 VARIABLE DESCRIPTION MGSF Mesh Growth Scale Factor Specifies the maximum mesh size that the vol ume 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 values 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 5 16 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_ CONTROL_MESH_MOV ICFD ICFD_CONTROL_MESH_MOV Purpose With this keyword the user can choos
57. e 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 DEV February 19 2013 revision 2754 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 E 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 DEV February 19 2013 revision 2754 ICFD BOUNDARY PRESCRIBED TEMP ITCFD ICFD_BOUNDARY_PRESCRIBED_TEMP Purpose Impose a fluid temperature 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 ile 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
58. e 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 DEV February 19 2013 revision 2754 MESH_SURFACE_ ELEMENT MESH MESH_ SURFACE ELEMENT Purpose Specify a set of surface elements quadrilateral or triangular 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 differ ent mesh sizes inside the volume see card MESH_SIZE Card 1 1 2 3 4 5 6 7 8 Variable EID PID NI 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_SURFACE _ELEMENTS cards PID Part ID A unique part identification number N1 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 trian gular face N3 N4 In 2 d N2 N3 N4 LS DYNA DEV February 19 2013 revision 2754 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
59. e 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 category 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 tetrahedral element meshes from bounding surface meshes is one of these tools Another is a new data output mecha nism for a limited set of variables from the solvers in this manual This mechanism is accessed through LSO keyword cards The CESE solver is a compressible flow solver based upon the Conservation Element Solution Ele ment 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 many non traditional features including a unified treatment of space and time the introduction of separate conservation elements CE and solution elements SE and a novel shock capturing strategy without using 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 interaction
60. e the type of algorithm for mesh movement Card 1 1 2 3 4 8 Variable MMSH PM RELTOL ABSTOL ieee _ITER 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 el ement sizes as stiffness EQ 3 mesh uses a Laplacian smoothing with stiffness 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 sys tem RELTOL Relative tolerance to use as a stopping criterion for the ball vertex method iterative linear solver conjugate gradient solver with diagonal scaling pre conditioner ABSTOL Absolute tolerance measure for the size of mesh displacement changes to use as a stopping criterion for the ball vertex iterative linear solver LS DYNA DEV February 19 2013 revision 2754 5 17 ICFD ICFD ICFD_CONTROL_OUTPUT 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 lev
61. e to also define SIDCURR This can be useful in circuits where various flow paths are possible for the current in order to force the entire current to go through SIDCURR LS DYNA DEV February 19 2013 revision 2754 4 7 EM EM EM_CIRCUIT_ROGO EM_CIRCUIT_ROGO Purpose Define Rogowsky coils to measure a global current vs time through a segment set or a node set Card Format 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 Anode set allows computing only a surface current An ASCII file rogo_xxx with xxx representing the rogold is generated for each EM_CIRCUIT_ROGO card giving the value of the current vs time 4 8 EM LS DYNA DEV February 19 2013 revision 2754 EM_CONTACT EM EM_CONTACT Purpose Define an electromagnetic contact between two sets of parts Used with the EM _CONTACT_RESISTANCE card Card 1 2 3 4 5 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 E
62. ection case IPULSE 0 LS DYNA DEV February 19 2013 revision 2754 7 3 STOCHASTIC STOCHASTIC STOCHASTIC_SPRAY_PARTICLES VARIABLE IDFUEL RHOP TIP PMASS PRTRTE STRINJ DURINJ XORIG YORIG ZORIG SMR VELINJ DRNOZ DTHNOZ TILTXY DESCRIPTION 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 C2H50H EQ 6 Gasoline C8H18 EQ 7 Jet A C12H23 EQ 8 Kerosene C12H23 EQ 9 Methanol CH3OH EQ 10 N_dodecane C12H26 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 center 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 counter clockwise from there 7 4 STOCHASTIC LS DYNA DEV February 19 2013 revision 2754 STOCHASTIC_SPRAY _PARTICLES STOCHASTIC VARIABLE DESCRIPTION TILTXZ Inclination angle in degrees of the injection in the x z plane where
63. ectromag netism behavior as EQ 0 EQ 2 Conductor carrying a source In these conductors the eddy cur rent problem is solved which gives the actual current density Typically 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 corre spond to the workpiece SIGMA Initial electrical conductivity of the material EOSID ID of the EOS to be used for the electrical conductivity see FEM_EOS card MUREL Relative permeability Is the ratio of the permeability of a specific medium to the permeability of free space u 4 Ho EOSMU ID of the EOS to be used to define the behavior of u by an equation of state Note if EOSMU is defined MUREL will be used for the initial value on ly LS DYNA DEV February 19 2013 revision 2754 4 25 EM EM EM OUTPUT 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 MESH MEM TIMING Type I I I I I I I Default 0 0 0 0 0 0 0 VARIABLE DESCRIPTION MATS Level of matrix assembly output to the screen EQ 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 matr
64. el 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 con vergence 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 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 5 18 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_CONTROL_PARTITION ICFD ICFD_CONTROL_PARTITION Purpose This keyword changes the default option for the partition in MPP thus it is only valid in MPP Card 1 Variable Type Default PTECH VARIABLE PTECH DESCRIPTION Indicates the type of partition the library Metis is used partition along the axis with higher aspect ratio partition along X axis partition along Y axis partition along Z axis
65. gar 2 UO ie GO ee eee AAA A T A 4 9 EM _ CONTACT RESISTANCE nta ia 4 10 ENL CONTRO id 4 12 EM CONTROL CONTACT us a a 4 13 ENM CONTROLE TIMESTER aia dni 4 14 EM EOS BURGESS sorry A aan 4 15 TENEROS MESDON cute bes od pd 4 18 PE EOS PERMEABILITY on id 4 21 SEMI OS TEE A EI Ll ci T 4 22 EM EXTERNAL FIELD oro AAA 4 23 ENC WA O iio 4 24 EM MAT OU 0 aii is ii id 4 25 UL OUTPUT asada 4 26 PEM ROTATION AXIS erriren e chiawhi EAN RAET RTN E AE REEE 4 28 TEN SOLVER BEN condir n E O ER 4 29 EM SOLVER BE NIMIA Tic ida 4 31 EM SOLVER FEM A 4 32 EM SOLVER FEMBEM sissien nrar r Aria 4 33 PACED BOUNDARY CONI HEA Tao is 5 3 CFB BOUNDARY FREESLIP cosita o o 5 4 ICED BOUNDARY PS Los a as 5 5 ICFD BOUNDARY NON SUID latinas 5 6 ICFD BOUNDARY PRESCRIBED MOVEMEBES Husos ir s 5 7 ICFD BOUNDARY PRESCRIBED WE Luisa rata nin 5 8 TED BOUNDARY PRESCRIBED PRE ua 5 9 ICFD BOUNDARY FLUX TEMP src Aire 5 10 AICED BOUNDARY PRESCRIBED TEMP sc EE ae EE 5 11 ICED CONTROL ADAP T oi ri 5 12 ACP CUNTROL ALAR SIE sind 5 13 ICED CONTROL POL A A 5 14 ACE DO CONTROL LOAD estra poda 5 15 KFD CONTROL MESA sti A ARS 5 16 ICE LUNTROL MESA MUS A IA AA 5 17 ICFD CONTROL OUTPUT asia cian 5 18 PTE CONTROL PARTITION us sida 5 19 0 4 TABLE OF CONTENTS LS DYNA DEV February 19 2013 revision 2754 TABLE OF CONTENTS ICPD CONTROL SURENES H uscar dat 5 20 ICEDICONTROL SPLUD 0 ii 5 21 E PDOCONTROL TIME rere eee re
66. he file name for this database is icfdavg dat with the different averaged variable values cop 1ed in a ASCII format 5 26 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_DEFINE_POINT ICFD ICFD_DEFINE_POINT Purpose This keyword defines a point in space that could be used for multiple purposes Card 1 1 2 3 4 8 Variable POID X x 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 LS DYNA DEV February 19 2013 revision 2754 5 27 ICFD ICFD ICFD_INITIAL 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 none 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 5 28 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_MAT ICFD ICFD_MAT Purpose Specify physical properties for the fluid material Card 1 1 2 3 4 5 6 8 Variable MID FLG RO VIS ST THD Type I I F F F F Default none 1 none 0 0
67. horic 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 FILE 1 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 DEV February 19 2013 revision 2754 CHEMISTRY_CONTROL_1D CHEMISTRY CHEMISTRY_CONTROL_1D 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 initializa tion of the CESE_INITIAL_CHEMISTRY cards Card 1 1 2 3 4 5 6 7 8 Variable ID XYZD DETDIR Type I F I Default none none none Card 2 1 Variable FILE1 Type A Card 3 1 Variable FILE2 Type A Card 4 1 Variable FILE3 Type A LS DYNA DEV February 19 2013 revision 2754 3 5 CHEMISTRY CHEMISTRY CHEMISTRY_CONTROL_1D 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
68. ifier defined by a CESE_EOS_card Remarks 1 flows Since material coefficients are only used in viscous flows the MID can be left blank for invisid LS DYNA DEV February 19 2013 revision 2754 2 33 CESE 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 keyword 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_OD 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 option is used then an 80 character string is read as a title from the first card of that keyword s input At pre sent LS DYNA does not make use of the title Inclusion of titles gives greater clarity to input decks LS DYNA DEV February 19 2013 revision 2754 3 1 CHEMISTRY CHEMISTRY CHEMISTRY_COMPOSITION CHEMISTRY_COMPOSITION Purpose Provides a general way to specify a chemical com
69. imilar to the one used by the keyword SET_NODE_LIST and serves a similar purpose MESH MESH MESH The keyword MESH is used to create a mesh that will be used in the analysis So far only tetrahe dral or triangular in 2 d elements can be generated The keyword cards in this section are defined 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 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 DEV February 19 2013 revision 2754 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 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 DEV February 19 2013 revision 2754 MESH_BL_SYM MESH MESH_BL_SYM Purpose Specify the part ID that will have symmetry conditions for the boundary layer
70. 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 aniso LS DYNA DEV February 19 2013 revision 2754 1 1 INTRODUCTION INTRODUCTION INTRODUCTION tropic 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 propaga tion 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 Forming bending or welding induced heating ring expansions and so forth The EM module allows the introduction of a source 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 me chanics 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
71. 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 DEV February 19 2013 revision 2754 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 maximum 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 CFD LS DYNA DEV February 19 2013 revision 2754 ICFD_CONTROL_ADAPT_SIZE ICFD ICFD_CONTROL_ADAPT_SIZE Purpose This keyword controls the re meshing of elements taking into account the element quality and distortion overwriting the default which only checks for inverted elements
72. ix assembly output to the messag file EQ 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 SOLS Level of solver output on the screen EQ 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 EQ 0 No output EQ 1 Global information at each FEM iteration EQ 2 Detailed information at each FEM iteration 4 26 EM LS DYNA DEV February 19 2013 revision 2754 EM_OUTPUT EM VARIABLE DESCRIPTION MESH Controls the output of the mesh data to the d3hsp file EQ 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 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 LS DYNA DEV February 19 2013 revision 2754 4 27 EM EM EM_ROTATION_AXIS 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
73. ndary 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 l 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 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD BOUNDARY PRESCRIBED PRE ITCFD ICFD_BOUNDARY_PRESCRIBED_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 l 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 th
74. nus sign in front of the segment set id 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 set ID for input voltage when CIRCTYP EQ 2 12 or CIRCTYP EQ 3 It is considered to be oriented as going into the structural mesh irrespective of the orientation of the segment Segment set ID for output voltage when CIRCTYP EQ 2 12 or CIRCTYP EQ 3 It is considered to be oriented as going out of the structur al mesh irrespective of the orientation of the segment Part ID associated to the Source Circuit Only needed if Option SOURCE 4 6 EM LS DYNA DEV February 19 2013 revision 2754 EM_CIRCUIT EM Correspondence between circuit type and card entries CIRCTYP 1 2 3 11 12 imp curr imp volt R L C F A t0 F A t0 LCID M M R L C VO M VA t0 M M SIDCURR M O M M O SIDVIN M 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 Remarks 1 When defining a circuit with an imposed current Type 1 or 11 and in cases of a closed loop ge ometry torus SIDVIN and SIDVOUT cannot be defined and thus only SIDCURR is neces sary 2 When defining a circuit with an imposed tension Type 2 or 12 1t is possibl
75. odel Card 1 1 2 3 4 5 6 7 8 Variable EOSID Pap Pig yap Aliq Map Hig Psatvap 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 Doy density of the saturated vapor Pig density of the saturated liquid vap sound speed of the saturated vapor tig sound speed of the saturated liquid Hap dynamic viscosity of the vapor Hig dynamic viscosity of the liquid Psatvap pressure of the 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 DEV February 19 2013 revision 2754 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 flu 1d solver Card 1 1 2 3 5 6 7 8 Variable EOSID Cv Cp Type I F F Default none TS 1004 5 VARIABLE DESCRIPTION EOSID Equation of state identifier Cv 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
76. position via a list of species mole num bers in the context of a chemkin database model REQ 1 1 2 3 4 5 6 7 8 Variable ID Type I Default none REQ Lora 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 DEV February 19 2013 revision 2754 CHEMISTRY_CONTROL_0D CHEMISTRY CHEMISTRY_CONTROL_0D 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 0e 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 1 Variable FILE1 Type A Card 4 1 Variable FILE2 Type A LS DYNA DEV February 19 2013 revision 2754 3 3 CHEMISTRY CHEMISTRY CHEMISTRY_CONTROL_0D Card 5 1 Variable FILE3 Type A VARIABLE DESCRIPTION ID Identifier for this OD computation COMPID Chemical composition identifier of composition to use SOLTYP Type of OD calculation EQ 1 Isoc
77. s for UUS relative to meter i e 1 e 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 conductivity defined in EM_MAT card Omat at room temperature o 0 OBurgess 9 Omat OBurgess Oroom Remarks 1 The Burgess model gives the electrical resistivity vs temperature and density for the solid phase liquid 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 com puted by 1 V 3 ate Om Omno 7 vo Da Vo Vo e IfT lt 0m solid phase model The solid phase electrical resistivity corresponds to the Meadon model ns C1 C0 C 1 where 0 is the temperature V is the specific volume and Vo is the reference specific volume zero pressure solid phase In 1 the volume dependence is given by fe G 2y 1 if EXPON 1 most materials 2a fe G ute if EXPON 1 tungsten 2b f GO if EXPON 0 stainless steel 2c with Y Yo ot 7 3 4 16 EM LS DYNA DEV February 19 2013 revision 2754 EM_EOS_BURGESS EM e If T gt On liquid phase model n
78. ses 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 8 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 8 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 8 Variable XORIG YORIG ZORIG SMR VELINJ DRNOZ DTHNOZ Type F F F F F F F 7 2 STOCHASTIC LS DYNA DEV February 19 2013 revision 2754 STOCHASTIC_SPRAY_ PARTICLES STOCHASTIC Card 4 6 Card N 1 2 3 4 5 6 7 8 Variable TILTXY TILTXZ CONE DCONE ANOZ AMPO Type F F F F F F Default VARIABLE DESCRIPTION INJDIST 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 IBRKUP Type of particle breakup model EQ 0 off no breakup EQ 1 TAP EQ 2 KHRT ICOLLDE Turn collision modeling on or off IEVAP Turn evaporation on or off IPULSE Type of injection EQ 0 continuous injection EQ 1 sine wave EQ 2 square wave LIMPRT Upper limit on the number of parent particles modeled in this spray This is not used with the continuous inj
79. set ID ESID U V W X y Z velocity components respectively RHO density P pressure T temperature Remarks 1 Usually only two of p P amp T are needed to be specified besides the velocity If all three are given only p and P will be used 2 The priority of this card is higher than CESE_INITIAL i e if an element is assigned a initial value by this card CESE_INITIAL will not apply to it anymore 2 24 CESE LS DYNA DEV February 19 2013 revision 2754 CESE_INITIAL_CHEMISTRY CHEMISTRY CESE_INITIAL_CHEMISTRY Purpose Initializes the chemistry and fluid state in every element of the CESE mesh that has not al ready been initialized by one of the other CESE_INITIAL_CHEMISTRY cards This is only used when chemistry is being solved with the CESE solver Card 1 1 2 El 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 UIC VIC WIC RHOIC PIC TIC HIC Type F F F F F F F Default none none none none none none none 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
80. t in order to save computational time as well as memory The electromagnetism is solved in 2D on a given cross section of the part defined by a segment set with asymmetry 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 compat ible with the symmetry i e each node in the part needs to be the child of a parent node on the seg ment set by a rotation around the axis Only the conductor parts with a EM_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 DEV February 19 2013 revision 2754 4 3 EM EM EM BOUNDARY EM_ BOUNDARY Purpose Define some boundary conditions for the electromagnetism problems Card 1 1 2 3 4 5 6 7 8 Variable SSID BIYPE LIYPE
81. t minimal size of the element 4 14 EM LS DYNA DEV February 19 2013 revision 2754 EM_EOS_BURGESS KM EM_EOS BURGESS Purpose Define the parameters for a Burgess model giving the electrical conductivity as as a function of the temperature and the density see T J Burgess Electrical resistivity model of metals 4th 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 Variable EOSID VO GAMMA THETA LF Cl 09 C3 Type I F F F F F F F Default none none none none none none none none Card 2 1 2 3 4 5 6 7 8 Variable C4 K EXPON LGTUNIT TIMUNIT TEMUNI 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 VO Reference specific volume Vo UUS GAMMAO Reference Gruneisen value yo no units THETA Reference melting temperature Omo in eV BUS LF Latent heat of fusion Lr in kJoule mol BUS Cl C1 constant BUS C2 C2 constant no units C3 C3 constant no units LS DYNA DEV February 19 2013 revision 2754 4 15 EM EM EM_EOS_BURGESS VARIABLE DESCRIPTION C4 C4 constant no units K Parameter k no units EXPON Exponent in equations 2 see remarks LGTUNIT Length unit
82. the surface elements that define the volume mesh PID 3 4 N H PARIDAS H fav fav PID 3 SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS ssss ICFD PART VOL SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS PART ID 5 is defined by the surfaces that enclose it ICFD PART VOL A N A O O NS pid secid mid 5 1 1 A E ESE ESRAS A E a a de pidl pid2 pid3 pid4 pid5 pid6 pid7 pid8 1 2 3 4 LS DYNA DEV February 19 2013 revision 2754 5 33 ICFD ICFD ICFD_SECTION ICFD_SECTION Purpose Define a section for the incompressible flow solver Card 1 1 2 3 4 5 6 7 8 Variable SID Type I Default none VARIABLE DESCRIPTION SID Section identifier 5 34 ICFD LS DYNA DEV February 19 2013 revision 2754 ICFD_SET_NODE ICFD 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 7 8 Variable SID PID Type I I Default none none Card 2 1 2 3 4 5 6 7 8 Variable NID1 NID2 NID3 NID4 NIDS5 NID6 NID7 NID8 Type I I I I I I I I Default none none none none none none none none VARIABLE DESCRIPTION SID Set ID PID Associated Part ID NIDI Node ID n Remarks 1 The convention is the s
83. to 0 2d In this last case only EOSID C1 C2 C3 TEMUNIT TIMUNIT and LGTUNIT need to be de fined with Y Y0 3 LS DYNA DEV February 19 2013 revision 2754 4 19 EM EM EM EOS MEADON The following table reports some sets of parameters given by Burgess in his paper 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 C2 0 113 0 131 0 170 0 465 0 233 0 330 C3 1 145 1 191 1 178 1 226 1 210 0 4133 EXPON 1 1 1 1 1 0 4 20 EM LS DYNA DEV February 19 2013 revision 2754 EM_ EOS PERMEABILITY KM EM_EOS_PERMEABILITY Purpose Define the parameters for the behavior of a material s permeability Card 1 1 2 3 4 5 6 7 8 Variable EOSID EOSTYPE LCID Type I I I Default none none none VARIABLE DESCRIPTION EOSID ID of the EM_EOS EOSTYPE Define the type of EOS EQ 1 Permeability defined by a H function of B curve H B u EQ 2 Permeability defined by a B function of H curve B uH LCID Load curve ID LS DYNA DEV February 19 2013 revision 2754 4 21 EM EM EM_EOS_TABULATED1 EM_EOS_TABULATED1 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
84. ue of the density for this boundary condition will be SF_RHO 2 10 CESE LS DYNA DEV February 19 2013 revision 2754 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 of the computation al 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 DEV February 19 2013 revision 2754 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 N1 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 2 This boundary condition has the same effect as a solid wall boundary condition for inviscid flows 2 12 CESE LS DYNA DEV February 19 2013 revision 2754 CESE BOUNDARY SOLID_ WALL CESE CESE_BOUNDARY_SOLID_ WALL_OPTION1_OP
85. urpose Define the fluid gas properties in a viscous flow for the CESE solver Card 1 1 2 3 4 5 6 7 8 Variable MID CIV 62 Pmd Type I F F F 1 458E Default none 6 110 4 0 72 VARIABLE DESCRIPTION MID Material identifier C1 C2 Two coefficients in the Sutherland s formula for viscosity i e u C T 3 T C2 where C and C are constants for a given gas For example for air at moderate temperatures C 1 458x 10 kg m s K 2 Co 110 4 K Prnd The Prandtl Number used to determine the coefficient of thermal conduc tivity Itis approximately constant for most gases For air at standard con ditions Prnd 0 72 Remarks 1 These material coefficients are only used to calculate the viscosity in viscous flows so for in visid flows this material card is not needed 2 As with other solvers in LS DYNA the user is responsible for unit consistency For example if dimensionless variables are used C and C2 should be replaced by the corresponding dimension less ones 2 32 CESE LS DYNA DEV February 19 2013 revision 2754 CESE_PART CESE CESE_PART Purpose Define CESE solver parts 1 e connect CESE material and EOS information Card 1 1 2 3 8 Variable PID MID EOSID Type I I I Default none none none VARIABLE DESCRIPTION PID Part identifier MID Material identifier defined by a CESE_MAT_ card EOSID Equation of state ident
86. ver iterations for the ball vertex linear sys tem RELTOL Relative tolerance to use as a stopping criterion for the iterative linear solv er 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 iterative linear solver LS DYNA DEV February 19 2013 revision 2754 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 Variable IFRAME IFLOW IGEOM Type I I I Default 0 0 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 giv en boundary conditions EQ 2 two dimensional 2D problem EQ 3 three dimensional 3D problem EQ 101 2D axisymmetric Remarks 1 Ifthe 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
87. y reer er ee her O 5 22 ICFD CONTROL TURBULENCE actos ia 5 23 KS PL DATABASE e e o ES 5 25 ICFD DATABASE AVERAGE scoreann eE REE EESAN OE EEREN Seelupaawiaectaslestancastieabotaees 5 26 PRD DEFINEN Doorriaaaana cts 5 27 ICFD INITIA Lerner i o re r E E da 5 28 A T ionia a EE ATE E OE NORE ENEAS EEN 5 29 ICED PAR Ti E EAE 5 31 E na a E E eer ee ere 5 32 ICFD SECTION soria A A RETEA AN 5 34 ACPO SEC Ne eden 5 35 MESA usado iaa 6 1 MESA BL diras 6 2 NESH BL SV oia ia 6 3 MESH EMBEDSHELL ci RA 6 4 ME NT ee il 6 5 MESHOSIZE cli 6 6 NESH SHE SHAPE attend 6 7 MESH SURFACE ELEMENT nirera r ea RE EAE N ETE RAV EREE E T 6 9 MESH SURPACE NODE roe E E O 6 10 MESH VOLUME cura aiii 6 11 MESA VOLUME ELEMENT aso 6 12 MESA VOLUME NODE csm 6 13 SIMI SEL VOLUME PART 0 a e 6 14 STOCHASTIC podios iio 7 1 STOCHASTICSP PRA Y SPAR TIOU ES ts isa 7 2 STOCHASTIC TES PARTICLES di 7 6 A O 8 1 LSO DOMAIN a 8 2 TESCO POINTE SE pi ii 8 4 LSO TIME SEQUENCE ona aii 8 5 LS DYNA DEV February 19 2013 revision 2754 0 5 TABLE OF CONTENTS 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 ar

Download Pdf Manuals

Related Search

Related Contents

www.WorkshopManuals.co.uk    Advisor - Benutzerhandbuch - OLAP Business Solutions GmbH  Communique de Presse  ADLER NEWS No.3  Hama 00113719 gaming control  キモチイイ SH。WER  New Machine Inspection and Installation PROCEDURE for  EL-200S DMX Manual  

Copyright © All rights reserved. Failed to retrieve file