LSDYNA動力學(xué)分析炸藥在土中爆炸GUI及APDL例子

LSDYNA動力學(xué)分析炸藥在土中爆炸GUI及APDL例子LSDYNA動力學(xué)分析炸藥在土中爆炸GUI及APDL例子/LSDYNA動力學(xué)分析炸藥在土中爆炸GUI及APDL例子第九章LS-DYNA動力學(xué)分析第四節(jié)LS-DYNA范例解析——炸藥在土壤內(nèi)部爆炸分析1．問題描述如圖9-1所示，條形炸藥設(shè)置在混凝土板下方一定距離處的土壤介質(zhì)中，引爆炸藥，試分析條形炸藥爆炸后土壤的鼓包運動與混凝土板的運動過程。

圖9-1炸藥與混凝土板的位置示意圖

2．建模

數(shù)值模型由炸藥、土壤層、空氣和混凝土板4部分組成，其中炸藥、土壤和空氣3種材料采用歐拉網(wǎng)格建模，單元使用多物質(zhì)ALE算法，混凝土板采用拉格朗日網(wǎng)格建模，混凝土板與空氣和土壤材料間采用耦合算法。由于是條形炸藥，在中心線起爆條件下，不考慮端部效應(yīng)時，可以將模型簡化成平面對稱問題。為了方便建模，采用單層實體網(wǎng)格建模。數(shù)值模型采用cm-g-?s單位制，具體模型尺寸見圖9-2所示。

3．分析步驟

（1）進(jìn)入ANSYS界面

1）啟動ANSYS，彈出ANSYS12.0Launcher窗口。

2）在Launcher對話框的SimulationEnvironment下拉列表框中選擇ANSYS，在License下拉列表框中選擇ANSYS/LS-DYNA。

圖9-2模型尺寸圖/mm

3）在對話框的WorkingDirectory文本框中輸入E:\explosion_undergro-

und作為工作目錄（假設(shè)工作目錄為E盤），在JobName文本框中輸入explosion_underground作為工作文件名，其他選項用默認(rèn)值。單擊Run按鈕，運行ANSYS程序，進(jìn)入ANSYS的操作界面。

（2）選擇單元類型

1）選擇菜單MainMenu：Preprocessor→ElementType→Add/Edit/Delete命令，彈出ElementTypes對話框。

2）單擊Add按鈕，出現(xiàn)LibraryofElementTypes對話框，在LibraryofElementTypes下拉列表中選擇LS-DYNAExplicit和3DSOLID164，單擊OK按鈕關(guān)閉對話框。

3）單擊ElementTypes對話框上的Close按鈕，關(guān)閉對話框。

（3）定義材料性能參數(shù)

1）選擇菜單MainMenu：Preprocessor→MaterialProps→MaterialModels命令，彈出DefineMaterialModelBehavior對話框。

2）選擇菜單LS-DYNA→EquationofState→Gruneisen→Null命令，彈出NullPropertiesforMaterialNumber1對話框。在文本框中，輸入以下數(shù)據(jù)：DENS=0.99821；C=1.647；S1=1.921；S2=-0.096；GAMAO=0.35。其余選項采用默認(rèn)值。輸入完成后單擊OK按鈕關(guān)閉對話框。

3）在DefineMaterialModelBehavior對話框中單擊Material→NewModel命令，彈出DefineMaterialID對話框，在文本框中輸入2，單擊OK按鈕關(guān)閉對話框。

4）選擇菜單LS-DYNA→Nonlinear→Inelastic→KinematicHardening→PlasticKinematic命令，彈出PlasticKinematicPropertiesforMaterialNumber2對話框，在文本框中輸入材料本構(gòu)參數(shù)（由于ANSYS前處理所帶的材料模型中沒有我們所需要的土壤模型，暫用隨動硬化塑性材料模型代替，后面將在K文件中對材料模型參數(shù)進(jìn)行修改），單擊OK按鈕關(guān)閉對話框。

5）在DefineMaterialModelBehavior對話框中單擊Material→NewModel命令，彈出DefineMaterialID對話框，在文本框中輸入3，單擊OK按鈕關(guān)閉對話框。

