ANSYS圆柱状压电模态分析，求助!后天答辩啊……
各位大侠，我做了一个圆柱形的压电分析，在顶面施加面压力，但是在分析模态分析时，它的形变方式不是我施加压力的方向，而是沿侧面变形了，小弟附上程序与图片望各位大侠帮帮忙啊，小弟不胜感激……程序如下：
/batch,list
/filname, Magnetostriction
/title, Static Analysis ofMagnetostriction
/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 !显示一阶频率 拜托啦……{:aa22 你不加力的结果和你加力的结果是一样的，明白？
页:
[1]