*HEADING
: TESTING TOTAL DSA FORMULATION FOR AXISYMMETRIC ELEMENTS (CAX_)
**
** o AXIAL PULLING OF A BAR
** o ELASTIC 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 = -2000.
young = 100000.
nu    = 0.3
thick = 1.0
shape = 0.0
*DESIGN PARAMETER
young,shape
*DSA CONTROLS, FORMULATION=TOTAL
*NODE,NSET=A
1001, 0.0, 1.0
1011, 0.0, 1.0
1021, 0.0, 1.0
1031, 0.0, 1.0
1041, 0.0, 1.0
1051, 0.0, 1.0
1061, 0.0, 1.0
1071, 0.0, 1.0
1081, 0.0, 1.0
1091, 0.0, 1.0
1101, 0.0, 1.0
1111, 0.0, 1.0
1121, 0.0, 1.0
1131, 0.0, 1.0
1141, 0.0, 1.0
1151, 0.0, 1.0
1161, 0.0, 1.0
1171, 0.0, 1.0
*NODE,NSET=B
1084, 0.0, 0.5
1094, 0.0, 0.5
1104, 0.0, 0.5
1114, 0.0, 0.5
1125, 0.0, 0.5
1135, 0.0, 0.5
1145, 0.0, 0.5
1155, 0.0, 0.5
1162, 0.0, 0.5
1172, 0.0, 0.5
*NODE,NSET=C
1002, 0.0, 0.0
1012, 0.0, 0.0
1022, 0.0, 0.0
1032, 0.0, 0.0
1062, 0.0, 0.0
1072, 0.0, 0.0
1082, 0.0, 0.0
1092, 0.0, 0.0
1122, 0.0, 0.0
1132, 0.0, 0.0
1142, 0.0, 0.0
1152, 0.0, 0.0
1163, 0.0, 0.0
1173, 0.0, 0.0
*NODE,NSET=D
1024, 0.5, 1.0
1034, 0.5, 1.0
1064, 0.5, 1.0
1074, 0.5, 1.0
2086, 0.5, 1.0
2096, 0.5, 1.0
2106, 0.5, 1.0
1128, 0.5, 1.0
1138, 0.5, 1.0
1148, 0.5, 1.0
1158, 0.5, 1.0
1164, 0.5, 1.0
1174, 0.5, 1.0
*NSET,NSET=DEQ
2086,
2096,
2106,
2116,
*NODE,NSET=E
1086, 0.5, 0.5
1096, 0.5, 0.5
1106, 0.5, 0.5
1165, 0.5, 0.5
1175, 0.5, 0.5
*NODE,NSET=F
1023, 0.5, 0.0
1033, 0.5, 0.0
1063, 0.5, 0.0
1073, 0.5, 0.0
1085, 0.5, 0.0
1095, 0.5, 0.0
1126, 0.5, 0.0
1136, 0.5, 0.0
1146, 0.5, 0.0
1156, 0.5, 0.0
1166, 0.5, 0.0
1176, 0.5, 0.0
*NODE,NSET=G
2003, 1.0, 1.0
2013, 1.0, 1.0
2024, 1.0, 1.0
2034, 1.0, 1.0
2064, 1.0, 1.0
2074, 1.0, 1.0
2083, 1.0, 1.0
2093, 1.0, 1.0
1124, 1.0, 1.0
1134, 1.0, 1.0
1144, 1.0, 1.0
1154, 1.0, 1.0
1167, 1.0, 1.0
1177, 1.0, 1.0
*NSET,NSET=GEQ
2083,
2093,
1124,
1134,
1144,
1154,
*NODE,NSET=H
2085, 1.0, 0.5
2095, 1.0, 0.5
2105, 1.0, 0.5
2115, 1.0, 0.5
1127, 1.0, 0.5
1137, 1.0, 0.5
1147, 1.0, 0.5
1157, 1.0, 0.5
1168, 1.0, 0.5
1178, 1.0, 0.5
*NODE,NSET=I
1003, 1.0, 0.0
1013, 1.0, 0.0
2023, 1.0, 0.0
2033, 1.0, 0.0
2063, 1.0, 0.0
2073, 1.0, 0.0
1083, 1.0, 0.0
1093, 1.0, 0.0
1113, 1.0, 0.0
1123, 1.0, 0.0
1133, 1.0, 0.0
1143, 1.0, 0.0
1153, 1.0, 0.0
1169, 1.0, 0.0
1179, 1.0, 0.0
*PARAMETER SHAPE VARIATION, PARAMETER=shape
G, 1., 1., 0.
*NSET,NSET=FIXED
C,
*NSET,NSET=FIXED1
A,B
*NSET,NSET=FIXED2
F,I
**
** ELEMENT DEFINITIONS
**
*ELEMENT,ELSET=CGAX3,TYPE=CGAX3
100,1001,1002,1003
200,1001,1003,2003
*ELEMENT,ELSET=CGAX3H,TYPE=CGAX3H
101,1011,1012,1013
201,1011,1013,2013
*ELEMENT,ELSET=CGAX4,TYPE=CGAX4
102,1021,1022,1023,1024
202,1024,1023,2023,2024
*ELEMENT,ELSET=CGAX4H,TYPE=CGAX4H
103,1031,1032,1033,1034
203,1034,1033,2033,2034
*ELEMENT,ELSET=CGAX4R,TYPE=CGAX4R
106,1061,1062,1063,1064
206,1064,1063,2063,2064
*ELEMENT,ELSET=CGAX4RH,TYPE=CGAX4RH
107,1071,1072,1073,1074
207,1074,1073,2073,2074
*ELEMENT,ELSET=CGAX6,TYPE=CGAX6
108,1081,1082,1083,1084,1085,1086
208,1081,1083,2083,1086,2085,2086
*ELEMENT,ELSET=CGAX6H,TYPE=CGAX6H
109,1091,1092,1093,1094,1095,1096
209,1091,1093,2093,1096,2095,2096
*ELEMENT,ELSET=CGAX8,TYPE=CGAX8
112,1121,1122,1123,1124,1125,1126,1127,1128
*ELEMENT,ELSET=CGAX8H,TYPE=CGAX8H
113,1131,1132,1133,1134,1135,1136,1137,1138
*ELEMENT,ELSET=CGAX8R,TYPE=CGAX8R
114,1141,1142,1143,1144,1145,1146,1147,1148
*ELEMENT,ELSET=CGAX8RH,TYPE=CGAX8RH
115,1151,1152,1153,1154,1155,1156,1157,1158
*ELEMENT, ELSET=CGAX6M, TYPE=CGAX6M
116,1161,1163,1169,1162,1166,1165
216,1161,1169,1167,1165,1168,1164
*ELEMENT, ELSET=CGAX6MH, TYPE=CGAX6MH
117,1171,1173,1179,1172,1176,1175
217,1171,1179,1177,1175,1178,1174
*ELSET,ELSET=ALL
CGAX3,CGAX3H,CGAX4,CGAX4H,CGAX4R,CGAX4RH,CGAX6,CGAX6H,
CGAX8,CGAX8H,CGAX8R,CGAX8RH,CGAX6M,CGAX6MH
*ELSET,ELSET=LOADP3
200,201,208,209,210,211,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
<thick>,
**
*MATERIAL,NAME=MAT
*ELASTIC
<young>,<nu>
**
*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
1086,1096,1165,1175
*NSET,NSET=DEQU,UNSORTED
2086,2096,1164,1174
*NSET,NSET=FU,UNSORTED
1023,1033,1063,1073,1085,
1095,1126,1136,1146,1156,
1166,1176
*NSET,NSET=DU,UNSORTED
1024,1034,1064,1074,2086,
2096,1128,1138,1148,1158,
1164,1174
*NSET,NSET=HU,UNSORTED
2085,2095,1127,1137,1147,1157,
1168,1178
*NSET,NSET=GEQU,UNSORTED
2083,2093,1124,1134,1144,1154,
1167,1177
*NSET,NSET=IU,UNSORTED
1003,1013,2023,2033,2063,2073,
1083,1093,1123,1133,1143,1153,
1169,1179
*NSET,NSET=GU,UNSORTED
2003,2013,2024,2034,2064,2074,
2083,2093,1124,1134,1144,1154,
1167,1177
*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.
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 1: GEOMETRICALLY NONLINEAR; TOTAL DSA
*STATIC
0.1, 0.4
**
*DLOAD
LOADP3,P3,<load>
LOADP4,P4,<load>
**
*CLOAD
G,5,5,<load>
*OUTPUT,FIELD,FREQ=2
*NODE OUTPUT,NSET=G
U,
*NODE OUTPUT,NSET=C
RF,
*ELEMENT OUTPUT
LE
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE, NSET=G
U,
*NODE RESPONSE, NSET=C
RF,
*ELEMENT RESPONSE
LE
*END STEP