*HEADING
** TESTING INCREMENTAL DSA FORMULATION FOR AXISYMMETRIC ELEMENTS WITH TWIST (CGAX_)
**
** o AXIAL PULLING OF A BAR
** o HYPERELASTIC MATERIAL
**
** Node/element numbering scheme:
**    element test |  second and subsequent node prefix label | element label
**        cgax3                  100_/200_                         100/200
**        cgax3h                 101_/201_                         101/201
**        cgax4                  102_/202_                         102/202
**        cgax4h                 103_/203_                         103/203
**        cgax4r                 106_/206_                         106/206
**        cgax4rh                107_/207_                         107/207
**        cgax6                  108_/208_                         108/208
**        cgax6h                 109_/209_                         109/209
**        cgax8                  112_                              112
**        cgax8h                 113_                              113
**        cgax8r                 114_                              114
**        cgax8rh                115_                              115
**        cgax6m                 116_/216_                         116/216
**        cgax6mh                117_/217_                         117/217
**
** 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 = -20.
C10   = 80.
C01   = 20.
D1    = 0.0001
thick = 1.0
shape = 0.0
*DESIGN PARAMETER
C10,shape
*DSA CONTROLS, FORMULATION=INCREMENTAL,SIZING FREQUENCY=10
shape,CD,1e-06
*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
*PARAMETER SHAPE VARIATION, PARAMETER=shape
G, 1., 1., 0.
*NSET,NSET=FIXED
C,
*NSET,NSET=FIXED1
A,B
*NSET,NSET=FIXED2
F,I
**
*NSET,NSET=OUT
G,
**
** ELEMENT DEFINITIONS
**
*ELEMENT,ELSET=CGAX3,TYPE=CGAX3
*ELEMENT,ELSET=CGAX3H,TYPE=CGAX3H
*ELEMENT,ELSET=CGAX4,TYPE=CGAX4
*ELEMENT,ELSET=CGAX4H,TYPE=CGAX4H
*ELEMENT,ELSET=CGAX4R,TYPE=CGAX4R
*ELEMENT,ELSET=CGAX4RH,TYPE=CGAX4RH
*ELEMENT,ELSET=CGAX6,TYPE=CGAX6
*ELEMENT,ELSET=CGAX6H,TYPE=CGAX6H
*ELEMENT,ELSET=CGAX8,TYPE=CGAX8
*ELEMENT,ELSET=CGAX8H,TYPE=CGAX8H
*ELEMENT,ELSET=CGAX8R,TYPE=CGAX8R
*ELEMENT,ELSET=CGAX8RH,TYPE=CGAX8RH
*ELEMENT, ELSET=CGAX6M, TYPE=CGAX6M
*ELEMENT, ELSET=CGAX6MH, TYPE=CGAX6MH
*ELSET,ELSET=ALL
CGAX3,CGAX3H,CGAX4,CGAX4H,CGAX4R,CGAX4RH,CGAX6,CGAX6H,
CGAX8,CGAX8H,CGAX8R,CGAX8RH,CGAX6M,CGAX6MH
*ELSET,ELSET=LOADP3
*ELSET,ELSET=LOADP4
**
*SOLID SECTION,ELSET=ALL,MATERIAL=MAT
<thick>,
**
*MATERIAL,NAME=MAT
*HYPERELASTIC,N=1
<C10>,<C01>,<D1>
**
*BOUNDARY
FIXED,1,2
FIXED1,1,1
FIXED2,2,2
FIXED,5,5
FIXED1,5,5
**
** 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.
**
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 1
*STATIC
0.1, 0.4
**
*DLOAD
LOADP3,P3,<load>
LOADP4,P4,<load>
**
*OUTPUT,FIELD,FREQ=2
*NODE OUTPUT,NSET=G
U,
*ELEMENT OUTPUT
LE
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE,NSET=G
U,
*ELEMENT RESPONSE
LE
*END STEP