shell181单元，分层叠合板，请大神给看看对不对

1203193201

finish
/clear
/prep7
k,1,
k,2,1,
k,3,1,1
k,4,,1
a,1,2,3,4
alls
wpcsys,,0
wprota,0,90,0  !绕x轴旋转
*do,i,1,9
wpave,0,0.1*i
asbw,all
*enddo
wpcsys,,0
wprota,0,0,90  
*do,i,1,9
wpave,0.1*i
asbw,all
*enddo

ET,1,shell181
MP,EX,1,3.25E10
MP,PRXY,1,0.2
MP,DENS,1,2500
ET,2,shell181
MP,EX,2,0.1502E10
MP,PRXY,2,0.2
MP,DENS,2,2400
keyopt,1,8,1
keyopt,2,8,1
sectype,1,shell
secdata,0.05,1
secdata,0.05,2
alls
aatt,,,,,1
mshape,0,2D
mshkey,1
allsel
lesize,all,1
amesh,all
eplot
/eshape,1
nsel,s,loc,y,0
d,all,all

save
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

我搞了个小例子，用181单元建立了个分层的板，最后模态分析和我用的solid45单元建立的模型算出来高频有一定的差距，前几阶频率很近，不知道命令流写的对不对，希望大神指教，181单元我不会用，看了半天help就是不确定写的对不！
谢谢，论坛兄弟姐妹！

1203193201

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

这是实体单元的模型，我对比了下

1203193201

壳单元的
sectype,1,shell
secdata,0.05,1
secdata,0.05,2
alls
改为
sectype,1,shell
secdata,0.3,1
secdata,0.3,2
alls

实体单元的
kgen,2,1,4,1,,,-0.05
kgen,2,1,4,1,,,0.05
改为
kgen,2,1,4,1,,,-0.3
kgen,2,1,4,1,,,0.3
差别就比较大了

1203193201

都没人回答啊！