** TEST OF REBAR: THERMAL LOAD, TEMP AND FV MATERIAL DEPENDENCE
**
** ELEMENTS TESTED: M3D4, C3D8R, S4R
** 
**                ELASTIC, FV2 DEPENDENT BASE
**                ELASTIC, TEMP AND FV1 DEPENDENT REBAR
**
** STEP 1: FV1=0, FV2=0, DELTA-T=100
** STEP 2: FV1=1, FV2=0, DELTA-T=100
** STEP 3: FV1=1, FV2=1, DELTA-T=100
** 
** ORIENTATION OF REBAR: ALONG X-DIRECTION
**  
** ANALYTIC RESULTS:
**
**           SOLID     REBAR       REBAR           REBAR
**          STRESS    STRESS   ELASTIC STRAIN  THERMAL STRAIN
**
** STEP 1:   900.0    450.00      -.015            .015
** STEP 2:   900.0    600.00      -.015            .015    
** STEP 3:  1200.0    600.00      -.015            .015    
**
*HEADING
 THERMAL LOADING OF REBAR, MATERIAL TESTS
*PREPRINT,MODEL=YES
**
**
** TEST CASE 1)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: C3D8R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
11,  0., 0.,  0.
12,  1., 0.,  0.
13,  1., 1.,  0.
14,  0., 1.,  0.
15,  0., 0.,  1.
16,  1., 0.,  1.
17,  1., 1.,  1.
18,  0., 1.,  1.
1224,           0.,           1.,          0.5
1225,           0.,           0.,          0.5
1226,           1.,           0.,          0.5
1227,           1.,           1.,          0.5
*NSET,NSET=BOT
11,12,15,16
*NSET,NSET=LEFT
11,14,15,18
*NSET,NSET=ALLZ
11,12,13,14
*NSET,NSET=RIGHT
12,13,16,17
*ELEMENT,TYPE=C3D8R,ELSET=SOLID1
11, 11,12,13,14,15,16,17,18
*SOLID SECTION,ELSET=SOLID1,MATERIAL=SOLID1
*ORIENTATION,NAME=REBAR1,DEFINITION=OFFSET TO NODES
2,4
3
*ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-1
1221,1224,1225,1226,1227
*SURFACE SECTION,ELSET=SFM3D4R-1
*REBAR LAYER,ORIENTATION=REBAR1
REB1,0.500000,1.000000,0.000000,REB1,90.000000,1
*EMBEDDED ELEMENT,HOST ELSET=SOLID1
1221,
**REBAR,ELEMENT=CONTINUUM,MATERIAL=REB1,NAME=REB1
**11, .5, 1.0, 90., .5, 1, 1 
*MATERIAL,NAME=SOLID1
*ELASTIC,DEPENDENCIES=2
6.0E4, .0, 0., 0., 0.
8.0E4, .0, 0., 0., 1.
*EXPANSION
1.5E-4,
*MATERIAL,NAME=REB1
*ELASTIC,DEPENDENCIES=2
2.0E4, ,   0., 0., 0.
3.0E4, , 100., 0., 0.
3.5E4, ,   0., 1., 0.
4.0E4, , 100., 1., 0.
2.0E4, ,   0., 0., 1.
3.0E4, , 100., 0., 1.
3.5E4, ,   0., 1., 1.
4.0E4, , 100., 1., 1.
*EXPANSION
1.5E-4,
**
**
** TEST CASE 2)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=ALLZ
111,  2., 0.
112,  3., 0.
113,  3., 1.
114,  2., 1.
*NSET,NSET=BOT
111,112
*NSET,NSET=LEFT
111,114
*NSET,NSET=RIGHT
112,113
*ELEMENT,TYPE=M3D4,ELSET=MEMB1
111, 111,112,113,114
*MEMBRANE SECTION,ELSET=MEMB1,MATERIAL=MEMB1
1.0,
*REBAR LAYER
REB2,0.5,1.0,0,REB2,0,1
*MATERIAL,NAME=MEMB1
*ELASTIC,DEPENDENCIES=2
6.0E4, .0, 0., 0., 0.
8.0E4, .0, 0., 0., 1.
*EXPANSION
1.5E-4,
*MATERIAL,NAME=REB2
*ELASTIC,DEPENDENCIES=2
2.0E4, ,   0., 0., 0.
3.0E4, , 100., 0., 0.
3.5E4, ,   0., 1., 0.
4.0E4, , 100., 1., 0.
2.0E4, ,   0., 0., 1.
3.0E4, , 100., 0., 1.
3.5E4, ,   0., 1., 1.
4.0E4, , 100., 1., 1.
*EXPANSION
1.5E-4,
**
**
** TEST CASE 3)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=ALLZ
211,  4., 0.
212,  5., 0.
213,  5., 1.
214,  4., 1.
*NSET,NSET=BOT
211,212
*NSET,NSET=LEFT
211,214
*NSET,NSET=RIGHT
212,213
*ELEMENT,TYPE=S4R,ELSET=SHELL1
211, 211,212,213,214
*SHELL SECTION,ELSET=SHELL1,MATERIAL=SHELL1
1.0,
*REBAR LAYER
REB3,0.5,1.0,0,REB3,0,1
*MATERIAL,NAME=SHELL1
*ELASTIC,DEPENDENCIES=2
6.0E4, .0, 0., 0., 0.
8.0E4, .0, 0., 0., 1.
*EXPANSION
1.5E-4,
*MATERIAL,NAME=REB3
*ELASTIC,DEPENDENCIES=2
2.0E4, ,   0., 0., 0.
3.0E4, , 100., 0., 0.
3.5E4, ,   0., 1., 0.
4.0E4, , 100., 1., 0.
2.0E4, ,   0., 0., 1.
3.0E4, , 100., 0., 1.
3.5E4, ,   0., 1., 1.
4.0E4, , 100., 1., 1.
*EXPANSION
1.5E-4,
**
**
**
** TEST CASE 4)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=ALLZ
221,  6., 0.
222,  7., 0.
223,  7., 1.
224,  6., 1.
*NSET,NSET=BOT
221,222
*NSET,NSET=LEFT
221,224
*NSET,NSET=RIGHT
222,223
*ELEMENT,TYPE=S4,ELSET=SHELL2
221, 221,222,223,224
*SHELL SECTION,ELSET=SHELL2,MATERIAL=SHELL2
1.0,
*REBAR LAYER
REB4,0.5,1.0,0,REB4,0,1
*MATERIAL,NAME=SHELL2
*ELASTIC,DEPENDENCIES=2
6.0E4, .0, 0., 0., 0.
8.0E4, .0, 0., 0., 1.
*EXPANSION
1.5E-4,
*MATERIAL,NAME=REB4
*ELASTIC,DEPENDENCIES=2
2.0E4, ,   0., 0., 0.
3.0E4, , 100., 0., 0.
3.5E4, ,   0., 1., 0.
4.0E4, , 100., 1., 0.
2.0E4, ,   0., 0., 1.
3.0E4, , 100., 0., 1.
3.5E4, ,   0., 1., 1.
4.0E4, , 100., 1., 1.
*EXPANSION
1.5E-4,
*ELSET,ELSET=EALL
SOLID1,MEMB1,SHELL1,SHELL2
*NSET,NSET=ALLN
NALL,ALLZ
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALLN, 0.
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=1
ALLN, 0.
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=2
ALLN, 0.
****************
** LOAD CASES **
****************
*STEP,NLGEOM
*STATIC
1.,1.
*RESTART,WRITE,FREQ=99
*NODE PRINT,FREQ=0
*BOUNDARY
ALLN,1,2
ALLZ,3
*TEMPERATURE
ALLN, 100.
**
*EL PRINT,ELSET=EALL,FREQ=99
S,
**
** TEST CASE 1)
**
*EL PRINT,ELSET=SOLID1,REBAR=REB1,FREQ=99
S,EE,THE,SENER,IVOL,RBFOR
*EL FILE,ELSET=SOLID1,REBAR=REB1,FREQ=99
S,E,ENER,IVOL,RBFOR
*EL FILE,ELSET=SFM3D4R-1,REBAR=REB1,FREQ=99
S,E,ENER,IVOL,RBFOR
**
** TEST CASE 2)
**
*EL PRINT,ELSET=MEMB1,REBAR=REB2,FREQ=99
S,EE,THE,SENER,IVOL,RBFOR
*EL FILE,ELSET=MEMB1,REBAR=REB2,FREQ=99
S,EE,THE,ENER,IVOL,RBFOR
**
** TEST CASE 3)
**
*EL PRINT,ELSET=SHELL1,REBAR=REB3,FREQ=99
S,EE,THE,SENER,IVOL
*EL FILE,ELSET=SHELL1,REBAR=REB3,FREQ=99
S,EE,THE,ENER,IVOL
**
** TEST CASE 4)
**
*EL PRINT,ELSET=SHELL2,REBAR=REB4,FREQ=99
S,EE,THE,SENER,IVOL
*EL FILE,ELSET=SHELL2,REBAR=REB4,FREQ=99
S,EE,THE,ENER,IVOL
*END STEP
*STEP,NLGEOM
*STATIC
1.,1.
*FIELD,VARIABLE=1
ALLN, 1.
*END STEP
*STEP,NLGEOM
*STATIC
1.,1.
*FIELD,VARIABLE=2
ALLN, 1.
*END STEP