*HEADING
BOUSSINESQ PROBLEM, 12 CAX8R + 4 CINAX5R & CAX8 only.
*preprint,model=yes,hist=yes
**************************************************
*part,name=analytic
*node,nset=all
*element,type=b21,elset=dummy
*beamsection,material=one,elset=dummy,sec=circ
*endpart
************************************************
*part,name=cax8cinax5
*NODE
*NGEN
*NSET,NSET=LHS,GENERATE
*ELEMENT,TYPE=CAX8R,ELSET=ALL
*ELGEN,ELSET=ALL
*ELEMENT,TYPE=CINAX5R,ELSET=ALL
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
*****************************************
*part,name=cax8
*NODE
**26,0.,-6.
**65,2.8333333,-2.8333333
**66,4.2426407,-4.2426407
**105,4.,0.
**106,6.,0.
*NGEN
*NSET,NSET=LHS,GENERATE
*NSET,NSET=FULLYFIXED,GEN
*ELEMENT,TYPE=CAX8R,ELSET=ALL
*ELGEN,ELSET=ALL
*SOLID SECTION,ELSET=ALL,MATERIAL=ONE
*endpart
*****************************************
*assembly,name=all
*instance,name=analytic,part=analytic
*endinstance
*instance,name=cax8cinax5,part=cax8cinax5
*endinstance
*instance,name=cax8,part=cax8
*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