** TEST OF REBAR: UNIAXIAL TENSION LOADING
**                ELASTIC BASE
**                ELASTIC REBAR
**                CPS4R
**
**                CHANGE ORIENTATION AND DIRECTION OF REBAR:
**
** 1)  0  DEGREE ORIENTATION, PARALLEL TO EDGE 1, 1/3 DISTANCE FROM EDGE1  
** 2)  45 DEGREE ORIENTATION, PARALLEL TO EDGE 1, 1/3 DISTANCE FROM EDGE1  
** 3)  90 DEGREE ORIENTATION, PARALLEL TO EDGE 1, 1/3 DISTANCE FROM EDGE1  
** 4)  0  DEGREE ORIENTATION, ALONG EDGE 1  
** 5)  0  DEGREE ORIENTATION, ALONG EDGE 2
** 6)  SAME AS TEST 1, BUT SKEW REBAR
** 7)  SAME AS TEST 2, BUT SKEW REBAR
** 8)  SAME AS TEST 3, BUT SKEW REBAR
** 9)  SAME AS TEST 4, BUT SKEW REBAR
** 10) SAME AS TEST 5, BUT SKEW REBAR
**  
** ANALYTIC RESULTS:
**       TEST CASE  /  REBAR STRESS   REBAR LOG STRAIN
**         1, 6          975.804         0.0487902
**         2, 7          499.800         0.0249900
**         3, 8            0.0           0.0
**         4, 9          975.804         0.0487902
**       **5,10          -38.279        -1.91393E-03
**
**    **Results are consistent with underlying element
**      deformation, but are not compared with an exact
**      analytic solution.
**  
*HEADING
 TEST OF REBAR ORIENTATIONS IN CPS4R ELEMENTS
