Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations KootK on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

Base Vibration excitation of Cantilever Beam

Status
Not open for further replies.

KSJTam

Mechanical
Nov 7, 2011
4
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, , , , , , , , , ,
 
Replies continue below

Recommended for you

I have also same problem ,,,

My problem is how to excite the base beam?,,,what is the post-processor procedure in ANSYS to get induced voltage in Piezoelectric material?,,
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor