*HEADING
STEADY STATE DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX3E]
INCLUDES RAYLEIGH DAMPING
*NODE
 1, .005, 0.0
 5, .025, 0.0
 21, .005, .02
 25, .025, .02
*NGEN
1,5
21,25
1,21,5
2,22,5
3,23,5
4,24,5
5,25,5
*ELEMENT,TYPE=CAX3E,ELSET=ALL
1,1,2,6
2,2,7,6
*ELGEN,ELSET=ALL
1,4,1,2,4,5,8
2,4,1,2,4,5,8
**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=PZT,ELSET=ALL,ORIENT=RECT
*MATERIAL,NAME=PZT
*ELASTIC,TYPE=ORTHO
13.90E10,7.428E10,11.54E10,7.784E10,7.428E10,13.90E10,2.564E10,2.564E10,
2.564E10, 
*PIEZOELECTRIC,TYPE=S
0.,0.,0.,12.71,0.,0.,-5.207,15.08
-5.207,0.,0.,0.,0.,0.,0.,0.
0.,12.71
*DIELECTRIC,TYPE=ANISO
6.752E-9,0.,5.872E-9,0.,0.,6.752E-9
*DENSITY
5555., 
*DAMPING,ALPHA=.01,BETA=.00001
*NSET,NSET=ELECTROD,GEN
22,25
*EQUATION
2, 
ELECTROD,9,1.0,21,9,-1.0
*NSET,NSET=TOPCEN
23, 
*NSET,NSET=BOTCEN
3, 
*NSET,NSET=BOTTOM,GEN
1,5
***   Node 3 has been transformed 90 degrees, so DOF 1 is constrained
*BOUNDARY
1,2
2,2
3,1
4,2
5,2
BOTTOM,9
*TRANSFORM,NSET=TOPCEN
0,1,0,-1,0,0
*TRANSFORM,NSET=BOTCEN
0,1,0,-1,0,0
*RESTART,WRITE,FREQ=1
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos, SIM=NO
45, , 
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*MODAL DAMPING,RAYLEIGH
1,45,.01,.00001
*ELSET,ELSET=TOPELE,GEN
26,32,2
*DLOAD
TOPELE,P2,1.0E6
*ELSET,ELSET=ELOUT
1,9,17,25
*EL PRINT
S,PHS
E,PHE
EPG,PHEPG
EFLX,PHEFL
*EL FILE,ELSET=ELOUT
S,PHS,E,PHE,EPG,PHEPG,EFLX,PHEFL
*NODE PRINT,NSET=TOPCEN,GLO=YES
U,PU
EPOT,PHPOT
*NODE PRINT,NSET=TOPCEN,GLO=NO
U,PU
EPOT,PHPOT
*NODE PRINT,NSET=BOTCEN,GLO=NO
RF,PRF
RCHG,
*NODE PRINT,NSET=BOTCEN,GLO=YES
RF,PRF
RCHG,
*NODE FILE,NSET=TOPCEN,GLO=NO
U,PU,EPOT,PHPOT
*NODE FILE,NSET=TOPCEN,GLO=YES
U,PU,EPOT,PHPOT
*NODE FILE,NSET=BOTCEN,GLO=NO
RF,PRF,RCHG
*NODE FILE,NSET=BOTCEN,GLO=YES
RF,PRF,RCHG
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DLOAD
TOPELE,P2,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49, 
*DLOAD
TOPELE,P2,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49, 
*MODAL DAMPING,RAYLEIGH
1,45,.01,.00001
*SELECT EIGENMODES, GENERATE
1,45,1
*DLOAD,LOADCASE=2
TOPELE,P2,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49, 
*DLOAD,LOADCASE=2
TOPELE,P2,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49, 
*DLOAD,LOADCASE=2
TOPELE,P2,1.0E6
*END STEP