1．問題描述
如圖9-1所示，條形炸藥設置在混凝土板下方一定距離處的土壤介質中，引爆炸藥，試分析條形炸藥爆炸后土壤的鼓包運動及混凝土板的運動過程。
圖9-1炸藥與混凝土板的位置示意圖
2．建模 數值模型由炸藥、土壤層、空氣和混凝土板4部分組成，其中炸藥、土壤和空氣3種材料採用歐拉網格建模，單元使用多物質ALE演算法，混凝土板採用拉格朗日網格建模，混凝土板與空氣和土壤材料間採用耦合演算法。由於是條形炸藥，在中心線起爆條件下，不考慮端部效應時，可以將模型簡化成平面對稱問題。為了方便建模，採用單層實體網格建模。數值模型採用cm-g-?s單位制，具體模型尺寸見圖9-2所示。 3．分析步驟 （1）進入ANSYS界面
1）啟動ANSYS，彈出ANSYS12.0 Launcher窗口。
2）在Launcher對話框的Simulation Environment下拉列表框中選擇ANSYS，在License下拉列表框中選擇ANSYS/LS-DYNA。 
圖9-2模型尺寸圖/mm
3）在File Management對話框的Working Directory文本框中輸入E:explosion_undergro- und作為工作目錄（假設工作目錄為E盤），在Job Name文本框中輸入explosion_underground作為工作文件名，其他選項用默認值。單擊Run按鈕，運行ANSYS程序，進入ANSYS的操作界面。
（2）選擇單元類型
1）選擇菜單Main Menu：Preprocessor→Element Type→Add/Edit/Delete命令，彈出Element Types對話框。 2）單擊Add按鈕，出現Library of Element Types對話框，在Library of Element Types下拉列表中選擇LS-DYNA Explicit和3D SOLID164，單擊OK按鈕關閉對話框。 3）單擊Element Types對話框上的Close按鈕，關閉對話框。 （3）定義材料性能參數 1）選擇菜單Main Menu：Preprocessor→Material Props→Material Models命令，彈出Define Material Model Behavior對話框。 2）選擇菜單LS-DYNA→Equation of State→Gruneisen→Null命令，彈出Null Properties for Material Number 1對話框。在文本框中，輸入以下數據：DENS=0.99821；C=1.647；S1=1.921；S2=-0.096；GAMAO=0.35。其餘選項採用默認值。輸入完成後單擊OK按鈕關閉對話框。 3）在Define Material Model Behavior對話框中單擊Material→New Model命令，彈出Define Material ID對話框，在文本框中輸入2，單擊OK按鈕關閉對話框。 4）選擇菜單LS-DYNA→Nonlinear→Inelastic→Kinematic Hardening→Plastic Kinematic命令，彈出Plastic Kinematic Properties for Material Number 2對話框，在文本框中輸入材料本構參數（由於ANSYS前處理所帶的材料模型中沒有我們所需要的土壤模型，暫用隨動硬化塑性材料模型代替，後面將在K文件中對材料模型參數進行修改），單擊OK按鈕關閉對話框。
5）在Define Material Model Behavior對話框中單擊Material→New Model命令，彈出Define Material ID對話框，在文本框中輸入3，單擊OK按鈕關閉對話框。
6）選擇菜單LS-DYNA→Equation of Static→Gruneisen→Null命令，彈出Null Properties for Material Number 3對話框。在文本框中，輸入以下數據：DENS=0.001252；C=0.03437； GAMAO=1.4。其餘選項採用默認值，輸入完成後單擊OK按鈕關閉對話框。 7）在Define Material Model Behavior對話框中單擊Material→New Model命令，出現Define Material ID對話框，在文本框中輸入4，單擊OK按鈕關閉對話框。 8）選擇菜單LS-DYNA→Nonlinear→Inelastic→Kinematic Hardening→Plastic Kinematic命令，彈出Plastic Kinematic Properties for Material Number 4對話框，在文本框中輸入以下數據：DENS=2.65；EX=40；NUXY=0.3；Yield Stress=1.0E-03；Tangent Modulus=4.0E-02；Hardening Parm=0.5；Failure Strain=0.8。其餘選項採用默認值，輸入完成後單擊OK按鈕關閉對話框。 9）在Define Material Model Behavior對話框中選擇Material→Exit命令，關閉對話框。此時共定義了4種材料。 （4）創建幾何模型、劃分網格 1）選擇菜單Utility Menu:WorkPlane→Change Active CS to→Working Plane命令。 2）選擇菜單Main Menu：Preprocessor→Modeling→Create→Volumes→Block→By Dimensions命令，彈出Create Block by Dimensions對話框，輸入模型體（包括炸藥、土壤和空氣域）的坐標，如圖9-3所示。
圖9-3 Create Block by Dimensions對話框
3）選擇菜單Utility Menu:WorkPlane→Offset WP by Increments命令，彈出Offset WP面板，將按鈕 下方的角度設置滑動條移到最右，滑動條上方的數值變為90，單擊按鈕 ，將視窗中的工作平面坐標繞X軸旋轉90度。 4）在Offset WP面板中的X，Y，Z Offsets文本框中輸入0，0，-6，單擊Apply按鈕，工作平面向Z軸負向移動6cm。 5）選擇菜單Main Menu：Preprocessor→Modeling→Operate→Booleans→Divide→Volu by WrkPlane命令，彈出Divide Vol by Wrk…拾取菜單，拾取視圖中的體，單擊OK按鈕，分割體。 6）按照步驟4）和步驟5）的方法分別將工作平面移至全局坐標點的（0，6，0）和（0，78，0），對體進行分割。 7）將工作平面沿WX方向平移6cm，單擊Offset WP對話框中的 ，工作平面繞WY軸旋轉90度，選擇菜單Main Menu：Preprocessor→Modeling→Operate→Booleans→Divide→Volu by WrkPlane命令，在彈出的Divide Vol by Wr…拾取菜單面板中單擊Pick All按鈕，將視圖中的4個體分割為8個體。 8）選擇菜單Utility Menu:WorkPlane→Align WP with→Global Cartesian命令，使工作平面與全局坐標系重合。
9）選擇菜單Main Menu：Preprocessor→Modeling→Create→Volumes→Block→By Dimensions命令，在彈出Create Block by Dimensions文本框中輸入混凝土板模型體的坐標，如圖9-4所示。
圖9-4混凝土板的模型尺寸 10）創建完成後的實體全模型如圖9-5所示。
圖9-5 實體全模型示意圖
（5）劃分網格 1）選擇菜單Main Menu：Preprocessor→Meshing→MeshTool命令，彈出MeshTool對話框，單擊Elements Attributes選擇欄下側的Set按鈕，彈出Meshing Attributes對話框，在[TYPE]下拉列表框中選擇1 SOLID164，在[MAT]下拉列表框中選擇1，其他欄用默認值，單擊OK按鈕，確認選擇，關閉對話框。 2）選擇菜單Utility Menu:Select→Entities...命令，彈出Select Entities對話框，將選項依次設置為Volumes、By Num/Pick、From Full，單擊Apply按鈕，拾取炸藥、土壤和空氣域的體。 3）將Select Entities對話框中的選項依次設置成Areas、Attached to 、Volumes、From Full，單擊Apply按鈕。改變Select Entities對話框中的設置，依次設置為Lines、Attached to、Areas、From Full，單擊Apply按鈕後，再單擊Plot按鈕。 4）在Mesh Tool對話框中，單擊Lines右側的Set按鈕，彈出Element Size on…拾取菜單，單擊Pick All按鈕，在彈出的對話框的SIZE文本框中輸入1.5（即指定單元長度為1.5cm），單擊OK按鈕，確認選擇。
5）在MeshTool對話框Mesh右側的下拉列表框中選定Volumes，並選取Shape選項的Hex和Mapped兩個按鈕。單擊Mesh按鈕，彈出Mesh Volumes拾取菜單，拾取視圖中的炸藥體，如圖9-6所示，單擊OK按鈕，對炸藥體進行映射網格劃分。
圖9-6炸藥域放大 6）用上述同樣方法改變[MAT]號分別對圖9-7和圖9-8中的體進行映射網格劃分。
圖9-7土壤域
圖9-8空氣域 7）依次單擊Select Entities對話框中的Invert、Plot按鈕，視圖中顯示出混凝土板的輪廓線圖。 8）在Mesh Toll對話框中，單擊Lines右側的Set按鈕，彈出Element Size on…拾取菜單，單擊Pick All按鈕，在彈出的對話框的SIZE文本框中輸入2（即指定單元長度為2cm），單擊OK按鈕，確認選擇，關閉對話框。 9）設置[MAT]為4，對混凝土板進行映射網格劃分。 10）選擇菜單Utility Menu:Select→Everything命令，然後選擇Utility Menu:Plot→Elements命令，視圖中顯示所有材料模型的有限元網格（由於混凝土板的網格與空氣網格疊加在一起，看上去好像只有一種顏色），如圖9-9所示。
圖9-9模型網格（局部）
（6）創建PART 選擇菜單Main Menu：Preprocessor→LS-DYNA Options→Parts Options命令，彈出Parts Data Writen for LS-DYNA對話框。在Options選項組中激活Create all parts單選按鈕，彈出圖9-10所示的PART信息框。檢查對話框中所有的PART定義是否正確，如有定義錯誤，需要重新劃分網格，並用Updata parts選項更新PART定義，直至PART定義正確為止。 
圖9-10EDPART Command對話框
（7）施加約束 1）選擇菜單Main Menu：Preprocessor→LS-DYNA Options→Constraints→Apply→On Areas命令，彈出Apply U，ROT on…拾取菜單。拾取如圖9-11中所有法線方向與Z軸方向相一致的18個面，單擊OK按鈕，彈出Apply U，ROT on Areas對話框。在Lab2 DOFs to be constrained列表框中選擇UZ，單擊OK按鈕，關閉對話框。 2）拾取如圖9-11中所示法線方向與X軸方向相重合的5個面（靠炸藥模型一側），單擊OK按鈕，在Apply U，ROT on Areas對話框的Lab2 DOFs to be constrained列表框中選擇UX，單擊OK按鈕，關閉對話框。 
圖9-11用於約束Z向位移的18個面
（8）設定分析選項 1）選擇菜單Main Menu：Solution→Analysis Options→Energy Options命令，彈出Energy Options對話框，激活Hourglass Energy，Sliding Interface和Rayliegh Energy 3個單選鈕，取消Stonewall Energy單選鈕，單擊OK按鈕，關閉對話框。 2）選擇菜單Main Menu：Solution→Analysis Options→CPU Limit命令，彈出CPU Limit對話框，使用默認值，單擊OK按鈕，關閉對話框。 3）選擇菜單Main Menu：Solution→Analysis Options→Bulk Viscosity命令，彈出Bulk Viscosity對話框，使用默認值，單擊OK按鈕，關閉對話框。 （9）設置求解時間和時間步控制 1）選擇菜單Main Menu：Solution→Time Controls→Solution Time命令，彈出Solution Time for LS-DYNA Explicit對話框，在Terminate at Time文本框中輸入100，單擊OK按鈕確認輸入。 2）選擇菜單Main Menu：Solution→Time Controls→Time Step Ctrls命令，彈出Spectify Time Step Scaling for LS-DYNA Explicit對話框，在Time step scale factor文本框中輸入0.6，單擊OK按鈕，關閉對話框。 （10）設置輸出類型和時間間隔 1）選擇菜單Main Menu：Solution→Output Controls→Output File Types命令，彈出Specify Output File Types for LS-DYNA Solver對話框，在File Options下拉列表框中選擇Add，在Produce output for...下拉列表框中選擇LS-DYNA，單擊OK按鈕，關閉對話框。 2）選擇菜單Main Menu：Solution→Output Controls→File Output Freq→Time Step Size命令，彈出Specify File Output Frequency對話框，使用默認值，單擊OK按鈕，關閉對話框。 （11）輸出K文件 1）選擇菜單Utility Menu：Select→Everything命令。 2）選擇菜單Main Menu：Solution→Write Jobname.k命令，彈出Input files to be Written for LS-DYNA對話框，在Write results files for...下拉列表框中選擇LS-DYNA，在Write input files to...文本框中輸入explosion_undergroud.k，單擊OK按鈕，程序將在工作目錄下生成explosion_undergroud.k文件。 （12）編輯修改explosion_undergroud.k文件 1）用記事本打開工作目錄下的explosion_undergroud.k文件。 2）對照後面所附的K文件修改explosion_undergroud.k文件里的關鍵字和參數，內容包括： ①將用於控制單元演算法的*SECTION_SOLID關鍵字修改為*SECTION_SOLID_ALE關鍵字，用於材料1、2和3的單元演算法定義，材料4的單元演算法仍採用*SECTION_SOLID定義； ②添加用於控制ALE演算法的*CONTROL_ALE關鍵字； ③添加關鍵字*CONSTRAINED_LAGRANGE_IN_SOLID和*SET_PART_LIST； ④添加關鍵字*ALE_MULTI-MATERIAL_GROUP； ⑤添加用於起爆點設置的*INITIAL_DETONATION關鍵字； ⑥用土壤泡沫材料模型*MAT_SOIL_AND_FOAM關鍵字替代原先所假定的土壤材料模型*MAT_PLASTIC_KINEMATIC； ⑦用炸藥材料模型*MAT_HIGH_EXPLOSIVE_BURN關鍵字和狀態方程*EOS_JWL關鍵字代替所假定的炸藥材料模型和狀態方程，修改後的材料模型關鍵字和狀態方程關鍵字分別為：炸藥——*MAT_HIGH_EXPLOSIVE_BURN和*EOS_JWL，添加本構模型中的個別參數。 3）將修改後的explosion_undergroud.k文件與K文件對照，確認無誤後存檔。 （13）求解 1）在Launcher對話框的Simulation Environment下拉列表框中選擇LS-DYNA Solver，並在Analysis Type下面的單選框中選擇Typical LS-DYNA Analysis。 2）在File Management對話框的Working Directory中選擇LS-DYNA的輸入文件存放的路徑，並在Keyword Input File欄中選定explosion_undergroud.k。 3）單擊Run按鈕，程序將調用LS-DYNA 971求解器開始求解。 4．後處理 （1）顯示土壤鼓包運動過程 1）運行LS-PREPOST程序，選擇菜單File：Open→Binary Plot命令，打開工作目錄下的d3plot文件，程序將讀入結果文件。 2）選擇菜單Misc.→Reflect→Reflect about YZ plane命令，將模型做YZ平面對稱鏡像，如圖9-12所示。 
圖9-12對稱鏡像 3）選擇主菜單區域Fcomp→Misc→Volume fraction mat#2，單擊Apply按鈕，然後單擊動畫控制區域中的動畫播放按鈕，程序將在圖形對話框中連續動態地顯示炸藥爆炸后土壤的鼓包運動過程，圖9-13a～d給出了幾個不同時刻爆炸場的形態。 （2）輸出土壤鼓包運動的速度及位移時程曲線 1）運行LS-PREPOST，讀入d3plot結果數據文件。
圖9-13爆炸后土壤的鼓包及混凝土板運動過程 2）選擇主菜單SelPar，在Part Selection中選擇H4（只顯示Part4）。 3）選擇主單History→Nodal→Y-velocity命令，在圖形對話框中選擇H4對稱面上最下方的一個節點（節點號為89606），如圖9-14所示。單擊Plot按鈕，將彈出Plot Window-1對話框，輸出混凝土板頂點的Y方向速度時程曲線，如圖9-15所示。
圖9-14節點89606 4）選擇菜單History→Nodal→Y-displacement命令，在圖形對話框中仍選擇節點89606，單擊Plot按鈕，將彈出Plot Window-1對話框，輸出混凝土板頂點的Y方向位移時程曲線，如圖9-16所示。
圖9-15鼓包運動的速度時程曲線
圖9-16鼓包運動的位移時程曲線 5．命令流文件 /PREP7 ！進入前處理器 ET,1,SOLID164 ！定義單元類型 MP,DENS,1,0.99821 ！定義材料密度 MP,EX,1,0 ！定義彈性模量 MP,NUXY,1,0 ！定義泊松比 MP,DENS,2,1.80 MP,EX,2, MP,NUXY,2, MP,DENS,3,0.001252 MP,EX,3,0 MP,NUXY,3,0 MP,DENS,4,2.65 MP,EX,4,40 MP,NUXY,4,0.3 CSYS,4 ！激活前面定義的數據 BLOCK,0,255,-250,140,0,1.5, ！創建模型體 wpro,,-90.000000, ！將工作平面繞X軸順時針旋轉90度 wpoff,0,0,-6 ！將工作平面向Z軸負向移動6cm VSBW,1 ！用工作平面分離體 WPSTYLE,,,,,,,,0 ！重新設置工作平面 wpoff,0,0,12 ！將工作平面向Z軸正向移動6cm FLST,2,2,6,ORDE,2 FITEM,2,2 FITEM,2,-3 VSBW,P51X ！用工作平面分離體 wpoff,0,0,72 ！將工作平面向Z軸正向移動78cm FLST,2,3,6,ORDE,3 FITEM,2,1 FITEM,2,-2 FITEM,2,4 VSBW,P51X ！用工作平面分離體 WPSTYLE,,,,,,,,0 ！重新設置工作平面 wpoff,6,0,0 ！將工作平面沿WX方向平移6cm wpro,,,-90.000000 ！將工作平面繞WY軸順時針旋轉90度 FLST,2,4,6,ORDE,2 FITEM,2,2 FITEM,2,-5 VSBW,P51X ！用工作平面分離體 WPSTYLE,,,,,,,,0 ！重新設置工作平面 WPCSYS,-1,0 ！將當前坐標系X-Y平面定義為工作平面 BLOCK,0,250,78,90,0,1.5, ！創建混凝土板 MAT,1 ！對炸藥模型進行網格劃分 REAL, ESYS,0 SECNUM, FLST,5,8,6,ORDE,3 FITEM,5,1 FITEM,5,6 FITEM,5,-12 VSEL,S,,,P51X ASLV,S LSLA,S FLST,5,59,4,ORDE,14 FITEM,5,1 FITEM,5,3 FITEM,5,8 FITEM,5,-14 FITEM,5,16 FITEM,5,-17 FITEM,5,19 FITEM,5,-22 FITEM,5,24 FITEM,5,-25 FITEM,5,27 FITEM,5,-33 FITEM,5,35 FITEM,5,-69 CM,_Y,LINE LSEL,,,,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1,1.5,,,,,,,1 ！設置單元長度為1.5 MSHAPE,0,3D MSHKEY,1 CM,_Y,VOLU VSEL,,,,9 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 MAT,2 ！對土壤模型進行網格劃分 REAL, ESYS,0 SECNUM, FLST,5,5,6,ORDE,4 FITEM,5,1 FITEM,5,6 FITEM,5,10 FITEM,5,-12 CM,_Y,VOLU VSEL,,,,P51X CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 MAT,3 ！對空氣模型進行網格劃分 REAL, ESYS,0 SECNUM, FLST,5,2,6,ORDE,2 FITEM,5,7 FITEM,5,-8 CM,_Y,VOLU VSEL,,,,P51X CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 MAT,4 ！對混泥土板模型進行網格劃分 REAL, ESYS,0 SECNUM, VSEL,S,,,2 ASLV,S LSLA,S EPLOT ！顯示模型材料的有限元網格 FLST,5,12,4,ORDE,10 FITEM,5,2 FITEM,5,4 FITEM,5,-7 FITEM,5,15 FITEM,5,18 FITEM,5,23 FITEM,5,26 FITEM,5,34 FITEM,5,70 FITEM,5,-71 CM,_Y,LINE LSEL,,,,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1,2,,,,,,,1 ！設置單元長度為2 CM,_Y,VOLU VSEL,,,,2 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 EDPART,CREATE ！創建PART ALLSEL,ALL ！選擇所有實體 FLST,2,18,5,ORDE,14 FITEM,2,2 FITEM,2,-5 FITEM,2,28 FITEM,2,-29 FITEM,2,32 FITEM,2,-33 FITEM,2,37 FITEM,2,-38 FITEM,2,40 FITEM,2,42 FITEM,2,45 FITEM,2,-48 FITEM,2,50 FITEM,2,-51 DA,P51X,UZ, ！對模型施加Z方向的約束 FLST,2,5,5,ORDE,5 FITEM,2,9 FITEM,2,11 FITEM,2,15 FITEM,2,19 FITEM,2,25 DA,P51X,UX, ！對靠近炸藥模型一側施加X方向的約束 EDENERGY,1,0,1,1 ！設定分析選項 EDCPU,0, EDBVIS,1.5,0.06, TIME,100, ！設置求解時間 EDCTS,0,0.6, ！設置時間步 EDOPT,ADD,blank,LSDYNA ！設置輸出類型 EDRST,0, EDHTIME,0, ！設置時間間隔 EDDUMP,0, ALLSEL,ALL ！選擇所有元素 EDWRITE,LSDYNA,explosion_undergroud,k ！生成K文件 /EXIT ！退出ANSYS 6．K文件 *KEYWORD *TITLE EXPLOSION_UNDERGROUND.K $ *DATABASE_FORMAT 0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ NODE DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *NODE 1 6.000000000E+00 6.000000000E+00 0.000000000E+00 0 0 2 6.000000000E+00 6.000000000E+00 1.500000000E+00 0 0 3 6.000000000E+00-6.000000000E+00 0.000000000E+00 0 0 ................................................................. 91366 2.480000000E+02 8.400000000E+01 1.500000000E+00 0 0 91367 2.480000000E+02 8.600000000E+01 1.500000000E+00 0 0 91368 2.480000000E+02 8.800000000E+01 1.500000000E+00 0 0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ SECTION DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *SECTION_SOLID_ALE 1,11 *SECTION_SOLID_ALE 2,11 *SECTION_SOLID_ALE 3,11 *SECTION_SOLID 4,0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ ALE DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *ALE_MULTI-MATERIAL_GROUP 1,1 2,1 3,1 $ *CONSTRAINED_LAGRANGE_IN_SOLID 1,2,0,0,0,5,3,0 0,0,0.15 *SET_PART_LIST 1 4 *SET_PART_LIST 2 2 3 $ *CONTROL_ALE 2,1,2,-1.0000000,0.0000000,0.0000000,0.0000000 0.0000000,0.0000000,0.0000000 $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ MATERIAL DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *MAT_HIGH_EXPLOSIVE_BURN 1,1.93,0.993,0.370,0.000E+00 *EOS_JWL 1,3.71,7.430E-02,4.15,0.950,0.300,7.000E-02,1.00 $ *INITIAL_DETONATION 1,0.00,0.00,0.00,0.00 *INITIAL_DETONATION 1,0.00,0.00,2.0,0.00 $ *MAT_SOIL_AND_FOAM 2,1.80,1.601E-04,1.328E+02,3.300E-03,1.310E-07,0.1232,0.0 0.0,0.0 0.0,0.050,0.090,0.11,0.15,0.19,0.21,0.22 0.25,0.30 0.0,3.420E-02,4.530E-02,6.760E-02,1.270E-01,2.080E-01,2.710E-01,3.920E-01 5.660E-01,1.230E+00 $ *MAT_NULL 3,0.125E-02,0.00,0.00,0.00,0.00,0.00,0.00 *EOS_GRUNEISEN 3,0.344,0.00,0.00,0.00,1.40,0.00,0.00 0.00 $ *MAT_PLASTIC_KINEMATIC 4,2.6500000,4.000E-01,0.3000000,1.000E-03,4.000E-02,0.5000000 0.0000000,0.0000000,0.8000000 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ PARTS DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ *PART Part 1 for Mat 1 and Elem Type 1 1 1 1 1 0 0 0 $ *PART Part 2 for Mat 2 and Elem Type 2 2 2 2 0 0 0 0 $ *PART Part 3 for Mat 3 and Elem Type 3 3 3 3 3 0 0 0 $ *PART Part 4 for Mat 4 and Elem Type 4 4 4 4 0 0 0 0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ ELEMENT DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *ELEMENT_SOLID 1 1 2 1 10 18 22 69 76 34 2 1 18 10 9 17 34 76 75 35 3 1 17 9 8 16 35 75 74 36 ......................................................................... 45118 4 90484 90485 89614 89615 91358 91363 90621 90622 45119 4 90485 90486 89613 89614 91363 91368 90620 90621 45120 4 90486 89742 89612 89613 91368 90618 90613 90620 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ COORDINATE SYSTEMS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ LOAD DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ RIGID BOUNDRIES $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ BOUNDARY DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *SET_NODE_LIST 1 0.000 0.000 0.000 0.000 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 20 21 22 31 32 33 ........................................................................ 91347 91348 91349 91350 91351 91352 91353 91354 91355 91356 91357 91358 91359 91360 91361 91362 91363 91364 91365 91366 91367 91368 *BOUNDARY_SPC_SET 1 0 0 0 1 0 0 0 *SET_NODE_LIST 2 0.000 0.000 0.000 0.000 19 23 24 25 26 27 28 29 30 55 59 60 61 62 63 64 ........................................................................ 89478 89479 89480 89481 89605 89606 89607 89608 89609 89610 89611 90487 90619 90744 90745 90746 90747 90748 *BOUNDARY_SPC_SET 2 0 1 0 1 0 0 0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ CONTACT DEFINITIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ CONTROL OPTIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *CONTROL_CPU 0.00 *CONTROL_ENERGY 2 1 2 2 *CONTROL_SHELL 20.0 1 -1 1 2 2 1 *CONTROL_TIMESTEP 0.0000 0.6000 0 0.00 0.00 *CONTROL_TERMINATION 0.300E+04 0 0.00000 0.00000 0.00000 $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ TIME HISTORY $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *DATABASE_BINARY_D3PLOT 20.00 *DATABASE_BINARY_D3THDT 4.000 $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ DATABASE OPTIONS $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ *DATABASE_EXTENT_BINARY 0 0 3 1 0 0 0 0 0 0 4 0 0 0 *END
推薦閱讀：
