*HEADING
: AQUA Buoyancy [R3D3]

 RIGID RECTANGULAR BOX:    HEEL THEN TRIM 

 A rigid 3D box is subject to AQUA buoyancy load type PB.
 The dimensions of the box are:
   
 Length = 100m
 Width  =  30m
 Detph  =  10m
 Draft  =   5m

In Step 1: only one side of the box is loaded. The expected reaction
force in the 2-direction is analytically determined to be 1.2569 E7 N.

In following steps: all of the box is loaded.
The buoyancy force is always  rho g Length Width Draft = 1.508E8 N,
i.e. RF3= -1.508E8. The body is then given 60 degreed of heel,
followed by 0, 10, and 20 degrees of trim. The reference values
for moment listed here agree within one quarter of 1% with values
produced by an independent code.

Heel: 0 degrees.   Trim: 0 degrees. (Side only.)  RF2= -1.2569E7.
Heel: 0 degrees.   Trim: 0 degrees.      Ref: RF3 = -1.508E8
Heel: 60 degrees.  Trim: 0 degrees.      Ref: RM1 = -5.167E8  RM2 = 0.
Heel: 60 degrees.  Trim: 10 degrees.     Ref: RM1 = -3.708E8, RM2 = -8.189E8 
Heel: 60 degrees.  Trim: 0 - 20 degrees. Ref: RM1 =  4.922E7, RM2 = -1.539E9

*******************
** R3D3 MODEL
*******************
*NODE,NSET=NR3D3
 101,   -50.0, -15.0, 10.0
 102,    50.0, -15.0, 10.0
 103,   -50.0, -15.0,  0.0
 104,    50.0, -15.0,  0.0
 105,   -50.0,  15.0,  0.0
 106,    50.0,  15.0,  0.0
 107,   -50.0,  15.0, 10.0
 108,    50.0,  15.0, 10.0
*NODE,NSET=QA_TEST
 109, 0.0, 0.0, 5.0
*ELEMENT,ELSET=ER3D3,TYPE=R3D3
 101,   101, 103, 102
 102,   102, 103, 104
 103,   103, 105, 104
 104,   104, 105, 106
 105,   105, 107, 106
 106,   106, 107, 108
 107,   107, 101, 108
 108,   108, 101, 102
 109,   102, 104, 106
 110,   102, 106, 108
 111,   101, 105, 103
 112,   101, 107, 105
*RIGID BODY, ELSET=ER3D3,REF=109
 1.0, 
*AQUA
 -200.0, 5.0, 9.8, 1026.05
*ELSET,ELSET=EALL
 ER3D3, 
*ELSET,ELSET=SIDE
101,102,
**************************************
*STEP,NLGEOM
 Heel: 0 degrees. Trim: 0 degrees.  Side only.  RF2= -1.2569E7.
*STATIC
** load on side only.
*DLOAD
 SIDE, PB, 1.0
*BOUNDARY
 QA_TEST, 1, 6
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
U3,UR1, UR2
RF1,RF2,RF3
RM1, RM2
*END STEP
*STEP,NLGEOM
 Heel: 0 degrees.  Trim: 0 degrees. Ref: RF3 = -1.508E8
*STATIC
*DLOAD,OP=NEW
 EALL, PB, 1.0
*BOUNDARY
 QA_TEST, 1, 6
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
U,
RF3, RM1, RM2
*END STEP
**
*STEP,NLGEOM
 Heel: 60 degrees.  Trim: 0 degrees. Ref: RM1 = -5.167E8  RM2 = 0.
*STATIC
*BOUNDARY
 QA_TEST, 4, , 1.0472
*END STEP
**
*STEP,NLGEOM
 Heel: 60 degrees.  Trim: 10 degrees.  Ref: RM1 = -3.708E8, RM2 = -8.189E8 
*STATIC
 1.0,  1.0
*BOUNDARY
 QA_TEST, 5, ,  0.08725
 QA_TEST, 6, ,  0.15115
*END STEP
*STEP,NLGEOM
 Heel: 60 degrees.  Trim: 0 - 20 degrees.  Ref: RM1 =  4.922E7, RM2 = -1.539E9
*STATIC
 1.0,  1.0
*BOUNDARY
 QA_TEST, 5, ,  0.1745
 QA_TEST, 6, ,  0.3023
*END STEP