MATLAB运用——计算船舶浮心（切面法）

已于 2022-03-13 11:29:15 修改
阅读量3.6k 收藏 29
点赞数 9
分类专栏： MATLAB学习 文章标签： matlab 开发语言 算法 3d
于 2022-03-13 11:25:04 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/go_bananas/article/details/123455457
版权

在上一次的运用中，我们计算了水花号的质心。
在这一篇文章中，我们的任务是计算在任意水线下，浪花号的浮心。

问题：

计算浪花号在任意水线下的浮心。

分析

  1. 船舶的浮心就是它排开水的部分的水的质心。所以计算浮心的问题最终还是回归到计算质心上来。
  2. 我们可以利用我们平时计算三维物体体积的积分的方式计算浮心：
    简单的说就是（原谅我的灵魂画图）
    在这里插入图片描述
  • 我们把船切成一个一个平面（紫色的线）
  • 把【平面的面积X区间长度】近似地看作一个区间的排开水的体积（红色区域）
  • 把每个区域的的水的质心算出来（橙色的点）
  • 最后算出整体的质心就是船的浮心了（蓝色的点）

解决过程（完整代码在最后）：

构建船体曲线

这次构造船体曲线和以往有些不同：
我们需要使用多项式拟合：原因是我们需要得到拟合图像的确切的公式。
有了这个公式我们才能计算出截面图形（定X变换Y得到Z）
同以前一样：我们把数据导入：
在这里插入图片描述
选择多项式拟合（我选择的是x，y各3次方）：
在这里插入图片描述
我们就可以得到图像：
在这里插入图片描述
最关键的是得到拟合的公式：
在这里插入图片描述
我们把这个公式写入matlab里面：我们就可以由x，y得到z了。

function z=caculate(x,y)       
       p00 =       24.74;
       p10 =     -0.7374;
       p01 =  -1.124e-16; 
       p20 =    0.002366; 
       p11 =    1.44e-18; 
       p02 =     0.01779; 
       p30 =  -1.795e-06; 
       p21 =  -3.016e-21; 
       p12 =  -1.137e-05; 
       p03 =  -1.453e-20; 
        z = p00 + p10*x + p01*y + p20*x^2 + p11*x*y + p02*y^2 + p30*x^3 + p21*x^2*y 
                    + p12*x*y^2 + p03*y^3;
end

定义各项数据：（所有代码都是MATLAB下运行）

船体与水的相关：

clear all
boat.L = 400;          %船的长度
boat.W = 100;          %船的宽度
boat.HB = boat.W / 2;  %船的半宽
boat.D = 60;           %船的深度
boat.HD = boat.D / 2;  %船的半宽

max_area = boat.D * boat.W     %沿X轴的切面（Y_Z平面切面）的最大面积
max_volume = max_area * boat.L %立方体的体积
density_water = 28.312;          %水的密度 （千克/立方英尺）
max_mass = max_volume * density_water  %排开水的最大质量

微元相关

dy = 1;   % meters（微元_Y轴）
dz = 1;   % meters（微元_X轴）

mesh.ys = -boat.HB:dy:boat.HB;        % meters  Y轴切片
mesh.zs = -boat.HD:dz:boat.HD;        % meters  Z轴切片

[mesh.ygrid,mesh.zgrid] = meshgrid(mesh.ys,mesh.zs);

total_area = boat.W * boat.D       %平方英尺
mesh.dA = total_area / numel(mesh.ygrid)  %微元面积

浮心相关：

%XMass,YMass,ZMass记录的是每个切片的浮心（2D）
ZMass=[];
XMass=[];
YMass=[];
MMass=[];%每个浮心点的权重（质量）
tot_mass=0;%总质量

切片

%对每个x所对应的点切片：
for i=2 : 11   %对每个x积分
    ZZ=[];%记录该X下的Z轴坐标
    for j=1 : 101 %讨论y从[-60,60]
        t=caculate(i*36,mesh.ys(1,j));
        ZZ=[ZZ t];
    end
    y = mesh.ys;
    figure(i);%绘制不同的图
    
    hull = mesh.zgrid > ZZ%在所计算的Z的值之上的就是船的截面
    redmap = [1,1,1;1,0,0];
    colormap(redmap);
    image(mesh.ys,mesh.zs,flipud(hull),'AlphaData',0.5);

