齊魯醫(yī)學(xué)LS-DYNA動(dòng)力學(xué)分析-炸藥在土中爆炸GUI及APDL例子要點(diǎn)

1、精品醫(yī)學(xué)文檔 精品醫(yī)學(xué)文檔 19/19精品醫(yī)學(xué)文檔 第九章 LS-DYNA動(dòng)力學(xué)分析第四節(jié) LS-DYNA范例解析炸藥在土壤內(nèi)部爆炸分析1問題描述如圖9-1所示，條形炸藥設(shè)置在混凝土板下方一定距離處的土壤介質(zhì)中，引爆炸藥，試分析條形炸藥爆炸后土壤的鼓包運(yùn)動(dòng)及混凝土板的運(yùn)動(dòng)過程。圖9-1炸藥與混凝土板的位置示意圖2建模數(shù)值模型由炸藥、土壤層、空氣和混凝土板4部分組成，其中炸藥、土壤和空氣3種材料采用歐拉網(wǎng)格建模，單元使用多物質(zhì)ALE算法，混凝土板采用拉格朗日網(wǎng)格建模，混凝土板與空氣和土壤材料間采用耦合算法。由于是條形炸藥，在中心線起爆條件下，不考慮端部效應(yīng)時(shí)，可以將模型簡化成平面對稱問題。為了方

2、便建模，采用單層實(shí)體網(wǎng)格建模。數(shù)值模型采用cm-g-?s單位制，具體模型尺寸見圖9-2所示。3分析步驟（1）進(jìn)入ANSYS界面1）啟動(dòng)ANSYS，彈出ANSYS12.0 Launcher窗口。2）在Launcher對話框的Simulation Environment下拉列表框中選擇ANSYS，在License下拉列表框中選擇ANSYS/LS-DYNA。圖9-2模型尺寸圖/mm3）在File Management對話框的Working Directory文本框中輸入E:explosion_undergro-und作為工作目錄（假設(shè)工作目錄為E盤），在Job Name文本框中輸入explosion

3、_underground作為工作文件名，其他選項(xiàng)用默認(rèn)值。單擊Run按鈕，運(yùn)行ANSYS程序，進(jìn)入ANSYS的操作界面。（2）選擇單元類型1）選擇菜單Main Menu：PreprocessorElement TypeAdd/Edit/Delete命令，彈出Element Types對話框。2）單擊Add按鈕，出現(xiàn)Library of Element Types對話框，在Library of Element Types下拉列表中選擇LS-DYNA Explicit和3D SOLID164，單擊OK按鈕關(guān)閉對話框。3）單擊Element Types對話框上的Close按鈕，關(guān)閉對話框。（3）定義

4、材料性能參數(shù)1）選擇菜單Main Menu：PreprocessorMaterial PropsMaterial Models命令，彈出Define Material Model Behavior對話框。2）選擇菜單LS-DYNAEquation of StateGruneisenNull命令，彈出Null Properties for Material Number 1對話框。在文本框中，輸入以下數(shù)據(jù)：DENS=0.99821；C=1.647；S1=1.921；S2=-0.096；GAMAO=0.35。其余選項(xiàng)采用默認(rèn)值。輸入完成后單擊OK按鈕關(guān)閉對話框。3）在Define Material

5、 Model Behavior對話框中單擊MaterialNew Model命令，彈出Define Material ID對話框，在文本框中輸入2，單擊OK按鈕關(guān)閉對話框。4）選擇菜單LS-DYNANonlinearInelasticKinematic HardeningPlastic Kinematic命令，彈出Plastic Kinematic Properties for Material Number 2對話框，在文本框中輸入材料本構(gòu)參數(shù)（由于ANSYS前處理所帶的材料模型中沒有我們所需要的土壤模型，暫用隨動(dòng)硬化塑性材料模型代替，后面將在K文件中對材料模型參數(shù)進(jìn)行修改），單擊OK按鈕關(guān)

6、閉對話框。5）在Define Material Model Behavior對話框中單擊MaterialNew Model命令，彈出Define Material ID對話框，在文本框中輸入3，單擊OK按鈕關(guān)閉對話框。6）選擇菜單LS-DYNAEquation of StaticGruneisenNull命令，彈出Null Properties for Material Number 3對話框。在文本框中，輸入以下數(shù)據(jù)：DENS=0.001252；C=0.03437； GAMAO=1.4。其余選項(xiàng)采用默認(rèn)值，輸入完成后單擊OK按鈕關(guān)閉對話框。7）在Define Material Model B

