*KEYWORD $========================== *CONTROL_ALE ==================================== *CONTROL_ALE $# dct nadv meth afac bfac cfac dfac efac 0 1 -2 -1.000000 $# start end aafac vfact prit ebc pref nsidebc 0.100 $# ncpl nbkt imascl checkr beamin mmgpref pdifmx dtmufac $# optimpp ialedr bndflx minmas 1 $=========================== S-ALE Mesh ===================================== *ALE_STRUCTURED_MESH $ mshid pid nbid ebid ityp nparts 1 1001 100001 100001 $ nptx npty nptz nid0 lcsid 1001 1002 1003 $------------------------------------------------------------------------------- *ALE_STRUCTURED_MESH_CONTROL_POINTS 1001 $ i1 x2 1 500. 25 24500. *ALE_STRUCTURED_MESH_CONTROL_POINTS 1002 $ i1 x2 1 1500. 21 21500. *ALE_STRUCTURED_MESH_CONTROL_POINTS 1003 $ i1 x2 1 0. 61 60000. $=========================== S-ALE Materials ================================= *ALE_STRUCTURED_MULTI-MATERIAL_GROUP $# mmgname mid eosid pref tnt 1 1 soil 5 air 2 2 vacuum 7 $--------------------( t,mm,s,MPa )--------------------------------------------- *MAT_HIGH_EXPLOSIVE_BURN $# mid ro d pcj beta k g sigy 1 1.63e-9 6930e3 2.10000E4 2.0 9500 4500 10 *EOS_JWL $# eosid a b r1 r2 omeg e0 vo 1 3.74000E53.750000E3 4.15 0.9 0.357.000000E3 1.0 $------------------------------------------------------------------------------- *MAT_SOIL_AND_FOAM $# mid ro g bulk a0 a1 a2 pc 5 2.2e-9 869 952 1e-5 8.31e-1 1.0e-6 -0.001 $# vcr ref $# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8 0.00E+00-1.500E-03-4.300E-03-1.010E-02-3.050E-02-5.130E-02-7.260E-02-9.430E-02 $# eps9 eps10 -1.740E-01 -0.5 $# p1 p2 p3 p4 p5 p6 p7 p8 0.00E+00 6.43E-01 1.84E+00 4.33E+00 1.31E+01 2.20E+01 3.12E+01 4.06E+01 $# p9 p10 7.55E+01 2.37E+02 $------------------------------------------------------------------------------- *MAT_NULL $# mid ro pc mu terod cerod ym pr 21.1684e-12 *EOS_IDEAL_GAS $# eosid cv0 cp0 cl cq t0 v0 vc0 2 718.0e6 1005.2e6 298.0 1.0 $# adiab $------------------------------------------------------------------------------- *MAT_VACUUM $# mid rho 7 1e-18 $======================== S-ALE Volume Filling =============================== *ALE_STRUCTURED_MESH_VOLUME_FILLING $ mshid to 1 air $ geometry ALL *ALE_STRUCTURED_MESH_VOLUME_FILLING $ mshid to 1 soil $ geometry in/out NID1 NID2 PLANE 99997 99998 *NODE 99997 0.0000000e+00 0.0000000e+00 32000.0 99998 0.0000000e+00 0.0000000e+00 31000.0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $ mshid to 1 TNT $ geometry in/out boxid BOXCOR 1 *DEFINE_BOX 1 11500. 13500. 10500. 12500. 32000. 33000. $====================== S-ALE Boundary Condition ============================== $ not recommending to use both no inflow/outflow bc and non-reflecting. $ here both listed for illustration purpose only. *BOUNDARY_SALE_MESH_FACE $ option mshid -x +x -y +y -z +z NOFLOW 1 1 1 1 1 1 1 NONREFL 1 1 1 1 1 1 1 $====================== S-ALE Initial Condition =============================== *SET_MULTI-MATERIAL_GROUP_LIST 1 tnt *INITIAL_DETONATION $# pid x y z lt - mmgset 1001 12500 11500 32500 0.01 1 $=============================================================================== *END $=============================================================================== *MAT_SOIL_AND_FOAM $# mid ro g bulk a0 a1 a2 pc 4 2.4e-9 869 952 0.4563 41.541000 -0.0519 -0.585 $# vcr ref $# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8 0.00E+00-1.500E-03-4.300E-03-1.010E-02-3.050E-02-5.130E-02-7.260E-02-9.430E-02 $# eps9 eps10 -1.740E-01-2.080E-01 $# p1 p2 p3 p4 p5 p6 p7 p8 0.0000E+001.7070E+013.7213E+015.9745E+011.1351E+021.7121E+022.4290E+023.7161E+02 $# p9 p10 2.1696E+033.3184E+03 *MAT_SOIL_AND_FOAM $# mid ro g bulk a0 a1 a2 pc 3 2.6e-9 1.1300E4 1.24000E4 1.73 41.541000 -0.0519 -3.0 $# vcr ref $# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8 0.00E+00-1.500E-03-4.300E-03-1.010E-02-3.050E-02-5.130E-02-7.260E-02-9.430E-02 $# eps9 eps10 -1.740E-01-2.080E-01 $# p1 p2 p3 p4 p5 p6 p7 p8 0.00E+00 22.581786 49.228296 79.036248 150.16888226.495312321.338816491.605504 $# p9 p10 2870.145024389.89926 *MAT_SOIL_AND_FOAM $# mid ro g bulk a0 a1 a2 pc 51 2.2e-9 869 952 1e-5 8.31e-1 1.0e-6 -0.001 $# vcr ref $# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8 0.00E+00-1.500E-03-4.300E-03-1.010E-02-3.050E-02-5.130E-02-7.260E-02-9.430E-02 $# eps9 eps10 -1.740E-01-2.080E-01 $# p1 p2 p3 p4 p5 p6 p7 p8 0.00E+008.5351E+001.8606E+012.9872E+015.6758E+018.5607E+011.2145E+021.8580E+02 $# p9 p10 1.0848E+031.6592E+03 $=============================================================================== *END