问题如标题：关于一个dyna中，定义节点位移为0，但是爆破后还是移动了，为什么
命令流和修改部分K文件如下
建模
finish
/clear
!单位统一cm-g-us
!(1)设置工程属性
/filname,smooth blast
/prep7
/view,,1,1,1 !设置图方式形显示
!(2)定义变量、单元类型、材料与实常数
L1=200 !炮孔长度
D1=4 !炮孔直径
R1=D1/2 !炮孔半径
N1=3 !炮孔个数
L2=50 !堵塞长度
L3=20 !单个炸药长度
D3=2 !炸药直径
R3=D3/2 !炸药半径
N3=5 !单孔炸药个数
l4=80 !炮孔孔距
l5=30 !炮孔距建基面的距离
l6=80 !炮孔预留距离
ei=10 !网格划分
l7=50
!******************************
ii=2 !影响倍数
x1=ii*l4 !l4炮孔孔距
x2=(ii+1)*l4
x3=(ii+2)*l4
x4=(2*ii+2)*l4
y1=x1
y2=l2+y1 !l2堵塞长度
y3=l1+y2 !l1炮孔长度
y4=y3+y1
z1=x1+l4
z2=z1+l5
z3=z2+l6
z4=z3+l7
!**********************************
*do,i,1,3
ET,i,SOLID164 !炸药参数
R,i
MP,EX,i,1.85
MP,NUXY,i,0
MP,DENS,i,1.85
*enddo
!***************************
ET,16,SOLID164 !炮孔内空气单元
R,16
MP,DENS,16,0.0012
MP,EX,16,0
MP,NUXY,16,0
tb,eos,16,,,2,1
!**********************************
et,17,solid164 !被爆结构单元
r,17 !
mp,ex,17,2.5e-2 !弹性模量
mp,dens,17,2.35 !密度
mp,nuxy,17,0.25 !泊松比
!****************************
!****************************
!岩石构造
!(3)实体模型建立
k,1,0,0,0 !最大体
k,2,0,y4,0
k,3,0,y4,z3
k,4,0,y1,z3
k,5,0,y1,z1
k,6,0,0,z1
a,1,2,3,4,5,6
vext,1,,,x4
*do,i,1,N1 !N1 !定义关键点
k,1000+3*i,x1-l4+i*l4,y3,z2
*enddo
*do,i,1,N1 !N1 !定义关键点
k,2000+3*i,x1-l4+i*l4,y4,z2
*enddo
wprota,,90
!矩形体
wprota,45
*do,i,1,n1
kwpave,2000+3*i
rprism,,(y4-y1),4,2*ei
*enddo
wprota,-45
!孔洞
*do,i,1,n1
kwpave,1000+3*i
cyl4,,,0,0,r1,360,l1
*enddo
VSEL,S,,,2,7
VPTN,ALL
numcmp,volu
VSEL,ALL
VPTN,ALL
numcmp,volu
VDELE,1,3,1
numcmp,volu
!孔洞2
*do,i,1,n1
kwpave,1000+3*i
cyl4,,,0,0,r1,360,-y1
wpoffs,,,l1
cyl4,,,0,0,r1,360,l2
*enddo
*do,i,1,n1
kwpave,1000+3*i
cyl4,,,0,0,r3,360,-y1
wpoffs,,,l1
cyl4,,,0,0,r3,360,l2
*enddo
VSEL,ALL
vptn,all
numcmp,volu
VGLUE,ALL
!分割矩形体
vsel,s,,,1
VSEL,INVE
WPROTA,,-90
kwpave,1000+3
vsbw,all
*do,i,1,n1
kwpave,1000+3*i
*do,j,1,3
wprota,,,45
vsbw,all
*enddo
wprota,,,-45*3
*enddo
!分割整体
vsel,all
kwpave,5
vsbw,all
wprota,,90
vsbw,all
kwpave,1
wprota,,-90
*do,i,1,N1 !N1
k,1000+3*i,x1-l4+i*l4,y3,z2
*enddo
kwpave,1003
vsbw,all
wpoffs,,,ei
vsbw,all
wpoffs,,,-2*ei
vsbw,all
wprota,,,90
*do,i,1,n1
kwpave,1000+3*i
vsbw,all
wpoffs,,,ei
vsbw,all
wpoffs,,,-2*ei
vsbw,all
*enddo
kwpave,1
wprota,,,-90
VSEL,ALL
wprota,,90
kwpave,1003
vsbw,all
*do,i,1,n3*2
wpoffs,,,l3
vsbw,all
*enddo
wprota,,-90
numcmp,volu
vsel,all
CM,guodu1,volu
vatt,17,17,17
!****************************
vsel,all
vsel,inve
*do,i,1,N1 !N1 !炸药
wprota,,90
kwpave,1000+3*i
vsel,all
vsel,inve
*do,j,1,n3
wpoffs,,,l3
cyl4,,,0,0,r3,360,l3
wpoffs,,,l3
*enddo
wprota,,-90
vsbw,all
*do,K,1,3
wprota,,,45
vsbw,all
*enddo
wprota,,,-45*3
wprota,,90
numcmp,volu
f=i
vatt,f,f,f
wprota,,-90
*enddo
numcmp,volu
cmsel,s,guodu1,volu
vsel,inve
cm,zhayao,volu
!****************************
vsel,all
vsel,inve
*do,i,1,N1 !N1 !炮孔空气
wprota,,90
kwpave,1000+3*i
*do,j,1,n3
cyl4,,,0,0,r3,360,l3
cyl4,,,r3,0,r1,360,l3
wpoffs,,,l3
cyl4,,,r3,0,r1,360,l3
wpoffs,,,l3
*enddo
wprota,,-90
*enddo
numcmp,volu
vglue,all
!****************************
kwpave,1003 !炮孔分割
vsbw,all
*do,i,1,n1
kwpave,1000+3*i
*do,j,1,3
wprota,,,45
vsbw,all
*enddo
wprota,,,-45*3
*enddo
numcmp,volu !169-528
cm,paokong,volu
vatt,16,16,16
vsel,all
!****************************
!****************************
!网格划分
kwpave,1
vsel,s,loc,x,x1-ei,x1+ei
vsel,a,loc,x,x2-ei,x2+ei
vsel,a,loc,x,x3-ei,x3+ei
vsel,r,loc,z,z2-ei,z2+ei
esize,,1
aslv,s
lsla,s
lsel,r,tan1,x,0
lsel,r,tan1,z,0
lesize,all,ei
aslv,s
lsla,s
lsel,r,radius,,r1
lesize,all,,,6
aslv,s
lsla,s
lsel,r,radius,,r3
lesize,all,,,6
vmesh,all
vsel,inve
esize,ei
vmesh,all
!****************************
!建立Part
allsel,all
edpart,create
!****************************
!定义接触面
!****************************
!施加约束
nsel,s,loc,z,0-1e-3,0+1e-3
d,all,all,0
nsel,s,loc,y,0-1e-3,0+1e-3
nsel,a,loc,y,y4-1e-3,y4+1e-3
d,all,uy,0
nsel,s,loc,x,0-1e-3,0+1e-3
nsel,a,loc,x,x4-1e-3,x4+1e-3
d,all,ux,0
asel,s,loc,x,x4-1e-3,x4+1e-3
asel,a,loc,x,0-1e-3,0+1e-3
cm,ffsbj,node
ednb,add,ffsbj,0,0
allssel,all
!****************************
finish
/solu
!生成K文件
edgcale,1,1,1,1
edale,add,all,-1,0,0,0,0,0,1e20,1
time,1500 !设置分析时间
edrst,100 !文件输出步数
edhtime,1 !时间历程文件输出步数
edenergy,1,1,1,1 !能量选项
edopt,add,blank,taurus !设置结果文件类型
edout,glstat
edstart,0,0,7,d3dump
edopt,add,blank,both
edwrite,both,'ex001','k',''
finish
/exit
K文件
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ SECTION DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID_ale
1 11
*SECTION_SOLID_ale
2 11
*SECTION_SOLID_ale
3 11
*SECTION_SOLID_ale
4 11
*SECTION_SOLID_ale
5 11
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ MATERIAL DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_high_explosive_burn
1,1.2,0.55,0.15,0,
*eos_jwl
1,7.41,0.18,5.56,1.65,0.35,0.036,1
*initial_detonation
*MAT_high_explosive_burn
2,1.2,0.55,0.15,0,
*eos_jwl
2,7.41,0.18,5.56,1.65,0.35,0.036,1
*MAT_high_explosive_burn
3,1.2,0.55,0.15,0,
*eos_jwl
3,7.41,0.18,5.56,1.65,0.35,0.036,1
*MAT_NULL
16,0.120E-02,0,0,0,0,0,0
*EOS_LINEAR_POLYNOMIAL
16,0,0,0,0.4,0.4,0,0
0.00,0.00
*MAT_soil_and_foam
17,2.5,1.6e-4,2.5e-4,3.3e-11,0,0,0
0,0
0,0.05,0.09,0.11,0.15,0.19,0.21,0.22
0.25,0.30
0,4.5e-6,5e-6,6.7e-6,1.3e-5,2.1e-5,2.7e-5,3.9e-5
5.5e-5,1.2e-4
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ PARTS DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part 1 for Mat 17 and Elem Type 17
1 1 17 0 0 0 0
$
*PART
Part 2 for Mat 1 and Elem Type 1
2 2 1 1 0 0 0
$
*PART
Part 3 for Mat 2 and Elem Type 2
3 3 2 2 0 0 0
$
*PART
Part 4 for Mat 3 and Elem Type 3
4 4 3 3 0 0 0
$
*PART
Part 5 for Mat 16 and Elem Type 16
5 5 16 16 0 0 0
$
$
$ CONTROL OPTIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_ENERGY
2 2 2 2
*CONTROL_SHELL
20.0 1 -1 1 2 2 1
*CONTROL_ALE
2 1 2 -1.00 0.00 0.00 0.00 0.00
0.00 0.100E+21 1.00 1.00 0.00 0
*ale_multi-material_group
1,1
2,1
3,1
4,1
5,1
*CONTROL_TIMESTEP
0.0000 1.0000 0 0.00 0.00
*CONTROL_TERMINATION
2000,0,0,0,0
*DATABASE_HISTORY_NODE
1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ TIME HISTORY $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_GLSTAT
2000
*DATABASE_BINARY_D3PLOT
10
*DATABASE_BINARY_D3THDT
2000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ DATABASE OPTIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
0 0 3 1 0 0 0 0
0 0 4 0 0 0
*END
