*HEADING
: STEADY STATE DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX4E]
RESTART ANALYSIS OF PPZOFRQ1 USING STEADY STATE DYNAMICS, DIRECT
*RESTART,READ,STEP=1,INC=1,WRITE,FREQ=999
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT=REALONLY
15915.49, 
*ELSET,ELSET=TOPELE,GEN
13,16
*DECHARGE,REAL,OP=NEW
TOPELE,ES3,1.0E-4
*EL PRINT
S,PHS
E,PHE
EPG,PHEPG
EFLX,PHEFL
*ELSET,ELSET=ELOUT
13, 
*EL FILE,ELSET=ELOUT
S,PHS,E,PHE,EPG,PHEPG,EFLX,PHEFL,
*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
*NODE FILE,NSET=BOTCEN
RF,PRF,RCHG,PHCHG
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DECHARGE,IMAGINARY,OP=MOD
TOPELE,ES3,1.0E-4
*END STEP