*HEADING
: STEADY STATE DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX4E]
RESTART ANALYSIS USING FREQUENCIES AND MODES FROM PPZOFRQ1 FOR MODAL STEPS
*RESTART,READ,WRITE,FREQ=999
**,STEP=1,INC=1,WRITE,FREQ=999
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*ELSET,ELSET=TOPELE,GEN
13,16
*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
*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, 
*DLOAD
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT=REALONLY
15915.49, 
*DLOAD
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos, SIM=NO
45, , 
*NODE FILE,F=0
*EL PRINT,F=0
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49, 
*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
*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,SUBSPACE PROJECTION,REAL ONLY
15915.49, 
*DLOAD
TOPELE,P3,1.0E6
*END STEP