6）選擇菜單LS-DYNA→EquationofStatic→Gruneisen→Null命令，彈出NullPropertiesforMaterialNumber3對話框。在文本框中，輸入以下數(shù)據(jù)：DENS=0.001252；C=0.03437；GAMAO=1.4。其余選項采用默認(rèn)值，輸入完成后單擊OK按鈕關(guān)閉對話框。

7）在DefineMaterialModelBehavior對話框中單擊Material→NewModel命令，出現(xiàn)DefineMaterialID對話框，在文本框中輸入4，單擊OK按鈕關(guān)閉對話框。

8）選擇菜單LS-DYNA→Nonlinear→Inelastic→KinematicHardening→PlasticKinematic命令，彈出PlasticKinematicPropertiesforMaterialNumber4對話框，在文本框中輸入以下數(shù)據(jù)：DENS=2.65；EX=40；NUXY=0.3；YieldStress=1.0E-03；TangentModulus=4.0E-02；HardeningParm=0.5；FailureStrain=0.8。其余選項采用默認(rèn)值，輸入完成后單擊OK按鈕關(guān)閉對話框。

9）在DefineMaterialModelBehavior對話框中選擇Material→Exit命令，關(guān)閉對話框。此時共定義了4種材料。

（4）創(chuàng)建幾何模型、劃分網(wǎng)格

1）選擇菜單UtilityMenu:WorkPlane→ChangeActiveCSto→WorkingPlane命令。

2）選擇菜單MainMenu：Preprocessor→Modeling→Create→Volumes→Block→ByDimensions命令，彈出CreateBlockbyDimensions對話框，輸入模型體（包括炸藥、土壤和空氣域）的坐標(biāo)，如圖9-3所示。

圖9-3

CreateBlockbyDimensions對話框

3）選擇菜單UtilityMenu:WorkPlane→OffsetWPbyIncrements命令，彈出OffsetWP面板，將按鈕下方的角度設(shè)置滑動條移到最右，滑動條上方的數(shù)值變?yōu)?0，單擊按鈕，將視窗中的工作平面坐標(biāo)繞X軸旋轉(zhuǎn)90度。

4）在OffsetWP面板中的X，Y，ZOffsets文本框中輸入0，0，-6，單擊Apply按鈕，工作平面向Z軸負(fù)向移動6cm。

5）選擇菜單MainMenu：Preprocessor→Modeling→Operate→Booleans→Divide→VolubyWrkPlane命令，彈出DivideVolbyWrk…拾取菜單，拾取視圖中的體，單擊OK按鈕，分割體。

6）按照步驟4）和步驟5）的方法分別將工作平面移至全局坐標(biāo)點的（0，6，0）和（0，78，0），對體進(jìn)行分割。

7）將工作平面沿WX方向平移6cm，單擊OffsetWP對話框中的，工作平面繞WY軸旋轉(zhuǎn)90度，選擇菜單MainMenu：Preprocessor→Modeling→Operate→Booleans→Divide→VolubyWrkPlane命令，在彈出的DivideVolbyWr…拾取菜單面板中單擊PickAll按鈕，將視圖中的4個體分割為8個體。

8）選擇菜單UtilityMenu:WorkPlane→AlignWPwith→GlobalCartesian命令，使工作平面與全局坐標(biāo)系重合。

9）選擇菜單MainMenu：Preprocessor→Modeling→Create→Volumes→Block→ByDimensions命令，在彈出CreateBlockbyDimensions文本框中輸入混凝土板模型體的坐標(biāo)，如圖9-4所示。

圖9-4混凝土板的模型尺寸

10）創(chuàng)建完成后的實體全模型如圖9-5所示。

圖9-5實體全模型示意圖

（5）劃分網(wǎng)格

1）選擇菜單MainMenu：Preprocessor→Meshing→MeshTool命令，彈出MeshTool對話框，單擊ElementsAttributes選擇欄下側(cè)的Set按鈕，彈出MeshingAttributes對話框，在[TYPE]下拉列表框中選擇1SOLID164，在[MAT]下拉列表框中選擇1，其他欄用默認(rèn)值，單擊OK按鈕，確認(rèn)選擇，關(guān)閉對話框。

