*HEADING
: STEADY STATE DYNAMIC ANALYSIS OF PIEZOELECTRIC MATERIALS [CAX4E]
INCLUDES RAYLEIGH DAMPING
*NODE
*NGEN
*ELEMENT,TYPE=CAX4E,ELSET=ALL
*ELGEN,ELSET=ALL
**local orientation matching global system
*ORIENTATION,NAME=RECT
*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
*DIELECTRIC,TYPE=ANISO
*DENSITY
5555.,
*DAMPING,ALPHA=.01,BETA=.00001
*NSET,NSET=ELECTROD,GEN
*EQUATION
ELECTROD,9,1.0,21,9,-1.0
*NSET,NSET=TOPCEN
*NSET,NSET=BOTCEN
*NSET,NSET=BOTTOM,GEN
***   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
*TRANSFORM,NSET=BOTCEN
*RESTART,WRITE,FREQ=999
*STEP,PERTURBATION
*FREQUENCY, eigensolver=lanczos
45, ,
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49,
*MODAL DAMPING,RAYLEIGH
*ELSET,ELSET=TOPELE,GEN
*DLOAD
TOPELE,P3,1.0E6
*ELSET,ELSET=ELOUT
*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,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49,
*DLOAD
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS
15915.49,
*MODAL DAMPING,RAYLEIGH
*SELECT EIGENMODES, GENERATE
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
15915.49,
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
15915.49,
*DLOAD,LOADCASE=2
TOPELE,P3,1.0E6
*END STEP