7、ehavior對話框中單擊MaterialNew Model命令，出現(xiàn)Define Material ID對話框，在文本框中輸入4，單擊OK按鈕關(guān)閉對話框。8）選擇菜單LS-DYNANonlinearInelasticKinematic HardeningPlastic Kinematic命令，彈出Plastic Kinematic Properties for Material Number 4對話框，在文本框中輸入以下數(shù)據(jù)：DENS=2.65；EX=40；NUXY=0.3；Yield Stress=1.0E-03；Tangent Modulus=4.0E-02；Hardening Parm

8、=0.5；Failure Strain=0.8。其余選項(xiàng)采用默認(rèn)值，輸入完成后單擊OK按鈕關(guān)閉對話框。9）在Define Material Model Behavior對話框中選擇MaterialExit命令，關(guān)閉對話框。此時(shí)共定義了4種材料。（4）創(chuàng)建幾何模型、劃分網(wǎng)格1）選擇菜單Utility Menu:WorkPlaneChange Active CS toWorking Plane命令。2）選擇菜單Main Menu：PreprocessorModelingCreateVolumesBlockBy Dimensions命令，彈出Create Block by Dimensions對話框

9、，輸入模型體（包括炸藥、土壤和空氣域）的坐標(biāo)，如圖9-3所示。圖9-3 Create Block by Dimensions對話框3）選擇菜單Utility Menu:WorkPlaneOffset WP by Increments命令，彈出Offset WP面板，將按鈕 下方的角度設(shè)置滑動(dòng)條移到最右，滑動(dòng)條上方的數(shù)值變?yōu)?0，單擊按鈕 ，將視窗中的工作平面坐標(biāo)繞X軸旋轉(zhuǎn)90度。4）在Offset WP面板中的X，Y，Z Offsets文本框中輸入0，0，-6，單擊Apply按鈕，工作平面向Z軸負(fù)向移動(dòng)6cm。5）選擇菜單Main Menu：PreprocessorModelingOperat

10、eBooleansDivideVolu by WrkPlane命令，彈出Divide Vol by Wrk拾取菜單，拾取視圖中的體，單擊OK按鈕，分割體。6）按照步驟4）和步驟5）的方法分別將工作平面移至全局坐標(biāo)點(diǎn)的（0，6，0）和（0，78，0），對體進(jìn)行分割。7）將工作平面沿WX方向平移6cm，單擊Offset WP對話框中的 ，工作平面繞WY軸旋轉(zhuǎn)90度，選擇菜單Main Menu：PreprocessorModelingOperateBooleansDivideVolu by WrkPlane命令，在彈出的Divide Vol by Wr拾取菜單面板中單擊Pick All按鈕，將視圖中

11、的4個(gè)體分割為8個(gè)體。8）選擇菜單Utility Menu:WorkPlaneAlign WP withGlobal Cartesian命令，使工作平面與全局坐標(biāo)系重合。9）選擇菜單Main Menu：PreprocessorModelingCreateVolumesBlockBy Dimensions命令，在彈出Create Block by Dimensions文本框中輸入混凝土板模型體的坐標(biāo)，如圖9-4所示。圖9-4混凝土板的模型尺寸10）創(chuàng)建完成后的實(shí)體全模型如圖9-5所示。圖9-5 實(shí)體全模型示意圖（5）劃分網(wǎng)格1）選擇菜單Main Menu：PreprocessorMeshingM

12、eshTool命令，彈出MeshTool對話框，單擊Elements Attributes選擇欄下側(cè)的Set按鈕，彈出Meshing Attributes對話框，在TYPE下拉列表框中選擇1 SOLID164，在MAT下拉列表框中選擇1，其他欄用默認(rèn)值，單擊OK按鈕，確認(rèn)選擇，關(guān)閉對話框。2）選擇菜單Utility Menu:SelectEntities.命令，彈出Select Entities對話框，將選項(xiàng)依次設(shè)置為Volumes、By Num/Pick、From Full，單擊Apply按鈕，拾取炸藥、土壤和空氣域的體。3）將Select Entities對話框中的選項(xiàng)依次設(shè)置成Areas

