*HEADING
: FREQUENCY EXTRACTION OF PIEZOELECTRIC MATERIALS [CAX8E]
*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=CAX8E,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., 
*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
*SURFACE, NAME=BOT
1,S1
2,S1
*STEP,PERTURBATION
*FREQUENCY         , eigensolver=lanczos
10, , 
*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
*END STEP