$ $ $****************************************************************************** $ $ UNITS: $ length -> millimeters $ time -> milliseconds $ mass -> grams $ force -> Newtons $ stress -> F/A = Newtons/mm^2 = MPa $ $ $****************************************************************************** $ $ $ Properties for Air $ *MAT_NULL $ MID RO PC MU TEROD CEROD YM PR 100, 1.29e-6, 0.0, 0.0, 0.0, 0.0 $ $========1=========2=========3=========4=========5=========6=========7=========8 $ EOS CARDS $========1=========2=========3=========4=========5=========6=========7=========8 $ $ Properties for Air $ atmosphric pressure 0.1 MPa = 1 bar $ *EOS_Linear_Polynomial $ EOSID C0 C1 C2 C3 C4 C5 C6 100 , 0.0, 0.0, 0.0, 0.0, 0.4, 0.4, 0.0 $ e0 v0 0.25, 1.0 $ $ ........................................................................... $ *MAT_Vacuum $ MID RHO 140, 1.29E-10 $ $ ........................................................................ $ $ Properties for TNT $ LS-DYNA User Manual V971 Page $ Units: LS-DYNA - Present (g-mm-ms) $ Density g/cm^3 = 10^(-3) g/mm^3 $ Pressure Mbar = 10^5 MPa $ Speed cm/us = 10^4 mm/ms $ ........................................................................ $ $ Properties TNT Equivalent *MAT_HIGH_EXPLOSIVE_BURN $ MID RO D PCJ BETA 8, 1.631E-3, 0.67174E4, 0.18503E5, 0.0 $ $========1=========2=========3=========4=========5=========6=========7=========8 $ EOS CARDS $========1=========2=========3=========4=========5=========6=========7=========8 $ $ TNT Equivalent *EOS_JWLB $ EOSID A1 A2 A3 A4 A5 8, 490.07e5, 56.868e5, 0.82426e5, 0.00093e5 $ R1 R2 R3 R4 R5 40.713, 9.6754, 2.4335, 0.15564 $ AL1 AL2 AL3 AL4 AL5 0.00, 11.468 $ BL1 BL2 BL3 BL4 BL5 1098.0, -6.5011 $ RL1 RL2 RL3 RL4 RL5 15.614, 2.1593 $ C OMEGA E V0 0.0071e5, 0.30270, 0.06656e5, 1.0 $ $ $ ........................................................................ $ Properties for COMPOSITION B (1.717) $ ........................................................................ $ $ MAT_HIGH_EXPLOSIVE_BURN Parameters: $ • RO = 1.717 g/cm^3 = 1.717E-3 g/mm^3 $ • D = 790 m/s = 7.980e3 mm/msec (detonation velocity) $ • PCJ = 29.5 GPa = 2.95E4 MPa (Chapman-Jouget pressure) $ • BETA = 2.0 (Beta burn flag = 0 -> Beta + programmed burn) $ $ EOS_JWL Parameters $ • A = 524.2 GPa = 5.242E5 MPa $ • B = 7.678 GPa = 7.678E3 MPa $ • R1 = 4.2 $ • R2 = 1.10 $ • OMEG = 0.34 $ • E0 = 8.5 GPa = 8.5E3 MPa $ • V0 = 1.0 $ $ *MAT_HIGH_EXPLOSIVE_BURN $ MID RO D PCJ BETA 1080, 1.717E-3, 7.98e3, 2.95E4, 0.0 $ $========1=========2=========3=========4=========5=========6=========7=========8 $ EOS CARDS $========1=========2=========3=========4=========5=========6=========7=========8 $ ........................................................................ $ COMPOSITION B (1.717) $ ........................................................................ *EOS_JWL $ EOSID A B R1 R2 OMEG E0 V0 1080 , 5.242E5, 7.678E3, 4.2, 1.10, 0.34, 8.5E3, 1.0 $ $ $ $ ........................................................................ $ Elastic Steel $ ........................................................................ *MAT_ELASTIC $ MID RO E PR N COUPLE M ALIAS/RE 1200, 7.85E-3, 205.0e3, 0.29 $ $ ........................................................................ $ Rigid Steel $ ........................................................................ *MAT_RIGID $ MID RO E PR N COUPLE M ALIAS/RE 1220, 7.85E-3, 205.0e3, 0.29 $ CMO CON1 CON2 1.0, 6.0, 7.0 $ LCO or A1 A2 A3 V1 V2 V3 0.0 $ $ ........................................................................ $ MAT005 Norway Sandy Soil $ ........................................................................ $ $ $ bulk modulus derived from L&S : K=density*c**2=2641*4634**2=56700. $ shear modulus derived from L&S : G=76.9 ( density=1674 ) $ G=37347.0 ( density=2641 ) $ Poisson ratio computed from both for MAT_016 $ PR=(3k-2G)/(6k+2G)=0.23 $ $ $ $ *MAT_SOIL_AND_FOAM $ MID RHO G Bulk a0 a1 a2 PC 30015, 1.37E-03, 37347., 56.7e3, 0.00, 0.00, 0.65, -6.9e-3 $ VCR , REF 0.000E+00, 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ EPS1, EPS2, EPS3, EPS4, EPS5, EPS6, EPS7, EPS8 0.000, 0.113, 0.176, 0.247, 0.296, 0.352, 0.395, 0.433 $ EPS9, EPS10 0.467, 0.470 $ P1, P2, P3, P4, P5, P6, P7, P8, 0.00, 14.98, 29.15, 58.18, 98.10, 179.44, 289.44, 450.20 $ P9, P10 650.66 , 800.0 $ $ ----------------------------------------------------------------------------- $ MAT010 Norway Sand $ ----------------------------------------------------------------------------- $ $ MAT005 Sigy = sqrt[3*(A0 + A1*P + A2*P^2] $ MAT010 Sigy = Sig0 + Eh*Eps + (AA1 + AA2*P)*max(p,0.0) $ $ If A0=A1=0 and Sig0=Eh=AA2=0 $ Then AA1 = sqrt(3*A2) $ $ Rafael Density as Norway Sandy Soil *MAT_ELASTIC_PLASTIC_HYDRO_SPALL $ MID RHO G SIGY EH PC FS CHARL 30011, 1.37E-03, 37347., 0.0, 0.0, -6.9e-3, 0.0, 0.0 $ A1 A2 SPALL 1.39642, 0.0, 1.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ EPS1, EPS2, EPS3, EPS4, EPS5, EPS6, EPS7, EPS8 $ EPS9, EPS10, EPS11, EPS12, EPS13, EPS14, EPS15, EPS16 $ ES1, ES2, ES3, ES4, ES5, ES6, ES7, ES8, $ ES9, ES10, ES11, ES12, ES13, ES14, ES15, ES16 $ *MAT_ELASTIC_PLASTIC_HYDRO_SPALL $ MID RHO G SIGY EH PC FS CHARL 30021, 1.37E-03, 37347., 0.0, 0.0, -6.9e-3, 0.0, 0.0 $ A1 A2 SPALL 1.39642, 0.0, 2.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ EPS1, EPS2, EPS3, EPS4, EPS5, EPS6, EPS7, EPS8 $ EPS9, EPS10, EPS11, EPS12, EPS13, EPS14, EPS15, EPS16 $ ES1, ES2, ES3, ES4, ES5, ES6, ES7, ES8, $ ES9, ES10, ES11, ES12, ES13, ES14, ES15, ES16 $ *MAT_ELASTIC_PLASTIC_HYDRO_SPALL $ MID RHO G SIGY EH PC FS CHARL 30022, 1.37E-03, 37347., 0.0, 0.0, -6.9e-3, 0.0, 0.0 $ A1 A2 SPALL 1.39642, 0.0, 3.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ EPS1, EPS2, EPS3, EPS4, EPS5, EPS6, EPS7, EPS8 $ EPS9, EPS10, EPS11, EPS12, EPS13, EPS14, EPS15, EPS16 $ ES1, ES2, ES3, ES4, ES5, ES6, ES7, ES8, $ ES9, ES10, ES11, ES12, ES13, ES14, ES15, ES16 $ *MAT_PSEUDO_TENSOR $ Poisson ratio computed as (3K-2G)/(6K+2G) $ not clear why it is needed as law is deviatoric $ MID RHO G PR 30012, 1.37E-03, 37347., 0.23 $ SIGF, A0, A1, A2, A0F, A1F, B1, PER 6.9e-3, 0,0,0,0,0,0 $ ER $ X1-X8 0.0, 1000.,10000. $ X9-X16 $ YS1-YS8 0.0, 1396.42, 13964.2 $ YS9-YS16 $ $ $ *MAT_SOIL_CONCRETE $ MID RHO G K , LCPV, LCYP 30010, 1.37E-03, 37347., 56.7e3, 781, 782 $ PC OUT B FAIL -6.9e-3, 0 , 0 , 0 $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ $ ----------------------------------------------------------------------------- $ EOS Norway Sand $ ----------------------------------------------------------------------------- $ *EOS_TABULATED_COMPACTION $ EOSID GAMA E0 V0 LCC LCT 30010, 0.0, 0.0, 1.0 $ EV1 EV2 EV3 EV4 EV5 0.000, -0.113, -0.176, -0.247, -0.296 $ EV6 EV7 EV8 EV9 EV10 -0.352, -0.395, -0.433, -0.467, -0.470 $ P1 P2 P3 P4 P5 0.00, 14.98, 29.15, 58.18, 98.10, $ P6 P7 P8 P9 P10 179.44, 289.44, 450.20, 650.66, 800.0 $ T1 T2 T3 T4 T5 0.0, 0.0, 0.0, 0.0, 0.0 $ T6 T7 T8 T9 T10 0.0, 0.0, 0.0, 0.0, 0.0 $ $ K1 K2 K3 K4 K5 56.7e3, 56.7e3, 56.7e3, 56.7e3, 56.7e3 $ K6 K7 K8 K9 K10 56.7e3, 56.7e3, 56.7e3, 56.7e3, 56.7e3 $ $========1=========2=========3=========4=========5=========6=========7=========8 $ $ *DEFINE_CURVE 781,0,0,0,0,0,0,0 0.000, 0.0 0.113, 14.98 0.176, 29.15 0.247, 58.18 0.296, 98.10 0.352, 179.44 0.395, 289.44 0.433, 450.20 0.467, 650.66 0.470, 800.00 *DEFINE_CURVE 782,0,0,0,0,0,0,0 0.000, 0.0 1000., 1396.42 10000., 13964.2 $ $ $========1=========2=========3=========4=========5=========6=========7=========8 $ Hourglass Cards $========1=========2=========3=========4=========5=========6=========7=========8 $ *Hourglass $ HGID IHQ QM IBQ Q1 Q2 QB/VDC QW 10, 1, 1.0e-4, 1, 1.5e-1, 6.0e-2 20, 3, 1.0e-1, 1, 1.5e-1, 6.0e-2 $ $ $