求压电分析的例子
求压电分析的例子 压电分析的例子 /POST1*GET,QT,NODE,NTOP,RF,CHRG ! GET TOTAL CHARGE ON TOP ELECTRODE
CP=ABS(QT) ! CAPACITANCE CP=Q/V, WHERE V=1V
EPZ0=8.854E-12 ! FREE SPACE PERMITTIVITY
PI=3.1415 ! PI CONSTANT
C=EP33*EPZ0*PI*A**2/T ! ANALYTICAL CAPACITANCE
/COM, 2-D CAPACITANCE (ANALYTICAL) =%C%, F
/COM, 2-D CAPACITANCE (ANSYS) = %CP%, F
FINI
/COM, ---------------------------------------------------------------
/COM, FINITE ELEMENT MODEL OF RLC-CIRCUIT
/COM, ---------------------------------------------------------------
/PREP7
DDELE,NTOP,VOLT ! DELETE VOLTAGE LOAD ON TOP ELECTRODE
ET,2,CIRCU94,0 ! DEFINE A RESISTOR
R=3000 ! RESISTANCE, OHM
R,1,R
N,1
TYPE,2
REAL,1
E,1,NTOP
ET,3,CIRCU94,1 ! DEFINE AN INDUCTOR
L=15 ! INDUCTANCE, H
R,2,L
N,2
TYPE,3
REAL,2
E,2,1
ET,4,CIRCU94,4 ! DEFINE A VOLTAGE SOURCE
V=1 ! VOLTAGE LOAD, V
R,3,V
N,3
TYPE,4
REAL,3
E,2,NBOT,3
FINI
/SOLU
ANTYPE,TRANS ! TRANSIENT ANALYSIS
NSUB,100 ! NUMBER OF TIME STEPS
TIME,2E-3 ! ANALYSIS TIME, S
TINTP,,0.25,0.5,0.5 ! INTEGRATION PARAMETERS FOR A PIEZOELECTRIC ANALYSIS
OUTRES,ALL,ALL
SOLVE
FINI
/GOPR
/COM, ANALYTICAL SOLUTION:
/COM, V_C = 1-EXP(-D*T)*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
/POST26
NUMVAR,20
NSOL,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(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)*
ADD,10,3,9,,V_C_EXACT ! 1-EXP(-D*T)*
! PRINT AND PLOT ANSYS AND EXACT VOLTAGE ACROSS THE PZT CAPACITOR
PRVAR,2,10
PLVAR,2,10
/NOPR
*DIM,VCE,ARRAY,5 ! EXACT SOLUTION FOR RESULTS TABLE
*DIM,VCA,ARRAY,5 ! ANSYS SOLUTION FOR RESULTS TABLE
*GET,VCE(1),VARI,10,RTIME,0.18E-3
*GET,VCE(2),VARI,10,RTIME,0.40E-3
*GET,VCE(3),VARI,10,RTIME,0.88E-3
*GET,VCE(4),VARI,10,RTIME,0.13E-2
*GET,VCE(5),VARI,10,RTIME,0.186E-2
*GET,VCA(1),VARI,2,RTIME,0.18E-3
*GET,VCA(2),VARI,2,RTIME,0.40E-3
*GET,VCA(3),VARI,2,RTIME,0.88E-3
*GET,VCA(4),VARI,2,RTIME,0.13E-2
*GET,VCA(5),VARI,2,RTIME,0.186E-2
*DIM,LABEL,CHAR,5
*DIM,VALUE,,5,3
LABEL(1) = '1.8E-2s','4.0E-2s','8.8E-2s','1.3E-1s','1.86E-1s'
*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 for t @ ...
/COM
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,' ',F12.4,' ',F12.4,' ',1F5.3)
/COM,----------------------------------------------------------------
/OUT
FINI
页:
[1]