贴出来大家看的清楚点。
方法1:利用mass21单元刚性区域法
/prep7
block,0,1,0,1,0,2
k,9,0.5,0.5,2.5
mp,ex,1,2e10
mp,prxy,1,0.2
mp,prxy,1,0.3
r,2,1e-6
et,1,45
et,2,21
keyopt,2,3,0
lesize,all,0.2
vmesh,all
ksel,s,,,9
type,2
real,2
kmesh,all
allsel
nsel,s,loc,z,2,3
NPLOT
CERIG,node(0.5,0.5,2.5),ALL,ALL, , , ,
allsel
/SOLU
f,node(0.5,0.5,2.5),my,100e3
FINISH
/SOL
nsel,s,loc,z,0
d,all,all
allsel
solve
方法2:利用mpc184单元
/prep7
block,0,1,0,1,0,2
mp,ex,1,2e10
mp,prxy,1,0.2
mp,prxy,1,0.3
et,1,45
et,2,184
keyopt,2,1,1
lesize,all,0.2
vmesh,all
n,1000,0.5,0.5,2.5
type,2
mat,2
*do,i,1,36
e,1000,36+i
*enddo
allsel
allsel
/SOLU
f,node(0.5,0.5,2.5),my,100e3
FINISH
/SOL
nsel,s,loc,z,0
d,all,all
allsel
solve
方法3:使用rbe3命令
/prep7
block,0,1,0,1,0,2
k,9,0.5,0.5,2.5
mp,ex,1,2e10
mp,prxy,1,0.2
mp,prxy,1,0.3
r,2,1e-6
et,1,45
et,2,21
keyopt,2,3,0
lesize,all,0.2
vmesh,all
ksel,s,,,9
type,2
real,2
kmesh,all
allsel
*dim,sla,array,36
*do,i,1,36
sla(i)=i+36
*enddo
*dim,sla2,array,36
*do,i,1,36
sla2(i)=i+36
*enddo
allsel
rbe3,node(0.5,0.5,2.5),all,sla,sla2
allsel
/SOLU
f,node(0.5,0.5,2.5),my,100e3
FINISH
/SOL
nsel,s,loc,z,0
d,all,all
allsel
solve |