求ansys模拟压电陶瓷片的代码，或详细的操作步骤~非常感谢~~~！！！！！

FINISH/UNITS,SI/CLEAR/FILNAME,gl2mmsunshang/CONFIG,FSPLIT,750/CONFIG,NRES,4000/PREP7 !进入前处理器/TITLE,sunshangET,1,SOLID45 !定义单元类型ET,2,SOLID5 !PZT-5AKEYOPT,2,1,3 !设置单元关键选项MP,EX,1,206E11MP,DENS,1,7850 !定义1号单元材料属性MP,PRXY,1,03MP,DENS,2,7750 !定义压电片相关参数PZT-5AEMUNIT,EPZRO,885E-12MP,PERX,2,916 !定义介电常数 k11=916,k33=830MP,PERY,2,916 !相对介电常数k11(Z轴极化)MP,PERZ,2,830 !k33!定义压电应力系数e33=158 !单位为c/m2e31=-54e15=123!定义各向异性弹性刚度c11=121E10 !单位为N/m2c12=754E10c13=752E10c33=111E10c44=211E10c66=226E10 TB,PIEZ,2 !创建压电应力系数TBDATA,3,e31 !e31TBDATA,6,e31 !e31TBDATA,9,e33 !e33TBDATA,14,e15 !e15TBDATA,16,e15 !e15TB,ANEL,2,,,0 !创建各向异性弹性刚度TBDATA,1,c11,c12,c13 !C11,C12,C13TBDATA,7,c11,c13 !C11,C13TBDATA,12,c33 !C33TBDATA,16,c66 !C66TBDATA,19,c44 !C44TBDATA,21,c44 !C44BLOCK,0,06,0,002,0,002 !生成实体模型BLOCK,0,0008,0006,0014,002,00203BLOCK,0296,0304,0006,0014,002,00203BLOCK,0504,0506,0,002,0018,002 !损伤尺寸宽2mm，深2mm，长度20mmVGLUE,1,2,3VSBV,5,4AGEN,2,25,,,0008 !产生面来切割体AGEN,2,25,,,0296AGEN,2,25,,,0304AGEN,2,25,,,0504AGEN,2,25,,,0506AGEN,2,25,,,0018AGEN,2,25,,,0286AGEN,2,25,,,0314AGEN,2,25,,,0494AGEN,2,25,,,0516AGEN,2,27,,,,0006AGEN,2,27,,,,0014AGEN,2,1,,,,,0018VSBA,ALL,ALLVSEL,S,VOLU,,4,66 !赋予属性VATT,1,,1VSEL,S,VOLU,,2,3VATT,2,,2VSEL,ALLLSEL,S,LENGTH,,0018 !选择线并指定网格划分大小LESIZE,ALL,,,6,,,,,1LSEL,S,LENGTH,,0008LESIZE,ALL,,,4,,,,,1LSEL,S,LENGTH,,001LESIZE,ALL,,,4,,,,,1LSEL,S,LENGTH,,0268LESIZE,ALL,,,22,,,,,1LSEL,S,LENGTH,,0180LESIZE,ALL,,,20,,,,,1LSEL,S,LENGTH,,0002LESIZE,ALL,,,2,,,,,1LSEL,S,LENGTH,,0084LESIZE,ALL,,,12,,,,,1LSEL,S,LENGTH,,0006LESIZE,ALL,,,3,,,,,1LSEL,S,LENGTH,,00003LESIZE,ALL,,,3,,,,,1VSEL,ALL !网格划分MSHKEY,1MSHAPE,0,3DVMESH,ALLNSEL,S,LOC,X,002 !施加约束NSEL,R,LOC,Z,0D,ALL,UX,0,,,,UY,UZNSEL,S,LOC,X,0587NSEL,R,LOC,Z,0D,ALL,UY,0D,ALL,UZ,0SET,V1,0NSEL,S,LOC,Z,002 !耦合第一个压电片的电压NSEL,R,LOC,X,0,0008NSEL,R,LOC,Y,0006,0014CP,1,VOLT,ALL !耦合电压GET,N1,NODE,,NUM,MIND,N1,VOLT,V1 !施加电压荷载NSEL,S,LOC,Z,00203NSEL,R,LOC,X,0,0008NSEL,R,LOC,Y,0006,0014CP,2,VOLT,ALLGET,N2,NODE,,NUM,MINNSEL,S,LOC,Z,002 !耦合第二个压电片的电压NSEL,R,LOC,X,0296,0304NSEL,R,LOC,Y,0006,0014CP,3,VOLT,ALLGET,N3,NODE,,NUM,MIND,N3,VOLT,V1 !施加电压荷载NSEL,S,LOC,Z,00203NSEL,R,LOC,X,0296,0304NSEL,R,LOC,Y,0006,0014CP,4,VOLT,ALLGET,N4,NODE,,NUM,MINALLSEL,ALLFINISHSAVESET,pi,31415926 !第一个信号（汉宁窗调制正弦波）SET,f,801e3 !80KHZSET,T0,10/f !80KHZ电压信号 （窗口宽度为10周期）DIM,V2,table,800,1,1 !定义一个数组t1=1E-10 !起始时间t2=T0/200 !时间间隔t3=T0 !终止时间t4=4T0DO,i,1,800,1 Ti=t1+it2 V2(i,0)=Ti IF,Ti,LE,t3,THEN V2(i,1)=05(1-COS(2piTi/T0))SIN(2pifTi) ELSE V2(i,1)=0 ENDIFENDDOFINISH!进行模态分析/SOLUANTYPE,2MODOPT,LANB,50EQSLV,SPARMXPAND,50, , ,0LUMPM,0PSTRES,0MODOPT,LANB,50,0,10000000000000000, ,OFFSOLVEFINI!提取第一阶模态频率GET,FREQ1,MODE,1,FREQ!提取第二阶模态频率GET,FREQ2,MODE,2,FREQ/POST1/VIEW,1,1,1,1/VUP,,Z/TYPE,,4DO,I,1,50SET,LIST,IPLDISP,I !显示变形几何形状ENDDOFINISH/SOLU !进入求解器ALLSEL,ALLANTYPE,TRANS !指定分析类型TRNOPT,FULL !指定瞬态分析方法NLGEOM,OFF !-----关闭大变形选项TINTP,,025,05,05 !压电瞬态分析中，瞬态积分参数EQSLV,SPAROUTPR,ALL,ALL !控制结果的输出OUTRES,ALL,ALL !保存所有的数据到结果文件中TIMINT,ON,ALL !使用瞬态效应!DAMPRATIO=002 !阻尼比!F1=2pi13549 !第一阶圆频率!F2=2PI18578 !第二阶圆频率!ALPHAD,2DAMPRATIOF1F2/(F1F2) !Raileigh阻尼!BETAD,2DAMPRATIO/(F1F2)!MP,DAMP,1,000022/(231416f/10) !假想钢的阻尼!MP,DAMP,2,00152/(231416f/10) !假想PZT的阻尼time,t4 ! 第一载荷步nsubst,800ALLSEL,ALLD,N2,VOLT,%V2%KBC,0 !指定为斜坡加载方式SOLVEFINISH/POST26NSOL,2,N2,VOLT,,V2 !定义变量NSOL,3,N4,VOLT,,V4/AXLAB,X,Time ( s ) /AXLAB,Y,Voltage ( V )PLVAR,2/WAIT,3PLVAR,3/GRTYP,2/GROPT,ASCAL,ON XVAR,1PLVAR,2,3GET,NUM_VAR,VARI,0, NSETS ! 将时间变量长度值赋给变量NUM_VARDIM,T_U,ARRAY,NUM_VAR,3 ! 定义保存时间历程结果的矩阵T_U，有num_var行，4列VGET,T_U(1,1,1),1 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）VGET,T_U(1,2,1),2 ! V2的时间历程结果VGET,T_U(1,3,1),3 ! V4CFOPEN,T_U,txt ! 定义数据结果文件，T_UtxtVWRITE,T_U(1,1,1),T_U(1,2,1),T_U(1,3,1) ! 将时间历程结果数据写到数据结果文件中(31F169) ! 数据的输出格式GET,NUM_VAR,VARI,0, NSETS ! 将时间变量长度值赋给变量NUM_VARDIM,T,ARRAY,NUM_VAR,1 ! 定义保存时间历程结果的矩阵T_U，有num_var行，4列VGET,T(1,1,1),1 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）CFOPEN,T,txt ! 定义数据结果文件，T_UtxtVWRITE,T(1,1,1) ! 将时间历程结果数据写到数据结果文件中(31F169) ! 数据的输出格式GET,NUM_VAR,VARI,0, NSETS ! 将时间变量长度值赋给变量NUM_VARDIM,v2,ARRAY,NUM_VAR,1 ! 定义保存时间历程结果的矩阵T_U，有num_var行，4列VGET,v2(1,1,1),2 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）CFOPEN,v2,txt ! 定义数据结果文件，T_UtxtVWRITE,v2(1,1,1) ! 将时间历程结果数据写到数据结果文件中(31F169)GET,NUM_VAR,VARI,0, NSETS ! 将时间变量长度值赋给变量NUM_VARDIM,v4,ARRAY,NUM_VAR,1 ! 定义保存时间历程结果的矩阵T_U，有num_var行，4列VGET,v4(1,1,1),3 ! 矩阵的第一列保存第一个时间历程变量（默认为TIME）CFOPEN,v4,txt ! 定义数据结果文件，T_UtxtVWRITE,v4(1,1,1) ! 将时间历程结果数据写到数据结果文件中(31F169)