*HEADING
: AQUA Buoyancy [R3D4]

 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

*******************
** R3D4 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=REF
 109, 0.0, 0.0, 5.0
*ELEMENT,ELSET=ER3D4,TYPE=R3D4
 101,    101, 103, 104, 102
 102,    103, 105, 106, 104
 103,    105, 107, 108, 106
 104,    107, 101, 102, 108
 105,    102, 104, 106, 108
 106,    107, 105, 103, 101
*RIGID BODY, ELSET=ER3D4,REF=109
 1.0, 
*AQUA
 -200.0, 5.0, 9.8, 1026.05
*ELSET,ELSET=EALL
 ER3D4,
*ELSET,ELSET=SIDE
101,
**************************************
*STEP,NLGEOM
 Heel: 0 degrees. Trim: 0 degrees.  Side only.  RF2= -1.2569E7.
*STATIC
** load on side only.
*DLOAD
 SIDE, PB, 1.0
*BOUNDARY
 REF, 1, 6
*NODE PRINT,NSET=REF
 U3, UR1, UR2, RF1, RF2, RF3, RM1, RM2
*OUTPUT,FIELD
*NODE OUTPUT,NSET=REF
U,
RF,
*END STEP
*STEP,NLGEOM
 Heel: 0 degrees.  Trim: 0 degrees. Ref: RF3 = -1.508E8
*STATIC
*DLOAD,OP=NEW
 EALL, PB, 1.0
*BOUNDARY
 REF, 1, 6
*NODE PRINT,NSET=REF
 U3, UR1, UR2, RF3, RM1, RM2
*OUTPUT,FIELD
*NODE OUTPUT,NSET=REF
U,
RF,
*END STEP
**
*STEP,NLGEOM
 Heel: 60 degrees.  Trim: 0 degrees. Ref: RM1 = -5.167E8  RM2 = 0.
*STATIC
*BOUNDARY
 REF, 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
 REF, 5, ,  0.08725
 REF, 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
 REF, 5, ,  0.1745
 REF, 6, ,  0.3023
*END STEP