我们就可以得到各个截面的图像了：
在这里插入图片描述
在这里插入图片描述
等等共11个图像

描述水线

   %描述水线 y=kx+d
    d = 0.05; 
    theta = 10;
    y = mesh.ys;
    z = tand(theta) .* y + d;%每个水线对应的方程
	water = mesh.zgrid < z;
	bluemap = [1,1,1; 0,0,1];
    plotMatrix(water,mesh,bluemap);%画出水线

在这里插入图片描述

计算水线和船体的公共面积

%计算水线与切面的公共面积
    sub_region = hull & water;%&位与运算（都是1才得1）
    purplemap = [1,1,1; 1,0,1];
    plotMatrix(sub_region,mesh,purplemap);%绘制图像

运行后得到图像：
在这里插入图片描述
等11个图像

分别计算每个区域的质心（2D）

 %计算浮心2D
    COB = centerOfMass2(sub_region,mesh);

function M = matrixSum(A)
    % matrixSum: returns total of all elements in the matrix
    % A: matrix
    % returns: scalar
    
    % normally sum(m) computes the sums of the columns
    % selecting m(:) flattens the matrix and computes the sum of all elements
    % see https://stackoverflow.com/questions/1721987/what-are-the-ways-to-sum-matrix-elements-in-matlab
    M = sum(A(:));
end

function COM = centerOfMass2(masses,mesh)
    % centerOfMass2: computes center of mass in 2D
    % masses: matrix of masses
    % mesh: structure containing ygrid and zgrid
    % returns: Vector [ycom,zcom]
    M = matrixSum(masses);
    ycom = matrixSum(masses .* mesh.ygrid) / M;
    zcom = matrixSum(masses .* mesh.zgrid) / M;
    COM = [ycom,zcom];
end

再根据2D的质心计算浮心（3D）

ANS=centerOfMass3(XMass,YMass,ZMass,MMass,tot_mass);%计算三维浮心

function ANS=centerOfMass3(XMass,YMass,ZMass,MMass,m)
x=0;
y=0;
z=0;
    for i=1:10
        x=x+XMass(1,i)*MMass(1,i)/m;
        y=y+YMass(1,i)*MMass(1,i)/m;
        z=z+YMass(1,i)*MMass(1,i)/m;
    end
    ANS=[x,y,z];
end

最后得到浪花号在z=tan(10)y+0.05的水线下的浮心是：
在这里插入图片描述

完整代码：（.mlx）

