急！急！对压电材料pzt-4的模型建立，请问有没有问题

tzhdacy

对压电材料pzt-4的模型建立，请问有没有问题
/prep7
!
/units,mks
!*******************************
ET,1,45              !  element type 45
ET,2,solid5          !  element type pzt
ET,3,42              !  2-D structural solid
!**************material 1 : PZT5H***********
MP,DENS,1,7500
MP,PERX,1,1700
MP,PERZ,1,1470
TB,PIEZ,1
TBDATA,3,-6.55
TBDATA,6,-6.55
TBDATA,9,23.3
TBDATA,14,17.0
TBDATA,16,17.0
TB,ANEL,1
TBDATA,1,12.6E10,7.95E10,8.41E10
TBDATA,7,12.6E10,8.41E10
TBDATA,12,11.7E10
TBDATA,16,2.35E10
TBDATA,19,2.30E10
TBDATA,21,2.30E10
MP,MURX,1,1
MP,KXX,1,1
!==============================
L     =2e-2
H     =2e-3
B     =2e-2
SIZE  =1e-3
!==============================
K,1,0,-H/2,-B/2
K,2,0,H/2,-B/2
K,3,0,H/2,B/2
K,4,0,-H/2,B/2
K,5,L,-H/2,-B/2
K,6,L,H/2,-B/2
K,7,L,H/2,B/2
K,8,L,-H/2,B/2
V,1,2,3,4,5,6,7,8
!==============================
LESIZE,1,,,4
LESIZE,3,,,4
LESIZE,6,,,4
LESIZE,10,,,4
LESIZE,2,,,6
LESIZE,5,,,6
LESIZE,8,,,6
LESIZE,9,,,6
LESIZE,4,,,6
LESIZE,7,,,6
LESIZE,11,,,6
LESIZE,12,,,6
VMESH,ALL
DSYM,SYMM,X
NSEL,S,LOC,x,0                ! SYMMETRY BOUNDARY CONSTRAINTS
NSEL,S,LOC,Y,0
DSYM,SYMM,Y
NSEL,ALL
FINISH
/SOLU

另外，要对其作HARMONIC分析，请指点电压加入时候应该注意事项。因为明天要交课题，请各位大哥帮忙

tzhdacy

斑竹，帮忙给点建议