*HEADING
BOUSSINESQ PROBLEM, 12 CAX8R + 4 CINAX5R & CAX8 only.
*preprint,model=yes,hist=yes
**************************************************
*part,name=analytic
*node,nset=all
1, 0, -0.5
2, 0, -1.0
3, 0, -2.0
4, 0, -3.0
5, 0, -4.0
*element,type=b21,elset=dummy
1, 1,2
2, 2,3
3, 3,4
4, 4,5
*beamsection,material=one,elset=dummy,sec=circ
1.0e-1,
*endpart
************************************************
*part,name=cax8cinax5
*NODE
1,0.,0.
3,.5,0.
2,0.,-.5
5,.375,-.875
7,.875,-.375
9,0.,-1.5
12,1.5,0.
10,.75,-1.25
11,1.25,-.75
4,0.,-1.
8,1.,0.
6,.75,-.75
21,0.,-2.
25,0.,-4.
101,2.,0.
105,4.,0.
61,1.5,-1.5
65,2.8333333,-2.8333333
27,0.,-8.
67,5.666666,-5.666666
107,8.,0.
*NGEN
21,25
61,65
101,105
21,61,10
22,62,20
23,63,10
24,64,20
25,65,10
27,67,20
61,101,10
62,102,20
63,103,10
64,104,20
65,105,10
67,107,20
*NSET,NSET=LHS,GENERATE
1,2,1
4,9,5
21,25,1
27,27,1
*ELEMENT,TYPE=CAX8R,ELSET=ALL
1,4,6,8,1,5,7,3,2
2,21,41,6,4,31,10,5,9
3,41,61,81,6,51,71,11,10
4,81,101,8,6,91,12,7,11
5,23,43,41,21,33,42,31,22
*ELGEN,ELSET=ALL
5,4,20,1,2,2,4
*ELEMENT,TYPE=CINAX5R,ELSET=ALL
13,45,25,27,47,35
14,65,45,47,67,55
15,85,65,67,87,75
16,105,85,87,107,95
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
*****************************************
*part,name=cax8
*NODE
1,0.,0.
3,.5,0.
2,0.,-.5
5,.375,-.875
7,.875,-.375
9,0.,-1.5
12,1.5,0.
10,.75,-1.25
11,1.25,-.75
4,0.,-1.
8,1.,0.
6,.75,-.75
21,0.,-2.
**26,0.,-6.
27,0.,-8.
61,1.5,-1.5
**65,2.8333333,-2.8333333
**66,4.2426407,-4.2426407
67,5.666666,-5.666666
101,2.,0.
**105,4.,0.
**106,6.,0.
107,8.,0.
*NGEN
21,27
61,67
101,107
21,61,10
22,62,10
23,63,10
24,64,10
25,65,10
26,66,10
27,67,10
61,101,10
62,102,10
63,103,10
64,104,10
65,105,10
66,106,10
67,107,10
*NSET,NSET=LHS,GENERATE
1,2,1
4,9,5
21,25,1
27,27,1
*NSET,NSET=FULLYFIXED,GEN
27,107,10
*ELEMENT,TYPE=CAX8R,ELSET=ALL
1,4,6,8,1,5,7,3,2
2,21,41,6,4,31,10,5,9
3,41,61,81,6,51,71,11,10
4,81,101,8,6,91,12,7,11
5,23,43,41,21,33,42,31,22
*ELGEN,ELSET=ALL
5,4,20,1,3,2,4
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
*****************************************
*assembly,name=all
*instance,name=analytic,part=analytic
0,0,0,
*endinstance
*instance,name=cax8cinax5,part=cax8cinax5
0,0,0
*endinstance
*instance,name=cax8,part=cax8
0,10,0
*endinstance
*endassembly
*****************************************
*MATERIAL,NAME=ONE
*ELASTIC
1.,.1
**********************************  
*BOUNDARY
cax8cinax5.LHS, 1,1,0.0
cax8.LHS,       1,1,0.0
cax8.fullyfixed,1,2,0.0
analytic.all,1,
*****************************************
*STEP
*STATIC
*RESTART,WRITE
*CLOAD
cax8cinax5.1,2,-1.
cax8.1,      2,-1.
*BOUNDARY
analytic.1,2,2,-.9803
analytic.2,2,2,-.4902
analytic.3,2,2,-.2451
analytic.4,2,2,-.1634
analytic.5,2,2,-.1225
*NODE FILE, NSET=cax8cinax5.LHS
U, 
*NODE PRINT,NSET=cax8cinax5.LHS
COORD,U
*NODE FILE,NSET=cax8.LHS
U,
*NODE PRINT,NSET=cax8.LHS
COORD,U
*NODE PRINT,NSET=analytic.all
COORD,U
*OUTPUT,FIELD
*NODE OUTPUT
U,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=cax8cinax5.LHS
U,
*NODE OUTPUT,NSET=cax8.LHS
U,
*NODE OUTPUT,NSET=analytic.all
U,
*END STEP