*HEADING
BOUSSINESQ PROBLEM, 12 CAX4 + 4 CINAX4 & CAX4 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=cax4cinax4
*NODE
1,0.,0.
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.833333333,-2.833333333
*NGEN
21,25,2
61,65,2
101,105,2
21,61,20
23,63,20
25,65,20
61,101,20
63,103,20
65,105,20
*NSET,NSET=INTER,GENERATE
25,105,20
*NCOPY,OLD SET=INTER,CHANGE NUMBER=2,POLE,NEW SET=FAR
1, 
*NSET,NSET=LHS,GENERATE
1,4,3
21,27,2
*ELEMENT,TYPE=CAX4,ELSET=ALL
1,4,6,8,1
2,21,41,6,4
3,41,61,81,6
4,81,101,8,6
5,23,43,41,21
*ELGEN,ELSET=ALL
5,4,20,1,2,2,4
*ELEMENT,TYPE=CINAX4,ELSET=ALL
13,45,25,27,47
14,65,45,47,67
15,85,65,67,87
16,105,85,87,107
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
***************************************************
*part,name=cax4
*NODE
1,0.,0.
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.833333333,-2.833333333
*NGEN
21,25,2
61,65,2
101,105,2
21,61,20
23,63,20
25,65,20
61,101,20
63,103,20
65,105,20
*NSET,NSET=INTER,GENERATE
25,105,20
*NCOPY,OLD SET=INTER,CHANGE NUMBER=2,POLE,NEW SET=FAR
1, 
*NSET,NSET=LHS,GENERATE
1,4,3
21,27,2
*ELEMENT,TYPE=CAX4,ELSET=ALL
1,4,6,8,1
2,21,41,6,4
3,41,61,81,6
4,81,101,8,6
5,23,43,41,21
*ELGEN,ELSET=ALL
5,4,20,1,2,2,4
*ELEMENT,TYPE=CAX4,ELSET=ALL
13,45,25,27,47
14,65,45,47,67
15,85,65,67,87
16,105,85,87,107
*NSET,Nset=fullyfixed
27,47,67,87,107
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
*****************************************
*assembly,name=all
*instance,name=analytic,part=analytic
0,0,0,
*endinstance
*instance,name=cax4cinax4,part=cax4cinax4
0,0,0
*endinstance
*instance,name=cax4,part=cax4
0,10,0
*endinstance
*endassembly
*****************************************
*MATERIAL,NAME=ONE
*ELASTIC
1.,.1
**********************************
*BOUNDARY
cax4cinax4.LHS,1
cax4.LHS,1
cax4.fullyfixed,1,2,0.0
analytic.all,1,
**********************************
*STEP
*STATIC
*RESTART,WRITE
*CLOAD
cax4cinax4.1,2,-1.
cax4.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=cax4cinax4.LHS
U,
*NODE PRINT,NSET=cax4cinax4.LHS
COORD,U
*NODE FILE,NSET=cax4.LHS
U,
*NODE PRINT,NSET=cax4.LHS
COORD,U
*NODE PRINT,NSET=analytic.all
COORD,U
*OUTPUT,FIELD
*NODE OUTPUT
U,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=cax4cinax4.LHS
U,
*NODE OUTPUT,NSET=cax4.LHS
U,
*NODE OUTPUT,NSET=analytic.all
U,
*END STEP