Ansys壓電分析

1、/POST1*GET,QT,NODE,NTOP,RF,CHRG ! GET TOTAL CHARGE ON TOP ELECTRODE CP=ABS(QT) ! CAPACITANCE CP=Q/V, WHERE V=1VEPZ0= ! FREE SPACE PERMITTIVITYPI= ! PI CONSTANTC=EP33*EPZ0*PI*A*2/T ! ANALYTICAL CAPACITANCE/COM, 2-D CAPACITANCE (ANALYTICAL) =%C%, F/COM, 2-D CAPACITANCE (ANSYS) = %CP%, FFINI/COM, -/COM

2、, FINITE ELEMENT MODEL OF RLC-CIRCUIT/COM, -/PREP7DDELE,NTOP,VOLT ! DELETE VOLTAGE LOAD ON TOP ELECTRODE ET,2,CIRCU94,0 ! DEFINE A RESISTORR=3000 ! RESISTANCE, OHMR,1,RN,1TYPE,2REAL,1E,1,NTOPET,3,CIRCU94,1 ! DEFINE AN INDUCTORL=15 ! INDUCTANCE, HR,2,LN,2TYPE,3REAL,2E,2,1ET,4,CIRCU94,4 ! DEFINE A VOL

3、TAGE SOURCEV=1 ! VOLTAGE LOAD, VR,3,VN,3TYPE,4REAL,3E,2,NBOT,3FINI/SOLUANTYPE,TRANS ! TRANSIENT ANALYSISNSUB,100 ! NUMBER OF TIME STEPSTIME,2E-3 ! ANALYSIS TIME, STINTP, ! INTEGRATION PARAMETERS FOR A PIEZOELECTRIC ANALYSISOUTRES,ALL,ALLSOLVEFINI/GOPR/COM, ANALYTICAL SOLUTION:/COM, V_C = 1-EXP(-D*T)

4、*COS(B*T)-D/B*EXP(-D*T)*SIN(B*T)/COM, WHERE:D=R/(2*L)B=SQRT(1/(L*C)-D*2)/NOPR/POST26NUMVAR,20NSOL,2,NTOP,VOLT,V_C_ANSYS! DERIVE EXACT SOLUTION *DIM,WORK1,ARRAY,100 *DIM,WORK2,ARRAY,100 FILLDATA,3,1 ! 1 EXP,4,1,-D,-1 ! -EXP(-D*T) PROD,5,1,B ! B*T VGET,WORK1(1),5 *VFUN,WORK2(1),COS,WORK1(1) VPUT,WORK2

5、(1),6 ! COS(B*T)*VFUN,WORK2(1),SIN,WORK1(1)VPUT,WORK2(1),7 ! SIN(B*T)ADD,8,6,7,D/B ! COS(B*T) + D/B*SIN(B*T)PROD,9,4,8 ! -EXP(-D*T)*COS(B*T) + D/B*SIN(B*T)! PRINT AND PLOT ANSYS AND EXACT VOLTAGE ACROSS THE PZT CAPACITOR PRVAR,2,10PLVAR,2,10/NOPR*DIM,VCE,ARRAY,5 ! EXACT SOLUTION FOR RESULTS TABLE *D

6、IM,VCA,ARRAY,5 ! ANSYS SOLUTION FOR RESULTS TABLE *GET,VCE(1),VARI,10,RTIME, *GET,VCE(2),VARI,10,RTIME, *GET,VCE(3),VARI,10,RTIME, *GET,VCE(4),VARI,10,RTIME, *GET,VCE(5),VARI,10,RTIME, *GET,VCA(1),VARI,2,RTIME, *GET,VCA(2),VARI,2,RTIME, *GET,VCA(3),VARI,2,RTIME, *GET,VCA(4),VARI,2,RTIME, *GET,VCA(5)

7、,VARI,2,RTIME,*DIM,LABEL,CHAR,5*DIM,VALUE,5,3 LABEL(1) = '','','','','' *VFILL,VALUE(1,1),DATA,VCE(1),VCE(2),VCE(3),VCE(4),VCE(5) *VFILL,VALUE(1,2),DATA,VCA(1),VCA(2),VCA(3),VCA(4),VCA(5) *VFILL,VALUE(1,3),DATA,ABS(VCA(1)/VCE(1),ABS(VCA(2)/VCE(2),ABS(VCA(3)/VCE(3), ABS(VCA(4)/VCE(4),ABS(VCA(5)/VCE(5)/OUT,vm237,vrt/COM/COM,- VM237 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO /COM,/COM, VC