*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER          
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, 3-D MODEL WITH
TWENTY-NODE BRICK ELEMENTS
ONE SECTION WITH AXISYMMETRIC BOUNDARY CONDITIONS
*NODE,NSET=START,SYSTEM=C
 1, 0.0001,0.,.0155
 5, .0110,0., .0155
 9, .0125,0.,.0155
 41, 0.0001,0.,.0125
 45, .0110,0., .0125
 49, .0125,0.,.0125
 125, .0110,0.,0.0
 129, .0125,0.,0.0
*NODE,NSET=END,SYSTEM=C
 1001, 0.0001,10.,.0155
 1005, .0110,10., .0155
 1009, .0125,10.,.0155
 1041, 0.0001,10.,.0125
 1045, .0110,10., .0125
 1049, .0125,10.,.0125
 1125, .0110,10.,0.0
 1129, .0125,10.,0.0
*NGEN,NSET=L1
1,5
5,9
*NGEN,NSET=L2
41,45
45,49
*NFILL,NSET=START
L1,L2,4,10
*NGEN,NSET=L3
1001,1005
1005,1009
*NGEN,NSET=L4
1041,1045
1045,1049
*NFILL,NSET=END
L3,L4,4,10
*NSET,NSET=L5,GENERATE
45,49
*NGEN,NSET=L6
125,129
*NFILL,NSET=START
L5,L6,8,10
*NSET,NSET=L7,GENERATE
1045,1049
*NGEN,NSET=L8
1125,1129
*NFILL,NSET=END
L7,L8,8,10
*NFILL,NSET=ALL
START,END,2,500
*ELEMENT,TYPE=C3D20E
 1, 1,3,23,21,1001,1003,1023,1021,2,13,22,11,1002,1013,1022,
    1011,501,503,523,521
*ELGEN
1,4,2,1,2,20,4
7,2,20,2
9,4,20,1
8,2,20,5
13,4,20,1
*ELSET,ELSET=PZT4,GENERATE
9,16
*ELSET,ELSET=BRASS,GENERATE
1,8
*SOLID SECTION,MATERIAL=PZT4,ELSET=PZT4,ORIENTATION=PZTDIR
*SOLID SECTION,MATERIAL=BRASS,ELSET=BRASS,ORIENTATION=PZTDIR
*ORIENTATION,NAME=PZTDIR,SYSTEM=USER
*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., 
*PIEZOELECTRIC,TYPE=S
0.,0.,0.,0.,0.,0.,0.,0.
0.,0.,0.,0.,0.,0.,0.,0.,
0.,0.
*NSET,NSET=MID,GENERATE
501,629
*NSET,NSET=IN,GENERATE
45,125,10
545,625,10
1045,1125,10
*NSET,NSET=SYM,GENERATE
125,129
625,629
1125,1129
*NSET,NSET=ASYM,GENERATE
1,41,10
501,541,10
1001,1041,10
*TRANSFORM, NSET=END
.98481,.17365,0.,-.17365,.98481,0.
*TRANSFORM,NSET=MID
.99619,.08716,0.,-.08716,.99619,0.
*BOUNDARY
SYM,3
ASYM,1,2
IN,9
START,2
END,2
MID,2
*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
107,549
*NODE FILE,NSET=OUTPUT1
U,EPOT
*NSET,NSET=OUTPUT2
1,45
*NODE FILE,NSET=OUTPUT2
RF,RCHG
*OUTPUT,FIELD,VARIABLE=PRESELECT,OP=NEW,FREQUENCY=9999
*OUTPUT,FIELD,OP=NEW
*ELEMENT OUTPUT,ELSET=OUTPUT
S,E,EFLX,EPG
*OUTPUT,HISTORY,OP=NEW
*ELEMENT OUTPUT,ELSET=OUTPUT
S,E,EFLX,EPG
*OUTPUT,FIELD,OP=NEW
*NODE OUTPUT,NSET=OUTPUT1
U,EPOT
*OUTPUT,HISTORY,OP=NEW
*NODE OUTPUT,NSET=OUTPUT1
U,EPOT
*OUTPUT,FIELD,OP=NEW
*NODE OUTPUT,NSET=OUTPUT2
RF,RCHG
*OUTPUT,HISTORY,OP=NEW
*NODE OUTPUT,NSET=OUTPUT2
RF,RCHG
*END STEP