*HEADING
: FREQUENCY EXTRACTION OF PIEZOELECTRIC MATERIALS [CAX4E]
*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=CAX4E,ELSET=ALL
1,1,2,7,6
*ELGEN,ELSET=ALL
1,4,1,1,4,5,4
**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., 
*NSET,NSET=ELECTROD,GEN
22,25
*EQUATION
2, 
ELECTROD,9,1.0,21,9,-1.0
*NSET,NSET=BOTTOM,GEN
1,5
*BOUNDARY
BOTTOM,2
BOTTOM,9
*RESTART,WRITE,FREQ=999
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos
10, , 
*ELSET,ELSET=TOPELE_O
13, 
*EL PRINT,MODE=1,ELSET=TOPELE_O
S,E
EPG,EFLX
*EL FILE,LASTMODE=1,ELSET=TOPELE_O
S,E,EFLX,EPG
*NSET,NSET=TOPCEN
23, 
*NODE PRINT,NSET=TOPCEN,MODE=2
U,EPOT
*NODE FILE,NSET=TOPCEN,MODE=2
U,EPOT
*NSET,NSET=B3
3, 
*NODE FILE,NSET=B3,LASTMODE=1
RF,RCHG
*END STEP