一种船舶敞水航行辅助决策方法技术

本发明专利技术属于船舶结构安全性领域，具体涉及一种船舶敞水航行辅助决策方法。包括建立数据库、拟合图像、对图像进行插值、将图像拟合出的曲面离散成若干点，求出船舶在此航速，浪向角下对应的发生警报的危险剖面的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值分别为My0,Mz0,Mx0,Fz0；在数据库中筛选出满足My1

【技术实现步骤摘要】
一种船舶敞水航行辅助决策方法
本专利技术属于船舶结构安全性领域，具体涉及一种船舶敞水航行辅助决策方法。
技术介绍
结构的安全性是保证船舶正常营运的基础，在船舶的结构设计中，设计人员需要参考有关船级社的结构设计规范或准则，结合相关的理论方法，确定船舶的各种设计载荷，进而进行有关的结构设计与强度评估。但是，对于实际海洋环境中航行的船舶来说，结构受到的外载荷具有很强的随机性，这些随机因素很难通过规范或准则中的方法进行完全准确地预报。因而，尽管船级社已经给出了船舶结构设计的规范或准则，设计人员也据此给出尽可能完善的结构设计方案，但船舶仍存在着一定程度结构失效的可能性。为了降低船舶结构安全的风险，一方面可以在设计建造阶段合理确定结构设防水平，不断完善船体结构强度分析与评估方法，另一方面是对船舶使用过程加强管理和监控。通过实时监测船体结构应力，并根据得到的监测数据对结构状态进行评估，当结构应力较大时发出报警，并给操船人员提供可供选择的辅助航行决策建议，则可以有效降低结构的安全风险。以往在船舶敞水航行时并没有明确的辅助决策方法，现根据三维频域线性水动力理论计算出危险剖面剖面载荷的百一值，并根据筛选百一值确定船舶航速及浪向角的变化情况，以此实现船舶敞水航行的辅助决策的方法。
技术实现思路
本专利技术的目的在于提供一种船舶敞水航行辅助决策方法，其能在船舶结构应力较大时发出报警，并给操船人员提供可供选择的辅助航行决策建议。一种船舶敞水航行辅助决策方法，具体包括如下步骤：步骤1、建立数据库；步骤2、拟合图像；步骤3、输入此时海况对应的参数值：海浪有义波高Hs和平均跨零周期Ts；步骤4、对图像进行插值，对拟合出的不同海浪有义波高，平均跨零周期下的三维图像进行插值，得出对应此波高，平均跨零周期下的三维图像，即航速V，浪向角β，同时得出垂向弯矩My，水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图及关系方程；步骤5、将图像拟合出的曲面离散成若干点，得到细化后的航速V，浪向角β对应的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图；步骤6、输入此时船的航速V0及浪向角β0；步骤7、求出船舶在此航速，浪向角下对应的发生警报的危险剖面的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值分别为My0,Mz0,Mx0,Fz0；步骤8、在数据库中筛选出满足My1<My0此条件所对应的(V,β)的集合{A0}，其次再从集合{A0}中筛选出满足Mz1<Mz0此条件所对应的(V,β)的集合{A1}，然后再从集合{A1}中筛选出满足Mx1<Mx0此条件所对应的(V,β)的集合{A2}，最后再从集合{A2}中筛选出满足Fz1<Fz0此条件所对应的(V,β)的集合{A3}；步骤9、通过筛选集合{A3}选出与所给航速V0及浪向角β0最接近的新的航速V1及浪向角β1，即使的值最小；步骤10、确定新的航速V1及浪向角β1后，在此有义波高，平均跨零周期下，求出危险剖面的载荷包括垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz，判断危险剖面上极值应力值是否小于许用应力值，若小于，结束，若极值应力不小于许用应力值，则继续从集合{A3}中筛选航速V1及浪向角β1，直到危险剖面上极值应力应小于许用应力值并且的值最小。所述一种船舶敞水航行辅助决策方法，步骤1具体包括基于三维频域线性水动力理论对船舶进行短期分析，计算出不同海况不同航速不同浪向角条件下的各个危险剖面垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值。所述一种船舶敞水航行辅助决策方法，步骤2具体包括以航速为X轴，浪向角为Y轴，各个危险剖面剖面载荷的百一值为Z轴，进行三维绘图及拟合图像，即描绘出某一确定海况下，某个危险剖面剖面载荷百一值与航速，浪向角的三维关系，拟合图像并确定出图像方程。本专利技术的有益效果在于：本专利技术根据三维频域线性水动力理论计算出危险剖面剖面载荷的百一值，并根据筛选百一值确定船舶航速及浪向角的变化情况，在已知当时海况波浪的波高及平均跨零周期，船舶此时航行的航速，浪向角的条件下，迅速找出满足小于此海况下的各个百一值所对应的新的航速及浪向角，通过筛选选出与所给航速及浪向角最接近的新的航速及浪向角。本专利技术在船体结构应力监测点发生报警时，能根据当时的海况、航向、航速，对船体载荷进行总体计算分析与评估，并快速作出如何改变航向、改变航速以及船舶姿态等操作指令。附图说明图1为监测对象某工程船有限元模型；图2为有义波高Hs＝4m，平均跨零周期Ts＝9.5s垂向弯矩My的百一值；图3为有义波高Hs＝4m，平均跨零周期Ts＝13.5s垂向弯矩My的百一值；图4为有义波高Hs＝6m，平均跨零周期Ts＝9.5s垂向弯矩My的百一值；图5为有义波高Hs＝6m，平均跨零周期Ts＝13.5s垂向弯矩My的百一值；图6为有义波高Hs＝5m，平均跨零周期Ts＝11.5s垂向弯矩My的百一值；图7为本专利技术的总过程流程图；图8为有义波高Hs＝5m，平均跨零周期Ts＝11.5s，航速V＝7kn，浪向角β＝130°各计算工况设计波参数；图9为工况1极值应力计算结果；图10为工况2极值应力计算结果；图11为工况3极值应力计算结果；图12为工况4极值应力计算结果；图13为工况5极值应力计算结果；图14为工况6极值应力计算结果；图15为工况7极值应力计算结果；图16为工况8极值应力计算结果。具体实施方式下面结合附图对本专利技术做进一步说明。如附图1所示，为监测对象某工程船有限元模型；主尺度为总长：122.5m；型宽：22.32m；型深：11.80m；设计吃水：7.85m，以有义波高Hs＝5m，平均跨零周期Ts＝11.5s，航速V＝10m/s，浪向角130°船中剖面发生报警情况为例说明本方法，附图7为本专利技术的总过程流程图，具体包括如下步骤。1、建立数据库，即基于三维频域线性水动力理论对某工程船进行规则波响应计算及短期分析，计算出不同海况航速为0，3kn，6kn，9kn，12kn，15kn，18kn，浪向角为0°,15°,30°,45°,60°,75°,90°,105°,120°,135°,150°,165°,180°条件下的各个危险剖面垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值。2、以航速为X轴，浪向角为Y轴，各个危险剖面剖面载荷的百一值为Z轴，进行三维绘图及拟合图像，通过筛选，得到需要的危险剖面三维图像，现已获知海浪有义波高Hs＝5m，平均跨零周期Ts＝11.5s，根据拟合出的不同海浪有义波高，平均跨零周期下的三维图像进行插值，得出对应此波高，平均跨零周期下的三维图像，即航速V，浪向角β与垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图及关系方程。以绘制垂向弯矩My的百一值的新图像为例，即有义波高Hs＝4m，平均跨零周期Ts＝9.5s的图2，有义波高Hs＝4m，平均跨零周期Ts＝13.5s的图3，有义波高Hs＝6m，平均跨零周期Ts＝9.5s的图4，有义波高Hs＝6m，平均跨零周期Ts＝13.5s的图5，再通过此四幅图像，进行插值，得到新图像，如图6。同理可绘制水平弯矩Mz，扭矩Mx，垂本文档来自技高网...

