*HEADING
EC3KERS2: PRESSURE AND HYDROSTATIC PRESSURE FOR 3-D SOLIDS  [C3D20RE]
*RESTART,WRITE
*NODE
  1, 0.0, 0.0, 0.0
  2, 2.0, 0.0, 0.0
  3, 2.0, 2.0, 0.0
  4, 0.0, 2.0, 0.0
  5, 0.0, 0.0, 1.0
  6, 2.0, 0.0, 1.0
  7, 2.0, 2.0, 1.0
  8, 0.0, 2.0, 1.0
  9, 1.0, 0.0, 0.0
 10, 2.0, 1.0, 0.0
 11, 1.0, 2.0, 0.0
 12, 0.0, 1.0, 0.0
 13, 1.0, 0.0, 1.0
 14, 2.0, 1.0, 1.0
 15, 1.0, 2.0, 1.0
 16, 0.0, 1.0, 1.0
 17, 0.0, 0.0, 0.5
 18, 2.0, 0.0, 0.5
 19, 2.0, 2.0, 0.5
 20, 0.0, 2.0, 0.5
*ELEMENT,TYPE=C3D20RE, ELSET=EALL
 1, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
 16,17,18,19,20
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*surface,NAME=SIDE
1,S2
*MATERIAL,NAME=A1
*ELASTIC,TYPE=ISOTROPIC
 30.0E6,0.3
*PIEZOELECTRIC,TYPE=S
0.,0.,0.,0.,0.,0.,0.,0.
0.,0.,0.,0.,0.,0.,0.,0.
0.,0.
*DIELECTRIC,TYPE=ISO
1.0E-3, 
*STEP,PERTURBATION
*STATIC
*BOUNDARY
 1,1,3
 2,2
 4,3
 5,1
 1,9
 2,9
 3,9
 4,9
 9,9
10,9
11,9
12,9
**ADDITIONAL B.C.'S ADDED TO AVOID SINGULARITIES
 9,2
12,3
17,1
20,1,,-1.73334E-04
20,2,,-6.99995E-05
20,3,,-6.6665E-06
*DLOAD
 1, P1, 1000.0
 1, P2, 1000.0
 1, P3, 1000.0
 1, P4, 1000.0
 1, P5, 1000.0
 1, P6, 1000.0
*DECHARGE
 1, ES1, 1000.0
 1, ES2, 1000.0
 1, ES3, 1000.0
 1, ES4, 1000.0
 1, ES5, 1000.0
 1, ES6, 1000.0
*CLOAD
 2,1,  -500.
 3,1,  -166.6667
 3,2,  -166.6667
 3,3,   333.33333
 4,1,  -166.6667
 4,2,  -500.0
 5,2,   166.666
 5,3,  -333.33333
 6,1,   166.6667
 6,2,   500.0
 7,1,   500.
 7,2,   500.
 7,3,   333.33333
 8,1,   500.0
 8,2,   166.6667
 9,1,   2000.
 9,2,   1333.33333
 9,3,   666.6667
 10,1,  1333.33333
 10,2,  666.6667
 10,3, -666.6667
 11,1,  666.6667
 11,2,  1333.33333
 11,3, -666.6667
 12,1,  1333.33333
 12,2,  2000.
 12,3,  666.6667
 13,1, -666.6667
 13,2, -1333.33333
 13,3,  666.6667
 14,1, -1333.33333
 14,2, -2000.
 14,3, -666.6667
 15,1, -2000.
 15,2, -1333.33333
 15,3, -666.6667
 16,1, -1333.33333
 16,2, -666.6667
 16,3,  666.66667
 17,1,  666.66667
 17,2,  666.66667
 17,3,  1333.33333
 18,1,  666.66667
 18,2, -666.66667
 19,1, -666.66667
 19,2, -666.66667
 19,3, -1333.33333
 20,1, -666.66667
 20,2,  666.66667
*CECHARGE
1,,666.667
2,,666.667
3,,666.667
4,,666.667
9,,-2000.
10,,-2000.
11,,-2000.
12,,-2000.
5,,333.333
6,,333.333
7,,333.333
8,,333.333
13,,-666.667
14,,-666.667
15,,-666.667
16,,-666.667
17,,-1333.333
18,,-1333.333
19,,-1333.333
20,,-1333.333
*EL PRINT
NFORC,
*SECTION PRINT,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
 ,1.,1,1.
 ,0.,1.,1., ,1.,0.,1
SOAREA,SOF
SOCF,
SOE, 
*SECTION FILE,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL
 ,1.,1,1.
 ,0.,1.,1., ,1.,0.,1.
SOAREA,SOF
SOCF,
SOE, 
*EL FILE
 S,COORD,E,
 LOADS,
 CHRGS,
 EFLX,EPG
 ENER,
 ELEN,
*ENERGY PRINT
*ENERGY FILE
*NODE FILE
 U,COORD,EPOT,RCHG,RF
*NODE PRINT
 U,COORD,EPOT
 RCHG,RF
*EL PRINT
 S,COORD,
 E,
 LOADS,
 CHRGS,
 EFLX,EPG
 ENER,
 ELEN,
*END STEP