ls-dyna_database

 LS-DYNA Database Binary Output Files Revised October 2012 Copyright ©, 1989-2012 LIVERMORE SOFTWARE TECHNOLOGY CORPORATION All Rights Reserved ®LS-DYNA DATABASE Support Address Livermore Software Technology Corporation 7374 Las Positas Road Livermore, California 94551 Tel: 925 449 2500 Fax: 925 449 2507 Email: support@lstc.com Copyright © 1989-2012 by Livermore Software Technology Corporation All rights Reserved LS-DYNA, LS-OPT and LS-PREPOST are registered trademarks of Livermore Software Technology Corporation. ii LS-DYNA DATABASE TABLE OF CONTENTS INTRODUCTION...........................................................................................................................1 STATE DATABASE (d3plot and d3part).....................................................................................2 CONTROL DATA...............................................................................................................3 MATERIAL TYPE DATA..................................................................................................7 FLUID MATERIAL ID DATA...........................................................................................8 SMOOTH PARTICLE HYDRODYNAMICS ELEMENT DATA FLAGS.......................9 PARTICLE DATA (NPEFG > 0)......................................................................................10 GEOMETRY DATA.........................................................................................................11 USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS............12 EXTRA 2 NODE CONNECTIVITY ARRAY FOR 10 NODE TETRAHEDRON ELEMENTS (ONLY IF NEL8 < 0)...................................................................................13 SMOOTH PARTICLE HYDRODYNAMICS NODE AND MATERIAL LIST..............13 RIGID ROAD SURFACE DATA.....................................................................................14 HEADER, PART & CONTACT INTERFACE TITLES..................................................15 DESCRIPTION OF BINARY FILE TYPES.....................................................................16 EXTRA DATA TYPES (OUTPUT FOR MULTI-SOLVER ANALYSIS)......................17 STATE DATA...................................................................................................................24 ELEMENT DELETION OPTION.....................................................................................31 SMOOTH PARTICLE HYDRODYNNAMICS NODE/ELEMENT STATE DATA......32 PARTICLE STATE DATA (NPEFG > 0).........................................................................33 ROAD SURFACE MOTION............................................................................................33 EXTRA DATA (MULTI-SOLVER ANALYSIS).............................................................34 END OF FILE MARKER..................................................................................................38 TIME HISTORY DATABASE (d3thdt)......................................................................................39 CONTROL DATA.............................................................................................................39 SMOOTH PARTICLE HYDRODYNAMICS ELEMENT DATA FLAGS.....................42 GEOMETRY DATA.........................................................................................................43 USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS............44 TIME HISTORY DATA....................................................................................................46 INTERFACE FORCE DATABASE............................................................................................56 CONTROL DATA.............................................................................................................56 GEOMETRY DATA.........................................................................................................58 USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS............59 STATE DATA...................................................................................................................61 FSIFOR FILE OUTPUT:...................................................................................................62 iii LS-DYNA DATABASE BLSTFOR FILE (NV2D=16 OR 7) OUTPUT:.................................................................62 CRACK FILE (d3crck)................................................................................................................64 DYNAIN BINARY FILE FORMAT (dynain.bin)......................................................................68 EXTRA DATA TYPE DEFINITIONS (NCFDV1 = 67108864).................................................69 4 LS-DYNA DATABASE INTRODUCTION The purpose of this information is to give guidance on how to access and read the various databases. The databases are written as word addressable fixed length binary files. The actual length depends on the amount of data saved, but will always be a multiple of 512 words (4 or 8 bytes each). Since it is likely that the database cannot be contained in a single file of length, FEMLEN, the data will spread over several files known as a family of files. Having a set of files enables them to be handled more easily than a single very large file. The root name for a family is the name of the first file member. Successive member names are compiled by appending a two or three digit number to the root name starting with 01, 02, through to 99, then 100 and ending with 999. Root names are limited to 75 characters. The original reason for a family of files was because the hard disks used for dyna3d runs could not cope with a single contiguous files large enough to contain all the data output. Subsequently, it has been found that splitting the output into separate files allows some unwanted data to be discarded and aids the copying, saving and movement of the data. Very large files can be impossible to transfer reliably over networks. Total output data can amount to several gigabytes or more depending on the model size. For ls-dyna runs with mesh adaptivity, the root name has a two letter appendage for each adapted mesh. Starting from ‘aa’ through ‘az’, then ‘ba’ through ‘bz’ and continuing up to ‘zz’, this gives a maximum of 676 possible adaptions. For example if the root name is ‘d3plot’ the subsequent files related to the original mesh are ‘d3plot01’, d3plot02, …, after adaption the new mesh and undeformed geometry is put in ‘d3plotaa’ and the subsequent files for the new mesh are: ‘d3plotaa01’, d3plotaa02, … The next adapted mesh will be in file ‘d3plotab’ and so on. A set of files at a particular adaption can be read separately by giving the root name with appendage, as the base file name. For example command: ‘lsprepost d3plot’ will read in all the file with root name ‘d3plot’ While: ‘lsprepost d3plotab’ will read in only files have ‘d3plotab’ in the name. LSPREPOST will read the binary databases separately or combined. Eg: lsprepost d3plot, lsprepost d3thdt, lsprepost iffname, lsprepost d3plot h=d3thdt f=iffname 1 Three databases are discussed, these are: 1. 2. 3. State Database (G=ptf, default name d3plot and d3part) Time History Database (F=thf, default name d3thdt) Interface Force Database (S=iff, no default name, typically: intfor) LS-DYNA DATABASE The file length used is set in the ls-dyna run as the default size of 7x512x512 words. The size can be changed on the command line with the ‘x=factor’ parameter giving a size of: factorx512x512 words. If the initial data or state data is larger than the given file length, the data will automatically split across files. This condition is not desirable because it is not clear whether any non root file can be discarded without destroying the continuity of the data. Ls-dyna checks before writing to a file, to ensure that there is room left in the file to contain the data at a particular state time. If not, it closes the current family member and starts writing the state data in the next file member. The files are written with a block size of 512 words, and if the data does not complete the last block it is padded out. This means that files cannot be concatenated and read together. The word size is 4 bytes for the single precision version of ls-dyna and 8 bytes for the double precision version, unless 32bit ieee format is defined, see *DATABASE_FORMAT, IBINARY. FILE GENERAL STRUCTURE The root file starts with a control words section, followed by node coordinates, then element connectivity for solids, thick shells, beams, and shells. Next are lists to reference the sequential internal numbering to the users number. State data is output next always starting with the time word. Data is of fixed length through the file members except where the mesh is adapted. The length of each area can be calculated from the information in the control words. The first file at adaption is like the root file in structure, so the new control words are used to recalculate the size of the subsequent data. The root file contains the initial data and also state data if there is room to write it. Further state data is written to the family members and each file will start with the time word provided data from the previous state did not overflow onto the file. If the disk address, DA, of the data being written exceeds the maximum file length, then data is written into file number int(DA/FAMLEN) at location DA-FAMLEN*int(DA/FAMLEN). If the state length is greater than the remaining length, the disk address is increased to start the writing at the beginning of the next file. STATE DATABASE (d3plot and d3part) There are three sections in this database. The first contains 64 words of control information plus extensions. The second contains geometric information including the nodal coordinates and element connectivities and user numbering lists. The third section contains the 2 LS-DYNA DATABASE results of the analysis at sequential output intervals. The output at a given time, called a state, contains a time word, global variables such as total energies and momenta for the whole model and each material (part), node data consisting of displacements, velocities, accelerations, and optionally temperatures, and finally element data that can include stresses and strains at integration points, and element deletion flags. The control data provides information about what is in the file and is used to calculate the various data length. There are two other state database files, namely: d3drfl and d3part, these are similar to d3plot but contain less data. The dynamic relaxation file, d3drfl, provides the state at the end of the DR process, while d3part is state output for a reduced number of parts in the model. CONTROL DATA VALUE Title Run time 10 1 DISK #WORDS ADDRESS 0 10 11 DESCRIPTION Model identification time in seconds since 00:00:00 UTC, January 1, 1970 d3plot=1 1=d3plot, 2=d3drlf, 3=d3thdt, 4=intfor, 5=d3part 6=blstfor, 7=d3cpm, 8=d3ale, 11=d3eigv, 12=d3mode, 13=d3iter, 21=d3ssd, 22=d3spcm, 23=d3psd, 24=d3rms, 25=d3ftg, 26=d3acs If > 1000, File type=INUM-1000 all external(users) numbers (Node, Element, Material and Rigid Surface Nodes) will be written in I8 format. Length of arbitrary numbering array = NARBS * 8 bytes for single precision files. Source version Version 1 1 12 13 14 ls-dyna version *1000000 + svn number Release number in character*4 form 50 for R5.0 511c for R5.1.1c Code version, floating number, eg 960.0 it is used to distinguish the floating point format, like cray, ieee, and dpieee NDIM 1 15 Number of dimensions (2 or 3). If 5 or 7 then an array of material types is read (MATTYP=1), element connectivities are unpacked and NDIM=3. If 4 then element connectivies are unpacked in the DYNA3D 3 Release number 1 INUM (File type) 1 LS-DYNA DATABASE database and NDIM is reset to 3. If >5 then state data contains movement of rigid road surface. NUMNP ICODE 1 1 16 17 Number of nodal points Flag to identify finite element code=2: old DYNA3D, code=6: NIKE3D, LS-DYNA/3D, LS-NIKE3D database NGLBV IT IU IV IA NEL8 1 1 1 1 1 1 18 19 20 21 22 23 Number of global variable to be read with each state NUMRW=number of rigid walls. NUMRBS=number of rigid body sets. = 6 + 6 * (NUMMAT8 + NUMMAT2 + NUMMAT4 + NUMATT+NUMRBS) + NUMRW * N N = 1 for DYNA3D and LS-DYNA3D N = 4 for LS-DYNA >= version 971 Flag for temperatures = 0, none, = 1, read in a temperature for each node = 2, read temperature for each node and heat flux for each node. = 3, read thermal shell middle temperature, thermal shell inner temperature, thermal shell outer temperature, and heat flux for each node. Solid node temperatures are repeated +=10, read mass scaling value for each node Flag for current geometry (=1 or 0) Flag for velocities (=1 or 0) Flag for accelerations (=1 or 0) Number of 8 node solid elements If NEL8 < 0, 2 extra nodes are output for ten node solids. Array is 2 * abs(NEL8), and is read after the arbitrary numbering arrays. NUMMAT8 BLANK BLANK NV3D NEL2 NUMMAT2 1 1 1 1 1 1 24 25 26 27 28 29 Number of materials used by the 8 node solids Insert zero Insert zero Number of values in database for each solid element. =7+NEIPH If NV3D is 8 * (7+NEIPH), each solid element has values at each Gauss point. Number of 2 node one-dimensional elements Number of materials used by the 2 node 1D elements 4 NV1D NEL4 NUMMAT4 NV2D NEIPH 1 1 1 1 1 30 31 32 33 34 LS-DYNA DATABASE Number of values in database for each 1D element = 6 + BEAMIP * 5 Number of four node two-dimensional elements Number of materials used by the 4 node 2D elements Number of values in database for each 2D element Are: MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS)+8 *IOSHL(3)+4*IOSHL(4)+12*ISTRN Number of additional values per solid element to be written in the type 6 database,NV3D=7+NEIPH, Actual number of history variables=NEIPH-6*ISTRN Number of additional values per integration point to be written into the type 6 database for shell elements. Number of integration points dumped for each shell. The magnitude of MAXINT must be greater than or equal to 3. if MAXINT>=0, then MDLOPT=0 MAXINT=MAXINT elseif MAXINT<0, then MDLOPT=1 MAXINT=abs(MAXINT) elseif MAXINT<10,000, then MDLOPT=2 MAXINT=abs(MAXINT)-10,000 endif MDLOPT controls the element deletion table (see below). The increase in state lengths allows deletion by nodes or elements. Element deletion flag (not standard) =xxx1 Solids deleted =xx1x Beams deleted =x1xx Shells deleted =1xxx Thick Shells deleted (Not used in LS-DYNA) Number of SPH Nodes Number of SPH materials Additional storage required for arbitrary node and element numbering in type 6 database =0 Sequential numbering =(10+NUMNP+ NEL8+NEL2+NEL4+NELT). 5 NEIPS MAXINT EDLOPT NMSPH NGPSPH NARBS 1 1 35 36 1 1 1 1 37 37 38 39 LS-DYNA DATABASE NELT NUMMATT NV3DT IOSHL(1) IOSHL(2) IOSHL(3) IOSHL(4) IALEMAT NCFDV1 NCFDV2 NADAPT NMMAT NUMFLUID 1 1 1 1 1 1 1 1 1 1 1 1 1 40 41 42 43 44 45 46 47 48 49 50 51 52 Number of 8 node thick shell elements. MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS)+ 12*ISTRN Number of materials used for the 8 node thick shell element. Number of values in database for each thick shell 6 stress components flag, if 1000 =1 else =0 Plastic strain flag, if 1000 =1 else =0 Shell force resultants flag, if 1000 =1 else =0 Shell thickness, energy+2 others, if 1000 =1 else =0 Size of array containing solid element parts numbers used as ALE material Bit flags for CFD nodal values. If = 67108864, then state contains CFD extra data – see below (ls980 version) Further bit flags for CFD nodal values. If extra data, then value equals number of data domains (ls980 version) Number of adapted element to parent pairs (not implemented) Total number of materials – not set in LS-DYNA/3D Total number of ALE fluid groups. Fluid density and volume fractions output as history variables, and a flag for the dominant group. If negative multi-material species mass for each group is also output. Order is: rho, vf1, … vfn, dvf flag, m1, … mn. Density is at position 8 after the location for plastic strain. Any element material history variables are written before the Ale variables, and the six element strains components after these if ISTRN=1. Invariant node numbering fore shell and solid elements See INN in card *CONTROL_ACCURACY Number of particle method data sets. Not used. Rate of change of temperature per node flag. INN NPEFG NVEFG IDTDT 1 1 1 1 53 54 55 56 6 NWORD WORDS 1 6 57 58-63 LS-DYNA DATABASE An array of dT/dt values of length NUMNP. Array is written after node temperature arrays. Additional number of control words. Used by D3THDT and INTFOR The value of ISTRN must be computed, it is not output in the control data ISTRN can only be computed as follows and if NV2D > 0. If NV2D-MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS)+8*IOSHL(3)+4*IOSHL(4) > 1 Then ISTRN = 1, else ISTRN = 0 If ISTRN=1, and NEIPH>=6, last the 6 additional values are the six strain components. Or NELT > 0 If NV3DT-MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS) > 1 Then ISTRN = 1, else ISTRN = 0 MATERIAL TYPE DATA The material section contains the material type numbers. This section is skipped if MATTYP is zero. This data is required because those shell elements that are in a rigid body have no element data output in the state data section. The normal length of the shell element state data is: NEL4 * NV2D, when the MATTYP flag is set the length is: (NEL4 – NUMRBE) * NV2D. When reading the shell element data, the material number must be check against IRBRTYP list to find the element’s material type. If the type = 20, then all the values for the element to zero. This option is set in *DATABASE_EXTENT_BINARY, with DCOMP=2 VALUE NUMRBE NUMMAT IRBTYP 1 1 NUMMAT Number of rigid body shell elements. Number of materials in the database. Material type numbers LENGTH DESCRIPTION 7 LS-DYNA DATABASE FLUID MATERIAL ID DATA The fluid material section contains the material numbers for solid elements that are used to define an Euler grid or Arbitrary Lagrangian Euler mesh. This section is skipped if IALEMAT is zero. VALUE FLUIDID LENGTH IALEMAT DESCRIPTION Fluid material number used in solid element mesh 8 LS-DYNA DATABASE SMOOTH PARTICLE HYDRODYNAMICS ELEMENT DATA FLAGS This section is only output if NMSPH > 0. The section is a list of flags to indicate what SPH data is output for each SPH node/element. The first number is the length in words for this array, currently = 11. SPH elements are centered at nodes, and cover a spherical volume defined by the radius of influence. They do not have a connection with other SPH elements. They should be displayed as a dot or a spherical surface, with radius scaling to reduce the size and enable each element to be distinguishable. As follows: If any value of isphfg(2) through isphfg(11) = 0, then the particular data item is not output for the particle. To calculated the size of data add the isphfg values from isphfg(2) through isphfg(11) plus one. One value is always output which is the material number as a floating point number for each particle. If this value is negative then the particle has been deleted from the model. Note: it is possible a SPH element could be deleted, or be non active in the initial states, and become active in later states. Full output for each particle is: mat#, radius, pressure, {sx, sy, sz, sxy, syz, sxz} ps, rho, ie, nn, {ex, ey, ez, exy, eyz, exz}, mass, hv1 … hvn. NUM_SPH_VARS = 1 + sum of isphfg(i), i=2 to isphfg(1) Hence, total size is 20 + the total number of history variables. When a particle is deleted from the model, data is still output for it because the length of data must always be the same for each state. 9 isphfg(1) = 11 - length of sph flags array isphfg(2) = 1 - radius of influence isphfg(3) = 1 - pressure in particle isphfg(4) = 6 - 6 true stress components isphfg(5) = 1 - plastic strain, > 0.0 if effective stress exceeds yield strength isphfg(6) = 1 - density of particle material isphfg(7) = 1 - internal energy (strain) isphfg(8) = 1 - number of neighbors affecting particle isphfg(9) = 6 - 6 true strain components isphfg(10)=1 - mass of element (>= ls971) isphfg(11)=1 - max number of sph history variables. LS-DYNA DATABASE PARTICLE DATA (NPEFG > 0) Control block If NPEFG > 0 airbag particles are output The first three digits of NPEFG are the number of airbags in the database = NPARTGAS NPARTGAS = NPRFG % 1000 SUBVER = NPEFG / 1000 In the extended control block: The first four words in the block are: 1. 2. 3. 4. 5. NLIST = NGEOM + NVAR + NSTGEOM NLIST words of output for variables listed to define the type of each variable, =1 for integer and 2= for floating point 2 * NLIST words of variable names (8 bytes per name or 16 for double precision output). NGEOM NVAR NPART number of geometry variables number of state variables number of particles number of state geometry variables number of chambers NSTGEOM NCHAMBER If SUBVER == 4 10 LS-DYNA DATABASE GEOMETRY DATA The geometry section contains the nodal coordinates and the element connectivities. The ordering of the nodal points is the same as the ordering of the nodal data in the state data that follows. If NDIM=3 the connectivities are assumed to be packed with 3 integers per word, if NDIM>3, then connectivities are not pack, (the default for LS-DYNA, LS-DYNA3D and LS-NIKE3D. The order of the elements are 3, 2, and 1 dimensional elements if the database is ICODE=2 or 6. VALUE X(3,1) IX8(9,1) If NEL8 < 0 IXT(9,1) IX2(6,1) LENGTH DESCRIPTION NDIM*NUMNP Array of nodal coordinates X1,Y1,Z1, X2,Y2,Z2, X3,Y3,Z3, ... ,Xn,Yn,Zn 9*NEL8 2*abs(NEL8) 9*NELT 6*NEL2 Connectivity and material number for each 8 node solid element. Extra nodes for ten node solids. Connectivity and material number for each 8 node thick shell element. Connectivity, orientation node, two null entries, and the material number for each 2 node beam element. For some beam types the last two number contain the beam type and length to width ratio * 100 and length to height ratio * 100 type = ix2(5,*) & 0x3F width = 0.01 * length / (ix2(5,*)>>6 height = 0.01 * length / ix2(6,*) Third node (orientation) may be > 1e9 Contain flag 1e9 to indicate a spot weld. IX4(5,1) 5*NEL4 Connectivity and material number for each 4 node shell element Note the node numbers are the LS-DYNA internal numbers for nodes, these will be the same as the user’s numbers if NARBS = 0, otherwise, the arbitrary number lists are used to find the user’s numbers, similarly, for element numbers and material numbers. 11 LS-DYNA DATABASE USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS Skip this section if NARBS (disk address 39) is zero. The user node and element numbers must be in ascending order. It assumed that if this option is used all node and element data anywhere in the databases is in ascending order in relation to the user numbering. The total length of the data in this data is equal to: NARBS=10+NUMNP+NEL8+NEL2+NEL4+NELT, if sequential numbering is used for the materials/parts. For arbitrary material numbering (NSORT < 0), the total length is increased by 6+NUMMAT8+NUMMAT4+NUMMAT2+NUMMATT. Material numbers are not in ascending order. VALUE NSORT LENGTH 1 DESCRIPTION Pointer to arbitrary node numbers in LS-DYNA source code. If < 0, it flags that arbitrary material identification numbers are also used. Pointer to arbitrary solid element numbers in LS-DYNA source code: =NSORT+NUMNP Pointer to arbitrary beam element numbers in LS-DYNA source code: =NSRH+NEL8 Pointer to arbitrary shell element numbers in LS-DYNA source code: =NSRB+NEL2 Pointer to arbitrary thick shell element numbers in LS-DYNA source code: =NSRS+NEL4 Number of nodal points Number of 8 node solid elements Number of 2 node beam elements Number of 4 node shell elements Number of 8 node thick shell elements Pointer to an array in the LS-DYNA source code that list the material ID’s in ascending order. Pointer to an array in the LS-DYNA source code that gives the material ID’s in the actual order that they are defined in the user input. 12 NSRH NSRB NSRS NSRT NSORTD NSRHD NSRBD NSRSD NSRTD NSRMA NSRMU 1 1 1 1 1 1 1 1 1 1 1 VALUE NSRMP LENGTH 1 DESCRIPTION LS-DYNA DATABASE Pointer to an array in the LS-DYNA source code that gives the location of a member in the array originating at NSRMU for each member in the array starting at NSRMA. Total number of materials Total number of nodal rigid body constraint sets Total number of materials Array of user defined node numbers Array of user defined solid element numbers Array of user defined beam element numbers Array of user defined shell element numbers Array of user defined thick shell numbers Ordered array of user defined material ID’s Unordered array of user material ID’s Cross reference array NSRTM NUMRBS NMMAT NUSERN NUSERH NUSERB NUSERS NUSERT NORDER NSRMU NSRMP 1 1 1 NSORTD NSORTH NSORTB NSORTS NSORTT NMMAT NMMAT NMMAT EXTRA 2 NODE CONNECTIVITY ARRAY FOR 10 NODE TETRAHEDRON ELEMENTS (only if NEL8 < 0) List of extra nodes for each 10 node tetrahedron element, 2 * abs(NEL8). Any 8 node solids have these two nodes set to zero. ADAPTED ELEMENT PARENT LIST (not implemented) List of element id pairs for H-type shell element adaptivity. Length of data is 2 * NADAPT, pairs are element number and element parent number SMOOTH PARTICLE HYDRODYNAMICS NODE AND MATERIAL LIST If NMSPH > 0 Length of data PARTICLE GEOMETRY DATA (NPEFG > 0) 13 List of sph node and its material number 2 * NUMSPH LS-DYNA DATABASE NPARTGAS blocks of NGEOM data to describe the geometry for each airbag: 1. 2. 3. 4. 5. RIGID ROAD SURFACE DATA If NDIM > 5 NNODE NSEG NSURF MOTION NODEID Number of nodes in road surface Total number of 4 noded road surface segments Number of road surfaces Flag to indicate motion data is output for each state NNODE list of IDs first particle ID for the airbag number of particles in the airbag ID for the airbag number of gas mixtures in the airbag number of chambers If NGEOM == 5 SURFNODE XYZ Coordinate for each node Lists of 4 noded segments for each surface SURFID Surface ID Number SURFNSEG Number of segments in surface SURFSEGS SURFNSEG of 4 node ids for each segment Length of data = 4 + NNODE + 3 * NNODE + NSURF * (2 + 4 * SURFNSEG) 14 LS-DYNA DATABASE HEADER, PART & CONTACT INTERFACE TITLES At the end of the first binary files, eg d3plot, the part and model titles are appended. If the model input includes *DATABASE_BINARY_D3PROP, all the d3prop part data is included. At the end of the first interface force file, titles and contact id are appended. This extra data is written at the end of the following files: d3plot, d3part and intfor files, and the header and part titles are written directly after the EOF (= -999999.0) marker. Header output ------------------------------------ NTYPE 1 entity type = 90000 HEAD 18 Header title (72 characters) For the interface force file (intfor), header and contact titles are written at the end of first file after the EOF (= -999999.0) marker Part title output Value Length Description ------------------------------- NTYPE 1 entity type = 90001 NUMPROP 1 number of parts For NUMPROP parts: IDP 1 part id PTITLE 18 Part title (72 characters) For the interface force file (intfor), header and contact titles are written at the end of first file after the EOF (= -999999.0) marker. Contact title output ------------------------------------ NTYPE 1 entity type = 90002 NUMCON 1 number of contacts For NUMCON contacts: IDC 1 contact id CTITLE 18 Contact title (72 characters) Header output ------------------------------------ NTYPE 1 entity type = 90000 HEAD 18 Header title (72 characters) The d3prop data is written to the d3plot file only if it is requested. 15 LS-DYNA DATABASE D3PROP output Values Length Description ------------------------------- NTYPE 1 entity type = 900100 NLINE 1 number of keyword lines For NLINE keyword lines: KEYWORD 20 keyword line (80 characters) DESCRIPTION OF BINARY FILE TYPES Control word 11 File type: 1=d3plot plot file of model and state data 2=d3drlf plot file of model and state data from a dynamic relaxation analysis 3=d3thdt time history plot file for a set of nodes and elements 4=intfor plot file of contact interfaces 5=d3part plot file of model and state data for a set of parts 6=blstfor plot file for a blast wave analysis 7=d3cpm 8=d3ale plot file for ale fluid-structure interface or fsifor 11=d3eigv plot file for an eigen value analysis 12=d3mode 13=d3iter 21=d3ssd plot file for steady state dynamic response. 22=d3spcm plot file for response spectrum analysis. 23=d3psd plot file for power spectral density of response, in random vibration. 24=d3rms plot file for root mean square of response, in random vibration. 25=d3ftg plot file for random fatigue analysis. 26=d3acs plot file for frequency domain acoustic FEM analysis 16 LS-DYNA DATABASE EXTRA DATA TYPES (Output for Multi-Solver Analysis) If NCFDV1 = 67108864, then NCFDV2 will be the number of additional datasets from different solver-mesh combinations that are included in the d3plot file. One of each of the solver-mesh combinations listed below can be among the NCFDV2 datasets. Currently defined solver-mesh combinations follow. For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. solver and domain ID: PFEM_IF number of volume vars output: nvolvar_pfem first volume variable ID: ID 1 ... last volume variable ID: ID nvolvar_pfem number of PFEM parts nPFEM_parts first internal part ID: partID 1 ... last internal part ID: partID nPFEM_parts first user part ID: user_partID 1 ... last user part ID: user_partID nPFEM_parts solver and domain ID: PFEM_IF_SURFACE number of surface vars output: nsurfvar_pfem first surface variable ID: ID 1 ... last surface variable ID: ID nsurfvar_pfem number of PFEM parts nPFEM_surfparts first internal part ID: partID 1 ... last internal part ID: partID nPFEM_surfparts first user part ID: user_partID 1 ... last user part ID: user_partID nPFEM_surfparts For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. solver and domain ID: CESE number of volume vars output: nvolvar_cese first volume variable ID: ID 1 ... last volume variable ID: ID nvolvar_cese number of CESE parts nCESE_parts first internal part ID: partID 1 ... last internal part ID: partID nCESE_parts For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. 17 LS-DYNA DATABASE first user part ID: user_partID 1 ... last user part ID: user_partID nCESE_parts solver and domain ID: CESE_SURFACE number of surface vars output: nsurfvar_cese first surface variable ID: ID 1 ... last surface variable ID: ID nsurfvar_cese number of CESE parts nCESE_surfparts first internal part ID: partID 1 ... last internal part ID: partID nCESE_surfparts first user part ID: user_partID 1 ... last user part ID: user_partID nCESE_surfparts For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. solver and domain ID: EM number of volume vars output: nvolvar_EM first volume variable ID: ID 1 ... last volume variable ID: ID nvolvar_EM number of EM parts nEM_parts first internal part ID: partID 1 ... last internal part ID: partID nEM_parts first user part ID: user_partID 1 ... last user part ID: user_partID nEM_parts For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. solver and domain ID: EM_SURFACE number of surface vars output: nsurfvar_EM first surface variable ID: ID 1 ... last surface variable ID: ID nsurfvar_EM number of EM parts nEM_surfparts first internal part ID: partID 1 ... last internal part ID: partID nEM_surfparts first user part ID: user_partID 1 ... last user part ID: user_partID nEM_surfparts For the following domain, the mesh can be completely different for each output state, so no mesh is output in this control block. 18 LS-DYNA DATABASE solver and domain ID: CESE_CFD_NODE size of each variable component: numnp_cese number of nodes: numnp_cese number of elements: numelh_cese user node numbers: nodes_cese_cfd(numnp_cese) array of nodal coordinates: x_cese_cfd(3, numnp_cese) element connectivity: ix8_cese_cfd(9, numelh_cese) number of output vars: nv_cese_cfd_node first variable ID: ID 1 ... last variable ID: ID nv_cese_cfd_node number of CESE parts ncese_parts first internal part ID: partID 1 ... last internal part ID: partID ncese_parts first user part ID: user_partID 1 ... last user part ID: user_partID ncese_parts user element number: for the first CESE element ... user element number: for the last CESE element In this domain, the variables are defined at the element centroid. solver and domain ID: CESE_CFD_ELEMENT size of each variable component: numelh_cese number of nodes: numnp_cese number of elements: numelh_cese user node numbers: nodes_cese_cfd(numnp_cese) array of nodal coordinates: x_cese_cfd(3, numnp_cese) element connectivity: ix8_ins_cfd(9, numelh_cese) number of output vars: nv_cese_cfd_ele first variable ID: ID 1 ... last variable ID: ID nv_cese_cfd_ele number of CESE parts ncese_parts first internal part ID: partID 1 ... last internal part ID: partID ncese_parts first user part ID: user_partID 1 ... last user part ID: user_partID ncese_parts user element number: for the first CESE element ... user element number: for the last CESE element In this domain, the variables are defined at the mesh nodes. solver and domain ID: CESE_CFD_ELEMENT_TS size of each variable component: 4*numelh_cese number of nodes: numnp_cese number of elements: numelh_cese user node numbers: nodes_cese_cfd(numnp_cese) array of nodal coordinates: x_cese_cfd(3, numnp_cese) element connectivity: ix8_ins_cfd(9, numelh_cese) In this domain, the variables are defined by Taylor series expanded around the element centroid. 19 LS-DYNA DATABASE number of output vars: nv_cese_cfd_ts first variable ID: ID 1 ... last variable ID: ID nv_cese_cfd_ts number of CESE parts ncese_parts first internal part ID: partID 1 ... last internal part ID: partID ncese_parts first user part ID: user_partID 1 ... last user part ID: user_partID ncese_parts user element number: for the first CESE element ... user element number: for the last CESE element In this domain, the variables are defined on structural solid elements. solver and domain ID: EM_FEMSTER_SOLID_INTEG_PTS size of each variable component: nip_solid_em * numelh number of nodes: numnp number of elements: numelh number of integration points: nip_solid_em number of output vars: nv_em_solid_integ first variable ID: ID 1 ... last variable ID: ID nv_em_solid_integ In this domain, the variables are defined on structural thick shell elements. solver and domain ID: EM_FEMSTER_TSHELL_INTEG_PTS size of each variable component: nip_tshell_em * numelt number of nodes: numnp number of elements: numelt number of integration points: nip_tshell_em number of output vars: nv_em_tshell_integ first variable ID: ID 1 ... last variable ID: ID nv_em_tshell_integ In this domain, the variables are defined on structural thin shell elements. solver and domain ID: EM_FEMSTER_SHELL_INTEG_PTS size of each variable component: nip_shell_em * numels number of nodes: numnp number of elements: numels number of integration points: nip_shell_em number of output vars: nv_em_shell_integ first variable ID: ID 1 ... last variable ID: ID nv_em_shell_integ In this domain, the variables are defined at the centroids of structural solid elements. solver and domain ID: EM_FEMSTER_SOLID_CENTROID size of each variable component: numelh number of nodes: numnp 20 LS-DYNA DATABASE number of elements: numelh number of output vars: nv_em_solid_cent first variable ID: ID 1 ... last variable ID: ID nv_em_solid_cent In this domain, the variables are defined at the centroids of structural thick shell elements. solver and domain ID: EM_FEMSTER_TSHELL_CENTROID size of each variable component: numelt number of nodes: numnp number of elements: numelt number of output vars: nv_em_tshell_cent first variable ID: ID 1 ... last variable ID: ID nv_em_tshell_cent solver and domain ID: EM_FEMSTER_SHELL_CENTROID size of each variable component: numels number of nodes: numnp number of elements: numels number of output vars: nv_em_shell_cent first variable ID: ID 1 ... last variable ID: ID nv_em_shell_cent In this domain, the variables are defined at the mesh nodes. solver and domain ID: EM_FEMSTER_AIR size of each variable component: nip_air_em * numelh_air_em number of nodes: numnp_air_em number of elements: numelh_air_em number of integration points: nip_air_em user node numbers: nodes_air_em(numnp_air_em) array of nodal coordinates: x_air_em(3,numnp_air_em) element connectivity: ix8_air_em(8, numelh_air_em) number of output vars: nv_em_air_integ first variable ID: ID 1 ... last variable ID: ID nv_em_air_integ In this domain, the variables are defined at the centroids of structural thin shell elements. solver and domain ID: RECT_AIR_EM_NODE size of each variable component: nx_rect_air_em * ny_rect_air_em * nz_rect_air_em number of x nodes: nx_rect_air_em number of y nodes: ny_rect_air_em number of z nodes: nz_rect_air_em minimum x coordinate: xmin_rect_air_em minimum y coordinate: ymin_rect_air_em minimum z coordinate: zmin_rect_air_em maximum x coordinate: xmax_rect_air_em maximum y coordinate: ymax_rect_air_em maximum z coordinate: zmax_rect_air_em number of output vars: nv_em_air_nd first variable ID: ID 1 ... In this domain, the variables are defined at the nodes of the implied rectangular mesh. 21 LS-DYNA DATABASE last variable ID: ID nv_em_air_nd In this domain, the variables are defined on faces of structural elements. solver and domain ID: EM_FEMSTER_BEM size of each variable component: nip_bem_em * nfaces_bem_em number of nodes: numnp_bem_em number of elements: nfaces_bem_em number of integration points: nip_bem_em number of BEM parts: em_numPartBem flag for BEM mesh: nBEMflag (first bit =0 if no motion,=1 if motion) (second bit =0 if no edge domain,=1 if edge domain) (third bit =0 if no node domain,=1 if node domain) user node numbers: nodes_bem_em(numnp_bem_em) array of nodal coordinates: x_bem_em(3, numnp_bem_em) element connectivity: ix4_bem_em(5, nfaces_bem_em) if (second bit(nBEMflag) = 1) number of edges nedges_bem_em number of edge domains nedgedomain_bem_em node edge connectivity edgex2_bem_em(3, nedges_bem_em) (internal node1,internal node2,partId) number of edges per domain numEdgesPerDomain(nedgedomain_bem_em) edge domain list edgeDomainList(sum(numEdgePerDomain(i))) endif if (third bit(nBEMflag) = 1) node element connectivity nodex4_bem_em(5, nfaces_bem_em) number of node domains nnodedomain_bem_em node domain array nodeDomain(numnp_bem_em) endif number of output vars: nv_em_bem_integ first variable ID: ID 1 ... last variable ID: ID nv_em_bem_integ number of output vars: n_prtcl_vars first variable ID: ID 1 ... last variable ID: ID n_prtcl_vars In this domain, the variables are defined at the particle positions. solver and domain ID: STOCHASTIC_PARTICLES Notes: The variable IDs are grouped into three groups: 1) D3PL_FIRST_SCALAR_ID <= ID < D3PL_FIRST_VECTOR_ID are scalar variables 2) D3PL_FIRST_VECTOR_ID <= ID < D3PL_FIRST_TENSOR_ID are vector variables (3 components per entry) 3) D3PL_FIRST_TENSOR_ID <= ID < D3PL_END_IDS 22 LS-DYNA DATABASE are symmetric tensor variables (6 component per entry) When a number of integration points are specified, it is assumed that they are distributed at the Gauss points of the given element type based upon how many are output. That is, for shell or face elements, 4 output points would imply the 2x2 Gauss points are used, while 9 output points would imply the 3x3 Gauss points are used, and so forth. Similarly, for volume elements, 8 output points would imply the 2x2x2 Gauss points are used, while 27 output points would imply the 3x3x3 Gauss points are used, and so forth. 23 LS-DYNA DATABASE STATE DATA The state data has three parts: VALUE TIME GLOBAL LENGTH 1 NGLBV DESCRIPTION Time word Global variables for this state. LS-DYNA Global Variables: KE, IE, TE, X, Y, and Z velocity IE for each material KE for each material X, Y, and Z velocity for mat 1 ... X, Y, and Z velocity for mat n Mass for each material Force for each rigid wall Xyz position of wall (ls971) = 6 + 7 * (NUMMAT8 + NUMMAT2 + NUMMAT4 + NUMMATT+NUMRBS) + N * NUMRW, N=1, for ls-dyna(ls971) N=4 Total nodal values for state. FOR LS-DYNA3D and LS-DYNA IT=1, node temperatures only, N=0 IT=2, node temperature and node flux, N=2 IT=3, 3 temperature per node and node flux, N=3 Temperature for shell node at inner, middle and outer layer, inner array, middle array, outer array. IT/10=1, mass scaling value at node. N+=1 =((IT+N)+NDIM*(IU+IV+IA))*NUMNP where IT=temperature flag, IU=cooordinates flag, IV=velocities flag, and IA=accelerations flag. Bit flag: NCFDV1, bits from right to left eg, Pressure, Resultant Vorticity, and Density NCFDV1=2+32+1024=1058 2 Pressure 3 X Vorticity 4 Y Vorticity 5 Z Vorticity 6 Resultant Vorticity 7 Enstrophy 8 Helicity 9 Stream Function 24 • • • Time word and global data Node data Element data for solids, shell, and beams, respectively NODEDATA CFDDATA NND CFD LS-DYNA DATABASE 10 Enthalpy 11 Density 12 Turbulent KE 13 Dissipation 14-20 Eddy Viscosity Bit flag: NCFDV2 2-11 Species 1 through 10 ELEMDATA ENN Total element data for state. =NEL8*NV3D+NELT*NV3DT+NEL2*NV1D+ NEL4*NV2D+NMSPH*NUM_SPH_VARS The organization of the element data for each element type is described below. The data for the solid elements (7 values/element) is printed first, followed by the data for the beam elements (6 values/element), and then the data for the shell elements (typical 33 or 45 values/element depending on whether the strains are included). This state data is repeated for each state in the database. Element data is defined at the integration points within the element. Contour and fringe plots require that the data be extrapolated to the nodal points. In LS_PREPOST the element values are averaged at the nodes. Element strains are not output by default, these are only output for solids, shell, and thick shell when *DATABASE_EXTENT_BINARY, STRFLG=1 25 LS-DYNA DATABASE SOLID ELEMENTS – 8 node Hexahedron, other solid elements like wedge, pyramid, and tetrahedron are identified by repeated final connectivities. Eg tetrahedron = 1,2,3,4,4,4,4,4 The database for solid elements consists of 7+NEIPH values per element. NEIPH extra values are defined if and only if NEIPH is greater than zero or if the model is an ALE analysis. If strain components are output, then the last 6 neiph values are true strains: ex, ey, ez, exy, eyz, exz, in the global system. They are: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10 . Sigma-x (true stress in the global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable First extra value (if NEIPH>0) Second extra value (if NEIPH >1) Etc. until NEIPH extra values are defined if ISTRN=1 7+NEIPH-5. Epsilon-x 7+NEIPH-4. Epsilon-y 7+NEIPH-3. Epsilon-z 7+NEIPH-2. Epsilon-xy 7+NEIPH-1. Epsilon-yz 7+NEIPH. Epsilon-zx 26 LS-DYNA DATABASE For thick shell elements the database contains NV3DT = MAXINT * (6 * IOSHL(1) + IOSHL(2) + NEIPS) +12 * ISTRN values per element. Three sets of global stresses are always put into the database for each thick shell and are located at the mid surface, the inner integration point surface, and the outer integration point surface, respectively. If one integration point is used the same through the thickness stress state is outputted three times. If two integration points are used then the mid surface value is taken as the average value. The inner values of the stress are always set to the values at the innermost integration point and likewise for outer values. If the integration point does not lie at the center, ie, an even number of integration points through the thickness, a value is computed that is an average of the two integration point nearest the mid surface. The IOSHL flags indicate which shell element data is included which is suppressed. The flags are set in ls-dyna by *DATABASE_EXTENT_BINARY, SIGFLG, EPSFLG, RLFLG, and ENGFLG The ordering of the data follows: 1. 2. 3. 4. 5. 6. 7. *. 8. 9. 10. 11. 12. 13. *. 15. 16. 17. 18. 19. 20. Sigma-x (mid surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for mid surface Sigma-x (inner surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Define NEIPS additional history values here for inner surface Sigma-x (outer surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx 27 14. Effective plastic strain or material dependent variable LS-DYNA DATABASE 21. *. *. If MAXINT > 3 then define an additional (MAXINT-3 )* (6 * IOSHL(1) +1*IOSHL(2)+NEIPS) quantities here. For beam elements the database contains NV1D=6 values per element. They are: 1. 2. 3. 4. 5. 6. 1. 2. 3. 4. 5. Axial force S shear resultant T shear resultant S bending moment T bending moment Torsional resultant RS shear stress TR shear stress Axial stress Plastic strain Axial strain Effective plastic strain or material dependent variable Define NEIPS additional history values here for outer surface If ISTRN=1, then define strain components Epsilon (x, y, z, xy, yz, zx) here for inner surface and outer surface If there are values output at beam integration points, then NV1D = 6 + 5 * BEAMIP BEAMIP is set in *DATABASE_EXTENT_BINARY 28 For shell elements the database contains NV2D values, where: LS-DYNA DATABASE NV2D=MAXINT* (6*IOSHL(1) + 1*IOSHL(2) + NEIPS) +8*IOSHL(3) + 4*IOSHL(4) + 12*ISTRN values per deformable element. If MATTYP=1 and IRBTYP(I)=20, where I=internal element number, then the material is rigid and the compressed database contains no data for the element. If the minimum value of MAXINT is used, i.e., 3, the stresses are typically located at the mid surface, the inner surface, and the outer surface, respectively. If one integration point is used the stress is written three times. If two integration points are used then the mid surface value is taken as the average value. The inner values of the stress are always set to the values at the innermost integration point and likewise for outer values. If no integration point lies at the center, i.e., an even number of integration points through the thickness, a value is computed that is an average of the two integration point lying nearest the mid surface. The ordering of the data follows: 1. 2. 3. 4. 5. 6. 7. *. 8. 9. 10. 11. 12. 13. *. Sigma-x (mid surface true stress in global system) Sigma-y Sigma-z Sigma-x Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for mid surface Sigma-x (inner surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Define NEIPS additional history values here for inner surface 14. Effective plastic strain or material dependent variable 29 LS-DYNA DATABASE 15. 16. 17. 18. 19. 20. 21. Sigma-x (outer surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable *. Define NEIPS additional history values here for outer surface If MAXINT >3 then define an additional (MAXINT-3 )* (6*IOSHL(1) + 1*IOSHL(2) + 8*IOSHL(3) + 4*IOSHL(4) + NEIPS) quantities here 22. Bending moment-Mx (local shell coordinate system) 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. Bending moment-My Bending moment-Mxy Shear resultant-Qx Shear resultant-Qy Normal resultant-Nx Normal resultant-Ny Normal resultant-Nxy Thickness Element dependent variable Element dependent variable 33. Internal energy (if and only if ISTRN=0) The following quantities are expected if and only if ISTRN=1 33. eps-x (inner surface strain in global system) 34. 35. 36. 37. 38. 39. 40. 41. 42. 43. 44. 45. eps-y eps-z eps-xy eps-yz eps-zx eps-x (outer surface strain in global system) eps-y eps-z eps-xy eps-yz eps-zx Internal energy (if and only if NV2D>=45) 30 LS-DYNA DATABASE ELEMENT DELETION OPTION Skip this section if the word MAXINT is greater than or equal to zero, (MDLOPT>=0). If MDLOPT=1, then the list is equal to the number of nodal points (NUMNP) and contains a one if the node is visible and a zero if the node is not visible, (only used in vec-dyna3d). If MDLOPT=2, then the list equals the total number of elements (NEL8 + NELT + NEL4 + NEL2), in this order, and each value is set to the element material number or =0, if the element is deleted. All these numbers are output as floating point values and not integers. 31 LS-DYNA DATABASE SMOOTH PARTICLE HYDRODYNNAMICS NODE/ELEMENT STATE DATA This section is only output if NMSPH>0 For each SPH node the follow values are output: NUM_SPH_DATA = 1 + ∑ isphfg(i), i=2:10 Length of data = NUM_SPH_DATA * NUMSPH Material number, if <=0 then element is deleted. Currently isphfg(1) = 10, ie number of sph data flags, this could be changed in the future. If isphfg(2) =1, radius of particle influence If isphfg(3) =1, pressure in particle If isphfg(4) =6, stress components for particle, sx, sy, sz, sxy, syz, sxz If isphfg(5) =1, plastic strain for particle If isphfg(6) =1, density of particle material If isphfg(7) =1, internal energy of particle If isphfg(8) =1, number of particle neighbors If isphfg(9) =6, strain components for particle, ex, ey, ez, exy, eyz, exz If isphfg(10)=1, mass of element (ls971) Note: it is possible a SPH element could be deleted, or be none active in the initial states, and become active in later states. 32 PARTICLE STATE DATA (NPEFG > 0) STATE DATA LS-DYNA DATABASE NPARTGAS blocks of NSTGEOM data to describe the state geometry for each bag: 1. 2. PARTICLE DATA NVAR words of data output for each particle: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. ROAD SURFACE MOTION If NDIM > 5 output rigid body displacement, dx, dy, dz and velocity, vx, vy, vz of each road surface. Length of data = 6 * NSURF 33 gas ID chamber ID leakage flag, 0 active, -1 fabric, -2 vent hole, -3 mistracked mass radius spin energy translational energy distance from particle to nearest segment x position y position z position x velocity y velocity z velocity number of active particles current bag volume LS-DYNA DATABASE EXTRA DATA (Multi-Solver Analysis) If NCFDV1 = 67108864, then the state data includes NCFDV2 additional datasets from solver-mesh combinations specified after the \"User material, node, and element identification numbers\" for the structural mesh. State data of the first solver-mesh combination ... State data of the last (NCFDV2-th) solver-mesh combination When the state data comes from the PFEM_IF domain, then the mesh is output first, followed by the data. Currently, the mesh is entirely tetrahedral, but we anticipate users will also specify mixed meshes in the near future: size of each volume variable component: nnpvol_pfem number of volume nodes: nnpvol_pfem number of tetrahedral elements: ntet_pfem number of pyramid elements: npyr_pfem number of wedge elements: nwdg_pfem number of hexahedral elements: nhex_pfem user volume node numbers: volnodes_pfem(nnpvol_pfem) array of volume nodal coordinates: xvol_pfem(3, nnpvol_pfem) tetrahedral element connectivity: ix4_pfem(5, ntet_pfem) pyramid element connectivity: ix5_pfem(6, npyr_pfem) wedge element connectivity: ix6_pfem(7, nwdg_pfem) hexahedral element connectivity: ix8_pfem(9, nhex_pfem) data for 1st volume variable (size is nnpvol_pfem) ... data for nvolvar_pfem-th volume variable (size is nnpvol_pfem) Notes: (1) the first four entries of ix4_pfem(5, ntet_pfem) for each element are the volume node numbers, while the last entry is the volume part number. (2) the first five entries of ix5_pfem(6, npyr_pfem) for each element are the volume node numbers, while the last entry is the volume part number. (3) the first six entries of ix6_pfem(7, nwdg_pfem) for each element are the volume node numbers, while the last entry is the volume part number. (4) the first eight entries of ix8_pfem(9, nhex_pfem) for each element are the volume node numbers, while the last entry is the volume part number. 34 LS-DYNA DATABASE When the state data comes from the PFEM_IF_SURFACE domain, then the surface mesh is output first, followed by the data. Currently, the surface mesh is entirely triangular, but we anticipate users will also specify mixed triangle-quadrilateral meshes in the near future: size of each surface variable component: nnpsurf_pfem number of surface nodes: nnpsurf_pfem number of surface elements: nelsurf_pfem user surface node numbers: surfnodes_pfem(nnpsurf_pfem) surface element connectivity: ixsurf_pfem(5, nelsurf_pfem) data for 1st surface variable (size is nnpsurf_pfem) ... data for nsurfvar_pfem-th surface variable (size is nnpsurf_pfem) Notes: (1) the first four entries of ixsurf_pfem(5, nelsurf_pfem) for each surface element are the surface node numbers, while the last entry is the surface part number. The 3rd and 4th node numbers are the same for triangles. When the state data comes from the CESE domain, then the mesh is output first, followed by the data. Currently, the mesh is entirely tetrahedral, but we anticipate users will also specify mixed meshes in the near future: size of each volume variable component: nele_cese number of volume nodes: nnpvol_cese number of tetrahedral elements: ntet_cese number of pyramid elements: npyr_cese number of wedge elements: nwdg_cese number of hexahedral elements: nhex_cese user volume node numbers: volnodes_cese(nnpvol_cese) array of volume nodal coordinates: xvol_cese(3, nnpvol_cese) tetrahedral element connectivity: ix4_cese(5, ntet_cese) pyramid element connectivity: ix5_cese(6, npyr_cese) wedge element connectivity: ix6_cese(7, nwdg_cese) hexahedral element connectivity: ix8_cese(9, nhex_cese) data for 1st volume variable (size is nele_cese) ... data for nvolvar_cese-th volume variable (size is nele_cese) Notes: (1) nele_cese = ntet_cese + npyr_cese + nwdg_cese + nhex_cese (2) the first four entries of ix4_cese(5, ntet_cese) for each element are the volume node numbers, while the last entry is the volume part number. (3) the first five entries of ix5_cese(6, npyr_cese) for each element are the volume node numbers, while the last entry is 35 LS-DYNA DATABASE the volume part number. (4) the first six entries of ix6_cese(7, nwdg_cese) for each element are the volume node numbers, while the last entry is the volume part number. (5) the first eight entries of ix8_cese(9, nhex_cese) for each element are the volume node numbers, while the last entry is the volume part number. When the state data comes from the CESE_SURFACE domain, then the surface mesh is output first, followed by the data. Currently, the surface mesh is entirely triangular, but we anticipate users will also specify mixed triangle-quadrilateral meshes in the near future: size of each surface variable component: nnpsurf_cese number of surface nodes: nnpsurf_cese number of surface elements: nelsurf_cese user surface node numbers: surfnodes_cese(nnpsurf_cese) surface element connectivity: ixsurf_cese(5, nelsurf_cese) data for 1st surface variable (size is nnpsurf_cese) ... data for nsurfvar_cese-th surface variable (size is nnpsurf_cese) Notes: (1) the first four entries of ixsurf_cese(5, nelsurf_cese) for each surface element are the surface node numbers, while the last entry is the surface part number. The 3rd and 4th node numbers are the same for triangles. When the state data comes from the EM domain, then the mesh is output first, followed by the data. Currently, the mesh is entirely tetrahedral, but we anticipate users will also specify mixed meshes in the near future: size of each volume variable component: nnpvol_EM number of volume nodes: nnpvol_EM number of tetrahedral elements: ntet_EM number of pyramid elements: npyr_EM number of wedge elements: nwdg_EM number of hexahedral elements: nhex_EM user volume node numbers: volnodes_EM(nnpvol_EM) array of volume nodal coordinates: xvol_EM(3, nnpvol_EM) tetrahedral element connectivity: ix4_EM(5, ntet_EM) pyramid element connectivity: ix5_EM(6, npyr_EM) wedge element connectivity: ix6_EM(7, nwdg_EM) hexahedral element connectivity: ix8_EM(9, nhex_EM) data for 1st volume variable (size is nnpvol_EM) ... data for nvolvar_EM-th volume variable (size is nnpvol_EM) 36 LS-DYNA DATABASE Notes: (1) the first four entries of ix4_EM(5, ntet_EM) for each element are the volume node numbers, while the last entry is the volume part number. (2) the first five entries of ix5_EM(6, npyr_EM) for each element are the volume node numbers, while the last entry is the volume part number. (3) the first six entries of ix6_EM(7, nwdg_EM) for each element are the volume node numbers, while the last entry is the volume part number. (4) the first eight entries of ix8_EM(9, nhex_EM) for each element are the volume node numbers, while the last entry is the volume part number. When the state data comes from the EM_SURFACE domain, then the surface mesh is output first, followed by the data. Currently, the surface mesh is entirely triangular, but we anticipate users will also specify mixed triangle-quadrilateral meshes in the near future: size of each surface variable component: nnpsurf_EM number of surface nodes: nnpsurf_EM number of surface elements: nelsurf_EM user surface node numbers: surfnodes_EM(nnpsurf_EM) surface element connectivity: ixsurf_EM(5, nelsurf_EM) data for 1st surface variable (size is nnpsurf_EM) ... data for nsurfvar_EM-th surface variable (size is nnpsurf_EM) Notes: (1) the first four entries of ixsurf_EM(5, nelsurf_EM) for each surface element are the surface node numbers, while the last entry is the surface part number. The 3rd and 4th node numbers are the same for triangles. When the state data comes from the STOCHASTIC_PARTICLES domain, then the size of each variable component: n_particles array of particle positions: x_particles(3, n_particles) data for 1st output variable ... data for n_prtcl_vars-th output variable There will always be at least the following two variables output for each particle domain: PARTICLE_SIZES and PARTICLE_VELOCITIES. That is, n_prtcl_vars >= 2. For each particle, both the position and velocity are a 3-component vector. 37 LS-DYNA DATABASE END OF FILE MARKER Value = -999999.0 (a floating point number) 38 LS-DYNA DATABASE TIME HISTORY DATABASE (d3thdt) There are three sections in the LS-DYNA time history database. The first used to contain 144 words of control information, but now depends upon the number of node and elements the user defines in LS-DYNA. The second contains geometric information including the nodal coordinates and element connectivities. The third section contains the results of the analysis at sequential output intervals for a subset of solids, beams, and shells. The output at a given time, called a state, contains a time word, global variables such as total energies and momenta, nodal data consisting of accelerations, velocities, and displacements, and finally element data is written that may include stresses and strains at integration points. The control information provides information on what is in the file and which database is contained. CONTROL DATA VALUE Title Run time File type 10 1 1 DISK #WORDS ADDRESS 0 10 11 DESCRIPTION Model identification time in seconds since 00:00:00 UTC, January 1, 1970 d3thdt=3 1=d3plot, 2=d3drlf, 3=d3thdt, 4=intfor, 5=d3part 6=blstfor, 7=d3cpm, 8=d3ale, 11=d3eigv, 12=d3mode, 13=d3iter, 21=d3ssd, 22=d3spcm, 23=d3psd, 24=d3rms, 25=d3ftg, 26=d3acs Source version Version NDIM 1 1 1 12 13 14 15 ls-dyna version *1000000 + svn number Release number in character*4 form 50 for R5.0 511c for R5.1.1c Code version, a real number, not integer Number of dimensions (2 or 3) is set to 4 if element connectivies are unpacked in the LS-DYNA/3D database and NDIM is reset to 3. NUMNP ICODE NGLBV 39 1 1 1 16 17 18 Number of nodal points Flag to identify finite element code =2 old DYNA3D, NIKE3D database =6 new LS-NIKE3D, LS-DYNA/3D database Number of global variables to be read in each state Release number 1 LS-DYNA DATABASE IT IU IV IA NEL8 NUMMAT8 NDS NST NV3D NEL2 NUMMAT2 NV1D NEL4 NUMMAT4 NV2D NEIPH NEIPS MAXINT NMSPH NGPSPH NARBS BLANK IOSHL(1) IOSHL(2) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 1 1 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 43 44 Flag for temperatures =0 none, =1 read in a temperature for each node Flag for current geometry (=1) Flag for velocities (=1) Flag for accelerations (=1) Number of 8 node solid elements Number of materials used by the 8 node solids Number of node blocks for plotting Number of element blocks for plotting. =NSTH + NSTB + NSTS +NSTT Number of values in database for each solid element Number of 2 node one-dimensional elements Number of materials used by the 2 node 1D elements Number of values in database for each 1D element Number of four node two-dimensional elements Number of materials used by the 4 node 2D elements Number of values in database for each 2D element Number of additional values per solid element to be written in the type 6 database =NEIPH-6*ISTRN Number of additional values per integration point to be written into the type 6 database for shell elements Number of integration points dumped for each shell element. Number of SPH Nodes Number of SPH materials Additional storage required for arbitrary node and element numbering in type 6 database Unused space stress components flag (=1000 yes) Strain components, ISTRN (=1000 yes) 40 IOSHL(3) IOSHL(4) BLANK NCFDV1 NCFDV2 BLANK NSTP IFLAGD NSTH NSTB NSTS NSTT NDSB NSTHB NSTBB NSTSB NSTTB NSTPB 1 1 1 1 1 8 1 1 1 1 1 1 45 46 47 48 49 50 58 59 60 61 62 63 64 LS-DYNA DATABASE Shell force resultants (=1000 yes) Shell thickness, energy + 2 others (=1000 yes) Unused space, ignore value Bit flags for CFD nodal values Further bit flags for CFD nodal values Unused space Number of SPH element blocks Number of node blocks + 1000 flag Number of solid element blocks Number of beam element blocks Number of shell element blocks Number of thick shell element block Node blocks for which time histories are output. The locations 2n-1, where n=1 through NDS correspond to the first node in the block and locations 2n correspond to the last node in the block. 2*NDS 2*NSTH 64+2*NDS Solid element blocks start and end numbers These are defined in a similar manner to the nodal time history blocks 2*NSTB 64+2*NDS Beam element block start and end numbers +2*NSTH 2*NSTS 64+2*NDS Shell element block start and end numbers +2*NSTH +2*NSTB 2*NSTT 2*NSTP 64+2*NDS Thick shell element block start and end numbers +2*NSTH +2*NSTB +2*NSTS 64+2*NDS SPH element block start and end numbers +2*NSTH +2*NSTB +2*NSTS +2*NSTT ISTRN can only be computed as follows and if NSTS > 0. If NV2D-MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS)+8*IOSHL(3)+4*IOSHL(4) > 10 Then ISTRN = 1, else ISTRN = 0 Or NSTT > 0 If NV3DT-MAXINT*(6*IOSHL(1)+IOSHL(2)+NEIPS) > 10 Then ISTRN = 1, else ISTRN = 0 41 LS-DYNA DATABASE SMOOTH PARTICLE HYDRODYNAMICS ELEMENT DATA FLAGS This section is only output if NMSPH > 0. The section is a list of flags to indicate what SPH data is output for each SPH node/element. The first number is the length in words for this array, currently = 10. SPH elements are centered at nodes, and cover a spherical volume defined by the radius of influence. They do not have a connectivity with other SPH elements. They should be displayed as a dot or a spherical surface, with radius scaling to reduce the size and enable each element to be distinguishable. As follows: If the value of isphfg(2-10) = 0, then the particular data item is not output for the particle. To calculated the size of data add the isphfg values from isphfg(2) through isphfg(10) and add one. One value is always output which is the material number as a floating point number for each particle. If this value is negative then the particle has been deleted from the model. Full output for each particle is: mat#, radius, pressure, {sx, sy, sz, sxy, syz, sxz} ps, rho, ie, nn, {ex, ey, ez, exy, eyz, exz}, mass. Hence total size is 20. When a particle is deleted from the model, data is still output for it because the length of data must always be the same for each state. isphfg(1) = 10 - length of sph flags array isphfg(2) = 1 - radius of influence isphfg(3) = 1 - pressure in particle isphfg(4) = 6 - 6 true stress components isphfg(5) = 1 - plastic strain, > 0.0 if effective stress exceeds yield strength isphfg(6) = 1 - density of particle material isphfg(7) = 1 - internal energy (strain) isphfg(8) = 1 - number of neighbors affecting particle isphfg(9) = 6 - 6 true strain components isphfg(10)=1 - mass of element 42 LS-DYNA DATABASE GEOMETRY DATA The geometry section contains the nodal coordinates and the element connectivities. The ordering of the nodal points is assumed to be the same as the ordering of the nodal data in the state data that follows. The connectivities are assumed to be packed with 3 integers per word unless NDIM is set to 4 as in the new LS-DYNA/3D, LS-NIKE3D databases. The order of the elements are 3, 2, and 1 dimensional elements if the database is ICODE=2 or 6. VALUE X(3,1) IX8(9,1) IXT(9,1) IX2(6,1) IX4(5,1) LENGTH NDIM*NUMNP 9*NEL8 9*NELT 6*NEL2 5*NEL4 DESCRIPTION Array of nodal coordinates X1,Y1,Z1, X2,Y2,Z2,X3,Y3,Z3, ... ,Xn,Yn,Zn Connectivity and material number for each 8 node solid element Connectivity and material number for each 8 node thick shell element Connectivity, orientation node, two null entries, and the material number for each 2 node beam element Connectivity and material number for each 4 node shell element 43 LS-DYNA DATABASE USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS Skip this section if NARBS (disk address 39) is zero. The user node and element numbers must be in ascending order. It is assume that if this option is used all the node and element data in the databases is in ascending order in relation to the user numbering. The total length of the data in this database is equal to NARBS=10+NUMNP+NEL8+NEL2+NEL4+NELT if sequential numbering is used of the materials. For arbitrary material numbering the total length is increased by 6+NUMMAT8+NUMMAT4+NUMMAT2+NUMMATT. Material numbers are not in ascending order. VALUE NSORT LENGTH 1 DESCRIPTION Pointer to arbitrary node numbers in LS-DYNA/3D source code, If < 0, it flags that arbitrary material identification numbers are also used. Pointer to arbitrary solid element numbers in LS-DYNA source code: =NSORT+NUMNP Pointer to arbitrary beam element numbers in LS-DYNA source code: =NSRH+NEL8 Pointer to arbitrary shell element numbers in LS-DYNA source code: =NSRB+NEL2 Pointer to arbitrary thick shell element numbers in LS-DYNA source code: =NSRS+NEL4 Number of nodal points Number of 8 node solid elements Number of 2 node beam elements Number of 4 node shell elements Number of 8 node thick shell elements Pointer to an array in the LS-DYNA source code that list the material ID’s in ascending order. Pointer to an array in the LS-DYNA source code that gives the material ID’s in the actual order that they are defined in the user input 44 NSRH NSRB NSRS NSRT NSORTD NSRHD NSRBD NSRSD NSRTD NSRMA NSRMU 1 1 1 1 1 1 1 1 1 1 1 VALUE NSRMP LENGTH 1 DESCRIPTION LS-DYNA DATABASE Pointer to an array in the LS-DYNA source code that gives the location of a member in the array originating at NSRMU for each member in the array starting at NSRMA. Total number of materials Total number of nodal rigid body constraint sets. Total number of materials Array of user defined node numbers Array of user defined solid element numbers Array of user defined beam element numbers Array of user defined shell element numbers Array of user defined solid shell numbers Ordered array of user defined material ID’s Unordered array of user material ID’s Cross reference array NSRTM NUMRBS NMMAT NUSERN NUSERH NUSERB NUSERS NUSERT NORDER NSRMU NSRMP 1 1 1 NSORTD NSORTH NSORTB NSORTS NSORTT NMMAT NMMAT NMMAT 45 LS-DYNA DATABASE TIME HISTORY DATA The time database contains the following data: VALUE TIME GLOBAL LENGTH 1 NGLBV DESCRIPTION Time word Global variables for this state. LS-DYNA Global Variables: KE, IE, TE, X, Y, and Z velocity IE for each material KE for each material X, Y, and Z velocity for mat 1 ... X, Y, and Z velocity for mat n Mass for each material Force for each rigid wall = 6 + 7 * (NUMMAT8+ NUMMAT2 + NUMMAT4 + NUMMATT + NUMRBS) + N*NUMRW, N=1 or N=4 (ls971) • • • • Time word Node data Node data for solids, thick shells, and shells, respectively Element data for solids, thick shells, beams, and shells, respectively SKIP THE FOLLOWING DATA IF THE NUMBER OF NODE BLOCKS FOR PLOTTING IS ZERO (VALUE NUMDS AT DISK ADDRESS 25) TIME NODEDATA 1 NND Time word Total nodal values for state where NLN=10*TNODS where TNODS is the number of nodes put into database. The database contains TNODS vectors each with up to 10 components: temperature (if IT=1); x, y, and z coordinates; x, y, and z velocities; and x, y, and z accelerations. 46 CFDDATA LS-DYNA DATABASE CFD Bit flag: NCFDV1, bits from right to left Eg Pressure, Resultant Vorticity, and Density NCFDV1=2+32+1024=1058 14 Pressure 15 X Vorticity 16 Y Vorticity 17 Z Vorticity 18 Resultant Vorticity 19 Enstrophy 20 Helicity 21 Stream Function 22 Enthalpy 23 Density 24 Turbulent KE 25 Dissipation 14-20 Eddy Viscosity Bit flag: NCFDV2 2-11 Species 1 through 10 Count number of bits on * NUMNP 47 LS-DYNA DATABASE SKIP THE FOLLOWING DATA IF THE NUMBER OF ELEMENT BLOCKS FOR IS ZERO (VALUE NUMDS AT DISK ADDRESS 26) ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR SOLID ELEMENTS **** TIME SOLIDDATA ENV Total nodal values for solid elements where ENV=56*TBELM where THELM is the total number of solid elements to be put into the database. The data contains THELM vectors each with 56 components ordered as follows: 8 connectivities: x,y,z coordinates for each of the 8 nodes; and, lastly, x,y,z velocities for each of the 8 nodes. 1 Time word For solid elements the database contains (7+NEIPH-6*ISTRN) values per element. One set of global stresses are always put into the database for each solid element followed by NEIPH history values. Only data for elements defined in the time history blocks is output. The ordering of the data follows: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10 . Sigma-x (true stress in the global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable First extra value (if NEIPH>0) Second extra value (if NEIPH >1) Etc. until NEIPH extra values are defined if ISTRN=1 7+NEIPH-5. Epsilon-x 7+NEIPH-4. Epsilon-y 7+NEIPH-3. Epsilon-z 7+NEIPH-2. Epsilon-xy 7+NEIPH-1. Epsilon-yz 7+NEIPH. Epsilon-zx 48 LS-DYNA DATABASE 49 LS-DYNA DATABASE ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR THICK SHELLS **** TIME TSHELLDATA ENV Total nodal values for thick shell elements where ENV=56*TBSEL where TBSEL is the total number of thick shell elements in the database. The data contains TBSEL vectors each with 56 components ordered as follows: 8 connectivities: x,y,z coordinates for each of the 8 nodes; and, lastly, x,y,z velocities for each of the 8 nodes. 1 Time word ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR SHELL ELEMENTS **** TIME SHELLDATA 1 ENVS Time word Total nodal values for shell elements where ENVS=28*TSELM where TSELM is the total number of shell elements in the database. The data contains TSELM vectors each with 28 components ordered as follows: 4 connectivities: x,y,z coordinates for each of the 4 nodes; and, lastly, x,y,z velocities for each of the 4 nodes. ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR THICK SHELLS **** For thick shell elements the database contains NV3DT values per element. Three sets of global stresses are always put into the database for each thick shell and are located at the mid surface, the inner surface, and the outer surface, respectively. If one integration point is used the single state is written three times. If two integration points are used then the mid surface value is taken as the average value. The inner values of the stress are always set to the values at the innermost integration point and likewise for outer values. If no integration point lies at the center, i. e. an even number of integration points through the thickness, a value is computed that is an average of the two integration point lying nearest the mid surface. Only data for elements defined in the time history blocks is output. The ordering of the data follows: 1. 2. 3. Sigma-x (mid surface true stress in global system) Sigma-y Sigma-z 50 4. 5. 6. 7. *. 8. 9. 10. 11. 12. 13. 14. *. 15. 16. 17. 18. 19. 20. 21. *. 21. Sigma-xy Sigma-yz Sigma-zx LS-DYNA DATABASE Effective plastic strain or material dependent variable Define NEIPS additional history values here for midsurface Sigma-x (inner surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for inner surface Sigma-x (outer surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for outer surface Effective plastic strain or material dependent variable *. Define NEIPS additional history values here for outer surface If MAXINT >3 then define an additional (MAXINT-3 )* (6*IOSHL(1) + 1*IOSHL(2) + NEIPS) quantities here *. If ISTRN=1, then define strain components Epsilon (x, y, z, xy, yz, zx) here for inner surface and outer surface 51 LS-DYNA DATABASE ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR BEAM ELEMENTS **** TIME BEAMDATA For beam elements the database contains NV1D=6 values per element. They are: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. Axial force S shear resultant T shear resultant S bending moment T bending moment Torsional resultant RS shear stress TR shear stress Axial stress Plastic strain Axial strain BESV Total element values for beam elements where BESV=NV1D*TBELM. The data contains TBELM vectors each with NV2D values ordered as described below. 1 Time word If there are values output at beam integration points, then NV1D = 6 + 5 * BEAMIP BEAMIP is set in *DATABASE_EXTENT_BINARY 52 LS-DYNA DATABASE ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR SHELL ELEMENTS **** For shell elements the database contains NV2D values per element. If the minimum value of MAXINT is 3, then the stresses are typically located at the mid surface, the inner surface, and the outer surface, respectively. If one integration point is used the stress is written three times. If two integration points are used then the mid surface value is taken as the average value. The inner values of the stress are always set to the values at the innermost integration point and likewise for outer values. If no integration point lies at the center, i. e. an even number of integration points through the thickness, a value is computed that is an average of the two integration point lying nearest the mid surface. Only data for elements defined in the time history blocks is output. The ordering of the data follows: 1. 2. 3. 4. 5. 6. 7. *. 8. 9. 10. 11. 12. 13. 14. *. 15. 16. 17. 18. 19. 20. 21. Sigma-x (mid surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for midsurface Sigma-x (inner surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable Define NEIPS additional history values here for inner surface Sigma-x (outer surface true stress in global system) Sigma-y Sigma-z Sigma-xy Sigma-yz Sigma-zx Effective plastic strain or material dependent variable *. Define NEIPS additional history values here for outer surface If MAXINT >3 then define an additional (MAXINT-3 )* (6*IOSHL(1) + 1*IOSHL(2) + 8*IOSHL(3) + 4*IOSHL(4) + NEIPS) quantities here 53 LS-DYNA DATABASE 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. Bending moment-mx (local shell coordinate system) Bending moment-my Bending moment-mxy Shear resultant-qx Shear resultant-qy Normal resultant-nx Normal resultant-ny Normal resultant-nxy Thickness Element dependent variable Element dependent variable 33. Internal energy (if and only if ISTRN=0) The following quantities are expected if and only if ISTRN=1 33. eps-x (inner surface strain in global system) 34. 35. 36. 37. 38. 39. 40. 41. 42. 43. 44. 45. eps-y eps-z eps-xy eps-yz eps-zx eps-x (outer surface strain in global system) eps-y eps-z eps-xy eps-yz eps-zx Internal energy (if and only if ISTRN=1) 54 LS-DYNA DATABASE ****SKIP THE FOLLOWING IF THERE IS NO DATA FOR SPH ELEMENTS **** TIME SPHDATA 1 SPHV Time word Data for each sph element according to the sph flags SPHV=NSTP*NUM_SPH_DATA 55 LS-DYNA DATABASE INTERFACE FORCE DATABASE There are three sections in the interface force database. The first contains 64 words of control information. The second contains geometric information, i.e. the nodal coordinates and segment connectivities for each segment contained in the master and slave surface definitions. The third section contains the results of the analysis at sequential output intervals. The output at a given time is called a state. The state contains a time word, global variables such as total energies and momenta, nodal data consisting of accelerations, velocities, and displacements, and finally segment data is written that include the pressure and shear stress acting on each segment and nodal forces for each node that defines the segment. The control information that follows provides information as to what is in the file and which database is being processed. CONTROL DATA VALUE Title Run time File type 10 1 1 0 10 11 Problem identification time in seconds since 00:00:00 UTC, January 1, 1970 intfor=4 1=d3plot, 2=d3drlf, 3=d3thdt, 4=intfor, 5=d3part 6=blstfor, 7=d3cpm, 8=d3ale, 11=d3eigv, 12=d3mode, 13=d3iter, 21=d3ssd, 22=d3spcm, 23=d3psd, 24=d3rms, 25=d3ftg, 26=d3acs Source version Version NDIM NUMNP ICODE NGLBV BLANK IU 1 1 1 1 1 1 1 1 14 15 16 17 18 19 20 12 13 ls-dyna version *1000000 + svn number Release number in character*4 form 50 for R5.0 511c for R5.1.1c Code version Insert 4 for LS-DYNA/3D database Number of nodal points Insert 6 for LS-DYNA/3D database Number of global variable to be read Insert zero Flag for current geometry (=1) 56 Release number 1 DISK #WORDS ADDRESS DESCRIPTION IV BLANK BLANK BLANK BLANK BLANK BLANK BLANK BLANK BLANK NUMSG NUMMAT4 NV2D BLANK NARBS BLANK NPRESU NSHEAR NFORCE NGAPC 1 1 1 1 1 1 1 1 1 1 1 1 1 5 1 17 1 1 1 1 21 22 23 24 25 26 27 28 29 30 31 32 33 34 39 40 57 58 59 60 LS-DYNA DATABASE Flag for velocities (default=1) Insert zero Insert zero Insert zero Insert zero Insert zero Insert zero Insert zero Insert zero Insert zero Total number of slave and master segments in sliding interface definitions. = 2 times the number of sliding interfaces. = 16, 17, 21 or 23 and see below values If NV2D is negative then the file is FSIFOR for an ALE model Unused space Additional storage required for arbitrary node and element numbering in type 6 database This number equals the sum of (10+ NUMNP+NEL8+NEL2+NEL4+ NELT) Unused space 1, 2, or 3 (default is 3) default is 1 default is 1 default is 1 If *DATABASE_EXTENT_BINARY is included in the model input the following values apply to the state output: NV2D = max(0,NPRESU) + max(0,NSHEAR)*3+max(0,NFORCE)*12 + max(0,NGAPC)*5 NGLBV: output option for global variables EQ.-1 no (NGLBV=0) EQ.1 on NVELO: output option for nodal velocities 57 LS-DYNA DATABASE EQ.-1 no (IV=0) EQ.1 yes NPRESU: output option for pressures EQ.-1 no EQ.1 output normal interface pressure only EQ.2 output normal interface pressure and peak pressure EQ.3 output normal interface pressure, peak pressure and time to peak pressure NSHEAR: output option for maximum interface shear stress, shear stress in r-direction and s-direction EQ.-1 no EQ.1 yes NFORCE: output option for X-, Y- and Z-force at all nodes EQ.-1 no EQ.1 yes NGAPC: output option for contact gap at all nodes and surface energy density EQ.-1 no EQ.1 yes GEOMETRY DATA The geometry section contains the nodal coordinates and the element connectivities. The ordering of the nodal points is assumed to be the same as the ordering of the nodal data in the state data that follows. VALUE X(3,1) IX4(5,1) LENGTH NDIM*NUMNP 5*NUMSG DESCRIPTION Array of nodal coordinates X1,Y1,Z1, X2,Y2,Z2,X3,Y3,Z3, ... ,Xn,Yn,Zn Connectivity and identification number for each 3 or 4 node interface segment. For sliding interface n the identification number in 2n-1 for the slave surface and 2n for the master surface. 58 LS-DYNA DATABASE USER MATERIAL, NODE, AND ELEMENT IDENTIFICATION NUMBERS Skip this section if NARBS (disk address 39) is zero. The user node and element numbers must be in ascending order. It is assumed that if this option is used all node and element data anywhere in the databases is in ascending order based on user numbering. The total length of the data in this database is equal to NARBS=10+NUMNP+NEL8+NEL2+NEL4+NELT if sequential numbering is used of the materials. For arbitrary material numbering the total length is increased by 6+NUMMAT8+NUMMAT4+NUMMAT2+NUMMATT. Material numbers are not in ascending order. VALUE NSORT LENGTH 1 DESCRIPTION Pointer to arbitrary node numbers in LS-DYNA source code, If < 0, it flags that arbitrary material identification numbers are also used. Pointer to arbitrary solid element numbers in LS-DYNA source code: =NSORT+NUMNP Pointer to arbitrary beam element numbers in LS-DYNA source code: =NSRH+NEL8 Pointer to arbitrary shell element numbers in LS-DYNA source code: =NSRB+NEL2 Pointer to arbitrary thick shell element numbers in LS-DYNA source code: =NSRS+NEL4 Number of nodal points Number of 8 node solid elements Number of 2 node beam elements Number of 4 node shell elements Number of 8 node thick shell elements Pointer to an array in the LS-DYNA source code that list the contact ID’s in ascending order. Pointer to an array in the LS-DYNA source code that gives the contact ID’s in the actual order that they are defined in the user input. NSRH NSRB NSRS NSRT NSORTD NSRHD NSRBD NSRSD NSRTD NSRMA NSRMU 1 1 1 1 1 1 1 1 1 1 1 59 LS-DYNA DATABASE VALUE NSRMP LENGTH 1 DESCRIPTION Pointer to an array in the LS-DYNA source code that gives the location of a member in the array originating at NSRMU for each member in the array starting at NSRMA. Total number of materials Total number of nodal rigid body constraint sets Total number of materials Array of user defined node numbers Array of user defined solid element numbers Array of user defined beam element numbers Array of user defined shell element numbers Array of user defined thick shell numbers Ordered array of user defined contact ID’s Unordered array of user contact ID’s Cross reference array NSRTM NUMRBS NMMAT NUSERN NUSERH NUSERB NUSERS NUSERT NORDER NSRMU NSRMP 1 1 1 NSORTD NSORTH NSORTB NSORTS NSORTT NMMAT NMMAT NMMAT 60 STATE DATA VALUE TIME GLOBAL NODEDATA LENGTH 1 NGLBV NND LS-DYNA DATABASE The state data for the interface forces have three parts: • • • Time word and global data Node data Force data for sliding interface segments DESCRIPTION Time word Global variables for this state Total nodal values for state where NND=(IT+NDIM*(IU+IV))*NUMNP LS-DYNA/3D writes 6 values per node, i.e., the three coordinates and the translational velocities. The Data is put into the database as two vectors: first X(3,NUMNP) and then V(3,NUMNP), respectively. Data for sliding interface segments where the quantity ENN = (16,17,21, or 23)*NV2D. The organization of the segment data is described below. SEGMDATA ENN This state data is repeated for each state in the database. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. Normal interface pressure acting on segment Maximum interface shear stress acting on segment Shear stress in local r-direction of segment Shear stress in local s-direction of segment X force at node n1 of segment Y force at node n1 of segment Z force at node n1 of segment X force at node n2 of segment Y force at node n2 of segment Z force at node n2 of segment X force at node n3 of segment Y force at node n3 of segment 61 For each sliding interface segment the database contains NV2D=16, 17, 21 or 23 values per segment. The data order is: LS-DYNA DATABASE 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. Z force at node n3 of segment X force at node n4 of segment Y force at node n4 of segment Z force at node n4 of segment contact gap at node n1 contact gap at node n2 contact gap at node n3 contact gap at node n4 surface energy density of segment peak pressure time to peak pressure If contact gap at node n1 = -1.0, then no values are set. Similarly, if surface energy density is = -1.0, no value is set. If contact gap is >= 999.0 ignore the value, this means that the interfaces are not in discernible contact. Note: original interface force files have 16 variables, while recent ones, from ls-dyna ls970 and ls971, have up to 23. Check size of NV2D. If *database_extent_intfor is include in the model input then NV2D is according to the formula above. FSIFOR file output: 1. Normal interface pressure acting on segment 2. 3. 4. 5. 6. 7. 8. BLSTFOR file (NV2D=16 or 7) output: 1. effective (combined incident and reflected) pressure applied to the segment 2. reflected wave (relevant only for BLAST=4) =-1: segment is below ground level and not exposed to blast = 0 : segment has not been subjected to blast waves = 1 : segment has been subjected to the initial incident wave 62 X force on segment Y force on segment Z force on segment relative interface velocity X interface velocity Y interface velocity Z interface velocity LS-DYNA DATABASE = 2 : segment has been subjected to the ground reflected wave = 3 : segment resides in the Mach stem region (Since known by geometrical considerations this value is fixed at time t=0. Thus, any non-zero reflected pressure on this segment is due to the Mach wave.) 3. incident pressure 4. mass density of air 5. global x-velocity of blast wind 6. global y-velocity of blast wind 7. global z-velocity of blast wind 63 LS-DYNA DATABASE CRACK FILE (d3crck) The crack file is created in LS-DYNA when the Winfrith Concrete material model is used for solid elements. This model allows up to three orthogonal crack planes to develop each with an origin at the center of the element. The plane is assumed to project to the surface of the solid, and can be represented by drawing a line on any of the six solid element faces where it emerges. This cutting line can be found by considering the intersection of each plane and each solid face. The crack file is written as a Fortran unformatted binary file, and each record in the file has a start record mark and an end record mark, each of which is 4 bytes for a single precision run and 8 bytes for a double precision run. The data is output for a state and contains: VALUE TIME 1 1 State Time word NC, Number of sets of crack data Sets of data for each crack NUMBER OF CRACKS CRACK DATA Data for each crack contains: Element ID, Flag for each crack plane, Normal vector for each crack plane and width for each crack plane. Crack plane flags are as follows: 0 = no crack, 1 = cracked, but no sustaining tensile load, 2 = cracked but closed up, and 3 = fully cracked with no tensile strength. Shown below is some C coding from LS-PREPOST to illustrate how the crack data is used and the crack lines established. The important coding is picked out in bold type. Also, after the crack data was read in, the crack flags (3 words) where stored as bits in one word. NG and MAT are the Group number for elements and the Material ID we assign in LS-PREPOST. /* elm id, 3 2bit flags 0 -> 2, group no., mat no. */ int id, pflag, ng, mat; float abc[9]; /* plane normals */ float cp[3]; /* crack width planes 1, 2 and 3 */ /* check brick element face against crack plane (nx,ny,nz) positioned * at element center (xe,ye,ze) */ int CrackPlane(float xe, float ye, float ze, float nx, float ny, float nz, float xyz[][3], float xs[2], float ys[2], float zs[2]) #WORDS DESCRIPTION 16 * NC 64 { int i, k; float x0, y0, z0, x1, y1, z1; float xp1, yp1, zp1, xp2, yp2, zp2, pn, dn, t; k = 0; x0 = xyz[3][0]; y0 = xyz[3][1]; z0 = xyz[3][2]; for (i=0; i<4; i++) { /* parametric line clip algorithm */ xp1 = x0 - xe; yp1 = y0 - ye; zp1 = z0 - ze; x1 = xyz[i][0]; y1 = xyz[i][1]; z1 = xyz[i][2]; xp2 = x1 - x0; yp2 = y1 - y0; zp2 = z1 - z0; x0 = x1; y0 = y1; z0 = z1; pn = nx*xp1 + ny*yp1 + nz*zp1; dn = nx*xp2 + ny*yp2 + nz*zp2; if (fabs(dn) < 1.0e-15) continue; t = -pn / dn; if (t < 0.0 || t > 1.0) continue; t = t - 1.0; xs[k] = x1 + t * xp2; ys[k] = y1 + t * yp2; zs[k] = z1 + t * zp2; k++; if (k > 1) break; } return k; } void SetCrackWidth(float v) { min_crack_width = MAX(0.0, v); } void DrawCracks(int ist, float *bg_color) { int i, k, m, n, nc, nd, kd, ip, ic, is; int id, facecode; int etype, nface, pflag, flag; unsigned int j; float xc, yc, zc, a, b, c; float xyz[24][3]; float xi[2], yi[2], zi[2]; int shrink, count, ns[2]; float dx, dy, dz, ds, d; LS-DYNA DATABASE 65 LS-DYNA DATABASE NDCOOR *nod; int ng=0; float rd, gn, bu; nod = node; nod--; rd = 1.0 - bg_color[0]; gn = 1.0 - bg_color[1]; bu = 1.0 - bg_color[2]; glDisable(GL_LIGHTING); glColor3f(rd, gn, bu); glLineWidth(2.0); GetCrackData(ist); glBegin(GL_LINES); nc = cstate[ist].nc; for (n=0; n>2) & FACEBITS1; if (j < BIT30 && facecode > 0) { id = j; ng = crack[n].ng; if (!part[ng].active) continue; dx = part[ng].dscale[0]; dy = part[ng].dscale[1]; dz = part[ng].dscale[2]; ds = part[ng].dscale[3]; shrink = part[ng].shrink_mode; etype = (active_list[k].akey) & 0x3; nface = FACE_NF[etype]; xc = yc = zc = 0.0; ic = 0; for (i=0; ixyz[0]; xyz[k][1] = (disp_state+kd)->xyz[1]; xyz[k][2] = (disp_state+kd)->xyz[2]; } else { GetScaledNodalCoord(kd, nod, disp_state, dx,dy,dz, xyz[k]); } xc += xyz[k][0]; yc += xyz[k][1]; zc += xyz[k][2]; ic += 1; 66 LS-DYNA DATABASE } } if (ic == 0) continue; d = 1.0 / (float)ic; xc *= d; yc *= d; zc *= d; if (shrink) { for (i=0; i<24; i++) { xyz[i][0] = xc + (xyz[i][0] - xc) * shrink_factor; xyz[i][1] = yc + (xyz[i][1] - yc) * shrink_factor; xyz[i][2] = zc + (xyz[i][2] - zc) * shrink_factor; } } for (i=0; i>is) & 3; if (flag == 0) continue; if ((min_crack_width < 0.5 && crack[n].cp[ip] >= min_crack_width) ||(min_crack_width >= 0.5 && flag == 3)) { a = crack[n].abc[ic]; b = crack[n].abc[ic+1]; c = crack[n].abc[ic+2]; count = CrackPlane(xc, yc, zc, a, b, c, &xyz[k], xi, yi, zi); if (count > 1) { glVertex3f(xi[0], yi[0], zi[0]); glVertex3f(xi[1], yi[1], zi[1]); } } } } } } } glEnd(); glLineWidth(1.0); } 67 LS-DYNA DATABASE DYNAIN BINARY FILE FORMAT (dynain.bin) /* Discription of Dynain binary format: * In first 100 words (integers) * head[0] = location of nodal data * head[1] = number of nodes * head[2] = location of solid element connectivities * head[3] = number of solid elements * head[4] = location of shell element connectivities + thicknesses * head[5] = number of shell elements * head[6] = location of adaptive constraints * head[7] = number of adaptive constraints * head[8] = location of initial stresses for solid elements * head[9] = number of initial stress states defined for solids * head[10] = location of initial stresses for shell elements * head[11] = number of initial stress states defined for shells * head[12] = location of initial strains for shell elements * head[13] = number of initial strains states defined for shells * head[14] = location of boundar spc's * head[15] = number of boundary spc's * head[16] = location of local coordinate systems by nodes * head[17] = number of local coordinate systems by nodes * head[18] = location of local coordinate systems by vector * head[19] = number of local coordinate systems by vector * head[20] = location of initial stress states for beams * head[21] = number of initial stress states for beams * head[22] = location of thick shell element connectivities * head[23] = number of thick shell elements * head[24] = location of initial stresses for thick shell elements * head[25] = number of initial stress states defined for thick shells * head[26] = location of beam element connectivities * head[27] = number of beam elements * head[28] = location of initial strains for solid elements * head[29] = number of initial strain states defined for solids */ 68 LS-DYNA DATABASE EXTRA DATA TYPE DEFINITIONS (NCFDV1 = 67108864) #ifndef _HAVE_D3PLOT #define _HAVE_D3PLOT 1 #define D3PL_FIRST_SCALAR_ID 0 #define D3PL_FIRST_VECTOR_ID 1000 #define D3PL_FIRST_TENSOR_ID 2000 #define D3PL_END_IDS 3000 /* scalar variable names */ enum { D3PL_Pressure_INS=0, D3PL_Temperature_INS, D3PL_Enstrophy_INS, D3PL_Helicity_INS, D3PL_Stream_function_INS, D3PL_Enthalpy_INS, D3PL_Turbulent_KE_INS, D3PL_Turbulent_eps_INS, D3PL_Eddy_Viscosity_INS, D3PL_Density_INS, D3PL_VolFractSpec1_INS, D3PL_VolFractSpec2_INS, D3PL_VolFractSpec3_INS, D3PL_VolFractSpec4_INS, D3PL_VolFractSpec5_INS, D3PL_VolFractSpec6_INS, D3PL_VolFractSpec7_INS, D3PL_VolFractSpec8_INS, D3PL_VolFractSpec9_INS, D3PL_VolFractSpec10_INS, D3PL_Density_CESE, D3PL_Pressure_CESE, D3PL_Temperature_CESE, D3PL_Total_energy_CESE, D3PL_Internal_energy_CESE, D3PL_Enthalpy_CESE, D3PL_Entropy_CESE, D3PL_Stream_function_CESE, D3PL_Density_TS_CESE, D3PL_Total_energy_TS_CESE, D3PL_Temperature_radflow, D3PL_Intensity_radflow, D3PL_Scalar_potential, D3PL_Electrical_conductivity, D3PL_Ohm_heating_power_FEM, D3PL_Ohm_heating_power_BEM, D3PL_Temperature_PFEM, D3PL_Pressure_PFEM, D3PL_K_PFEM, D3PL_eps_PFEM, D3PL_particle_size, D3PL_particle_temperature, D3PL_particle_cnt_child_particles, D3PL_Vorticity_PFEM, D3PL_Cp_PFEM, D3PL_Qc_PFEM, D3PL_Shear_PFEM, D3PL_void_fraction_CESE, D3PL_schlieren_number_CESE, 69 LS-DYNA DATABASE D3PL_LEVELSET_PFEM }; /* vector variable names */ enum { D3PL_Velocity_INS=1000, D3PL_Vorticity_INS, D3PL_Velocity_CESE, D3PL_Vorticity_CESE, D3PL_Momentum_CESE, D3PL_Momentum_TS_CESE, D3PL_E_field_radflow, D3PL_H_field_radflow, D3PL_Current_density_FEM, D3PL_Electric_field_FEM, D3PL_Magnetic_field_FEM, D3PL_Lorentz_force_FEM, D3PL_Vector_potential_FEM, D3PL_Current_density_BEM, D3PL_Electric_field_BEM, D3PL_Magnetic_field_BEM, D3PL_Lorentz_force_BEM, D3PL_Vector_potential_BEM, D3PL_Surface_current, D3PL_Surface_magnetic_field, D3PL_Surface_Lorentz_force, D3PL_Velocity_PFEM, D3PL_Vorticity_vect_PFEM, D3PL_particle_velocity }; /* symmetric tensor variable names */ enum { D3PL_INS_VELOCITY_GRAD=2000 }; typedef struct _d3pnt { char * name; int id; } D3PLOT_NAME_TABLE; /* Identifiers for solver-mesh combinations */ enum { FEM_Q1Q0_INS_CFD=0, CESE_CFD_NODE, CESE_CFD_ELEMENT, CESE_CFD_ELEMENT_TS, RADFLOW_FULL, RADFLOW_NODE, EM_FEMSTER_SOLID_INTEG_PTS, EM_FEMSTER_TSHELL_INTEG_PTS, EM_FEMSTER_SHELL_INTEG_PTS, EM_FEMSTER_SOLID_CENTROID, EM_FEMSTER_TSHELL_CENTROID, EM_FEMSTER_SHELL_CENTROID, EM_FEMSTER_AIR, RECT_AIR_EM_NODE, 70 LS-DYNA DATABASE EM_FEMSTER_BEM, PFEM_IF, PFEM_IF_SURFACE, STOCHASTIC_PARTICLES, CESE, CESE_SURFACE, EM, EM_SURFACE }; static D3PLOT_NAME_TABLE d3plot_solver_name[] = { {\"Incompressible FEM CFD\",FEM_Q1Q0_INS_CFD}, {\"CESE CFD node\",CESE_CFD_NODE}, {\"CESE CFD element\",CESE_CFD_ELEMENT}, {\"CESE CFD taylor series\",CESE_CFD_ELEMENT_TS}, {\"Radiation transport (w/groups)\",RADFLOW_FULL}, {\"Radiation transport\",RADFLOW_NODE}, {\"EM solid integ. pts\",EM_FEMSTER_SOLID_INTEG_PTS}, {\"EM tshell integ. pts\",EM_FEMSTER_TSHELL_INTEG_PTS}, {\"EM shell integ. pts\",EM_FEMSTER_SHELL_INTEG_PTS}, {\"EM solid centroid\",EM_FEMSTER_SOLID_CENTROID}, {\"EM tshell centroid\",EM_FEMSTER_TSHELL_CENTROID}, {\"EM shell centroid\",EM_FEMSTER_SHELL_CENTROID}, {\"EM air\",EM_FEMSTER_AIR}, {\"EM air - rectangular grid\",RECT_AIR_EM_NODE}, {\"EM BEM\",EM_FEMSTER_BEM}, {\"Incompressible CFD\",PFEM_IF}, {\"Incomp. CFD surfaces\",PFEM_IF_SURFACE}, {\"Stochastic particles\",STOCHASTIC_PARTICLES}, {\"CESE compressible CFD\",CESE}, {\"Comp. CFD surfaces\",CESE_SURFACE}, {\"EM nodes\",EM}, {\"EM surface nodes\",EM_SURFACE} }; static D3PLOT_NAME_TABLE d3plot_et_name[] = { {\"Pressure\",D3PL_Pressure_INS}, {\"Temperature\",D3PL_Temperature_INS}, {\"Enstrophy\",D3PL_Enstrophy_INS}, {\"Helicity\",D3PL_Helicity_INS}, {\"Stream function\",D3PL_Stream_function_INS}, {\"Enthalpy\",D3PL_Enthalpy_INS}, {\"Turbulent KE\",D3PL_Turbulent_KE_INS}, {\"Turbulent eps\",D3PL_Turbulent_eps_INS}, {\"Eddy Viscosity\",D3PL_Eddy_Viscosity_INS}, {\"Density\",D3PL_Density_INS}, {\"Volume fraction-1\",D3PL_VolFractSpec1_INS}, {\"Volume fraction-2\",D3PL_VolFractSpec2_INS}, {\"Volume fraction-3\",D3PL_VolFractSpec3_INS}, {\"Volume fraction-4\",D3PL_VolFractSpec4_INS}, {\"Volume fraction-5\",D3PL_VolFractSpec5_INS}, {\"Volume fraction-6\",D3PL_VolFractSpec6_INS}, {\"Volume fraction-7\",D3PL_VolFractSpec7_INS}, {\"Volume fraction-8\",D3PL_VolFractSpec8_INS}, {\"Volume fraction-9\",D3PL_VolFractSpec9_INS}, {\"Volume fraction-10\",D3PL_VolFractSpec10_INS}, {\"Fluid_velocity\",D3PL_Velocity_INS}, {\"Vorticity\",D3PL_Vorticity_INS}, {\"grad(velocity)\",D3PL_INS_VELOCITY_GRAD} {\"Density\",D3PL_Density_CESE}, {\"Pressure\",D3PL_Pressure_CESE}, {\"Temperature\",D3PL_Temperature_CESE}, {\"Total energy\",D3PL_Total_energy_CESE}, 71 LS-DYNA DATABASE {\"Enthalpy\",D3PL_Enthalpy_CESE}, {\"Entropy\",D3PL_Entropy_CESE}, {\"Stream function\",D3PL_Stream_function_CESE}, {\"Void fraction\",D3PL_void_fraction_CESE}, {\"Schlieren_number\",D3PL_schlieren_number_CESE}, {\"Density Taylor series\",D3PL_Density_TS_CESE}, {\"Total energy Taylor series\",D3PL_Total_energy_TS_CESE}, {\"Fluid_velocity\",D3PL_Velocity_CESE}, {\"Vorticity\",D3PL_Vorticity_CESE}, {\"Momentum\",D3PL_Momentum_CESE}, {\"Momentum Taylor series\",D3PL_Momentum_TS_CESE}, {\"Temperature radflow\",D3PL_Temperature_radflow}, {\"Intensity radflow\",D3PL_Intensity_radflow}, {\"E-field radflow\",D3PL_E_field_radflow}, {\"H-field radflow\",D3PL_H_field_radflow}, {\"Scalar potential\",D3PL_Scalar_potential}, {\"Electrical conductivity\",D3PL_Electrical_conductivity}, {\"Ohm heating power FEM\",D3PL_Ohm_heating_power_FEM}, {\"Ohm heating power BEM\",D3PL_Ohm_heating_power_BEM}, {\"Current density FEM\",D3PL_Current_density_FEM}, {\"Electric field FEM\",D3PL_Electric_field_FEM}, {\"Magnetic field FEM\",D3PL_Magnetic_field_FEM}, {\"Lorentz force FEM\",D3PL_Lorentz_force_FEM}, {\"Vector potential FEM\",D3PL_Vector_potential_FEM}, {\"Current density BEM\",D3PL_Current_density_BEM}, {\"Electric field BEM\",D3PL_Electric_field_BEM}, {\"Magnetic field BEM\",D3PL_Magnetic_field_BEM}, {\"Lorentz force BEM\",D3PL_Lorentz_force_BEM}, {\"Vector potential BEM\",D3PL_Vector_potential_BEM}, {\"Surface current\",D3PL_Surface_current}, {\"Surface magnetic field\",D3PL_Surface_magnetic_field}, {\"Surface Lorentz force\",D3PL_Surface_Lorentz_force}, {\"Fluid velocity\",D3PL_Velocity_PFEM}, {\"Fluid temperature\",D3PL_Temperature_PFEM}, {\"Fluid pressure\",D3PL_Pressure_PFEM}, {\"Fluid vortcity\",D3PL_Vorticity_PFEM}, {\"Fluid pressure\",D3PL_Pressure_PFEM}, {\"Turbulent K.E.\",D3PL_K_PFEM}, {\"Turbulent eps.\",D3PL_eps_PFEM}, {\"Particle size\",D3PL_particle_size}, {\"Particle velocity\",D3PL_particle_velocity}, {\"Particle temperature\",D3PL_particle_temperature}, {\"# of child particles\",D3PL_particle_cnt_child_particles} {\"Pressure Coefficient\",D3PL_Cp_PFEM}, {\"Q Criterion\",D3PL_Qc_PFEM}, {\"Surface Shear\",D3PL_Shear_PFEM}, {\"Level Set\",D3PL_LEVELSET_PFEM} }; #endif 72

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文