ANSYS圆柱状压电模态分析，求助!后天答辩啊……

lantianshi1hao

各位大侠，我做了一个圆柱形的压电分析，在顶面施加面压力，但是在分析模态分析时，它的形变方式不是我施加压力的方向，而是沿侧面变形了，小弟附上程序与图片望各位大侠帮帮忙啊，小弟不胜感激……
程序如下：
/batch,list
/filname, Magnetostriction
/title, Static Analysis of  Magnetostriction
/NOPR
/COM,
/PREP7
D=35E-3
H1=1E-3

CYLIND,0,D/2,0,H1,0,360

! PZT-5H材料特性参数
MP, DENS, 1, 7750         !密度，Kg/m3
MP, PERX, 1, 919   !介电常数
MP, PERY, 1, 919
MP, PERZ, 1, 826
TB, ANEL, 1               !弹性劲度系数，N/m2
TBDATA, 1, 12.035E10, 7.5179E10, 7.509E10     
TBDATA, 7, 12.035E10, 7.509E10  
TBDATA, 12, 11.087E10
TBDATA, 16, 2.2584E10  
TBDATA, 19, 2.1053E10
TBDATA, 21, 2.1053E10   

TB, PIEZ, 1               !压电（应力）常数，C/m2
TBDATA, 3, -5.3512
TBDATA, 6, -5.3512
TBDATA, 9, 15.784
TBDATA, 14, 12.295
TBDATA, 16, 12.295

ET, 1, SOLID5,3

VSEL, S, VOLU, ,1
MAT, 1 $ TYPE, 1              !压电片网络划分

VSEL, ALL
wprota,,90  $ vsbw,all                      ! 切分圆柱体
wprota,,,90  $ vsbw,all                     ! 再切分圆柱体
mshape,0,3d  $ mshkey,1               ! 定义单元形状、网格划分类型
lsel,s,loc,z,-H1-H2  $ lesize,all,,,8             ! 每条圆周线网格划分数为 6
lsel,s,length,,2*H1+2*H2  $ lesize,all,,,8    ! 每条柱面侧线网格划分数为 8
vmesh,all
wprota,,,-90
wprota,,-90

NSEL, S, LOC, Z, 0
CP,1,VOLT, ALL
*GET, N_BOT, NODE, 0, NUM, MIN
D,N_BOT,VOLT,0

NSEL,S,LOC,Z,H1
CP,2,VOLT,ALL
*GET,N_TOP,NODE,0,NUM,MIN

NSEL,ALL
FINI
/SOLU                           !进行模态分析
ANTYPE, MODAL
MODOPT, LANB, 3
MXPAND, 3

NSEL,S,LOC,z,0
D,ALL,UZ,0
ASEL, S, LOC, Z,H1
SFA, ALL, 1, PRES, 15000

ALLSEL, ALL
SOLVE
FINISH
/POST1
SET, 1, 1               
PLDISP, 2                   !显示一阶频率

lantianshi1hao

拜托啦……

love5211109

你不加力的结果和你加力的结果是一样的，明白？