13、、Attached to 、Volumes、From Full，單擊Apply按鈕。改變Select Entities對話框中的設(shè)置，依次設(shè)置為Lines、Attached to、Areas、From Full，單擊Apply按鈕后，再單擊Plot按鈕。4）在Mesh Tool對話框中，單擊Lines右側(cè)的Set按鈕，彈出Element Size on拾取菜單，單擊Pick All按鈕，在彈出的對話框的SIZE文本框中輸入1.5（即指定單元長度為1.5cm），單擊OK按鈕，確認(rèn)選擇。5）在MeshTool對話框Mesh右側(cè)的下拉列表框中選定Volumes，并選取Shape選項(xiàng)的Hex和Mapp

14、ed兩個(gè)按鈕。單擊Mesh按鈕，彈出Mesh Volumes拾取菜單，拾取視圖中的炸藥體，如圖9-6所示，單擊OK按鈕，對炸藥體進(jìn)行映射網(wǎng)格劃分。圖9-6炸藥域放大6）用上述同樣方法改變MAT號分別對圖9-7和圖9-8中的體進(jìn)行映射網(wǎng)格劃分。圖9-7土壤域圖9-8空氣域7）依次單擊Select Entities對話框中的Invert、Plot按鈕，視圖中顯示出混凝土板的輪廓線圖。8）在Mesh Toll對話框中，單擊Lines右側(cè)的Set按鈕，彈出Element Size on拾取菜單，單擊Pick All按鈕，在彈出的對話框的SIZE文本框中輸入2（即指定單元長度為2cm），單擊OK按鈕，確

15、認(rèn)選擇，關(guān)閉對話框。9）設(shè)置MAT為4，對混凝土板進(jìn)行映射網(wǎng)格劃分。10）選擇菜單Utility Menu:SelectEverything命令，然后選擇Utility Menu:PlotElements命令，視圖中顯示所有材料模型的有限元網(wǎng)格（由于混凝土板的網(wǎng)格與空氣網(wǎng)格疊加在一起，看上去好像只有一種顏色），如圖9-9所示。圖9-9模型網(wǎng)格（局部）（6）創(chuàng)建PART選擇菜單Main Menu：PreprocessorLS-DYNA OptionsParts Options命令，彈出Parts Data Writen for LS-DYNA對話框。在Options選項(xiàng)組中激活Create al

16、l parts單選按鈕，彈出圖9-10所示的PART信息框。檢查對話框中所有的PART定義是否正確，如有定義錯(cuò)誤，需要重新劃分網(wǎng)格，并用Updata parts選項(xiàng)更新PART定義，直至PART定義正確為止。圖9-10EDPART Command對話框（7）施加約束1）選擇菜單Main Menu：PreprocessorLS-DYNA OptionsConstraintsApplyOn Areas命令，彈出Apply U，ROT on拾取菜單。拾取如圖9-11中所有法線方向與Z軸方向相一致的18個(gè)面，單擊OK按鈕，彈出Apply U，ROT on Areas對話框。在Lab2 DOFs to

17、be constrained列表框中選擇UZ，單擊OK按鈕，關(guān)閉對話框。2）拾取如圖9-11中所示法線方向與X軸方向相重合的5個(gè)面（靠炸藥模型一側(cè)），單擊OK按鈕，在Apply U，ROT on Areas對話框的Lab2 DOFs to be constrained列表框中選擇UX，單擊OK按鈕，關(guān)閉對話框。圖9-11用于約束Z向位移的18個(gè)面（8）設(shè)定分析選項(xiàng)1）選擇菜單Main Menu：SolutionAnalysis OptionsEnergy Options命令，彈出Energy Options對話框，激活Hourglass Energy，Sliding Interface和Ray

18、liegh Energy 3個(gè)單選鈕，取消Stonewall Energy單選鈕，單擊OK按鈕，關(guān)閉對話框。2）選擇菜單Main Menu：SolutionAnalysis OptionsCPU Limit命令，彈出CPU Limit對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。3）選擇菜單Main Menu：SolutionAnalysis OptionsBulk Viscosity命令，彈出Bulk Viscosity對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。（9）設(shè)置求解時(shí)間和時(shí)間步控制1）選擇菜單Main Menu：SolutionTime ControlsSolution Tim

19、e命令，彈出Solution Time for LS-DYNA Explicit對話框，在Terminate at Time文本框中輸入100，單擊OK按鈕確認(rèn)輸入。2）選擇菜單Main Menu：SolutionTime ControlsTime Step Ctrls命令，彈出Spectify Time Step Scaling for LS-DYNA Explicit對話框，在Time step scale factor文本框中輸入0.6，單擊OK按鈕，關(guān)閉對話框。（10）設(shè)置輸出類型和時(shí)間間隔1）選擇菜單Main Menu：SolutionOutput ControlsOutput Fi

