*HEADING DSA MODIFIED NAFEMS BENCHMARK EXAMPLE 3DNLG-1: ELEMENT S4R 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. 101, 0., 20., 0. 125, 60., 20., 0. 149, 120., 20., 30. 173, 180., 20., 30. *SYSTEM *NGEN 1,25,1 25,49,1 49,73,1 101,125,1 125,149,1 149,173,1 *ELEMENT, TYPE = S4R 1, 1, 2, 102, 101 *ELGEN, ELSET=CANTE 1, 72,1 *ELSET,ELSET=ELTIPA 24,25 *SHELL SECTION, ELSET=CANTE, POISSON=0.3, MATERIAL=ELAS, ORIENT=OR ,11 *MATERIAL, NAME=ELAS *ELASTIC 2.0E5, *Plastic 4000.0, 0.0 8000.0, 0.05 9000.0, 0.10 *ORIENTATION, NAME=OR, SYSTEM=RECT 1,1,0,-1,1,0 3,0 *NSET,NSET=END 1,101 *BOUNDARY END,ENCASTRE *RESTART,WRITE,OVERLAY **_________________________ *Parameter T=1.7 P=50000.0 P2=P/2.0 P3=0.25*P2 Theta=0.0 *Design Parameter T,Theta,P ***Dsa Controls,Formulation=Total *Nset, Nset=ND 73 *Elset, Elset=ED 25 *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 126,0.0, 0.0, 3.00 127,0.0, 0.0, 6.00 128,0.0, 0.0, 9.00 129,0.0, 0.0, 12.00 130,0.0, 0.0, 15.00 131,0.0, 0.0, 18.00 132,0.0, 0.0, 21.00 133,0.0, 0.0, 24.00 134,0.0, 0.0, 27.00 135,0.0, 0.0, 30.00 136,0.0, 0.0, 33.00 137,0.0, 0.0, 36.00 138,0.0, 0.0, 39.00 139,0.0, 0.0, 42.00 140,0.0, 0.0, 45.00 141,0.0, 0.0, 48.00 142,0.0, 0.0, 51.00 143,0.0, 0.0, 54.00 144,0.0, 0.0, 57.00 145,0.0, 0.0, 60.00 146,0.0, 0.0, 63.00 147,0.0, 0.0, 66.00 148,0.0, 0.0, 69.00 149,0.0, 0.0, 72.00 150,0.0, 0.0, 72.00 151,0.0, 0.0, 72.00 152,0.0, 0.0, 72.00 153,0.0, 0.0, 72.00 154,0.0, 0.0, 72.00 155,0.0, 0.0, 72.00 156,0.0, 0.0, 72.00 157,0.0, 0.0, 72.00 158,0.0, 0.0, 72.00 159,0.0, 0.0, 72.00 160,0.0, 0.0, 72.00 161,0.0, 0.0, 72.00 162,0.0, 0.0, 72.00 163,0.0, 0.0, 72.00 164,0.0, 0.0, 72.00 165,0.0, 0.0, 72.00 166,0.0, 0.0, 72.00 167,0.0, 0.0, 72.00 168,0.0, 0.0, 72.00 169,0.0, 0.0, 72.00 170,0.0, 0.0, 72.00 171,0.0, 0.0, 72.00 172,0.0, 0.0, 72.00 173,0.0, 0.0, 72.00 **_____LOAD__25%___1______ *STEP,NLGEOM,INC=4000,DSA=YES *STATIC 0.01,1,,0.1 *CLOAD 73,1, 73,2, 173,1, 173,2, *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, Parameter=Field, Field=Displacement 1e-07 *ENDSTEP **_____LOAD__100%___2______ *STEP *STATIC 0.01,1,,0.1 *CLOAD 73,1, 73,2, 173,1, 173,2, *ENDSTEP **_____UNLOAD__75%___3______ *STEP *STATIC 0.01,1,,0.1 *CLOAD 73,1, 73,2, 173,1, 173,2, *ENDSTEP **____UNLOAD___100%___4_ *STEP *STATIC 0.01,1.0,,0.1 *CLOAD,OP=NEW *ENDSTEP