*HEADING
: DLOAD, DECHARGE, AND TEMPERATURE FOR PIEZO. ELEMENTS  [CPS6E]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1    , 0.        , 0.
3    , 7.        , 0.
7    , 0.        , 7.
9    , 7.        , 7.
*NGEN,NSET=BOTTOM
1    , 3
*NGEN,NSET=TOP
7    , 9
*NGEN,NSET=LEFT
1    , 7    , 3
*NGEN,NSET=RIGHT
3    , 9    , 3
*NGEN
4    , 6
*NSET,NSET=ALL,GENERATE
1    , 9
*NSET,NSET=SIDES
 1,2,3,4,5,6,7,8,9
*ELEMENT,TYPE=CPS6E,ELSET=ELEMS
1    , 1    , 3    , 9    , 2    , 6    , 5
2    , 1    , 9    , 7    , 5    , 8    , 4
*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.,
12.710,0.
*DIELECTRIC,TYPE=ISO
5.872E-9, 
*EXPANSION
0.0001, 
*DENSITY
5.E-5, 
*INITIAL CONDITION,TYPE=TEMPERATURE
ALL  , -10.
*BOUNDARY
SIDES,1,2
SIDES,9
*STEP,PERTURBATION
*STATIC
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*NODE FILE
 RF, 
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*DLOAD
ELEMS,BX,1000.
*EL PRINT
 LOADS,
*EL FILE
 LOADS,
*NODE PRINT
 RF, 
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
ELEMS,BY,2000.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
ELEMS,CENT,3000.,0.,0.,0.,0.,1.,0.
*END STEP
*STEP,PERTURBATION
*STATIC
*DLOAD
 ELEMS,CENTRIF,6.E7,0.,0.,0.,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,1000.
*EL PRINT
 CHRGS,
*EL FILE
 CHRGS,
*NODE FILE
RCHG, 
*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
*STEP,PERTURBATION
*STATIC
*TEMPERATURE
ALL,40.
*EL FILE
 S, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*NODE FILE
 RF,U
*OUTPUT,FIELD
*NODE OUTPUT
RF,U
*EL PRINT
 S, 
*NODE PRINT
 RF,U
*END STEP
*STEP,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
LEFT,1
BOTTOM,2
ALL,9
*TEMPERATURE,OP=MOD
ALL,40.
*END STEP