*HEADING
 o TESTING INCREMENTAL DSA FORMULATION FOR PLANE STRAIN ELEMENTS (CPE_)
**
** o AXIAL PULLING OF A BAR
** o HYPERELASTIC MATERIAL
**
** Node/element numbering scheme:
**    element test |  second and subsequent node prefix label | element label
**        cpe3                  100_/200_                         100/200
**        cpe3h                 101_/201_                         101/201
**        cpe4                  102_/202_                         102/202
**        cpe4h                 103_/203_                         103/203
**        cpe4i                 104_/204_                         104/204
**        cpe4ih                105_/205_                         105/205
**        cpe4r                 106_/206_                         106/206
**        cpe4rh                107_/207_                         107/207
**        cpe6                  108_/208_                         108/208
**        cpe6h                 109_/209_                         109/209
**        cpe6m                 110_/210_                         110/210
**        cpe6mh                111_/211_                         111/211
**        cpe8                  112_                              112
**        cpe8h                 113_                              113
**        cpe8r                 114_                              114
**        cpe8rh                115_                              115
**
** 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.
C10   = 80.
C01   = 20.
D1    = 0.0001
thick = 1.0
shape = 0.0
*DESIGN PARAMETER
C10,shape
*DSA CONTROLS, FORMULATION=INCREMENTAL
C10,CD,1.2e-05
SHAPE,CD,1.5e-06
*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
*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
*NODE,NSET=C
1002, 0.0, 0.0
1012, 0.0, 0.0
1022, 0.0, 0.0
1032, 0.0, 0.0
1042, 0.0, 0.0
1052, 0.0, 0.0
1062, 0.0, 0.0
1072, 0.0, 0.0
1082, 0.0, 0.0
1092, 0.0, 0.0
1102, 0.0, 0.0
1112, 0.0, 0.0
1122, 0.0, 0.0
1132, 0.0, 0.0
1142, 0.0, 0.0
1152, 0.0, 0.0
*NODE,NSET=D
1024, 0.5, 1.0
1034, 0.5, 1.0
1044, 0.5, 1.0
1054, 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
2116, 0.5, 1.0
1128, 0.5, 1.0
1138, 0.5, 1.0
1148, 0.5, 1.0
1158, 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
1116, 0.5, 0.5
*NODE,NSET=F
1023, 0.5, 0.0
1033, 0.5, 0.0
1043, 0.5, 0.0
1053, 0.5, 0.0
1063, 0.5, 0.0
1073, 0.5, 0.0
1085, 0.5, 0.0
1095, 0.5, 0.0
1105, 0.5, 0.0
1115, 0.5, 0.0
1126, 0.5, 0.0
1136, 0.5, 0.0
1146, 0.5, 0.0
1156, 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
2044, 1.0, 1.0
2054, 1.0, 1.0
2064, 1.0, 1.0
2074, 1.0, 1.0
2083, 1.0, 1.0
2093, 1.0, 1.0
2103, 1.0, 1.0
2113, 1.0, 1.0
1124, 1.0, 1.0
1134, 1.0, 1.0
1144, 1.0, 1.0
1154, 1.0, 1.0
*NSET,NSET=GEQ
2083,
2093,
2103,
2113,
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
*NODE,NSET=I
1003, 1.0, 0.0
1013, 1.0, 0.0
2023, 1.0, 0.0
2033, 1.0, 0.0
2043, 1.0, 0.0
2053, 1.0, 0.0
2063, 1.0, 0.0
2073, 1.0, 0.0
1083, 1.0, 0.0
1093, 1.0, 0.0
1103, 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
*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=CPE3,TYPE=CPE3
100,1001,1002,1003
200,1001,1003,2003
*ELEMENT,ELSET=CPE3H,TYPE=CPE3H
101,1011,1012,1013
201,1011,1013,2013
*ELEMENT,ELSET=CPE4,TYPE=CPE4
102,1021,1022,1023,1024
202,1024,1023,2023,2024
*ELEMENT,ELSET=CPE4H,TYPE=CPE4H
103,1031,1032,1033,1034
203,1034,1033,2033,2034
*ELEMENT,ELSET=CPE4I,TYPE=CPE4I
104,1041,1042,1043,1044
204,1044,1043,2043,2044
*ELEMENT,ELSET=CPE4IH,TYPE=CPE4IH
105,1051,1052,1053,1054
205,1054,1053,2053,2054
*ELEMENT,ELSET=CPE4R,TYPE=CPE4R
106,1061,1062,1063,1064
206,1064,1063,2063,2064
*ELEMENT,ELSET=CPE4RH,TYPE=CPE4RH
107,1071,1072,1073,1074
207,1074,1073,2073,2074
*ELEMENT,ELSET=CPE6,TYPE=CPE6
108,1081,1082,1083,1084,1085,1086
208,1081,1083,2083,1086,2085,2086
*ELEMENT,ELSET=CPE6H,TYPE=CPE6H
109,1091,1092,1093,1094,1095,1096
209,1091,1093,2093,1096,2095,2096
*ELEMENT,ELSET=CPE6M,TYPE=CPE6M
110,1101,1102,1103,1104,1105,1106
210,1101,1103,2103,1106,2105,2106
*ELEMENT,ELSET=CPE6MH,TYPE=CPE6MH
111,1111,1112,1113,1114,1115,1116
211,1111,1113,2113,1116,2115,2116
*ELEMENT,ELSET=CPE8,TYPE=CPE8
112,1121,1122,1123,1124,1125,1126,1127,1128
*ELEMENT,ELSET=CPE8H,TYPE=CPE8H
113,1131,1132,1133,1134,1135,1136,1137,1138
*ELEMENT,ELSET=CPE8R,TYPE=CPE8R
114,1141,1142,1143,1144,1145,1146,1147,1148
*ELEMENT,ELSET=CPE8RH,TYPE=CPE8RH
115,1151,1152,1153,1154,1155,1156,1157,1158
*ELSET,ELSET=ALL
CPE3,CPE3H,CPE4,CPE4H,CPE4I,CPE4IH,CPE4R,CPE4RH,CPE6,CPE6H,CPE6M,CPE6MH
CPE8,CPE8H,CPE8R,CPE8RH
*ELSET,ELSET=LOADP3
200,201,208,209,210,211
*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
*HYPERELASTIC,N=1
<C10>,<C01>,<D1>
**
*BOUNDARY
FIXED,1,2
FIXED1,1,1
FIXED2,2,2
**
** KINEMATIC CONSTRAINTS (TO MAINTAIN UNIFORM DEFORMATION)
**
*NSET,NSET=EU,UNSORTED
1086,1096,1106,1116
*NSET,NSET=DEQU,UNSORTED
2086,2096,2106,2116
*NSET,NSET=FU,UNSORTED
1023,1033,1043,1053,1063,1073,1085,
1095,1105,1115,1126,1136,1146,1156
*NSET,NSET=DU,UNSORTED
1024,1034,1044,1054,1064,1074,2086,
2096,2106,2116,1128,1138,1148,1158
*NSET,NSET=HU,UNSORTED
2085,2095,2105,2115,1127,1137,1147,1157
*NSET,NSET=GEQU,UNSORTED
2083,2093,2103,2113,1124,1134,1144,1154
*NSET,NSET=IU,UNSORTED
1003,1013,2023,2033,2043,2053,2063,2073,
1083,1093,1103,1113,1123,1133,1143,1153
*NSET,NSET=GU,UNSORTED
2003,2013,2024,2034,2044,2054,2064,2074,
2083,2093,2103,2113,1124,1134,1144,1154
*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.
**
*DIAGNOSTICS,NONHYBRID=WARNING
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 1
*STATIC
0.1, 0.4
**
*DLOAD
LOADP3,P3,<load>
LOADP4,P4,<load>
**
*OUTPUT,FIELD,FREQ=0
*END STEP
**
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 2
*STATIC
0.1, 0.4
*OUTPUT,FIELD,FREQ=0
*END STEP
**
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 3
*STATIC
0.1, 0.4
*DLOAD
LOADP3,P3,-20.0
LOADP4,P4,-20.0
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE,NSET=G
U
*ELEMENT RESPONSE
S,LE,SINV,SP
*END STEP
*STEP,NLGEOM,PERT,DSA,UNSYMM=YES
PERT STEP
*STATIC
*DLOAD
LOADP3,P3,-10.0
LOADP4,P4,-10.0
*DESIGN RESPONSE,FREQ=2
*NODE RESPONSE,NSET=G
U,
*ELEMENT RESPONSE
S,E,SINV,SP
*END STEP