*heading
First order cylindrical surface elements with rebar embedded in 
cylindrical solid elements, 180 degree segment modelled.
*preprint, echo=yes, model=yes, history=yes
*node, nset=nall
1, 1.0, 0.0, 0.0
2,0.70711,0.70711,0
3, 0.0, 1.0, 0.0
7, 2.0, 0.0, 0.0
8,1.4142,1.4142,0
9, 0.0, 2.0, 0.0
21,1.0,0,1.0
22,0.70711,0.70711,1.0
23,0,1.0,1.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
105,-1.4142,1.4142,0
106,-1.4142,1.4142,1.0
107,-0.70711,0.70711,1.0
103,-1,0,1
100,-1,0,0
101,-2,0,0
102,-2,0,1
201,1.25,0,1
202,1.25,0,0
203,0,1.25,0
204,0,1.25,1
205,0.88388,0.88388,0
206,0.88388,0.88388,1
207,-1.25,0,0
208,-1.25,0,1
209,-0.88388,0.88388,0
210,-0.88388,0.88388,1
*element, type=ccl12,elset=eall
31,3,9,29,23,1,7,27,21,2,8,28,22
32,100,101,102,103,3,9,29,23,104,105,106,107
*element,type=sfmcl6,elset=mem3d
1,201,202,203,204,205,206
2,204,203,207,208,209,210
*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,21,27
*nset,nset=end2
100,101,102,103
*nset,nset=face
1,2,3,7,8,9,104,105,100,101,202,203,205,207,209
*nset,nset=top
9,29
*nset,nset=outer
7,27,8,28,9,29,105,106,101,102
*nset,nset=inner
1,21,2,22,3,23,104,107,100,103
*boundary
end1,2,2
end2,2,2
inner,1,3
*embedded element,host elset=eall,roundoff=1.0e-3
mem3d
*step, nlgeom
*static
1, 1.0
*dload
eall,p4,-5.0E3
*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