这是用多项式拟合出来的函数图像数据
X = [ 0.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 36.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 72.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 108.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 144.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 180.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 216.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 252.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 288.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 324.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 360.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 396.000000 432.000000 ];
Y = [ 0.000000 -27.750000 -27.750000 -17.000000 -16.000000 -9.750000 -7.500000 7.500000 9.750000 16.000000 17.000000 27.750000 27.750000 -41.750000 -41.750000 -39.375000 -36.500000 -32.500000 -29.500000 -24.375000 -16.500000 -9.750000 9.750000 16.500000 24.375000 29.500000 32.500000 36.500000 39.375000 41.750000 41.750000 -51.250000 -50.875000 -50.875000 -49.750000 -46.250000 -41.250000 -36.250000 -27.500000 -16.250000 -6.250000 6.250000 16.250000 27.500000 36.250000 41.250000 46.250000 49.750000 50.875000 50.875000 51.250000 -57.250000 -57.250000 -55.250000 -54.000000 -51.875000 -48.625000 -44.125000 -35.250000 -25.500000 -9.250000 9.250000 25.500000 35.250000 44.125000 48.625000 51.875000 54.000000 55.250000 57.250000 57.250000 -59.750000 -58.750000 -56.750000 -56.250000 -56.250000 -53.500000 -51.125000 -41.750000 -27.500000 -14.750000 14.750000 27.500000 41.750000 51.125000 53.500000 56.250000 56.250000 56.750000 58.750000 59.750000 -58.750000 -58.250000 -58.000000 -56.500000 -56.500000 -55.250000 -49.250000 -43.000000 -28.750000 -17.500000 17.500000 28.750000 43.000000 49.250000 55.250000 56.500000 56.500000 58.000000 58.250000 58.750000 -58.500000 -58.500000 -57.250000 -57.000000 -56.750000 -53.875000 -50.000000 -46.500000 -16.250000 -15.250000 15.250000 16.250000 46.500000 50.000000 53.875000 56.750000 57.000000 57.250000 58.500000 58.500000 -55.750000 -55.750000 -55.000000 -54.750000 -53.500000 -49.000000 -43.500000 -33.000000 -19.500000 -10.750000 10.750000 19.500000 33.000000 43.500000 49.000000 53.500000 54.750000 55.000000 55.750000 55.750000 -50.500000 -50.500000 -50.000000 -48.500000 -48.000000 -45.500000 -33.500000 -21.000000 -11.500000 -8.750000 8.750000 11.500000 21.000000 33.500000 45.500000 48.000000 48.500000 50.000000 50.500000 50.500000 -48.000000 -48.000000 -47.500000 -45.875000 -43.000000 -34.000000 -16.875000 -10.500000 -6.500000 -4.500000 4.500000 6.500000 10.500000 16.875000 34.000000 43.000000 45.875000 47.500000 48.000000 48.000000 -42.625000 -41.375000 -41.000000 -41.000000 -33.500000 -11.500000 11.500000 33.500000 41.000000 41.000000 41.375000 42.625000 0.000000 ];
Z = [ 24.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -6.000000 0.000000 6.000000 12.000000 24.000000 18.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 24.000000 18.000000 12.000000 18.000000 24.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 18.000000 24.000000 12.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 18.000000 24.000000 12.000000 6.000000 0.000000 24.000000 18.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 18.000000 24.000000 0.000000 6.000000 12.000000 12.000000 6.000000 0.000000 24.000000 18.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 18.000000 24.000000 0.000000 6.000000 12.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 18.000000 24.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 18.000000 24.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 18.000000 24.000000 24.000000 18.000000 12.000000 6.000000 0.000000 -6.000000 -12.000000 -18.000000 -24.000000 -30.000000 -30.000000 -24.000000 -18.000000 -12.000000 -6.000000 0.000000 6.000000 12.000000 24.000000 18.000000 6.000000 12.000000 24.000000 18.000000 0.000000 -6.000000 -6.000000 0.000000 24.000000 18.000000 12.000000 6.000000 24.000000 ];
[a,b]=createFit(X,Y,Z)
axis('equal');%使得各个坐标轴长度相等
zlim([-40,40]);%限定Z轴范围[-40,40]
计算浮心(单位：英尺)：
clear all
boat.L = 400;          %船的长度
boat.W = 100;          %船的宽度
boat.HB = boat.W / 2;  %船的半宽
boat.D = 60;           %船的深度
boat.HD = boat.D / 2;  %船的半宽

max_area = boat.D * boat.W     %沿X轴的切面（Y_Z平面切面）的最大面积
max_volume = max_area * boat.L %立方体的体积
density_water = 28.312;          %水的密度 （千克/立方英尺）
max_mass = max_volume * density_water  %排开水的最大质量
从Y,Z轴切片（取微元）：
dy = 1;   % meters（微元_Y轴）
dz = 1;   % meters（微元_X轴）

mesh.ys = -boat.HB:dy:boat.HB;        % meters  Y轴切片
mesh.zs = -boat.HD:dz:boat.HD;        % meters  Z轴切片

[mesh.ygrid,mesh.zgrid] = meshgrid(mesh.ys,mesh.zs);

