*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER      
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, 3-D MODEL WITH
EIGHT-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,2
5,9,2
*NGEN,NSET=L2
41,45,2
45,49,2
*NFILL,NSET=START
L1,L2,4,10
*NGEN,NSET=L3
1001,1005,2
1005,1009,2
*NGEN,NSET=L4
1041,1045,2
1045,1049,2
*NFILL,NSET=END
L3,L4,4,10
*NSET,NSET=L5,GENERATE
45,49,2
*NGEN,NSET=L6
125,129,2
*NFILL,NSET=START
L5,L6,8,10
*NSET,NSET=L7,GENERATE
1045,1049,2
*NGEN,NSET=L8
1125,1129,2
*NFILL,NSET=END
L7,L8,8,10
*NFILL,NSET=ALL
START,END,1,1000
*ELEMENT,TYPE=C3D8
1,1,3,23,21,1001,1003,1023,1021
*ELEMENT,TYPE=C3D8E
9,45,47,67,65,1045,1047,1067,1065
*ELGEN
1,4,2,1,2,20,4
9,4,20,1
9,2,2,4
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., 
*NSET,NSET=IN,GENERATE
45,125,10
1045,1125,10
*NSET,NSET=SYM,GENERATE
125,129
1125,1129
*NSET,NSET=ASYM,GENERATE
1,41,10
1001,1041,10
*TRANSFORM, NSET=END
.98481,.17365,0.,-.17365,.98481,0.
*BOUNDARY
SYM,3
ASYM,1,2
IN,9
START,2
END,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,13
*EL FILE,ELSET=OUTPUT
S,E,EFLX,EPG
*NSET,NSET=OUTPUT1
1047,1049
*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,FIELD,OP=NEW
*NODE OUTPUT,NSET=OUTPUT1
U,EPOT
*OUTPUT,FIELD,OP=NEW
*NODE OUTPUT,NSET=OUTPUT2
RF,RCHG
*END STEP