*HEADING
: TESTING TOTAL DSA FORMULATION FOR GENERALIZED PLANE STRAIN ELEMENTS (CPEG_)
**
** 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
**        cpeg3                  100_/200_                         100/200
**        cpeg3h                 101_/201_                         101/201
**        cpeg4                  102_/202_                         102/202
**        cpeg4h                 103_/203_                         103/203
**        cpeg4i                 104_/204_                         104/204
**        cpeg4ih                105_/205_                         105/205
**        cpeg4r                 106_/206_                         106/206
**        cpeg4rh                107_/207_                         107/207
**        cpeg6                  108_/208_                         108/208
**        cpeg6h                 109_/209_                         109/209
**        cpeg8                  112_                              112
**        cpeg8h                 113_                              113
**        cpeg8r                 114_                              114
**        cpeg8rh                115_                              115
**        cpeg6m                 116_                              116/216
**        cpeg6mh                117_                              117/217
**
** The "_" is replaced by the node number for an element.  All models
** are two elements except for the 10 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
theta1= 0.0
theta2= 0.0
shape = 0.0
*DESIGN PARAMETER
nu,shape
*DSA CONTROLS, FORMULATION=TOTAL
*NODE,NSET=A
10001, 0.0, 1.0
10101, 0.0, 1.0
10201, 0.0, 1.0
10301, 0.0, 1.0
10401, 0.0, 1.0
10501, 0.0, 1.0
10601, 0.0, 1.0
10701, 0.0, 1.0
10801, 0.0, 1.0
10901, 0.0, 1.0
11201, 0.0, 1.0
11301, 0.0, 1.0
11401, 0.0, 1.0
11501, 0.0, 1.0
11601, 0.0, 1.0
11701, 0.0, 1.0
*NODE,NSET=B
10804, 0.0, 0.5
10904, 0.0, 0.5
11205, 0.0, 0.5
11305, 0.0, 0.5
11405, 0.0, 0.5
11505, 0.0, 0.5
11602, 0.0, 0.5
11702, 0.0, 0.5
*NODE,NSET=C
10002, 0.0, 0.0
10102, 0.0, 0.0
10202, 0.0, 0.0
10302, 0.0, 0.0
10402, 0.0, 0.0
10502, 0.0, 0.0
10602, 0.0, 0.0
10702, 0.0, 0.0
10802, 0.0, 0.0
10902, 0.0, 0.0
11202, 0.0, 0.0
11302, 0.0, 0.0
11402, 0.0, 0.0
11502, 0.0, 0.0
11603, 0.0, 0.0
11703, 0.0, 0.0
*NODE,NSET=D
10204, 0.5, 1.0
10304, 0.5, 1.0
10404, 0.5, 1.0
10504, 0.5, 1.0
10604, 0.5, 1.0
10704, 0.5, 1.0
20806, 0.5, 1.0
20906, 0.5, 1.0
11208, 0.5, 1.0
11308, 0.5, 1.0
11408, 0.5, 1.0
11508, 0.5, 1.0
11604, 0.5, 1.0
11704, 0.5, 1.0
*NSET,NSET=DEQ
20806,
20906,
*NODE,NSET=E
10806, 0.5, 0.5
10906, 0.5, 0.5
11605, 0.5, 0.5
11705, 0.5, 0.5
*NODE,NSET=F
10203, 0.5, 0.0
10303, 0.5, 0.0
10403, 0.5, 0.0
10503, 0.5, 0.0
10603, 0.5, 0.0
10703, 0.5, 0.0
10805, 0.5, 0.0
10905, 0.5, 0.0
11206, 0.5, 0.0
11306, 0.5, 0.0
11406, 0.5, 0.0
11506, 0.5, 0.0
11606, 0.5, 0.0
11706, 0.5, 0.0
*NODE,NSET=G
20003, 1.0, 1.0
20103, 1.0, 1.0
20204, 1.0, 1.0
20304, 1.0, 1.0
20404, 1.0, 1.0
20504, 1.0, 1.0
20604, 1.0, 1.0
20704, 1.0, 1.0
20803, 1.0, 1.0
20903, 1.0, 1.0
11204, 1.0, 1.0
11304, 1.0, 1.0
11404, 1.0, 1.0
11504, 1.0, 1.0
11607, 1.0, 1.0
11707, 1.0, 1.0
*NSET,NSET=GEQ
20803,
20903,
11204,
11304,
11404,
11504,
*NODE,NSET=H
20805, 1.0, 0.5
20905, 1.0, 0.5
11207, 1.0, 0.5
11307, 1.0, 0.5
11407, 1.0, 0.5
11507, 1.0, 0.5
11608, 1.0, 0.5
11708, 1.0, 0.5
*NODE,NSET=I
10003, 1.0, 0.0
10103, 1.0, 0.0
20203, 1.0, 0.0
20303, 1.0, 0.0
20403, 1.0, 0.0
20503, 1.0, 0.0
20603, 1.0, 0.0
20703, 1.0, 0.0
10803, 1.0, 0.0
10903, 1.0, 0.0
11203, 1.0, 0.0
11303, 1.0, 0.0
11403, 1.0, 0.0
11503, 1.0, 0.0
11609, 1.0, 0.0
11709, 1.0, 0.0
** 
** EXTRA NODES:
*NODE,NSET=REFND
1, -0.5,-0.5
*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=CPEG3,TYPE=CPEG3
100,10001,10002,10003
200,10001,10003,20003
*ELEMENT,ELSET=CPEG3H,TYPE=CPEG3H
101,10101,10102,10103
201,10101,10103,20103
*ELEMENT,ELSET=CPEG4,TYPE=CPEG4
102,10201,10202,10203,10204
202,10204,10203,20203,20204
*ELEMENT,ELSET=CPEG4H,TYPE=CPEG4H
103,10301,10302,10303,10304
203,10304,10303,20303,20304
*ELEMENT,ELSET=CPEG4I,TYPE=CPEG4I
104,10401,10402,10403,10404
204,10404,10403,20403,20404
*ELEMENT,ELSET=CPEG4IH,TYPE=CPEG4IH
105,10501,10502,10503,10504
205,10504,10503,20503,20504
*ELEMENT,ELSET=CPEG4R,TYPE=CPEG4R
106,10601,10602,10603,10604
206,10604,10603,20603,20604
*ELEMENT,ELSET=CPEG4RH,TYPE=CPEG4RH
107,10701,10702,10703,10704
207,10704,10703,20703,20704
*ELEMENT,ELSET=CPEG6,TYPE=CPEG6
108,10801,10802,10803,10804,10805,10806
208,10801,10803,20803,10806,20805,20806
*ELEMENT,ELSET=CPEG6H,TYPE=CPEG6H
109,10901,10902,10903,10904,10905,10906
209,10901,10903,20903,10906,20905,20906
*ELEMENT,ELSET=CPEG8,TYPE=CPEG8
112,11201,11202,11203,11204,11205,11206,11207,11208
*ELEMENT,ELSET=CPEG8H,TYPE=CPEG8H
113,11301,11302,11303,11304,11305,11306,11307,11308
*ELEMENT,ELSET=CPEG8R,TYPE=CPEG8R
114,11401,11402,11403,11404,11405,11406,11407,11408
*ELEMENT,ELSET=CPEG8RH,TYPE=CPEG8RH
115,11501,11502,11503,11504,11505,11506,11507,11508
*ELEMENT,ELSET=CPEG6M,TYPE=CPEG6M
116, 11601, 11603, 11609, 11602, 11606, 11605
216, 11601, 11609, 11607, 11605, 11608, 11604
*ELEMENT,ELSET=CPEG6MH,TYPE=CPEG6MH
117, 11701, 11703, 11709, 11702, 11706, 11705
217, 11701, 11709, 11707, 11705, 11708, 11704
*ELSET,ELSET=ALL
CPEG3,CPEG3H,CPEG4,CPEG4H,CPEG4I,CPEG4IH,CPEG4R,CPEG4RH,CPEG6,CPEG6H,
CPEG8,CPEG8H,CPEG8R,CPEG8RH,CPEG6M,CPEG6MH
*ELSET,ELSET=LOADP3
200,201,208,209,116,117,216,217
*ELSET,ELSET=LOADP4
102,202,103,203,104,204,105,205,106,206,107,207,
112,113,114,115
**
*SOLID SECTION,ELSET=ALL,MATERIAL=MAT,REFNODE=REFND
<thick>,<theta1>,<theta2>
**
*MATERIAL,NAME=MAT
*ELASTIC
<young>,<nu>
**
*BOUNDARY
REFND,3,5
FIXED,1,2
FIXED1,1,1
FIXED2,2,2
**
** KINEMATIC CONSTRAINTS (TO MAINTAIN UNIFORM DEFORMATION)
**
*NSET,NSET=EU,UNSORTED
10806,10906,11605,11705
*NSET,NSET=DEQU,UNSORTED
20806,20906,11604,11704
*NSET,NSET=FU,UNSORTED
10203,10303,10403,10503,10603,10703,10805,
10905,11206,11306,11406,11506,11606,11706
*NSET,NSET=DU,UNSORTED
10204,10304,10404,10504,10604,10704,20806,
20906,11208,11308,11408,11508,11604,11704
*NSET,NSET=HU,UNSORTED
20805,20905,11207,11307,11407,11507,11608,
11708
*NSET,NSET=GEQU,UNSORTED
20803,20903,11204,11304,11404,11504,11607,
11707
*NSET,NSET=IU,UNSORTED
10003,10103,20203,20303,20403,20503,20603,20703,
10803,10903,11203,11303,11403,11503,11609,11709
*NSET,NSET=GU,UNSORTED
20003,20103,20204,20304,20404,20504,20604,20704,
20803,20903,11204,11304,11404,11504,11607,11707
*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,E
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE,NSET=G
U,
*NODE RESPONSE,NSET=C
RF,
*ELEMENT RESPONSE
S,E
*END STEP