*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-1: ELEMENT S4
 ELASTIC LARGE DEFLECTION RESPONSE OF Z-SHAPED CANTILEVER UNDER END LOAD
*SYSTEM 
0,0,0,   1,1,0  
 -1,1,0 
*NODE   
1,    0.,  0., 0.     
25,   60.,  0., 0.    
49, 120.,  0., 30. 
73, 180.,  0., 30. 

101,   0., 20., 0.
125,  60., 20., 0.
149, 120., 20., 30.
173, 180., 20., 30. 
*NGEN   
1,25,1   
25,49,1  
49,73,1
101,125,1   
125,149,1  
149,173,1
*ELEMENT, TYPE = S4   
1, 1, 2, 102, 101
*ELGEN, ELSET=CANTE     
1, 72,1
*ELSET,ELSET=ELTIPA
24,25
*SHELL GENERAL SECTION, ELSET=CANTE, POISSON=0.3, MATERIAL=ELAS, ORIENT=OR
1.7,
*MATERIAL, NAME=ELAS
*ELASTIC
2.0E5,
*ORIENTATION, NAME=OR, SYSTEM=RECT
1,1,0,-1,1,0
3,0
*NSET,NSET=END
1,101  
*BOUNDARY       
END,ENCASTRE 
*RESTART,WRITE,OVERLAY     
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,52.25
173,3,52.25
*CONTROLS, PARAMETERS=FIELD

,0.0,
*NSET,NSET=NOUT
73,173  
*NODE FILE,NSET=NOUT
CF   ,    
U    ,    
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NOUT
CF,
U,
*EL FILE,ELSET=ELTIPA,POSITION=AVERAGED AT NODES,FREQ=1     
SF  , 
*NODE PRINT,F=0
*EL PRINT,F=0
*NODE PRINT,NSET=NOUT,F=99
CF3  ,      ,     
U3   ,   
*nodeprint, nset=END, f=99
RF, 
*EL PRINT,POSITION=AVERAGED AT NODES,FREQ=99     
SM1, 
*ENDSTEP
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,631.5     
173,3,631.5 
*ENDSTEP
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,2000       
173,3,2000
*ENDSTEP