type('119.inp')
*HEADING
 VERIFICATION OF ABAQUS ELECTRICAL REACTION CHARGE OUTPUT TO MATLAB (RCHG, RECORD KEY 119)
*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
*ELEMENT,TYPE=C3D8E, ELSET=EALL
 1, 1,2,3,4,5,6,7,8
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*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, 
*BOUNDARY
 1,1,3
 2,2
 4,3
 5,1
 1,9
 2,9
 3,9
 4,9
*SURFACE, NAME=SIDE1485
 1,S6
*STEP,PERTURBATION
*STATIC
*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,  1500.00
 3,1,   500.00
 3,2,   500.00
 3,3, -1000.00
 4,1,   500.00
 4,2,  1500.00
 5,2,  -500.00
 5,3,  1000.00
 6,1,  -500.00
 6,2, -1500.00
 7,1, -1500.00
 7,2, -1500.00
 7,3, -1000.00
 8,1, -1500.00
 8,2,  -500.00
*CECHARGE
 1,,-2000.
 2,,-2000.
 3,,-2000.
 4,,-2000.
 5,,-1000.
 6,,-1000.
 7,,-1000.
 8,,-1000.
*FILE FORMAT, ASCII
*NODE FILE
 RCHG
*END STEP