20、le Types命令，彈出Specify Output File Types for LS-DYNA Solver對話框，在File Options下拉列表框中選擇Add，在Produce output for.下拉列表框中選擇LS-DYNA，單擊OK按鈕，關(guān)閉對話框。2）選擇菜單Main Menu：SolutionOutput ControlsFile Output FreqTime Step Size命令，彈出Specify File Output Frequency對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。（11）輸出K文件1）選擇菜單Utility Menu：SelectEvery

21、thing命令。2）選擇菜單Main Menu：SolutionWrite Jobname.k命令，彈出Input files to be Written for LS-DYNA對話框，在Write results files for.下拉列表框中選擇LS-DYNA，在Write input files to.文本框中輸入explosion_undergroud.k，單擊OK按鈕，程序?qū)⒃诠ぷ髂夸浵律蒭xplosion_undergroud.k文件。（12）編輯修改explosion_undergroud.k文件1）用記事本打開工作目錄下的explosion_undergroud.k文件。2

22、）對照后面所附的K文件修改explosion_undergroud.k文件里的關(guān)鍵字和參數(shù)，內(nèi)容包括：將用于控制單元算法的*SECTION_SOLID關(guān)鍵字修改為*SECTION_SOLID_ALE關(guān)鍵字，用于材料1、2和3的單元算法定義，材料4的單元算法仍采用*SECTION_SOLID定義；添加用于控制ALE算法的*CONTROL_ALE關(guān)鍵字；添加關(guān)鍵字*CONSTRAINED_LAGRANGE_IN_SOLID和*SET_PART_LIST；添加關(guān)鍵字*ALE_MULTI-MATERIAL_GROUP；添加用于起爆點(diǎn)設(shè)置的*INITIAL_DETONATION關(guān)鍵字；用土壤泡沫材料模型

23、*MAT_SOIL_AND_FOAM關(guān)鍵字替代原先所假定的土壤材料模型*MAT_PLASTIC_KINEMATIC；用炸藥材料模型*MAT_HIGH_EXPLOSIVE_BURN關(guān)鍵字和狀態(tài)方程*EOS_JWL關(guān)鍵字代替所假定的炸藥材料模型和狀態(tài)方程，修改后的材料模型關(guān)鍵字和狀態(tài)方程關(guān)鍵字分別為：炸藥*MAT_HIGH_EXPLOSIVE_BURN和*EOS_JWL，添加本構(gòu)模型中的個(gè)別參數(shù)。3）將修改后的explosion_undergroud.k文件與K文件對照，確認(rèn)無誤后存盤。（13）求解1）在Launcher對話框的Simulation Environment下拉列表框中選擇LS-DY

24、NA Solver，并在Analysis Type下面的單選框中選擇Typical LS-DYNA Analysis。2）在File Management對話框的Working Directory中選擇LS-DYNA的輸入文件存放的路徑，并在Keyword Input File欄中選定explosion_undergroud.k。3）單擊Run按鈕，程序?qū)⒄{(diào)用LS-DYNA 971求解器開始求解。4后處理（1）顯示土壤鼓包運(yùn)動(dòng)過程1）運(yùn)行LS-PREPOST程序，選擇菜單File：OpenBinary Plot命令，打開工作目錄下的d3plot文件，程序?qū)⒆x入結(jié)果文件。2）選擇菜單Misc.R

25、eflectReflect about YZ plane命令，將模型做YZ平面對稱鏡像，如圖9-12所示。圖9-12對稱鏡像3）選擇主菜單區(qū)域FcompMiscVolume fraction mat#2，單擊Apply按鈕，然后單擊動(dòng)畫控制區(qū)域中的動(dòng)畫播放按鈕，程序?qū)⒃趫D形對話框中連續(xù)動(dòng)態(tài)地顯示炸藥爆炸后土壤的鼓包運(yùn)動(dòng)過程，圖9-13ad給出了幾個(gè)不同時(shí)刻爆炸場的形態(tài)。（2）輸出土壤鼓包運(yùn)動(dòng)的速度及位移時(shí)程曲線1）運(yùn)行LS-PREPOST，讀入d3plot結(jié)果數(shù)據(jù)文件。圖9-13爆炸后土壤的鼓包及混凝土板運(yùn)動(dòng)過程2）選擇主菜單SelPar，在Part Selection中選擇H4（只顯示Par