total_area = boat.W * boat.D       %平方英尺
mesh.dA = total_area / numel(mesh.ygrid)  %微元面积
XMass,YMass,ZMass记录的是每个切片的浮心（2D）
ZMass=[];
XMass=[];
YMass=[];
MMass=[];%每个浮心点的权重（质量）
tot_mass=0;%总质量
对每个x所对应的点切片：
for i=2 : 11   %对每个x积分
    ZZ=[];%记录该X下的Z轴坐标
    for j=1 : 101 %讨论y从[-60,60]
        t=caculate(i*36,mesh.ys(1,j));
        ZZ=[ZZ t];
    end
    y = mesh.ys;
    figure(i);%绘制不同的图
    
    hull = mesh.zgrid > ZZ%在所计算的Z的值之上的就是船的截面
    redmap = [1,1,1;1,0,0];
    colormap(redmap);
    image(mesh.ys,mesh.zs,flipud(hull),'AlphaData',0.5);
    
    %描述水线 y=kx+d
    d = 0.05; 
    theta = 10;
    y = mesh.ys;
    z = tand(theta) .* y + d;%每个水线对应的方程
    
   
    water = mesh.zgrid < z;
    bluemap = [1,1,1; 0,0,1];
    plotMatrix(water,mesh,bluemap);
     %计算水线与切面的公共面积
    
     sub_region = hull & water;%&位与运算（都是1才得1）
    purplemap = [1,1,1; 1,0,1];
    plotMatrix(sub_region,mesh,purplemap);%绘制图像
    
    %计算浮心2D
    COB = centerOfMass2(sub_region,mesh);
    ZMass=[ZMass COB(1,1)];
    YMass=[YMass COB(1,2)];
    XMass=[XMass (i-0.5)*36];
    MMass=[MMass matrixSum(sub_region)];
    tot_mass=tot_mass+matrixSum(sub_region);
end
ANS=centerOfMass3(XMass,YMass,ZMass,MMass,tot_mass);%计算三维浮心
ANS

function ANS=centerOfMass3(XMass,YMass,ZMass,MMass,m)
x=0;
y=0;
z=0;
    for i=1:10
        x=x+XMass(1,i)*MMass(1,i)/m;
        y=y+YMass(1,i)*MMass(1,i)/m;
        z=z+YMass(1,i)*MMass(1,i)/m;
    end
    ANS=[x,y,z];
end

function M = matrixSum(A)
    % matrixSum: returns total of all elements in the matrix
    % A: matrix
    % returns: scalar
    
    % normally sum(m) computes the sums of the columns
    % selecting m(:) flattens the matrix and computes the sum of all elements
    % see https://stackoverflow.com/questions/1721987/what-are-the-ways-to-sum-matrix-elements-in-matlab
    M = sum(A(:));
end

function COM = centerOfMass2(masses,mesh)
    % centerOfMass2: computes center of mass in 2D
    % masses: matrix of masses
    % mesh: structure containing ygrid and zgrid
    % returns: Vector [ycom,zcom]
    M = matrixSum(masses);
    ycom = matrixSum(masses .* mesh.ygrid) / M;
    zcom = matrixSum(masses .* mesh.zgrid) / M;
    COM = [ycom,zcom];
end
function plotMatrix(A,mesh,cmap)
    % plotMatrix: plots a matrix using image
    % A: matrix
    % mesh: srtruct containing ys and zs
    % cmap: Colormap
    colormap(cmap);
    image(mesh.ys,mesh.zs,flipud(A),'AlphaData',0.5);
end

function [fitresult, gof] = createFit(X, Y, Z)
%CREATEFIT(X,Y,Z)
%  Create a fit.
%
%  Data for 'untitled fit 1' fit:
%      X Input : X
%      Y Input : Y
%      Z Output: Z
%  Output:
%      fitresult : a fit object representing the fit.
%      gof : structure with goodness-of fit info.
%
%  另请参阅 FIT, CFIT, SFIT.

%  由 MATLAB 于 12-Mar-2022 16:36:02 自动生成


%% Fit: 'untitled fit 1'.
[xData, yData, zData] = prepareSurfaceData( X, Y, Z );

% Set up fittype and options.
ft = fittype( 'poly33' );

% Fit model to data.
[fitresult, gof] = fit( [xData, yData], zData, ft, 'Normalize', 'on' );

% Plot fit with data.
figure( 'Name', 'untitled fit 1' );
h = plot( fitresult, [xData, yData], zData );
legend( h, 'untitled fit 1', 'Z vs. X, Y', 'Location', 'NorthEast', 'Interpreter', 'none' );
% Label axes
xlabel( 'X', 'Interpreter', 'none' );
ylabel( 'Y', 'Interpreter', 'none' );
zlabel( 'Z', 'Interpreter', 'none' );
grid on
view( 10.1, 31.6 );
end
多项式函数拟合部分：
Linear model Poly33:
     f(x,y) = p00 + p10*x + p01*y + p20*x^2 + p11*x*y + p02*y^2 + p30*x^3 + p21*x^2*y 
                    + p12*x*y^2 + p03*y^3
