*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER 
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, REFINED MODEL WITH
EIGHT-NODE AXISYMMETRIC PIEZOELECTRIC ELEMENTS
*NODE
 17, .0110,.0
 25, .0125, 0.0
 1617, .0110,.0125
 1625, .0125,.0125
 2417, .0110, .0155
 2425, .0125,.0155 
 1601, 0.0,.0125
 2401, 0.0,.0155
*NGEN,NSET=SYM
17,25
*NGEN,NSET=JUNCTION
1617,1625 
*NGEN,NSET=TOPB
2401,2417
*NGEN,NSET=TOPP
2417,2425
*NGEN,NSET=ASYM
1601,2401,100
*NGEN,NSET=RIGHT
1617,2417,100
*NFILL
SYM,JUNCTION,16,100
JUNCTION,TOPP,8,100
ASYM,RIGHT,16,1
*ELEMENT,TYPE=CAX8RE
 1, 17,19,219,217,18,119,218,117
*ELEMENT,TYPE=CAX8R
 1001, 1601,1603,1803,1801,1602,1703,1802,1701
*ELGEN,ELSET=PZT4
1,4,2,1,8,200,4
*ELGEN,ELSET=BRASS
1001,12,2,1,4,200,12
**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
13.90E10,7.428E10,11.54E10,7.784E10,7.428E10,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
17,1617,100
*BOUNDARY
SYM,2
ASYM,1
IN,9
*RESTART,WRITE
**
**
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos, SIM=NO
5, ,100000.,10. 
*EL PRINT,FREQUENCY=0
*NODE PRINT
U,EPOT 
RF,RCHG
*ELSET,ELSET=OUTPUT
1,1001
*EL FILE,ELSET=OUTPUT
S,E,EFLX,EPG
*NSET,NSET=OUTPUT1
218,1819
*NODE FILE,NSET=OUTPUT1
U,EPOT
*NSET,NSET=OUTPUT2
17,1601
*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