*HEADING
:  STATIC ANALYSIS OF PIEZOELECTRIC MATERIALS   [CAX4E]
*NODE
 1, .005, 0.0
 5, .025, 0.0
 21, .005, .02
 25, .025, .02
*ELEMENT,TYPE=CAX4E,ELSET=ALL
1,1,5,25,21
**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.710
*DIELECTRIC,TYPE=ANISO
6.752E-9,0.,5.872E-9,0.,0.,6.752E-9
*DENSITY
5555., 
** CREATE NSET FOR DEPENDENT POTENTIALS ON TOP ELECTRODE
*NSET,NSET=ELECTDEP
25, 
** CONSTRAIN ALL POTENTIALS ON TOP ELECTRODE TO HAVE SAME VALUES
*EQUATION
2, 
ELECTDEP,9,1.0,21,9,-1.0
*NSET,NSET=BOTTOM
1,5
*BOUNDARY
BOTTOM,2
** BOTTOM ELECTRODE ALL HAVE ZERO POTENTIAL
BOTTOM,9
*RESTART,WRITE,FREQ=999
*STEP,PERTURBATION
Step with distributed mechanical load
*STATIC
*ELSET,ELSET=TOP
1, 
*DLOAD
TOP,P3,1.0E6
*EL PRINT
S,E
EFLX,EPG
LOADS,CHRGS
*NODE PRINT
U,EPOT
RF,RCHG
*ELSET,ELSET=ELOUT
1, 
*EL FILE,ELSET=ELOUT
S,E,EFLX,EPG
LOADS,CHRGS
*NSET,NSET=NOUT
1,5,21,25
*NODE FILE,NSET=NOUT
U,EPOT
RF,RCHG
*ENERGY PRINT
*ENERGY FILE
*END STEP
*STEP,PERTURBATION
Step with distributed electrical charge load
*STATIC
*DECHARGE
TOP,ES3,1.0E-3
*END STEP
*STEP,PERTURBATION
Prescribed displacements and potentials, approximates solution to Step 1
*STATIC
*BOUNDARY
1,1,1,1.11E-08
5,1,1,5.55E-08
21,1,1,1.11E-08
25,1,1,5.55E-08
21,2,2,-1.65E-07
25,2,2,-1.65E-07
21,9,,-502
25,9,,-502
*END STEP
*STEP,PERTURBATION
Prescribed displacements and potentials, approximates solution to Step 2
*STATIC
*BOUNDARY
1,1,1,5.34E-08
5,1,1,2.67E-07
21,1,1,5.34E-08
25,1,1,2.67E-07
21,2,2,-5.02E-07
25,2,2,-5.02E-07
21,9,,1737
25,9,,1737
*END STEP