*heading
First order cylindrical surface elements with rebar embedded in 
cylindrical solid elements, 180 degree segment modelled.
*preprint, echo=yes, model=yes, history=yes
*RESTART, WRITE
**
***************************************************
** Node and element definition for structure 1    *
***************************************************
**
*node, nset=nall1
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=eall1
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=mem3d1
1,201,202,203,204,205,206
2,204,203,207,208,209,210
**
***************************************************
** Node and element definition for structure 2    *
***************************************************
**
*node, nset=nall2
1001, 1.0, 0.0, 0.0
1002,0.70711,0.70711,0
1003, 0.0, 1.0, 0.0
1007, 2.0, 0.0, 0.0
1008,1.4142,1.4142,0
1009, 0.0, 2.0, 0.0
1021,1.0,0,1.0
1022,0.70711,0.70711,1.0
1023,0,1.0,1.0
1027,2.0,0,1.0
1028,1.4142,1.4142,1.0
1029,0,2.0,1.0
1104,-0.70711,0.70711,0
1105,-1.4142,1.4142,0
1106,-1.4142,1.4142,1.0
1107,-0.70711,0.70711,1.0
1103,-1,0,1
1100,-1,0,0
1101,-2,0,0
1102,-2,0,1
1201,1.25,0,1
1202,1.25,0,0
1203,0,1.25,0
1204,0,1.25,1
1205,0.88388,0.88388,0
1206,0.88388,0.88388,1
1207,-1.25,0,0
1208,-1.25,0,1
1209,-0.88388,0.88388,0
1210,-0.88388,0.88388,1
*element, type=ccl12,elset=eall2
1031,1003,1009,1029,1023,1001,1007,1027,1021,1002,1008,1028,1022
1032,1100,1101,1102,1103,1003,1009,1029,1023,1104,1105,1106,1107
*element,type=sfmcl6,elset=mem3d2
1001,1201,1202,1203,1204,1205,1206
1002,1204,1203,1207,1208,1209,1210
**
***************************************************
** Sections and material properties               *
***************************************************
**
*surface section,elset=mem3d1
*REBAR LAYER,ORI=REBAR
REBAR,0.3,0.3,0,MAT2,90.0,1
*surface section,elset=mem3d2
*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=eall1, material=st1
*solid section, elset=eall2, material=st1
*material, name=st1
*elastic
200.0E9,0.3
**
******************************************************
** Boundary conditions                               *
******************************************************
*transform, type=c, nset=nall1
0.0,0.0,0.0, 0.0,0.0,1.0
*transform, type=c, nset=nall2
0.0,0.0,0.0, 0.0,0.0,1.0
*nset,nset=end1
1,7,21,27,1001,1007,1021,1027
*nset,nset=end2
100,101,102,103,1100,1101,1102,1103
*nset,nset=inner
1,21,2,22,3,23,104,107,100,103,1001,1021,1002,1022,1003,1023,
1104,1107,1100,1103
*boundary
end1,2,2
end2,2,2
inner,1,3
*embedded element,host elset=eall1,roundoff=1.0e-3
mem3d1
*embedded element,host elset=eall2,roundoff=1.0e-3
mem3d2
**
*******************************************************
** History data                                       *
*******************************************************
*step, nlgeom
*static
1, 1.0
*dload
eall1,p4,-5.0E3
eall2,p4,-2.5E3
*el print,elset=eall1
s
e
*el print,elset=mem3d1,rebar
s
e
*node print,total=yes,nset=nall1
u
rf
*el print,elset=eall2
s
e
*el print,elset=mem3d2,rebar
s
e
*node print,total=yes,nset=nall2
u
rf
*output,field
*element output,elset=eall1
s
e
*element output,elset=mem3d1,rebar
s
e
*node output, nset=nall1
u
rf
*element output,elset=eall2
s
e
*element output,elset=mem3d2,rebar
s
e
*node output, nset=nall2
u
rf
*end step
*step, nlgeom
*static
1, 1.0
*dload
eall1,p4,-5.0E3
eall2,p4,-5.0E3
*el print,elset=eall1
s
e
*el print,elset=mem3d1,rebar
s
e
*node print,total=yes,nset=nall1
u
rf
*el print,elset=eall2
s
e
*el print,elset=mem3d2,rebar
s
e
*node print,total=yes,nset=nall2
u
rf
*output,field
*element output,elset=eall1
s
e
*element output,elset=mem3d1,rebar
s
e
*node output, nset=nall1
u
rf
*element output,elset=eall2
s
e
*element output,elset=mem3d2,rebar
s
e
*node output, nset=nall2
u
rf
*end step