*HEADING
 : NAFEMS BENCHMARK EXAMPLE 3DNLG-10
 ELASTO-PLASTIC BEHAVIOUR OF A STIFFENED CYLINDRICAL PANEL UNDER COMPRESSIVE END LOAD
*NODE,NSET=PANEL,SYSTEM=C
**
** Transform panel nodes to a cylindrical system.
**
*TRANSFORM,NSET=PANEL,TYPE=C
**
** Define a new co-ordinate system to use in the definition
** of the stiffener.
**
*SYSTEM
**
** Define stiffener nodes in rectangular co-ordinates.
**
*NODE,NSET=STIFF,SYSTEM=R
**
** Transform the stiffener nodes.
**
*TRANSFORM,NSET=STIFF,TYPE=R
*ELEMENT,TYPE=STRI65,ELSET=PANEL
*ELGEN,ELSET=PANEL
*ELEMENT,TYPE=STRI65,ELSET=STIFF
*ELGEN,ELSET=STIFF
*NSET,NSET=EDGEA,GEN
*NSET,NSET=ENDAP
*NSET,NSET=ENDAS
*NSET,NSET=SYMMP
*NSET,NSET=SYMMS
*NSET,NSET=NODE3
*SHELL SECTION, ELSET=PANEL, MATERIAL=MAT1
1.0,9
*SHELL SECTION, ELSET=STIFF, MATERIAL=MAT1
1.0,9
*MATERIAL, NAME=MAT1
*ELASTIC, TYPE=ISO
210000., 0.3
*PLASTIC
350,0.0
*BOUNDARY
EDGEA,2,2,0
EDGEA,4,4,0
EDGEA,6,6,0
ENDAP,1,2,0
ENDAP,4,4,0
ENDAP,6,6,0
ENDAS,2,4,0
SYMMP,3,5,0
SYMMS,1,1,0
SYMMS,5,6,0
*ELSET, ELSET=LOAD, GENERATE
*SURFACE, NAME=SURF_LOAD
 LOAD, E2
**
**
** Initially this job is linear - until plasticity
** occurs.
** In step 1 apply load using a normal static
** step.
** In step to apply the remainder of the load
** using the RIKS algorithm
**
**
**
*RESTART,WRITE,F=10
*STEP, AMPLITUDE=RAMP, INC=10000, NLGEOM
*STATIC,RIKS
0.05, 1.,,0.1,,5,3,0.2
*MONITOR,NODE=5,DOF=3
*CONTROLS,PARAM=FIELD,FIELD=DISPLACEMENT
.01,1.
*CONTROLS,PARAM=FIELD,FIELD=ROTATION
.01,1.
*DSLOAD, CONSTANT RESULTANT=YES, FOLLOWER=NO
 SURF_LOAD, EDLD, 191., 0.,0.,1.
*OUTPUT, HISTORY
*NODE OUTPUT, NSET=NODE3
 U3,
*NODE PRINT,F=0
*EL PRINT,F=0
*NODE FILE,F=1,NSET=NODE3
U,
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NODE3
U,
*NODE FILE,F=1,NSET=ENDAP
CF,
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=ENDAP
CF,
*NODE FILE,F=1,NSET=ENDAS
CF,
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=ENDAS
CF,
*ENDSTEP