mehdikasaee
Electrical
Hi , In this simulation I try to connect a resistor to a piezoelectric crystal and measure the voltage input which is caused because of acceleration , I have written following codes but these codes does not work!!! the first problem is ansys does not solve these codes
units,MKS
/prep7
lb=50e-3
wb=6e-3
tb=0.3e-3
lp=13e-3
wp=6e-3
tp=1e-3
lm=2.3e-3
wm=6e-3
tm=6e-3
/prep7
locall,11
local,12,,,,,,180,,
et,1,solid45
et,2,solid5,3
et,3,circu94,0
CSYS,11
block,0,lb,0,wb,0,tb
block,0,lp,0,wp,tb,tb+tp
block,lb-lm,lb,0,wm,tb,tb+tm
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
EMUNIT, MKS, ! Define Piezoelectric #PZT PSI-5H4E) Material Model
MPTEMP, 1, 25, , , , ,
MPDATA,DENS,1,,7700
TB,PIEZ,1,,,0 ! Piezoelectric matrix
TBDATA, 1, 14.6913
TBDATA, 4, -2.8012
TBDATA, 7, -2.8012
TBDATA, 11, 9.8568
TBDATA, 18, 9.8568
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Clamped Dielectric Constant
EMUNIT, EPZRO, 8.85E-12
MPData, PERX, 1 , ,701
MPData, PERY, 1 , ,829
MPData, PERZ, 1 , ,829
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/com Pz26 X-polarized
/com Stiffness
TB, ANEL, 1 , 1 , 0
TBTEMP,25
TBDATA, 1, 1.2264E+11 , 9.9905E+10 , 9.9905E+10
TBDATA, 7, 1.6800E+11 , 1.1035E+11
TBDATA, 12, 1.6800E+11
TBDATA, 16, 3.0126E+10
TBDATA, 19, 2.8825E+10
TBDATA, 21, 3.0126E+10
!*
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
MPDATA,DENS,2,,2330 ! Density of Shim
MPDATA,EX,2,,1.11E11 ! Elastic Modulus of Shim
MPDATA,PRXY,2,,0.23 ! Poissons Ratio of Shim
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!alternative!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
!mp,dens,2,2330
!mp,ex,2,1.11e11
!mp,nuxy,2,0.23
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Define Mass #Tungsten Alloy) Material Model
MPDATA,DENS,3,,17000 ! Density of Mass
MPDATA,EX,3,,4E11 ! Elastic Modulus of Mass
!MPDATA,Ey,3,,4E11
!MPDATA,Ez,3,,4E11
MPDATA,PRXY,3,,0.28 ! Poissons Ratio of Mass
VSEL,All
VGLUE,ALL ! Combine all Volumes
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! DEFINING ELEMENT TYPE, AND mATERIAL PROPERTIES
Vsel,s,loc,z,0,tb,,
Vsel,r,loc,x,0,lb,,
nsel,r,loc,y,0,wb,,
Vatt, 2, , 1, 11
vsel,s,loc,z,tb,tb+tp
vsel,r,loc,x,0,lp,,
nsel,r,loc,y,0,wp,,
vatt, 1 ,, 2, 11
vsel,s,loc,z,tb,tb+tm,,
vsel,r,loc,x,lb-lm,lm,,
nsel,r,loc,y,0,wm,,
vatt,3 ,, 1, 11
esize,4e-4,0,
Vsel,all
Vsweep,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Node!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
n,n1,lb/2,wb/2,2*tm
n,n2,lb/2,wb/2,0
n,n3,Lp/2,wb/2,Tb+tp
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! electrode!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
NSEL,S,LOC,Z,0,0,tb, ,
NSEL,R,LOC,X,0,lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,1,VOLT,ALL ! Couple nodes of bottom
*GET,ELEC_BOT,node,0,num,min
NSEL,S,LOC,Z,Tb+Tp,Tb+Tp, ,
NSEL,R,LOC,X,0,Lp, ,
NSEL,R,LOC,Y,0,Wp, ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!REsistores!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
R,1,200
type,3
Real,1
e,7377,7378
R,2,200
type,3
real,2
e,7377,7379
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Fixing the end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ASEL,S,LOC,Z,0,Tb+TP,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
ASEL,S,LOC,Z,0,Tb+TP,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UY,0,,
ASEL,S,LOC,Z,0,Tb+TP,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UZ,0,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!solution!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOL
ANTYPE,HARMIC ,NEW ! Harmonic Analysis
HROPT,FULL
EQSLV,SPARSE ! Sparse Solver
KBC,1
D,ELEC_bot,VOLT,0, ! Setting Voltage on Middle to zero
ACEL,0,0,2.25, ! Acceleration
HARFRQ,104,106, ! Frequency Range
NSUBST,3, !Substeps
!*
ALPHAD,0,
BETAD,0,
DMPRAT,0.02, ! Constant Damping Ratio
!*
ALLSEL,ALL
/REPLOT
/STATUS,SOLU
SOLVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
units,MKS
/prep7
lb=50e-3
wb=6e-3
tb=0.3e-3
lp=13e-3
wp=6e-3
tp=1e-3
lm=2.3e-3
wm=6e-3
tm=6e-3
/prep7
locall,11
local,12,,,,,,180,,
et,1,solid45
et,2,solid5,3
et,3,circu94,0
CSYS,11
block,0,lb,0,wb,0,tb
block,0,lp,0,wp,tb,tb+tp
block,lb-lm,lb,0,wm,tb,tb+tm
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
EMUNIT, MKS, ! Define Piezoelectric #PZT PSI-5H4E) Material Model
MPTEMP, 1, 25, , , , ,
MPDATA,DENS,1,,7700
TB,PIEZ,1,,,0 ! Piezoelectric matrix
TBDATA, 1, 14.6913
TBDATA, 4, -2.8012
TBDATA, 7, -2.8012
TBDATA, 11, 9.8568
TBDATA, 18, 9.8568
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Clamped Dielectric Constant
EMUNIT, EPZRO, 8.85E-12
MPData, PERX, 1 , ,701
MPData, PERY, 1 , ,829
MPData, PERZ, 1 , ,829
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/com Pz26 X-polarized
/com Stiffness
TB, ANEL, 1 , 1 , 0
TBTEMP,25
TBDATA, 1, 1.2264E+11 , 9.9905E+10 , 9.9905E+10
TBDATA, 7, 1.6800E+11 , 1.1035E+11
TBDATA, 12, 1.6800E+11
TBDATA, 16, 3.0126E+10
TBDATA, 19, 2.8825E+10
TBDATA, 21, 3.0126E+10
!*
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
MPDATA,DENS,2,,2330 ! Density of Shim
MPDATA,EX,2,,1.11E11 ! Elastic Modulus of Shim
MPDATA,PRXY,2,,0.23 ! Poissons Ratio of Shim
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!alternative!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
!mp,dens,2,2330
!mp,ex,2,1.11e11
!mp,nuxy,2,0.23
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MPTEMP, 1, 25, , , , , ! Define Mass #Tungsten Alloy) Material Model
MPDATA,DENS,3,,17000 ! Density of Mass
MPDATA,EX,3,,4E11 ! Elastic Modulus of Mass
!MPDATA,Ey,3,,4E11
!MPDATA,Ez,3,,4E11
MPDATA,PRXY,3,,0.28 ! Poissons Ratio of Mass
VSEL,All
VGLUE,ALL ! Combine all Volumes
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! DEFINING ELEMENT TYPE, AND mATERIAL PROPERTIES
Vsel,s,loc,z,0,tb,,
Vsel,r,loc,x,0,lb,,
nsel,r,loc,y,0,wb,,
Vatt, 2, , 1, 11
vsel,s,loc,z,tb,tb+tp
vsel,r,loc,x,0,lp,,
nsel,r,loc,y,0,wp,,
vatt, 1 ,, 2, 11
vsel,s,loc,z,tb,tb+tm,,
vsel,r,loc,x,lb-lm,lm,,
nsel,r,loc,y,0,wm,,
vatt,3 ,, 1, 11
esize,4e-4,0,
Vsel,all
Vsweep,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Node!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
n,n1,lb/2,wb/2,2*tm
n,n2,lb/2,wb/2,0
n,n3,Lp/2,wb/2,Tb+tp
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! electrode!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
NSEL,S,LOC,Z,0,0,tb, ,
NSEL,R,LOC,X,0,lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,1,VOLT,ALL ! Couple nodes of bottom
*GET,ELEC_BOT,node,0,num,min
NSEL,S,LOC,Z,Tb+Tp,Tb+Tp, ,
NSEL,R,LOC,X,0,Lp, ,
NSEL,R,LOC,Y,0,Wp, ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!REsistores!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
R,1,200
type,3
Real,1
e,7377,7378
R,2,200
type,3
real,2
e,7377,7379
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Fixing the end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ASEL,S,LOC,Z,0,Tb+TP,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
ASEL,S,LOC,Z,0,Tb+TP,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UY,0,,
ASEL,S,LOC,Z,0,Tb+TP,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UZ,0,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!solution!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOL
ANTYPE,HARMIC ,NEW ! Harmonic Analysis
HROPT,FULL
EQSLV,SPARSE ! Sparse Solver
KBC,1
D,ELEC_bot,VOLT,0, ! Setting Voltage on Middle to zero
ACEL,0,0,2.25, ! Acceleration
HARFRQ,104,106, ! Frequency Range
NSUBST,3, !Substeps
!*
ALPHAD,0,
BETAD,0,
DMPRAT,0.02, ! Constant Damping Ratio
!*
ALLSEL,ALL
/REPLOT
/STATUS,SOLU
SOLVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!