*HEADING
 SIMPLE SHEAR -- C3D4 ELEMENT TYPE TESTED TO 300 PERCENT SHEAR
**
**Analytical solution provided by Dienes (1979):
**S11=4*mu*(cos(2*beta)*ln(cos(beta))+beta*sin(2*beta)-sin^2(beta))
**S12=2*mu*cos(2*beta)*(2*beta-2*tan(2*beta)*ln(cos(beta)-tan(beta))
**where tan(beta)=gamma
**
*NODE,NSET=NALL
  101,0.,0.
  102,1.,0.
  103,0.,1.
  104,1.,1.
  105,0.,0.,1.
  106,1.,0.,1.
  107,0.,1.,1.
  108,1.,1.,1.
  201,0.,2.
  202,1.,2.
  203,0.,3.
  204,1.,3.
  205,0.,2.,1.
  206,1.,2.,1.
  207,0.,3.,1.
  208,1.,3.,1.
*NSET,NSET=FIX
101,102,105,106,202,206,204,208
*NSET,NSET=XSHEAR
103,104,107,108
*NSET,NSET=YSHEAR
201,203,205,207
*NSET,NSET=NOUT
103,203
*NSET,NSET=QA_TEST
NALL,
*ELEMENT,TYPE=C3D4,ELSET=EALL
1,101,103,105,102
2,107,105,103,108
3,104,102,108,103
4,106,108,102,105
5,102,103,105,108
6,201,203,205,202
7,207,205,203,208
8,204,202,208,203
9,206,208,202,205
10,202,203,205,208
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT
*MATERIAL,NAME=MAT
*DENSITY
1.346E-4, 
*ELASTIC
1.,0.
*BOUNDARY
FIX,1,3
XSHEAR,2,3
YSHEAR,1,1
YSHEAR,3,3
*RESTART,WRITE,NUMBER INTERVAL=2
*STEP
*DYNAMIC,EXPLICIT
,1.5
*BOUNDARY,TYPE=VELOCITY
XSHEAR,1,1,4.0
YSHEAR,2,2,4.0
***OUTPUT,HISTORY,TIME INTERVAL=.0015
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL,DT
***ELEMENT OUTPUT,ELSET=EALL
**S,LE
***NODE OUTPUT,NSET=NOUT
**U,
*FILE OUTPUT,NUMBER INTERVAL=1,TIMEMARKS=YES
*EL FILE
SP,LE,NE,
*NODE FILE
U, 
***********************
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,
************************ 
*END STEP