【技术保护点】
1.一种船舶敞水航行辅助决策方法，其特征在于，具体包括如下步骤：步骤1、建立数据库；步骤2、拟合图像；步骤3、输入此时海况对应的参数值：海浪有义波高Hs和平均跨零周期Ts；步骤4、对图像进行插值，对拟合出的不同海浪有义波高，平均跨零周期下的三维图像进行插值，得出对应此波高，平均跨零周期下的三维图像，即航速V，浪向角β，同时得出垂向弯矩My，水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图及关系方程；步骤5、将图像拟合出的曲面离散成若干点，得到细化后的航速V，浪向角β对应的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图；步骤6、输入此时船的航速V0及浪向角β0；步骤7、求出船舶在此航速，浪向角下对应的发生警报的危险剖面的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值My0,Mz0,Mx0和Fz0；步骤8、在数据库中筛选出满足My1

【技术特征摘要】
1.一种船舶敞水航行辅助决策方法，其特征在于，具体包括如下步骤：步骤1、建立数据库；步骤2、拟合图像；步骤3、输入此时海况对应的参数值：海浪有义波高Hs和平均跨零周期Ts；步骤4、对图像进行插值，对拟合出的不同海浪有义波高，平均跨零周期下的三维图像进行插值，得出对应此波高，平均跨零周期下的三维图像，即航速V，浪向角β，同时得出垂向弯矩My，水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图及关系方程；步骤5、将图像拟合出的曲面离散成若干点，得到细化后的航速V，浪向角β对应的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值在此海况下的关系图；步骤6、输入此时船的航速V0及浪向角β0；步骤7、求出船舶在此航速，浪向角下对应的发生警报的危险剖面的垂向弯矩My,水平弯矩Mz,扭矩Mx,垂向剪力Fz的百一值My0,Mz0,Mx0和Fz0；步骤8、在数据库中筛选出满足My1<My0此条件所对应的(V,β)的集合{A0}，其次再从集合{A0}中筛选出满足Mz1<Mz0此条件所对应的(V,β)的集合{A1}，然后再从集合{A1}中筛选出满足Mx1<Mx0此条件所对应的(V,β...

【专利技术属性】
技术研发人员：刘剑飞，李磊，冯国庆，冯玮，康有为，邓燕飞，王波，田宇，呼文佳，陈思均，
申请(专利权)人：哈尔滨工程大学，中集海洋工程有限公司，
类型：发明
国别省市：黑龙江,23