26、t4）。3）選擇主單HistoryNodalY-velocity命令，在圖形對話框中選擇H4對稱面上最下方的一個(gè)節(jié)點(diǎn)（節(jié)點(diǎn)號為89606），如圖9-14所示。單擊Plot按鈕，將彈出Plot Window-1對話框，輸出混凝土板頂點(diǎn)的Y方向速度時(shí)程曲線，如圖9-15所示。圖9-14節(jié)點(diǎn)896064）選擇菜單HistoryNodalY-displacement命令，在圖形對話框中仍選擇節(jié)點(diǎn)89606，單擊Plot按鈕，將彈出Plot Window-1對話框，輸出混凝土板頂點(diǎn)的Y方向位移時(shí)程曲線，如圖9-16所示。圖9-15鼓包運(yùn)動(dòng)的速度時(shí)程曲線圖9-16鼓包運(yùn)動(dòng)的位移時(shí)程曲線5命令流文件/PRE

27、P7！進(jìn)入前處理器ET,1,SOLID164！定義單元類型MP,DENS,1,0.99821！定義材料密度MP,EX,1,0！定義彈性模量MP,NUXY,1,0！定義泊松比MP,DENS,2,1.80MP,EX,2,MP,NUXY,2,MP,DENS,3,0.001252MP,EX,3,0MP,NUXY,3,0MP,DENS,4,2.65MP,EX,4,40MP,NUXY,4,0.3CSYS,4！激活前面定義的數(shù)據(jù)BLOCK,0,255,-250,140,0,1.5,！創(chuàng)建模型體wpro,-90.000000,！將工作平面繞X軸順時(shí)針旋轉(zhuǎn)90度wpoff,0,0,-6！將工作平面向Z軸負(fù)向移動(dòng)

28、6cmVSBW,1！用工作平面分離體WPSTYLE,0！重新設(shè)置工作平面wpoff,0,0,12！將工作平面向Z軸正向移動(dòng)6cmFLST,2,2,6,ORDE,2FITEM,2,2FITEM,2,-3VSBW,P51X！用工作平面分離體wpoff,0,0,72！將工作平面向Z軸正向移動(dòng)78cmFLST,2,3,6,ORDE,3FITEM,2,1FITEM,2,-2FITEM,2,4VSBW,P51X！用工作平面分離體WPSTYLE,0！重新設(shè)置工作平面wpoff,6,0,0！將工作平面沿WX方向平移6cmwpro,-90.000000！將工作平面繞WY軸順時(shí)針旋轉(zhuǎn)90度FLST,2,4,6,O

29、RDE,2FITEM,2,2FITEM,2,-5VSBW,P51X！用工作平面分離體WPSTYLE,0！重新設(shè)置工作平面WPCSYS,-1,0！將當(dāng)前坐標(biāo)系X-Y平面定義為工作平面BLOCK,0,250,78,90,0,1.5,！創(chuàng)建混凝土板MAT,1！對炸藥模型進(jìn)行網(wǎng)格劃分REAL,ESYS,0SECNUM,FLST,5,8,6,ORDE,3FITEM,5,1FITEM,5,6FITEM,5,-12VSEL,S,P51XASLV,SLSLA,SFLST,5,59,4,ORDE,14FITEM,5,1FITEM,5,3FITEM,5,8FITEM,5,-14FITEM,5,16FITEM,5,

30、-17FITEM,5,19FITEM,5,-22FITEM,5,24FITEM,5,-25FITEM,5,27FITEM,5,-33FITEM,5,35FITEM,5,-69CM,_Y,LINELSEL,P51XCM,_Y1,LINECMSEL,_YLESIZE,_Y1,1.5,1！設(shè)置單元長度為1.5MSHAPE,0,3DMSHKEY,1CM,_Y,VOLUVSEL,9CM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2MAT,2！對土壤模型進(jìn)行網(wǎng)格劃分REAL,ESYS,0SECNUM,FLST,5,5

31、,6,ORDE,4FITEM,5,1FITEM,5,6FITEM,5,10FITEM,5,-12CM,_Y,VOLUVSEL,P51XCM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2MAT,3！對空氣模型進(jìn)行網(wǎng)格劃分REAL,ESYS,0SECNUM,FLST,5,2,6,ORDE,2FITEM,5,7FITEM,5,-8CM,_Y,VOLUVSEL,P51XCM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2MAT

