MATLAB运用——设计船舶模型

已于 2022-04-12 00:07:41 修改
阅读量4.8k 收藏 56
点赞数 13
分类专栏： MATLAB学习 文章标签： matlab 材质
于 2022-04-12 00:05:37 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/go_bananas/article/details/124112440
版权

这次我们的目标如下：
在这里插入图片描述
可以看到项目里面最大的挑战是控制船的倾翻角度。

设计步骤：

在这里插入图片描述
下面我们将按照这一设计步骤一步一步地来完成。

计算重心，计算空载曲线

把这两步合并在一起的原因是：船的重心其实就等于空载曲线下面排开水的部分的质心，所以当我们算出船的空载曲线时，船的重心也就出来了。

可以看我的这篇文章来了解浮心的算法。
我们计划的船型是一个100X30X20(cm)的船。
为了方便计算与拼装，我们不妨设船舶的每一片切面的形状都是相同的，只是大小发生变化。
于是我们设计出船舶的切片的标准形状：
在这里插入图片描述
你可以认为船的每一片切面都是由这个“标准切面”等比例缩放而来的
而这个等比例缩放在x轴满足下面的分布：
在这里插入图片描述

有了这几点前提我们的空载曲线就很好算了：
首先定义出水平面：

%这条是空载水线z=-4
    d = -4;
    theta = 0;
    y = mesh.ys;
    z = tand(theta) .* y + d;%每个水线对应的方程

为什么空载曲线是z=-4?
其实是我们利用二分的思想一次次试出来的，
我们计划船舶的重量是5800g左右,排开水的重量随排水线的高度单调增加，所以我们二分排水线高度，当排开水的体积等一重量时，这条水线就是我们的空载水线。

下面是计算浮心和重量的代码：

    %这条是空载水线z=-4
    d = -4;
    theta = 0;
    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);%绘制图像
    mass=matrixSum(sub_region)
    COB = centerOfMass2(sub_region,mesh); 
    V=0;
    %计算浮心2D
    for i=1:30 %X从0到100切30次每次长度为3
        Zjudge=Judge(i*3)
        bili=(10-Zjudge)/20.0%相关的缩放比例
        if i>=5 && i<=25
            bili=1;
        end
        yyy=COB(1);%标准面的重心的y坐标
        zzz=COB(2);%标准面的重心的x坐标
        ZMass(end+1)=Zjudge+(zzz+10)*bili
        YMass(end+1)=yyy*bili
        XMass(end+1)=(i*3)
        MMass(end+1)=mass*bili
        tot_mass=tot_mass+mass*bili;%用来算总质心的权重比
        V=V+mass*3*bili;%总体积
    end
   ANS=centerOfMass3(XMass,YMass,ZMass,MMass,tot_mass);%计算三维浮心
   M=V*1%总重量

在这里插入图片描述

算出来在这条水线下的排开水的体积等于4275g(还会加上电机电池与两个载重)与我们船的质量5800g相近，我们就得出了船的空载水线z=-4；
但是算质心时，我们应该再加上电机和电池的重心和载重物品的重心：
于是：

算出裸船的重心在（48，0，-5.6）处。
再加上两个电机一个电池(800g),重心在(65,0,-6)
   ANS(1)=ANS(1,1)*M/(M+800)+65*800/(M+800)
   ANs(2)=ANS(1,2)*M/(M+800)+0*800/(M+800)
   Ans(3)=ANS(1,3)*M/(M+800)+(-6)*800/(M+800)
   ANS
   M=M+800;
   
再加上载重（800g）,重心再在（50，0，6）
   ANS(1)=ANS(1)*M/(M+800)+50*800/(M+800)
   ANs(2)=ANS(2)*M/(M+800)+0*800/(M+800)
   Ans(3)=ANS(3)*M/(M+800)+(6)*800/(M+800)
   ANS
   M=M+800
得出船在正常情况载重为800g时的重心：
（48.4，0，-5.6）、
船重：5800g

在这里插入图片描述

计算130度水线的浮心