Coefficients (with 95% confidence bounds):
       p00 =       24.74  (18.55, 30.93)
       p10 =     -0.7374  (-0.8701, -0.6048)
       p01 =  -1.124e-16  (-0.1365, 0.1365)
       p20 =    0.002366  (0.001669, 0.003062)
       p11 =    1.44e-18  (-0.001667, 0.001667)
       p02 =     0.01779  (0.01487, 0.02072)
       p30 =  -1.795e-06  (-2.821e-06, -7.681e-07)
       p21 =  -3.016e-21  (-3.618e-06, 3.618e-06)
       p12 =  -1.137e-05  (-2.312e-05, 3.759e-07)
       p03 =  -1.453e-20  (-3.92e-05, 3.92e-05)
 
function z=caculate(x,y)       
       p00 =       24.74;
       p10 =     -0.7374;
       p01 =  -1.124e-16; 
       p20 =    0.002366; 
       p11 =    1.44e-18; 
       p02 =     0.01779; 
       p30 =  -1.795e-06; 
       p21 =  -3.016e-21; 
       p12 =  -1.137e-05; 
       p03 =  -1.453e-20; 
        z = p00 + p10*x + p01*y + p20*x^2 + p11*x*y + p02*y^2 + p30*x^3 + p21*x^2*y 
                    + p12*x*y^2 + p03*y^3;
end

可待提升的空间：

  1. 多项式拟合可以使用更高的次数，结果会更精确。
  2. 船的切片可以切得更密集一些，这样进度会更高。下面是制造更多切面的代码（C++）
