CN104764413B - 海洋结构物甲板片焊接变形测量方法 - Google Patents
海洋结构物甲板片焊接变形测量方法 Download PDFInfo
- Publication number
- CN104764413B CN104764413B CN201410811276.8A CN201410811276A CN104764413B CN 104764413 B CN104764413 B CN 104764413B CN 201410811276 A CN201410811276 A CN 201410811276A CN 104764413 B CN104764413 B CN 104764413B
- Authority
- CN
- China
- Prior art keywords
- msub
- deck plate
- girder
- point
- mrow
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了海洋结构物甲板片焊接变形测量新方法,它包括以下步骤:在CAD三维软件中建立海洋结构物甲板片的三维模型,并对甲板片的主梁外轮廓直线信息进行提取;使用三维激光扫描系统对焊接后的甲板片进行扫描,得到实际焊接后甲板片的点云数据;将甲板片的三维模型的轮廓信息进行坐标变换,统一到点云数据的坐标系下;搜索甲板片主梁外轮廓周围区域内的点获得特征邻域点云集合;对特征邻域点云集合进行分析计算,提取出实际的甲板片主梁外轮廓直线;计算出主梁两个中线交点的坐标;求取所有需要测量的甲板片主梁交点位置坐标;对比实际主梁两个中线交点坐标与理论主梁两个中线交点坐标的偏差。采用本方法测量精度和效率可以大幅提高。
Description
技术领域
本发明涉及一种焊接变形测量方法,尤其涉及海洋结构物甲板片的焊接变形测量方法。
背景技术
海洋结构物甲板片焊接后往往会产生变形,因此就需要在焊接后对甲板片进行测量,检验出焊接质量是否合格。当前多数海洋结构物甲板片焊接变形测量采用全站仪对甲板片主梁的交点进行逐点测量,然后与设计的交点坐标进行比对检验。该方法效率较低,而且容易产生测量误差。
发明内容
本发明的目的在于克服已有技术的缺点,提供一种测量精度和效率可以大幅提高的海洋结构物甲板片焊接变形测量方法。
为了达到上述目的,本发明采用的技术方案是:
本发明的海洋结构物甲板片焊接变形测量方法,它包括以下步骤:
(1)在CAD三维软件中建立海洋结构物甲板片的三维模型,并对甲板片的主梁外轮廓直线信息进行提取;
(2)使用三维激光扫描系统对焊接后的甲板片进行扫描,得到实际焊接后甲板片的点云数据;
(3)将甲板片的三维模型的轮廓信息进行坐标变换,统一到点云数据的坐标系下;
(4)搜索甲板片主梁外轮廓周围区域内的点获得特征邻域点云集合;
(5)对特征邻域点云集合进行分析计算,提取出实际的甲板片主梁外轮廓直线;
(6)通过计算得到实际的甲板片的主梁中线,从而计算出主梁两个中线交点的坐标;
(7)重复所述的步骤(1)至步骤(6),求取所有需要测量的甲板片主梁交点位置坐标;
(8)对比实际主梁两个中线交点坐标与理论主梁两个中线交点坐标的偏差,检验出焊接质量是否合格。
本发明的优点在于:采用三维激光扫描系统扫描甲板片,可一次性提取出需要测量的交点信息,相比于以往的全站仪逐点测量,本方法的测量精度和效率可以大幅提高;因为海洋结构物甲板片焊接变形测量主要测量的是主梁焊接交点,所以在CAD模型辅助下处理点云数据,将减少其他不必要的数据计算工作,提高效率。
附图说明
图1是本发明的海洋结构物甲板片焊接变形测量方法采用的坐标计算三维示意图;
图2是发明的海洋结构物甲板片焊接变形测量方法采用的特征信息周围区域点云搜索示意图。
图3是本发明的海洋结构物甲板片焊接变形测量方法采用的轮廓线中线及其交点示意图。
具体实施方式
下面结合具体实施例对本发明进行详细描述。
本发明是海洋结构物甲板片焊接变形测量方法,它包括以下步骤:
(1)在CAD三维软件中建立海洋结构物甲板片焊接后的理论三维模型,并对甲板片的主梁外轮廓直线信息进行提取,甲板片的主梁的结构形式为工字钢焊接而成的钢架结构;
(2)使用三维激光扫描系统对甲板片进行扫描,得到实际焊接后甲板片的点云数据,然后对点云数据进行优化处理,得到优化的实际焊接后甲板片点云数据;优化处理通常包括去噪等处理。三维激光扫描系统可以在市场购得,三维激光扫描系统的工作原理为:
三维激光扫描系统是采用脉冲激光测距的方法,通过激光脉冲二极管发射激光脉冲,射向被测物体,激光碰到被测物体的表面后反射回来,由系统内的探测装置接收并记录。根据激光束发射和接收的时间差T,计算扫描点到仪器的斜距S,同时获取此激光束的水平方向角度α和垂直方向角度θ,计算出扫描点用空间坐标系中X方向、Y方向、Z方向3个矢量数据表示。见图1。
计算扫描点空间坐标时,坐标原点为系统内部测距激光发射的几何中心,Z轴为仪器的竖轴,Y轴为扫描默认起始方位,X轴与Y、Z轴垂直构成空间坐标系。通过计算可得扫描点P的坐标(Xp,Yp,Zp),即
式(1)中距离是利用激光发射和接收的时间差T来计算的,为:
式(2)中C0为光速。
(3)将提取的理论焊接后甲板片主梁的轮廓线信息进行坐标变换,统一到点云数据的坐标系下;
坐标变换的具体算法如下:
取三维模型轮廓信息坐标系的原点为P0和在X轴、Y轴上各取一点P1和P2,点云数据坐标系的原点Q0,取X轴、Y轴上各一点Q1和Q2。
以P0点做为原点,P0到P1方向做为X轴方向,用上述点构造理论模型坐标系,构建得到的单位向量:
Z轴方向确定是通过X轴方向向量e1与P1到P2的方向矢量的叉乘,构建出的Z轴的单位向量为:
再以e1、e3的叉乘计算得到Y轴方向的单位向量:
e2=e1×e3 (5)
同上,通过点云数据中的点Q0、Q1、Q2进行点云数据中三个坐标轴的计算,得到的单位向量:
e′2=e′1×e′3 (8)
在坐标变换完成后,两个坐标系基本可以达到重合的状态,若M和M’为同一向量在这两个坐标系中的坐标,则:
[e′1 e′2 e′3]M′=[e1 e2 e3]T×M (9)
因为是一个正交矩阵,则匹配矩阵里的旋转R为:
R=[e′1 e′2 e′]×[e1 e2 e3]T (10)
M′=([e′1 e′2 e′3]T)-1×[e1 e2 e3]T×M (11)
匹配矩阵的平移矩阵T为:
T=Q0-P0 (12)
完整的匹配矩阵为:
(4)以经过坐标变换的理论焊接后甲板片主梁的轮廓线为轴线,对每一个轮廓特征线周围设定区域的领域点云数据进行搜索;
邻域点云数据搜索具体步骤为:
步骤一、将理论三维模型的轮廓特征信息整合到点云数据的坐标系后,通过人工移动使点云数据和理论三维模型的轮廓线位置基本重合;
步骤二、以经过坐标变换的理论焊接后甲板片主梁的轮廓特征线为轴线,生成圆柱体,然后用球体对特征周围形成的圆柱体进行空间填充,用一系列球体所包含的点集近似代替圆柱体内部点;
步骤三、使用k-邻域点云搜索方法对轮廓线周围进行邻域搜索;例如计算某空间点p的k-邻域的方法是求出其余n-1个点到点p的空间距离di,i∈[1,n-1],然后对其进行从小到大排序,取前k个点作为该空间点的k-邻域。该方法详见2013年5月南京师范大学硕士学位论文“影响特征线辅助下的三维激光点云建筑物建模”一文。
(5)将搜索到的邻域点云数据进行轮廓线提取,获得实际焊接后甲板片主梁的轮廓线;
轮廓线提取方法具体如下:
由于主梁轮廓是甲板片的上面与其轴向两个侧面的分界线,在分界线处点云法向或曲率存在突变。此处为了降低其他因素的干扰,步骤一、先将点云数据中具有法向或曲率突变的点舍去;步骤二、根据甲板片上面与其轴向两个侧面内点云法向指向不同特性进行区域分类,然后得到上表面和与其相交的各面两区域点云;步骤三、利用最小二乘法拟合两区域空间平面P1、P2,数学表达式为:
式中a1、b1、c1、d1、a2、b2、c2、d2分别为空间平面P1、P2表达式中的系数,x、y、z分别为两面上点的坐标。
联立两空间平面方程,求取两空间平面的交线来表示实际焊接后甲板片主梁的轮廓线。相交空间直线方程数学表示:
式(15)中,(x0,y0,z0)为空间直线上任意一点,为空间直线的方向向量。
(6)每条主梁都可获得两条轮廓线,计算两轮廓线中线,然后计算出中线的交点。
如图3所示,具体步骤如下;
第一步,计算出L1与L4交点D1,L3与L4交点D2。
第二步,计算出空间平行线L1与L3的中间线L2的表达式,即:
式(16)中,(x3,y3,z3)为点D1和点D2的中点坐标,为空间直线L1的方向向量。
第三步,重复以上两步,计算出点D3坐标,然后计算出空间平行线L4与L6的中间线L5的表达式。
第四步,将计算得到的两条中间线L2和L5进行联立,解方程后得出交点坐标。
(7)重复所述的步骤(1)至步骤(6),求取所有需要测量的甲板片主梁交点位置坐标;
(8)将实际焊接后甲板片主梁各交点位置的坐标与理论坐标进行比对,分析出偏差是否在控制范围内,检验焊接后甲板片是否合格。
本方法中使用三维激光扫描系统,对甲板片进行扫描,一次性获取甲板片主梁交点;通过CAD软件三维建模,在三维模型的辅助下对点云数据进行处理,能够准确快速的获得甲板片主梁的三维信息,减少对其他部分不必要的分析计算,提高测量的精度和效率。
Claims (3)
1.海洋结构物甲板片焊接变形测量方法,其特征在于,它包括以下步骤:
(1)在CAD三维软件中建立海洋结构物甲板片的三维模型,并对甲板片的主梁外轮廓直线信息进行提取;
(2)使用三维激光扫描系统对焊接后的甲板片进行扫描,得到实际焊接后甲板片的点云数据;
(3)将甲板片的三维模型的轮廓信息进行坐标变换,统一到点云数据的坐标系下;
(4)搜索甲板片主梁外轮廓周围区域内的点获得特征邻域点云集合;
(5)对特征邻域点云集合进行分析计算,提取出实际的甲板片主梁外轮廓直线;
(6)通过计算得到实际的甲板片的主梁中线,从而计算出主梁两个中线交点的坐标;
(7)重复所述的步骤(1)至步骤(6),求取所有需要测量的甲板片主梁交点位置坐标;
(8)对比实际主梁两个中线交点坐标与理论主梁两个中线交点坐标的偏差,检验出焊接质量是否合格。
2.根据权利要求1所述的海洋结构物甲板片焊接变形测量方法,其特征在于:所述的步骤(4)中邻域点云数据搜索具体步骤为:
步骤一、将理论三维模型的轮廓特征信息整合到点云数据的坐标系后,通过人工移动使点云数据和理论三维模型的轮廓线位置基本重合;
步骤二、以经过坐标变换的理论焊接后甲板片主梁的轮廓特征线为轴线,生成圆柱体,然后用球体对特征周围形成的圆柱体进行空间填充,用一系列球体所包含的点集近似代替圆柱体内部点;
步骤三、使用k-邻域点云搜索方法对轮廓线周围进行邻域搜索。
3.根据权利要求1或2所述的海洋结构物甲板片焊接变形测量方法,其特征在于:所述的步骤(5)中特征提取方法具体如下:
步骤一、先将点云数据中具有法向或曲率突变的点舍去;
步骤二、根据甲板片上面与其轴向两个侧面内点云法向指向不同特性进行区域分类,然后得到上表面和与其相交的各面两区域点云;
步骤三、利用最小二乘法拟合两区域空间平面,数学表达式为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mi>z</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>:</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mi>z</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中a1、b1、c1、d1、a2、b2、c2、d2分别为空间平面P1、P2表达式中的系数,x、y、z分别为两面上点的坐标;
联立两空间平面方程,求取两空间平面的交线来表示实际焊接后甲板片主梁的轮廓线;相交空间直线方程数学表示:
<mrow>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
</mfrac>
</mrow>
式中,(x0,y0,z0)为空间直线上任意一点,(θ1,θ2,θ3)为空间直线的方向向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410811276.8A CN104764413B (zh) | 2014-12-22 | 2014-12-22 | 海洋结构物甲板片焊接变形测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410811276.8A CN104764413B (zh) | 2014-12-22 | 2014-12-22 | 海洋结构物甲板片焊接变形测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104764413A CN104764413A (zh) | 2015-07-08 |
CN104764413B true CN104764413B (zh) | 2017-10-27 |
Family
ID=53646393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410811276.8A Active CN104764413B (zh) | 2014-12-22 | 2014-12-22 | 海洋结构物甲板片焊接变形测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104764413B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105783782B (zh) * | 2016-05-12 | 2018-06-15 | 电子科技大学 | 表面曲率突变光学轮廓测量方法 |
CN107393006B (zh) * | 2017-07-26 | 2021-02-12 | 河海大学 | 一种衡量隧道整体变形的方法 |
CN108180878A (zh) * | 2017-11-21 | 2018-06-19 | 博迈科海洋工程股份有限公司 | 一种大型模块结构中甲板片安装误差控制方法 |
CN110426001A (zh) * | 2019-08-30 | 2019-11-08 | 四川大学 | 一种基于三维激光扫描的危岩体旋转位移监测方法 |
CN110425995B (zh) * | 2019-08-30 | 2020-04-14 | 四川大学 | 一种基于点云平均域矢量算法的滑坡变形监测方法 |
CN111811442A (zh) * | 2020-07-14 | 2020-10-23 | 上海船舶工艺研究所(中国船舶工业集团公司第十一研究所) | 一种船舶大面积甲板平面度快速测量与计算方法 |
CN112478082B (zh) * | 2020-12-11 | 2022-07-08 | 天津博迈科海洋工程有限公司 | 大型海洋油气模块多层甲板建造精度控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288123A (zh) * | 2011-06-27 | 2011-12-21 | 清华大学 | 一种用于焊接变形三维非接触式测量的精确定位方法 |
CN102607446A (zh) * | 2011-10-21 | 2012-07-25 | 中建工业设备安装有限公司 | 基于网格测量钢结构焊接变形的方法 |
CN103940356A (zh) * | 2014-02-27 | 2014-07-23 | 山东交通学院 | 一种基于三维激光扫描技术的建筑物整体变形监测方法 |
CN104006757A (zh) * | 2014-06-13 | 2014-08-27 | 苏州西博三维科技有限公司 | 焊接变形实时检测系统及其检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007298343A (ja) * | 2006-04-28 | 2007-11-15 | Tokyo Institute Of Technology | 溶接変形および残留応力の評価方法 |
KR100910090B1 (ko) * | 2007-05-31 | 2009-07-30 | 주식회사 한진중공업 | 용접이면의 변형량 측정장치 |
-
2014
- 2014-12-22 CN CN201410811276.8A patent/CN104764413B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288123A (zh) * | 2011-06-27 | 2011-12-21 | 清华大学 | 一种用于焊接变形三维非接触式测量的精确定位方法 |
CN102607446A (zh) * | 2011-10-21 | 2012-07-25 | 中建工业设备安装有限公司 | 基于网格测量钢结构焊接变形的方法 |
CN103940356A (zh) * | 2014-02-27 | 2014-07-23 | 山东交通学院 | 一种基于三维激光扫描技术的建筑物整体变形监测方法 |
CN104006757A (zh) * | 2014-06-13 | 2014-08-27 | 苏州西博三维科技有限公司 | 焊接变形实时检测系统及其检测方法 |
Non-Patent Citations (5)
Title |
---|
大型船体焊接变形仿真技术探究及其应用;李靖;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20100715(第07(2010)期);C036-35 * |
应用三维激光扫描法测量板材的焊接变形;何洪文等;《焊接学报》;20111231;第32卷(第12期);9-13 * |
影像特征线辅助下的三维激光点云建筑物建模;王波;《中国优秀硕士学位论文全文数据库基础科学辑》;20140215(第02(2014)期);正文第42-45页 * |
点云平面拟合在三维激光扫描仪变形监测中的应用;蔡来良等;《测绘科学》;20100930;第35卷(第05期);231-232 * |
焊接变形测量方法研究进展;刘建伟等;《热加工工艺》;20140930;第43卷(第17期);5-9、14 * |
Also Published As
Publication number | Publication date |
---|---|
CN104764413A (zh) | 2015-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104764413B (zh) | 海洋结构物甲板片焊接变形测量方法 | |
CN104316016B (zh) | 一种整体叶盘复杂曲面叶片的纵向测量方法 | |
Yu et al. | Adaptive sampling method for inspection planning on CMM for free-form surfaces | |
CN103453849B (zh) | 多光学传感器协同的复杂曲面零件三维测量方法与系统 | |
CN104574502A (zh) | 基于钢结构模型的激光断面特征识别方法 | |
CN110059879B (zh) | 车身三坐标测量的自动规划方法 | |
CN108645339A (zh) | 一种生物发电厂料垛点云数据采集及体积计算方法 | |
CN105066915A (zh) | 模具曲面加工误差和表面粗糙度在机检测装置及检测方法 | |
CN103884291B (zh) | 基于nurbs参数曲面的建筑物表面柔性变形监测方法 | |
JP6186305B2 (ja) | たわみ推定装置、及びプログラム | |
CN107622530B (zh) | 一种高效鲁棒的三角网切割方法 | |
CN103438826A (zh) | 激光与视觉相结合的钢板的三维测量系统及方法 | |
CN108332685B (zh) | 一种编码结构光三维测量方法 | |
CN104359415B (zh) | 一种水火弯板角变形量测量方法及系统 | |
CN110220476A (zh) | 一种基于三维激光扫描的航空零件自动快速测量方法 | |
JP2010286267A (ja) | 対象物検出装置及び対象物検出方法 | |
CN105387842A (zh) | 基于感知驱动的自航式海底地形地貌测绘系统及测绘方法 | |
CN109141266A (zh) | 一种钢结构测量方法及系统 | |
CN111540001A (zh) | 航空发动机涡轮叶片气膜孔轴线方向检测方法 | |
CN112733428A (zh) | 一种光学测量的扫描姿态与覆盖路径规划方法 | |
CN104596464B (zh) | 基于径向扫描结合环向扫描的摆臂轮廓面形检测方法 | |
CN107292039B (zh) | 一种基于小波聚类的uuv巡岸轮廓构建方法 | |
CN108446496A (zh) | 一种航海模拟器中实时岸线碰撞检测方法和装置 | |
CN105651202A (zh) | 一种用于测量矿山体积的三维扫描方法及装置 | |
CN106407603B (zh) | 一种虚拟齿轮测量中心上回转体工件快速碰撞检测的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |