CN111336957B - 一种基于数学建模的浮船坞挠度精确测量方法 - Google Patents
一种基于数学建模的浮船坞挠度精确测量方法 Download PDFInfo
- Publication number
- CN111336957B CN111336957B CN202010155392.4A CN202010155392A CN111336957B CN 111336957 B CN111336957 B CN 111336957B CN 202010155392 A CN202010155392 A CN 202010155392A CN 111336957 B CN111336957 B CN 111336957B
- Authority
- CN
- China
- Prior art keywords
- dock
- measuring
- deflection
- point
- fore
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B13/00—Measuring arrangements characterised by the use of fluids
- G01B13/24—Measuring arrangements characterised by the use of fluids for measuring the deformation in a solid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
- G01M5/005—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种基于数学建模的浮船坞挠度精确测量方法,包括以下步骤:在浮船坞左右坞墙内,沿艏、艉方向各均匀安装数个垂直于船坞甲板的测量管,且测量管的底部用连通管相互连通,向测量管内注入水或其他液体,在每根测量管下部的同一水平位置各安装一个液位传感器,用于测量各个测量管内的液柱高度;将包括首末测量点在内的所有测量点的数据全部纳入处理范围,通过最小二乘法和拉格朗日插值法相结合的方法拟合计算出各测量点未变形的液位高度;通过有限的各测量点挠度值,拟合出船坞艏艉方向各连续点的挠度函数模型。本发明通过对有限的测量数据,拟合出近似的挠度数学函数模型,进而可以预测出船坞纵向任意点的挠度情况。
Description
技术领域
本发明涉及浮船坞挠度的测量方法,尤其是涉及一种基于数学建模的浮船坞挠度精确测量方法
背景技术
浮船坞是一种通过其底部压载水舱的水量进行压载和泵排来控制其浮沉的修船平台,是修造船业的重要设施工具。在船只进出坞的时候,浮船坞需要下沉和上浮操作,中间过程中坞体会产生一定的形变,称为挠度。浮船坞挠度的实时监测对浮船坞的安全操作非常重要。
目前常用的几种浮船坞挠度测量数据处理方法,一方面对数据的利用率不够,会对首末的挠度数据进行理想处理,定性为不发生形变;另一方面只能测量设有测量点处的挠度值,无法得知未设置测量点地方的挠度值,进而无法精确控制各个水仓的进水和出水,导致船坞应力受力不均,容易导致安全事故。
发明内容
针对现有技术不足,本发明提供了一种基于数学建模的浮船坞挠度精确测量方法。
本发明解决上述技术问题所采用的技术方案为:
一种基于数学建模的浮船坞挠度精确测量方法,其包括以下步骤:在浮船坞左右坞墙内,沿艏、艉方向各均匀安装数个垂直于船坞甲板的测量管,且测量管的底部用连通管相互连通,向测量管内注入水或其他液体,在每根测量管下部的同一水平位置各安装一个液位传感器,用于测量各个测量管内的液柱高度;将包括首末测量点在内的所有测量点的数据全部纳入处理范围,通过最小二乘法和拉格朗日插值法相结合的方法拟合计算出各测量点未变形的液位高度;通过有限的各测量点挠度值,拟合出船坞艏艉方向各连续点的挠度函数模型,根据挠度函数模型可以计算浮船坞艏、艉方向任意点的挠度值。
进一步地,所述的基于数学建模的浮船坞挠度精确测量方法,包括以下具体步骤:
(1)将船坞方向设为X轴,连通管液位高度设为Y轴,挠度值设为Z轴,同时艏、艉方向设有n点测量管;
(2)将船坞长度设为L,则坞艏X坐标为0,坞艉X坐标为L,各测量点与坞艏的距离为各测量点的X坐标值,记为X1,X2,X3…,Xn;
(3)记n个连通管液位高度由艏至艉分别为Y1,Y2,Y3…,Yn,则得到一组坐标数据(X1,Y1),(X2,Y2),(X3,Y3)…,(Xn,Yn);
(4)当浮船坞不存在变形时,测量点位置和测量管液位高度为一次线性关系,通过最小二乘法和拉格朗日插值法可分别拟出两条直线方程Y=a0+a1X和Y=b0+b1X,两者相结合,系数均值处理得到的新拟合函数Y=(a0+b0)/2+(a1+b1)/2X,记为Y=c0+c1X,则各测量点的挠度值为:Z1=Y1-(c0+c1X1),Z2=Y2-(c0+c1X2),Z3=Y3-(c0+c1X3)…,Zn=Yn-(c0+c1Xn);
(5)根据步骤(4)可得到一组测量点与挠度值的新数据(X1,Z1),(X2,Z2),(X3,Z3)…,(Xn,Zn),通过以上n组数据拟合出挠度与测量点的近似多项式组成的数学函数,即可以计算出船坞艏艉方向任意一点的挠度值,将具有n组数据的测量数值,利用拉格朗日插值法可拟合出(n-1)次的多项式Z=f(X)=Z=e0+e1X+e2X2+…+e(n-1)X(n-1),其中X为浮船坞艏艉测量管的连通管方向上任意一点距坞艏的距离,Z为该点对应的挠度值。
进一步地,在浮船坞艏、艉方向设置所述测量管的数量为3-7个。
与现有技术相比,本发明具备的有益效果:其一,优化了对现有技术中首末两个测量点的数据处理,不再单纯将此两点作为参考点使用;其二,将最小二乘法和拉格朗日插值法两者相结合,优势叠加,优于单独使用一种数值处理方法;其三,通过对有限的测量数据,拟合出近似的挠度数学函数模型,进而可以预测出船坞纵向任意点的挠度情况。
附图说明
图1为采用本发明挠度精确测量方法的浮船坞处于正常水平无变形状态的结构示意图;
图2为采用本发明挠度精确测量方法的浮船坞处于艏倾甲板变形状态的结构示意图。
其中,1、浮船坞,11、船坞甲板,12、坞艏,13、坞艉,21、测量管,22、连通管,23、未变形水位高度拟合线。
具体实施方式
下面结合附图和具体实施例对本发明做进一步地说明。
一种基于数学建模的浮船坞挠度精确测量方法,其包括以下步骤:在浮船坞1左右坞墙内,沿坞艏12、坞艉13方向各均匀安装数个垂直于船坞甲板11的测量管21,且测量管21的底部用连通管22相互连通,向测量管21内注入大约1/2测量管高度的水或其他液体,在每根测量管21下部的同一水平位置各安装一个液位传感器,用于测量各个测量管21内的液柱高度;在正常情况下,如图1所示,如果浮船坞甲板没有变形,则各测量管21内的液柱高度基本相同,此时各测量管21内的液柱高度记为未变形水位高度拟合线23;而如果坞体发生变形,会致使各测量管21的相对位置发生变化,如图2所示,从而引起各测量管21内的液体高度发生变化,通过对测量点及测量管21内液体高度的数据进行处理,从而获得对应测量点的挠度值。对应测量点的挠度值的测量方法为:将包括首末测量点在内的所有测量点的数据全部纳入处理范围,通过最小二乘法和拉格朗日插值法相结合的方法拟合计算出各测量点未变形的液位高度;通过有限的各测量点挠度值,拟合出船坞艏艉方向各连续点的挠度函数模型,根据挠度函数模型可以计算浮船坞艏、艉方向任意点的挠度值。
进一步地,所述的基于数学建模的浮船坞挠度精确测量方法,包括以下具体步骤:
(1)将船坞方向设为X轴,连通管液位高度设为Y轴,挠度值设为Z轴,同时在本实施例中坞艏、艉方向设有7点测量管,在其他实施例中还可以选择不同数量的、间隔相同的测量管,若为了增加测量精度,可选择设置更多的测量点;
(2)将船坞长度设为L,如图1、2所示,则坞艏12的X坐标为0,坞艉13的X坐标为L,各测量点与坞艏的距离为各测量点的X坐标值,记为X1、X2、X3、X4、X5、X6、X7;
(3)记7个连通管液位高度由艏至艉分别为Y1、Y2、Y3、Y4、Y5、Y6、Y7,则得到一组坐标数据(X1,Y1)、(X2,Y2)、(X3,Y3)、(X4,Y4)、(X5,Y5)、(X6,Y6)、(X7,Y7);
(4)当浮船坞不存在变形时,测量点位置和测量管21液位高度应为一次线性关系,变形导致的实测数据分布于这条线的附近,通过最小二乘法和拉格朗日插值法可分别拟出两条直线方程Y=a0+a1X和Y=b0+b1X,两者相结合,系数均值处理得到的新拟合函数Y=(a0+b0)/2+(a1+b1)/2X,记为Y=c0+c1X,为最接近无变形时的直线方程,则各测量点的挠度值为:Z1=Y1-(c0+c1X1)、Z2=Y2-(c0+c1X2)、Z3=Y3-(c0+c1X3)、Z4=Y4-(c0+c1X4)、Z5=Y5-(c0+c1X5)、Z6=Y6-
(c0+c1X6)、Z7=Y7-(c0+c1X7),关于最小二乘法和朗格朗日插值法的原理作为公开原理,在此就不再赘述;
(5)根据步骤(4)可得到一组测量点与挠度值的新数据(X1,Z1)、(X2,Z2)、(X3,Z3)、(X4,Z4)、(X5,Z5)、(X6,Z6)、(X7,Z7),通过以上7组数据拟合出挠度与测量点的近似多项式组成的数学函数,即可以计算出船坞艏艉方向任意一点的挠度值,将本实施例具有的7组数据测量数值,利用拉格朗日插值法可拟合出6次的多项式Z=f(X)=Z=e0+e1X+e2X2+e3X3+e4X4+e5X5+e6X6,其中X为浮船坞艏艉测量管的连通管方向上任意一点距坞艏的距离,Z为该点对应的挠度值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于数学建模的浮船坞挠度精确测量方法,其特征在于,包括以下步骤:在浮船坞左右坞墙内,沿艏、艉方向各均匀安装数个垂直于船坞甲板的测量管,且测量管的底部用连通管相互连通,向测量管内注入水或其他液体,在每根测量管下部的同一水平位置各安装一个液位传感器,用于测量各个测量管内的液柱高度;将包括首末测量点在内的所有测量点的数据全部纳入处理范围,通过最小二乘法和拉格朗日插值法相结合的方法拟合计算出各测量点未变形的液位高度;通过有限的各测量点挠度值,拟合出船坞艏艉方向各连续点的挠度函数模型;
所述的基于数学建模的浮船坞挠度精确测量方法,具体步骤为:
(1)将船坞方向设为X轴,连通管液位高度设为Y轴,挠度值设为Z轴,同时艏、艉方向设有n点测量管;
(2)将船坞长度设为L,则坞艏X坐标为0,坞艉X坐标为L,各测量点与坞艏的距离为各测量点的X坐标值,记为X1,X2,X3…,Xn;
(3)记n个连通管液位高度由艏至艉分别为Y1,Y2,Y3…,Yn,则得到一组坐标数据(X1,Y1) ,(X2,Y2) ,(X3,Y3)…,(Xn,Yn);
(4)当浮船坞不存在变形时,测量点位置和测量管液位高度为一次线性关系,通过最小二乘法和拉格朗日插值法可分别拟出两条直线方程Y=a0+a1X和Y=b0+b1X,两者相结合,系数均值处理得到的新拟合函数Y=(a0+b0)/2+(a1+b1)/2X,记为Y=c0+c1X,则各测量点的挠度值为:Z1=Y1-(c0+c1X1),Z2=Y2-(c0+c1X2),Z3=Y3-(c0+c1X3)…,Zn=Yn-(c0+c1Xn);
(5)根据步骤(4)可得到一组测量点与挠度值的新数据(X1,Z1),(X2,Z2) ,(X3,Z3)…,(Xn,Zn),通过以上n组数据拟合出挠度与测量点的近似多项式组成的数学函数,即可以计算出船坞艏艉方向任意一点的挠度值,将具有n组数据的测量数值,利用拉格朗日插值法可拟合出(n-1)次的多项式Z=f(X)=Z=e0+e1X+e2X2+…+e(n-1)X(n-1),其中X为浮船坞艏艉测量管的连通管方向上任意一点距坞艏的距离,Z为该点对应的挠度值。
2.根据权利要求1所述的一种基于数学建模的浮船坞挠度精确测量方法,其特征在于在浮船坞艏、艉方向设置所述测量管的数量为3-7个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010155392.4A CN111336957B (zh) | 2020-03-09 | 2020-03-09 | 一种基于数学建模的浮船坞挠度精确测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010155392.4A CN111336957B (zh) | 2020-03-09 | 2020-03-09 | 一种基于数学建模的浮船坞挠度精确测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111336957A CN111336957A (zh) | 2020-06-26 |
CN111336957B true CN111336957B (zh) | 2022-03-11 |
Family
ID=71182240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010155392.4A Active CN111336957B (zh) | 2020-03-09 | 2020-03-09 | 一种基于数学建模的浮船坞挠度精确测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111336957B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112284325A (zh) * | 2020-11-09 | 2021-01-29 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种升船机承船厢挠度动态测试方法 |
CN112623150B (zh) * | 2020-12-09 | 2022-12-20 | 大连中远海运重工有限公司 | 油船在漂浮状态下甲板整体换新的方法 |
CN114861252A (zh) * | 2021-02-04 | 2022-08-05 | 西安电子科技大学青岛计算技术研究院 | 一种基于bim+gis技术的智慧水库管理系统及其实现方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB773160A (en) * | 1955-01-31 | 1957-04-24 | Thomas Mercer Ltd | Improvements relating to pneumatic gauging devices for effecting linear measurements |
WO2009140233A2 (en) * | 2008-05-12 | 2009-11-19 | Smith David Q | Floating dock deflection management systems |
CN201569421U (zh) * | 2009-12-17 | 2010-09-01 | 杨光祥 | 一种浮船坞挠度监测装置 |
CN203100703U (zh) * | 2012-12-25 | 2013-07-31 | 上海中远船务工程有限公司 | 浮船坞挠度测试装置 |
CN203745170U (zh) * | 2014-03-17 | 2014-07-30 | 重庆交通大学 | 桥梁静态挠度监测装置 |
CN103996038A (zh) * | 2014-05-01 | 2014-08-20 | 朱毅 | 一种调整手指静脉识别图像倾斜角度的方法 |
CN103995492A (zh) * | 2014-04-30 | 2014-08-20 | 中国人民解放军总后勤部军事交通运输研究所 | 浮船坞升沉自动控制系统 |
CN104139836A (zh) * | 2014-07-24 | 2014-11-12 | 江苏驷博电气有限公司 | 浮船坞扰度测量装置 |
CN204096060U (zh) * | 2014-07-24 | 2015-01-14 | 江苏驷博电气有限公司 | 浮船坞扰度测量装置 |
CN108716904A (zh) * | 2018-04-09 | 2018-10-30 | 水利部南京水利水文自动化研究所 | 基于有限测斜仪测点测值的坝体挠度获取方法 |
CN109186531A (zh) * | 2018-09-26 | 2019-01-11 | 招商局重工(深圳)有限公司 | 一种集中水箱式挠度测量系统 |
-
2020
- 2020-03-09 CN CN202010155392.4A patent/CN111336957B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB773160A (en) * | 1955-01-31 | 1957-04-24 | Thomas Mercer Ltd | Improvements relating to pneumatic gauging devices for effecting linear measurements |
WO2009140233A2 (en) * | 2008-05-12 | 2009-11-19 | Smith David Q | Floating dock deflection management systems |
CN201569421U (zh) * | 2009-12-17 | 2010-09-01 | 杨光祥 | 一种浮船坞挠度监测装置 |
CN203100703U (zh) * | 2012-12-25 | 2013-07-31 | 上海中远船务工程有限公司 | 浮船坞挠度测试装置 |
CN203745170U (zh) * | 2014-03-17 | 2014-07-30 | 重庆交通大学 | 桥梁静态挠度监测装置 |
CN103995492A (zh) * | 2014-04-30 | 2014-08-20 | 中国人民解放军总后勤部军事交通运输研究所 | 浮船坞升沉自动控制系统 |
CN103996038A (zh) * | 2014-05-01 | 2014-08-20 | 朱毅 | 一种调整手指静脉识别图像倾斜角度的方法 |
CN104139836A (zh) * | 2014-07-24 | 2014-11-12 | 江苏驷博电气有限公司 | 浮船坞扰度测量装置 |
CN204096060U (zh) * | 2014-07-24 | 2015-01-14 | 江苏驷博电气有限公司 | 浮船坞扰度测量装置 |
CN108716904A (zh) * | 2018-04-09 | 2018-10-30 | 水利部南京水利水文自动化研究所 | 基于有限测斜仪测点测值的坝体挠度获取方法 |
CN109186531A (zh) * | 2018-09-26 | 2019-01-11 | 招商局重工(深圳)有限公司 | 一种集中水箱式挠度测量系统 |
Non-Patent Citations (4)
Title |
---|
一种求最佳直线的新方法;张弢 等;《沈阳工业大学学报》;20020630;第24卷(第3期);262-264 * |
智能压差式液体密度计的设计与实现;杨云涛;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20140315;C030-5 * |
浮船坞实时数据监控系统的设计;唐家和 等;《青岛远洋船员学院学报》;20030515(第2期);45-49 * |
浮船坞挠度自动监测系统设计;李志;《船舶工程》;20100915;第32卷(第S1期);58-60 * |
Also Published As
Publication number | Publication date |
---|---|
CN111336957A (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111336957B (zh) | 一种基于数学建模的浮船坞挠度精确测量方法 | |
CN102285431A (zh) | 通航船舶吃水量实时测量装置及测量方法 | |
CN114636383B (zh) | 一种沉管隧道管节施工过程动态形变测量方法 | |
CN104392138B (zh) | 一种船舶破舱进水应急决策辅助方法 | |
CN113720296A (zh) | 一种沉管水下形变监测方法 | |
CN111024178A (zh) | 一种浮式海洋平台压载监测系统及舱中液体体积计算方法 | |
CN110523941A (zh) | 连铸粘结漏钢多级风险控制方法及控制装置 | |
EP0372631B1 (en) | A method and system for measuring a vertical density profile of a fluid | |
KR101403347B1 (ko) | 압력강하 및 전열성능을 고려한 증기발생기 관막음률 환산 장치 | |
AU2018344386B2 (en) | Method for determining the geometry of a defect and for determining a load limit | |
EP3532712B1 (en) | System and method to determine a flow rate of a washing liquid at an exhaust gas inlet of an exhaust gas cleaning system of a marine vessel | |
CN110228571B (zh) | 一种耙吸挖泥船船舶吃水计算优化方法 | |
CN115592294A (zh) | 焊材在管道焊接的外环焊缝全位置成形质量评价方法 | |
CN105625484B (zh) | 一种大型物品深水测控系统和方法 | |
KR20170126624A (ko) | 선박의 평형계측 장치 및 평형계측 방법 | |
CN101354279A (zh) | 铁水罐车重量计算液面的方法 | |
CN104006938A (zh) | 一种横管冷却器泄漏点的检出方法 | |
CN106647599A (zh) | 一种中间包对中检测及控制系统与方法 | |
KR20160031987A (ko) | 연료 질량 및 연료 밀도를 결정하기 위한 방법 | |
CN114954854A (zh) | 水下无人航行器水舱容积建模系统及建模方法 | |
CN112325815B (zh) | 一种集装箱船加长改装的三维检测方法 | |
CN112033347A (zh) | 一种水下遥控载具及其姿态检测方法 | |
KR20120057762A (ko) | 탕면부 유동 측정장치 및 측정방법 | |
CN112180435A (zh) | 辅船拖带的震源的位置监控方法及装置 | |
CN216049593U (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 |