我们只要把上面的水线的代码稍作改动，使得它的角度为130度排开水的体积是船的重量即可，经一番测试，我们得出水线为：

    d = 9;
    theta = 50;%取>的部分故反过来为50度
    y = mesh.ys;
    z = tand(theta) .* y + d;%每个水线对应的方程

时，截取的船截面为：
在这里插入图片描述
计算这种水线下的浮心为：
在这里插入图片描述
于是我们得出了130度时的浮心位置。

计算扶正力矩：

给出计算图：
在这里插入图片描述
会发现在130度的情况下船舶是具有扶正能力的，且在140度时不具备扶正能力。

后续

接下来就是按照这一模型进行soildwork建模，造出实物来验算了。

完整版代码如下：

计算浮心(单位：厘米)：
clear all
boat.L = 100;          %船的长度(厘米)
boat.W = 30;          %船的宽度
boat.HB = boat.W / 2;  %船的半宽
boat.D = 20;           %船的深度
boat.HD = boat.D / 2;  %船的半宽

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

mesh.ys = -1.0*(-boat.HB:dy:boat.HB);        % meters  Y轴切片
mesh.zs = -1.0*(-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所对应的点切片：
    ZZ=[];
    for i=0 : 30 %讨论y从[-10,10]
        t=caculate(i);
        ZZ=[ZZ t];
    end
    y = mesh.ys;
    figure(1);%绘制不同的图
    %axis xy;
    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
    %这条是空载水线z=-4
    d = -4;
    theta = 0;
    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);%绘制图像
    mass=matrixSum(sub_region)
    COB = centerOfMass2(sub_region,mesh); 
    V=0;
    %计算浮心2D
    for i=1:30 %X从0到100切30次每次长度为3
        Zjudge=Judge(i*3)
        bili=(10-Zjudge)/20.0%相关的缩放比例
        if i>=5 && i<=25
            bili=1;
        end
        yyy=COB(1);%标准面的重心的y坐标
        zzz=COB(2);%标准面的重心的x坐标
        ZMass(end+1)=Zjudge+(zzz+10)*bili
        YMass(end+1)=yyy*bili
        XMass(end+1)=(i*3)
        MMass(end+1)=mass*bili
        tot_mass=tot_mass+mass*bili;%用来算总质心的权重比
        V=V+mass*3*bili;%总体积
    end
   ANS=centerOfMass3(XMass,YMass,ZMass,MMass,tot_mass);%计算三维浮心
   M=V*1%总重量
算出裸船的重心在（48，0，-5.6）处。
再加上两个电机一个电池(800g),重心在(65,0,-6)
   ANS(1)=ANS(1,1)*M/(M+800)+65*800/(M+800)
   ANs(2)=ANS(1,2)*M/(M+800)+0*800/(M+800)
   Ans(3)=ANS(1,3)*M/(M+800)+(-6)*800/(M+800)
   ANS
   M=M+800;
   
再加上载重（800g）,重心再在（50，0，6）
   ANS(1)=ANS(1)*M/(M+800)+50*800/(M+800)
   ANs(2)=ANS(2)*M/(M+800)+0*800/(M+800)
   Ans(3)=ANS(3)*M/(M+800)+(6)*800/(M+800)
   ANS
   M=M+800
得出船在正常情况载重为800g时的重心：
（48.4，0，-5.6）、
船重：5800g

function z=Judge(zj)
    if zj<15
       z=-4.0/3*zj+10; 
    end
    if zj>75
        z=4.0/3*(zj-75)-10;
    end
    if zj>=15 && zj<=75
        z=-10;
    end
end

