*HEADING
 SIMPLE SHEAR -- M3D3 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
201,1.5,0.
202,2.5,0.
203,2.5,1.
204,1.5,1.
205,1.5,2.
206,2.5,2.
207,2.5,3.
208,1.5,3.
*NSET,NSET=FIX
201,202,206,207
*NSET,NSET=XSHEAR
203,204
*NSET,NSET=YSHEAR
205,208
*NSET,NSET=NOUT
204,205
*NSET,NSET=QA_TEST_NALL
FIX,XSHEAR,YSHEAR,NOUT
*ELEMENT,TYPE=M3D3,ELSET=EALL
1, 201,202,203
2, 201,203,204
3, 205,206,207
4, 205,207,208
*ELSET,ELSET=QA_TEST_EALL
EALL,
*MEMBRANE SECTION,ELSET=EALL,MATERIAL=MAT, POISSON=0.
0.01, 
*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
LE,
*NODE FILE
U, 
*****************************************
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*ELEMENT OUTPUT
S,
*ELEMENT OUTPUT,ELSET=QA_TEST_EALL
LE,
*NODE OUTPUT, NSET=QA_TEST_NALL
U, 
*OUTPUT,HISTORY,TIME INTERVAL=1.5
*ENERGY OUTPUT, VAR=ALL
*****************************************
*END STEP