*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
*NODE,NSET=B
*NODE,NSET=C
*NODE,NSET=D
*NSET,NSET=DEQ
*NODE,NSET=E
*NODE,NSET=F
*NODE,NSET=G
*NSET,NSET=GEQ
*NODE,NSET=H
*NODE,NSET=I
**
** EXTRA NODES:
*NODE,NSET=REFND
*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
*ELEMENT,ELSET=CPEG3H,TYPE=CPEG3H
*ELEMENT,ELSET=CPEG4,TYPE=CPEG4
*ELEMENT,ELSET=CPEG4H,TYPE=CPEG4H
*ELEMENT,ELSET=CPEG4I,TYPE=CPEG4I
*ELEMENT,ELSET=CPEG4IH,TYPE=CPEG4IH
*ELEMENT,ELSET=CPEG4R,TYPE=CPEG4R
*ELEMENT,ELSET=CPEG4RH,TYPE=CPEG4RH
*ELEMENT,ELSET=CPEG6,TYPE=CPEG6
*ELEMENT,ELSET=CPEG6H,TYPE=CPEG6H
*ELEMENT,ELSET=CPEG8,TYPE=CPEG8
*ELEMENT,ELSET=CPEG8H,TYPE=CPEG8H
*ELEMENT,ELSET=CPEG8R,TYPE=CPEG8R
*ELEMENT,ELSET=CPEG8RH,TYPE=CPEG8RH
*ELEMENT,ELSET=CPEG6M,TYPE=CPEG6M
*ELEMENT,ELSET=CPEG6MH,TYPE=CPEG6MH
*ELSET,ELSET=ALL
CPEG3,CPEG3H,CPEG4,CPEG4H,CPEG4I,CPEG4IH,CPEG4R,CPEG4RH,CPEG6,CPEG6H,
CPEG8,CPEG8H,CPEG8R,CPEG8RH,CPEG6M,CPEG6MH
*ELSET,ELSET=LOADP3
*ELSET,ELSET=LOADP4
**
*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
*NSET,NSET=DEQU,UNSORTED
*NSET,NSET=FU,UNSORTED
*NSET,NSET=DU,UNSORTED
*NSET,NSET=HU,UNSORTED
*NSET,NSET=GEQU,UNSORTED
*NSET,NSET=IU,UNSORTED
*NSET,NSET=GU,UNSORTED
*EQUATION
EU,1,1.,DEQU,1,-1.
FU,1,1.,DU,1,-1.
HU,1,1.,GEQU,1,-1.
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