*HEADING
: MODAL DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX4E]
RESTART ANALYSIS USING FREQUENCIES AND MODES FROM PPZOFRQ1
** DEFINE SINUSOIDAL AMPLITUDE FUNCTION
*AMPLITUDE,NAME=SIN,DEF=PERIODIC
1,100000.,0.
0.,1.
*RESTART,READ,STEP=1,INC=1,WRITE,FREQ=999
*STEP,PERTURBATION
*MODAL DYNAMIC
.01E-4,1.6E-4
*ELSET,ELSET=TOPELE,GEN
13,16
*DLOAD,AMP=SIN
TOPELE,P3,1.0E6
*EL PRINT,FREQ=1000
S,E
EPG,EFLX
*EL FILE,FREQ=1000
S,E,EFLX,EPG
*NODE PRINT,NSET=TOPCEN
U, 
*NODE FILE,NSET=TOPCEN
U,EPOT
*NODE PRINT,NSET=TOPCEN,FREQ=10
EPOT, 
*NODE FILE,NSET=BOTTOM,FREQ=1000
RF,RCHG
*NODE PRINT,NSET=BOTTOM,FREQ=1000
RF,RCHG
*END STEP