CN103486984A - 一种风洞内型面同轴度的检测方法 - Google Patents
一种风洞内型面同轴度的检测方法 Download PDFInfo
- Publication number
- CN103486984A CN103486984A CN201310476290.2A CN201310476290A CN103486984A CN 103486984 A CN103486984 A CN 103486984A CN 201310476290 A CN201310476290 A CN 201310476290A CN 103486984 A CN103486984 A CN 103486984A
- Authority
- CN
- China
- Prior art keywords
- wind
- tunnel
- profile
- point
- coordinate
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 claims description 16
- 230000008602 contraction Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Excavating Of Shafts Or Tunnels (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种风洞内型面同轴度的检测方法,包括以下步骤:(a)用三维激光扫描仪扫描风洞内型面,获取风洞内型面点云;(b)将扫描所得到的点云数据的坐标统一到同一坐标系内;(c)沿风洞轴线方向对内型面的点云数据进行等间距剖切,在每个剖切面上可获得风洞内型面点云与剖切面的相交线;(d)将风洞内型面点云与剖切面的相交线上的点进行等间距处理,并计算每个剖切面上等间距处理后点集的质心点坐标;(e)根据内型面点云与剖切面相交线上等间距点集的质心点拟合出一条直线,该直线为基准轴线;(f)计算风洞内各站内型面点云与剖切面相交线上等间距点集的质心点到对应基准轴线的距离。其解决了风洞内型面同轴度误差的全面检测问题。
Description
技术领域
本发明涉及一种同轴度的检测方法,尤其涉及一种用于风洞内型面同轴度的检测方法。
背景技术
风洞内型面的质量和精度是影响风洞内模拟风场品质的重要因素,风洞各部段的型面同轴度偏差会影响空气流场的特性参数,导致风洞实验结果不准确,因此,在风洞工程建设中,对于风洞内型面的同轴度提出了很高的要求。目前,由于风洞内型面尺寸大,形状复杂多样,且传统测量仪器和测量方法只能对少数离散单点进行测量,无法获取型面的整体三维空间信息,不能对内型面同轴度进行有效的全面检测。
根据《产品几何技术规范(GPS)几何公差形状、方向、位置和跳动公差标注》(GB/T1182-2008)中关于同轴度误差的定义,同轴度是指被测圆柱面轴线对基准轴线不共轴的程度。根据定义,用与基准轴线同轴的圆柱面来包容被测物体的实测轴线,在被测长度内,最小包容圆柱面的直径即为同轴度误差值。风洞内型面的截面形状有矩形、八边形,与定义中所述圆柱面不同,且风洞尺寸较大,工业领域的同轴度检测方法无法应用于风洞内型面的检测。
近几年,三维激光扫描技术不断发展,已逐渐应用于地形测绘、基坑监测、工业测量等领域,成为继GPS技术之后的又一项测绘技术革命。三维激光扫描能快速获取目标物体表面大量采样点的空间位置坐标,得到一个描述三维实体的点集合,即点云。因此,三维激光扫描技术的出现和发展为解决风洞内型面同轴度检测问题提供了新的思路。
发明内容
为了解决上述问题,本发明提供了一种风洞内型面同轴度的检测方法,其利用三维激光扫描仪能够高效率、高精度获取风洞内型面大量空间点坐标信息的优势,将实际内型面的空间几何信息以点云数据的方式复制到计算机中,通过分析内型面点云数据解决了风洞内型面同轴度误差的高精度全面检测问题。
为了实现以上目的,本发明提供了一种风洞内型面同轴度的检测方法,风洞中设置有至少两个站,风洞内型面同轴度的检测方法包括以下步骤:(a)用三维激光扫描仪沿风洞轴线方向扫描风洞内型面,获取风洞内型面的点云;(b)将扫描所得到的点云数据的坐标通过坐标变换统一到同一坐标系内;(c)沿风洞轴线方向对内型面的点云数据进行等间距剖切,在每个剖切面上可获得风洞内型面点云与剖切面的相交线,该相交线由一系列离散线构成;(d)将风洞内型面点云与剖切面的相交线上的离散点进行等间距处理,并计算每个剖切面上等间距处理后的点集的质心点坐标;(e)根据内型面点云与剖切面相交线上等间距点集的质心点拟合出一条直线,该直线为基准轴线;(f)计算风洞内各站内型面点云与剖切面相交线上等间距点集的质心点到对应基准轴线的距离,其中最大距离值的两倍即为同轴度误差值。
优选地,步骤(a)具体过程为将三维激光扫描仪置于风洞内部,使三维激光扫描仪的激光发射器发射的激光扫描风洞内型面,获取风洞内型面的点云数据通过三维激光扫描仪接收器接收。
优选地,在步骤(a)中,测量风洞内型面的方法具体为在风洞中的相邻两站之间布设至少三个标靶球;每站可通视区域设置至少三个标靶板,并采用全站型电子速测仪测量标靶板中心的坐标。
优选地,在步骤(b)中统一点云数据坐标系的方法具体为:根据标靶球和标靶板的点云数据以及标靶板中心的坐标,将各站扫描所得点云数据的坐标通过点云数据配准和坐标变换统一到同一坐标系内。
优选地,在步骤(e)中,拟合基准轴线的具体方法为:对于风洞内与风洞收缩段同轴的区域,以收缩段区域的内型面点云与剖切面的相交线上等间距点集的质心点拟合出的直线为基准轴线;对于除与收缩段同轴的区域以外的其他区域,以同轴各区域内型面点云与剖切面相交线上等间距点集的质心点拟合出的直线为基准轴线。
本发明相较于现有技术所具有的优势在于:利用三维激光扫描仪能够高效率、高精度获取内型面大量空间点坐标信息的优势,将实际内型面的空间几何信息以点云数据的方式复制到计算机中,通过分析内型面点云数据解决了风洞内型面同轴度误差的高精度全面检测问题。该方法与传统检测方法相比,能获取风洞内型面的完整空间几何信息,检测结果精度高。
附图说明
图1为三维激光扫描仪扫描获取的风洞内型面点云变换到同一坐标系中的示意图;
图2为剖切后获得的内型面点云与剖切面的相交线示意图;
图3为内型面的点云与一个剖切面的相交线上的离散点集及该离散点集的质心点示意图。
具体实施方式
下面,通过示例性的实施方式对本发明进行具体描述。然而应当理解,在没有进一步叙述的情况下,一个实施方式中的元件、结构和特征也可以有益地结合到其他实施方式中。
本发明公开了一种风洞内型面同轴度的检测方法,风洞中设置有至少两个站,该站的划分是根据三维激光扫描仪的有效扫描范围以及风洞内可通视范围综合考虑划分的,其中风洞内型面同轴度的检测方法包括以下步骤:(a)用三维激光扫描仪沿风洞轴线方向扫描风洞内型面,获取风洞内型面的点云数据,具体过程为将三维激光扫描仪置于风洞内部,使三维激光扫描仪的激光发射器发射的激光沿风洞轴线方向扫描风洞内型面,获取风洞内型面的点云数据通过三维激光扫描仪接收器接收。其中优选地,测量风洞内型面的方法具体为在风洞中的相邻两站之间布设至少三个标靶球,标靶球的优选布设方法为标靶球尽量共面,但不要共线;每站可通视区域设置至少三个标靶板,并采用全站型电子速测仪测量标靶板中心的坐标,该全站型电子速测仪的测角精度优选为±0.5";
(b)将扫描所得到的点云数据的坐标通过坐标变换统一到同一坐标系内;上述优选测量风洞内型面方法的实施方式中,具体为:根据标靶球和标靶板的点云数据以及标靶板中心的坐标,将各站扫描所得点云数据的坐标通过点云数据配准和坐标变换统一到同一坐标系内;如图1所示,即为上述三维激光扫描仪扫描获取的风洞内型面点云变换到同一坐标系后的示意图。
(c)沿风洞轴线方向对内型面的点云数据进行等间距剖切,在每个剖切面上可获得风洞内型面点云与剖切面的相交线,如图2所示,该相交线由一系列离散线构成;
(d)通过插值计算,将风洞内型面点云与剖切面的相交线上的离散点进行等间距处理,并计算每个剖切面上等间距处理后的离散点集的质心点坐标。具体算法如下:设某一剖切面相交线上的离散线包含n个点,第i个点的坐标为(xi,yi,zi),质量为ρi,则该离散线内点集的质心坐标(xo,yo,zo)为:
由于点等间距分布,设各点的质量相同,由上式可得,质心点坐标即为各点坐标的平均值;如图3所示,即为内型面的点云与一个剖切面的相交线上的离散点集及该离散点集的质心点示意图,图中所示的中心点即为质心点。
(e)根据内型面点云与剖切面相交线上等间距点集的质心点拟合出一条直线,该直线为基准轴线,具体为:对于风洞内与收缩段同轴的区域,以收缩段区域的内型面点云与剖切面的相交线上等间距点集的质心点拟合出的直线为基准轴线,对于除与收缩段同轴的区域以外的其他区域,以同轴各区域内型面点云与剖切面相交线上等间距点集的质心点拟合出的直线为基准轴线;拟合基准轴线的具体方法如下:
首先根据各剖切面上求得的质心点的x坐标和y坐标采用最小二乘法拟合出空间直线在xy平面上的投影方程:
a1x+b1y+d1=0
根据x坐标和z坐标拟合出空间直线在xz平面上的投影方程:
a2x+c2z+d2=0
则空间直线即轴线的一般方程为:
由一般方程可求出参数式方程:
式中,t为参数,(x0,y0,z0)为直线上任取一点的坐标;
(f)计算各站内型面点云与剖切面相交线上等间距点集的质心点到对应基准轴线的距离,其中最大距离值的两倍即为同轴度误差值。具体方法如下:
设某一剖切平面上的质心点坐标为(xj,yj,zj),则过该点与基准轴线垂直的平面方程为:
b1c2(x-xj)-a1c2(y-yj)-a2b1(z-zj)=0
将基准轴线的参数方程代入平面方程可求得基准轴线与垂直平面的交点坐标(xk,yk,zk),利用空间两点的距离公式即可求出点(xj,yj,zj)与点(xk,yk,zk)的距离,即点(xj,yj,zj)到基准轴线的距离。按照该方法求出同轴各段所有剖切面上的质心点到基准轴线的距离,其中距离最大值dmax即为与基准轴线同轴且包容所有质心点(即实测轴线)的最小包容圆柱面的半径,直径2dmax即为同轴度偏差值。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
Claims (5)
1.一种风洞内型面同轴度的检测方法,风洞中设置有至少两个站,其特征在于,风洞内型面同轴度的检测方法包括以下步骤:
(a)用三维激光扫描仪沿风洞轴线方向扫描风洞内型面,获取风洞内型面的点云;
(b)将扫描所得到的点云数据的坐标通过坐标变换统一到同一坐标系内;
(c)沿风洞轴线方向对内型面的点云数据进行等间距剖切,在每个剖切面上可获得风洞内型面点云与剖切面的相交线,该相交线由一系列离散点构成;
(d)将风洞内型面点云与剖切面的相交线上的离散点进行等间距处理,并计算每个剖切面上等间距处理后点集的质心点坐标;
(e)根据内型面点云与剖切面相交线上等间距点集的质心点拟合出一条直线,该直线为基准轴线;
(f)根据插值计算,计算风洞内各站内型面点云与剖切面相交线上等间距点集的质心点到对应基准轴线的距离,其中最大距离值的两倍即为同轴度误差值。
2.根据权利要求1所述的风洞内型面同轴度的检测方法,其特征在于,步骤(a)具体过程为将三维激光扫描仪置于风洞内部,使三维激光扫描仪的激光发射器发射的激光沿风洞轴线方向扫描风洞内型面,获取风洞内型面的点云数据通过三维激光扫描仪接收器接收。
3.根据权利要求1或2所述的风洞内型面同轴度的检测方法,其特征在于,在步骤(a)中,测量风洞内型面的方法具体为在风洞中的相邻两站之间布设至少三个标靶球;每站可通视区域设置至少三个标靶板,并采用全站型电子速测仪测量标靶板中心的坐标。
4.根据权利要求3所述的风洞内型面同轴度的检测方法,其特征在于,在步骤(b)中统一点云数据坐标系的方法具体为:根据标靶球和标靶板的点云数据以及标靶板中心的坐标,将各站扫描所得点云数据的坐标通过点云数据配准和坐标变换统一到同一坐标系内。
5.根据权利要求1所述的风洞内型面同轴度的检测方法,其特征在于,在步骤(e)中,拟合基准轴线的具体方法为:对于风洞内与收缩段同轴的区域,以收缩段区域的内型面点云与剖切面的相交线上等间距点集的质心点拟合出的直线为基准轴线,对于除与收缩段同轴的区域以外的其他区域,以同轴各区域内型面点云与剖切面相交线上等间距点集的质心点拟合出的直线为基准轴线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310476290.2A CN103486984B (zh) | 2013-10-12 | 2013-10-12 | 一种风洞内型面同轴度的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310476290.2A CN103486984B (zh) | 2013-10-12 | 2013-10-12 | 一种风洞内型面同轴度的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103486984A true CN103486984A (zh) | 2014-01-01 |
CN103486984B CN103486984B (zh) | 2015-11-25 |
Family
ID=49827400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310476290.2A Expired - Fee Related CN103486984B (zh) | 2013-10-12 | 2013-10-12 | 一种风洞内型面同轴度的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103486984B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104374224A (zh) * | 2014-11-19 | 2015-02-25 | 金龙精密铜管集团股份有限公司 | 强化蒸发传热管 |
CN105910814A (zh) * | 2016-07-01 | 2016-08-31 | 成都市龙泉通惠实业有限责任公司 | 活塞体通孔检测装置 |
CN106392609A (zh) * | 2016-08-30 | 2017-02-15 | 武汉冶钢结构有限责任公司 | 大型风洞安装同轴度调整辅助装置及其调整方法 |
CN106441165A (zh) * | 2016-08-30 | 2017-02-22 | 哈尔滨锅炉厂有限责任公司 | 一种补偿计算容器同轴度的检测方法 |
CN107167090A (zh) * | 2017-03-13 | 2017-09-15 | 深圳市速腾聚创科技有限公司 | 车辆外廓尺寸测量方法及系统 |
CN107270881A (zh) * | 2017-07-12 | 2017-10-20 | 中国人民解放军63926部队 | 一种风洞体内型面施工控制点高效测设方法 |
CN113418495A (zh) * | 2021-08-03 | 2021-09-21 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种风洞喷管同轴度测量方法 |
CN114842039A (zh) * | 2022-04-11 | 2022-08-02 | 中国工程物理研究院机械制造工艺研究所 | 一种含回转体微结构金刚石压砧的同轴度误差计算方法 |
CN117168365A (zh) * | 2023-09-05 | 2023-12-05 | 天津大学 | 一种电机定子同轴度测量方法及装置 |
CN118089598A (zh) * | 2024-04-17 | 2024-05-28 | 四川吉埃智能科技有限公司 | 一种基于三维检测的管道接长方法、装置、设备及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2618972A (en) * | 1948-03-12 | 1952-11-25 | Leduc Rene | Testing apparatus for scale models of aerodynamic ducts |
CN101825444A (zh) * | 2010-04-09 | 2010-09-08 | 上海辉格科技发展有限公司 | 基于面结构光的车载路谱测试系统 |
CN102393183A (zh) * | 2011-04-19 | 2012-03-28 | 程效军 | 基于控制网的海量点云快速配准方法 |
CN102446354A (zh) * | 2011-08-29 | 2012-05-09 | 北京建筑工程学院 | 一种高精度多源地面激光点云的整体配准方法 |
CN102692190A (zh) * | 2011-03-24 | 2012-09-26 | 北京林业大学 | 一种基于数字高程模型的树木体积计测方法 |
CN102798412A (zh) * | 2012-07-31 | 2012-11-28 | 同济大学 | 一种基于三维激光扫描评定隧道钻爆施工质量的方法 |
-
2013
- 2013-10-12 CN CN201310476290.2A patent/CN103486984B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2618972A (en) * | 1948-03-12 | 1952-11-25 | Leduc Rene | Testing apparatus for scale models of aerodynamic ducts |
CN101825444A (zh) * | 2010-04-09 | 2010-09-08 | 上海辉格科技发展有限公司 | 基于面结构光的车载路谱测试系统 |
CN102692190A (zh) * | 2011-03-24 | 2012-09-26 | 北京林业大学 | 一种基于数字高程模型的树木体积计测方法 |
CN102393183A (zh) * | 2011-04-19 | 2012-03-28 | 程效军 | 基于控制网的海量点云快速配准方法 |
CN102446354A (zh) * | 2011-08-29 | 2012-05-09 | 北京建筑工程学院 | 一种高精度多源地面激光点云的整体配准方法 |
CN102798412A (zh) * | 2012-07-31 | 2012-11-28 | 同济大学 | 一种基于三维激光扫描评定隧道钻爆施工质量的方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104374224A (zh) * | 2014-11-19 | 2015-02-25 | 金龙精密铜管集团股份有限公司 | 强化蒸发传热管 |
CN105910814A (zh) * | 2016-07-01 | 2016-08-31 | 成都市龙泉通惠实业有限责任公司 | 活塞体通孔检测装置 |
CN106392609A (zh) * | 2016-08-30 | 2017-02-15 | 武汉冶钢结构有限责任公司 | 大型风洞安装同轴度调整辅助装置及其调整方法 |
CN106441165A (zh) * | 2016-08-30 | 2017-02-22 | 哈尔滨锅炉厂有限责任公司 | 一种补偿计算容器同轴度的检测方法 |
CN107167090A (zh) * | 2017-03-13 | 2017-09-15 | 深圳市速腾聚创科技有限公司 | 车辆外廓尺寸测量方法及系统 |
CN107270881A (zh) * | 2017-07-12 | 2017-10-20 | 中国人民解放军63926部队 | 一种风洞体内型面施工控制点高效测设方法 |
CN107270881B (zh) * | 2017-07-12 | 2019-07-05 | 中国人民解放军63926部队 | 一种风洞体内型面施工控制点高效测设方法 |
CN113418495A (zh) * | 2021-08-03 | 2021-09-21 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种风洞喷管同轴度测量方法 |
CN114842039A (zh) * | 2022-04-11 | 2022-08-02 | 中国工程物理研究院机械制造工艺研究所 | 一种含回转体微结构金刚石压砧的同轴度误差计算方法 |
CN117168365A (zh) * | 2023-09-05 | 2023-12-05 | 天津大学 | 一种电机定子同轴度测量方法及装置 |
CN118089598A (zh) * | 2024-04-17 | 2024-05-28 | 四川吉埃智能科技有限公司 | 一种基于三维检测的管道接长方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103486984B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103486984B (zh) | 一种风洞内型面同轴度的检测方法 | |
CN107167790B (zh) | 一种基于标定场的激光雷达两步标定方法 | |
CN103605135B (zh) | 一种基于断面剖分的道路特征提取方法 | |
CN103453849B (zh) | 多光学传感器协同的复杂曲面零件三维测量方法与系统 | |
CN106597417A (zh) | 一种远距离扫描激光雷达测量误差的修正方法 | |
CN103542820B (zh) | 一种检测风洞内表面平整度的方法 | |
CN105136054A (zh) | 基于地面三维激光扫描的构筑物精细变形监测方法及系统 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN109141266B (zh) | 一种钢结构测量方法及系统 | |
CN108871268B (zh) | 一种基于激光点云的隧道超欠挖数值计算方法 | |
CN104007444A (zh) | 一种基于中心投影的地面激光雷达反射强度图像的生成方法 | |
CN108120394A (zh) | 柔性化的高速列车车头曲面质量检测方法 | |
CN102749623B (zh) | 一种基于靶标的高精度遥感卫星地面采样距离测试方法 | |
CN106568365A (zh) | 一种球面孔系复合位置度误差的检测与评定方法 | |
CN112857315B (zh) | 基于三维激光扫描的钢管柱垂直度测量方法 | |
CN111765902A (zh) | 一种基于多棱锥标靶的激光点云精度评估方法 | |
CN104729529A (zh) | 地形图测量系统误差判断的方法和系统 | |
CN106646413B (zh) | 一种雷达组网垂线交叉融合定位方法及误差解算方法 | |
CN107389293B (zh) | 一种检测风洞收缩段内型面曲面误差的方法 | |
CN107424185A (zh) | 基于点云数据处理技术的圆锥体结构特征参数检测方法 | |
Pan et al. | Assessment method of slope excavation quality based on point cloud data | |
CN106197274A (zh) | 一种数字化检测具有曲面的模具的刻线位置度的方法 | |
CN108332784B (zh) | 一种距离测量验证方法 | |
CN106323587A (zh) | 机翼风洞试验模型弹性变形的单目视频高精度测量方法 | |
CN106772332B (zh) | 一种旋转式距离测量验证方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 |