** TEST OF REBAR: INITIAL STRESSES
**
** ELEMENTS TESTED: M3D4, C3D8R, S4R
** 
**             CASE 1-3: INITIAL STRESSES IN REBAR (-100.0)
**                  
**             CASE 4-6: PRESTRESS HOLD
**
** ORIENTATION OF REBAR: ALONG X-DIRECTION
**  
** ANALYTIC RESULTS:
**               SOLID STRESS   REBAR STRESS   REBAR STRAIN
** CASE   1-3:     -2.5.           5.0         -1.25E-03
** CASE   4-6:     -5.0           10.0            0.0
**
*HEADING
 UNIAXIAL TEST OF REBAR, INITIAL STRESS, PRESTRESS HOLD
*PREPRINT,MODEL=YES
**
**
** TEST CASE 1)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: C3D8R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALLS
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.
1234,           0.,           1.,          0.5
1235,           0.,           0.,          0.5
1236,           1.,           0.,          0.5
1237,           1.,           1.,          0.5
*NSET,NSET=BOT
11,12,15,16
*NSET,NSET=LEFT
11,14,15,18
*NSET,NSET=BACK
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
1231,1234,1235,1236,1237
*SURFACE SECTION,ELSET=SFM3D4R-1
*REBAR LAYER,ORIENTATION=REBAR1
REBAR_SET-1,0.500000,1.000000,0.000000,REBAR_MATERIAL-1,90.000000,1
*EMBEDDED ELEMENT
1231,
**REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR_MATERIAL-1,NAME=REBAR_SET-1
**11, .5, 1.0, 90., .5, 1, 1 
*MATERIAL,NAME=SOLID1
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-1
*ELASTIC
2.E4, 
**
**
** TEST CASE 2)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
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
REBAR_SET-2,0.5,1.0,0,REBAR_MATERIAL-2,0,1 
*MATERIAL,NAME=MEMB1
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-2
*ELASTIC
2.E4, 
**
**
** TEST CASE 3)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
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
REBAR_SET-3,0.5,1.0,0,REBAR_MATERIAL-3,0,1 
*MATERIAL,NAME=SHELL1
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-3
*ELASTIC
2.E4, 
**
** TEST CASE 4)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: C3D8R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALLS
21,  0., 2.,  0.
22,  1., 2.,  0.
23,  1., 3.,  0.
24,  0., 3.,  0.
25,  0., 2.,  1.
26,  1., 2.,  1.
27,  1., 3.,  1.
28,  0., 3.,  1.
2234,           0.,           3.,          0.5
2235,           0.,           2.,          0.5
2236,           1.,           2.,          0.5
2237,           1.,           3.,          0.5
*NSET,NSET=BOT
21,22,25,26
*NSET,NSET=LEFT
21,24,25,28
*NSET,NSET=BACK
21,22,23,24
*NSET,NSET=RIGHT
22,23,26,27
*ELEMENT,TYPE=C3D8R,ELSET=SOLID2
21, 21,22,23,24,25,26,27,28
*SOLID SECTION,ELSET=SOLID2,MATERIAL=SOLID2
*ELEMENT,TYPE=SFM3D4R,ELSET=SFM3D4R-2
2231,2234,2235,2236,2237
*SURFACE SECTION,ELSET=SFM3D4R-2
*REBAR LAYER,ORIENTATION=REBAR1
REBAR_SET-4,0.500000,1.000000,0.000000,REBAR_MATERIAL-4,90.000000,1
*EMBEDDED ELEMENT
2231,
**REBAR,ELEMENT=CONTINUUM,MATERIAL=REBAR_MATERIAL-4,NAME=REBAR_SET-4
**21, .5, 1.0, 90., .5, 1, 1 
*MATERIAL,NAME=SOLID2
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-4
*ELASTIC
2.E4, 
**
**
** TEST CASE 5)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
121,  2., 2.
122,  3., 2.
123,  3., 3.
124,  2., 3.
*NSET,NSET=BOT
121,122
*NSET,NSET=LEFT
121,124
*NSET,NSET=RIGHT
122,123
*ELEMENT,TYPE=M3D4,ELSET=MEMB2
121, 121,122,123,124
*MEMBRANE SECTION,ELSET=MEMB2,MATERIAL=MEMB2
1.0,
*REBAR LAYER
REBAR_SET-5,0.5,1.0,0,REBAR_MATERIAL-5,0,1 
*MATERIAL,NAME=MEMB2
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-5
*ELASTIC
2.E4, 
**
**
** TEST CASE 6)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
221,  4., 2.
222,  5., 2.
223,  5., 3.
224,  4., 3.
*NSET,NSET=BOT
221,222
*NSET,NSET=LEFT
221,224
*NSET,NSET=RIGHT
222,223
*ELEMENT,TYPE=S4R,ELSET=SHELL2
221, 221,222,223,224
*SHELL SECTION,ELSET=SHELL2,MATERIAL=SHELL2
1.0,
*REBAR LAYER
REBAR_SET-6,0.5,1.0,0,REBAR_MATERIAL-6,0,1
*MATERIAL,NAME=SHELL2
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-6
*ELASTIC
2.E4, 
**
** TEST CASE 7)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
231,  6., 2.
232,  7., 2.
233,  7., 3.
234,  6., 3.
*NSET,NSET=BOT
231,232
*NSET,NSET=LEFT
231,234
*NSET,NSET=RIGHT
232,233
*ELEMENT,TYPE=S4,ELSET=SHELL3
231, 231,232,233,234
*SHELL SECTION,ELSET=SHELL3,MATERIAL=SHELL3
1.0,
*REBAR LAYER
REBAR_SET-7,0.5,1.0,0,REBAR_MATERIAL-7,0,1
*MATERIAL,NAME=SHELL3
*ELASTIC
1.E4, 0.25
*MATERIAL,NAME=REBAR_MATERIAL-7
*ELASTIC
2.E4, 
*ELSET,ELSET=ALL_ELEMENTS-2
SOLID2,MEMB2,SHELL2,SHELL3,SFM3D4R-2
*ELSET,ELSET=ALL_ELEMENTS
SOLID1,MEMB1,SHELL1
SOLID2,MEMB2,SHELL2,SHELL3,SFM3D4R-1,SFM3D4R-2
*INITIAL CONDITIONS,TYPE=STRESS,REBAR
**ALL_ELEMENTS, REBAR_SET-1, 10.
ALL_ELEMENTS, REBAR_SET-2, 10.
ALL_ELEMENTS, REBAR_SET-3, 10.
ALL_ELEMENTS, REBAR_SET-4, 10.
ALL_ELEMENTS, REBAR_SET-5, 10.
ALL_ELEMENTS, REBAR_SET-6, 10.
ALL_ELEMENTS, REBAR_SET-7, 10.
*INITIAL CONDITIONS,TYPE=STRESS,USER
**
****************
** LOAD CASES **
****************
*STEP,NLGEOM
*STATIC
1.,1.
*RESTART,WRITE,FREQ=99
*NODE PRINT,NSET=RIGHT,FREQ=99
U, 
*BOUNDARY
NALL,3
LEFT, 1
BOT,  2
BACK, 3
*PRESTRESS HOLD
ALL_ELEMENTS-2, REBAR_SET-4
ALL_ELEMENTS-2, REBAR_SET-5
ALL_ELEMENTS-2, REBAR_SET-6
ALL_ELEMENTS-2, REBAR_SET-7
*ELPRINT,ELSET=ALL_ELEMENTS,FREQ=99
S, 
**
**
** TEST CASE 1)
**
*EL PRINT,ELSET=SOLID1,REBAR=REBAR_SET-1,FREQ=99
S,E,SENER,IVOL,RBFOR
*EL FILE,ELSET=SOLID1,REBAR=REBAR_SET-1,FREQ=99
S,E,ENER,IVOL,RBFOR,NE
*EL FILE,ELSET=SFM3D4R-1,REBAR=REBAR_SET-1,FREQ=99
S,E,ENER,IVOL,RBFOR,NE
**
** TEST CASE 2)
**
*EL PRINT,ELSET=MEMB1,REBAR=REBAR_SET-2,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB1,REBAR=REBAR_SET-2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 3)
**
*EL PRINT,ELSET=SHELL1,REBAR=REBAR_SET-3,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL1,REBAR=REBAR_SET-3,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 4)
**
*EL PRINT,ELSET=SOLID2,REBAR=REBAR_SET-4,FREQ=99
S,E,SENER,IVOL,RBFOR
*EL FILE,ELSET=SOLID2,REBAR=REBAR_SET-4,FREQ=99
S,E,ENER,IVOL,RBFOR
*EL FILE,ELSET=SFM3D4R-2,REBAR=REBAR_SET-4,FREQ=99
S,E,ENER,IVOL,RBFOR
**
** TEST CASE 5)
**
*EL PRINT,ELSET=MEMB2,REBAR=REBAR_SET-5,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB2,REBAR=REBAR_SET-5,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 6)
**
*EL PRINT,ELSET=SHELL2,REBAR=REBAR_SET-6,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL2,REBAR=REBAR_SET-6,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
**
** TEST CASE 6)
**
*EL PRINT,ELSET=SHELL3,REBAR=REBAR_SET-7,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL3,REBAR=REBAR_SET-7,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
*END STEP