2）選擇菜單UtilityMenu:Select→Entities...命令，彈出SelectEntities對話框，將選項依次設(shè)置為Volumes、ByNum/Pick、FromFull，單擊Apply按鈕，拾取炸藥、土壤和空氣域的體。

3）將SelectEntities對話框中的選項依次設(shè)置成Areas、Attachedto、Volumes、FromFull，單擊Apply按鈕。改變SelectEntities對話框中的設(shè)置，依次設(shè)置為Lines、Attachedto、Areas、FromFull，單擊Apply按鈕后，再單擊Plot按鈕。

4）在MeshTool對話框中，單擊Lines右側(cè)的Set按鈕，彈出ElementSizeon…拾取菜單，單擊PickAll按鈕，在彈出的對話框的SIZE文本框中輸入1.5（即指定單元長度為1.5cm），單擊OK按鈕，確認(rèn)選擇。

5）在MeshTool對話框Mesh右側(cè)的下拉列表框中選定Volumes，并選取Shape選項的Hex和Mapped兩個按鈕。單擊Mesh按鈕，彈出MeshVolumes拾取菜單，拾取視圖中的炸藥體，如圖9-6所示，單擊OK按鈕，對炸藥體進(jìn)行映射網(wǎng)格劃分。

圖9-6炸藥域放大

6）用上述同樣方法改變[MAT]號分別對圖9-7和圖9-8中的體進(jìn)行映射網(wǎng)格劃分。

圖9-7土壤域

圖9-8空氣域

7）依次單擊SelectEntities對話框中的Invert、Plot按鈕，視圖中顯示出混凝土板的輪廓線圖。

8）在MeshToll對話框中，單擊Lines右側(cè)的Set按鈕，彈出ElementSizeon…拾取菜單，單擊PickAll按鈕，在彈出的對話框的SIZE文本框中輸入2（即指定單元長度為2cm），單擊OK按鈕，確認(rèn)選擇，關(guān)閉對話框。

9）設(shè)置[MAT]為4，對混凝土板進(jìn)行映射網(wǎng)格劃分。

10）選擇菜單UtilityMenu:Select→Everything命令，然后選擇UtilityMenu:Plot→Elements命令，視圖中顯示所有材料模型的有限元網(wǎng)格（由于混凝土板的網(wǎng)格與空氣網(wǎng)格疊加在一起，看上去好像只有一種顏色），如圖9-9所示。

圖9-9模型網(wǎng)格（局部）

（6）創(chuàng)建PART

選擇菜單MainMenu：Preprocessor→LS-DYNAOptions→PartsOptions命令，彈出PartsDataWritenforLS-DYNA對話框。在Options選項組中激活Createallparts單選按鈕，彈出圖9-10所示的PART信息框。檢查對話框中所有的PART定義是否正確，如有定義錯誤，需要重新劃分網(wǎng)格，并用Updataparts選項更新PART定義，直至PART定義正確為止。

圖9-10EDPARTCommand對話框

（7）施加約束

1）選擇菜單MainMenu：Preprocessor→LS-DYNAOptions→Constraints→Apply→OnAreas命令，彈出ApplyU，ROTon…拾取菜單。拾取如圖9-11中所有法線方向與Z軸方向相一致的18個面，單擊OK按鈕，彈出ApplyU，ROTonAreas對話框。在Lab2DOFstobeconstrained列表框中選擇UZ，單擊OK按鈕，關(guān)閉對話框。

2）拾取如圖9-11中所示法線方向與X軸方向相重合的5個面（靠炸藥模型一側(cè)），單擊OK按鈕，在ApplyU，ROTonAreas對話框的Lab2DOFstobeconstrained列表框中選擇UX，單擊OK按鈕，關(guān)閉對話框。

圖9-11用于約束Z向位移的18個面

（8）設(shè)定分析選項

1）選擇菜單MainMenu：Solution→AnalysisOptions→EnergyOptions命令，彈出EnergyOptions對話框，激活HourglassEnergy，SlidingInterface和RaylieghEnergy3個單選鈕，取消StonewallEnergy單選鈕，單擊OK按鈕，關(guān)閉對話框。

