Stiffness matrix and translation

[mancellin]

I'm not actually sure if this is a bug or a feature that I don't understand:

from meshmagick.mmio import load_mesh
from meshmagick.mesh import Mesh
from meshmagick.hydrostatics import Hydrostatics

reference_mesh = Mesh(*load_mesh("SEAREV.hst", 'hst'))

for shift in range(0, 100, 20):
    shifted_mesh = reference_mesh.copy()
    shifted_mesh.translate_x(shift)
    print shift
    print Hydrostatics(shifted_mesh, cog=(shift, 0, 0)).hs_data['stiffness_matrix'][2, 2]
    print Hydrostatics(shifted_mesh, cog=(shift, 0, 0)).hs_data['Iyy']

The inertia matrix is invariant when the body is translated in the x direction.
I would have expected the same to be true for the stiffness matrix.
Are the rotation dofs defined around the same axis in both cases ?

[mancellin]

Actually it seems that the x and y coordinates given in the cog argument have no influence on the results.

[ryancoe]

Looking at this a bit more, here's my current understanding:

stiffness: I think there was a bit of a problem here, which I've submitted a fix for with #22 (note however that what I've done is not yet propagated back to the Hydrostatics class; if you think it looks good we can finish doing this).

hydrostatics: these look like they are currently calculated relative to a point on a vertical line passing through origin but with the correct vertical position:

This section uses the x and y location of the origin (xi, yi = polyverts[0, :2])

meshmagick/meshmagick/hydrostatics.py

Lines 605 to 621 in 0d214f9

	xi, yi = polyverts[0, :2]
	for (xii, yii) in polyverts[1:, :2]:
	dx = xii - xi
	dy = yii - yi
	px = xi + xii
	py = yi + yii
	# a = xi * xi + xii * xii
	sigma0 += dy * px
	sigma1 += dy * (px * px - xi * xii)
	sigma2 += dx * (py * py - yi * yii)
	# sigma3 += dy * (py * a + 2 * px * (xi * yi + xii * yii))
	sigma3 += dy * (py * px * px + yi * xi * xi + yii * xii * xii)
	sigma4 += dy * (xi * xi + xii * xii) * px
	sigma5 += dx * (yi * yi + yii * yii) * py
	xi, yi = xii, yii

This section uses the vertical location of the cog (self.zg)

meshmagick/meshmagick/hydrostatics.py

Lines 653 to 655 in 0d214f9

	a = self.zg - zb # BG
	gm_x = transversal_metacentric_radius - a
	gm_y = longitudinal_metacentric_radius - a

Thus, it might be appropriate to change this method so that it calculates based on COG, right?

meshmagick/meshmagick/hydrostatics.py

Line 599 in 0d214f9

polyverts = clipper.clipped_crown_mesh.vertices[polygon]

This line would be changed to polyverts = clipper.clipped_crown_mesh.vertices[polygon] - self._gravity_center

[ryancoe]

related to #1

[frongere]

I just released v3.0 which should fix these issues.

Note however that I have totally refactored the hydrostatics module using only functions with no class. The code is far more simple and the number of lines of code has lowered.

The stiffness matrix and metacentric stuffs are now expressed on a line passing though COG/buoyancy centre.

Positions such as buoyancy centre are expressed in the initial frame. It allows to directly compare results to other commercial software such as Maat Hydro.