CN102607484B - 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 - Google Patents
一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 Download PDFInfo
- Publication number
- CN102607484B CN102607484B CN 201210083667 CN201210083667A CN102607484B CN 102607484 B CN102607484 B CN 102607484B CN 201210083667 CN201210083667 CN 201210083667 CN 201210083667 A CN201210083667 A CN 201210083667A CN 102607484 B CN102607484 B CN 102607484B
- Authority
- CN
- China
- Prior art keywords
- point
- sigma
- curve
- delta
- segment
- 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
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明提出一种变间距自适应搜寻测量点的飞机壁板法矢确定方法,依据壁板理论数学模型利用等曲率积分法自适应搜寻测量点,对测量点插值出两条三次B样曲线以确定曲面法矢。与现有方法相比,本发明方法利用了壁板曲面的理论数学模型提供的有益数学参数,采用等曲率积分法搜寻测量点可保证随曲线曲率的增大密化测量点,随曲线曲率的减小疏化测量点,通过优化测量点位的分布提高了壁板曲面的法矢求解精度,求解结果可用于引导自动钻铆设备对壁板进行制孔与铆接,并可用于分析估计壁板变形,作为检测壁板外形的依据,对于保障飞机壁板加工质量有重要意义。
Description
技术领域
本发明涉及测量检测技术领域,具体为一种变间距自适应搜寻测量点的飞机壁板法矢确定方法。
背景技术
自动钻铆技术主要应用于飞机壁板类部件的加工。其中,壁板钻孔法向精度直接影响飞机的气动外形乃至部件结构的安全性能,其技术关键在于如何利用壁板上有限个测量点确定钻孔点处的法矢。现有的法矢求解方法一般是在待测点周围选择一些临近测量点,然后根据这些点的测量值在其附近的拓扑形状采用平面、曲线或曲面拟合来计算,如Ruey-Tsung Lee等在文献Calculation of the unit normal vector using thecross-curve moving mask method for probe radius compensation of a freeform surfacemeasurement(Measurement,pp469-478,2010(43))中提出的9点3×3非标准二次贝齐尔曲面法及5点、9点非标准二次贝齐尔曲线插值法,Daoshan OuYang等在文献On thenormal vector estimation for point cloud data from smooth surface(Computer-AidedDesign,pp1071-1079,2005(37))中提出的点云估方法,易传云等在文献“数字化曲面的法矢求解”(华中科技大学学报,pp49-51,2002,30(8))中提出的累加弦长三次参数样条法等,这些曲线、曲面拟合法可以获得比较好的拟合效果,对任意空间曲面法矢求解精度高,普适性较强,但对壁板类曲面法矢的求解计算较为复杂,且上述方法求解法矢需要获得实际曲面上的大量测量点。求解效率方面,如秦现生等在文献“大型壁板数控钻铆的三点快速调平算法”(航空学报,pp1455-1460,2007,28(6))中提出的3点快速调平法数据采集量小,计算效率高,适用于求解大尺寸小曲率曲面法矢,但求解大曲度壁板法矢精度偏低。
此外,上述方法测量点的选取固定,测量点间距不随壁板曲面的变化而改变,并且测量点的选取规则往往依据经验,与壁板曲面的理论数学模型脱离联系,使理论数学模型在引导测量曲面法矢方面失去应用作用。虽然在实际生产过程中壁板变形导致其实际外形与理论数学模型存在偏差,但仍可将理论数学模型包含的相关数学特征信息作为搜寻测量点的依据,提高法矢求解精度。
发明内容
要解决的技术问题
为解决现有技术存在的问题,本发明提出一种变间距自适应搜寻测量点的飞机壁板法矢确定方法,依据壁板理论数学模型利用等曲率积分法自适应搜寻测量点,对测量点插值出两条三次B样曲线以确定曲面法矢。
技术方案
本发明的技术方案为:
所述一种变间距自适应搜寻测量点的飞机壁板法矢确定方法,其特征在于:包括以下步骤:
步骤1:根据壁板理论数学模型提取待测点P0的空间坐标、P0点所在壁板表面的曲面方程r(u,w)、以及过P0点的展向曲线方程rp(u)和过P0点的弦向曲线方程rp(w);
步骤2:以P0点为起点,沿P0周围的M、N、E、F四个方向按照步长h分别进行采样,其中M、N向为P0两侧沿着展向向外的方向,E、F向为P0两侧沿着弦向向外的方向;采样得到有序点集其中n为采样点数,且60mm≤nh≤100mm;
其中,ΔMi-1MiMi+1为三角形Mi-1MiMi+1的面积,ΔNi-1NiNi+1为三角形Ni-1NiNi+1的面积,ΔEi-1EiEi+1为三角形Ei-1EiEi+1的面积,ΔFi-1FiFi+1为三角形Fi-1FiFi+1的面积;
步骤5:取曲率积分值ΔA=0.1ω,ω为区域影响因子,其中当壁板处于机身、机翼区域时,ω取0.3~1,当壁板处于机翼前缘、尾翼前缘、整流罩、翼身融合处区域时,ω取1~2.5;
计算得到测量点在步骤2中采样得到的有序点集中的点位I1,I2,I3,I4,I5,I6,则得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
计算得到测量点在步骤2中采样得到的有序点集中的点位I1,I2,I3,I4,I5,I6,则得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
步骤6:采用测量仪器测量P0点以及步骤5中得到的六个测量点对应实际壁板上的空间位置,得到P0点以及六个测量点的实际测量值:P′0,P′1,P′2,P′3,P′4,P′5,P′6;
分别计算P′0点处沿r′p(u)和r′p(w)两个方向的切矢P′(u)和P′(w),得到P′0点处法矢为:
有益效果
与现有方法相比,本发明方法利用了壁板曲面的理论数学模型提供的有益数学参数,采用等曲率积分法搜寻测量点可保证随曲线曲率的增大密化测量点,随曲线曲率的减小疏化测量点,通过优化测量点位的分布提高了壁板曲面的法矢求解精度,求解结果可用于引导自动钻铆设备对壁板进行制孔与铆接,并可用于分析估计壁板变形,作为检测壁板外形的依据,对于保障飞机壁板加工质量有重要意义。以某型飞机机身及机翼壁板曲面作为分析对象实例验证表明,采样本发明方法求解法矢有如下优点:
1、法矢求解精度较高,选取相同的测量点数时,对于计算小曲率平缓区域(如机身)上的待测点,等间距、变间距选取测量点两种方法求解精度相近;但对于计算大曲率区域(如机翼前缘、整流罩处)上的待测点,本发明方法精度降低幅度小,求解精度提高30%-55%,完全满足实际规定的精度要求。
2、计算效率较高,本发明方法在配置Intel(R)Core(TM)2处理器2.93GHz、内存2GB的PC机上运用MATLAB编程仿真时,平均计算时间为0.828s,满足实际使用要求。
附图说明
图1:本发明方法中机翼壁板待测钻孔点及参数曲线示意图;
图2:本发明方法中的对参数曲线均匀弧长采样原理图;
图3:本发明方法中的等曲率积分法搜寻测量点原理图;
图4:本发明方法中的3次B样条插值求解法矢原理图;
具体实施方式
下面结合具体实施例描述本发明:
本实施例为测量某型飞机机翼前缘壁板上待测钻铆点P0的法向矢量,由于现代飞机壁板多为自由复杂曲面,为了既迅速快捷又精确地测量并计算待测点处的法矢,本实施例中采用本发明提出的方法来确定壁板曲面上待测点的法矢,包括以下步骤:
步骤1:如附图1所示,在CATIA V5三维环境中根据壁板理论数学模型提取待测点P0处空间坐标、P0点所在壁板表面的曲面方程r(u,w)、以及过P0点的展向曲线方程rp(u)和过P0点的弦向曲线方程rp(w)。
步骤2:如附图2所示,对参数曲线rp(u)、rp(w)均匀步长采样并求采样点数,即以待测点P0点为起点,沿P0周围的M、N、E、F四个方向按照步长h分别进行采样,步长h取0.1mm,其中M、N向为P0两侧沿着展向向外的方向,E、F向为P0两侧沿着弦向向外的方向;采样得到有序点集其中n为采样点数,且60mm≤nh≤100mm,得到四个有序点集的采样点数
其中,ΔMi-1MiMi+1为三角形Mi-1MiMi+1的面积,ΔNi-1NiNi+1为三角形Ni-1NiNi+1的面积,ΔEi-1EiEi+1为三角形Ei-1EiEi+1的面积,ΔFi-1FiFi+1为三角形Fi-1FiFi+1的面积。
步骤5:如附图3所示,以P0为起始点,分别沿M、N方向对参数曲线rp(u)上曲线段MnP0Nn的曲率求积分,分别沿E、F方向对参数曲线rp(w)上曲线段EnP0Fn的曲率求积分。取曲率积分值ΔA=0.1ω,ω为区域影响因子,其中当壁板处于机身、机翼区域时,ω取0.3~1,当壁板处于机翼前缘、尾翼前缘、整流罩、翼身融合处区域时,ω取1~2.5;本实施例中壁板处于机翼前缘,取ω=1.2,得到积分面积单元ΔA=0.12。
当时,说明壁板展向曲率变化大,则在曲线段MnP0Nn上提取四个测量点,在曲线段EnP0Fn上提取两个测量点,以提高法矢求解精度;由以下公式得到有序点集 中每块积分面积单元ΔA包含的采样点数,即测量点的点位I1,I2,I3,I4,I5,I6:
从而得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
当时,说明壁板弦向曲率变化大,则在曲线段MnP0Nn上提取两个测量点,在曲线段EnP0Fn上提取四个测量点,以提高法矢求解精度;由以下公式得到有序点集中每块积分面积单元ΔA包含的采样点数,即测量点的点位I1,I2,I3,I4,I5,I6:
从而得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
步骤6:采用激光跟踪仪或电涡流传感器测量P0点以及步骤5中得到的六个测量点对应实际壁板上的空间位置,得到P0点以及六个测量点的实际测量值:P′0,P′1,P′2,P′3,P′4,P′5,P′6;
当对P′1,P′0,P′2三次B样条插值得到一条三次B样条曲线r′p(u),对P′3,P′4,P′0,P′5,P′6三次B样条插值得到一条三次B样条曲线r′p(w);
分别计算P′0点处沿r′p(u)和r′p(w)两个方向的切矢P′(u)和P′(w),得到P′0点处法矢为:
Claims (1)
1.一种变间距自适应搜寻测量点的飞机壁板法矢确定方法,其特征在于:包括以下步骤:
步骤1:根据壁板理论数学模型提取待测点P0的空间坐标、P0点所在壁板表面的曲面方程r(u,w)、以及过P0点的展向曲线方程rp(u)和过P0点的弦向曲线方程rp(w);
步骤2:以P0点为起点,沿P0周围的M、N、E、F四个方向按照步长h分别进行采样,其中M、N向为P0两侧沿着展向向外的方向,E、F向为P0两侧沿着弦向向外的方向;采样得到有序点集 其中n为采样点数,且60mm≤nh≤100m;
其中,ΔMi-1MiMi+1为三角形Mi-1MiMi+1的面积,ΔNi-1NiNi+1为三角形Ni-1NiNi+1的面积,ΔEi-1EiEi+1为三角形Ei-1EiEi+1的面积,ΔFi-1FiFi+1为三角形Fi-1FiFi+1的面积;
步骤5:取曲率积分值ΔA=0.1ω,ω为区域影响因子,其中当壁板处于机身、机翼区域时,ω取0.3~1,当壁板处于机翼前缘、尾翼前缘、整流罩、翼身融合处区域时,ω取1~2.5;
计算得到测量点在步骤2中采样得到的有序点集中的点位I1,I2,I3,I4,I5,I6,则得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
计算得到测量点在步骤2中采样得到的有序点集中的点位I1,I2,I3,I4,I5,I6,则得到曲线段MnP0Nn和曲线段EnP0Fn上的测量点为:
步骤6:采用测量仪器测量P0点以及步骤5中得到的六个测量点对应实际壁板上的空间位置,得到P0点以及六个测量点的实际测量值:P0′,P1′,P2′,P3′,P4′,P5′,P6′;
分别计算P0′点处沿rp′(u)和rp′(w)两个方向的切矢P′(u)和P′(w),得到P0′点处法矢为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210083667 CN102607484B (zh) | 2012-03-27 | 2012-03-27 | 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210083667 CN102607484B (zh) | 2012-03-27 | 2012-03-27 | 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102607484A CN102607484A (zh) | 2012-07-25 |
CN102607484B true CN102607484B (zh) | 2013-12-25 |
Family
ID=46525091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210083667 Expired - Fee Related CN102607484B (zh) | 2012-03-27 | 2012-03-27 | 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102607484B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104034261B (zh) * | 2014-06-10 | 2016-12-07 | 深圳先进技术研究院 | 一种曲面法向测量装置和曲面法向测量方法 |
CN104019779B (zh) * | 2014-06-18 | 2016-09-21 | 清华大学 | 曲面法矢测量精度的优化方法 |
CN104006781B (zh) * | 2014-06-18 | 2017-01-04 | 清华大学 | 曲面法矢测量精度的计算方法 |
CN104697478B (zh) * | 2015-03-16 | 2016-12-28 | 厦门大学 | 基于恒定接触力的直纹面自适应测量方法 |
CN108803480B (zh) * | 2018-07-12 | 2021-04-30 | 深圳市汇川技术股份有限公司 | 轨迹间g2连续光顺转接方法、设备及计算机可读存储介质 |
CN109579680A (zh) * | 2018-11-22 | 2019-04-05 | 中国十七冶集团有限公司 | 一种临时门洞的位移监测方法 |
CN111664801A (zh) * | 2020-05-22 | 2020-09-15 | 成都飞机工业(集团)有限责任公司 | 一种基于高斯曲率的光学靶标点投影规划方法 |
CN111907727B (zh) * | 2020-06-17 | 2022-04-08 | 成都飞机工业(集团)有限责任公司 | 一种飞机机翼骨架外形测量点选取方法 |
CN111844043B (zh) * | 2020-07-31 | 2021-06-18 | 西北工业大学 | 一种机器人锪窝自适应控制***及控制方法 |
CN114266148A (zh) * | 2021-12-15 | 2022-04-01 | 成都飞机工业(集团)有限责任公司 | 一种用于航空结构件曲面群孔加工的patch设计方法 |
CN116680816B (zh) * | 2023-07-27 | 2023-11-10 | 成都飞机工业(集团)有限责任公司 | 一种飞机部件制孔法矢修正方法、装置、设备及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102248450A (zh) * | 2011-04-20 | 2011-11-23 | 上海交通大学 | 用于大曲率半径曲面法向矢量快速检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005044156A (ja) * | 2003-07-23 | 2005-02-17 | National Institute Of Advanced Industrial & Technology | 曲面の曲率演算方法及びプログラム |
-
2012
- 2012-03-27 CN CN 201210083667 patent/CN102607484B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102248450A (zh) * | 2011-04-20 | 2011-11-23 | 上海交通大学 | 用于大曲率半径曲面法向矢量快速检测方法 |
Non-Patent Citations (5)
Title |
---|
《面向曲面重构的切平面法矢方向调整算法》;刘大峰等;《机械科学与技术》;20080229;第27卷(第2期);第192-197页 * |
《飞机壁板自动钻铆法向量测量方法研究》;应高明等;《机床与液压》;20101231;第38卷(第23期);第1-4、8页 * |
JP特开2005-44156A 2005.02.17 |
刘大峰等.《面向曲面重构的切平面法矢方向调整算法》.《机械科学与技术》.2008,第27卷(第2期), |
应高明等.《飞机壁板自动钻铆法向量测量方法研究》.《机床与液压》.2010,第38卷(第23期), |
Also Published As
Publication number | Publication date |
---|---|
CN102607484A (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102607484B (zh) | 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 | |
CN102735204B (zh) | 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法 | |
Spentzos et al. | Investigation of three-dimensional dynamic stall using computational fluid dynamics | |
CN101672637B (zh) | 一种复杂曲面的数字化检测方法 | |
CN103486996B (zh) | 未知cad模型的航空发动机叶片型面测量方法 | |
Potsdam et al. | Unstructured mesh CFD aerodynamic analysis of the NREL Phase VI rotor | |
CN109766604B (zh) | 一种基于随机等几何分析的叶片高刚度设计方法 | |
CN103411574A (zh) | 航空发动机叶片型面三坐标测量方法 | |
CN103942837B (zh) | 基于序列线性规划的叶片点云模型截面曲线直接构造方法 | |
Bolzon et al. | Force measurements and wake surveys of a swept tubercled wing | |
CN105095540B (zh) | 一种井间相对位置关系的法面扫描解析方法 | |
CN103808286A (zh) | 一种基于全站仪的钢结构三维精度检测分析方法及其应用 | |
CN103900520A (zh) | 一种整体壁板板坯模型几何尺寸检测方法 | |
Castelli et al. | Numerical analysis of the influence of airfoil asymmetry on VAWT performance | |
CN103366056A (zh) | 基于有限元分析的机械零件圆柱度误差的测量方法 | |
CN109190232B (zh) | 一种飞机平尾区动能损失计算评估方法 | |
CN103177160B (zh) | 基于特征的飞机结构件自动装夹实现方法 | |
Venkatachari et al. | Assessment of transition modeling capability in OVERFLOW with emphasis on swept-wing configurations | |
CN105550383B (zh) | 一种非定常气动力测量试验***的设计方法 | |
Baker | Drag reduction of a blunt trailing-edge airfoil | |
CN105574221A (zh) | 一种改进cst翼型参数化方法 | |
CN102728658B (zh) | 一种用于确定叶片加工弯曲度误差的方法 | |
Mariens | Wing Shape multidisciplinary design optimization | |
CN112362290B (zh) | 一种机翼厚度公差对阻力系数影响的快速分析方法 | |
Wan et al. | Constructing process models of engine blade surfaces for their adaptive machining: an optimal approach |
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 |
Granted publication date: 20131225 Termination date: 20150327 |
|
EXPY | Termination of patent right or utility model |