*PREPRINT,MODEL=YES
**
**
** TEST CASE 1)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
** 
**
*NODE,NSET=NALL
11,  0., 0.
12,  1., 0.
13,  1., 1.
14,  0., 1. 
*NSET,NSET=LEFT
11,14
*NSET,NSET=BOT
11,12
*NSET,NSET=RIGHT
12,13
*NSET,NSET=TOP
13,14
*ELEMENT,TYPE=CPS4R,ELSET=BASE1
11, 11,12,13,14
*SOLID SECTION,ELSET=BASE1,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB1
11, .5, .5, 0., .333333, 1 
*ELSET,ELSET=EALL
BASE1,
**
** TEST CASE 2)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 45.
** 
**
*NODE,NSET=NALL
21,  2., 0.
22,  3., 0.
23,  3., 1.
24,  2., 1. 
*NSET,NSET=LEFT
21,24
*NSET,NSET=BOT
21,22
*NSET,NSET=RIGHT
22,23
*NSET,NSET=TOP
23,24
*ELEMENT,TYPE=CPS4R,ELSET=BASE2
21, 21,22,23,24
*SOLID SECTION,ELSET=BASE2,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB2
21, .5, .5, 45., .333333, 1 
*ELSET,ELSET=EALL
BASE2,
**
** TEST CASE 3)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 89.9
** 
**
*NODE,NSET=NALL
31,  4., 0.
32,  5., 0.
33,  5., 1.
34,  4., 1. 
*NSET,NSET=LEFT
31,34
*NSET,NSET=BOT
31,32
*NSET,NSET=RIGHT
32,33
*NSET,NSET=TOP
33,34
*ELEMENT,TYPE=CPS4R,ELSET=BASE3
31, 31,32,33,34
*SOLID SECTION,ELSET=BASE3,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB3
31, .5, .5, 89.9, .333333, 1 
*ELSET,ELSET=EALL
BASE3,
**
**
** TEST CASE 4)
** MESH: 1 ELEMENT 
** REBAR POSITION: ALONG EDGE 1
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0. 
** 
**
*NODE,NSET=NALL
41,  6., 0.
42,  7., 0.
43,  7., 1.
44,  6., 1. 
*NSET,NSET=LEFT
41,44
*NSET,NSET=BOT
41,42
*NSET,NSET=RIGHT
42,43
*NSET,NSET=TOP
43,44
*ELEMENT,TYPE=CPS4R,ELSET=BASE4
41, 41,42,43,44
*SOLID SECTION,ELSET=BASE4,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB4
41, .5, .5, 0.0, 0.0, 1 
*ELSET,ELSET=EALL
BASE4,
**
** TEST CASE 5)
** MESH: 1 ELEMENT 
** REBAR POSITION: ALONG EDGE 2
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0. 
** 
**
*NODE,NSET=NALL
51,  8., 0.
52,  9., 0.
53,  9., 1.
54,  8., 1. 
*NSET,NSET=LEFT
51,54
*NSET,NSET=BOT
51,52
*NSET,NSET=RIGHT
52,53
*NSET,NSET=TOP
53,54
*ELEMENT,TYPE=CPS4R,ELSET=BASE5
51,51,52,53,54
*SOLID SECTION,ELSET=BASE5,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB5
51, .5, .5, 0.0, 0.0, 2
*ELSET,ELSET=EALL
BASE5,
**
**
** TEST CASE 6)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0.
** 
**
*NODE,NSET=NALL
61,  0., 2.
62,  1., 2.
63,  1., 3.
64,  0., 3. 
*NSET,NSET=LEFT
61,64
*NSET,NSET=BOT
61,62
*NSET,NSET=RIGHT
62,63
*NSET,NSET=TOP
63,64
*ELEMENT,TYPE=CPS4R,ELSET=BASE6
61, 61,62,63,64
*SOLID SECTION,ELSET=BASE6,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB6,GEOMETRY=SKEW
61, .5, .5, 0.
0., .33333, 0., .666666
*ELSET,ELSET=EALL
BASE6,
**
** TEST CASE 7)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 45.
** 
**
*NODE,NSET=NALL
71,  2., 2.
72,  3., 2.
73,  3., 3.
74,  2., 3. 
*NSET,NSET=LEFT
71,74
*NSET,NSET=BOT
71,72
*NSET,NSET=RIGHT
72,73
*NSET,NSET=TOP
73,74
*ELEMENT,TYPE=CPS4R,ELSET=BASE7
71, 71,72,73,74
*SOLID SECTION,ELSET=BASE7,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB7,GEOMETRY=SKEW
71, .5, .5, 45.0
0., .33333, 0., .666666
*ELSET,ELSET=EALL
BASE7,
**
** TEST CASE 8)
** MESH: 1 ELEMENT 
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 89.9
** 
**
*NODE,NSET=NALL
81,  4., 2.
82,  5., 2.
83,  5., 3.
84,  4., 3. 
*NSET,NSET=LEFT
81,84
*NSET,NSET=BOT
81,82
*NSET,NSET=RIGHT
82,83
*NSET,NSET=TOP
83,84
*ELEMENT,TYPE=CPS4R,ELSET=BASE8
81, 81,82,83,84
*SOLID SECTION,ELSET=BASE8,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB8,GEOMETRY=SKEW
81, .5, .5, 89.9
0., .33333, 0., .666666
*ELSET,ELSET=EALL
BASE8,
**
**
** TEST CASE 9)
** MESH: 1 ELEMENT 
** REBAR POSITION: ALONG EDGE 1
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0. 
** 
**
*NODE,NSET=NALL
91,  6., 2.
92,  7., 2.
93,  7., 3.
94,  6., 3. 
*NSET,NSET=LEFT
91,94
*NSET,NSET=BOT
91,92
*NSET,NSET=RIGHT
92,93
*NSET,NSET=TOP
93,94
*ELEMENT,TYPE=CPS4R,ELSET=BASE9
91, 91,92,93,94
*SOLID SECTION,ELSET=BASE9,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB9,GEOMETRY=SKEW
91, .5, .5, 0.
0., 0.0001, 0., 1.0
*ELSET,ELSET=EALL
BASE9,
**
** TEST CASE 10)
** MESH: 1 ELEMENT 
** REBAR POSITION: ALONG EDGE 2
** ELEMENT TYPE: CPS4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** ORIENTATION ANGLE: 0. 
** 
**
*NODE,NSET=NALL
101,  8., 2.
102,  9., 2.
103,  9., 3.
104,  8., 3. 
*NSET,NSET=LEFT
101,104
*NSET,NSET=BOT
101,102
*NSET,NSET=RIGHT
102,103
*NSET,NSET=TOP
103,104
*ELEMENT,TYPE=CPS4R,ELSET=BASE10
101,101,102,103,104
*SOLID SECTION,ELSET=BASE10,MATERIAL=BASE
1.0, 
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB10,GEOMETRY=SKEW
101, .5, .5, 0.
1.0, 0., 0.00001, 0.
*ELSET,ELSET=EALL
BASE10,
**
**
**
*MATERIAL,NAME=BASE
*ELASTIC
1.E4, 0.3
*DENSITY
1.E-4, 
*MATERIAL,NAME=REBAR
*ELASTIC
2.E4, 
*DENSITY
1.E-4, 
**
**
**
****************
** LOAD CASES **
****************
*AMPLITUDE,NAME=SMOOTH, DEFINITION=SMOOTH STEP
0.0,0.0,.01,1.0
*STEP
*DYNAMIC,EXPLICIT
,.01
**
*BOUNDARY,AMPLITUDE=SMOOTH
LEFT,1
BOT,2
RIGHT,1,1,.05
TOP,2
**
*RESTART,TIMEMARKS=YES,WRITE,NUM=1
*FILE OUTPUT,NUMBER INTERVAL=2,TIMEMARKS=YES
*EL FILE
LE,S,MISES,PRESS,PEEQ,TEMP,ENER
ELEN, 
*NODE FILE
RF, 
*ENERGY FILE
***********************************************
*OUTPUT,HISTORY,TIME INTERVAL=0.01
*ENERGY OUTPUT, VAR=ALL
***********************************************
***OUTPUT,HISTORY,TIME INTERVAL=.00001
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLWK,ETOTAL,
**DT,
***ELEMENT OUTPUT,ELSET=EALL
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE1,REBAR=REB1
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE2,REBAR=REB2
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE3,REBAR=REB3
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE4,REBAR=REB4
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE5,REBAR=REB5
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE6,REBAR=REB6
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE7,REBAR=REB7
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE8,REBAR=REB8
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE9,REBAR=REB9
**S,LE,ENER
**ELEN,
***ELEMENT OUTPUT,ELSET=BASE10,REBAR=REB10
**S,LE,ENER
**ELEN,
***NODE OUTPUT,NSET=RIGHT
**RF,
*END STEP