*HEADING
ECA3EFDF: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CAX3E]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1,1.,0.
 3,4.,0.
 7,1.,3.
 9,4.,3.
*ELEMENT,TYPE=CAX3E,ELSET=ELEMS
 1,1,3,9
 2,1,9,7
*SOLID SECTION,MATERIAL=A1, ELSET=ELEMS
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 3.E6,0.3
*PIEZOELECTRIC,TYPE=S
15.08,-5.207,-5.207,0.,0.,0.,0.,0.
0.,12.710,0.,0.,0.,0.,0.,0.,
0.,0.
*DIELECTRIC,TYPE=ISO
5.872E-9, 
*DENSITY
5.E-5, 
*BOUNDARY
 ALL,1,2
 ALL,9
*STEP,PERTURBATION
*STATIC
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*DLOAD
 ELEMS,BR,1000.
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,BZ,2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,CENT,3000.,,,,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,GRAV,400.,0.,-1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP1,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP2,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,HP2,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,HP3,3000.,3.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P1,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P2,2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,P2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,P3,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
ELEMS,EBF,3000.
*EL PRINT
 CHRGS,
*OUTPUT,FIELD
*NODE OUTPUT
RCHG, 
*NODE PRINT
RCHG, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
1,ES1,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
1,ES2,2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
2,ES2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
2,ES3,3000.
*END STEP