*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER 
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, COARSE MODEL WITH
FOUR-NODE AXISYMMETRIC PIEZOELECTRIC ELEMENTS
*NODE
 1, .0110,.0
 2, .01175, 0.0
 3, .0125, 0.0
 13, .0110, .0125
 14, .01175,.0125
 15, .0125,.0125
 16, .0110, .0155
 17, .01175,.0155
 18, .0125,.0155
 19, 0.0,.0125
 21,.0073333,.0125
 22, 0.0,.0155
 24,.0073333,.0155
*NGEN
1,13,3
2,14,3
3,15,3
19,21
22,24
*ELEMENT,TYPE=CAX4E
 1, 1,2,5,4
*ELEMENT,TYPE=CAX4,ELSET=BRASS
 9, 13,14,17,16
10, 14,15,18,17
11, 19,20,23,22
12, 20,21,24,23
13, 21,13,16,24
*ELGEN,ELSET=PZT4
1,2,1,1,4,3,2
**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=PZT4,ELSET=PZT4,ORIENT=RECT
*SOLID SECTION,MATERIAL=BRASS,ELSET=BRASS
*MATERIAL,NAME=PZT4
*ELASTIC,TYPE=ORTHOTROPIC
11.54E10,7.428E10,13.90E10,7.428E10,7.784E10,13.90E10,2.564E10,2.564E10
2.564E10, 
*PIEZOELECTRIC,TYPE=S
15.08,-5.207,-5.207,0.,0.,0.,0.,0.
0.,12.710,0.,0.,0.,0.,0.,0.,
12.710,0.
*DIELECTRIC,TYPE=ANISO
5.872E-9,0.,6.752E-9,0.,0.,6.752E-9
*DENSITY
7500., 
*MATERIAL,NAME=BRASS
*ELASTIC,TYPE=ISOTROPIC
 10.4E10,0.37
*DENSITY
8500., 
*NSET,NSET=IN,GENERATE
1,13,3
*NSET,NSET=SYM
1,2,3
*NSET,NSET=ASYM
19,22
*BOUNDARY
SYM,2
ASYM,1
IN,9
*RESTART,WRITE
**
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos, SIM=NO
5, ,100000.,10. 
*EL PRINT
S,E
EFLX,EPG
*NODE PRINT
U,EPOT
RF,RCHG
*ELSET,ELSET=OUTPUT
1,9
*EL FILE,ELSET=OUTPUT
S,E,EFLX,EPG
*NSET,NSET=OUTPUT1
5,8
*NODE FILE,NSET=OUTPUT1
U,EPOT
*NSET,NSET=OUTPUT2
1,19
*NODE FILE,NSET=OUTPUT2
RF,RCHG
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=OUTPUT
S,E,EFLX,EPG
*OUTPUT,FIELD
*NODE OUTPUT
U,
*NODE OUTPUT,NSET=OUTPUT1
EPOT, 
*OUTPUT,FIELD
*NODE OUTPUT,NSET=OUTPUT2
RF,RCHG
*END STEP