首页/文章/ 详情

DYNA中爆炸算法

6月前浏览391

一、 load_blast 加载(优点:与实验较符合,且计算方便,缺点:只能是球形炸药,且不考虑流固耦合)

建好模型之后选择需要加载的面的壳单元,定义两条空载荷曲线,并加载在选择的单元上

进入K文件,修改加载的曲线为:

*DEFINE_CURVE

1

*DEFINE_CURVE

2

再添加关键字*load_blast

*LOAD_BLAST

20,0.00,0.00,58.00,0.00,4,2

0,0,0,0

4025b8a0-8a6c-417c-95e4-67fd8e4278be.png

a4770883-645b-401f-a182-77d5b0e000ba.png

4e2c7795-b9fc-42b6-8250-756f321e4878.png

同时需要修改

*LOAD_SHELL_SET

1,-2,1.00,0.00

718b85af-02a4-4189-a844-236ba3a7dd95.png

31574874-61a2-4b02-b730-bc4b5c991cc8.png


二、ale方法

(1)、*ALE_MULTI-MATERIAL_GROUP

6(空气part号),1

7(炸药part号),1

*INITIAL_VOLUME_FRACTION_GEOMETRY

1,1,1

6(炸药几何形状,球形),0,7

0,0,0(爆点坐标),1.64(炸药半径)

*CONSTRAINED_lagrange_in_SOLID

1(结构part号),6(空气part号),1,1,3,4,2,0

0,1e10,0.1,0,0.3,0


*INITIAL_DETONATION

7(炸药part号),0,0,0,0.0




(2)、*SECTION_SOLID_ale

         4         11

此行空白




(3)、材料参数(g-cm-um)

*MAT_NULL

         6 0.120E-02  0.00      0.00      0.00      0.00      0.00      0.00    

*EOS_LINEAR_POLYNOMIAL

         6  0.00      0.00      0.00      0.00     0.400     0.400      0.00    

  0.00      0.00   

*MAT_HIGH_EXPLOSIVE_BURN

7(炸药材料号),1.63,0.693,0.21,0.0,0,0,0

*EOS_JWL

7(状态方程号),3.7400002,0.0737,4.1500001,0.95,0.3,0.043,1.00



(4)、 *CONTROL_BULK_VISCOSITY

1.5,0.06,1



免责声明:

本页面/内容部分素材来源于互联网公 开 信 息,旨在传递更多信息,不代表本平台立场。

版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。

本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。

ACTUM爆炸材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-03-27
最近编辑:6月前
仿真支持爱好者
在仿真的路上越走越远
获赞 333粉丝 26文章 318课程 0
点赞
收藏
作者推荐

ABAQUS地震反应谱分析

ABAQUS反应谱法计算地震反应的简单实例问题描述:(木例的问题引用《有限元法及其应用》一书中陆新征博士ANSYS算例的问题)悬臂柱高12m,工字型截面(图1),密度7800kg/m3,EX=2.1ellPa,泊松比0.3,所有振型的阻尼比为2%,在3m高处行-集中质量160kg,在6m、9m、12m处分别为120kg的集中质量。反应谱按7度多遇地震,取地震影响系数为0.08,第一组,III类场地,卓越周期Tg=0.45s,几点说明:本例建模过程使用CAE:添加反应谱必须在inp中加关键词实现,CAE不支持反应谱;*Spectrum不可以在keywordeditor中添力口,keywordeditor不支持此关键词读入:ABAQUS的反应谱法计算过程以及后处理耍比ANSYS方便的多。操作过程为:(1)打开ABAQUS/CAE,点击createmodeldatabase0(2)进入Part模块,点击createpart,命名为column,3D、deformationwire。OK(3)Createlines:connected,分别输入0,0:0,3;0,6:0»9:0,12cOK。退出sketch。(4)进入property模块,creatematerial,name:steel»general-»density»massdensity:7800,mechanical—»elasticity—»elastic,young'smodulus:2.1el1»poisson'sratio:0.3.OK(5)Createsection,name:I,category:beam,type:beam.Continue,createprofile,name:I,shape:I,按图1尺寸输入界面尺寸,ok。在profilename选择I,materialname选择steelcOk0(6)Assignsection,选择全部,done,弹出的对话框选择section:I,Oh(7)Assignbeamorientation,选择全部,默认值,OK.(8)View—»partdisplayoptions,在弹出的对话框里勾选,renderbeamprofiles,以可视化梁截面形状。(9)添力口集中质量,Special—»inertia—»create,name:massl,type:pointmass/inertia>continue*选择0,3位置点,done»mass:160,ok。create,name:mass2,type:pointmass/inertia«continue*选择0,6;0»9:0,12位置点(按shift多选),done»mass:120»Ok,dismiss0(10)Assembly—»instancepart,instancetype选dependent(meshonpart),Ok.(11)Step—»createstep,name:step-1»proceduretype选freqency>continue,在basic选项卡中,eigensolver选择频率提取方法,木例选用lanczos法,numberofeigenvaluesrequest,选value,输入10。OkeMcreatestep>createstep,name:step-2,proceduretype选responsespectrum,continue*在basic选项卡中,excitations选择单向singledirection,sumations选择squarerootofthesumofsquares(SRSS)法,useresponsespectrum:sp(反应谱的name,后面再inp中添加),方向余弦(0,0,1),scalefactor:1.进入damping选项卡,阻尼使用吏接模态(directmodal),勾选directdampingdata,startmode:1,endmode:8,criticaldampingfraction:0.02。Ok(12)进入load模块,Load—»createboundarycondition,name:fixed,step选择initial*categorymechanical,typesdisplacement/rotation,continue,选择0,0点,done»勾选ul~ur3所有6个自由度。Ok。(13)进入mesh模块,object选择part,点seededgebynumber,选择所有杆,done*输入3,done,点assignelementtype,选择全部杆,done,默认B31,ok。点meshpart,yes。(14)进入job模块,name:demo-spc»source:model,continue,默认,0匕进入jobmanager,点击writeinput,在T.作目录生成demo-spc.inp文件。(15)进入ABAQUS工作目录,使用UltraEdit软件(或其他类似软件)打开demo-spc.inp,*Boundary关键词的后面加如下根据问题叙述确定的反应谱:*Spectrum,type=acceleration,name=sp0.1543,0.167,00.1915,0.25,00.2102,0.333,00.2241,0.444,00.25,0.5,00.3295,0.667,00.4843,1,00.5987,1.25,00.7868,1.667,01.0342,2.222,01.0342,10,00.3528,10000,0第一列为加速度,第二列为频率(图2)保存。(16)进入job模块,createjob.name:spc,source选择inputfile,inputfileselect:工作目录下的demo-spc.inp,continues,默认,ok,UAjobmanager,选择spc,submit,计算成功!(17)点击results进入后处理模块,可以看到最大位移为3.159cm,这与陆博士讲解的ANSYS结果3.1611cm基木一致。可以在看工.作目录下的spc.dat文件得到详细的频率和模态分析结果(图34)免责声明:本页面/内容部分素材来源于互联网公开信息,旨在传递更多信息,不代表本平台立场。版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