function ANS=centerOfMass3(XMass,YMass,ZMass,MMass,m)
x=0;
y=0;
z=0;
    for i=1:30
        x=x+XMass(1,i)*MMass(1,i)/m;
        y=y+YMass(1,i)*MMass(1,i)/m;
        z=z+ZMass(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 z=caculate(a)
    a=a-15;
    if a<=-11 
       z=-5.0/2.0*a-55.0/2;
    end
    if a>=11
       z=5.0/2.0*a-55.0/2;
    end
    if a>(-11) && a<11
      z=10.0/121*(a-11)*(a+11);
    end
end
go_bananas
关注
  • 13
    点赞
  • 踩
  • 56
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
无人船 无人艇路径跟踪控制 Matlab Simulink效果
QFkMbEA的博客
01-27 750
通过利用Matlab和Simulink的强大功能，我们可以基于Foseen船舶模型，采用LOS制导子系统和反步控制器实现路径跟踪控制，并利用ESO实现漂角和侧滑角的估计。在无人船和无人艇的控制过程中，反步控制器可以用于实现航向和速度的控制，使得船艇能够按照预定的路径稳定地行驶。总之，路径跟踪控制技术的研究和应用将进一步推动无人船和无人艇技术的发展，并为海洋科研、水下勘察和环境监测等领域提供更加可靠、高效的解决方案。相信随着技术的不断发展，无人船和无人艇的路径跟踪控制将取得更加重要的进展和应用。
船运动程序matlab_船舶_；船舶航向_PID控制_
10-02
MATLAB，基于PID船舶航向控制程序，船舶
船舶海上运动matlab仿真
12-20
根据平稳随机过程理论，根据船舶海上运动频谱，模拟了船舶的海上横摇、纵摇运动
基于matlab的船体三维模型仿真
weixin_47108207的博客
10-04 1954
船体数据来源来自于C++代码。在C++中，将船体数据导入到EXCEL中，然后再MATLAB中调用EXCEL的船体数据，进行三维建模。在建模过程中，因为船体数据不可能每个点的距离都是无限小，所以我们采用了拟合的办法，调用了MATLB的函数 griddata 。
使用模型预测控制对USV进行自主控制（Matlab代码实现）
m0_73907476的博客
04-08 1477
本文介绍了一种模型预控制（MPC）算法，旨在自动驾驶无人水面车辆（USV）驶向一组航路点。USV和扰动的建模已经简化，因为这项工作旨在证明概念：对于现实世界的实施，应该考虑更准确的建模。无人船(unmanned surface vehicles,USV)是一种船端无人操控的水面船舶,近年来受到了广泛关注。如何实现自主航行是USV面临的核心问题,而设计一种具有精确航迹控制能力的运动控制器是解决该问题的基础。[1]柳晨光. 基于预测控制的无人船运动控制方法研究[D]. 武汉理工大学.
【路径规划】基于RRT算法实现实际地图船舶航路避障规划附matlab代码
Matlab_dashi的博客
06-04 267
✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。????个人主页：Matlab科研工作室????个人信条：格物致知。⛄ 内容介绍一种基于RRT算法的无人船海面路径规划方法及系统,包括获取海图数据,得到海图整体信息;按照海图多边形的定义呈现海面情况,并将原始的矢量数据进行栅格化得到完整的海图...
船舶强度与结构设计大作业二matlab,船舶强度与结构设计题目
weixin_36328628的博客
03-19 1263
1、船底是船体等值梁的下缘，下列各选项载荷中属于船底板架受到的载荷是( ) A：总纵弯曲应力 B：货物横向载荷 C：首部冲击力 D：以上三项都是 参考答案：D 为了减小舱口角隅处应力集中，除了要求去除方角而采用圆弧形角隅，并在角隅处采用加复板或厚板进行加强，同时选用的钢材要采用 级钢； A：级或E/EH级钢 B：级或E/EH级钢 C：级或D/DH级钢 D：级或D/DH级钢 参考答案：A 下列各项重...
【MATLAB代码】基于模拟退化算法的船舶调度优化问题、船舶进港排班
运筹不帷幄的博客
03-18 5010
目录一、模型建立二、模型数据三、模拟退火算法编写四、运行结果 一、模型建立 符号定义： g——港口位置 ggg ∈\in∈ G,1,2,…nG,{1,2,…n}G,1,2,…n iii——泊位位置iii ∈\in∈ I,1,2,…nI, {1,2,…n}I,1,2,…n jjj——船舶位置jjj ∈\in∈ J,1,2,…nJ,{1,2,…n}J,1,2,…n TwjTwjTwj——船j的等待时间 TkjTkjTkj——船j的卸货时间（一个集装箱的卸货时间是15分钟） TajTajTaj——船j的到达时间
MATLAB编程——灰色预测模型源代码.rar
01-08
通过对原始数据序列X0进行准光滑性检验，得出：预测值，残差数列，后验差检验，方差比 ，计算小误差概率等，非常适合选数据处理题的数学建模同学，压缩包里几乎涵盖了所有灰色预测模型的代码。
船舶一阶野本（Nomoto）模型
09-04
matlab源代码，非simulink，野本模型，Z形试验，参数由辨识模型所得
SailingYachtModel.z​ip:在 Matlab/Simulink 中使用 3D 动画模拟帆船游艇模型。-matlab开发
06-01
使用 Simulink 3D 动画工具箱对具有 3D 动画的帆船（包括翻滚）进行 4DOF 模拟。 模型、符号和参数取自 Xiao 和 Jouffroy，“帆船的建模和非线性航向控制”，IEEE 海洋工程杂志，卷。 39，没有。 2、2014 年。 在 Matlab r2013a、Simulink 8.1、Simulink 3D Animation 7.0 中测试。
船舶运动一键建模软件
09-10
只要输入船舶具体参数就可以直接得到船舶模型
船舶运动数学模型
04-10
杨盐生的船舶运动数学模型，适用于船舶建模，关键词为nomoto,船舶建模，MMG
船舶运动控制Matlab
05-11
船舶运动控制matlab模拟仿真系统编程
Matlab模拟船舶操纵运动模型
06-29
1.版本：matlab2019a，不会运行可私信 2.领域：基础教程 3.内容：Matlab模拟船舶操纵运动模型 4.适合人群：本科，硕士等教研学习使用
matlab GUI课程设计——刮刮乐
05-25
平台：matlab appdesigner 功能：模拟刮刮乐活动。具备鼠标刮奖、增删奖项功能 文章链接：https://blog.csdn.net/gouqikong5855/article/details/124972392
matlab GUI课程设计——简易计算器
05-25
平台：matlab GUIDE 功能：演示了基本的四则运算操作。 文章链接：https://blog.csdn.net/gouqikong5855/article/details/124970542
matlab数据 —— sd.xlsx
06-27
matlab数据 —— sd.xlsx
matlab神经网络设计电机模型
最新发布
07-28
在使用MATLAB实现神经网络来设计电机模型时，首先需要将电机的相关特征和参数进行数据采集和整理，并作为输入数据。这些特征和参数可以包括电机的电压、电流、转速、负载等信息。 接下来，在MATLAB中选择适当的神经网络模型，如前馈神经网络或循环神经网络。然后，将之前采集的电机数据作为训练集，通过神经网络的训练算法对神经网络进行训练。这可以使用MATLAB中自带的神经网络工具箱来完成。 在训练过程中，首先需要对数据进行预处理，如归一化、标准化等操作，以便提高训练的效果。然后，使用训练函数对神经网络进行训练，不断调整神经网络的权重和阈值，使得模型能够尽可能准确地预测电机的行为。 训练完成后，可以使用已经训练好的神经网络来进行电机模型的预测。这可以通过将新的电机输入数据输入到神经网络中，并获取神经网络的输出结果，即电机的预测行为。 最后，需要对神经网络模型进行评估和优化。可以比较神经网络的预测结果与实际电机行为的差异，以评估神经网络模型的性能。如果需要，可以通过调整神经网络的结构和参数，如增加隐藏层、调整神经元个数等操作，进一步优化电机模型的预测精度。 综上所述，使用MATLAB进行神经网络设计电机模型，需要数据采集、网络选择、训练、预测和优化等步骤。通过这些步骤，能够帮助我们更好地理解和预测电机的行为，提高电机的性能与效能。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

go_bananas

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

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

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

打赏作者

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

抵扣说明：

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

余额充值