**
** TEST OF REBAR: THERMAL STRAIN FOR VARIOUS ELEMENTS WITH REBAR
**                ELASTIC BASE
**                ELASTIC REBAR
**
** 1)  CPS4R 
** 2)  CPE4R
** 3)  CAX4R
** 4)  C3D8R
** 5)  M3D4R
** 6)  S4R,   TEMPERATURE GRADIENT
** 7)  S4RS,  TEMPERATURE GRADIENT
** 8)  S4RSW, TEMPERATURE GRADIENT
** 9)  SAX1,  TEMPERATURE GRADIENT
** 10) S4R,   TEMPERATURE AT SECTION POINTS
** 11) S4RS,  TEMPERATURE AT SECTION POINTS
** 12) S4RSW, TEMPERATURE AT SECTION POINTS
** 13) SAX1,  TEMPERATURE AT SECTION POINTS
** 14) S4R,   TEMPERATURE GRADIENT, COMPOSITE
** 15) S4R,   TEMPERATURE AT SECTION POINTS, COMPOSITE
** 16) S4R,   TEMPERATURE GRADIENT, COMPOSITE, NODAL THICKNESS
** 17) S4R,   TEMPERATURE AT SECTION POINTS, COMPOSITE, NODAL THICKNESS
**
**  ANALYTIC RESULTS: THE STRESS IN THE REBAR AT THE END OF THE
**                    STEP SHOULD BE 160.0
**  
** TEMPERATURES IN THE LAYERED SHELLS BASE15 and BASE17 and REBAR should be:
** from SNEG-->SPOS: 20, 27.5, 35, 42.5, 50, 57.5, 65, 72.5, 80
** TEMPS in their REBARS should be 40
**
*HEADING
 DATA DECK TO VERIFY THERMAL STRAIN FOR DIFFERENT ELEMENTS WITH REBARS
