*Heading
 o TESTING DSA FOR BRICK ELEMENTS: C3D8
 o 2X2 SOLID; TWO ELASTIC MATERIALS
 o LINEAR GEOMETRY, TOTAL DSA FORMULATION
*DESIGN PARAMETER
E1,N1,E2,S2
*PARAMETER
E1=2.E5
N1=0.45
E2=3.E4
N2=0.1
S2=0.0
*Dsa controls, FORMULATION=TOTAL
N1, CD, 1.0e-05
*Node
      1,           5.,          -5.,           0.
      2,           5.,           0.,           0.
      3,           5.,           5.,           0.
      4,           5.,          -5.,           5.
      5,           5.,           0.,           5.
      6,           5.,           5.,           5.
      7,           5.,          -5.,          10.
      8,           5.,           0.,          10.
      9,           5.,           5.,          10.
     10,           0.,          -5.,           0.
     11,           0.,           0.,           0.
     12,           0.,           5.,           0.
     13,           0.,          -5.,           5.
     14,           0.,           0.,           5.
     15,           0.,           5.,           5.
     16,           0.,          -5.,          10.
     17,           0.,           0.,          10.
     18,           0.,           5.,          10.
     19,          -5.,          -5.,           0.
     20,          -5.,           0.,           0.
     21,          -5.,           5.,           0.
     22,          -5.,          -5.,           5.
     23,          -5.,           0.,           5.
     24,          -5.,           5.,           5.
     25,          -5.,          -5.,          10.
     26,          -5.,           0.,          10.
     27,          -5.,           5.,          10.
*NSET,NSET=OUT, GENERATE
1,  7,1
10,15,1
19,26,1
*Element, type=C3D8, Elset=E12
1, 10, 11, 14, 13,  1,  2,  5,  4
2, 11, 12, 15, 14,  2,  3,  6,  5
3, 13, 14, 17, 16,  4,  5,  8,  7
4, 14, 15, 18, 17,  5,  6,  9,  8
*ELEMENT, TYPE=C3D8, ELSET=E34
5, 19, 20, 23, 22, 10, 11, 14, 13
6, 20, 21, 24, 23, 11, 12, 15, 14
7, 22, 23, 26, 25, 13, 14, 17, 16
8, 23, 24, 27, 26, 14, 15, 18, 17
*NSET, NSET=TOPOUT
3, 9, 21, 27
*NSET, NSET=TOPIN
6, 12, 18, 24
*NSET, NSET=TOPCEN
15, 
*NSET,NSET=BOTTOM
1, 4, 7, 10, 13, 16, 
19, 22, 25
*SOLID SECTION, ELSET=E12, MATERIAL=HARD
*SOLID SECTION, ELSET=E34, MATERIAL=SOFT
**
*MATERIAL,NAME=SOFT
*ELASTIC
<E1>,<N1>
*MATERIAL,NAME=HARD
*ELASTIC
<E2>,<N2>
*PARAMETER SHAPE VARIATION, PARAMETER=S2
5,  0.1, 0.0, 0.0
14, 0.1, 0.0, 0.0
23, 0.1, 0.0, 0.0
**
*BOUNDARY 
BOTTOM,1,3,0.0
*STEP,DSA
   SMALL DISPLACEMENT ANALYSIS
*STATIC,DIRECT
1., 5.
**
*CLOAD
TOPOUT, 2, 1.E4
TOPIN,  2, 2.E4
TOPCEN, 2, 4.E4
**
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT,NSET=OUT
U, 
*ELEMENT OUTPUT
S, EVOL, ELEN
*DESIGN RESPONSE,FREQ=10
*NODE RESPONSE,NSET=OUT
U
*ELEMENT RESPONSE
S, EVOL, ELEN
*END STEP