*heading
Second order cylindrical surface elements with rebar embedded in 
triangular cylindrical solid elements, 180 degree segment modelled.
*preprint, echo=yes, model=yes, history=yes
*node, nset=nall
1, 1.0, 0.0, 0.5
2,0.70711,0.70711,0.5
3, 0.0, 1.0, 0.5
7, 2.0, 0.0, 0.0
8,1.4142,1.4142,0
9, 0.0, 2.0, 0.0
27,2.0,0,1.0
28,1.4142,1.4142,1.0
29,0,2.0,1.0
104,-0.70711,0.70711,0.5
105,-1.4142,1.4142,0
106,-1.4142,1.4142,1.0
100,-1,0,0.5
101,-2,0,0
102,-2,0,1
201,1.25,0,0.625
202,1.25,0,0.375
203,0,1.25,0.375
204,0,1.25,0.625
205,0.88388,0.88388,0.375
206,0.88388,0.88388,0.625
207,-1.25,0,0.375
208,-1.25,0,0.625
209,-0.88388,0.88388,0.375
210,-0.88388,0.88388,0.625
300,1.25,0,0.5
301,0,1.25,0.5
302,0.88388,0.88388,0.5
303,-1.25,0,0.5
304,-0.88388,0.88388,0.5
305,1.5,0,0.25
306,2,0,0.5
307,1.5,0,0.75
308,0,1.5,0.25
309,0,2,0.5
310,0,1.5,0.75
311,1.060655,1.060655,0.25
312,1.4142,1.4142,0.5
313,1.060655,1.060655,0.75
314,-1.5,0,0.25
315,-2,0,0.5
316,-1.5,0,0.75
317,-1.060655,1.060655,0.25
318,-1.4142,1.4142,0.5
319,-1.060655,1.060655,0.75
*element, type=ccl18,elset=eall
31,3,9,29,1,7,27,2,8,28,308,309,310,
305,306,307,311,312,313
32,100,101,102, 3,9,29,104,105,106,314,315,316,
308,309,310,317,318,319
*element,type=sfmcl9,elset=mem3d
1,201,202,203,204,300,205,301,206,302
2,204,203,207,208,301,209,303,210,304
*surface section,elset=mem3d
*REBAR LAYER,ORI=REBAR
REBAR,0.3,0.3,0,MAT2,90.0,1
*ORIENTATION, NAME=REBAR, DEFINITION=OFFSET TO NODES
2, 4
3
*MATERIAL,NAME=MAT2
*ELASTIC
6.0E6,
*solid section, elset=eall, material=st1
*material, name=st1
*elastic
3.0E6
**
*transform, type=c, nset=nall
0.0,0.0,0.0, 0.0,0.0,1.0
*nset,nset=end1
1,7,27,305,306,307,201,202,300
*nset,nset=end2
100,101,102,314,315,316,207,208,303
*nset,nset=face
7,27,8,28,9,29,105,106,101,102,306,309,312,
315,318
*nset,nset=top
9,29,309
*nset,nset=outer
7,27,8,28,9,29,105,106,101,102,306,309,312,
315,318
*nset,nset=inner
1,2,3,104,100
*boundary
end1,2,2
end2,2,2
inner,1,3
*embedded element,host elset=eall,roundoff=1.0e-3
mem3d
*step
*static
1, 1.0
*boundary
face,1,1,0.1
*el print,elset=eall
s
e
*el print,elset=mem3d,rebar
s
e
*node print,total=yes
u
rf
*el file,elset=eall
s
e
*el file,elset=mem3d,rebar
s
e
*node file
u
rf
*output,field
*element output,elset=eall
s
e
*element output,elset=mem3d,rebar
s
e
*node output
u
rf
*end step