CN107607080B - 一种棱柱状管道横截面形变测量计算方法 - Google Patents
一种棱柱状管道横截面形变测量计算方法 Download PDFInfo
- Publication number
- CN107607080B CN107607080B CN201710609214.2A CN201710609214A CN107607080B CN 107607080 B CN107607080 B CN 107607080B CN 201710609214 A CN201710609214 A CN 201710609214A CN 107607080 B CN107607080 B CN 107607080B
- Authority
- CN
- China
- Prior art keywords
- section
- cross
- profile
- mass center
- section profile
- 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.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种棱柱状管道横截面形变测量计算方法，根据棱柱状管道横截面轮廓点的坐标测量值，对此截面轮廓点的数据进行处理，拟合出此截面的轮廓，根据拟合出的轮廓与基准截面轮廓上某些特征的相对关系来反映截面的真实变形情况，对截面的变形进行了详细的分析计算，精确的描述截面的变形情况。本发明创新的提出了一种棱柱状管道截面形变的测量计算方法，根据基准轮廓与拟合出的轮廓的差异，提出一些算法，通过算法来判定截面的变形情况。本发明是一种高精度的棱柱状横截面管道形变测量计算方法，可推广应用到对类似的棱柱状棒料或工件截面形变测量的测量计算。
Description
技术领域
本发明涉及一种管道横截面形变测量计算领域，尤其是涉及一种棱柱状管道横截面形变测量计算方法。
背景技术
在实际的检测过程中，对棱柱状管道截面形变的测量计算方法较少，在管道或易变性工件的检测过程中，现有的棱柱状管道截面形变测量计算方法有测量实际截面的对边距与未变形时的对边距的差值、实际截面的对角距与未变形时的对角距的差值两种，这两种测量计算方法不精确，不能真实反映截面的变形情况，测量人员对于测量结果较模糊，根据上述两种测量计算方法得出的结论会对测量人员造成误导，对于高精度管道或工件的形变，迫切需要一种新的测量计算方法，真实全面的反映截面形变，为管道或工件的维护和更换提供依据。
发明内容
本发明要解决的问题是，针对现有技术的不足，提供一种棱柱状管道横截面形变测量计算方法。
为了解决上述问题，本发明采用的技术方案是：利用测量装置测出棱柱状管道某一截面轮廓上各测点的位置数据，根据棱柱状管道横截面轮廓点的坐标测量值，对此截面轮廓点的数据进行处理，拟合出此截面的轮廓，根据拟合出的轮廓与基准截面轮廓上某些特征的相对关系来反映截面的真实变形情况，对截面的变形进行了详细的分析和计算，精确的描述截面的变形情况。
所述一种棱柱状管道横截面形变测量计算方法，包括以下步骤：
1.利用测量装置测出棱柱状管道某一横截面轮廓点的坐标测量值，对这些坐标测量值进行处理，拟合出此截面的轮廓；
2.根据拟合出的截面轮廓，求出轮廓的质心以轮廓的质心为圆心，作轮廓的最小外接圆和最大内切圆，最小外接圆半径R2，最大内切圆半径R1，通过两圆的半径差来计算截面形变；
3.以基准截面轮廓质心O(0，0)为圆心，作拟合截面轮廓的最小外接圆和最大内切圆，最小外接圆半径R22，最大内切圆半径R11，通过两圆的半径差来计算截面形变；
4.求出拟合截面轮廓的质心与基准截面轮廓质心O(0，0)的距离Rj，通过两质心的距离Rj和拟合截面轮廓的质心和基准截面轮廓的质心O(0，0)连线与X轴正方向夹角为ωj来计算质心偏移距离和方向；
5.求出拟合截面轮廓各个顶点和拟合截面轮廓质心连线与X轴的夹角θ1i，基准截面轮廓的相应顶点和基准截面轮廓质心连线与X轴的夹角θ2i，求出两者的差值θi，用θi作为扭转计算参数AA，θi中的最大值max(θi)来计算横截面的扭转程度；
6.求出拟合截面轮廓各个顶点与拟合截面轮廓质心连线和基准截面轮廓的相应顶点与拟合截面轮廓质心连线的角度差αi，用αi作为扭转计算参数BB，用αi的最大值max(αi)来计算横截面的扭转程度；
7.求出拟合截面轮廓各个顶点与基准截面轮廓质心连线和基准截面轮廓的相应顶点与基准截面轮廓质心连线的角度差γi，用γi作为扭转计算参数C，以γi的最大值max(γi)来计算横截面的扭转程度；
8.求出拟合截面轮廓的各个边与X轴的夹角β1i，基准截面轮廓的相应边与X轴的夹角β2i，两者的差βi，用βi作为扭转计算参数D，以βi的最大值max(βi)来计算横截面边的扭转程度。
收益效果：与现有技术相比，本发明所具有的有益效果为：本发明提出的一种棱柱状管道截面形变测量计算方法，利用测量装置测出棱柱状管道某一截面轮廓上各测点的位置数据，根据棱柱状管道横截面各测点的位置数据，对此截面轮廓点的数据进行处理，拟合出此截面的轮廓，根据拟合出的截面轮廓与基准截面轮廓上某些特征的相对关系来反映截面的真实变形情况，对截面的变形进行了详细的分析和计算，精确、全面的描述截面的变形情况。
附图说明
图1是本发明拟合截面轮廓质心的最小外接圆和最大内切圆示意图；
图2是本发明基准截面轮廓质心最小外接圆和最大内切圆示意图；
图3是本发明拟合截面轮廓质心与基准截面轮廓质心的相对位置示意图；
图4是本发明扭转计算参数AA示意图。
图5是本发明扭转计算参数BB示意图。
图6是本发明扭转计算参数C示意图。
图7是本发明扭转计算参数D示意图。
具体实施方式
下面结合附图对本发明的实施例作详细说明，本实施例以本发明技术方案为前提，以棱柱状管道截面的形变测量计算方法为例，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如图1所示，测量装置测出多边形管道某一截面轮廓上各测点的位置数据，用数学算法对此截面各边的位置数据进行处理，得到数据处理完后拟合截面轮廓，计算出截面的实际质心点和六边形顶点坐标最小外接圆半径R2，最大内切圆半径R1，半径差值其中i为同一截面内顶点或边序号，i＝1,2,3,4,5,6。j为检测截面序号，j＝1,2,3…n,n为检测截面个数,其中：
如图2所示，六边形顶点坐标未变形截面的质心点O(0，0),最小外接圆半径R22，最大内切圆半径R11，半径差值其中i为同一截面内顶点或边序号，i＝1,2,3,4,5,6。j为检测截面序号，j＝1,2,3…n,n为检测截面个数,其中：
如图3所示，实际截面多边形的质心点未变形截面的质心点O(0，0)，两质心的距离Rj，实际截面六边形的质心和未变形截面的质心连线与X轴正方向夹角为ωj，其中i为同一截面内顶点或边序号，i＝1,2,3,4,5,6。j为检测截面序号，j＝1,2,3…n,n为检测截面个数,其中：
如图4所示，求出拟合截面轮廓某个顶点和拟合截面轮廓质心连线与X轴的夹角θ1i，基准截面轮廓的相应顶点和基准截面轮廓质心连线与X轴的夹角θ2i，求出两者的差值θi，用θi中的最大值max(θi)来计算横截面的扭转程度，其中i为截面内顶点或边序号，i＝1,2,3…6。
如图5所示，求出拟合截面轮廓各个顶点与拟合截面轮廓质心连线和基准截面轮廓的相应顶点与拟合截面轮廓质心连线的角度差的最大值max(αi)来计算横截面的扭转程度，其中i为截面内顶点或边序号，i＝1,2,3…6。
如图6所示，求出拟合截面轮廓各个顶点与基准截面轮廓质心连线和基准截面轮廓的相应顶点与基准截面轮廓质心连线的角度差max(γi)来计算横截面的扭转程度，其中i为截面内顶点或边序号，i＝1,2,3…6。
如图7所示，拟合截面轮廓的一边与X轴的夹角为β1i，基准截面轮廓对应边与X轴的夹角为β2i，两者的差值为βi，以βi的最大值max(βi)来计算横截面的扭转程度，其中：
βi＝β1-β2
本发明未叙述部分为现有技术。
Claims (1)
1.一种棱柱状管道截面形变的测量计算方法，其特征在于：所述测量计算方法包含三个截面变形计算参数和四个截面扭转变形计算参数，包括以下步骤：
1)利用测量装置测出棱柱状管道某一横截面轮廓点的坐标测量值，对这些坐标测量值进行处理，拟合出此截面的轮廓；
2)根据拟合出的截面轮廓，求出轮廓的质心以轮廓的质心为圆心，作截面轮廓多边形的最小外接圆和最大内切圆，最小外接圆半径R2，最大内切圆半径R1，通过两圆的半径差来计算截面形变，其中为检测截面拟合截面轮廓多边形的顶点坐标，为检测截面拟合截面轮廓多边形相邻顶点连成直线的斜率，为检测截面拟合截面轮廓多边形相邻顶点连成直线的截距，其中i为截面内顶点或边序号，i＝1,2,3…m，m为截面内顶点个数，j为检测截面序号，j＝1,2,3…n,n为检测截面个数，其中：
3)以基准截面轮廓质心O(0，0)为圆心，作截面轮廓多边形的最小外接圆和最大内切圆，最小外接圆半径R22，最大内切圆半径R11，通过两圆的半径差来计算截面形变，其中i为截面内顶点或边序号，i＝1,2,3…m，m为截面内顶点个数，j为检测截面序号，j＝1,2,3…n,n为检测截面个数，其中：
4)求出拟合截面轮廓的质心与基准截面轮廓质心O(0，0)的距离Rj，通过两质心的距离Rj和拟合截面轮廓的质心和基准截面轮廓的质心O(0，0)连线与X轴正方向夹角为ωj来计算质心偏移距离和方向，其中i为截面内顶点或边序号，i＝1,2,3…m，m为截面内顶点个数，j为检测截面序号，j＝1,2,3…n,n为检测截面个数，其中：
5)求出拟合截面轮廓各个顶点和拟合截面轮廓质心连线与X轴的夹角θ1i，基准截面轮廓的相应顶点和基准截面轮廓质心连线与X轴的夹角θ2i，求出两者的差值θi，用θi作为扭转计算参数AA，θi中的最大值max(θi)来计算横截面的扭转程度，其中i为截面内顶点序号，i＝1,2,3…m，m为截面内顶点个数；
6)求出拟合截面轮廓各个顶点与拟合截面轮廓质心连线和基准截面轮廓的相应顶点与拟合截面轮廓质心连线的角度差αi，用αi作为扭转计算参数BB，用αi的最大值max(αi)来计算横截面的扭转程度，其中i为截面内顶点序号，i＝1,2,3…m，m为截面内顶点个数；
7)求出拟合截面轮廓各个顶点与基准截面轮廓质心连线和基准截面轮廓的相应顶点与基准截面轮廓质心连线的角度差γi，用γi作为扭转计算参数C，以γi的最大值max(γi)来计算横截面的扭转程度，其中i为截面内顶点序号，i＝1,2,3…m，m为截面内顶点个数；
8)求出拟合截面轮廓的各个边与X轴的夹角β1i，基准截面轮廓的相应边与X轴的夹角β2i，两者的差βi，用βi作为扭转计算参数D，以βi的最大值max(βi)来计算横截面边的扭转程度，其中i为截面内边序号，i＝1,2,3…m，m为截面边或顶点个数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710609214.2A CN107607080B (zh) | 2017-07-24 | 2017-07-24 | 一种棱柱状管道横截面形变测量计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710609214.2A CN107607080B (zh) | 2017-07-24 | 2017-07-24 | 一种棱柱状管道横截面形变测量计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107607080A CN107607080A (zh) | 2018-01-19 |
| CN107607080B true CN107607080B (zh) | 2019-05-28 |
Family
ID=61059551
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710609214.2A Expired - Fee Related CN107607080B (zh) | 2017-07-24 | 2017-07-24 | 一种棱柱状管道横截面形变测量计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107607080B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108444431B (zh) * | 2018-03-19 | 2020-02-18 | 成都飞机工业（集团）有限责任公司 | 基于三个安装点的管件外形检测方法 |
| CN109373961B (zh) * | 2018-12-25 | 2020-03-24 | 浙江大学 | 一种六方轴直线度检测装置与数据处理方法 |
| CN109916366B (zh) * | 2019-03-25 | 2021-04-20 | 西安电子科技大学 | 一种实时监测管道形变与姿态的系统及方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08285572A (ja) * | 1995-04-14 | 1996-11-01 | Nippon Steel Corp | 薄板形状の評価方法 |
| JP2000002527A (ja) * | 1998-06-17 | 2000-01-07 | Daiwa Seiko Inc | 管状体のたわみ特性の指標を表示する測定方法、及びその測定方法を指標とした管状体、並びにそのたわみ特性を測定する測定装置 |
| WO2009042997A2 (en) * | 2007-09-28 | 2009-04-02 | Lam Research Corporation | Wafer bow metrology arrangements and methods thereof |
| CN101536016A (zh) * | 2006-09-18 | 2009-09-16 | 华盛顿大学 | 光学微层析成像的焦平面跟踪 |
| CN103018024A (zh) * | 2012-12-03 | 2013-04-03 | 奇瑞汽车股份有限公司 | 一种扭转刚度的评价方法 |
| CN204495559U (zh) * | 2015-04-07 | 2015-07-22 | 中国直升机设计研究所 | 一种弹性轴承扭转刚度测试装置 |
| CN106840022A (zh) * | 2016-07-25 | 2017-06-13 | 中国科学院长春光学精密机械与物理研究所 | 一种复杂曲面面形误差评价方法 |
-
2017
- 2017-07-24 CN CN201710609214.2A patent/CN107607080B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08285572A (ja) * | 1995-04-14 | 1996-11-01 | Nippon Steel Corp | 薄板形状の評価方法 |
| JP2000002527A (ja) * | 1998-06-17 | 2000-01-07 | Daiwa Seiko Inc | 管状体のたわみ特性の指標を表示する測定方法、及びその測定方法を指標とした管状体、並びにそのたわみ特性を測定する測定装置 |
| CN101536016A (zh) * | 2006-09-18 | 2009-09-16 | 华盛顿大学 | 光学微层析成像的焦平面跟踪 |
| WO2009042997A2 (en) * | 2007-09-28 | 2009-04-02 | Lam Research Corporation | Wafer bow metrology arrangements and methods thereof |
| CN103018024A (zh) * | 2012-12-03 | 2013-04-03 | 奇瑞汽车股份有限公司 | 一种扭转刚度的评价方法 |
| CN204495559U (zh) * | 2015-04-07 | 2015-07-22 | 中国直升机设计研究所 | 一种弹性轴承扭转刚度测试装置 |
| CN106840022A (zh) * | 2016-07-25 | 2017-06-13 | 中国科学院长春光学精密机械与物理研究所 | 一种复杂曲面面形误差评价方法 |
Non-Patent Citations (3)
| Title |
|---|
| Assessment of the integrity of pipes containing circumferential defects using the method of real elements;Markotchev, VM et al;《INTERNATIONAL JOURNAL OF PRESSURE VESSELS AND PIPING》;19991130;第76卷(第14-15期);第935-943页 |
| 管道弯曲塑性屈曲应变的评价;郑茂盛等;《焊管》;20151031;第38卷(第10期);第19-27、31页 |
| 高速卧式加工中心床身结构模态优化设计研究;周军等;《机械制造》;20110331;第49卷(第559期);第37-41页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107607080A (zh) | 2018-01-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107607080B (zh) | 一种棱柱状管道横截面形变测量计算方法 | |
| CN103886593A (zh) | 一种基于三维点云曲面圆孔检测方法 | |
| CN103292701A (zh) | 基于机器视觉的精密器件在线尺寸测量方法 | |
| CN102735204B (zh) | 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法 | |
| CN101196394A (zh) | 小段圆弧圆度的优化最小二乘评价方法 | |
| CN104050372A (zh) | 一种三维几何形状误差自动评定方法 | |
| TWI724696B (zh) | 工件孔洞量測方法 | |
| JP2015232513A (ja) | たわみ推定装置、及びプログラム | |
| CN107091627A (zh) | 一种多表面系统的综合测量与评估方法 | |
| CN110806571A (zh) | 一种多结构光传感器空间姿态标定件及其标定方法 | |
| Yang et al. | An accurate and fast method to inspect rail wear based on revised global registration | |
| CN102980532B (zh) | 采用三坐标测量仪拼接测量大口径非球面面形的方法 | |
| CN112017293A (zh) | 一种关于圆钢管几何初始缺陷的测量方法 | |
| CN103292760A (zh) | 一种薄壁叶片误差分析方法 | |
| CN107504915B (zh) | 埋地管道变形程度检测方法及埋地管道变形程度评价方法 | |
| CN105627935A (zh) | 产品料厚三维分析系统及方法 | |
| CN113689394A (zh) | 一种基于三维点云的隧道初支平整度检测方法 | |
| CN101013022B (zh) | 一种消除光学元件干涉采样数据中随机误差的方法 | |
| CN114781056B (zh) | 一种基于特征匹配的飞机整机外形测量方法 | |
| CN106197274A (zh) | 一种数字化检测具有曲面的模具的刻线位置度的方法 | |
| CN113240637B (zh) | 一种基于机器学习的墙面平整度信息化检测方法和系统 | |
| CN106546189A (zh) | 一种系列尺寸岩体结构面起伏幅度各向异性的评价方法 | |
| JP4013986B1 (ja) | 固定構造体の曲げ応力測定方法、記録媒体及びコンピュータ | |
| CN102728658B (zh) | 一种用于确定叶片加工弯曲度误差的方法 | |
| CN105004293A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190528 Termination date: 20210724 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |