*HEADING
 MODIFIED DSA NAFEMS BENCHMARK EXAMPLE 3DNLG-1: ELEMENT B33
 ELASTIC-PLASTIC 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. 
*System
*NGEN   
1,25,1   
25,49,1  
49,73,1 
*ELEMENT, TYPE = B33
1,1,2 
*ELGEN, ELSET=CANTE     
1, 72, 1
*BEAM SECTION, SECT=RECT ,ELSET=CANTE,POISSON=0.3,MATERIAL=ABAQUS
20,<T>
-1,1,0
5,11
*MATERIAL, NAME=ABAQUS
*ELASTIC
2.0E5,  0.3
*Plastic
4000.0, 0.0
8000.0, 0.05
9000.0, 0.10
*RESTART,WRITE,OVERLAY
**______________________
*Parameter
T=1.7
P=50000.0
P2=0.25*P
Theta=0.0
*Design Parameter
T,P,Theta
***Dsa Controls, Formulation=Total
*Nset, Nset=ND
73
*Elset, Elset=ED
25
*Nset, Nset=All, Generate
2, 73, 1
*Parameter Shape Variation, Parameter=Theta
26, 0.0, 0.0,  3.00 
27, 0.0, 0.0,  6.00 
28, 0.0, 0.0,  9.00 
29, 0.0, 0.0, 12.00 
30, 0.0, 0.0, 15.00 
31, 0.0, 0.0, 18.00 
32, 0.0, 0.0, 21.00 
33, 0.0, 0.0, 24.00 
34, 0.0, 0.0, 27.00 
35, 0.0, 0.0, 30.00 
36, 0.0, 0.0, 33.00 
37, 0.0, 0.0, 36.00 
38, 0.0, 0.0, 39.00 
39, 0.0, 0.0, 42.00 
40, 0.0, 0.0, 45.00 
41, 0.0, 0.0, 48.00 
42, 0.0, 0.0, 51.00 
43, 0.0, 0.0, 54.00 
44, 0.0, 0.0, 57.00 
45, 0.0, 0.0, 60.00 
46, 0.0, 0.0, 63.00 
47, 0.0, 0.0, 66.00 
48, 0.0, 0.0, 69.00 
49, 0.0, 0.0, 72.00 
50, 0.0, 0.0, 72.00
51, 0.0, 0.0, 72.00
52, 0.0, 0.0, 72.00
53, 0.0, 0.0, 72.00
54, 0.0, 0.0, 72.00
55, 0.0, 0.0, 72.00
56, 0.0, 0.0, 72.00
57, 0.0, 0.0, 72.00
58, 0.0, 0.0, 72.00
59, 0.0, 0.0, 72.00
60, 0.0, 0.0, 72.00
61, 0.0, 0.0, 72.00
62, 0.0, 0.0, 72.00
63, 0.0, 0.0, 72.00
64, 0.0, 0.0, 72.00
65, 0.0, 0.0, 72.00
66, 0.0, 0.0, 72.00
67, 0.0, 0.0, 72.00
68, 0.0, 0.0, 72.00
69, 0.0, 0.0, 72.00
70, 0.0, 0.0, 72.00
71, 0.0, 0.0, 72.00
72, 0.0, 0.0, 72.00
73, 0.0, 0.0, 72.00
**______________________     
*BOUNDARY       
1,ENCASTRE 
**_______LOAD_25%___1_
*STEP,NLGEOM,INC=4000,DSA     
*STATIC 
0.000625,1,,0.1  
*CLOAD  
73,1,<P2>
73,2,<P2>
*OUTPUT,FIELD, Freq=1
*NODE OUTPUT, Nset=ND
U,
*ELEMENT OUTPUT, Elset=ED
SF,
*NODE PRINT,F=0
*EL PRINT,F=10
PEEQ,
*DEsign Response, Freq=1
*Node Response, Nset=ND
U,
*Element Response,ELSET=ED
SF,
*Controls, Parameters=field, field=displacement
1e-05
*ENDSTEP
**_______LOAD_100%___2_
*STEP     
*STATIC 
0.000625,1,,0.1  
*CLOAD  
73,1,<P>
73,2,<P>
*Controls, Parameters=field, field=displacement
1e-07
*ENDSTEP
**____UNLOAD_75%_____3_
*STEP
*STATIC
0.01,1.0,,0.1
*CLOAD
73,1,<P2>
73,2,<P2>
*ENDSTEP
**____UNLOAD_100%_____4_
*STEP
*STATIC
0.01,1.0,,0.1
*CLOAD,OP=NEW
*ENDSTEP