*HEADING
Test of angular acceleration loading for C3D10E elements.
*RESTART,WRITE
*NODE,NSET=NALL
 1,9.,-1.,1.
 3,11.,-1.,1.
 19,9.,1.,1.
 21,11.,1.,1.
 7,9.,-1.,-1.
 9,11.,-1.,-1.
 25,9.,1.,-1.
 27,11.,1.,-1.
*NGEN,NSET=NALL
 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
*ELEMENT,TYPE=C3D10E,ELSET=EALL
 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=STEEL,ELSET=EALL,ORIENT=RECT
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,.33
*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
*DENSITY
2000., 
*BOUNDARY
NALL,1,3
*STEP,PERTURBATION
*STATIC
*DLOAD
EALL,ROTA,10.,0.,0.,0.,0,0,1
*NODE PRINT,TOTALS=YES
RF, 
*EL PRINT
EVOL, 
*EL PRINT,FREQ=0
*OUTPUT,FIELD
*NODE OUTPUT
RF, 
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT
*END STEP