#include<stdio.h>
#include<bits/stdc++.h>
using namespace std;
double XY_Plane[22][12]=//这是船的 XY平面 
{   //Y轴负半轴 
   0,0,0,0,0,0,0,0,0,0,0,0, 
   0,0,	0,	6.25,	9.25,	14.75,	17.5,	15.25,	10.75,	8.75,	4.5,	0,
   0,0,	9.75,	16.25,	25.5,	27.5,	28.75,	16.25,	19.5,	11.5,	6.5,	0,
   0,0,	16.5,	27.5,	35.25,	41.75,	43,	46.5,	33,	21,	10.5,	0,
   0,0,	24.375,	36.25,	44.125,	51.125,	49.25,	50,	43.5,	33.5,	16.875,	0,
   0,7.5,	29.5,	41.25,	48.625,	53.5,	55.25,	53.875,	49,	45.5,	34,	11.5,
   0,9.75,	32.5,	46.25,	51.875,	56.75,	58,	56.75,	53.5,	48,	43,	33.5,//DWL（水平面） XY_Plane[6][N]
   0,16,	36.5,	49.75,	54,	58.75,	58.25,	57,	54.75,	48.5,	45.875,	42.625,
   0,17,	39.375,	51.25,	55.25,	59.75,	58.75,	57.25,	55,	50,	47.5,	41.375,
   0,27.75,	41.75,	50.875,	57.25,	56.25,	56.5,	58.5,	55.75,	50.5,	48,	41,
   0,27.75,	41.75,	50.875,	57.25,	56.25,	56.5,	58.5,	55.75,	50.5,	48,	41,
   0,0,0,0,0,0,0,0,0,0,0,0,//----------------------------------->x轴 XY_Plane[11][N] 
   0,27.75,	41.75,	50.875,	57.25,	56.25,	56.5,	58.5,	55.75,	50.5,	48,	41,
   0,27.75,	41.75,	50.875,	57.25,	56.25,	56.5,	58.5,	55.75,	50.5,	48,	41,
   0,17,	39.375,	51.25,	55.25,	59.75,	58.75,	57.25,	55,	50,	47.5,	41.375,
   0,16,	36.5,	49.75,	54,	58.75,	58.25,	57,	54.75,	48.5,	45.875,	42.625,
   0,9.75,	32.5,	46.25,	51.875,	56.75,	58,	56.75,	53.5,	48,	43,	33.5,//DWL（水平面） XY_Plane[16][N]
   0,7.5,	29.5,	41.25,	48.625,	53.5,	55.25,	53.875,	49,	45.5,	34,	11.5,
   0,0,	24.375,	36.25,	44.125,	51.125,	49.25,	50,	43.5,	33.5,	16.875,	0,
   0,0,	16.5,	27.5,	35.25,	41.75,	43,	46.5,	33,	21,	10.5,	0,
   0,0,	9.75,	16.25,	25.5,	27.5,	28.75,	16.25,	19.5,	11.5,	6.5,	0,
   0,0,	0,	6.25,	9.25,	14.75,	17.5,	15.25,	10.75,	8.75,	4.5,	0,//21
   //Y轴正半轴 
};//每一行代表一个水线 
vector<double> XX,YY,ZZ;//存储算好的X,Y,Z数据  
double Z[22][12];//Z[i][j]表示在第一张图中第i行第j列的点的高度
void Insert_Num(){
   for(int i=1;i<=21;i++){
   	if(i==11)continue;
   	for(int j=1;j<11;j++){
   		if(XY_Plane[i][j]&&XY_Plane[i][j+1]){
   			XX.push_back( (j+0.5)*36);
   			YY.push_back( (XY_Plane[i][j]+XY_Plane[i][j+1])/2 );
   			ZZ.push_back( Z[i][j] );
   		}
   	}
   }
}
int main(){
   //高于DWL 
   for(int i=12,k=4;i<=15;i++,k--){
   	for(int j=1;j<=11;j++){
   		if(XY_Plane[i][j]==0)continue;
   		//Z为正值 
   		Z[i][j]=(k)*6;//Y轴正半轴 
   		Z[22-i][j]=Z[i][j];//Y轴负半轴
   	} 
   }
   //低于DWL
   for(int i=17;i<=21;i++){
   	for(int j=1;j<=11;j++){
   		if(XY_Plane[i][j]==0)continue;
   		//Z为负值 
   		Z[i][j]=-1*(i-16)*6;//Y轴正半轴 
   		Z[22-i][j]=Z[i][j];//Y轴负半轴 
   	}
   }
   //导出数据
   
   for(int j=1;j<=11;j++){//列 
   	for(int i=12;i<=21;i++){//行 
   		if(XY_Plane[i][j]){
   			XX.push_back(j*36);
   			YY.push_back(XY_Plane[i][j]);
   			ZZ.push_back(Z[i][j]);//Y轴正半轴 
   			
   			XX.push_back(j*36);
   			YY.push_back(-XY_Plane[i][j]);
   			ZZ.push_back(Z[i][j]);//Y轴负半轴  
   		}
   	}
   }
   //(0,0,288)最左顶点 
   XX.push_back(0); 
   YY.push_back(0);
   ZZ.push_back(24);
   //最右顶点
   XX.push_back(12*36);
   YY.push_back(0);
   ZZ.push_back(24); 
   
   Insert_Num(); 
   
   	//导出X轴
   printf("X = [ ");
   for(int i=0;i<XX.size();i++){
   	printf("%lf ",XX[i]);
   }
   printf("];\n");
   
   //导出Y轴 
   printf("Y = [ ");
   for(int i=0;i<YY.size();i++){
   	printf("%lf ",YY[i]);
   }
   printf("];\n");
   
   //导出Z轴 
   printf("Z = [ ");
   for(int i=0;i<ZZ.size();i++){
   	printf("%lf ",ZZ[i]);
   }
   printf("];\n");
}
go_bananas
关注
  • 9
    点赞
  • 踩
  • 29
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
一种通过maxsurf计算稳心的方法
wz199399的博客
06-08 759
最近在做浮体的稳性分析的时候，客户要求计算浮体的稳心。传统的方法是通过计算在倾斜角度10°时候的浮力的作用线与中线的交点来代稳心。
NX二次开发-UFUN导出STL函数UF_STD_put_solid_in_stl_file博客文章源代码
10-21
NX二次开发-UFUN导出STL函数UF_STD_put_solid_in_stl_file博客文章源代码
CuraEngine-master_STL_开源_切片_stl切片源代码_
09-30
开源的STL切片源码，适合大部分程序员代码研读
Matlab拟合曲面，计算质心浮心
northern_light_的博客
03-13 664
船体主要要素：主尺度、船型系数、尺度比代码实现x 长y 宽z 高 / 深dy = 0.01;% mesh.ys 还是一个struct数组 mesh . zs = 0 : dz : boat . D;meshgrid生成网格dA是是网格中一个单元格的面积，以平方米为单位。numel：返回数组 A 中的元素数目 n。
stlshow_stl分层_STL分层_stlmatlab_STL切片_stl分层
09-11
在matlab实现STL分层，自定义厚度，参考意义
stl读取.rar_STL文件_matlab_文件读取_读取stl文件
07-15
stl文件读取，用于读取stl文件，在进行其他操作
几何质心的定义和计算
04-20
几何质心的定义和计算
stl文件分层
01-23
可用于stl文件的快速分层，其分层速度相对比较快，应该可满足应用需求
MATLAB的Curve Fitting Tool 生成代码并被调用
热门推荐
liumangmao1314的博客
04-21 1万+
以MATLAB2014b为例： 打开Curve Fitting：应用程序>>Curve Fitting； 拟合：选择数据源>>X/Y/Zdata（没有的话运行下你想曲面拟合的代码即可） 生成代码：最上菜单栏>>CURVE FITTING TOOL>>文件>>Generate Code>>自动生成一个creatFit....
stlshow_stl分层_STL分层_stlmatlab_STL切片_stl分层.zip
10-10
stlshow_stl分层_STL分层_stlmatlab_STL切片_stl分层.zip
ThreeJS_STL_Loader
04-03
ThreeJS STL加载器 Projeto que carrega um modelo STL 第三方ThreeJS e STL： / 跨浏览器动画和。 Three.js美国 。 Código-fonte： STL以及3D印象深刻的格式。 Esse formato aproxima代表建模的超高级能力，而e信息代表了3D几何的卓越表现（无模型计算的CAD程序）。 Modelo 3D实用工具：ibudmen的 铁线虫 节点JS ExpressJS的 Nodemon :house_with_garden: :hotel: :department_store: :convenience_store: :Japanese_post_office: :house:
MATLAB "slice"切片功能实现四维显示
07-20
利用MATLAB 7.0的"slice"切片功能将三维表示的数据，通过对图形的线型、立面、色彩、渲染、光线、视角等的控制，可形象表现数据四维特性。
matlab代码设参数范围-NBUT-MADE1-Boat-Design:NBUT-MADE1-船设计
05-23
matlab 代码 设参数范围 NBUT 2020 MADE I Boat Design writer：zzzcd0x 思路 根据板材长度以及合理的长宽比等因素设定A和B的范围 第一步是粗筛，即枚举A的范围内每一个长度为0.6的区间的中点，计算中点位置的A、B、H对应的船体底面方程的重心、正浮吃水线、倾斜135°时的吃水线平面方程、浮心和复原力矩，最终判断该答案是否合理的条件是135°时复原力矩绝对值小于0.01 得到区间中点后手动选取几组理想的数据然后放入精筛，精筛中讲枚举选取的每一个中点对应的区间中的每一个值，筛选条件比粗筛多了长宽比、复原力矩绝对值大小要求小于0.001。 最终从精筛得到的数据中选取一组理想的作为最终数据 各个函数的作用 low_accuracy 粗筛主函数 high _accuracy 精筛主函数 getCOM 返回值为船体重心的z轴坐标 getCOB 返回值为储存135°时船体浮心的一维矩阵 getTOI 返回值为复原力矩 waterline 返回值为正浮吃水线高度 waterline135 返回值为135°时吃水线方程的截距 注意事项和优化方案 粗筛过程代码
遗传算法分层遗传算法.rar_GA_算法 C#_遗传算法 _遗传算法 分层_遗传算法’
07-14
遗传算法分层遗传算法
基于STL数据模型动态拓扑重构的快速切片算法
02-12
分析了现有金属激光立体成形(MLSF)切片算法,提出了一种基于STL模型动态拓扑重构的快速切片算法。根据STL模型中三角面片的几何信息和切片厚度,通过建立分组矩阵,减小了三角面片遍历的次数; 通过构建三角面片之间的局部动态拓扑关系,减小了切片平面与三角面片的求交计算次数; 并根据切片过程中大部分三角面片的毗邻关系不发生改变这一事实,提出了动态拓扑重构的算法,减小了切片过程中三角面片毗邻关系的查找次数,从而提高了切片算法的整体效率。在该算法的基础上,使用Visual C++和OpenGL开发了金属激光立体成形软件系统。
EXCELE计算船舶静水力模版
03-14
用EXCELE表格计算船舶静水力性能，包括浮心，排水体积，等
水下滑翔机的外形设计与水动力计算
11-15
水下滑翔机是一种不依靠外部推进装置的水下自主航行器， 它通过改变净浮力和浮心位置进行运动， 因此要求其具有优良的水动力性能。本文介绍实验室所研制的水下滑翔机的外形设计， 对水下滑翔机斜航运动与定常回转...
基于Solidworks二次开发的两栖车辆浮态设计及静稳性分析 (2009年)
05-26
以浮心计算为基础，通过对车辆的受力分析和运动分析，确定了车辆浮态的计算方法；利用SolidWorks的API函数进行二次开发，形成了计算静稳性的程序，并进行了实例计算。结果表明：与传统的作图法相比较，该方法计算...
鱼雷豚跳运动及初始转动角速度选取? (2015年)
05-22
附加质量、浮力、浮心、浸湿面积、浸水体积等变量随鱼雷出水姿态、出水过程不同而不同. 在构建鱼雷豚跳运动模型过程中, 充分考虑了各个变化物理量的导数项, 分析了鱼雷水动阻力系数与攻角的关系, 构建了运动模型, ...
matlab船舶六自由度代码
最新发布
04-28
以下是一个简单的Matlab船舶六自由度代码示例： ```matlab % 定义常数 g = 9.81; % 重力加速度 rho = 1025; % 海水密度 m = 100000; % 船舶质量 Iz = 20000000; % 船舶绕z轴转动惯量 B = 20; % 船宽 L = 100; % 船长 T = 10; % 船吃水深度 xg = 0; % 质心x坐标 yg = 0; % 质心y坐标 zg = 0; % 质心z坐标 xb = 0; % 浮心x坐标 yb = 0; % 浮心y坐标 zb = 0; % 浮心z坐标 % 定义初值 u0 = 0; % 初始速度 v0 = 0; % 初始横向速度 w0 = 0; % 初始垂向速度 p0 = 0; % 初始横摇速度 q0 = 0; % 初始纵摇速度 r0 = 0; % 初始航向角速度 phi0 = 0; % 初始横摇角 theta0 = 0; % 初始纵摇角 psi0 = 0; % 初始航向角 % 定义时间间隔和模拟时间 dt = 0.1; % 时间间隔 t_end = 100; % 模拟时间 % 初始化状态变量 u = u0; v = v0; w = w0; p = p0; q = q0; r = r0; phi = phi0; theta = theta0; psi = psi0; % 定义控制量 delta = 0; % 舵角 thrust = 0; % 推力 % 循环计算状态变量 for t=0:dt:t_end % 计算水动力力矩 X = 0; % 船舶横向力 Y = 0; % 船舶纵向力 Z = -rho*g*T*B*L; % 船舶垂向力 K = 0; % 船舶横摇力矩 M = 0; % 船舶纵摇力矩 N = 0; % 船舶航向力矩 % 计算控制力矩 K_delta = 0; % 舵面横摇力矩系数 N_delta = 0; % 舵面航向力矩系数 K_thr = 0; % 推力横摇力矩系数 N_thr = 0; % 推力航向力矩系数 K_control = K_delta*delta + K_thr*thrust; N_control = N_delta*delta + N_thr*thrust; % 计算运动方程 u_dot = r*v-q*w-X/m; v_dot = p*w-r*u+Y/m; w_dot = q*u-p*v+Z/m; p_dot = (K-K_control)/Iz; q_dot = (M)/Iz; r_dot = (N-N_control)/Iz; phi_dot = p + tan(theta)*(q*sin(phi)+r*cos(phi)); theta_dot = q*cos(phi)-r*sin(phi); psi_dot = (q*sin(phi)+r*cos(phi))/cos(theta); % 更新状态变量 u = u + u_dot*dt; v = v + v_dot*dt; w = w + w_dot*dt; p = p + p_dot*dt; q = q + q_dot*dt; r = r + r_dot*dt; phi = phi + phi_dot*dt; theta = theta + theta_dot*dt; psi = psi + psi_dot*dt; end ``` 这只是一个简单的示例，实际情况可能更加复杂。具体实现要根据船舶的具体情况和模型来确定。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

go_bananas

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值