2）選擇菜單MainMenu：Solution→AnalysisOptions→CPULimit命令，彈出CPULimit對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。

3）選擇菜單MainMenu：Solution→AnalysisOptions→BulkViscosity命令，彈出BulkViscosity對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。

（9）設(shè)置求解時間和時間步控制

1）選擇菜單MainMenu：Solution→TimeControls→SolutionTime命令，彈出SolutionTimeforLS-DYNAExplicit對話框，在TerminateatTime文本框中輸入100，單擊OK按鈕確認(rèn)輸入。

2）選擇菜單MainMenu：Solution→TimeControls→TimeStepCtrls命令，彈出SpectifyTimeStepScalingforLS-DYNAExplicit對話框，在Timestepscalefactor文本框中輸入0.6，單擊OK按鈕，關(guān)閉對話框。

（10）設(shè)置輸出類型和時間間隔

1）選擇菜單MainMenu：Solution→OutputControls→Output命令，彈出SpecifyOutputforLS-DYNASolver對話框，在下拉列表框中選擇Add，在Produceoutputfor...下拉列表框中選擇LS-DYNA，單擊OK按鈕，關(guān)閉對話框。

2）選擇菜單MainMenu：Solution→OutputControls→Freq→TimeStepSize命令，彈出SpecifyFrequency對話框，使用默認(rèn)值，單擊OK按鈕，關(guān)閉對話框。

（11）輸出K文件

1）選擇菜單UtilityMenu：Select→Everything命令。

2）選擇菜單MainMenu：Solution→WriteJobname.k命令，彈出InputfilestobeWrittenforLS-DYNA對話框，在Writeresultsfilesfor...下拉列表框中選擇LS-DYNA，在Writeinputfilesto...文本框中輸入explosion_undergroud.k，單擊OK按鈕，程序?qū)⒃诠ぷ髂夸浵律蒭xplosion_undergroud.k文件。

（12）編輯修改explosion_undergroud.k文件

1）用記事本打開工作目錄下的explosion_undergroud.k文件。

2）對照后面所附的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；

⑤添加用于起爆點設(shè)置的*INITIAL_DETONATION關(guān)鍵字；

⑥用土壤泡沫材料模型*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)模型中的個別參數(shù)。

3）將修改后的explosion_undergroud.k文件與K文件對照，確認(rèn)無誤后存盤。

（13）求解

1）在Launcher對話框的SimulationEnvironment下拉列表框中選擇LS-DYNASolver，并在AnalysisType下面的單選框中選擇TypicalLS-DYNAAnalysis。

2）在對話框的WorkingDirectory中選擇LS-DYNA的輸入文件存放的路徑，并在KeywordInputFile欄中選定explosion_undergroud.k。

3）單擊Run按鈕，程序?qū)⒄{(diào)用LS-DYNA971求解器開始求解。

4．后處理

（1）顯示土壤鼓包運動過程

1）運行LS-PREPOST程序，選擇菜單→BinaryPlot命令，打開工作目錄下的d3plot文件，程序?qū)⒆x入結(jié)果文件。

2）選擇菜單Misc.→Reflect→ReflectaboutYZplane命令，將模型做YZ平面對稱鏡像，如圖9-12所示。

圖9-12對稱鏡像

3）選擇主菜單區(qū)域Fcomp→Misc→Volumefractionmat#2，單擊Apply按鈕，然后單擊動畫控制區(qū)域中的動畫播放按鈕，程序?qū)⒃趫D形對話框中連續(xù)動態(tài)地顯示炸藥爆炸后土壤的鼓包運動過程，圖9-13a～d給出了幾個不同時刻爆炸場的形態(tài)。

（2）輸出土壤鼓包運動的速度與位移時程曲線

1）運行LS-PREPOST，讀入d3plot結(jié)果數(shù)據(jù)文件。

圖9-13爆炸后土壤的鼓包與混凝土板運動過程

2）選擇主菜單SelPar，在PartSelection中選擇H4（只顯示Part4）。

3）選擇主單History→Nodal→Y-velocity命令，在圖形對話框中選擇H4對稱面上最下方的一個節(jié)點（節(jié)點號為89606），如圖9-14所示。單擊Plot按鈕，將彈出PlotWindow-1對話框，輸出混凝土板頂點的Y方向速度時程曲線，如圖9-15所示。

