** TEST OF REBAR: FOR Nonlinear Analyses: (NLGEOM)
**                ELASTIC BASE
**                ELASTIC REBAR
**
** ORIENTATION OF REBAR: 60deg to Xaxis
**  
*HEADING
 UNIAXIAL TEST OF REBAR 
*PREPRINT,MODEL=YES
**
** TEST CASE 1)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
111,  0., 2.
112,  1., 2.
113,  1., 3.
114,  0., 3.
*NODAL THICKNESS
LEFT,1.0
RIGHT,0.5
*NSET,NSET=BOT
111,112
*NSET,NSET=LEFT
111,114
*NSET,NSET=RIGHT
112,113
*NSET,NSET=TOP
113,114
*ELEMENT,TYPE=M3D4,ELSET=MEMB1
111, 111,112,113,114
*MEMBRANE SECTION,ELSET=MEMB1,MATERIAL=ELAS,NODAL THICKNESS
1.0,
*REBAR LAYER
REB6,0.5,1.0,0,REBAR,60,1 
REB6_2,0.5,1.0,0,REBAR,-60,1 
** TEST CASE 2)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D4R
** 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
*NSET,NSET=TOP
123,124
*NODAL THICKNESS
LEFT,1.0
RIGHT,0.5
*ELEMENT,TYPE=M3D4R,ELSET=MEMB2
121, 121,122,123,124
*MEMBRANE SECTION,ELSET=MEMB2,MATERIAL=ELAS,NODAL THICKNESS
1.0,
*REBAR LAYER
REB7,0.5,1.0,0,REBAR,60,1 
REB7_2,0.5,1.0,0,REBAR,-60,1
**
** TEST CASE 3)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D8
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
131,  4.0, 2.
132,  5.0, 2.
133,  5.0, 3.
134,  4.0, 3.
135,  4.5, 2.
136,  5.0, 2.5
137,  4.5, 3.
138,  4.0, 2.5
*NSET,NSET=BOT
131,132,135
*NSET,NSET=LEFT
131,134,138
*NSET,NSET=RIGHT
132,133,136
*NSET,NSET=TOP
133,134,137
*NODAL THICKNESS
131,1.0
132,0.5
133,0.5
134,1.0
135,0.75
136,0.5
137,0.75
138,1.0
*ELEMENT,TYPE=M3D8,ELSET=MEMB3
131, 131,132,133,134,135,136,137,138
*MEMBRANE SECTION,ELSET=MEMB3,MATERIAL=ELAS, NODAL THICKNESS
1.0,
*REBAR LAYER
REB8,0.5,1.0,0,REBAR,60,1 
REB8_2,0.5,1.0,0,REBAR,-60,1  
** TEST CASE 4)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D8R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
141,  6.0, 2.
142,  7.0, 2.
143,  7.0, 3.
144,  6.0, 3.
145,  6.5, 2.
146,  7.0, 2.5
147,  6.5, 3.
148,  6.0, 2.5
*NODAL THICKNESS
141,1.0
142,0.5
143,0.5
144,1.0
145,0.75
146,0.5
147,0.75
148,1.0
*NSET,NSET=BOT
141,142,145
*NSET,NSET=LEFT
141,144,148
*NSET,NSET=RIGHT
142,143,146
*NSET,NSET=TOP
143,144,147
*ELEMENT,TYPE=M3D8R,ELSET=MEMB4
141, 141,142,143,144,145,146,147,148
*MEMBRANE SECTION,ELSET=MEMB4,MATERIAL=ELAS, NODAL THICKNESS
1.0,
*REBAR LAYER
REB9,0.5,1.0,0,REBAR,60,1 
REB9_2,0.5,1.0,0,REBAR,-60,1
** TEST CASE 5)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: M3D9R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
151,  8.0, 2.
152,  9.0, 2.
153,  9.0, 3.
154,  8.0, 3.
155,  8.5, 2.
156,  9.0, 2.5
157,  8.5, 3.
158,  8.0, 2.5
159,  8.5, 2.5
*NSET,NSET=BOT
151,152,155
*NSET,NSET=LEFT
151,154,158
*NSET,NSET=RIGHT
152,153,156
*NSET,NSET=TOP
153,154,157
*NODAL THICKNESS
151, 1.0
152, 0.5
153, 0.5
154, 1.0
155, 0.75
156, 0.5
157, 0.75
158, 1.0
159, 0.75
*ELEMENT,TYPE=M3D9R,ELSET=MEMB5
151, 151,152,153,154,155,156,157,158,159
*MEMBRANE SECTION,ELSET=MEMB5,MATERIAL=ELAS, NODAL THICKNESS
1.0,
*REBAR LAYER
REB10,0.5,1.0,0,REBAR,60,1 
REB10_2,0.5,1.0,0,REBAR,-60,1
**
**
**
** TEST CASE 6)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4R
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
221,  2., 4.
222,  3., 4.
223,  3., 5.
224,  2., 5.
*NODAL THICKNESS
LEFT,1.0
RIGHT,0.5
*NSET,NSET=BOT
221,222
*NSET,NSET=LEFT
221,224
*NSET,NSET=RIGHT
222,223
*NSET,NSET=TOP
223,224
*ELEMENT,TYPE=S4R,ELSET=SHELL1
221, 221,222,223,224
*SHELL SECTION,ELSET=SHELL1,MATERIAL=ELAS,NODAL THICKNESS
1.0,
*REBAR LAYER
REB12,0.5,1.0,0,REBAR,60,1 
REB12_2,0.5,1.0,0,REBAR,-60,1
** TEST CASE 6)
** MESH: 1 ELEMENT 
** REBAR POSITION: MIDSURFACE
** ELEMENT TYPE: S4
** BASE MATERIAL: ELASTIC
** REBAR MATERIAL: ELASTIC
** 
**
*NODE,NSET=NALL
231,  4., 4.
232,  5., 4.
233,  5., 5.
234,  4., 5.
*NODAL THICKNESS
LEFT,1.0
RIGHT,0.5
*NSET,NSET=BOT
231,232
*NSET,NSET=LEFT
231,234
*NSET,NSET=RIGHT
232,233
*NSET,NSET=TOP
233,234
*ELEMENT,TYPE=S4,ELSET=SHELL2
231, 231,232,233,234
*SHELL SECTION,ELSET=SHELL2,MATERIAL=ELAS,NODAL THICKNESS
1.0,
*REBAR LAYER
REB13,0.5,1.0,0,REBAR,60,1 
REB13_2,0.5,1.0,0,REBAR,-60,1 
**
**
**
*MATERIAL,NAME=ELAS
*ELASTIC
1.E4, .25
*MATERIAL,NAME=REBAR
*ELASTIC
2.E4, 
*ELSET,ELSET=EALL
MEMB1,MEMB2,MEMB3,MEMB4,MEMB5,SHELL1,SHELL2
****************
** LOAD CASES **
****************
*STEP, NLGEOM, INC=50
*STATIC
0.1,1.
*RESTART,WRITE
*NODE PRINT
RF, 
*BOUNDARY
NALL,3
LEFT, 1
BOT,  2
TOP,  2
RIGHT,1,1,1.0
**
**
** TEST CASE 1)
**
*EL PRINT,ELSET=MEMB1,FREQ=99
STH, 
*EL FILE,ELSET=MEMB1,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=MEMB1
STH,
*ELEMENT OUTPUT,ELSET=MEMB1,REBAR=REB6
S,E
*ELEMENT OUTPUT,ELSET=MEMB1,REBAR=REB6,POSITION=CENTROIDAL
S,E
*EL PRINT,ELSET=MEMB1,REBAR=REB6,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB1,REBAR=REB6,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB1,REBAR=REB6_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 2)
**
*EL PRINT,ELSET=MEMB2,FREQ=99
STH, 
*EL FILE,ELSET=MEMB2,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=MEMB2
STH,
*EL PRINT,ELSET=MEMB2,REBAR=REB7,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB2,REBAR=REB7,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB2,REBAR=REB7_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 3)
**
*EL PRINT,ELSET=MEMB3,FREQ=99
STH, 
*EL FILE,ELSET=MEMB3,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=MEMB3
STH,
*EL PRINT,ELSET=MEMB3,REBAR=REB8,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB3,REBAR=REB8,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB3,REBAR=REB8_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 4)
**
*EL PRINT,ELSET=MEMB4,FREQ=99
STH, 
*EL FILE,ELSET=MEMB4,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=MEMB4
STH,
*EL PRINT,ELSET=MEMB4,REBAR=REB9,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB4,REBAR=REB9,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB4,REBAR=REB9_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 5)
**
*EL PRINT,ELSET=MEMB5,FREQ=99
STH, 
*EL FILE,ELSET=MEMB5,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=MEMB5
STH,
*EL PRINT,ELSET=MEMB5,REBAR=REB10,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB5,REBAR=REB10,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=MEMB5,REBAR=REB10_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
**
** TEST CASE 6)
**
*EL PRINT,ELSET=SHELL1,FREQ=99
STH, 
*EL FILE,ELSET=SHELL1,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=SHELL1
STH,
*EL PRINT,ELSET=SHELL1,REBAR=REB12,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL1,REBAR=REB12,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL1,REBAR=REB12_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
** TEST CASE 7)
**
*EL PRINT,ELSET=SHELL2,FREQ=99
STH, 
*EL FILE,ELSET=SHELL2,FREQ=99
STH, 
*OUTPUT,FIELD,FREQ=99
*ELEMENT OUTPUT,ELSET=SHELL2
STH,
*EL PRINT,ELSET=SHELL2,REBAR=REB13,FREQ=99
S,E,SENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL2,REBAR=REB13,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
*EL FILE,ELSET=SHELL2,REBAR=REB13_2,FREQ=99
S,E,ENER,IVOL,RBFOR,RBANG,RBROT
**
*END STEP