*HEADING
 o TESTING TOTAL DSA FORMULATION FOR CONTINUUM 3D ELEMENTS (C3D_)
**
** Node/element numbering scheme:
**    element test  |   node labels   |   element labels
**        c3d4             10_               101-106
**        c3d4h            20_               201-206
**        c3d6             30_               301-302
**        c3d6h            40_               401-402
**        c3d8             50_               501-502
**        c3d8h            60_               601-602
**        c3d8i            70_               701-702
**        c3d8ih           80_               801-802
**        c3d8r            90_               901-902
**        c3d8rh          100_              1001-1002
**        c3d10           110_              1101-1106
**        c3d10h          120_              1201-1206
**        c3d10m          130_              1301-1306
**        c3d10mh         140_              1401-1406
**        c3d15           150_              1501-1502
**        c3d15h          160_              1601-1602
**        c3d20           170_              1701-1702
**        c3d20h          180_              1801-1802
**        c3d20r          190_              1901-1902
**        c3d20rh         200_              2001-2002
**
**        c3d15v          210_              2101-2102
**        c3d15vh         220_              2201-2202
**        c3d27           230_              2301-2302
**        c3d27h          240_              2401-2402
**        c3d27r          250_              2501-2502
**        c3d27rh         260_              2601-2602
**
**         | 3
**         |
**         |      
**               O-----------O----------O
**              /                      /|
**             /                      / |
**            O     FACE2 O          O  |
**           /                      /   O
**          /                      /    |
**         O-----------O----------O   FACE4
**         |                      |  O  |
**   FACE6 |                      |     O
**         |    FACE3             |    /
**         O           O          O   /  
**         | /2                   |  O 
**         |/                     | /
**         /                      |/
**         O-----------O----------O  ---- 1
**                  FACE1 
**
*PARAMETER
dload1 =  -.2
dload2 =  -.1
dload3 =  -.05
shape = 0.0
young = 100000.
nu    = 0.3
*DESIGN PARAMETER
nu,shape
*DSA CONTROLS, FORMULATION=TOTAL
nu,CD,1.e-03
shape,CD,1e-07
**
*NODE,NSET=C3D41
10101, 0.0, 1.0, 0.0
10102, 0.0, 0.0, 0.0
10103, 1.0, 0.0, 0.0
10104, 1.0, 1.0, 0.0
*NSET,NSET=C3D4
C3D41,
*NCOPY,OLD SET=C3D41,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D4
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D4H1
20101, 0.0, 1.0, 0.0
20102, 0.0, 0.0, 0.0
20103, 1.0, 0.0, 0.0
20104, 1.0, 1.0, 0.0
*NSET,NSET=C3D4H
C3D4H1,
*NCOPY,OLD SET=C3D4H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D4H
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D61
30101, 0.0, 1.0, 0.0
30102, 0.0, 0.0, 0.0
30103, 1.0, 0.0, 0.0
30104, 1.0, 1.0, 0.0
*NSET,NSET=C3D6
C3D61,
*NCOPY,OLD SET=C3D61,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D6
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D6H1
40101, 0.0, 1.0, 0.0
40102, 0.0, 0.0, 0.0
40103, 1.0, 0.0, 0.0
40104, 1.0, 1.0, 0.0
*NSET,NSET=C3D6H
C3D6H1,
*NCOPY,OLD SET=C3D6H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D6H
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D81
50101, 0.0, 1.0, 0.0
50102, 0.0, 0.0, 0.0
50103, 1.0, 0.0, 0.0
50104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8
C3D81,
*NCOPY,OLD SET=C3D81,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D81,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D8H1
60101, 0.0, 1.0, 0.0
60102, 0.0, 0.0, 0.0
60103, 1.0, 0.0, 0.0
60104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8H
C3D8H1,
*NCOPY,OLD SET=C3D8H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8H
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D8H1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8H
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D8I1
70101, 0.0, 1.0, 0.0
70102, 0.0, 0.0, 0.0
70103, 1.0, 0.0, 0.0
70104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8I
C3D8I1,
*NCOPY,OLD SET=C3D8I1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8I
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D8I1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8I
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D8IH1
80101, 0.0, 1.0, 0.0
80102, 0.0, 0.0, 0.0
80103, 1.0, 0.0, 0.0
80104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8IH
C3D8IH1,
*NCOPY,OLD SET=C3D8IH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8IH
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D8IH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8IH
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D8R1
90101, 0.0, 1.0, 0.0
90102, 0.0, 0.0, 0.0
90103, 1.0, 0.0, 0.0
90104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8R
C3D8R1,
*NCOPY,OLD SET=C3D8R1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8R
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D8R1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8R
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D8RH1
100101, 0.0, 1.0, 0.0
100102, 0.0, 0.0, 0.0
100103, 1.0, 0.0, 0.0
100104, 1.0, 1.0, 0.0
*NSET,NSET=C3D8RH
C3D8RH1,
*NCOPY,OLD SET=C3D8RH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D8RH
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D8RH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D8RH
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D101
110101, 0.0, 1.0, 0.0
110102, 0.0, 0.0, 0.0
110103, 1.0, 0.0, 0.0
110104, 1.0, 1.0, 0.0
110105, 0.0, 0.5, 0.0
110106, 0.5, 0.0, 0.0
110107, 1.0, 0.5, 0.0
110108, 0.5, 1.0, 0.0
110109, 0.5, 0.5, 0.0
*NSET,NSET=C3D10
C3D101,
*NCOPY,OLD SET=C3D101,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D10
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D101,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D10
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D10H1
120101, 0.0, 1.0, 0.0
120102, 0.0, 0.0, 0.0
120103, 1.0, 0.0, 0.0
120104, 1.0, 1.0, 0.0
120105, 0.0, 0.5, 0.0
120106, 0.5, 0.0, 0.0
120107, 1.0, 0.5, 0.0
120108, 0.5, 1.0, 0.0
120109, 0.5, 0.5, 0.0
*NSET,NSET=C3D10H
C3D10H1,
*NCOPY,OLD SET=C3D10H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D10H
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D10H1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D10H
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D10M1
130101, 0.0, 1.0, 0.0
130102, 0.0, 0.0, 0.0
130103, 1.0, 0.0, 0.0
130104, 1.0, 1.0, 0.0
130105, 0.0, 0.5, 0.0
130106, 0.5, 0.0, 0.0
130107, 1.0, 0.5, 0.0
130108, 0.5, 1.0, 0.0
130109, 0.5, 0.5, 0.0
*NSET,NSET=C3D10M
C3D10M1,
*NCOPY,OLD SET=C3D10M1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D10M
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D10M1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D10M
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D10MH1
140101, 0.0, 1.0, 0.0
140102, 0.0, 0.0, 0.0
140103, 1.0, 0.0, 0.0
140104, 1.0, 1.0, 0.0
140105, 0.0, 0.5, 0.0
140106, 0.5, 0.0, 0.0
140107, 1.0, 0.5, 0.0
140108, 0.5, 1.0, 0.0
140109, 0.5, 0.5, 0.0
*NSET,NSET=C3D10MH
C3D10MH1,
*NCOPY,OLD SET=C3D10MH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D10MH
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D10MH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D10MH
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D151
150101, 0.0, 1.0, 0.0
150102, 0.0, 0.0, 0.0
150103, 1.0, 0.0, 0.0
150104, 1.0, 1.0, 0.0
150105, 0.0, 0.5, 0.0
150106, 0.5, 0.0, 0.0
150107, 1.0, 0.5, 0.0
150108, 0.5, 1.0, 0.0
150109, 0.5, 0.5, 0.0
*NSET,NSET=C3D15
C3D151,
*NCOPY,OLD SET=C3D151,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D15
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D151,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D15
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D15H1
160101, 0.0, 1.0, 0.0
160102, 0.0, 0.0, 0.0
160103, 1.0, 0.0, 0.0
160104, 1.0, 1.0, 0.0
160105, 0.0, 0.5, 0.0
160106, 0.5, 0.0, 0.0
160107, 1.0, 0.5, 0.0
160108, 0.5, 1.0, 0.0
160109, 0.5, 0.5, 0.0
*NSET,NSET=C3D15H
C3D15H1,
*NCOPY,OLD SET=C3D15H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D15H
0.0, 0.0, 0.5
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D15H1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D15H
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D201
170101, 0.0, 1.0, 0.0
170102, 0.0, 0.0, 0.0
170103, 1.0, 0.0, 0.0
170104, 1.0, 1.0, 0.0
170105, 0.0, 0.5, 0.0
170106, 0.5, 0.0, 0.0
170107, 1.0, 0.5, 0.0
170108, 0.5, 1.0, 0.0
*NSET,NSET=C3D20
C3D201,
*NCOPY,OLD SET=C3D201,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D20
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D201,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D20
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D201,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D20
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D201,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D20
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D20H1
180101, 0.0, 1.0, 0.0
180102, 0.0, 0.0, 0.0
180103, 1.0, 0.0, 0.0
180104, 1.0, 1.0, 0.0
180105, 0.0, 0.5, 0.0
180106, 0.5, 0.0, 0.0
180107, 1.0, 0.5, 0.0
180108, 0.5, 1.0, 0.0
*NSET,NSET=C3D20H
C3D20H1,
*NCOPY,OLD SET=C3D20H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D20H
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20H1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D20H
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20H1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D20H
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20H1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D20H
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D20R1
190101, 0.00, 1.00, 0.00
190102, 0.00, 0.00, 0.00
190103, 1.00, 0.00, 0.00
190104, 1.00, 1.00, 0.00
190105, 0.00, 0.50, 0.00
190106, 0.50, 0.00, 0.00
190107, 1.00, 0.50, 0.00
190108, 0.50, 1.00, 0.00
*NSET,NSET=C3D20R
C3D20R1,
*NCOPY,OLD SET=C3D20R1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D20R
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20R1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D20R
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20R1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D20R
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20R1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D20R
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D20RH1
200101, 0.00, 1.00, 0.00
200102, 0.00, 0.00, 0.00
200103, 1.00, 0.00, 0.00
200104, 1.00, 1.00, 0.00
200105, 0.00, 0.50, 0.00
200106, 0.50, 0.00, 0.00
200107, 1.00, 0.50, 0.00
200108, 0.50, 1.00, 0.00
*NSET,NSET=C3D20RH
C3D20RH1,
*NCOPY,OLD SET=C3D20RH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D20RH
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20RH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D20RH
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20RH1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D20RH
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D20RH1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D20RH
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D15V1
210101, 0.00, 1.00, 0.00
210102, 0.00, 0.00, 0.00
210103, 1.00, 0.00, 0.00
210104, 1.00, 1.00, 0.00
210105, 0.00, 0.50, 0.00
210106, 0.50, 0.00, 0.00
210107, 1.00, 0.50, 0.00
210108, 0.50, 1.00, 0.00
210109, 0.50, 0.50, 0.00
*NSET,NSET=C3D15V
C3D15V1,
*NCOPY,OLD SET=C3D15V1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D15V
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D15V1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D15V
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D15VH1
220101, 0.00, 1.00, 0.00
220102, 0.00, 0.00, 0.00
220103, 1.00, 0.00, 0.00
220104, 1.00, 1.00, 0.00
220105, 0.00, 0.50, 0.00
220106, 0.50, 0.00, 0.00
220107, 1.00, 0.50, 0.00
220108, 0.50, 1.00, 0.00
220109, 0.50, 0.50, 0.00
*NSET,NSET=C3D15VH
C3D15VH1,
*NCOPY,OLD SET=C3D15VH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D15VH
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D15VH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D15VH
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D271
230101, 0.00, 1.00, 0.00
230102, 0.00, 0.00, 0.00
230103, 1.00, 0.00, 0.00
230104, 1.00, 1.00, 0.00
230105, 0.00, 0.50, 0.00
230106, 0.50, 0.00, 0.00
230107, 1.00, 0.50, 0.00
230108, 0.50, 1.00, 0.00
230109, 0.50, 0.50, 0.00
*NSET,NSET=C3D27
C3D271,
*NCOPY,OLD SET=C3D271,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D27
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D271,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D27
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D271,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D27
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D271,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D27
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D27H1
240101, 0.00, 1.00, 0.00
240102, 0.00, 0.00, 0.00
240103, 1.00, 0.00, 0.00
240104, 1.00, 1.00, 0.00
240105, 0.00, 0.50, 0.00
240106, 0.50, 0.00, 0.00
240107, 1.00, 0.50, 0.00
240108, 0.50, 1.00, 0.00
240109, 0.50, 0.50, 0.00
*NSET,NSET=C3D27H
C3D27H1,
*NCOPY,OLD SET=C3D27H1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D27H
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27H1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D27H
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27H1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D27H
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27H1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D27H
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D27R1
250101, 0.00, 1.00, 0.00
250102, 0.00, 0.00, 0.00
250103, 1.00, 0.00, 0.00
250104, 1.00, 1.00, 0.00
250105, 0.00, 0.50, 0.00
250106, 0.50, 0.00, 0.00
250107, 1.00, 0.50, 0.00
250108, 0.50, 1.00, 0.00
250109, 0.50, 0.50, 0.00
*NSET,NSET=C3D27R
C3D27R1,
*NCOPY,OLD SET=C3D27R1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D27R
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27R1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D27R
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27R1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D27R
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27R1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D27R
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NODE,NSET=C3D27RH1
260101, 0.00, 1.00, 0.00
260102, 0.00, 0.00, 0.00
260103, 1.00, 0.00, 0.00
260104, 1.00, 1.00, 0.00
260105, 0.00, 0.50, 0.00
260106, 0.50, 0.00, 0.00
260107, 1.00, 0.50, 0.00
260108, 0.50, 1.00, 0.00
260109, 0.50, 0.50, 0.00
*NSET,NSET=C3D27RH
C3D27RH1,
*NCOPY,OLD SET=C3D27RH1,CHANGE NUMBER=100,SHIFT,MULTIPLE=1,NEW SET=C3D27RH
0.0, 0.0, 0.25
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27RH1,CHANGE NUMBER=200,SHIFT,MULTIPLE=1,NEW SET=C3D27RH
0.0, 0.0, 0.50
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27RH1,CHANGE NUMBER=300,SHIFT,MULTIPLE=1,NEW SET=C3D27RH
0.0, 0.0, 0.75
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*NCOPY,OLD SET=C3D27RH1,CHANGE NUMBER=400,SHIFT,MULTIPLE=1,NEW SET=C3D27RH
0.0, 0.0, 1.00
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
********
*NSET,NSET=PAR1BOT,GENERATE
10101,260101,10000
*NSET,NSET=PAR2BOT,GENERATE
10102,260102,10000
*NSET,NSET=PAR3BOT,GENERATE
10103,260103,10000
*NSET,NSET=PAR4BOT,GENERATE
10104,260104,10000
*NSET,NSET=PAR5BOT,GENERATE
110105,260105,10000
*NSET,NSET=PAR6BOT,GENERATE
110106,260106,10000
*NSET,NSET=PAR7BOT,GENERATE
110107,260107,10000
*NSET,NSET=PAR8BOT,GENERATE
110108,260108,10000
*NSET,NSET=PAR9BOT,GENERATE
110109,160109,10000
210109,260109,10000
*NSET,NSET=PAR1MID,GENERATE
50201,160201,10000
170301,200301,10000
210101,220101,10000
230101,260101,10000
*NSET,NSET=PAR2MID,GENERATE
50202,160202,10000
170302,200302,10000
210202,220202,10000
230302,260302,10000
*NSET,NSET=PAR3MID,GENERATE
50203,160203,10000
170303,200303,10000
210203,220203,10000
230303,260303,10000
*NSET,NSET=PAR4MID,GENERATE
50204,160204,10000
170304,200304,10000
210204,220204,10000
230304,260304,10000
*NSET,NSET=PAR5MID,GENERATE
110205,160205,10000
170305,200305,10000
210205,220205,10000
230305,260305,10000
*NSET,NSET=PAR6MID,GENERATE
110206,160206,10000
170306,200306,10000
210206,220206,10000
230306,260306,10000
*NSET,NSET=PAR7MID,GENERATE
110207,160207,10000
170307,200307,10000
210207,220207,10000
230307,260307,10000
*NSET,NSET=PAR8MID,GENERATE
110208,160208,10000
170308,200308,10000
210208,220208,10000
230308,260308,10000
*NSET,NSET=PAR9MID,GENERATE
110209,160209,10000
210209,220209,10000
230309,260309,10000
*NSET,NSET=PAR1TOP,GENERATE
10201,40201,10000
50301,160301,10000
170501,200501,10000
210301,220301,10000
230501,260501,10000
*NSET,NSET=PAR2TOP,GENERATE
10202,40202,10000
50302,160302,10000
170502,200502,10000
210302,220302,10000
230502,260502,10000
*NSET,NSET=PAR3TOP,GENERATE
10203,40203,10000
50303,160303,10000
170503,200503,10000
210303,220303,10000
230503,260503,10000
*NSET,NSET=PAR4TOP,GENERATE
10204,40204,10000
50304,160304,10000
170504,200504,10000
210304,220304,10000
230504,260504,10000
*NSET,NSET=PAR5TOP,GENERATE
110305,160305,10000
170505,200505,10000
210305,220305,10000
230505,260505,10000
*NSET,NSET=PAR6TOP,GENERATE
110306,160306,10000
170506,200506,10000
210306,220306,10000
230506,260506,10000
*NSET,NSET=PAR7TOP,GENERATE
110307,160307,10000
170507,200507,10000
210307,220307,10000
230507,260507,10000
*NSET,NSET=PAR8TOP,GENERATE
110308,160308,10000
170508,200508,10000
210308,220308,10000
230508,260508,10000
*NSET,NSET=PAR9TOP,GENERATE
110309,160309,10000
210309,220309,10000
230509,260509,10000
**
*PARAMETER SHAPE VARIATION, PARAMETER=shape
**
PAR1BOT,   0.5,  0.0,  0.250
PAR2BOT,   0.0,  0.0,  0.000
PAR3BOT,   0.0, -0.5,  0.250
PAR4BOT,   1.0,  1.0,  0.500
PAR5BOT,   0.25, 0.0,  0.125
PAR6BOT,   0.0,  0.25, 0.125
PAR7BOT,   0.75, 0.75, 0.375
PAR8BOT,   0.75, 0.75, 0.375
PAR9BOT,   0.50, 0.5,  0.250
**
PAR1MID,  -0.5,  0.0, -0.250
PAR2MID,  -0.0,  0.0, -0.000
PAR3MID,   0.0, -0.5, -0.250
PAR4MID,  -1.0, -1.0, -0.500
PAR5MID,  -0.25, 0.0, -0.125
PAR6MID,   0.0, -0.25,-0.125
PAR7MID,  -0.75,-0.75,-0.375
PAR8MID,  -0.75,-0.75,-0.375
PAR9MID,  -0.50,-0.5, -0.250
**
PAR1TOP,   1.0,  0.0,  0.50
PAR2TOP,   0.0,  0.0,  0.00
PAR3TOP,   0.0,  1.0,  0.50
PAR4TOP,   2.0,  2.0,  1.00
PAR5TOP,   0.5,  0.0,  0.25
PAR6TOP,   0.0,  0.5,  0.25
PAR7TOP,   1.5,  1.5,  0.75
PAR8TOP,   1.5,  1.5,  0.75
PAR9TOP,   1.0,  1.0,  0.50
**
*NSET,NSET=AMID,GENERATE
50201,160201,10000
170301,200301,10000
210201,220201,10000
230301,260301,10000
*NSET,NSET=BMID,GENERATE
50202,160202,10000
170302,200302,10000
210202,220202,10000
230302,260302,10000
*NSET,NSET=CMID,GENERATE
50203,160203,10000
170303,200303,10000
210203,220203,10000
230303,260303,10000
*NSET,NSET=DMID,GENERATE
50204,160204,10000
170304,200304,10000
210204,220204,10000
230304,260304,10000
*NSET,NSET=ATOP,GENERATE
10201,40201,10000
50301,160301,10000
170501,200501,10000
210301,220301,10000
230501,260501,10000
*NSET,NSET=BTOP,GENERATE
10202,40202,10000
50302,160302,10000
170502,200502,10000
210302,220302,10000
230502,260502,10000
*NSET,NSET=CTOP,GENERATE
10203,40203,10000
50303,160303,10000
170503,200503,10000
210303,220303,10000
230503,260503,10000
*NSET,NSET=DTOP,GENERATE
10204,40204,10000
50304,160304,10000
170504,200504,10000
210304,220304,10000
230504,260504,10000
*NSET,NSET=FACE1,GENERATE
10101,10104,1
20101,20104,1
30101,30104,1
40101,40104,1
50101,50104,1
60101,60104,1
70101,70104,1
80101,80104,1
90101,90104,1
100101,100104,1
110101,110109,1
120101,120109,1
130101,130109,1
140101,140109,1
150101,150109,1
160101,160109,1
170101,170108,1
180101,180108,1
190101,190108,1
200101,200108,1
210101,210108,1
220101,220108,1
230101,230109,1
240101,240109,1
250101,250109,1
260101,260109,1
**
*NSET,NSET=NOUT
PAR4TOP,
**
*ELEMENT,ELSET=C3D4,TYPE=C3D4
101,10201,10102,10203,10202
102,10201,10203,10102,10103
103,10201,10103,10102,10101
104,10201,10203,10104,10204
105,10201,10104,10203,10103
106,10201,10104,10103,10101
*ELEMENT,ELSET=C3D4H,TYPE=C3D4H
201,20201,20102,20203,20202
202,20201,20203,20102,20103
203,20201,20103,20102,20101
204,20201,20203,20104,20204
205,20201,20104,20203,20103
206,20201,20104,20103,20101
*ELEMENT,ELSET=C3D6,TYPE=C3D6
301,30104,30101,30102,30204,30201,30202
302,30104,30102,30103,30204,30202,30203
*ELEMENT,ELSET=C3D6H,TYPE=C3D6H
401,40104,40101,40102,40204,40201,40202
402,40104,40102,40103,40204,40202,40203
*ELEMENT,ELSET=C3D8,TYPE=C3D8
501,50101,50102,50103,50104,50201,50202,50203,50204
*ELGEN,ELSET=C3D8
501,2,100,1
*ELEMENT,ELSET=C3D8H,TYPE=C3D8H
601,60101,60102,60103,60104,60201,60202,60203,60204
*ELGEN,ELSET=C3D8H
601,2,100,1
*ELEMENT,ELSET=C3D8I,TYPE=C3D8I
701,70101,70102,70103,70104,70201,70202,70203,70204
*ELGEN,ELSET=C3D8I
701,2,100,1
*ELEMENT,ELSET=C3D8IH,TYPE=C3D8IH
801,80101,80102,80103,80104,80201,80202,80203,80204
*ELGEN,ELSET=C3D8IH
801,2,100,1
*ELEMENT,ELSET=C3D8R,TYPE=C3D8R
901,90101,90102,90103,90104,90201,90202,90203,90204
*ELGEN,ELSET=C3D8R
901,2,100,1
*ELEMENT,ELSET=C3D8RH,TYPE=C3D8RH
1001,100101,100102,100103,100104,100201,100202,100203,100204
*ELGEN,ELSET=C3D8RH
1001,2,100,1
*ELEMENT,ELSET=C3D10,TYPE=C3D10
1101,110301,110102,110303,110302,110205,110206,110309,110305,110202,110306
1102,110301,110303,110102,110103,110309,110206,110205,110209,110203,110106
1103,110301,110103,110102,110101,110209,110106,110205,110201,110109,110105
1104,110301,110303,110104,110304,110309,110207,110208,110308,110307,110204
1105,110301,110104,110303,110103,110208,110207,110309,110209,110107,110203
1106,110301,110104,110103,110101,110208,110107,110209,110201,110108,110109
*ELEMENT,ELSET=C3D10H,TYPE=C3D10H
1201,120301,120102,120303,120302,120205,120206,120309,120305,120202,120306
1202,120301,120303,120102,120103,120309,120206,120205,120209,120203,120106
1203,120301,120103,120102,120101,120209,120106,120205,120201,120109,120105
1204,120301,120303,120104,120304,120309,120207,120208,120308,120307,120204
1205,120301,120104,120303,120103,120208,120207,120309,120209,120107,120203
1206,120301,120104,120103,120101,120208,120107,120209,120201,120108,120109
*ELEMENT,ELSET=C3D10M,TYPE=C3D10M
1301,130301,130102,130303,130302,130205,130206,130309,130305,130202,130306
1302,130301,130303,130102,130103,130309,130206,130205,130209,130203,130106
1303,130301,130103,130102,130101,130209,130106,130205,130201,130109,130105
1304,130301,130303,130104,130304,130309,130207,130208,130308,130307,130204
1305,130301,130104,130303,130103,130208,130207,130309,130209,130107,130203
1306,130301,130104,130103,130101,130208,130107,130209,130201,130108,130109
*ELEMENT,ELSET=C3D10MH,TYPE=C3D10MH
1401,140301,140102,140303,140302,140205,140206,140309,140305,140202,140306
1402,140301,140303,140102,140103,140309,140206,140205,140209,140203,140106
1403,140301,140103,140102,140101,140209,140106,140205,140201,140109,140105
1404,140301,140303,140104,140304,140309,140207,140208,140308,140307,140204
1405,140301,140104,140303,140103,140208,140207,140309,140209,140107,140203
1406,140301,140104,140103,140101,140208,140107,140209,140201,140108,140109
*ELEMENT,ELSET=C3D15,TYPE=C3D15
1501,150104,150101,150102,150304,150301,150302,150108,150105,150109,
150308,150305,150309,150204,150201,150202
1502,150104,150102,150103,150304,150302,150303,150109,150106,150107,
150309,150306,150307,150204,150202,150203
*ELEMENT,ELSET=C3D15H,TYPE=C3D15H
1601,160104,160101,160102,160304,160301,160302,160108,160105,160109,
160308,160305,160309,160204,160201,160202
1602,160104,160102,160103,160304,160302,160303,160109,160106,160107,
160309,160306,160307,160204,160202,160203
*ELEMENT,ELSET=C3D20,TYPE=C3D20
1701,170101,170102,170103,170104,170301,170302,170303,170304,
170105,170106,170107,170108,170305,170306,170307,170308,
170201,170202,170203,170204
1702,170301,170302,170303,170304,170501,170502,170503,170504,
170305,170306,170307,170308,170505,170506,170507,170508,
170401,170402,170403,170404
*ELEMENT,ELSET=C3D20H,TYPE=C3D20H
1801,180101,180102,180103,180104,180301,180302,180303,180304,
180105,180106,180107,180108,180305,180306,180307,180308,
180201,180202,180203,180204
1802,180301,180302,180303,180304,180501,180502,180503,180504,
180305,180306,180307,180308,180505,180506,180507,180508,
180401,180402,180403,180404
*ELEMENT,ELSET=C3D20R,TYPE=C3D20R
1901,190101,190102,190103,190104,190301,190302,190303,190304,
190105,190106,190107,190108,190305,190306,190307,190308,
190201,190202,190203,190204
1902,190301,190302,190303,190304,190501,190502,190503,190504,
190305,190306,190307,190308,190505,190506,190507,190508,
190401,190402,190403,190404
*ELEMENT,ELSET=C3D20RH,TYPE=C3D20RH
2001,200101,200102,200103,200104,200301,200302,200303,200304,
200105,200106,200107,200108,200305,200306,200307,200308,
200201,200202,200203,200204
2002,200301,200302,200303,200304,200501,200502,200503,200504,
200305,200306,200307,200308,200505,200506,200507,200508,
200401,200402,200403,200404
*ELEMENT,ELSET=C3D15V,TYPE=C3D15V
2101,210104,210101,210102,210304,210301,210302,210108,210105,210109,
210308,210305,210309,210204,210201,210202,210208,210205,210209
2102,210104,210102,210103,210304,210302,210303,210109,210106,210107,
210309,210306,210307,210204,210202,210203,210209,210206,210207
*ELEMENT,ELSET=C3D15VH,TYPE=C3D15VH
2201,220104,220101,220102,220304,220301,220302,220108,220105,220109,
220308,220305,220309,220204,220201,220202,220208,220205,220209
2202,220104,220102,220103,220304,220302,220303,220109,220106,220107,
220309,220306,220307,220204,220202,220203,220209,220206,220207
*ELEMENT,ELSET=C3D27,TYPE=C3D27
2301,230101,230102,230103,230104,230301,230302,230303,230304,
230105,230106,230107,230108,230305,230306,230307,230308,
230201,230202,230203,230204,230209,230109,230309,
230205,230206,230207,230208
2302,230301,230302,230303,230304,230501,230502,230503,230504,
230305,230306,230307,230308,230505,230506,230507,230508,
230401,230402,230403,230404,230409,230309,230509,
230405,230406,230407,230408
*ELEMENT,ELSET=C3D27H,TYPE=C3D27H
2401,240101,240102,240103,240104,240301,240302,240303,240304,
240105,240106,240107,240108,240305,240306,240307,240308,
240201,240202,240203,240204,240209,240109,240309,
240205,240206,240207,240208
2402,240301,240302,240303,240304,240501,240502,240503,240504,
240305,240306,240307,240308,240505,240506,240507,240508,
240401,240402,240403,240404,240409,240309,240509,
240405,240406,240407,240408
*ELEMENT,ELSET=C3D27R,TYPE=C3D27R
2501,250101,250102,250103,250104,250301,250302,250303,250304,
250105,250106,250107,250108,250305,250306,250307,250308,
250201,250202,250203,250204,250209,250109,250309,
250205,250206,250207,250208
2502,250301,250302,250303,250304,250501,250502,250503,250504,
250305,250306,250307,250308,250505,250506,250507,250508,
250401,250402,250403,250404,250409,250309,250509,
250405,250406,250407,250408
*ELEMENT,ELSET=C3D27RH,TYPE=C3D27RH
2601,260101,260102,260103,260104,260301,260302,260303,260304,
260105,260106,260107,260108,260305,260306,260307,260308,
260201,260202,260203,260204,260209,260109,260309,
260205,260206,260207,260208
2602,260301,260302,260303,260304,260501,260502,260503,260504,
260305,260306,260307,260308,260505,260506,260507,260508,
260401,260402,260403,260404,260409,260309,260509,
260405,260406,260407,260408
**
*ELSET,ELSET=ALL
C3D4,C3D4H,C3D6,C3D6H,C3D8,C3D8H,C3D8I,C3D8IH,C3D8R,C3D8RH,
C3D10,C3D10H,C3D10M,C3D10MH,C3D15,C3D15H,C3D20,C3D20H,C3D20R,C3D20RH,
C3D15V,C3D15VH,C3D27,C3D27H,C3D27R,C3D27RH
**
** DLOADS FOR MOST ELEMENTS:
*ELSET,ELSET=DLOAD2
104,204,301,302,401,402,502,602,702,802,902,1002,
1104,1204,1304,1404,1501,1502,1601,1602,1702,
1802,1902,2002,2102,2202,2302,2402,2502,2602
**
*ELSET,ELSET=DLOAD4
101,201,1101,1201,1301,1401
**
*ELSET,ELSET=EOUT,GENERATE
102,2602,100
**
*SOLID SECTION,ELSET=ALL,MATERIAL=MAT
**
*MATERIAL,NAME=MAT
*ELASTIC
<young>,<nu>
**
*BOUNDARY
FACE1,1,3
**
*STEP,DSA
STEP 1: GEOMETRICALLY LINEAR; TOTAL DSA
*STATIC
0.1, 0.4
**
*DLOAD
DLOAD2,P2, <dload1>
DLOAD4,P4, <dload1>
**
*BOUNDARY
ATOP,1,1,-.001
ATOP,2,2,-.001
BTOP,1,1, .001
BTOP,2,2,-.001
CTOP,1,1, .001
CTOP,2,2, .001
DTOP,1,1,-.001
DTOP,2,2, .001
*BOUNDARY
AMID,1,1,-.0005
AMID,2,2,-.0005
BMID,1,1, .0005
BMID,2,2,-.0005
CMID,1,1, .0005
CMID,2,2, .0005
DMID,1,1,-.0005
DMID,2,2, .0005
**
*OUTPUT,FIELD,FREQ=4
*NODE OUTPUT,NSET=NOUT
U,RF
*ELEMENT OUTPUT,ELSET=EOUT
S,E
*DESIGN RESPONSE,FREQ=4
*NODE RESPONSE,NSET=NOUT
U,RF
*ELEMENT RESPONSE,ELSET=EOUT
S,E
*END STEP
**
*STEP, DSA=NO
STEP 2: GEOMETRICALLY LINEAR; NO DSA
*STATIC
0.1, 0.2
**
*BOUNDARY,OP=NEW
FACE1,1,3
ATOP,1,1, .001
ATOP,2,2, .001
CTOP,1,1,-.001
CTOP,2,2,-.001
AMID,1,1,-.0005
AMID,2,2,-.0005
BMID,1,1, .0005
BMID,2,2,-.0005
CMID,1,1, .0005
CMID,2,2, .0005
DMID,1,1,-.0005
DMID,2,2, .0005
*DLOAD,OP=NEW
DLOAD2,P2, <dload2>
DLOAD4,P4, <dload2>
**
*OUTPUT,FIELD,FREQ=0
*END STEP
**
*STEP,NLGEOM,DSA,UNSYMM=YES
STEP 3: GEOMETRICALLY NONLINEAR; TOTAL DSA
*STATIC
0.1, 0.4
**
*BOUNDARY,FIXED
FACE1,1,3
ATOP,1,1, .001
ATOP,2,2, .001
CTOP,1,1,-.001
CTOP,2,2,-.001
AMID,1,1,-.0005
AMID,2,2,-.0005
BMID,1,1, .0005
BMID,2,2,-.0005
CMID,1,1, .0005
CMID,2,2, .0005
DMID,1,1,-.0005
DMID,2,2, .0005
*DLOAD,OP=NEW
DLOAD2,P2, <dload3>
DLOAD4,P4, <dload3>
**
*OUTPUT,FIELD,FREQ=4
*NODE OUTPUT,NSET=NOUT
U,RF
*ELEMENT OUTPUT,ELSET=EOUT
S,LE
*DESIGN RESPONSE,FREQ=4
*NODE RESPONSE,NSET=NOUT
U,RF
*ELEMENT RESPONSE,ELSET=EOUT
S,LE
*END STEP
*STEP,NLGEOM,DSA,UNSYMM=YES,PERT
STEP 4: PERTURBATION STEP 
*STATIC
*DLOAD
DLOAD2,P2,-0.025
DLOAD4,P2,-0.025
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NOUT
U,RF
*ELEMENT OUTPUT,ELSET=EOUT
S,E
*DESIGN RESPONSE,FREQ=4
*NODE RESPONSE,NSET=NOUT
U,RF
*ELEMENT RESPONSE,ELSET=EOUT
S,E
*END STEP