*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER 
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, COARSE MODEL WITH
EIGHT-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, .014
17, .01175,.014
18, .0125,.014
19, .0110, .0155
20, .01175,.0155
21, .0125,.0155
22, 0.0,.0125
25,.00825,.0125
26, 0.0,.014
27, .0055,.014
28, 0.0,.0155
31,.00825,.0155
*NGEN
1,13,3
2,14,6
3,15,3
22,25
28,31
*ELEMENT,TYPE=CAX8RE
1, 1,3,9,7,2,6,8,4
2, 7,9,15,13,8,12,14,10
*ELEMENT,TYPE=CAX8R
3, 13,15,21,19,14,18,20,16
4, 22,24,30,28,23,27,29,26
5, 24,13,19,30,25,16,31,27
*ELSET,ELSET=PZT4
1,2
*ELSET,ELSET=BRASS
3,4,5
**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
22,26,28
*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,4
*EL FILE,ELSET=OUTPUT
S,E,EFLX,EPG
*NSET,NSET=OUTPUT1
8,14
*NODE FILE,NSET=OUTPUT1
U,EPOT
*NSET,NSET=OUTPUT2
1,26
*NODE FILE,NSET=OUTPUT2
RF,RCHG
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=OUTPUT
S,E,EFLX,EPG
*OUTPUT,FIELD
*NODE OUTPUT,NSET=OUTPUT1
U,EPOT
*OUTPUT,FIELD
*NODE OUTPUT
U, 
*NODE OUTPUT,NSET=OUTPUT2
RF,RCHG
*END STEP