CN111591407A - 一种船体型线特征参数提取方法 - Google Patents

一种船体型线特征参数提取方法 Download PDF

Info

Publication number
CN111591407A
CN111591407A CN202010489918.2A CN202010489918A CN111591407A CN 111591407 A CN111591407 A CN 111591407A CN 202010489918 A CN202010489918 A CN 202010489918A CN 111591407 A CN111591407 A CN 111591407A
Authority
CN
China
Prior art keywords
discrete points
ship
line
area
parabola
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010489918.2A
Other languages
English (en)
Other versions
CN111591407B (zh
Inventor
沈冠之
陈京普
冯君
孙文愈
刘希洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
702th Research Institute of CSIC
Original Assignee
702th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 702th Research Institute of CSIC filed Critical 702th Research Institute of CSIC
Priority to CN202010489918.2A priority Critical patent/CN111591407B/zh
Publication of CN111591407A publication Critical patent/CN111591407A/zh
Application granted granted Critical
Publication of CN111591407B publication Critical patent/CN111591407B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B71/00Designing vessels; Predicting their performance

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种船体型线特征参数提取方法,涉及船舶技术领域,该方法通过对船体的船型文件进行处理,可以自动化提取得到船体型线在给定吃水深度下的特征参数,在进行特征提取时,根据各特征参数的参数类型采用不同方法提取,对尺度、位置等代数型参数直接提取,对面积、体积等积分型参数进行积分求解处理提取,对角度、曲率等微分型参数进行微分求解处理提取,从而可以提取得到船体的表征整体特征的全局参数、表征区域特征的局部参数以及表征曲面细节特征的细部参数三个层次的特征参数,为后续的数据挖掘提供了数据基础。

Description

一种船体型线特征参数提取方法
技术领域
本发明涉及船舶技术领域,尤其是一种船体型线特征参数提取方法。
背景技术
船体型线智能化设计过程中,需要从船体型线特征、流场特征、性能指标等大量的数据中提取有价值的信息,并开展下一步的结果分析和数据挖掘等工作。因此,需要对表达船体型线的特征参数进行合理的分析和提取,但船体型线特征复杂,目前较难开展深层的数据挖掘。
发明内容
本发明人针对上述问题及技术需求,提出了一种船体型线特征参数提取方法,该方法包括如下步骤:
获取船体的船型文件并获取吃水深度,船型文件中包括表征船体的结构的型线;
直接基于船型文件的数据提取船体的尺度参数;
对船型文件的数据进行积分求解处理提取船体的剖面面积和排水体积;
对船型文件的数据进行微分求解处理提取船体在吃水深度下的横剖面船型特征值,横剖面船型特征值用于表征横剖面UV度;
提取得到船体型线的特征参数,包括尺度参数、位置参数、剖面面积、排水体积以及横剖面船型特征值。
本发明的有益技术效果是:
本申请公开了一种船体型线特征参数提取方法,该方法通过对船体的船型文件进行处理,可以自动化提取得到船体型线在给定吃水深度下的特征参数,在进行特征提取时,根据各特征参数的参数类型采用不同方法提取,对尺度、位置等代数型参数直接提取,对面积、体积等积分型参数进行积分求解处理提取,对角度、曲率等微分型参数进行微分求解处理提取,从而可以提取得到船体的表征整体特征的全局参数、表征区域特征的局部参数以及表征曲面细节特征的细部参数三个层次的特征参数,为后续的数据挖掘提供了数据基础。
附图说明
图1是本申请的船体型线特征参数提取方法的方法流程图。
图2是船体的纵剖面线的示意图。
图3是船体的横剖面线的一种局部示意图。
图4是对曲线进行拟合抛物线时的示意图。
图5是船体的横剖面线的另一种局部示意图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。
本申请公开了一种船体型线特征参数提取方法,该方法包括如下步骤,请结合图1所示的流程图:
步骤S1,获取船体的船型文件并获取吃水深度T,船型文件中包括表征船体的结构的型线,常见的包括横剖面线、纵剖面线、设计水线等等。
步骤S2,直接基于船型文件的数据提取船体的尺度参数。
步骤S3,对船型文件的数据进行积分求解处理提取船体的剖面面积和排水体积。
步骤S4,对船型文件的数据进行微分求解处理提取船体在当前的吃水深度下的横剖面船型特征值,横剖面船型特征值用于表征横剖面UV度。
步骤S5,提取得到船体型线的特征参数,包括尺度参数、位置参数、剖面面积、排水体积以及横剖面船型特征值。
一、在上述步骤S2中,直接基于船型文件的数据提取到的主要包括一些代数型参数,具体的:
(1)提取得到的代数型参数主要是如上所述的船体结构的尺度参数,包括船体总长、垂线间长、水线长以及型宽,请结合图2所示的纵剖面线示意图,提取方法为,:
确定船体的船艏最前端至船艉最后端的水平距离为船体总长Loa。
确定船体的艏垂线与艉垂线之间的水平距离为垂线间长Lpp,艏垂线是通过设计水线与艏柱前缘的交点所作的垂直于设计水线面的垂线,艉垂线位于0站横剖面上。
确定设计水线最前端至船艉最后端的水平距离为水线长Lwl。
确定船体两侧型表面之间垂直于中线面的水平距离为型宽B,图中未示出,计算时采用中横剖面设计水线处的宽度。
(2)对于某些船型,这些船型的横剖面线或部分横剖面线具有较大的曲率,请结合图3,则提取得到的代数型参数还包括横剖面线的反曲点,比如常见的,一些船型的艉部的横剖面线具有较大的曲率。在提取每条横剖面线的反曲点时,采用逐步逼近法确定最逼近于横剖面线的逼近直线,如图3中虚线所示,遍历横剖面线上的所有点找到与逼近直线之间直线距离最远的点作为横剖面线的反曲点。
二、在上述步骤S3中,对船型文件的数据进行积分求解处理提取到的主要包括一些积分型参数,主要为面积参数和体积参数,常见的即为上述剖面面积和排水体积,在求解时,用抛物线拟合船体型线并计算抛物线包含面积得到剖面面积,然后进一步用抛物线拟合剖面面积并计算抛物线包含面积得到排水体积。
对于船体型线和剖面面积中的任意一种曲线,假设构成曲线的所有离散点共形成n个包络,每相邻两个离散点之间形成一个包络,若n=1,则直接使用梯形计算公式计算,但通常情况下这种情况较少,其余情况用抛物线拟合曲线并计算抛物线包含面积的方法为,请结合图4:
按照三个离散点为一组将离散点划分为若干组离散点,每组离散点中的三个离散点形成两个包络,n≥2。比如图4中,第一组离散点形成包括包络1和包括2,第二组离散点包括包络3和包络4,以此类推,相邻两组离散点会共用一个离散点。
基于每组离散点中的三个离散点拟合得到对应的子抛物线,子抛物线方程表示为y=ax2+bx+c。假设该组离散点中三个离散点的坐标分别为(x0,y0)、(x1,y1)和(x2,y2),则求解该子抛物线有:
Figure BDA0002520678550000031
Figure BDA0002520678550000032
Figure BDA0002520678550000033
引入变量令h1=x1-x0,h2=x2-x1,l=h1+h2=x2-x0,则求解上述公式可得子抛物线方程的系数为:
Figure BDA0002520678550000034
对于每组离散点得到的上述子抛物线,求解得到子抛物线的面积,子抛物线的面积包括该组离散点形成的两个包络的总面积,具体的:
Figure BDA0002520678550000041
时,同时计算三个离散点包含的面积得到该组离散点对应的子抛物线的面积:
Figure BDA0002520678550000042
Figure BDA0002520678550000043
或者
Figure BDA0002520678550000044
时,分别计算该组离散点形成的两个包络的面积两个包络的面积的求解公式为:
Figure BDA0002520678550000045
Figure BDA0002520678550000046
求和得到该组离散点对应的子抛物线的面积S=S1+S2
由此可以求解得到每组离散点对应的子抛物线的面积,若n为偶数,则在划分离散点时正好可以划分为n/2组,此时直接对每组离散点对应的子抛物线的面积进行求和得到抛物线包含面积,也即相应得到剖面面积或排水体积。
若n为奇数,则此时会存在最后一个离散点未被划分到任何一组离散点内,若仅将每组离散点对应的子抛物线的面积进行求和则最后一个包络的面积未计算到,也即图中的包络n的面积未计算到,此时将该最后一个离散点与前面的数据进行组合,也即未划分到任何一组离散点内的最后一个离散点及其最接近的前两个离散点拟合得到对应的子抛物线,假设以(x2,y2)表示未划分到任何一组离散点内的最后一个离散点的坐标、(x1,y1)表示(x2,y2)最接近的前一个离散点的坐标、(x0,y0)表示(x2,y2)次接近的前一个离散点的坐标时,则得到的子抛物线的计算公式同样如上所述,则利用
Figure BDA0002520678550000047
即能基于该子抛物线求解得到最后一个离散点及其最接近的一个离散点之间形成的包络的落单包络面积,也即图4中包络n的面积,对求解得到的落单包络面积以及其余每组离散点对应的子抛物线的面积进行求和得到抛物线包含面积,也即相应得到剖面面积或排水体积。
在上述计算过程中,进一步的,当n=2且
Figure BDA0002520678550000051
时,还可以直接利用两个梯形的面积之和求解得到子抛物线的面积也即
Figure BDA0002520678550000052
三、在上述步骤S4中,对船型文件的数据进行微分求解处理提取得到的主要包括一些微分型参数,主要是一些角度、曲率的参数,常见的即为上述提到的每条横剖面线的横剖面船型特征值,具体的,请结合图5,对于每条横剖面线:
取横剖面线在0.2T~0.8T范围内的型线段,其中T表示吃水深度。求得最逼近于该型线段的直线的直线斜率代表该剖面的特征,求解方式为:
假设最逼近于该型线段的直线的方程式为y=m·z+b,利用最小二乘方法求解直线斜率m的方程式为:
Figure BDA0002520678550000053
解得:
Figure BDA0002520678550000054
式中,
Figure BDA0002520678550000055
j=1,2,3,4,j表示水线序号,yj为zj处的型值,代入上式进一步可得:
Figure BDA0002520678550000056
因为m尚受B/T影响,因此确定横剖面线的横剖面船型特征值为:
Figure BDA0002520678550000057
其中B表示船体的型宽。
上述得到的横剖面船型特征值为按绝对型值计算得到的,进一步的,本申请对型值进行无量纲化得到无量纲参数
Figure BDA0002520678550000058
Figure BDA0002520678550000059
则基于无量纲参数将按绝对型值计算的横剖面船型特征值处理为按无量纲型值计算的横剖面船型特征值表示为
Figure BDA0002520678550000061
以上所述的仅是本申请的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。

Claims (10)

1.一种船体型线特征参数提取方法,其特征在于,所述方法包括:
获取船体的船型文件并获取吃水深度,所述船型文件中包括表征所述船体的结构的型线;
直接基于所述船型文件的数据提取所述船体的尺度参数;
对所述船型文件的数据进行积分求解处理提取所述船体的剖面面积和排水体积;
对所述船型文件的数据进行微分求解处理提取所述船体在所述吃水深度下的横剖面船型特征值,所述横剖面船型特征值用于表征横剖面UV度;
提取得到船体型线的特征参数,包括所述尺度参数、位置参数、剖面面积、排水体积以及横剖面船型特征值。
2.根据权利要求1所述的方法,其特征在于,直接基于所述船型文件的数据提取得到的所述尺度参数包括船体总长、垂线间长、水线长以及型宽。
3.根据权利要求2所述的方法,其特征在于,所述直接基于所述船型文件的数据提取所述船体的尺度参数,包括:
确定所述船体的船艏最前端至船艉最后端的水平距离为所述船体总长;
确定所述船体的艏垂线与艉垂线之间的水平距离为所述垂线间长,所述艏垂线是通过设计水线与艏柱前缘的交点所作的垂直于设计水线面的垂线,所述艉垂线位于0站横剖面上;
确定设计水线最前端至船艉最后端的水平距离为所述水线长;
确定所述船体两侧型表面之间垂直于中线面的水平距离为所述型宽。
4.根据权利要求1-3任一所述的方法,其特征在于,所述方法还包括:
直接基于所述船型文件的数据提取所述船体的横剖面线的反曲点,在提取每条横剖面线的反曲点时,采用逐步逼近法确定最逼近于所述横剖面线的逼近直线,遍历所述横剖面线上的所有点找到与所述逼近直线之间直线距离最远的点作为所述横剖面线的反曲点。
5.根据权利要求1所述的方法,其特征在于,所述对所述船型文件的数据进行积分求解处理提取所述船体的剖面面积和排水体积,包括:
用抛物线拟合船体型线并计算抛物线包含面积得到所述剖面面积;
用抛物线拟合所述剖面面积并计算抛物线包含面积得到所述排水体积。
6.根据权利要求5所述的方法,其特征在于,对于所述船体型线和所述剖面面积中的任意一种曲线,用抛物线拟合所述曲线并计算抛物线包含面积的方法为:
假设构成所述曲线的所有离散点共形成n个包络,每相邻两个离散点之间形成一个包络,则按照三个离散点为一组将离散点划分为若干组离散点,每组离散点中的三个离散点形成两个包络,n≥2;
基于每组离散点中的三个离散点拟合得到对应的子抛物线并求解得到所述子抛物线的面积,所述子抛物线的面积包括该组离散点形成的两个包络的总面积;
当n为偶数时,对每组离散点对应的子抛物线的面积进行求和得到所述抛物线包含面积;
当n为奇数时,基于未划分到任何一组离散点内的最后一个离散点及其最接近的前两个离散点拟合得到对应的子抛物线,并基于所述子抛物线求解得到所述最后一个离散点及其最接近的一个离散点之间形成的包络的落单包络面积,对求解得到的所述落单包络面积以及每组离散点对应的子抛物线的面积进行求和得到所述抛物线包含面积。
7.根据权利要求6所述的方法,其特征在于,
基于三个离散点拟合得到的对应的子抛物线表示为y=ax2+bx+c且:
Figure FDA0002520678540000021
其中,(x0,y0)、(x1,y1)和(x2,y2)分别为三个离散点的坐标值,h1=x1-x0,h2=x2-x1,l=h1+h2=x2-x0
当所述三个离散点属于同一组离散点且
Figure FDA0002520678540000022
时,同时计算三个离散点包含的面积得到该组离散点对应的子抛物线的面积:
Figure FDA0002520678540000023
当所述三个离散点属于同一组离散点且
Figure FDA0002520678540000031
或者
Figure FDA0002520678540000032
时,分别计算该组离散点形成的两个包络的面积并求和得到该组离散点对应的子抛物线的面积,两个包络的面积的求解公式为:
Figure FDA0002520678540000033
Figure FDA0002520678540000034
当所述三个离散点不属于同一组离散点,且(x2,y2)表示未划分到任何一组离散点内的最后一个离散点的坐标、(x1,y1)表示(x2,y2)最接近的前一个离散点的坐标、(x0,y0)表示(x2,y2)次接近的前一个离散点的坐标时,按照如下公式求解得到所述落单包络面积:
Figure FDA0002520678540000035
8.根据权利要求1所述的方法,其特征在于,所述对所述船型文件的数据进行微分求解处理提取所述船体在所述吃水深度下的横剖面船型特征值,包括对于每条横剖面线:
取所述横剖面线在0.2T~0.8T范围内的型线段,T表示所述吃水深度;
求得最逼近于所述型线段的直线的直线斜率;
确定所述横剖面线的横剖面船型特征值为
Figure FDA0002520678540000036
m为所述直线斜率,B表示所述船体的型宽。
9.根据权利要求8所述的方法,其特征在于,
求得的所述直线斜率的表达式为
Figure FDA0002520678540000037
其中,yj
Figure FDA0002520678540000038
处的型值;则求得的按绝对型值计算的所述横剖面船型特征值为:
Figure FDA0002520678540000039
10.根据权利要求9所述的方法,其特征在于,所述方法还包括:
对型值进行无量纲化得到无量纲参数
Figure FDA0002520678540000041
Figure FDA0002520678540000042
基于所述无量纲参数将按绝对型值计算的所述横剖面船型特征值处理为按无量纲型值计算的所述横剖面船型特征值表示为
Figure FDA0002520678540000043
CN202010489918.2A 2020-06-02 2020-06-02 一种船体型线特征参数提取方法 Active CN111591407B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010489918.2A CN111591407B (zh) 2020-06-02 2020-06-02 一种船体型线特征参数提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010489918.2A CN111591407B (zh) 2020-06-02 2020-06-02 一种船体型线特征参数提取方法

Publications (2)

Publication Number Publication Date
CN111591407A true CN111591407A (zh) 2020-08-28
CN111591407B CN111591407B (zh) 2021-07-27

Family

ID=72180937

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010489918.2A Active CN111591407B (zh) 2020-06-02 2020-06-02 一种船体型线特征参数提取方法

Country Status (1)

Country Link
CN (1) CN111591407B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115352598A (zh) * 2022-08-31 2022-11-18 大连海事大学 一种仰扫式船舶实载吨位计算方法及系统
CN115688273A (zh) * 2022-10-28 2023-02-03 武汉理工大学 一种长江中下游货船型线谱系生成方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB340000A (en) * 1929-06-14 1930-12-15 Erling Lars Kristofersen Improvements relating to screw propellers
SU1409518A1 (ru) * 1984-02-22 1988-07-15 Предприятие П/Я Г-4556 Носова оконечность судна ледового плавани
CN200942862Y (zh) * 2006-08-23 2007-09-05 国营武昌造船厂 一种大型客滚船船型
CN105916763A (zh) * 2013-12-20 2016-08-31 阿克北极科技公司 破冰船舶
CN109263799A (zh) * 2018-10-25 2019-01-25 中船黄埔文冲船舶有限公司 一种船体型线光顺方法
CN109263798A (zh) * 2018-10-25 2019-01-25 中船黄埔文冲船舶有限公司 一种艉柱不对称船舶船体型线放样方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB340000A (en) * 1929-06-14 1930-12-15 Erling Lars Kristofersen Improvements relating to screw propellers
SU1409518A1 (ru) * 1984-02-22 1988-07-15 Предприятие П/Я Г-4556 Носова оконечность судна ледового плавани
CN200942862Y (zh) * 2006-08-23 2007-09-05 国营武昌造船厂 一种大型客滚船船型
CN105916763A (zh) * 2013-12-20 2016-08-31 阿克北极科技公司 破冰船舶
CN109263799A (zh) * 2018-10-25 2019-01-25 中船黄埔文冲船舶有限公司 一种船体型线光顺方法
CN109263798A (zh) * 2018-10-25 2019-01-25 中船黄埔文冲船舶有限公司 一种艉柱不对称船舶船体型线放样方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张军等: "船舶前体横剖面UV度的系列变换及衡准", 《舰船科学》 *
盛振邦等: "《船舶原理》", 30 September 2003 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115352598A (zh) * 2022-08-31 2022-11-18 大连海事大学 一种仰扫式船舶实载吨位计算方法及系统
CN115688273A (zh) * 2022-10-28 2023-02-03 武汉理工大学 一种长江中下游货船型线谱系生成方法
CN115688273B (zh) * 2022-10-28 2023-08-15 武汉理工大学 一种长江中下游货船型线谱系生成方法

Also Published As

Publication number Publication date
CN111591407B (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
CN111591407B (zh) 一种船体型线特征参数提取方法
CN105701508B (zh) 基于多级卷积神经网络的全局‑局部优化模型及显著性检测算法
CN110599506B (zh) 一种复杂异形曲面机器人三维测量的点云分割方法
CN113450596A (zh) 基于船舶轨迹特征点提取的时空dp方法
CN102013102A (zh) 一种图像中的粘连细胞分割方法
KR102164137B1 (ko) 요트 선체 표면의 생성 및 파라메트릭 수정을 위한 새로운 설계 방법
CN108846844B (zh) 一种基于海天线的海面目标检测方法
CN112464366B (zh) 基于数据挖掘的自主水下航行器多保真度外形优化方法
CN104408260A (zh) 一种潮流能水轮机叶片翼型设计方法
CN105550691A (zh) 基于尺度空间的自适应山谷山脊线提取方法及系统
CN104156730B (zh) 一种基于骨架的抗噪声汉字特征提取方法
CN104137119B (zh) 图像处理装置及图像处理方法
Vernengo et al. Numerical investigation on the hydrodynamic performance of fast SWATHs with optimum canted struts arrangements
CN107330901A (zh) 一种基于骨架的物体构件分解方法
CN110580388B (zh) 一种基于众源轨迹数据的航道网络提取方法
CN115146412A (zh) 一种船型性能预报与优化方法及系统
CN113409332B (zh) 一种基于三维点云的建筑物平面分割方法
CN115100439B (zh) 一种用于水产养殖的水体异常检测方法
KR101227953B1 (ko) 선박의 선형 설계 방법 및 장치
CN110084779B (zh) 一种基于激光扫描的飞机厚蒙皮端面特征点提取与去噪方法
CN115457277A (zh) 一种智能化路面病害识别检测方法及系统
CN109241975B (zh) 一种基于字符中心点定位的车牌字符分割方法
CN102930558B (zh) 一种多特征融合的红外图像目标实时跟踪方法
CN104851105B (zh) 一种改进的基于分水岭变换的泡沫图像分割方法
CN102222220B (zh) 手背静脉模式纹理提取方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant