*HEADING
PPZOFRQ2: FREQUENCY EXTRACTION OF PIEZOELECTRIC MATERIALS [CAX8RE]
*NODE
 1, .005, 0.0
 5, .025, 0.0
 21, .005, .02
 25, .025, .02
*NGEN
1,5
21,25
1,21,5
2,22,5
3,23,5
4,24,5
5,25,5
*ELEMENT,TYPE=CAX8RE,ELSET=ALL
1,1,3,13,11,2,8,12,6
2,3,5,15,13,4,10,14,8
3,11,13,23,21,12,18,22,16
4,13,15,25,23,14,20,24,18
**local orientation matching global system 
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 1, 0.0
*SOLID SECTION,MATERIAL=PZT,ELSET=ALL,ORIENT=RECT
*MATERIAL,NAME=PZT
*ELASTIC,TYPE=ORTHO
13.90E10,7.428E10,11.54E10,7.784E10,7.428E10,13.90E10,2.564E10,2.564E10,
2.564E10, 
*PIEZOELECTRIC,TYPE=S
0.,0.,0.,12.71,0.,0.,-5.207,15.08
-5.207,0.,0.,0.,0.,0.,0.,0.
0.,12.71             
*DIELECTRIC,TYPE=ANISO
6.752E-9,0.,5.872E-9,0.,0.,6.752E-9
*DENSITY
5555., 
*DAMPING,ALPHA=.01,BETA=.00001
*NSET,NSET=ELECTROD,GEN
22,25
*EQUATION
2, 
ELECTROD,9,1.0,21,9,-1.0
*NSET,NSET=TOPCEN
23, 
*NSET,NSET=BOTCEN
3, 
*NSET,NSET=BOTTOM,GEN
1,5
***   Node 3 has been transformed 90 degrees, so DOF 1 is constrained
*BOUNDARY             
1,2
2,2
3,1
4,2
5,2
BOTTOM,9
*TRANSFORM,NSET=TOPCEN
0,1,0,-1,0,0
*TRANSFORM,NSET=BOTCEN
0,1,0,-1,0,0
*RESTART,WRITE
*surface, NAME=BOT
1,S1
2,S1
*STEP,PERTURBATION
*FREQUENCY         , eigensolver=lanczos, SIM=NO
45, ,,,2,15 
*EL PRINT,FREQ=0
,*NODE PRINT,tot=YES
U,EPOT
RF,
*SECTION PRINT, NAME=BOT, SURF=BOT
***SECTION FILE, NAME=BOT, SURF=BOT
**SOAREA,SOF
***MODAL FILE
*OUTPUT,FIELD,VARIABLE=PRESELECT
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*MODAL DAMPING,,RAYLEIGH
1,45,.01,.00001
*SELECT EIGENMODES, GENERATE
1,45,1
*ELSET,ELSET=TOPELE
3,4
*DLOAD
TOPELE,P3,1.0E6
*ELSET,ELSET=ELOUT
3, 
*EL PRINT
S,PHS
E,PHE
EPG,PHEPG
EFLX,PHEFL
***EL FILE,ELSET=ELOUT
**S,PHS,E,PHE,EPG,PHEPG,EFLX,PHEFL
*NODE PRINT,NSET=TOPCEN,GLO=YES
U,PU
EPOT,PHPOT
*NODE PRINT,NSET=TOPCEN,GLO=NO
U,PU
EPOT,PHPOT
*NODE PRINT,NSET=BOTCEN,GLO=NO
RF,PRF
RCHG
*NODE PRINT,NSET=BOTCEN,GLO=YES
RF,PRF
RCHG
*OUTPUT, FIELD
*NODE OUTPUT,NSET=TOPCEN
U,EPOT
*NODE OUTPUT,NSET=BOTCEN
RF,RCHG
***NODE FILE,NSET=TOPCEN,GLO=NO
**U,PU,EPOT,PHPOT
***NODE FILE,NS ET=TOPCEN,GLO=YES
**U,PU,EPOT,PHPOT
***NODE FILE,NSET=BOTCEN,GLO=NO
**RF,PRF,RCHG
***NODE FILE,NSET=BOTCEN,GLO=YES
**RF,PRF,RCHG
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DLOAD
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49, 
*DLOAD
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*MODAL DAMPING,RAYLEIGH
1,45,.01,.00001
*SELECT EIGENMODES, GENERATE
1,45,1
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49, 
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP