*HEADING
: MODAL DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX4E]
INCLUDES RAYLEIGH DAMPING
*NODE
*NGEN
*ELEMENT,TYPE=CAX4E,ELSET=ALL
*ELGEN,ELSET=ALL
**local orientation matching global system
*ORIENTATION,NAME=RECT
*SOLID SECTION,MATERIAL=PZT,ELSET=ALL,ORIENT=RECT
*MATERIAL,NAME=PZT
*ELASTIC,TYPE=ORTHO
13.90E10,7.428E10,11.54E10,7.784E10,7.428E10,13.90E10,2.564E10,2.564E10,
2.564E10,
*PIEZOELECTRIC,TYPE=S
*DIELECTRIC,TYPE=ANISO
*DENSITY
5555.,
*NSET,NSET=ELECTROD,GEN
*EQUATION
ELECTROD,9,1.0,21,9,-1.0
*NSET,NSET=BOTTOM,GEN
*BOUNDARY
BOTTOM,2
BOTTOM,9
*RESTART,WRITE,FREQ=999
** DEFINE SINUSOIDAL AMPLITUDE FUNCTION
*AMPLITUDE,NAME=SIN,DEF=PERIODIC
1,100000.,0.
0.,1.
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos
45, ,
*ELSET,ELSET=TOPELE
*EL PRINT,MODE=1,ELSET=TOPELE
S,E
EPG,EFLX
*EL FILE,LASTMODE=1,ELSET=TOPELE
S,E,EFLX,EPG
*NSET,NSET=TOPCEN
*NODE PRINT,NSET=TOPCEN,MODE=2
U,EPOT
*NODE FILE,NSET=TOPCEN,MODE=2
U,EPOT
*NSET,NSET=B3
*NODE FILE,NSET=B3,LASTMODE=1
RF,RCHG
*END STEP
*STEP,PERTURBATION
*MODAL DYNAMIC
.01E-4,1.6E-4
*MODAL DAMPING,RAYLEIGH
*ELSET,ELSET=TOPELE,GEN
*DLOAD,AMP=SIN
TOPELE,P3,1.0E6
*EL PRINT,FREQ=1000
S,E
EPG,EFLX
*EL FILE,FREQ=1000
S,E,EFLX,EPG
*NSET,NSET=TOPCEN
*NODE PRINT,NSET=TOPCEN
U,EPOT
*NODE FILE,NSET=TOPCEN
U,
*NODE FILE,NSET=TOPCEN,FREQ=10
EPOT,
*NODE FILE,NSET=BOTTOM,FREQ=1000
RF,RCHG
*END STEP