y coordinates of center of buoyancy

Hello,

I have plotted a ship that is partially immersed in water.

With the following code, I calculated the area of the ship :

ship=polyshape(x,y)
figure(1)
hold on 
plot(ship)
axis equal
plot(0,0,'*')
[xc,yc]=centroid(ship)
plot(centroid(ship),'x','Color','k')
rectangle=polyshape([-12 -12 12 12],[-6 6 6 -6])
plot(rectangle)
area(intersect(rectangle,coque_bateau))

Now, I need to find the y coordinate of the centroid of the immersed volume

I only know the total value of the immersed volume is Vim = 34,14 m^3

How can I find the centroid ? I am guessing I need to use again the intersect method, but I don't see how

I hope my question is clear enough, thank you in advance

Answer 1

[Are Mjaavatten]

Is this what you want?

x = [-8,-8,-6,0,6,8,8,-8];
y = [10,5,1,0,1,5,10,10];
ship=polyshape(x,y);
figure(1)
clf
hold on 
plot(ship)
axis equal
[xc,yc]=centroid(ship);
plot(xc,yc,'xk')
rectangle=polyshape([-12 -12 12 12],[-6 6 6 -6]);
plot(rectangle)
submerged = intersect(rectangle,ship);
[xs,ys] = centroid(submerged);
A = area(submerged);
plot(xs,ys,'ok')
fprintf('Area = %.2f\n',A);
fprintf('Centroid y = %.2f\n',ys)