*HEADING
Input deck for elbow test plan for S4 elements
*restart,write
************************
*** PARAMETER DEFINITION
************************
*PARAMETER
a                 = 10.0   # length a of the straight sections
bend_radius       = 4.0    # radius of the curved section (bend)
outer_pipe_radius = 0.5    # outer radius of the pipe
wall_thick        = 0.08   # wall thickness of the pipe
num_elem_cf       = 40     # number of elements around the pipe circumference
num_elem_s        = 25     # number of elements along the straight sections of the pipe
num_elem_c        = 25     # number of elements around the bend
disp              = 1.0    # displacement at the end of the pipe
young             = 28.1E6 # Young's modulus
poisson           = 0.0    # Poisson's ratio 
#
#    total number of elements along the length of the pipe
#
num_elem = num_elem_c + 2*num_elem_s 
#
#    distance between pipe C/L and midsurface
#
midsurf_dist = outer_pipe_radius - (wall_thick/2)
#
#    geometrical properties of the problem:
#
#    Points along pipe C/L
# 
x1 = a + bend_radius
y1 = 0.
x2 = a + bend_radius
y2 = a 
x3 = a 
y3 = a + bend_radius
x4 = 0.
y4 = a + bend_radius
#
#    Points on pipe midsurface
#
x10 = x1 - midsurf_dist 
y10 = 0.
x20 = x2 - midsurf_dist
y20 = a
x30 = a
y30 = x2 - midsurf_dist
x40 = 0.
y40 = x1 - midsurf_dist
#
#
#     node number control data
#
#     ninc_cf: increment in node numbers around pipe
#     ninc   : increment in node numbers along the pipe
#           
#
ninc_cf = 1
ninc    = (num_elem_cf + 1)*ninc_cf 
#
#     Layer 1
#
n1                  = 1
#
n1_first            = n1
n1_second           = n1_first  +  ninc_cf
n1_last             = n1_first  + (num_elem_cf - 1)*ninc_cf
n1_dummy            = n1_last   +  ninc_cf
n1_first_plus_inc   = n1_first  +  ninc
n1_second_plus_inc  = n1_second +  ninc
n1_last_plus_inc    = n1_last   +  ninc
#
#     Layer 2
#
n2_first= n1_first + num_elem_s*ninc 
n2_dummy= n2_first + num_elem_cf*ninc_cf 
#
#     Layer 3
#
n3_first= n2_first + num_elem_c*ninc
n3_dummy= n3_first + num_elem_cf*ninc_cf  
#
#     Layer 4
#
n4_first= n3_first + num_elem_s*ninc
n4_dummy= n4_first + num_elem_cf*ninc_cf 
#
#     bend angle
#
angle = 90.0/num_elem_c
#
#     element number control data
#
#     einc_cf: increment in element number around the pipe
#     einc:    increment in element number along the pipe
#     num_elem_c_min_1:  number of elements around bend and circ minus one 
#     num_elem_cf_min_1: number around circ plus node inc
#     e1_first: first element number in first layer
#     e1_last:  last element number in first layer
#
einc_cf  = 1
einc     = num_elem_cf*einc_cf
#
num_elem_c_min_1     = num_elem_c  - 1 
num_elem_cf_min_1    = num_elem_cf - 1
#
e1_first = 1
e1_last  = e1_first + num_elem_cf_min_1*einc_cf
#
# END OF PARAMETER SECTION
#
#*******************
#** Node definitions  
#*******************
#*** First straight section *****
*NODE
<n1_first>,<x10>,<y10>,0.
<n1_dummy>,<x10>,<y10>,0.
*NGEN,LINE=c,NSET=circ1
<n1_first>,<n1_dummy>,<ninc_cf>,,<x1>,<y1>,0.,0.,1.,0.
*NODE
<n2_first>  , <x20>,<y20>,0.
<n2_dummy>  , <x20>,<y20>,0.
<n3_first>  , <x30>,<y30>,0.
<n3_dummy>  , <x30>,<y30>,0.
<n4_first>  , <x40>,<y40>,0.
<n4_dummy>  , <x40>,<y40>,0.
*NGEN,LINE=c,NSET=circ2
<n2_first>,<n2_dummy>,<ninc_cf>,,<x2>,<y2>,0.,0.,1.,0.
*NFILL,NSET=straight1
circ1,circ2,<num_elem_s>,<ninc>
*** Curved section **************
*NGEN,LINE=c,nset=circ3
<n3_first>,<n3_dummy>,<ninc_cf>,,<x3>,<y3>,0.,-1.,0.,0.
*NGEN,LINE=c,nset=circ4
<n4_first>,<n4_dummy>,<ninc_cf>,,<x4>,<y4>,0.,-1.,0.,0.
*NCOPY,old set=circ2, change number=<ninc>,shift, multiple=<num_elem_c_min_1>,
newset=curved
0.,0.,0.
<a>,<a>,0.,<a>,<a>,1.,<angle>
*NFILL,nset=straight2
circ3,circ4,<num_elem_s>,<ninc>
*******************
** Element generation 
*ELEMENT,TYPE=S4
<e1_first>,<n1_first>,<n1_first_plus_inc>,<n1_second_plus_inc>,<n1_second>
*ELGEN,ELSET=all_elem1
<e1_first>,<num_elem_cf_min_1>,<ninc_cf>,<einc_cf>,<num_elem>,<ninc>,<einc>
*ELEMENT,TYPE=S4
<e1_last>,<n1_last>,<n1_last_plus_inc>,<n1_first_plus_inc>,
<n1_first>,
*ELGEN,ELSET=all_elem2
<e1_last>,<num_elem>,<ninc>,<einc>
*ELSET, ELSET=all_elem
all_elem1,all_elem2
**********************
** Material definition
**********************
*MATERIAL,NAME=pipe
*ELASTIC
<young>,<poisson>
**************************
** Element definition
**************************
*SHELL SECTION,ELSET=all_elem,MATERIAL=PIPE
<wall_thick>,
**********************
*NSET,NSET=boun,gen
<n1_second>,<n1_last>,<ninc_cf>
*NSET,NSET=reac
<n1_first>,
*EQUATION
2, 
boun,1,1,<n1_first>,1,-1
*BOUNDARY
circ1,2,2
circ1,3,3
circ4,1,1
circ4,3,3
*STEP, NLGEOM
Apply displacements at the ends of the structure
*STATIC, STABILIZE
1.,1.
*BOUNDARY
reac,1,1,<disp>
circ4,2,2,<disp>
*NODE FILE, freq=9999,nset=reac
rf,
*NODE PRINT, freq=9999,nset=reac
rf,
*NODE PRINT, freq=9999,nset=circ1
rf,
*NODE PRINT, freq=9999,nset=circ4
rf,
*EL FILE, freq=0
*EL PRINT, freq=0
*END STEP