$ $ Thanks to Paul Du Bois , author of this model $ kindly agrees to share the model with our ls-dyna users. $ *KEYWORD *TITLE EMI 180 deg units mm-ms-g-N-MPa $ $ $ *INCLUDE controls.key *INCLUDE outputs.key *INCLUDE materials.key $ $ $ $ $ detonation in the tip of the detonator $ *INITIAL_DETONATION $# pid x y z lt mmgid 1000001 000.000 000.0000 800.00 00.00 1 *SET_MULTI-MATERIAL_GROUP_LIST 1 compB $ $ $ $ below is the SALE specific input $ $ $ *NODE 600001 0.0000000 0.0000 000.0000 600002 0.0000000 0.0000 000.0000 600003 1.0000000 0.0000 000.0000 600004 0.0000000 1.0000 000.0000 $ 600011 0.0000000 0.0000 000.0000 600012 0.0000000 0.0000 850.0000 600013 0.0000000 0.0000 763.0000 600014 0.0000000 0.0000 800.0000 *DEFINE_COORDINATE_NODES 234 600002 600003 600004 600001 $ *ALE_STRUCTURED_MESH 1 1000001 1000001 1000001 1001 1002 1003 600001 234 *ALE_STRUCTURED_MESH_CONTROL_POINTS 1001 1 -0.00 46 1200.00 *ALE_STRUCTURED_MESH_CONTROL_POINTS 1002 1 0.000 46 1200.00 *ALE_STRUCTURED_MESH_CONTROL_POINTS 1003 1 -0.0000 74 1450.00 $ *ALE_STRUCTURED_MULTI-MATERIAL_GROUP $# mid eosid hgid steel 1200 10 0.0 air 100 100 0.1 soil 30021 30010 10 0.0 compB 1080 1080 0.0 vacuum 140 0.0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $# id - ammg 1 air $# geotype in_opt ALL 0 *ALE_STRUCTURED_MESH_VOLUME_FILLING $# id - ammg 1 steel $# geotype in_opt BOXCOR 0 1000005 *ALE_STRUCTURED_MESH_VOLUME_FILLING $# id - ammg 1 soil $# geotype in_opt CYLINDER 0 600011 600012 315. 315. *ALE_STRUCTURED_MESH_VOLUME_FILLING $# id - ammg 1 compB $# geotype in_opt CYLINDER 0 600013 600014 56.5 56.5 $# *DEFINE_BOX 1000005 -0.00 400.00 00.0 400.00 1050.00 1110.00 $ $ $ *BOUNDARY_SALE_MESH_FACE $# option mid -x +x -y +y -z +z SYM 1 1 0 1 0 0 0 FIXED 1 0 0 0 0 1 0 $ $ $ *END $ $ $ *ALE_STRUCTURED_MESH_VOLUME_FILLING $# id - ammg 1 rha $# geotype in_opt BOXCOR 0 1000005 $# *DEFINE_BOX 1000005 -0.00 250.00 36.0 106.00 -01.00 1.00 p2 11 1 p3 11 1 p4 11 1 p6 11 1 $ $ $ $$ *BOUNDARY_SPC_SET_ID $# cid heading 11 $# nsid cid dofx dofy dofz dofrx dofry dofrz 11 0 0 0 1 *BOUNDARY_SPC_SET_ID $# cid heading 12 $# nsid cid dofx dofy dofz dofrx dofry dofrz 12 0 0 0 1 $ *SET_NODE_GENERAL 11 BOX 11 *SET_NODE_GENERAL 12 BOX 12 $ *DEFINE_BOX $ boxid xmin xmax ymin ymax zmin zmax 11 -1. 380. -1.0 340.00 -1.00 0.5 *DEFINE_BOX $ boxid xmin xmax ymin ymax zmin zmax 12 -1. 380. -1.0 340.00 23.5 24.50