*HEADING
: TESTING TOTAL DSA FORMULATION FOR PLANE STRESS ELEMENTS (CPS_)
**
** o AXIAL PULLING OF A BAR
** o GEOMETRICALLY LINEAR AND NONLINEAR STEPS
** o DSA AND NON-DSA STEPS
** o ELASTIC MATERIAL
**
** Node/element numbering scheme:
**    element test |  second and subsequent node prefix label | element label
**        cps3                  100_/200_                         100/200
**        cps4                  102_/202_                         102/202
**        cps4i                 104_/204_                         104/204
**        cps4r                 106_/206_                         106/206
**        cps6                  108_/208_                         108/208
**        cps6m                 110_/210_                         110/210
**        cps8                  112_                              112
**        cps8r                 114_                              114
**
** The "_" is replaced by the node number for an element.  All models
** are two elements except for the 8 noded elements.  Nodes shared by 
** both elements will be labeled as belonging to the first element.
**
** Grid of nodes:
**
**       | 2
**       |
**
**       A O         D O         G O
**          
**                        
**
**       B O         E O         H O
**
**
**
**       C O         F O         I O  ---- 1
**
*PARAMETER
load = -10.
young = 100000.
nu    = 0.3
thick = 1.0
shape = 0.0
*DESIGN PARAMETER
nu,shape
*DSA CONTROLS, FORMULATION=TOTAL
*NODE,NSET=A
1001, 0.0, 1.0
1021, 0.0, 1.0
1041, 0.0, 1.0
1061, 0.0, 1.0
1081, 0.0, 1.0
1101, 0.0, 1.0
1121, 0.0, 1.0
1141, 0.0, 1.0
*NODE,NSET=B
1084, 0.0, 0.5
1104, 0.0, 0.5
1125, 0.0, 0.5
1145, 0.0, 0.5
*NODE,NSET=C
1002, 0.0, 0.0
1022, 0.0, 0.0
1042, 0.0, 0.0
1062, 0.0, 0.0
1082, 0.0, 0.0
1102, 0.0, 0.0
1122, 0.0, 0.0
1142, 0.0, 0.0
*NODE,NSET=D
1024, 0.5, 1.0
1044, 0.5, 1.0
1064, 0.5, 1.0
2086, 0.5, 1.0
2106, 0.5, 1.0
1128, 0.5, 1.0
1148, 0.5, 1.0
*NSET,NSET=DEQ
2086,
2106,
*NODE,NSET=E
1086, 0.5, 0.5
1106, 0.5, 0.5
*NODE,NSET=F
1023, 0.5, 0.0
1043, 0.5, 0.0
1063, 0.5, 0.0
1085, 0.5, 0.0
1105, 0.5, 0.0
1126, 0.5, 0.0
1146, 0.5, 0.0
*NODE,NSET=G
2003, 1.0, 1.0
2024, 1.0, 1.0
2044, 1.0, 1.0
2064, 1.0, 1.0
2083, 1.0, 1.0
2103, 1.0, 1.0
1124, 1.0, 1.0
1144, 1.0, 1.0
*NSET,NSET=GEQ
2083,
2103,
1124,
1144,
*NODE,NSET=H
2085, 1.0, 0.5
2105, 1.0, 0.5
1127, 1.0, 0.5
1147, 1.0, 0.5
*NODE,NSET=I
1003, 1.0, 0.0
2023, 1.0, 0.0
2043, 1.0, 0.0
2063, 1.0, 0.0
1083, 1.0, 0.0
1103, 1.0, 0.0
1123, 1.0, 0.0
1143, 1.0, 0.0
*PARAMETER SHAPE VARIATION, PARAMETER=shape
A, 0., 1., 0.
D, 0., 1., 0.
G, 1., 1., 0.
H, 1., 0., 0.
I, 1., 0., 0.
*NSET,NSET=FIXED
C,
*NSET,NSET=FIXED1
A,B
*NSET,NSET=FIXED2
F,I
**
** ELEMENT DEFINITIONS
**
*ELEMENT,ELSET=CPS3,TYPE=CPS3
100,1001,1002,1003
200,1001,1003,2003
*ELEMENT,ELSET=CPS4,TYPE=CPS4
102,1021,1022,1023,1024
202,1024,1023,2023,2024
*ELEMENT,ELSET=CPS4I,TYPE=CPS4I
104,1041,1042,1043,1044
204,1044,1043,2043,2044
*ELEMENT,ELSET=CPS4R,TYPE=CPS4R
106,1061,1062,1063,1064
206,1064,1063,2063,2064
*ELEMENT,ELSET=CPS6,TYPE=CPS6
108,1081,1082,1083,1084,1085,1086
208,1081,1083,2083,1086,2085,2086
*ELEMENT,ELSET=CPS6M,TYPE=CPS6M
110,1101,1102,1103,1104,1105,1106
210,1101,1103,2103,1106,2105,2106
*ELEMENT,ELSET=CPS8,TYPE=CPS8
112,1121,1122,1123,1124,1125,1126,1127,1128
*ELEMENT,ELSET=CPS8R,TYPE=CPS8R
114,1141,1142,1143,1144,1145,1146,1147,1148
*ELSET,ELSET=ALL
CPS3,CPS4,CPS4I,CPS4R,CPS6,CPS6M,CPS8,CPS8R
*ELSET,ELSET=LOADP3
200,208,210
*ELSET,ELSET=LOADP4
102,202,104,204,106,206,112,114
**
*SOLID SECTION,ELSET=ALL,MATERIAL=MAT
<thick>,
**
*MATERIAL,NAME=MAT
*ELASTIC
<young>,<nu>
**
*BOUNDARY
FIXED,1,2
FIXED1,1,1
FIXED2,2,2
**
** KINEMATIC CONSTRAINTS (TO MAINTAIN UNIFORM DEFORMATION)
**
*NSET,NSET=EU,UNSORTED
1086,1106
*NSET,NSET=DEQU,UNSORTED
2086,2106
*NSET,NSET=FU,UNSORTED
1023,1043,1063,1085,1105,1126,1146
*NSET,NSET=DU,UNSORTED
1024,1044,1064,2086,2106,1128,1148
*NSET,NSET=HU,UNSORTED
2085,2105,1127,1147
*NSET,NSET=GEQU,UNSORTED
2083,2103,1124,1144
*NSET,NSET=IU,UNSORTED
1003,2023,2043,2063,1083,1103,1123,1143
*NSET,NSET=GU,UNSORTED
2003,2024,2044,2064,2083,2103,1124,1144
*EQUATION
2, 
EU,1,1.,DEQU,1,-1.
2, 
FU,1,1.,DU,1,-1.
2, 
HU,1,1.,GEQU,1,-1.
2, 
IU,1,1.,GU,1,-1.
**
** FIRST DSA STEP
**
*STEP,DSA
STEP 1: GEOMETRICALLY LINEAR; TOTAL DSA
*STATIC
0.1, 0.4
**
*DLOAD
LOADP3,P3,<load>
LOADP4,P4,<load>
**
*OUTPUT,FIELD,FREQ=2
*NODE OUTPUT,NSET=G
U,
*NODE OUTPUT,NSET=C
RF,
*ELEMENT OUTPUT
S,E
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE, NSET=G
U,
*NODE RESPONSE, NSET=C
RF,
*ELEMENT RESPONSE
S,E
*END STEP
**
**  STATIC STEP WITH *NO* DSA
**
*STEP, DSA=NO
STEP 2: GEOMETRICALLY LINEAR; NO DSA
*STATIC
0.1, 0.4
*OUTPUT,FIELD,FREQ=2
*NODE OUTPUT,NSET=G
U,
*NODE OUTPUT,NSET=C
RF,
*ELEMENT OUTPUT
S,E
*END STEP
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 3: GEOMETRICALLY NONLINEAR; TOTAL DSA
*STATIC
0.1, 0.4
*OUTPUT,FIELD,FREQ=2
*NODE OUTPUT,NSET=G
U,
*NODE OUTPUT,NSET=C
RF,
*ELEMENT OUTPUT
S,LE
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE,NSET=G
U,
*NODE RESPONSE,NSET=C
RF,
*ELEMENT RESPONSE
S,LE
*END STEP