刚性区域相关

flypanzxl

做了一个筒子，想在筒子的上表面设置刚性区域，自由度ROTY，结果提示说不能，用面单元就可以，实体单元难道不能这么做？因为实体单元只有3个平动自由度？如要要实现，应该怎么做呢？
下附命令流：
fini
/clear
r1=880
r2=860
h1=2300
r3=988.5
r4=736.5
h2=60
r5=890
h3=20
h4=820
b=620
r6=100
h5=320
m=5.0288e9
fc=3.5934e5
fs=3.3874e4
/prep7
et,1,95
et,2,21,,,2
mp,ex,1,2.1e5
mp,prxy,1,0.26
/view,1,1,1,1
/ang,1,-120,zs,1
/replot
cyl4,,,r5,,,,h3
cyl4,,,r1,,r2,,h1+h2+h3
wpoffs,,,h3+h1
cyl4,,,r3,,r4,,h2
vsbv,1,2,,,keep
vsbv,3,2,,,keep
wpcsys,-1
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
wpcsys,-1
wpoffs,,,h3
vsbw,all
wpoffs,,,h1
vsbw,all
alls
wpcsys,-1
mshkey,1
mshape,0,3d
esize,100
vmesh,all
type,2
kmax=kpinqr(0,14)
k,kmax+1,,,h1+h2+h3
kmesh,kmax+1
nmax=node(0,0,h1+h2+h3)
nsel,s,loc,z,h1+h2+h3
cerig,nmax,all,roty