这是我的部分命令流：/prep7
et,1,plane42,,,1
et,2,solid5,3
MP,ex,1,9.8e10 !黄铜
MP,nuxy,1,.32
MP,dens,1,7500
MP,DENS,2,7600 !压电陶瓷密度
MP,PERX,2,7.124e-9 !X和Y轴向的介电常数
MP,PERZ,2,5.841E-9 !Z轴向的介电常数
TB,PIEZ,2 !压电电材料参数定义
TBDATA,16,10.5
TBDATA,14,10.5
TBDATA,3,-4.1
TBDATA,6,-4.1
TBDATA,9,14.1
TB,ANEL,2
TBDATA,1,13.2E10,7.1E10,7.3E10
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA, 21, 2.6E10
!压电层几何模型 (R1=5mm,R2=9mm, H1 =1mm)
R1=5e-3
R2=9e-3
H1=1e-3
H2=1e-3


! 建立几何模型
cylind,R1,0,H2,0,360 !铜片层
cylind,R1,R2,H2,H1+H2,0,360 !PZT
!进行网格划分
LESIZE,ALL,10e-4, , , ,1, , ,1,
TYPE,2
MAT,2
ALLSEL,ALL
VSEL,S,LOC,Z,0,H2
VSWEEP,ALL
ALLSEL,ALL
TYPE,1
MAT,1
ALLSEL,ALL
VSEL,S,LOC,Z,H2,H1+H2
VSWEEP,ALL

ALLSEL,ALL
NSEL,S,LOC,Z,0
D,ALL,UZ
ALLSEL,ALL
allsel,all
finish
/solu !瞬态分析
antype,trans
trnopt,full
tintp,,0.25,0.5,0.5
outpr,all,all
outres,all,all
kbc,1
nsel,s,loc,z,H2 !耦合电压
d,all,volt,0
cp,1,volt,all
nsel,s,loc,Z,H1+H2
cp,2,volt,all