**
** 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
*NSET,NSET=NALL
ALL1,
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=CPS4R,ELSET=BASE1
*SOLID SECTION,ELSET=BASE1, MATERIAL=BASE
1.0,
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB1
*ELSET,ELSET=EALL
BASE1,
**
** TEST CASE 2)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CPE4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL2
*NSET,NSET=NALL
ALL2,
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=CPE4R,ELSET=BASE2
*SOLID SECTION,ELSET=BASE2, MATERIAL=BASE
1.0,
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB2
*ELSET,ELSET=EALL
BASE2,
**
** TEST CASE 3)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: CAX4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL3
*NSET,NSET=NALL
ALL3,
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=CAX4R,ELSET=BASE3
*SOLID SECTION,ELSET=BASE3, MATERIAL=BASE
1.0,
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB3
*ELSET,ELSET=EALL
BASE3,
**
** TEST CASE 4)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: C3D8R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL4
*NSET,NSET=NALL
ALL4,
*NSET,NSET=TOP
*NSET,NSET=BOT
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=C3D8R,ELSET=BASE4
*SOLID SECTION,ELSET=BASE4,MATERIAL=BASE
*REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR,NAME=REB4
*ELSET,ELSET=EALL
BASE4,
**
** TEST CASE 5)
** REBAR POSITION: MIDPLANE
** ELEMENT TYPE: M3D4R
** TEMPERATURE: APPLIED AT NODE
**
*NODE,NSET=ALL5
*NSET,NSET=TEMPMEMB
ALL5,
*NSET,NSET=NALL
ALL5,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=M3D4R,ELSET=BASE5
*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
*NSET,NSET=NALL
ALL6,
*NSET,NSET=TEMPGRAD
ALL6,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE6
*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
*NSET,NSET=NALL
ALL7,
*NSET,NSET=TEMPGRAD
ALL7,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4RS,ELSET=BASE7
*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
*NSET,NSET=NALL
ALL8,
*NSET,NSET=TEMPGRAD
ALL8,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4RSW,ELSET=BASE8
*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
*NSET,NSET=NALL
ALL9,
*NSET,NSET=TEMPGRAD
ALL9,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=SAX1,ELSET=BASE9
*SHELL SECTION,ELSET=BASE9,MATERIAL=BASE
1.0,
*REBAR LAYER, ANGULAR SPACING
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
*NSET,NSET=NALL
ALL10,
*NSET,NSET=TEMPPOINTS
ALL10,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE10
*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
*NSET,NSET=NALL
ALL11,
*NSET,NSET=TEMPPOINTS
ALL11,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4RS,ELSET=BASE11
*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
*NSET,NSET=NALL
ALL12,
*NSET,NSET=TEMPPOINTS
ALL12,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4RSW,ELSET=BASE12
*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
*NSET,NSET=NALL
ALL13,
*NSET,NSET=TEMPPOINTS
ALL13,
*ELEMENT,TYPE=SAX1,ELSET=BASE13
*SHELL SECTION,ELSET=BASE13,MATERIAL=BASE,TEMPERATURE=3
1.0,
*REBAR LAYER, ANGULAR SPACING
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
*NSET,NSET=NALL
ALL14,
*NSET,NSET=TEMPGRAD
ALL14,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE14
*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
*NSET,NSET=NALL
ALL15,
*NSET,NSET=TEMPPOINTS-COMP
ALL15,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE15
*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
*NSET,NSET=NALL
ALL16,
*NSET,NSET=TEMPGRAD
ALL16,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE16
*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
**
** TEST CASE 17)
** REBAR POSITION: 1/3 FROM BOTTOM
** ELEMENT TYPE: S4R, COMPOSITE, NODAL THICKNESS
** TEMPERATURE:  APPLIED AT SECTION POINTS
**
*NODE,NSET=ALL17
*NSET,NSET=NALL
ALL17,
*NSET,NSET=TEMPPOINTS-COMP
ALL17,
*NSET,NSET=LEFT
*NSET,NSET=RIGHT
*ELEMENT,TYPE=S4R,ELSET=BASE17
*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
**
**
**
*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
TEMP
*ELEMENT OUTPUT,ELSET=TEMP_COMP,REBAR
TEMP
*END STEP