Hi, I've been trying to simulate the power output of a piezoelectric layer on a cantilever beam undergoing base excitation.
In the paper I'm trying to duplicate their ANSYS model has a base displacement of 0.1mm and uses a harmonic analysis to sweep the frequencies.
My question is, How do I apply a base displacement to my model? Should I use the D command and apply it to the fixed end only, or to all the nodes of my model?
I have attatched my entire code, however my problem lies in the /SOL section.
finish
/clear
/COM,Unimorph Cantilever Model: By Kenneth S.J. Tam
/FILNAME,Unimorph,0
/TITLE, Unimorph Power Analysis
*SET,L,0.0485
*SET,W,0.01
*SET,T,4E-4
*SET,PL,7E-3
*SET,PT,0.5E-4
*SET,PW,0.01
/PREP7
ET,1,SOLID226,1001 !Piezo
ET,2,SOLID186 !Structural
ET,3,CIRCU94,0
EMUNIT, MKS, ! Define Piezoelectric (PZT PSI-5H4E) Material Model 1
MPTEMP, 1, 0, , , , , ! PZT Density
MPDATA,DENS,1,,7600
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,0, , , , , ! Clamped Dielectric Constant
MPDATA,PERX,1,,1550
MPDATA,PERY,1,,1550
MPDATA,PERZ,1,,1390
TB, ANEL, 1,1,0 ! Anisotropic elastic matrix
TBTEMP,0
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, 0, , , , , ! Define Structral (Brass) Material Model 2
MPDATA,DENS,2,,8410
MPDATA,EX,2,,9.2E10
MPDATA,PRXY,2,,0.35
BLOCK,0,L,0,W,0,T, !Brass !Mass
BLOCK,0,PL,0,PW,T,T+PT, !PZT
VSEL,All
VGLUE,ALL ! Combine all Volumes
VSEL,S,LOC,Z,0,T, , ! Defining Material Properties of layers
VATT, 2, , 2, 0
VSEL,S,LOC,Z,T,T+PT, ,
VATT, 1, , 1, 0
MSHAPE,1,3D
MSHKEY,0
MSHMID,0
ESIZE,1E-3,0 ! Mesh
VSEL,ALL
VSWEEP,All
NSEL,S,LOC,Z,(T),(T), ,
NSEL,R,LOC,X,(0),(PL), ,
CP,1,VOLT,ALL ! Couple nodes of bottom
*GET,ELEC_BOT,node,0,num,min
NSEL,S,LOC,Z,(T+PT),(T+PT), ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
!*
R,1,400000, , , , , ! Define Resistor
RMORE, , , , , ,
RMORE, , , ,
!*
TYPE, 3
MAT, 0
REAL, 1
ESYS, 0
SECNUM,
TSHAP,LINE
!*
EN,99999,ELEC_TOP,ELEC_BOT ! Drawing Resistor
!*
ALLSEL,ALL
/REPLOT
ASEL,S,LOC,Z,0,T+PT,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
DA,ALL,UY,0,,
FINISH ! My problem is with this section
/SOL
ANTYPE,HARMIC,NEW ! Harmonic Analysis
HROPT,FULL
EQSLVE,SPARSE
KBC,1
D,ELEC_BOT,VOLT,0, !Bottom Electrode Zero
D,ALL,UZ,0.001,, !This displaces all the nodes
HARFRQ,100,140, !Frequency Sweep
NSUBST,30, !Substeps
!*
ALPHAD,0,
BETAD,0,
DMPRAT,0.01, !Damping
!*
ALLSEL,ALL
/REPLOT
/STATUS,SOLU
SOLVE
/POST26 ! Plot Power Graph
FILE,'Unimorph','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, , , , , , , , , ,
In the paper I'm trying to duplicate their ANSYS model has a base displacement of 0.1mm and uses a harmonic analysis to sweep the frequencies.
My question is, How do I apply a base displacement to my model? Should I use the D command and apply it to the fixed end only, or to all the nodes of my model?
I have attatched my entire code, however my problem lies in the /SOL section.
finish
/clear
/COM,Unimorph Cantilever Model: By Kenneth S.J. Tam
/FILNAME,Unimorph,0
/TITLE, Unimorph Power Analysis
*SET,L,0.0485
*SET,W,0.01
*SET,T,4E-4
*SET,PL,7E-3
*SET,PT,0.5E-4
*SET,PW,0.01
/PREP7
ET,1,SOLID226,1001 !Piezo
ET,2,SOLID186 !Structural
ET,3,CIRCU94,0
EMUNIT, MKS, ! Define Piezoelectric (PZT PSI-5H4E) Material Model 1
MPTEMP, 1, 0, , , , , ! PZT Density
MPDATA,DENS,1,,7600
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,0, , , , , ! Clamped Dielectric Constant
MPDATA,PERX,1,,1550
MPDATA,PERY,1,,1550
MPDATA,PERZ,1,,1390
TB, ANEL, 1,1,0 ! Anisotropic elastic matrix
TBTEMP,0
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, 0, , , , , ! Define Structral (Brass) Material Model 2
MPDATA,DENS,2,,8410
MPDATA,EX,2,,9.2E10
MPDATA,PRXY,2,,0.35
BLOCK,0,L,0,W,0,T, !Brass !Mass
BLOCK,0,PL,0,PW,T,T+PT, !PZT
VSEL,All
VGLUE,ALL ! Combine all Volumes
VSEL,S,LOC,Z,0,T, , ! Defining Material Properties of layers
VATT, 2, , 2, 0
VSEL,S,LOC,Z,T,T+PT, ,
VATT, 1, , 1, 0
MSHAPE,1,3D
MSHKEY,0
MSHMID,0
ESIZE,1E-3,0 ! Mesh
VSEL,ALL
VSWEEP,All
NSEL,S,LOC,Z,(T),(T), ,
NSEL,R,LOC,X,(0),(PL), ,
CP,1,VOLT,ALL ! Couple nodes of bottom
*GET,ELEC_BOT,node,0,num,min
NSEL,S,LOC,Z,(T+PT),(T+PT), ,
CP,2,VOLT,ALL ! Couple nodes of top
*GET,ELEC_TOP,node,0,num,min
!*
R,1,400000, , , , , ! Define Resistor
RMORE, , , , , ,
RMORE, , , ,
!*
TYPE, 3
MAT, 0
REAL, 1
ESYS, 0
SECNUM,
TSHAP,LINE
!*
EN,99999,ELEC_TOP,ELEC_BOT ! Drawing Resistor
!*
ALLSEL,ALL
/REPLOT
ASEL,S,LOC,Z,0,T+PT,, ! Fixing End of Shim
ASEL,R,LOC,X,0,0,,
DA,ALL,UX,0,,
DA,ALL,UY,0,,
FINISH ! My problem is with this section
/SOL
ANTYPE,HARMIC,NEW ! Harmonic Analysis
HROPT,FULL
EQSLVE,SPARSE
KBC,1
D,ELEC_BOT,VOLT,0, !Bottom Electrode Zero
D,ALL,UZ,0.001,, !This displaces all the nodes
HARFRQ,100,140, !Frequency Sweep
NSUBST,30, !Substeps
!*
ALPHAD,0,
BETAD,0,
DMPRAT,0.01, !Damping
!*
ALLSEL,ALL
/REPLOT
/STATUS,SOLU
SOLVE
/POST26 ! Plot Power Graph
FILE,'Unimorph','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, , , , , , , , , ,
