*HEADING
: STEADY STATE DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX8E]
RESTART ANALYSIS USING FREQUENCIES AND MODES FROM PPZOFRQ2 FOR MODAL STEPS
*RESTART,READ,STEP=1,INC=1,WRITE,FREQ=999
*SURFACE, NAME=BOTT
1,S1
2,S1
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*ELSET,ELSET=TOPELE
3,4
*DLOAD
TOPELE,P3,1.0E6
*EL PRINT
S,PHS
E,PHE
EPG,PHEPG
EFLX,PHEFL
*EL FILE
S,PHS,E,PHE,EPG,PHEPG,EFLX,PHEFL
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,E,EPG,EFLX
*NSET,NSET=TOPCEN
23, 
*NSET,NSET=BOTCEN
3, 
*NODE PRINT,NSET=TOPCEN
U,PU
EPOT,PHPOT
*NODE PRINT,NSET=BOTCEN
RF,PRF
RCHG,PHCHG
*NODE FILE,NSET=TOPCEN
U,PU,EPOT,PHPOT
*OUTPUT,FIELD
*NODE OUTPUT,NSET=TOPCEN
U,EPOT
*NODE FILE,NSET=BOTCEN
RF,PRF,RCHG,PHCHG
*OUTPUT,FIELD
*NODE OUTPUT,NSET=BOTCEN
RF,RCHG
*NODE PRINT,total=yes
RF,
*NODE FILE
RF,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*SECTION PRINT, NAME=BOTT, SURFACE=BOTT
SOF,SOAREA
*SECTION FILE, NAME=BOTT, SURFACE=BOTT
SOF,SOAREA
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DLOAD
TOPELE,P3,1.0E6
*NODE PRINT,total=yes
RF,
*NODE FILE
RF,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*SECTION PRINT, NAME=BOTT, SURFACE=BOTT
SOF,SOAREA
*SECTION FILE,  NAME=BOTT, SURFACE=BOTT
SOF,SOAREA
*END STEP