*HEADING
 NAFEMS BENCHMARK EXAMPLE 3DNLG-1: ELEMENT B31H
 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. 
*NGEN   
1,25,1   
25,49,1  
49,73,1 
*ELEMENT, TYPE = B31H
1,1,2 
*ELGEN, ELSET=CANTE     
1, 72, 1
*BEAM GENERAL SECTION, SECT=RECT ,ELSET=CANTE,POISSON=0.3
20,1.7  
-1,1,0   
2.0E5,  
*BOUNDARY       
1,ENCASTRE 
*RESTART,WRITE,OVERLAY     
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,104.5
*CONTROLS, PARAMETERS=FIELD

,0.0,
*NSET,NSET=N19
73   ,  
*NODE FILE,NSET=N19
CF   ,      ,     
U    ,    
*OUTPUT,FIELD
*NODE OUTPUT,NSET=N19
CF,,
U,
*EL FILE,POSITION=AVERAGED AT NODES,FREQ=1     
SF   ,    
*NODE PRINT,F=0
*EL PRINT,F=0
*NODE PRINT,NSET=N19,F=99
CF3  ,      ,     
U3   ,    
*EL PRINT,POSITION=AVERAGED AT NODES,FREQ=99     
SM1, 
*ENDSTEP
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,1263       
*ENDSTEP
*STEP,NLGEOM,INC=40     
*STATIC 
0.01,1  
*CLOAD  
73,3,4000       
*ENDSTEP