*HEADING
EC3KEFDF: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [C3D20E]
*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
*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=C3D20E, ELSET=EALL
1,  1,3,9,7,19,21,27,25,2,6,8,4,20,24,26,
    22,10,12,18,16
**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=A1, ELSET=EALL,ORIENT=RECT
*MATERIAL,NAME=A1
*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
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*DLOAD
 1,BX,1000.
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,BY,2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,BZ,3000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,GRAV,100.,1.,0.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,CENT,300.,-1000.,3.5,-3.5,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,CENTRIF,30.,-1000,3.5,-3.5,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P1,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P2,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P3,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P4,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P5,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,P6,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP1,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP2,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP3,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP4,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP5,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 1,HP6,3000.,0.,-7.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,EBF,1000.
*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,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES3,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES4,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES5,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DECHARGE
 1,ES6,1000.
*END STEP
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
 ALL,40.
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S, 
*OUTPUT,FIELD
*NODE OUTPUT
 RF, 
*EL PRINT
 S, 
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 NS1,1
 NS2,2
 NS3,3
 ALL,9
*TEMPERATURE,OP=MOD
 ALL,40.
*OUTPUT,FIELD
*NODE OUTPUT
 U, 
*EL PRINT,FREQ=0
*NODE PRINT
 U, 
*END STEP