圖9-14節(jié)點89606

4）選擇菜單History→Nodal→Y-displacement命令，在圖形對話框中仍選擇節(jié)點89606，單擊Plot按鈕，將彈出PlotWindow-1對話框，輸出混凝土板頂點的Y方向位移時程曲線，如圖9-16所示。

圖9-15鼓包運動的速度時程曲線

圖9-16鼓包運動的位移時程曲線5．命令流文件

/PREP7

！進(jìn)入前處理器

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

！激活前面定義的數(shù)據(jù)

BLOCK,0,255,-250,140,0,1.5,

！創(chuàng)建模型體

wpro,,-90.000000,

！將工作平面繞X軸順時針旋轉(zhuǎn)90度

wpoff,0,0,-6

！將工作平面向Z軸負(fù)向移動6cm

VSBW,1

！用工作平面分離體

WPSTYLE,,,,,,,,0

！重新設(shè)置工作平面

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

！重新設(shè)置工作平面

wpoff,6,0,0

！將工作平面沿WX方向平移6cm

wpro,,,-90.000000

！將工作平面繞WY軸順時針旋轉(zhuǎn)90度

FLST,2,4,6,ORDE,2

FITEM,2,2

FITEM,2,-5

VSBW,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,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

！設(shè)置單元長度為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

！對土壤模型進(jìn)行網(wǎng)格劃分

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

！對空氣模型進(jìn)行網(wǎng)格劃分

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

！對混泥土板模型進(jìn)行網(wǎng)格劃分

REAL,

ESYS,0

SECNUM,

VSEL,S,,,2

ASLV,S

LSLA,S

EPLOT

！顯示模型材料的有限元網(wǎng)格

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

！設(shè)置單元長度為2

CM,_Y,VOLU

VSEL,,,,2

CM,_Y1,VOLU

CHKMSH,'VOLU'

CMSEL,S,_Y

VMESH,_Y1

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

EDPART,CREATE

！創(chuàng)建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,

！對靠近炸藥模型一側(cè)施加X方向的約束

EDENERGY,1,0,1,1

！設(shè)定分析選項

EDCPU,0,

EDBVIS,1.5,0.06,

TIME,100,

！設(shè)置求解時間

EDCTS,0,0.6,

！設(shè)置時間步

EDOPT,ADD,blank,LSDYNA

！設(shè)置輸出類型

EDRST,0,

EDHTIME,0,

！設(shè)置時間間隔

EDDUMP,0,

ALLSEL,ALL

！選擇所有元素

EDWRITE,LSDYNA,explosion_undergroud,k

！生成K文件

/EXIT

！退出ANSYS

6．K文件

*KEYWORD

*TITLE

EXPLOSION_UNDERGROUND.K

$

*DATABASE_FORMAT

0

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

NODEDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*NODE

16.000000000E+006.000000000E+000.000000000E+00

0

0

26.000000000E+006.000000000E+001.500000000E+00

0

0

36.000000000E+00-6.000000000E+000.000000000E+00

0

0

913662.480000000E+028.400000000E+011.500000000E+00

0

0

913672.480000000E+028.600000000E+011.500000000E+00

0

0

913682.480000000E+028.800000000E+011.500000000E+00

0

0

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

SECTIONDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*SECTION_SOLID_ALE

1,11

*SECTION_SOLID_ALE

2,11

*SECTION_SOLID_ALE

3,11

*SECTION_SOLID

4,0

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

ALEDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*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

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

MATERIALDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*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

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

PARTSDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$

*PART

Part

1forMat

1andElemType

1

1

1

1

1

0

0

0

$

*PART

Part

2forMat

2andElemType

2

2

2

2

0

0

0

0

$

*PART

Part

3forMat

3andElemType

3

3

3

3

3

0

0

0

$

*PART

Part

4forMat

4andElemType

4

4

4

4

0

0

0

0

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

ELEMENTDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*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

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

COORDINATESYSTEMS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

LOADDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

RIGIDBOUNDRIES

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

BOUNDARYDEFINITIONS

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

*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