*HEADING
EC3AEFDF: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [C3D10E]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
 1,0.,0.,0.
 3,7.,0.,0.
 19,0.,7.,0.
 21,7.,7.,0.
 7,0.,0.,-7.
 9,7.,0.,-7.
 25,0.,7.,-7.
 27,7.,7.,-7.
*NGEN,NSET=ALL
 1,3
 1,7,3
 3,9,3
 7,9
 1,19,9
 3,21,9
 7,25,9
 9,27,9
 19,25,3
 21,27,3
 19,21
 25,27
 8,26,9
 2,20,9
 2,8,3
 10,16,3
 20,26,3
 12,18,3
*NSET,NSET=NS1
1    , 4    , 7    , 16   , 25   , 22   , 19   , 10
*NSET,NSET=NS2
1    , 2    , 3    , 6    , 9    , 8    , 7    , 4
*NSET,NSET=NS3
1    , 2    , 3    , 12   , 21   , 20   , 19   , 10
*ELEMENT,TYPE=C3D10E,ELSET=ELEMS
 1, 7,9,25,1, 8,17,16, 4,5,13
 2, 1,9,21,3, 5,15,11, 2,6,12
 3, 9,27,25,21, 18,26,17, 15,24,23
 4, 1,21,25,19, 11,23,13, 10,20,22
 5, 1,21,9,25, 11,15,5, 13,23,17
**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=ELMAT, ELSET=ELEMS,ORIENT=RECT
*MATERIAL,NAME=ELMAT
*ELASTIC
 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.,
 12.710,0.
*DIELECTRIC,TYPE=ANISO
 5.872E-9,0.,6.752E-9,0.,0.,6.752E-9
*EXPANSION
 .0001, 
*DENSITY
 10.0, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL,-10.
*BOUNDARY
 ALL,1,3
 ALL,9
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,BX,3000.
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
*OUTPUT,FIELD
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=ALL
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,BY,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,BZ,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,GRAV,300.,1.,0.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,CENT,300.,-1000.,3.5,-3.5,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,CENTRIF,30.,-1000,3.5,-3.5,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P1,3000.
 3,P1,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,P4,3000.
 4,P2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P2,3000.
 2,P2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,P3,3000.
 3,P2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 3,P3,3000.
 4,P3,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P4,3000.
 4,P4,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP1,3000.,0.,-7.
 3,HP1,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,HP4,3000.,0.,-7.
 4,HP2,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP2,3000.,0.,-7.
 2,HP2,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 2,HP3,3000.,0.,-7.
 3,HP2,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 3,HP3,3000.,0.,-7.
 4,HP3,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP4,3000.,0.,-7.
 4,HP4,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
ELEMS,EBF,3000.
*EL PRINT
 CHRGS,
*NODE PRINT
RCHG, 
*OUTPUT,FIELD
*NODE OUTPUT
RCHG, 
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ALL
RCHG, 
*ENERGY OUTPUT,VAR=ALL
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES1,3000.
 3,ES1,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 2,ES4,3000.
 4,ES2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES2,3000.
 2,ES2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 2,ES3,3000.
 3,ES2,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 3,ES3,3000.
 4,ES3,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES4,3000.
 4,ES4,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,40.
*EL PRINT
 S, 
*NODE PRINT
 RF, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S, 
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=ALL
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 NS1,1
 NS2,2
 NS3,3
 ALL,9
*TEMPERATURE,OP=MOD
 ALL,40.
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*OUTPUT,FIELD
*NODE OUTPUT
U, 
*ELEMENT OUTPUT
S, 
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=ALL
*END STEP