*HEADING
EIGENVALUE ANALYSIS OF CYLINDRICAL TRANSDUCER
PIEZOELECTRIC CERAMIC AND BRASS END CAPS, REFINED MODEL WITH
EIGHT-NODE AXISYMMETRIC PIEZOELECTRIC ELEMENTS
*NODE
*NGEN,NSET=SYM
*NGEN,NSET=JUNCTION
*NGEN,NSET=TOPB
*NGEN,NSET=TOPP
*NGEN,NSET=ASYM
*NGEN,NSET=RIGHT
*NFILL
SYM,JUNCTION,16,100
JUNCTION,TOPP,8,100
ASYM,RIGHT,16,1
*ELEMENT,TYPE=CAX8RE
*ELEMENT,TYPE=CAX8R
*ELGEN,ELSET=PZT4
*ELGEN,ELSET=BRASS
**local orientation matching global system
*ORIENTATION,NAME=RECT
*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
*DIELECTRIC,TYPE=ANISO
*DENSITY
7500.,
*MATERIAL,NAME=BRASS
*ELASTIC,TYPE=ISOTROPIC
 10.4E10,0.37
*DENSITY
8500.,
*NSET,NSET=IN,GENERATE
*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
*EL FILE,ELSET=OUTPUT
S,E,EFLX,EPG
*NSET,NSET=OUTPUT1
*NODE FILE,NSET=OUTPUT1
U,EPOT
*NSET,NSET=OUTPUT2
*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