-
1
- #1
KSJTam
Mechanical
- Nov 7, 2011
- 4
Hi, I am trying to create a cantilever piezoelectric bimorph. I have successfully created the structure and found the master nodes of the electrodes. I am trying to connect the bimorph with a parallel circuit. My load is just a single resistor. One end has to connect to the middle electrode and the other end has to connect with both the top and bottom electrodes. I'm not really sure how to connect the resistor. Should I couple the voltage for the top and bottom electrodes?
My code is here:
finish
/clear
/COM,Bimorph Model
/FILNAME, Bimorph,0
/TITLE, Bimorph Power Analysis
*SET,Lm,17.3E-3 !Length of mass
*SET,Hm,7.7E-3 !Height of mass
*SET,Wm,3.6E-3 !Width of mass
*SET,Lb,10.7E-3 !Length of beam
*SET,Wb,3.2E-3 !Width of beam
*SET,Tp,0.278E-3 !Pizeoelectric Thickness
*SET,Tsh,0.102E-3 !Shim Thickness
*SET,Rload,151E3 !Resistance
/UNITS,MKS ! MKS system (m, kg, s, °C#
/PREP7
local,11
local,12,,,,,,180,,
ET,1,SOLID45 !Set Element Types
ET,2,SOLID5,3
ET,3,CIRCU94,0
CSYS,11
BLOCK,0,Lb+Lm,0,Wb,0,Tsh, !Draw Shim
BLOCK,0,Lb,0,Wb,Tsh,Tsh+Tp, !Drawing Top Pizeoelectric Layer
BLOCK,0,Lb,0,Wb,0,-Tp, !Drawing Bottom Pizeoelectric Layer
BLOCK,Lb,Lb+Lm,#(Wb/2#-(Wm/2##,##Wb/2#+#Wm/2##,Tsh,Tsh+Hm !Drawing Mass
EMUNIT, MKS, ! Define Piezoelectric #PZT PSI-5H4E) Material Model
MPTEMP, 1, 25, , , , , ! PZT Density
MPDATA,DENS,1,,7500
TB,PIEZ,1,,,0 ! Piezoelectric matrix
TBDATA, 3, -12
TBDATA, 6, -12
TBDATA, 9, 18.2
TBDATA, 14, 21.9
TBDATA, 16, 21.9
MPTEMP, 1, 25, , , , , ! Clamped Dielectric Constant
MPDATA,PERX,1,,1550
MPDATA,PERY,1,,1550
MPDATA,PERZ,1,,1390
TB, ANEL, 1,1,0 ! Anisotropic elastic matrix
TBTEMP,25
TBDATA, 1, 16.9E+10 , 11.8E+10 ,10.9E+10
TBDATA, 7, 16.9E+10 , 10.9E+10
TBDATA, 12, 12.3E+10
TBDATA, 16, 2.5E+10
TBDATA, 19, 2.7E+10
TBDATA, 21, 2.7E+10
!*
MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
MPDATA,DENS,2,,8530 ! Density of Shim
MPDATA,EX,2,,1.11E11 ! Elastic Modulus of Shim
MPDATA,PRXY,2,,0.375 ! Poissons Ratio of Shim
MPTEMP, 1, 25, , , , , ! Define Mass #Tungsten Alloy) Material Model
MPDATA,DENS,3,,17000 ! Density of Mass
MPDATA,EX,3,,4.11E11 ! Elastic Modulus of Mass
MPDATA,PRXY,3,,0.28 ! Poissons Ratio of Mass
VSEL,All
VGLUE,ALL ! Combine all Volumes
VSEL,S,LOC,Z,0,Tsh, , ! Defining Element Type & Material Properties of layers
VSEL,R,LOC,X,0,Lb+Lm, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 2, , 1, 11 !Shim Layer
VSEL,S,LOC,Z,Tsh,Tsh+Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 1, , 2, 11 !Top Pizeoelectric Layer
VSEL,S,LOC,Z,0,-Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 1, , 2, 11 !Bottom Pizeoelectric Layer
VSEL,S,LOC,Z,Tsh,Tsh+Hm, ,
NSEL,R,LOC,X,Lb,Lb+Lm, ,
VATT, 3, , 1, 11 ! Mass
ESIZE,4E-4,0, ! Mesh Size
VSEL,ALL
VSWEEP,All !Sweep Volumes
NSEL,S,LOC,Z,-Tp,-Tp, ,
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,Tsh+Tp,Tsh+Tp, ,
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
NSEL,S,LOC,Z,0,0
NSEL,R,LOC,X,0,Lb, , ! Creating components
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_1,NODE
NSEL,S,LOC,Z,Tsh,Tsh
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_2,NODE
CMSEL,S,CM_1
CMSEL,A,CM_2
CP,3,VOLT,ALL ! Couple nodes of middle
*GET,ELEC_MID,node,0,num,min
ALLSEL,ALL
/REPLOT
R,1,Rload, , , , , ! Define Resistor
RMORE, , , , , ,
RMORE, , , ,
!*
TYPE, 3
MAT, 0
REAL, 1
ESYS, 11
SECNUM,
TSHAP,LINE
!*
EN,99999,ELEC_TOP,ELEC_MID,, ! Creating Resistor
!*
ALLSEL,ALL
/REPLOT
ASEL,S,LOC,Z,0,Tsh,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
ASEL,S,LOC,Z,0,Tsh,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UY,0,,
ASEL,S,LOC,Z,0,Tsh,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UZ,0,,
ALLSEL,ALL
/REPLOT
FINISH
/SOL
ANTYPE,HARMIC ,NEW ! Harmonic Analysis
HROPT,FULL
EQSLV,SPARSE ! Sparse Solver
KBC,1
D,ELEC_MID,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
/POST26 ! Plot Power Graph
FILE,'Bimorph','rst','.'
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
PLCPLX,0
PRCPLX,1
FILLDATA,191,,,,1,1
REALVAR,191,191
!*
!*
!*
FORCE,TOTAL
SHELL,TOP
LAYERP26,0
ESOL,2,99999, ,NMIS,1,PWR
!*
PLVAR,PWR, , , , , , , , , ,
My code is here:
finish
/clear
/COM,Bimorph Model
/FILNAME, Bimorph,0
/TITLE, Bimorph Power Analysis
*SET,Lm,17.3E-3 !Length of mass
*SET,Hm,7.7E-3 !Height of mass
*SET,Wm,3.6E-3 !Width of mass
*SET,Lb,10.7E-3 !Length of beam
*SET,Wb,3.2E-3 !Width of beam
*SET,Tp,0.278E-3 !Pizeoelectric Thickness
*SET,Tsh,0.102E-3 !Shim Thickness
*SET,Rload,151E3 !Resistance
/UNITS,MKS ! MKS system (m, kg, s, °C#
/PREP7
local,11
local,12,,,,,,180,,
ET,1,SOLID45 !Set Element Types
ET,2,SOLID5,3
ET,3,CIRCU94,0
CSYS,11
BLOCK,0,Lb+Lm,0,Wb,0,Tsh, !Draw Shim
BLOCK,0,Lb,0,Wb,Tsh,Tsh+Tp, !Drawing Top Pizeoelectric Layer
BLOCK,0,Lb,0,Wb,0,-Tp, !Drawing Bottom Pizeoelectric Layer
BLOCK,Lb,Lb+Lm,#(Wb/2#-(Wm/2##,##Wb/2#+#Wm/2##,Tsh,Tsh+Hm !Drawing Mass
EMUNIT, MKS, ! Define Piezoelectric #PZT PSI-5H4E) Material Model
MPTEMP, 1, 25, , , , , ! PZT Density
MPDATA,DENS,1,,7500
TB,PIEZ,1,,,0 ! Piezoelectric matrix
TBDATA, 3, -12
TBDATA, 6, -12
TBDATA, 9, 18.2
TBDATA, 14, 21.9
TBDATA, 16, 21.9
MPTEMP, 1, 25, , , , , ! Clamped Dielectric Constant
MPDATA,PERX,1,,1550
MPDATA,PERY,1,,1550
MPDATA,PERZ,1,,1390
TB, ANEL, 1,1,0 ! Anisotropic elastic matrix
TBTEMP,25
TBDATA, 1, 16.9E+10 , 11.8E+10 ,10.9E+10
TBDATA, 7, 16.9E+10 , 10.9E+10
TBDATA, 12, 12.3E+10
TBDATA, 16, 2.5E+10
TBDATA, 19, 2.7E+10
TBDATA, 21, 2.7E+10
!*
MPTEMP, 1, 25, , , , , ! Define Shim #Brass) Material Model
MPDATA,DENS,2,,8530 ! Density of Shim
MPDATA,EX,2,,1.11E11 ! Elastic Modulus of Shim
MPDATA,PRXY,2,,0.375 ! Poissons Ratio of Shim
MPTEMP, 1, 25, , , , , ! Define Mass #Tungsten Alloy) Material Model
MPDATA,DENS,3,,17000 ! Density of Mass
MPDATA,EX,3,,4.11E11 ! Elastic Modulus of Mass
MPDATA,PRXY,3,,0.28 ! Poissons Ratio of Mass
VSEL,All
VGLUE,ALL ! Combine all Volumes
VSEL,S,LOC,Z,0,Tsh, , ! Defining Element Type & Material Properties of layers
VSEL,R,LOC,X,0,Lb+Lm, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 2, , 1, 11 !Shim Layer
VSEL,S,LOC,Z,Tsh,Tsh+Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 1, , 2, 11 !Top Pizeoelectric Layer
VSEL,S,LOC,Z,0,-Tp, ,
VSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
VATT, 1, , 2, 11 !Bottom Pizeoelectric Layer
VSEL,S,LOC,Z,Tsh,Tsh+Hm, ,
NSEL,R,LOC,X,Lb,Lb+Lm, ,
VATT, 3, , 1, 11 ! Mass
ESIZE,4E-4,0, ! Mesh Size
VSEL,ALL
VSWEEP,All !Sweep Volumes
NSEL,S,LOC,Z,-Tp,-Tp, ,
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,Tsh+Tp,Tsh+Tp, ,
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
NSEL,S,LOC,Z,0,0
NSEL,R,LOC,X,0,Lb, , ! Creating components
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_1,NODE
NSEL,S,LOC,Z,Tsh,Tsh
NSEL,R,LOC,X,0,Lb, ,
NSEL,R,LOC,Y,0,Wb, ,
ESLN,S,1
CM,cm_2,NODE
CMSEL,S,CM_1
CMSEL,A,CM_2
CP,3,VOLT,ALL ! Couple nodes of middle
*GET,ELEC_MID,node,0,num,min
ALLSEL,ALL
/REPLOT
R,1,Rload, , , , , ! Define Resistor
RMORE, , , , , ,
RMORE, , , ,
!*
TYPE, 3
MAT, 0
REAL, 1
ESYS, 11
SECNUM,
TSHAP,LINE
!*
EN,99999,ELEC_TOP,ELEC_MID,, ! Creating Resistor
!*
ALLSEL,ALL
/REPLOT
ASEL,S,LOC,Z,0,Tsh,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
ASEL,S,LOC,Z,0,Tsh,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UY,0,,
ASEL,S,LOC,Z,0,Tsh,,
ASEL,R,LOC,X,0,0,,
DA,ALL,UZ,0,,
ALLSEL,ALL
/REPLOT
FINISH
/SOL
ANTYPE,HARMIC ,NEW ! Harmonic Analysis
HROPT,FULL
EQSLV,SPARSE ! Sparse Solver
KBC,1
D,ELEC_MID,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
/POST26 ! Plot Power Graph
FILE,'Bimorph','rst','.'
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
PLCPLX,0
PRCPLX,1
FILLDATA,191,,,,1,1
REALVAR,191,191
!*
!*
!*
FORCE,TOTAL
SHELL,TOP
LAYERP26,0
ESOL,2,99999, ,NMIS,1,PWR
!*
PLVAR,PWR, , , , , , , , , ,