隧道爆破LSDYNAAPDL命令流.docx
《隧道爆破LSDYNAAPDL命令流.docx》由会员分享，可在线阅读，更多相关《隧道爆破LSDYNAAPDL命令流.docx（11页珍藏版）》请在冰豆网上搜索。
隧道爆破LSDYNAAPDL命令流
!
（1）工作环境设置,进入前处理器
fini
/cle
/prep7
!
（2）定义单元类型
ET,1,SHELL163
ET,2,SOLID164
!
（3）定义材料模型
MP,DENS,1,2400
MP,EX,1,1.9e10
MP,NUXY,1,0.29
MP,DENS,2,1150
MP,EX,2,6.72E6
MP,NUXY,2,0.33
MP,DENS,3,2600
MP,EX,3,2.0E10
MP,NUXY,3,0.2
MP,DENS,4,1.2e-3
MP,EX,4,
MP,NUXY,4,
TB,EOS,4,,,2,1
!
（5）建立模型
K,1,0,1.98,0
K,2,3.43,0,0
K,3,-3.43,0,0
K,4,0,4.13,0
CIRCLE,1,3.39
CIRCLE,2,7.35
CIRCLE,3,7.35
CIRCLE,4,7.4
lcsl,all
lsel,s,line,,33,38,1
lsel,a,line,,28,29,1
lsel,a,line,,31,32,1
lsel,a,line,,47,48,1
lsel,a,line,,11
lsel,a,line,,21,23,1
lsel,a,line,,43
lsel,a,line,,45
lsel,a,line,,26
lsel,a,line,,8
lsel,a,line,,41,43,1
lsel,a,line,,3,4,1
lsel,a,line,,17
lsel,a,line,,20
ldel,all,,,1
!
创建关键点
K,101,-45,-45,0,
K,102,-9,-45,0,
K,103,9,-45,0,
K,104,45,-45,0,
K,105,45,-8,0,
K,106,9,-8,0,
K,107,-9,-8,0,
K,108,-45,-8,0,
K,109,-45,8,0,
K,110,-9,8,0,
K,111,9,8,0,
K,112,45,8,0,
K,113,45,45,0,
K,114,9,45,0,
K,115,-9,45,0,
K,116,-45,45,0,
LSTR,101,102
LSTR,102,103
LSTR,103,104
LSTR,104,105
LSTR,105,106
LSTR,106,107
LSTR,107,108
LSTR,108,109
LSTR,109,110
LSTR,110,111
LSTR,111,112
LSTR,112,113
LSTR,113,114
LSTR,114,115
LSTR,115,116
LSTR,101,108
LSTR,109,116
LSTR,102,107
LSTR,107,110
LSTR,110,115
LSTR,103,106
LSTR,106,111
LSTR,111,114
LSTR,105,112
allsel,all
a,101,104,113,116
asbl,all,all
adele,10
numcmp,area
lgen,2,1,17,1,,,-60,,,,
lgen,2,20,23,1,,,-60,,,,
lgen,2,26,,,,,-60,,,,
lgen,2,28,29,1,,,-60,,,,
lgen,2,1,17,1,,,-42,,,,
lgen,2,20,23,1,,,-42,,,,
lgen,2,26,,,,,-42,,,,
lgen,2,28,29,1,,,-42,,,,
LCOMB,19,44,0
LCOMB,25,39,0
LCOMB,18,46,0
LCOMB,30,40,0
NUMMRG,ALL
k,500,0,0,-42
circle,500,0.03
lsel,s,loc,z,-59,-61
a,5,9,25,29
asbl,all,all
allsel,all
lsel,s,loc,z,-41,-43
a,43,46,55,58
asbl,all,all
allsel,all
!
拖拉成体
asel,s,loc,z,-41,-43
vext,all,,,,,2
asel,s,loc,z,-59,-61
vext,all,,,,,18
asel,s,loc,z,-1,1
vext,all,,,,,-40
allsel,all
wpoff,0,0,-40.6
vsbw,1
wpoff,0,0,-0.8
vsbw,30
WPSTYLE,,,,,,,,0
vplot
!
粘结全体
vglue,all
!
划分网格
lsel,s,line,,83,86,1,
lsel,a,line,,225,229,2,
lsel,a,line,,230,
lsel,a,line,,89,
lsel,a,line,,237,238,1,
lsel,a,line,,235,
lsel,a,line,,39,40,1,
lsel,a,line,,44,46,2,
lesize,all,,,8,
lsel,all
lsel,s,line,,224,228,2,
lsel,a,line,,223,
lsel,a,line,,87,88,1,
lsel,a,line,,90,
lsel,a,line,,236
lesize,all,,,24,
lsel,all
lsel,s,line,,239,242,1,
lesize,all,,,32,
lsel,all
vsel,s,volume,,29
vatt,1,,2,0,
vsweep,all
vsel,s,volume,,31
vatt,2,,2,0,
vsweep,all
vsel,s,volume,,1
vatt,4,,2,0,
vsweep,all
vsel,all
vsel,s,volume,,30,36,3
vsel,a,volume,,38
vatt,1,,2,0,
esize,2.5,
mshape,0,3d
mshkey,1,
vmesh,all
lsel,s,line,,2,14,6
lsel,a,line,,29
lesize,all,,,18,
vsel,s,volume,,32,34,2
vsel,a,volume,,35,37,2
vatt,1,,2,0,
vsweep,all
allsel,all
AADD,20,29
!
中间网格划分
lsel,s,line,,95
lsel,a,line,,107
lsel,a,line,,123
lsel,a,line,,130
lesize,all,,,4
vsel,s,volume,,2,4,2
vsel,a,volume,,7,9,2
vatt,1,,2,0,
vsweep,all
allsel,all
vsel,s,volume,,3,5,2
vsel,a,volume,,6,8,2
vatt,1,,2,0,
vsweep,all
allsel,all
vsel,s,volume,,39,47,1
vatt,1,,2,0,
vsweep,all
allsel,all
lsel,s,line,,18,19,1
lsel,a,line,,24,27,3
lsel,a,line,,271,274,3
lsel,a,line,,276,279,3
lesize,all,1
vsel,s,volume,,48
vatt,1,,2,0,
vsweep,all
allsel,all
vsel,s,volume,,49
vatt,1,,2,0,
smrtsize,6
MSHAPE,1,3D
mshkey,0
vmesh,all
allsel,all
!
（4）定义实参数
*SET,_RC_SET,1,
R,1
RMODIF,1,1,5/6,,0.1,,,,
asel,s,area,,135,138,1
asel,a,area,,197,200,1
aatt,3,1,1,0,
amesh,all
allsel,all
!
建立part
edpart,create
!
施加约束
nsel,s,loc,y,-44.9,-45.1
d,all,uy,0.0
d,all,ux,0.0
d,all,uz,0.0
allsel,all
!
施加无反射边界条件
nsel,s,loc,y,-44.9,-45.1
nsel,a,loc,y,44.9,45.1
nsel,a,loc,z,-0.1,0.1
nsel,a,loc,z,-59.9,-60.1
nsel,a,loc,x,-44.9,-45.1
nsel,a,loc,x,44.9,45.1
cm,nonr,node
ednb,add,nonr,1,1
allsel,all
!
退出前处理器
finish
!
进入求解器
/solu
!
ALE算法控制选项
edgcale,1,1,1
edale,add,all,-1,0,0,0,0,0,1e20,1
!
设置能量控制选项
edenergy,1,1,1,1
!
设置求解时间
time,0.08
!
设置结果文件的类型和输出步数
edopt,add,blank,sdyna
edrst,1600
edhtime,1,
升级会员 