*HEADING
 SIMPLE SHEAR -- C3D8R 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=ODB_NSET
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.
*NSET,NSET=FIX
101,102,106,107, 109,110,114,115
*NSET,NSET=XSHEAR
103,104,111,112
*NSET,NSET=YSHEAR
105,108,113,116
*NSET,NSET=NOUT
104,105
*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
*MATERIAL,NAME=MAT
*DENSITY
1.346E-4, 
*ELASTIC
1.,0.
*ELSET,ELSET=EALL
BRICK, 
*BOUNDARY
FIX,1,3
XSHEAR,2,3
YSHEAR,1,1
YSHEAR,3,3
*RESTART,WRITE,NUMBER INTERVAL=2,TIMEMARKS=NO
**
** Sets added by addodb
**
*Elset,Elset=ODB_ESET
BRICK,
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, FIX, XSHEAR, YSHEAR
NOUT, 
**
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
BRICK,
**
*NSET, NSET=QA_TEST_ALLNODES
ODB_NSET, FIX, XSHEAR, YSHEAR
NOUT, ODB_NSET
**
*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=NO
*EL FILE
SP,LE,NE,
*NODE FILE
U, 
*OUTPUT, FIELD, OP=ADD, NUMBER INTERVAL=2, TIMEMARKS=NO
*ELEMENT OUTPUT
ELEN,
S,SP,LE,LEP,NE,NEP,PRESS,
*NODE OUTPUT
U, 
*OUTPUT, HISTORY, TIME INTERVAL=0.75
*ELEMENT OUTPUT, Elset=ODB_ESET
ELEN,
S,SP,LE,LEP,NE,NEP,PRESS,
*NODE OUTPUT, Nset=ODB_NSET
U,
*ENERGY OUTPUT,VARIABLE=ALL
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=NO
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS
SP,LE,NE,
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U, 
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=1.5
*END STEP