32、,4！對混泥土板模型進(jìn)行網(wǎng)格劃分REAL,ESYS,0SECNUM,VSEL,S,2ASLV,SLSLA,SEPLOT！顯示模型材料的有限元網(wǎng)格FLST,5,12,4,ORDE,10FITEM,5,2FITEM,5,4FITEM,5,-7FITEM,5,15FITEM,5,18FITEM,5,23FITEM,5,26FITEM,5,34FITEM,5,70FITEM,5,-71CM,_Y,LINELSEL,P51XCM,_Y1,LINECMSEL,_YLESIZE,_Y1,2,1！設(shè)置單元長度為2CM,_Y,VOLUVSEL,2CM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YV

33、MESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2EDPART,CREATE！創(chuàng)建PARTALLSEL,ALL！選擇所有實(shí)體FLST,2,18,5,ORDE,14FITEM,2,2FITEM,2,-5FITEM,2,28FITEM,2,-29FITEM,2,32FITEM,2,-33FITEM,2,37FITEM,2,-38FITEM,2,40FITEM,2,42FITEM,2,45FITEM,2,-48FITEM,2,50FITEM,2,-51DA,P51X,UZ,！對模型施加Z方向的約束FLST,2,5,5,ORDE,5FITEM,2,9FITEM,2,11FITE

34、M,2,15FITEM,2,19FITEM,2,25DA,P51X,UX,！對靠近炸藥模型一側(cè)施加X方向的約束EDENERGY,1,0,1,1！設(shè)定分析選項(xiàng)EDCPU,0,EDBVIS,1.5,0.06,TIME,100,！設(shè)置求解時(shí)間EDCTS,0,0.6,！設(shè)置時(shí)間步EDOPT,ADD,blank,LSDYNA！設(shè)置輸出類型EDRST,0,EDHTIME,0,！設(shè)置時(shí)間間隔EDDUMP,0,ALLSEL,ALL！選擇所有元素EDWRITE,LSDYNA,explosion_undergroud,k！生成K文件/EXIT！退出ANSYS6K文件*KEYWORD*TITLEEXPLOSION_

35、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

36、8.600000000E+01 1.500000000E+00 0 0 91368 2.480000000E+02 8.800000000E+01 1.500000000E+00 0 0$ SECTION DEFINITIONS $*SECTION_SOLID_ALE1,11*SECTION_SOLID_ALE2,11*SECTION_SOLID_ALE3,11*SECTION_SOLID4,0$ ALE DEFINITIONS $*ALE_MULTI-MATERIAL_GROUP1,12,13,1$*CONSTRAINED_LAGRANGE_IN_SOLID1,2,0,0,0,5,3,00,

37、0,0.15*SET_PART_LIST 1 4*SET_PART_LIST 2 2 3$*CONTROL_ALE2,1,2,-1.0000000,0.0000000,0.0000000,0.00000000.0000000,0.0000000,0.0000000$ MATERIAL DEFINITIONS $*MAT_HIGH_EXPLOSIVE_BURN1,1.93,0.993,0.370,0.000E+00*EOS_JWL1,3.71,7.430E-02,4.15,0.950,0.300,7.000E-02,1.00$*INITIAL_DETONATION1,0.00,0.00,0.00

38、,0.00*INITIAL_DETONATION1,0.00,0.00,2.0,0.00$*MAT_SOIL_AND_FOAM2,1.80,1.601E-04,1.328E+02,3.300E-03,1.310E-07,0.1232,0.00.0,0.00.0,0.050,0.090,0.11,0.15,0.19,0.21,0.220.25,0.300.0,3.420E-02,4.530E-02,6.760E-02,1.270E-01,2.080E-01,2.710E-01,3.920E-015.660E-01,1.230E+00$*MAT_NULL3,0.125E-02,0.00,0.00,

39、0.00,0.00,0.00,0.00*EOS_GRUNEISEN3,0.344,0.00,0.00,0.00,1.40,0.00,0.000.00$*MAT_PLASTIC_KINEMATIC4,2.6500000,4.000E-01,0.3000000,1.000E-03,4.000E-02,0.50000000.0000000,0.0000000,0.8000000$ PARTS DEFINITIONS $*PARTPart 1 for Mat 1 and Elem Type 1 1 1 1 1 0 0 0$*PARTPart 2 for Mat 2 and Elem Type 2 2 2 2 0 0 0 0$*PARTPart 3 for Ma