cylsh_postbuck.inp:Input file for the postbuckling analysis

*heading
Input file for the postbuckling analysis.
*node,system=c
1,5.0,0.0,0.0
240,5.0,-1.5,0.0
4801,5.0,0.0,2
5040,5.0,-1.5,2
*ngen,line=c,nset=bottom
1,240,1,,0.0,0.0,0.0,0.0,0.0,1.0
*ncopy,new set=top,old set=bottom,shift,change number=4800
0.0,0.0,2
0.0,0.0,0.0,0.0,0.0,1.0,0.0
*nfill
bottom,top,20,240
**specify the imperfection as a function of modeshape amplitude
*imperfection,file=cylsh_buck1,step=2
*include, input=max_round_input.dat
*element,type=s4
1,1,2,242,241
240,240,1,241,480
*elgen,elset=cylinder
1,239,1,1,20,240,240
240,1,,,20,240,240
*shell section, elset=cylinder, material=mat_1
0.01,
*material,name=mat_1
*elastic
30e+06,0.3
*nset,nset=ends
bottom,top
*nset,nset=node_out
5040,
*transform,type=c,nset=ends
0.0,0.0,0.0,0.0,0.0,1.0
*boundary
ends,1,2
ends,4,4
ends,6,6
1,3
**
*step,nlgeom
static preload for internal pressure
*static
1.0,1.0
*monitor,node=5040, dof=3
*dload
cylinder,p,0.0
*end step
**
*step,nlgeom,inc=50
postbuckling (riks) analysis
*static,riks
0.05,1.0,,0.05,,5040,3,-0.1
*monitor,node=5040,dof=3
*controls,parameter=field,field=global
1.e-5,
*cload
top,3,-49.29375
bottom,3,49.29375
*node file,freq=20,nset=top
u,
*output,field,freq=20
*node output
u,
*output,history,freq=1
*node output,nset=node_out
u,
*end step