*HEADING
 o 2X2 AXISYMMETRIC RING; TWO ELASTIC MATERIALS
 o LINEAR GEOMETRY, TOTAL DSA FORMULATION
 o IVOL, E, U - Design responses
**
*DSA CONTROLS, FORMULATION=TOTAL
*PARAMETER
y1 = 20.E6
y1t = 20.E5
nu1 = 0.2
y2 = 30.E6
y2t = 30.E5
nu2 = 0.3
n5c1 = 0.0
n5c2 = 0.0
*DESIGN PARAMETER
y1t,y2t,n5c2
*PARAMETER SHAPE VARIATION, PARAMETER=n5c1
5, 1., 0., 0.
*NODE,NSET=NALL
 1, 1.250, -0.5
 2, 2.000, -0.5
 3, 2.750, -0.5
 4, 1.250,  0.0
 5, 2.000,  0.0
 6, 2.750,  0.0
 7, 1.250,  0.5
 8, 2.000,  0.5
 9, 2.750,  0.5
*NSET, NSET=SL
1,4,7
*NSET, NSET=SR
3,6,10
*ELEMENT, TYPE=CAX4, ELSET=E1
1,  1,  2,  5,  4
*ELEMENT, TYPE=CAX4, ELSET=E2
2,  2,  3,  6,  5
*ELEMENT, TYPE=CAX4, ELSET=E3
3,  4,  5,  8,  7
*ELEMENT, TYPE=CAX4, ELSET=E4
4,  5,  6,  9,  8
*ELSET,ELSET=ALL
1,2,3,4
*SOLID SECTION, MATERIAL=HKS1, ELSET=E1
*SOLID SECTION, MATERIAL=HKS1, ELSET=E2
*SOLID SECTION, MATERIAL=HKS2, ELSET=E3
*SOLID SECTION, MATERIAL=HKS2, ELSET=E4
*MATERIAL, NAME=HKS1
*ELASTIC
<y1>,<nu1>,0.0
<y1t>,<nu1>,100.0
*MATERIAL, NAME=HKS2
*ELASTIC
<y2>,<nu2>,0.0
<y2t>,<nu2>,100.0
*BOUNDARY
 1, ENCASTRE
 2, 2,2,0.0
 3, 2,2,0.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL, 50.0
*STEP,DSA
  STRECHING STEP
*STATIC
0.1, 0.2
*CLOAD
8,2,5.0E6
7,2,2.5E6
9,2,2.5E6
*ELPRINT, ELSET=ALL, FREQUENCY=1,TOTALS=YES
S,
E,
*NODE PRINT, FREQUENCY=1,TOTALS=YES
U,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E, IVOL
*NODE OUTPUT
U, 
*DESIGN RESPONSE
*ELEMENT RESPONSE
E,IVOL
*NODE RESPONSE
U
*END STEP