*HEADING
: MODAL DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX8E]
RESTART ANALYSIS USING FREQUENCIES AND MODES FROM PPZOFRQ2
** DEFINE SINUSOIDAL AMPLITUDE FUNCTION
*AMPLITUDE,NAME=SIN,DEF=PERIODIC
1,100000.,0.
0.,1.
*RESTART,READ,STEP=1,INC=1,WRITE,FREQ=999
*SURFACE, NAME=BOTT
1,S1
2,S1
*STEP,PERTURBATION
*MODAL DYNAMIC
.01E-4,1.6E-4
*ELSET,ELSET=TOPELE
3,4
*DLOAD,AMP=SIN
TOPELE,P3,1.0E6
*EL PRINT,FREQ=1000
S,E
EPG,EFLX
*EL FILE,FREQ=1000
S,E,EFLX,EPG
*OUTPUT,FIELD,FREQ=1000
*ELEMENT OUTPUT
S,E,EFLX,EPG
*NSET,NSET=TOPCEN
23, 
*NODE PRINT,NSET=TOPCEN,FREQ=1000
U,EPOT
*NODE FILE,NSET=TOPCEN,FREQ=1000
U,EPOT
RF,
*OUTPUT,FIELD,FREQ=1000
*NODE OUTPUT,NSET=TOPCEN
U,EPOT
RF,
*NODE FILE,NSET=BOTTOM,FREQ=1000
RF,RCHG
*OUTPUT,FIELD,FREQ=1000
*NODE OUTPUT,NSET=BOTTOM
RF,RCHG
*NODE PRINT,totals=yes,FREQ=1000
U,EPOT
RF,
*NODE FILE,FREQ=1000
U,EPOT
RF,RCHG
*OUTPUT,FIELD,FREQ=1000
*NODE OUTPUT
U,EPOT
RF,RCHG
*SECTION PRINT, NAME=BOTT, SURFACE=BOTT,FREQ=1000
SOF,SOAREA
*SECTION FILE,  NAME=BOTT, SURFACE=BOTT,FREQ=1000
SOF,SOAREA
*END STEP