finish
/clear
/prep7
k,1,
k,2,1,
k,3,1,1
k,4,,1
kgen,2,1,4,1,,,-0.05
kgen,2,1,4,1,,,0.05
v,1,2,3,4,5,6,7,8
v,9,10,11,12,1,2,3,4
!et,2,mesh200,6
et,1,45
alls
wpcsys,,0
wprota,0,90,0 !绕x轴旋转
*do,i,1,9
wpave,0,0.1*i
vsbw,all
*enddo
wpcsys,,0
wprota,0,0,90
*do,i,1,9
wpave,0.1*i
vsbw,all
*enddo
MP,EX,1,3.25E10
MP,PRXY,1,0.2
MP,DENS,1,2500
vsel,s,loc,z,-1,0
vatt,1,,1,
MP,EX,2,0.1502E10
MP,PRXY,2,0.2
MP,DENS,2,2400
vsel,s,loc,z,0,1
vatt,2,,1,
alls
mshape,0,3D
mshkey,1
lesize,all,1
vmesh,all
!!nsel,s,loc,x,-0.1,0.01
!!nsel,a,loc,x,0.91,1.01
!!nsel,u,loc,y,0.01,0.99
!!nsel,r,loc,z,-0.01,0.01
!!d,all,all
nsel,s,loc,y,0
d,all,all
finish
allsel
/solu
antype,2
modopt,lanb,10,,,,1
mxpand,10,,,,yes
solve
finish
/solu
expass,on
mxpand,10
solve
finish
save
/post1
set,list
这是实体单元的模型,我对比了下 |