Bug Beamdyn: center of mass and moment of inertia in summary file are incorrect

[andrew-platt]

Bug description

The moment of inertia about the z-axis for a straight beam is reported by BeamDyn to be zero. This is incorrect indicating an issue in the calculations performed BD_ComputeElementMass routine. Further inspection of the routine also shows that the x_cm and y_cm terms may not be accounted for in calculating the center of gravity in an element (they may be buried in the p%uu0 terms.

To Reproduce

Run standalone BeamDyn with the 5 MW straight blade. Look at the summary file to see that

Blade mass moment of inertia: 
       1.11816E+07       0.00000E+00       0.00000E+00
       0.00000E+00       1.11816E+07       0.00000E+00
       0.00000E+00       0.00000E+00       0.00000E+00

Expected behavior

We expect to see some moment of inertia reported about the z-axis ( (3,3) term).

Notes

This does not affect internal calculations, but only the values reported in the summary file. These values are not used anywhere else in the code. This bug is present in all versions of BeamDyn to date.

[andrew-platt]

See issue #1337 for more complete writeup.

Closing this as duplicate.