*HEADING
Template file for pinched cylinder parametric study     
Principal parameters: Element type (eltype), Mesh density (m_density)
Tables: nnum for node numbers of main elements
        egen for data lines on *ELGEN			
*RESTART,WRITE
*NODE,NSET=CORNER
1, 300.,0.,0., 1.,0.,0.
21, 300.,90.,0., 1.,0.,0.
2001, 300.,0.,300., 1.,0.,0.
2021, 300.,90.,300., 1.,0.,0.
*NGEN,NSET=LOADB
1,21
*NGEN,NSET=BND3
2001,2021
*NGEN,NSET=BND2
1,2001,100
*NGEN,NSET=BND4
21,2021,100
*NFILL,NSET=ALL
BND2,BND4,20
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
0.,0.,0., 300.,0.,0.
0.,300.,0.

*PARAMETER
m_density=10
eltype='S3R'
# parameters toggle1 and toggle2 are used to switch between triangular
# and quadrilateral elements in the element connectivity definition
toggle1=''
toggle2='**'
*PARAMETER DEPENDENCE,TABLE=nnum,NUMBERV=9
1,5,405,401,3,205,403,201,5
1,3,203,201,2,103,202,101,10
1,2,102,101,20,20,20,20,20
*PARAMETER,TABLE=NNUM,DEPENDENT=(n1,n2,n3,n4,n5,n6,n7,n8),
INDEPENDENT=(M_DENSITY)
*PARAMETER
# make sure tabularly dependent parameters are converted to integers
# (in general, table evaluation produces reals)
in1=int(n1)
in2=int(n2)
in3=int(n3)
in4=int(n4)
in5=int(n5)
in6=int(n6)
in7=int(n7)
in8=int(n8)
*ELEMENT,TYPE=<eltype>
1,<in1>,<in2>,<in3>,<in4>,<in5>,<in6>,<in7>,<in8>
<toggle1>2,<in1>,<in3>,<in4>
*PARAMETER
il1=20//m_density
il2=il1*100
*ELGEN,ELSET=EALL
<toggle1>1,<m_density>,<il1>,2,<m_density>,<il2>,100
<toggle1>2,<m_density>,<il1>,2,<m_density>,<il2>,100
<toggle2>1,<m_density>,<il1>,1,<m_density>,<il2>,100
*ELSET,ELSET=PRINT
 1, 
*TRANSFORM,TYPE=C,NSET=ALL
.0,.0,.0,300,.0,.0
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT
3.0, 
*MATERIAL,NAME=MAT
*ELASTIC
3.0E6,0.3
*BOUNDARY
BND2,YSYMM
BND3,ZSYMM
BND4,YSYMM
LOADB,2
LOADB,6
LOADB,1
*STEP
*STATIC
*CLOAD
2001,1,-0.25
*EL PRINT,ELSET=PRINT
*EL FILE,ELSET=PRINT
S,E
*NODE FILE,NSET=CORNER
U,RF
*END STEP