*HEADING
 SIMPLE SHEAR -- MULITIPLE ELEMENT TYPES 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,0.,0.
102,1.0,0.,0.
103,1.0,1.,0.
104,0.0,1.,0.
105,0.0,2.,0.
106,1.0,2.,0.
107,1.0,3.,0.
108,0.0,3.,0.
109,0.0,0.,1.
110,1.0,0.,1.
111,1.0,1.,1.
112,0.0,1.,1.
113,0.0,2.,1.
114,1.0,2.,1.
115,1.0,3.,1.
116,0.0,3.,1.
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.
301,3.0,0.
302,4.0,0.
303,4.0,1.
304,3.0,1.
305,3.0,2.
306,4.0,2.
307,4.0,3.
308,3.0,3.
401,4.5,0.
402,5.5,0.
403,5.5,1.
404,4.5,1.
405,4.5,2.
406,5.5,2.
407,5.5,3.
408,4.5,3.
501,6.0,0.
502,7.0,0.
503,7.0,1.
504,6.0,1.
505,6.0,2.
506,7.0,2.
507,7.0,3.
508,6.0,3.
*NSET,NSET=FIX
101,102,106,107, 109,110,114,115, 201,202,206,207, 301,302,306,307,
401,402,406,407, 501,502,506,507
*NSET,NSET=XSHEAR
103,104,111,112, 203,204, 303,304, 403,404, 503,504
*NSET,NSET=YSHEAR
105,108,113,116, 205,208, 305,308, 405,408, 505,508
*NSET,NSET=NOUT,GENERATE
104,504,100
105,505,100
*NSET,NSET=QA_TEST
NALL,
*ELEMENT,TYPE=C3D8R,ELSET=BRICK
1, 101,102,103,104,109,110,111,112
6, 105,106,107,108,113,114,115,116
*SOLID SECTION,ELSET=BRICK,MATERIAL=MAT
*ELEMENT,TYPE=CPE4R,ELSET=PE
2, 201,202,203,204
7, 205,206,207,208
*SOLID SECTION,ELSET=PE,MATERIAL=MAT
*ELEMENT,TYPE=CPS4R,ELSET=PS
3, 301,302,303,304
8, 305,306,307,308
*SOLID SECTION,ELSET=PS,MATERIAL=MAT
*ELEMENT,TYPE=M3D4R,ELSET=MEMBRANE
4, 401,402,403,404
9, 405,406,407,408
*MEMBRANE SECTION,ELSET=MEMBRANE,MATERIAL=MAT
.01, 
*ELEMENT,TYPE=S4R,ELSET=SHELL
 5, 501,502,503,504
10, 505,506,507,508
*SHELL SECTION,ELSET=SHELL,MATERIAL=MAT,SECTION INTEGRATION=GAUSS
.01,
*MATERIAL,NAME=MAT
*DENSITY
1.346E-4, 
*ELASTIC
1.,0.
*ELSET,ELSET=EALL
BRICK,PE,PS,MEMBRANE,SHELL
*BOUNDARY
FIX,1,3
XSHEAR,2,3
YSHEAR,1,1
YSHEAR,3,3
*RESTART,WRITE,NUMBER INTERVAL=2,TIMEMARKS=NO
*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