CN113065259A - 一种壁板组件与骨架位姿协调方法 - Google Patents
一种壁板组件与骨架位姿协调方法 Download PDFInfo
- Publication number
- CN113065259A CN113065259A CN202110406038.9A CN202110406038A CN113065259A CN 113065259 A CN113065259 A CN 113065259A CN 202110406038 A CN202110406038 A CN 202110406038A CN 113065259 A CN113065259 A CN 113065259A
- Authority
- CN
- China
- Prior art keywords
- assembly
- wing box
- coordination
- pose
- skeleton
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Toys (AREA)
Abstract
Description
技术领域
本发明涉及飞机装配领域,尤其涉及一种改进的飞机翼盒壁板组件与骨架位姿协调方法。
背景技术
在飞机机翼翼盒数字化装配系统中,采用壁板调姿定位系统实现壁板与骨架的位姿协调与对合,首先基于最佳匹配原理计算位姿协调参数,然后通过多数控定位器协同运动,实现壁板位姿调整及其与翼盒骨架的对合操作。但大型飞机翼盒零组件尺寸大、几何外形复杂,极易受温度等影响产生变形,导致翼盒装配间隙的出现,且壁板与骨架之间的装配间隙呈多区域分布的特点,对合时如果壁板与翼盒骨架的相对位姿状态不佳,可能导致翼盒局部区域装配间隙较大。现有技术中,有关飞机翼盒壁板与骨架位姿协调算法中,未能充分考虑各种变形因素对翼盒装配的综合影响,无法精确计算出有效的装配位姿协调参数,导致最终装配效果不佳。
发明内容
本发明要解决的技术问题,在于提供一种壁板组件与骨架位姿协调方法,实现最优化位姿协调参数的精准、快速查找,提高装配效率和装配质量。
本发明提供了一种壁板组件与骨架位姿协调方法,包括:
步骤S1、计算飞机翼盒各个翼肋处装配间隙Sgap(j);
步骤S3、基于壁板与骨架关键特征点集三维协调误差的加权范数构建最小二乘目标函数Fobj,所述最小二乘目标函数Fobj为:
其中,ωi为翼盒上第i个关键特征点的权重值,S表示位姿协调六元组参数,S=(x,y,z,α,β,γ),x,y,z表示刚体位置,α,β,γ表示刚体姿态,n表示翼盒关键特征点总个数;
步骤S5、求解使翼盒壁板与骨架关键特征点集协调误差加权范数以及多区域装配间隙加权范数综合评价指标最小化的位姿协调参数Sopt=[xopt,yopt,zopt,αopt,βopt,γopt],作为最优位姿协调参数输出结果。
进一步的,所述装配间隙Sgap(j)的计算公式如下:
Sgap(j)=∫[Cpanel(x)-Crib(x)]dx;
其中,Cpanel为飞机翼盒翼肋处的壁板内形轮廓,Crib为对应翼肋处的外形轮廓。
进一步的,所述壁板内形轮廓Cpanel和翼肋外形轮廓Crib通过3D扫描仪器设备测量得到。
其中,所述ΔKC(T,g)panel是基于机翼装配系统热弹变形物理仿真模型求解得到的壁板与翼盒骨架位姿协调时由温度T和重力g误差源导致的壁板关键特征点偏差,所述ΔKC(T,g)skeleton是基于机翼装配系统热弹变形物理仿真模型求解得到的壁板与翼盒骨架位姿协调时由温度T和重力g误差源导致的骨架关键特征点偏差。
进一步的,所述步骤S5中求解操作具体为,以装配间隙容差为约束,融合非迭代SVD算法和非线性L-M算法进行求解。
本发明的优点在于:
本发明的综合飞机结构变形和装配间隙的翼盒壁板与骨架位姿协调方法,可以有效提升壁板与骨架之间的位姿协调效率,并使所有翼盒装配间隙的总和最小,从而为后续装配间隙抑制和装配偏差控制提供良好的前提条件和数据依据,优化装配工艺,提高飞机装配质量和效率。
附图说明
下面参照附图结合实施例对本发明作进一步的说明。
图1为本发明一种壁板组件与骨架位姿协调方法的执行流程图。
图2为飞机翼盒壁板与骨架位姿协调对装配间隙影响示意图。
具体实施方式
如图1和图2所示,本发明的一种壁板组件与骨架位姿协调方法,包括:
步骤S1、计算飞机翼盒各个翼肋处装配间隙Sgap(j);该装配间隙为机翼壁板内形和翼肋外形之间的层间间隙;较佳的,所述装配间隙Sgap(j)的计算公式如下:
Sgap(j)=∫[Cpanel(x)-Crib(x)]dx;
其中,Cpanel为飞机翼盒翼肋处的壁板内形轮廓,Crib为对应翼肋处的外形轮廓。较佳的,所述壁板内形轮廓Cpanel和翼肋外形轮廓Crib通过3D扫描仪器设备测量得到。在实际操作过程中可通过将测量误差以高斯噪声形式引入其内外形轮廓模型,然后通过该积分公式求解两条轮廓曲线之间的面积,即可得到某翼肋处的壁板内形和翼肋外形之间的装配间隙。如图2所示,骨架外形(即翼肋外形)与壁板内形之间均为装配间隙。
其中,所述ΔKC(T,g)panel是基于机翼装配系统热弹变形物理仿真模型求解得到的壁板与翼盒骨架位姿协调时由温度T和重力g误差源导致的壁板关键特征点偏差,所述ΔKC(T,g)skeleton是基于机翼装配系统热弹变形物理仿真模型求解得到的壁板与翼盒骨架位姿协调时由温度T和重力g误差源导致的骨架关键特征点偏差。通过上述的ΔKC(T,g)panel和ΔKC(T,g)skeleton补偿壁板与翼盒骨架关键特征点集(该关键特征点集由翼盒装配工艺给定)实际测量数据中的热弹变形误差,得到消除翼盒热弹变形的壁板与骨架关键特征点集测量数据和从而提高计算精度。由于温度和重力对复杂外形的飞机结构的影响比较复杂,所以这里采用有限元分析方法(即采用基于有限元分析软件自行建立的机翼装配系统热弹变形物理仿真模型)计算出温度和重力造成的非线性的变形偏差,然后在测量数据中去除这部分偏差,从而将温度和重力耦合条件下壁板与翼盒骨架关键特征点集的非线性映射关系线性化。
步骤S3、基于壁板与骨架关键特征点集三维协调误差的加权范数构建最小二乘目标函数Fobj,所述最小二乘目标函数Fobj为:
其中,ωi为翼盒上第i个关键特征点的权重值,S为位姿协调六元组参数,S=(x,y,z,α,β,γ),x,y,z表示刚体位置,α,β,γ表示刚体姿态,n表示翼盒关键特征点总个数;六元组参数是刚体位姿的一种表达方式,α,β,γ表示刚体姿态,即ZYX欧拉角。翼盒上的关键特征点是以“点对”形式存在的,每一翼盒上的关键特征点包括对应的壁板上的关键特征点以及骨架上关键特征点。
考虑工程实际中不同关键特征点的重要性不尽相同,为各关键特征点对的协调误差施加不同权重ωi,当只考虑各关键特征点的影响时,可将基于线性Singular ValueDecomposition(SVD)算法求解的位姿协调参数Ssvd=(xsvd,ysvd,zsvd,αsvd,βsvd,γsvd)作为非线性Levenberg-Marquard(L-M)算法的初值,可计算出局部最优位姿参数Slm=[xlm,ylm,zlm,αlm,βlm,γlm]使目标函数Fobj最小化,此处的Ssvd表示基于线性SingularValueDecomposition(SVD)算法求解的位姿协调参数,Slm表示基于非线性Levenberg-Marquard(L-M)算法计算的位姿协调参数;
步骤S5、求解使翼盒壁板与骨架关键特征点集协调误差加权范数以及多区域装配间隙加权范数综合评价指标最小化的位姿协调参数Sopt=[xopt,yopt,zopt,αopt,βopt,γopt],作为最优位姿协调参数输出结果。较佳的,所述步骤S5中求解操作具体为,以装配间隙容差为约束,融合非迭代SVD算法和非线性L-M算法进行求解。
飞机装配时通常基于待调姿组件和调姿基准对应关键特征点集的最佳匹配求解位姿协调六元组参数S=(x,y,z,α,β,γ),进而协调组件之间的位姿关系。
本发明的飞机翼盒壁板与骨架位姿协调方法中,基于飞机大部件装配间隙模型,综合考虑翼盒多区域装配间隙Sgap(j)的加权范数,考虑工程实际中不同区域的装配间隙对翼盒装配质量的影响不同,为各区域装配间隙施加不同权重υj,另外,为了平衡关键特征点集协调误差和翼盒装配间隙对位姿协调参数最优化求解收敛过程的影响,为两者施加不同权重,使其更加贴合实际装配对象和装配工艺,能够提供更加精准的壁板和骨架位姿最优协调方案,为壁板与骨架对合装配提供数据依据,提高装配效率和装配质量。
虽然以上描述了本发明的具体实施方式,但是熟悉本技术领域的技术人员应当理解,我们所描述的具体的实施例只是说明性的,而不是用于对本发明的范围的限定,熟悉本领域的技术人员在依照本发明的精神所作的等效的修饰以及变化,都应当涵盖在本发明的权利要求所保护的范围内。
Claims (5)
1.一种壁板组件与骨架位姿协调方法,其特征在于:包括:
步骤S1、计算飞机翼盒各个翼肋处装配间隙Sgap(j);
步骤S3、基于壁板与骨架关键特征点集三维协调误差的加权范数构建最小二乘目标函数Fobj,所述最小二乘目标函数Fobj为:
其中,ωi为翼盒上第i个关键特征点的权重值,S表示位姿协调六元组参数,S=(x,y,z,α,β,γ),x,y,z表示刚体位置,α,β,γ表示刚体姿态,n表示翼盒关键特征点总个数;
步骤S5、求解使翼盒壁板与骨架关键特征点集协调误差加权范数以及多区域装配间隙加权范数综合评价指标最小化的位姿协调参数Sopt=[xopt,yopt,zopt,αopt,βopt,γopt],作为最优位姿协调参数输出结果。
2.如权利要求1所述的一种壁板组件与骨架位姿协调方法,其特征在于:所述装配间隙Sgap(j)的计算公式如下:
Sgap(j)=∫[Cpanel(x)-Crib(x)]dx;
其中,Cpanel为飞机翼盒翼肋处的壁板内形轮廓,Crib为对应翼肋处的外形轮廓。
3.如权利要求2所述的一种壁板组件与骨架位姿协调方法,其特征在于:所述壁板内形轮廓Cpanel和翼肋外形轮廓Crib通过3D扫描仪器设备测量得到。
5.如权利要求1所述的一种壁板组件与骨架位姿协调方法,其特征在于:所述步骤S5中求解操作具体为,以装配间隙容差为约束,融合非迭代SVD算法和非线性L-M算法进行求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406038.9A CN113065259B (zh) | 2021-04-15 | 2021-04-15 | 一种壁板组件与骨架位姿协调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110406038.9A CN113065259B (zh) | 2021-04-15 | 2021-04-15 | 一种壁板组件与骨架位姿协调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113065259A true CN113065259A (zh) | 2021-07-02 |
CN113065259B CN113065259B (zh) | 2022-05-31 |
Family
ID=76566715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110406038.9A Active CN113065259B (zh) | 2021-04-15 | 2021-04-15 | 一种壁板组件与骨架位姿协调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113065259B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114019915A (zh) * | 2021-11-04 | 2022-02-08 | 山东大学 | 多移动机器人协同运动控制系统及其控制方法 |
CN114954997A (zh) * | 2022-07-14 | 2022-08-30 | 成都飞机工业(集团)有限责任公司 | 一种舱门装配阶差的控制方法、装置、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103979118A (zh) * | 2014-04-01 | 2014-08-13 | 浙江大学 | 一种机翼壁板数字化定位方法以及定位装置 |
CN106514343A (zh) * | 2016-12-14 | 2017-03-22 | 中航飞机股份有限公司西安飞机分公司 | 一种飞机壁板自动制孔铆接定位工装及定位方法 |
US20190005162A1 (en) * | 2017-06-29 | 2019-01-03 | The Boeing Company | Methods and apparatus to generate a synthetic point cloud of a spacecraft |
CN112193432A (zh) * | 2020-08-21 | 2021-01-08 | 浙江大学 | 飞机壁板的安装方法及其与骨架间的装配间隙的计算方法 |
CN112307561A (zh) * | 2020-10-30 | 2021-02-02 | 泉州装备制造研究所 | 大型装备上多体复杂结构装配用的位姿协调方法及系统 |
-
2021
- 2021-04-15 CN CN202110406038.9A patent/CN113065259B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103979118A (zh) * | 2014-04-01 | 2014-08-13 | 浙江大学 | 一种机翼壁板数字化定位方法以及定位装置 |
CN106514343A (zh) * | 2016-12-14 | 2017-03-22 | 中航飞机股份有限公司西安飞机分公司 | 一种飞机壁板自动制孔铆接定位工装及定位方法 |
US20190005162A1 (en) * | 2017-06-29 | 2019-01-03 | The Boeing Company | Methods and apparatus to generate a synthetic point cloud of a spacecraft |
CN112193432A (zh) * | 2020-08-21 | 2021-01-08 | 浙江大学 | 飞机壁板的安装方法及其与骨架间的装配间隙的计算方法 |
CN112307561A (zh) * | 2020-10-30 | 2021-02-02 | 泉州装备制造研究所 | 大型装备上多体复杂结构装配用的位姿协调方法及系统 |
Non-Patent Citations (1)
Title |
---|
陈磊: "数字化测量辅助的飞机翼身对接装配协调技术研究", 《中国博士学位论文工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114019915A (zh) * | 2021-11-04 | 2022-02-08 | 山东大学 | 多移动机器人协同运动控制系统及其控制方法 |
CN114954997A (zh) * | 2022-07-14 | 2022-08-30 | 成都飞机工业(集团)有限责任公司 | 一种舱门装配阶差的控制方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113065259B (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111964575B (zh) | 一种面向移动机器人铣削加工的数字孪生建模方法 | |
CN106709124B (zh) | 用于垫补弯曲体的系统和方法 | |
CN113065259B (zh) | 一种壁板组件与骨架位姿协调方法 | |
Peng et al. | Total differential methods based universal post processing algorithm considering geometric error for multi-axis NC machine tool | |
Maropoulos et al. | A new paradigm in large-scale assembly—research priorities in measurement assisted assembly | |
US10450053B2 (en) | Methods using predictive shimming to optimize part-to-part alignment | |
CN103164584B (zh) | 基于关键特性的协调准确度计算方法 | |
Wang et al. | An assembly gap control method based on posture alignment of wing panels in aircraft assembly | |
EP3382581B1 (en) | Method for designing a die surface | |
CN112697041B (zh) | 一种基于蒙特卡洛法的装配位姿测量精度预评估方法 | |
CN112380616B (zh) | 高复杂易形变航天舱段高精度数字孪生对接装配方法 | |
CN107545080B (zh) | 一种基于bim模型的钢桥虚拟预拼装方法 | |
Schleich et al. | A quantitative comparison of tolerance analysis approaches for rigid mechanical assemblies | |
CN112613118A (zh) | 火箭发动机内部不可测装配质量数字孪生建模与追溯方法 | |
CN116611166B (zh) | 一种考虑初始误差的飞机部件多层级装配偏差分析方法 | |
Wang et al. | 3D machining allowance analysis method for the large thin-walled aerospace component | |
CN103921954B (zh) | 基于三轴数控定位器的飞机壁板装配变形的数字化校正方法 | |
CN103950552B (zh) | 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法 | |
CN114969976B (zh) | 基于数字化实测数据的一体化结构虚拟装配方法 | |
Maset et al. | Procrustes analysis for the virtual trial assembly of large-size elements | |
Guo et al. | Analysis on quantifiable and controllable assembly technology for aeronautical thin-walled structures | |
CN114239170A (zh) | 一种飞机装配的误差补偿方法 | |
Guo et al. | Assembly error propagation modeling and coordination error chain construction for aircraft | |
Muelaner et al. | Integrated dimensional variation management in the digital factory | |
CN112307561B (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 |