LS-DYNA动力学分析--炸药在土中爆炸GUI及APDL例子要点.docx

第九章 LS-DYNA动力学分析第四节 LS-DYNA范例解析炸药在土壤内部爆炸分析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模型尺寸图/mm3在File Management对话框的Working Directory文本框中输入E:explosion_undergro-und作为工作目录假设工作目录为E盘，在Job Name文本框中输入explosion_underground作为工作文件名，其他选项用默认值。单击Run按钮，运行ANSYS程序，进入ANSYS的操作界面。2选择单元类型1选择菜单Main Menu：PreprocessorElement TypeAdd/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：PreprocessorMaterial PropsMaterial Models命令，弹出Define Material Model Behavior对话框。2选择菜单LS-DYNAEquation of StateGruneisenNull命令，弹出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对话框中单击MaterialNew Model命令，弹出Define Material ID对话框，在文本框中输入2，单击OK按钮关闭对话框。4选择菜单LS-DYNANonlinearInelasticKinematic HardeningPlastic Kinematic命令，弹出Plastic Kinematic Properties for Material Number 2对话框，在文本框中输入材料本构参数由于ANSYS前处理所带的材料模型中没有我们所需要的土壤模型，暂用随动硬化塑性材料模型代替，后面将在K文件中对材料模型参数进行修改，单击OK按钮关闭对话框。5在Define Material Model Behavior对话框中单击MaterialNew Model命令，弹出Define Material ID对话框，在文本框中输入3，单击OK按钮关闭对话框。6选择菜单LS-DYNAEquation of StaticGruneisenNull命令，弹出Null Properties for Material Number 3对话框。在文本框中，输入以下数据：DENS=0.001252；C=0.03437； GAMAO=1.4。其余选项采用默认值，输入完成后单击OK按钮关闭对话框。7在Define Material Model Behavior对话框中单击MaterialNew Model命令，出现Define Material ID对话框，在文本框中输入4，单击OK按钮关闭对话框。8选择菜单LS-DYNANonlinearInelasticKinematic HardeningPlastic 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对话框中选择MaterialExit命令，关闭对话框。此时共定义了4种材料。4创立几何模型、划分网格1选择菜单Utility Menu:WorkPlaneChange Active CS toWorking Plane命令。2选择菜单Main Menu：PreprocessorModelingCreateVolumesBlockBy Dimensions命令，弹出Create Block by Dimensions对话框，输入模型体包括炸药、土壤和空气域的坐标，如图9-3所示。图9-3 Create Block by Dimensions对话框3选择菜单Utility Menu:WorkPlaneOffset WP by Increments命令，弹出Offset WP面板，将按钮 下方的角度设置滑动条移到最右，滑动条上方的数值变为90，单击按钮 ，将视窗中的工作平面坐标绕X轴旋转90度。4在Offset WP面板中的X，Y，Z Offsets文本框中输入0，0，-6，单击Apply按钮，工作平面向Z轴负向移动6cm。5选择菜单Main Menu：PreprocessorModelingOperateBooleansDivideVolu 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：PreprocessorModelingOperateBooleansDivideVolu by WrkPlane命令，在弹出的Divide Vol by Wr拾取菜单面板中单击Pick All按钮，将视图中的4个体分割为8个体。8选择菜单Utility Menu:WorkPlaneAlign WP withGlobal Cartesian命令，使工作平面与全局坐标系重合。9选择菜单Main Menu：PreprocessorModelingCreateVolumesBlockBy Dimensions命令，在弹出Create Block by Dimensions文本框中输入混凝土板模型体的坐标，如图9-4所示。图9-4混凝土板的模型尺寸10创立完成后的实体全模型如图9-5所示。图9-5 实体全模型示意图5划分网格1选择菜单Main Menu：PreprocessorMeshingMeshTool命令，弹出MeshTool对话框，单击Elements Attributes选择栏下侧的Set按钮，弹出Meshing Attributes对话框，在TYPE下拉列表框中选择1 SOLID164，在MAT下拉列表框中选择1，其他栏用默认值，单击OK按钮，确认选择，关闭对话框。2选择菜单Utility Menu:SelectEntities.命令，弹出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:SelectEverything命令，然后选择Utility Menu:PlotElements命令，视图中显示所有材料模型的有限元网格由于混凝土板的网格与空气网格叠加在一起，看上去好似只有一种颜色，如图9-9所示。图9-9模型网格局部6创立PART选择菜单Main Menu：PreprocessorLS-DYNA OptionsParts 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：PreprocessorLS-DYNA OptionsConstraintsApplyOn 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：SolutionAnalysis OptionsEnergy Options命令，弹出Energy Options对话框，激活Hourglass Energy，Sliding Interface和Rayliegh Energy 3个单项选择钮，取消Stonewall Energy单项选择钮，单击OK按钮，关闭对话框。2选择菜单Main Menu：SolutionAnalysis OptionsCPU Limit命令，弹出CPU Limit对话框，使用默认值，单击OK按钮，关闭对话框。3选择菜单Main Menu：SolutionAnalysis OptionsBulk Viscosity命令，弹出Bulk Viscosity对话框，使用默认值，单击OK按钮，关闭对话框。9设置求解时间和时间步控制1选择菜单Main Menu：SolutionTime ControlsSolution Time命令，弹出Solution Time for LS-DYNA Explicit对话框，在Terminate at Time文本框中输入100，单击OK按钮确认输入。2选择菜单Main Menu：SolutionTime ControlsTime Step Ctrls命令，弹出Spectify Time Step Scaling for LS-DYNA Explicit对话框，在Time step scale factor文本框中输入0.6，单击OK按钮，关闭对话框。10设置输出类型和时间间隔1选择菜单Main Menu：SolutionOutput ControlsOutput File Types命令，弹出Specify Output File Types for LS-DYNA Solver对话框，在File Options下拉列表框中选择Add，在Produce output for.下拉列表框中选择LS-DYNA，单击OK按钮，关闭对话框。2选择菜单Main Menu：SolutionOutput ControlsFile Output FreqTime Step Size命令，弹出Specify File Output Frequency对话框，使用默认值，单击OK按钮，关闭对话框。11输出K文件1选择菜单Utility Menu：SelectEverything命令。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按钮，程序将在工作目录下生成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：OpenBinary Plot命令，翻开工作目录下的d3plot文件，程序将读入结果文件。2选择菜单Misc.ReflectReflect about YZ plane命令，将模型做YZ平面对称镜像，如图9-12所示。图9-12对称镜像3选择主菜单区域FcompMiscVolume fraction mat#2，单击Apply按钮，然后单击动画控制区域中的动画播放按钮，程序将在图形对话框中连续动态地显示炸药爆炸后土壤的鼓包运动过程，图9-13ad给出了几个不同时刻爆炸场的形态。2输出土壤鼓包运动的速度及位移时程曲线1运行LS-PREPOST，读入d3plot结果数据文件。图9-13爆炸后土壤的鼓包及混凝土板运动过程2选择主菜单SelPar，在Part Selection中选择H4只显示Part4。3选择主单HistoryNodalY-velocity命令，在图形对话框中选择H4对称面上最下方的一个节点节点号为89606，如图9-14所示。单击Plot按钮，将弹出Plot Window-1对话框，输出混凝土板顶点的Y方向速度时程曲线，如图9-15所示。图9-14节点896064选择菜单HistoryNodalY-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.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！激活前面定义的数据BLOCK,0,255,-250,140,0,1.5,！创立模型体wpro,-90.000000,！将工作平面绕X轴顺时针旋转90度wpoff,0,0,-6！将工作平面向Z轴负向移动6cmVSBW,1！用工作平面别离体WPSTYLE,0！重新设置工作平面wpoff,0,0,12！将工作平面向Z轴正向移动6cmFLST,2,2,6,ORDE,2FITEM,2,2FITEM,2,-3VSBW,P51X！用工作平面别离体wpoff,0,0,72！将工作平面向Z轴正向移动78cmFLST,2,3,6,ORDE,3FITEM,2,1FITEM,2,-2FITEM,2,4VSBW,P51X！用工作平面别离体WPSTYLE,0！重新设置工作平面wpoff,6,0,0！将工作平面沿WX方向平移6cmwpro,-90.000000！将工作平面绕WY轴顺时针旋转90度FLST,2,4,6,ORDE,2FITEM,2,2FITEM,2,-5VSBW,P51X！用工作平面别离体WPSTYLE,0！重新设置工作平面WPCSYS,-1,0！将当前坐标系X-Y平面定义为工作平面BLOCK,0,250,78,90,0,1.5,！创立混凝土板MAT,1！对炸药模型进行网格划分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,-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！设置单元长度为1.5MSHAPE,0,3DMSHKEY,1CM,_Y,VOLUVSEL,9CM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2MAT,2！对土壤模型进行网格划分REAL,ESYS,0SECNUM,FLST,5,5,6,ORDE,4FITEM,5,1FITEM,5,6FITEM,5,10FITEM,5,-12CM,_Y,VOLUVSEL,P51XCM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2MAT,3！对空气模型进行网格划分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,4！对混泥土板模型进行网格划分REAL,ESYS,0SECNUM,VSEL,S,2ASLV,SLSLA,SEPLOT！显示模型材料的有限元网格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！设置单元长度为2CM,_Y,VOLUVSEL,2CM,_Y1,VOLUCHKMSH,VOLUCMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2EDPART,CREATE！创立PARTALLSEL,ALL！选择所有实体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,11FITEM,2,15FITEM,2,19FITEM,2,25DA,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！退出ANSYS6K文件*KEYWORD*TITLEEXPLOSION_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_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,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,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,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 Mat 3 and Elem Type 3 3 3 3 3 0 0 0$*PARTPart 4 for Mat 4 and Elem Type 4 4 4 4 0 0 0 0$ ELEMENT DEFINITIONS $