**
** TEST CASE 1)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPS4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL1
11,0.,0.
12,1.,0.
13,1.,1.
14,0.,1.
*NSET,NSET=NALL
ALL1, 
*NSET,NSET=TOP
13,14
*NSET,NSET=BOT
11,12
*NSET,NSET=LEFT
11,14
*NSET,NSET=RIGHT
12,13
*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., .33333333, 1
*ELSET,ELSET=EALL
BASE1,  
**
** TEST CASE 2)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPE4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL2
21,1.5,0.
22,2.5,0.
23,2.5,1.
24,1.5,1.
*NSET,NSET=NALL
ALL2, 
*NSET,NSET=TOP
23,24
*NSET,NSET=BOT
21,22
*NSET,NSET=LEFT
21,24
*NSET,NSET=RIGHT
22,23
*ELEMENT,TYPE=CPE4R,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, 0., .33333333, 1 
*ELSET,ELSET=EALL
BASE2, 
**
** TEST CASE 3)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CAX4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL3
31,3.,0.
32,4.,0.
33,4.,1.
34,3.,1.
*NSET,NSET=NALL
ALL3, 
*NSET,NSET=TOP
33,34
*NSET,NSET=BOT
31,32
*NSET,NSET=LEFT
31,34
*NSET,NSET=RIGHT
32,33
*ELEMENT,TYPE=CAX4R,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, 0., .33333333, 1, ,1.5
*ELSET,ELSET=EALL
BASE3, 
**
** TEST CASE 4)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: C3D8R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL4
41,  4.5, 0.,  0.
42,  5.5, 0.,  0.
43,  5.5, 0., -1.
44,  4.5, 0., -1.
45,  4.5, 1.,  0.
46,  5.5, 1.,  0.
47,  5.5, 1., -1.
48,  4.5, 1., -1.
*NSET,NSET=NALL
ALL4, 
*NSET,NSET=TOP
45,46,47,48
*NSET,NSET=BOT
41,42,43,44
*NSET,NSET=LEFT
41,44,45,48
*NSET,NSET=RIGHT
42,43,46,47
*ELEMENT,TYPE=C3D8R,ELSET=BASE4
41,41,42,43,44,45,46,47,48
*SOLID SECTION,ELSET=BASE4,MATERIAL=BASE
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB4
41, .5, .5, 90., .33333333, 1, 1
*ELSET,ELSET=EALL
BASE4, 
**
** TEST CASE 5)
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: M3D4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL5
51,  6.0, 0.,  0.
52,  7.0, 0.,  0.
53,  7.0, 0., -1.
54,  6.0, 0., -1.
*NSET,NSET=TEMPMEMB
ALL5, 
*NSET,NSET=NALL
ALL5, 
*NSET,NSET=LEFT
51,54
*NSET,NSET=RIGHT
52,53
*ELEMENT,TYPE=M3D4R,ELSET=BASE5
51,51,52,53,54
*MEMBRANE SECTION,ELSET=BASE5,MATERIAL=BASE
1.0,
*REBAR LAYER
REB5,0.5,0.5,,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE5, 
**
** TEST CASE 6)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R
** TEMPERATURE: APPLIED AS GRADIENT
**
*NODE,NSET=ALL6
61,  0.0, 3.,  0.
62,  1.0, 3.,  0.
63,  1.0, 3., -1.
64,  0.0, 3., -1.
*NSET,NSET=NALL
ALL6, 
*NSET,NSET=TEMPGRAD
ALL6, 
*NSET,NSET=LEFT
61,64
*NSET,NSET=RIGHT
62,63
*ELEMENT,TYPE=S4R,ELSET=BASE6
61,61,62,63,64
*SHELL SECTION,ELSET=BASE6,MATERIAL=BASE
1.0,
*REBAR LAYER
REB6,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE6, 
**
** TEST CASE 7)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4RS
** TEMPERATURE: APPLIED AS GRADIENT
**
*NODE,NSET=ALL7
71,  2.0, 3.,  0.
72,  3.0, 3.,  0.
73,  3.0, 3., -1.
74,  2.0, 3., -1.
*NSET,NSET=NALL
ALL7, 
*NSET,NSET=TEMPGRAD
ALL7, 
*NSET,NSET=LEFT
71,74
*NSET,NSET=RIGHT
72,73
*ELEMENT,TYPE=S4RS,ELSET=BASE7
71,71,72,73,74
*SHELL SECTION,ELSET=BASE7,MATERIAL=BASE
1.0,
*REBAR LAYER
REB7,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE7, 
**
** TEST CASE 8)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4RSW
** TEMPERATURE: APPLIED AS GRADIENT
**
*NODE,NSET=ALL8
81,  4.0, 3.,  0.
82,  5.0, 3.,  0.
83,  5.0, 3., -1.
84,  4.0, 3., -1.
*NSET,NSET=NALL
ALL8, 
*NSET,NSET=TEMPGRAD
ALL8, 
*NSET,NSET=LEFT
81,84
*NSET,NSET=RIGHT
82,83
*ELEMENT,TYPE=S4RSW,ELSET=BASE8
81,81,82,83,84
*SHELL SECTION,ELSET=BASE8,MATERIAL=BASE
1.0,
*REBAR LAYER
REB8,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE8, 
**
** TEST CASE 9)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: SAX1
** TEMPERATURE: APPLIED AS GRADIENT
**
*NODE,NSET=ALL9
91,  6.0, 3.,  0.
92,  7.0, 3.,  0.
*NSET,NSET=NALL
ALL9, 
*NSET,NSET=TEMPGRAD
ALL9, 
*NSET,NSET=LEFT
91, 
*NSET,NSET=RIGHT
92, 
*ELEMENT,TYPE=SAX1,ELSET=BASE9
91,91,92
*SHELL SECTION,ELSET=BASE9,MATERIAL=BASE
1.0,
*REBAR LAYER, GEOMETRY=ANGULAR
REB9,0.5,6.37,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE9, 
**
**
** TEST CASE 10)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R
** TEMPERATURE: APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL10
101,  0.0, 5.,  0.
102,  1.0, 5.,  0.
103,  1.0, 5., -1.
104,  0.0, 5., -1.
*NSET,NSET=NALL
ALL10, 
*NSET,NSET=TEMPPOINTS
ALL10, 
*NSET,NSET=LEFT
101,104
*NSET,NSET=RIGHT
102,103
*ELEMENT,TYPE=S4R,ELSET=BASE10
101,101,102,103,104
*SHELL SECTION,ELSET=BASE10,MATERIAL=BASE,TEMPERATURE=3
1.0,
*REBAR LAYER
REB10,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE10,
**
** TEST CASE 11)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4RS
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL11
111,  2.0, 5.,  0.
112,  3.0, 5.,  0.
113,  3.0, 5., -1.
114,  2.0, 5., -1.
*NSET,NSET=NALL
ALL11, 
*NSET,NSET=TEMPPOINTS
ALL11, 
*NSET,NSET=LEFT
111,114
*NSET,NSET=RIGHT
112,113
*ELEMENT,TYPE=S4RS,ELSET=BASE11
111,111,112,113,114
*SHELL SECTION,ELSET=BASE11,MATERIAL=BASE,TEMPERATURE=3
1.0,
*REBAR LAYER
REB11,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE11,
**
** TEST CASE 12)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4RSW
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL12
121,  4.0, 5.,  0.
122,  5.0, 5.,  0.
123,  5.0, 5., -1.
124,  4.0, 5., -1.
*NSET,NSET=NALL
ALL12, 
*NSET,NSET=TEMPPOINTS
ALL12, 
*NSET,NSET=LEFT
121,124
*NSET,NSET=RIGHT
122,123
*ELEMENT,TYPE=S4RSW,ELSET=BASE12
121,121,122,123,124
*SHELL SECTION,ELSET=BASE12,MATERIAL=BASE,TEMPERATURE=3
1.0,
*REBAR LAYER
REB12,0.5,0.5,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE12,
**
** TEST CASE 13)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: SAX1
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL13
131,  6.0, 5.,  0.
132,  7.0, 5.,  0.
*NSET,NSET=NALL
ALL13, 
*NSET,NSET=TEMPPOINTS
ALL13, 
*ELEMENT,TYPE=SAX1,ELSET=BASE13
131,131,132
*SHELL SECTION,ELSET=BASE13,MATERIAL=BASE,TEMPERATURE=3
1.0,
*REBAR LAYER, GEOMETRY=ANGULAR
REB13,0.5,6.37,-.1666667,REBAR,0.0,1  
*ELSET,ELSET=EALL
BASE13,
**
** TEST CASE 14)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R, COMPOSITE
** TEMPERATURE:  APPLIED AS GRADIENT
**
*NODE,NSET=ALL14
141,  0.0, 7.,  0.
142,  1.0, 7.,  0.
143,  1.0, 7., -1.
144,  0.0, 7., -1.
*NSET,NSET=NALL
ALL14, 
*NSET,NSET=TEMPGRAD
ALL14, 
*NSET,NSET=LEFT
141,144
*NSET,NSET=RIGHT
142,143
*ELEMENT,TYPE=S4R,ELSET=BASE14
141,141,142,143,144
*SHELL SECTION,ELSET=BASE14,COMPOSITE
.25, 3, BASE
.50, 5, BASE
.25, 3, BASE
*REBAR LAYER
REB14,0.5,0.5,-.1666667,REBAR,0.0,1 
*ELSET,ELSET=EALL
BASE14,
**
** TEST CASE 15)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R, COMPOSITE
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL15
151,  2.0, 7.,  0.
152,  3.0, 7.,  0.
153,  3.0, 7., -1.
154,  2.0, 7., -1.
*NSET,NSET=NALL
ALL15, 
*NSET,NSET=TEMPPOINTS-COMP
ALL15, 
*NSET,NSET=LEFT
151,154
*NSET,NSET=RIGHT
152,153
*ELEMENT,TYPE=S4R,ELSET=BASE15
151,151,152,153,154
*SHELL SECTION,ELSET=BASE15,COMPOSITE,TEMPERATURE=3
.25, 3, BASE
.50, 5, BASE
.25, 3, BASE
*REBAR LAYER
REB15,0.5,0.5,-.1666667,REBAR,0.0,1 
*ELSET,ELSET=EALL
BASE15,
**
** TEST CASE 16)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R, COMPOSITE, NODAL THICKNESS
** TEMPERATURE:  APPLIED AS GRADIENT
**
*NODE,NSET=ALL16
161,  4.0, 7.,  0.
162,  5.0, 7.,  0.
163,  5.0, 7., -1.
164,  4.0, 7., -1.
*NSET,NSET=NALL
ALL16, 
*NSET,NSET=TEMPGRAD
ALL16, 
*NSET,NSET=LEFT
161,164
*NSET,NSET=RIGHT
162,163
*ELEMENT,TYPE=S4R,ELSET=BASE16
161,161,162,163,164
*SHELL SECTION,ELSET=BASE16,COMPOSITE,NODAL THICKNESS
.125, 3, BASE
.250, 5, BASE
.125, 3, BASE
*REBAR LAYER
REB16,0.5,0.5,-.08333333,REBAR,0.0,1 
*ELSET,ELSET=EALL
BASE16,
*NODAL THICKNESS
161, .5
164, .5
162, 1.5
163, 1.5
**
** TEST CASE 17)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R, COMPOSITE, NODAL THICKNESS
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL17
171,  6.0, 7.,  0.
172,  7.0, 7.,  0.
173,  7.0, 7., -1.
174,  6.0, 7., -1.
*NSET,NSET=NALL
ALL17, 
*NSET,NSET=TEMPPOINTS-COMP
ALL17, 
*NSET,NSET=LEFT
171,174
*NSET,NSET=RIGHT
172,173
*ELEMENT,TYPE=S4R,ELSET=BASE17
171,171,172,173,174
*SHELL SECTION,ELSET=BASE17,COMPOSITE,NODAL THICKNESS,TEMPERATURE=3
.125, 3, BASE
.250, 5, BASE
.125, 3, BASE
*REBAR LAYER
REB17,0.5,0.5,-.08333333, REBAR, 0.0,1 
*ELSET,ELSET=EALL
BASE17,
*NODAL THICKNESS
171, .5
174, .5
172, 1.5
173, 1.5
**
**
** 
**
*MATERIAL,NAME=BASE
*ELASTIC
1.E4, 0.495
*DENSITY
1.E-3, 
*EXPANSION
1.E-4, 
*MATERIAL,NAME=REBAR
*ELASTIC
2.E4, 
*DENSITY
1.E-3, 
*EXPANSION
2.E-4, 
********************************
**    THERMAL LOADING STEP    **
******************************** 
*INITIAL CONDITION,TYPE=TEMPERATURE
TOP,0.
BOT,0.
TEMPMEMB, 0.0
TEMPPOINTS,0.,0.,0.
TEMPPOINTS-COMP,0.,0,0.0,0.,0.0,0.0,0.0
TEMPGRAD,0.,0.
**TOP,0.
**BOT,0.
**TEMPPOINTS, 0.,0.,0.
**TEMPPOINTS-COMP, 0.,0.,0.,0.,0.,0.,0.
**TEMPGRAD,   0.,0.
*AMPLITUDE,NAME=RAMP
0.0,0.0,.05,1.0
*RESTART,TIMEMARKS=YES,WRITE,NUMBER=1
**
**
*ELSET,ELSET=TEMP_COMP
BASE15, BASE17
*STEP
*DYNAMIC,EXPLICIT
,.05
**
*TEMPERATURE,AMP=RAMP
TOP,80.
BOT,20.
TEMPMEMB, 40.
TEMPPOINTS,20.,50.,80.
TEMPPOINTS-COMP,20.,27.5,35.0,50.,65.0,72.5,80.0
TEMPGRAD,50.,60.
**
*BOUNDARY
NALL, 1,6
**
**
*FILE OUTPUT,NUMBER INTERVAL=1,TIMEMARKS=YES
*EL FILE
LE,S,MISES,PRESS,TEMP,ENER
ELEN, 
*NODE FILE
RF, 
*ENERGY FILE
***************************
**REBAR OUTPUT AVAILABLE **
** for M3D4R ELEMENTS    **
***************************
***
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT, REBAR=REB5
LE,S,TEMP,ENER
************
***OUTPUT,HISTORY,TIME INTERVAL=0.001
***ELEMENT OUTPUT,ELSET=BASE5,REBAR=REB5
**S11,LE11,RBFOR
***************************
*OUTPUT,FIELD,NUMBER INTERVAL = 1 
*ELEMENT OUTPUT
LE,S,MISES,PRESS,TEMP,ENER
ELEN
********************************
*******TEMP IN COMPOSITES*******
*OUTPUT,FIELD,NUMBER INTERVAL=2
*ELEMENT OUTPUT,ELSET=TEMP_COMP
1,2,3,4,5,6,7,8,9,9,10,11,12
TEMP
*ELEMENT OUTPUT,ELSET=TEMP_COMP,REBAR
TEMP
*END STEP