*heading
Second 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
211,1.25,0,0.5
212,0,1.25,0.5
213,0.88388,0.88388,0.5
214,-1.25,0,0.5
215,-0.88388,0.88388,0.5
301,1.5,0,0
302,2,0,0.5
303,1.5,0,1
304,1,0,0.5
305,1.060655,1.060655,0
306,1.4142,1.4142,0.5
307,1.060655,1.060655,1.0
308,0.70711,0.70711,0.5
309,0,1.5,0
310,0,2,0.5
311,0,1.5,1
312,0,1,0.5
313,-1.060655,1.060655,0
314,-1.4142,1.4142,0.5
315,-1.060655,1.060655,1
316,-0.70711,0.70711,0.5
317,-1.5,0,0
318,-2,0,0.5
319,-1.5,0,1
320,-1,0,0.5
*element, type=ccl24,elset=eall
31,3,9,29,23,1,7,27,21,2,8,28,22,309,310,311,312,301,302,303,304,
305,306,307,308
32,100,101,102,103,3,9,29,23,104,105,106,107,317, 318,319,320,
309,310,311,312,313,314,315,316
*element,type=sfmcl9,elset=mem3d
1,201,202,203,204,211,205,212,206,213
2,204,203,207,208,212,209,214,210,215
*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,301,302,303,304
*nset,nset=end2
100,101,102,103,317,318,319,320
*nset,nset=face
1,2,3,7,8,9,104,105,100,101,202,203,205,207,
209,301,305,309,313,317
*nset,nset=top
9,29,310
*nset,nset=outer
7,27,302,8,28,306,9,29,310,105,106,
314,101,102,318
*nset,nset=inner
1,21,304,2,22,308,3,23,312,104,107,316,
100,103,320
*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