CN103950552A - 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法 - Google Patents

基于六轴数控定位器的飞机壁板装配变形的数字化校正方法 Download PDF

Info

Publication number
CN103950552A
CN103950552A CN201410171705.XA CN201410171705A CN103950552A CN 103950552 A CN103950552 A CN 103950552A CN 201410171705 A CN201410171705 A CN 201410171705A CN 103950552 A CN103950552 A CN 103950552A
Authority
CN
China
Prior art keywords
check point
shaft
error
machining test
aircraft
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
Application number
CN201410171705.XA
Other languages
English (en)
Other versions
CN103950552B (zh
Inventor
毕运波
严伟苗
柯映林
屠晓伟
姜丽萍
沈立恒
朱宇
邢宏文
周庆慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Shanghai Aircraft Manufacturing Co Ltd
Original Assignee
Zhejiang University ZJU
Shanghai Aircraft Manufacturing Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, Shanghai Aircraft Manufacturing Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201410171705.XA priority Critical patent/CN103950552B/zh
Publication of CN103950552A publication Critical patent/CN103950552A/zh
Application granted granted Critical
Publication of CN103950552B publication Critical patent/CN103950552B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了一种基于六轴数控定位器的飞机壁板装配变形的数字化校正方法。本发明的数字化校正方法中,通过偏最小二乘回归反演建模方法,建立六轴数控定位器运动参数和检测点的位置误差数据之间的关系得到数字化校正模型,实现了大型飞机壁板装配变形的数字化校正,不仅有效降低了大型飞机壁板因装配变形引起的装配应力,同时保证了机身段装配中各个壁板的高效、高精度调姿和对接,最终提升了飞机大部件的装配质量。本发明的数字化校正方法通过六轴数控定位器的协调运动,成功解决了大型飞机壁板装配变形校正和准确定位问题,有效降低了大型飞机壁板因装配变形引起的装配应力,提升飞机大部件的装配质量。

Description

基于六轴数控定位器的飞机壁板装配变形的数字化校正方法
技术领域
本发明涉及飞机装配技术领域,尤其涉及一种基于六轴数控定位器的飞机壁板装配变形的数字化校正方法。
背景技术
飞机装配作为飞机制造环节中极其重要的一环,在很大程度上决定了飞机的最终质量、制造成本和交货周期,是整个飞机制造过程中的关键和核心技术。大型飞机一般由多个机身段对接装配而成,而各个机身段又由若干壁板组装拼接而成。壁板作为现代大型飞机的重要组件之一,既是构成飞机气动外形的重要组成部分,同时也是机身、机翼等的主要承力构件。壁板装配是将蒙皮、长桁、隔框、角片等薄壁类零件按照设计和技术要求进行定位、制孔并通过以铆接为主的手段进行连接而成,是飞机装配中极为重要的环节,但零件特性和装配方式往往造成其自身刚度、强度相对不足。
虽然航空制造企业正逐步采用整体壁板代替组装壁板,减少壁板所含零件数量,降低壁板整体重量,并在一定程度上提高了壁板的强度和刚度,提升了气动表面与外形的装配质量,但由于大型飞机壁板表面轮廓为复杂的空间自由曲面,面积较大,在装配过程中不利于承受集中载荷,同时加上定位误差、制孔、铆接、插螺栓、强迫装配以及自身重量、残余应力等实际装配因素的影响,往往导致壁板局部刚度过低,变形量超过容差范围,造成部件间交点不协调,使飞机最终的实际装配外形与理论外形存在较大偏差,影响整机的气动外形,并对后续工序产生不良影响。
虽然我国航空制造企业多通过绷带等工具进行强迫装配,或增加修配和精加工等工序的方式进行补救,达到总体精度要求,但这样势必会增加企业的生产成本,延长飞机的装配周期。因此,在大型飞机壁板装配过程中如何有效控制和减小壁板装配变形是目前我国航空工业亟需解决和攻克的重要技术难题之一。
发明内容
针对当前大型飞机壁板在装配中存在的变形问题,本发明提供了一种基于六轴数控定位器的飞机壁板装配变形的数字化校正方法。
一种基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,包括:
(1)在飞机壁板各个隔框上均匀布置若干个检测点;
(2)根据六轴数控定位器的数量、以及各个六轴数控定位器在X、Y、Z方向的移动量阈值和转动量阈值确定正交表,正交表的每一行表示六轴数控定位器的平移量和转动量;
(3)将所述正交表中的每一行数据作为载荷样本,将各个载荷样本加载至飞机壁板的理论有限元模型,获得各个载荷样本作用下飞机壁板变形有限元模型;
(4)针对任意一个载荷样本,利用对应的变形有限元模型计算当前载荷样本作用下各个检测点的位置误差和各个工艺球头球心的位置及转动误差;
(5)根据所有载荷样本作用下,各个检测点的位置误差和各个工艺球头球心的位置及转动误差,采用偏最小二乘回归反演建模方法,建立飞机壁板装配变形的数字化校正模型;
(6)获取各个检测点的实际位置误差,并将各个检测点的实际位置误差代入所述的数字化校正模型,计算得到各个六轴数控定位器的校形数据;
(7)根据所述的校形数据,对六轴数控定位器运动进行位置调整,完成大型飞机壁板的装配变形校正。
所述步骤(1)中的布置检测点时避免所有检测点位于同一直线上,尽量保证面分布。本发明中在飞机壁板的每个隔框上提取相同数量的检测点。本发明中检测点的个数为10~40。
根据六轴数控定位器的数量、以及各个六轴数控定位器在在X、Y、Z方向的移动量阈值和转动量阈值确定正交表,正交表的每一行表示各个六轴数控定位器的平移量和转动量;
所述步骤(2)中各个六轴数控定位器在在X、Y、Z方向的移动量阈值和转动量阈值取决于六轴数控定位器本身,各个六轴数控定位器的移动量阈值和转动量阈值在制作后就已经固定。根据各个六轴数控定位器的个数查表确定合适大小的正交表(正交表的列数和行数)。行数就是仿真试验样本数,直接取决于六轴数控定位器的个数。对于列数,对于六轴数控定位器,每个六轴数控定位器具有6个自由度,因此需要保证正交表的列数大于或等于6w,w为六轴数控定位器的个数。
本发明的数字化校正方法中,针对当前大型飞机壁板在装配中存在的变形问题,通过建立六轴数控定位器运动参数(位置及转动误差)和检测点的位置误差数据之间的反演计算得到数字化校正模型,实现了大型飞机壁板装配变形的数字化校正,不仅有效降低了大型飞机壁板因装配变形引起的装配应力,同时保证了机身段装配中各个壁板的高效、高精度调姿和对接,最终提升了飞机大部件的装配质量。
所述步骤(3)包括以下步骤:
(3-1)从理论有限元模型获取工艺球头的球心和各个检测点的理论坐标;
(3-2)将正交表中每一行数据作为一个载荷样本,将各个载荷样本施加至工艺球头的球心,并通过有限元模拟得到壁板装配变形的变形有限元模型。
飞机壁板的理论有限元模型为飞机壁板的有限元设计模型。一个载荷样本实际上包含一次仿真实验中各个六轴数控定位器在X、Y和Z方向上的平移量和转动量。将各个载荷样本施加至工艺球头的球心实际上是以载荷样本中各个六轴数控定位器的平移量和转动量为增量,改变各个轴数控定位器对应的工艺球头球心的坐标。有限元模拟通过有限元模拟仿真软件实现,首先输入飞机壁板的理论有限元模型,然后将载荷样本施输入有限元模拟仿真软件,即可得到该载荷样本作用下变形的变形有限元模型。
所述步骤(4)通过以下步骤计算各个检测点的位置误差和各个工艺球头球心的位置及转动误差:
(4-1)从变形有限元模型中获取各个检测点和工艺球头的球心的实际坐标;
(4-2)根据各个检测点的理论坐标和实际坐标计算得到姿态变换矩阵和位置平移向量;
(4-3)根据姿态变换矩阵和位置平移向量计算各个检测点的位置误差,以及各个工艺球头的球心的位置及转动误差,其中:
第i个检测点的位置误差为ei,根据公式:
ei=X'Ki-(RXKi+t)
计算得到,其中X'Ki为第i个检测点的实际坐标,XKi为第i个检测点的理论坐标值,R为姿态变换矩阵,t为位置平移向量,i=1,2,……,v,v为检测点的个数;
第j个工艺球头的球心的位置及转动误差fj为:
fj=[fj1,fj2,fj3,fj4,fj5,fj6],
其中,[fj1,fj2,fj3]=X'Bj-(RXBj+t),X'Bj为第j个工艺球头的球心的实际坐标,XBj为第j个工艺球头的球心的理论坐标值,[fj4,fj5,fj6]当前载荷样本对应的第j个六轴数控定位器的转动量,j=1,2,……,w,w为六轴数控定位器的个数。
所述步骤(5)中包括:
(5-1)将所有载荷样本作用下的各个检测点的位置误差组合得到检测点误差矩阵,将所有载荷样本作用下的各个工艺球头的球心的位置及转动误差组合得到球心误差矩阵;
(5-2)对检测点误差矩阵和球心误差矩阵进行标准化处理,并根据标准处理后的检测点误差矩阵和球心误差矩阵,采用偏最小二乘回归建模方法得到飞机壁板装配变形的数字化校正模型:
η=Aε+ηconst
其中,η为六轴数控定位器的校形数据,ε为各个检测点的位置误差,A为系数矩阵,ηconst为常数项。
以每个载荷样本作用下的各个检测点的位置误差作为检测点位置矩阵中的一行,从而得到检测点矩阵。以每个载荷样本作用下的各个球心的位置即转动误差作为球心误差矩阵中的一行,从而得到球心误差矩阵。
步骤(5-2)中标准化处理即单位化处理,将检测点矩阵和球心误差矩阵转化为单位矩阵,从而消除因单位不同引起的计算误差。
假设正交表的行数为n(即仿真试验的样本数为n),列数为6w,检测点的个数为v,因此构建的检测点误差矩阵为n×p(p=3v),球心误差矩阵大小为n×q(q=6w),检测点误差矩阵记为M=[m1,...,mp]n×p,球心误差矩阵记为N=[n1,...,nq]n×q
对检测点误差矩阵M和球心误差矩阵N进行标准化处理,得到标准化处理后的检测点误差矩阵和球心误差矩阵。其中,标准化处理后的检测点误差矩阵为E0=[E01,...,E0p]n×p,标准化处理后球心误差矩阵为F0=[F01,...,E0q]n×q
在偏最小二乘回归建模过程中,以标准化处理后的检测点误差矩阵E0为自变量,以标准化处理后的球心误差矩阵F0为因变量。记E0和F0的第一个主轴分别为w1和c1,而t1和u1分别为E0和F0的第一个主成分,且有t1=E0w1,u1=F0c1,并求解以下优化问题:
maxCov(t1,u1)=(E0w1)T·(F0c1)
s . t . : w 1 T w 1 = 1 c 1 T c 1 = 1
maxCov(t1,u1)表示对Cov(t1,u1)取最大值,s.t.为subject to的缩写,表示后面内容为约束条件。
引入拉格朗日乘子λ1和λ2,并记:
s = w 1 T E 0 T F 0 c 1 - λ 1 ( w 1 T w 1 - 1 ) - λ 2 ( c 1 T c 1 - 1 ) ,
对s分别求关于w1、c1、λ1和λ2的偏导数,可以推得:
θ 1 = 2 λ 1 = 2 λ 2 = w 1 T E 0 T F 0 c 1 , E 0 T F 0 F 0 T E 0 w 1 = θ 1 2 w 1 , F 0 T E 0 E 0 T F 0 c 1 = θ 1 2 c 1 ,
可见,w1是矩阵E0的特征向量,且对应的特征值为。θ1是目标函数值,它要求取最大值,所以,w1是对应于矩阵E0最大特征值的单位特征向量;同理,c1也应是对应于矩阵F0最大特征值的单位特征向量。
求得第一个轴w1和c1后,即可得到成分t1和u1,然后,分别求E0和F0对t1、u1的回归方程:
E 0 = t 1 p 1 T = t 1 E 0 t 1 T | | t 1 | | 2 + E 1
F 0 = t 1 r 1 T = t 1 F 0 t 1 T | | t 1 | | 2 + F 1 ,
其中,E1、F1分别是以上两式的残差矩阵。
用残差矩阵E1和F1取代E0和F0,然后求它们的第二个轴w2和c2以及第二个成分t2和u2,如此计算下去,如果检测点误差矩阵M的秩为λ,则会有下式成立:
E 0 = Σ i = 1 λ t i p i T ,
F 0 = Σ i = 1 λ t i r i T + F λ .
令ηi=F0i,εj=E0j,上式可进一步可表示为:
η=Aε+ηconst
该式即为壁板装配变形的数字化校正模型,其中η表示六轴数控定位器在X、Y、Z方向上的校形量(包括平移量和转动量),ε表示检测点的位置误差,A为系数矩阵,ηconst为常数项。其中,校形量η的大小为q×1,系数矩阵A的大小为q×p,检测点的位置误差ε的大小为p×1,常数项ηconst的大小为q×1。
所述步骤(6)中各个检测点的实际位置误差采用激光跟踪仪测量系统扫描六轴数控定位器支撑下飞机壁板得到。
实际位置误差实际上为装配时,飞机壁板中各个检测点的坐标(位置)相对于理论模型的位置的偏差,直接通过激光跟踪仪测量系统扫描测量得到,简单,且易于实现。
所述步骤(7)具体如下:
以校形数据为增量,调整六轴数控定位器的位置。
为实现飞机壁板装配变形的数字化校正,因此得到校形数据后,用校形数据调整六轴数控定位器的位置,从而消除变形引起的装配误差。
所述的装配变形数字化校正方法在完成大型飞机壁板的装配变形校正后对校正结果进行检测,具体如下:
获取校正完成后各个检测点的位置误差,将各个检测点的位置误差数据与设定的容差进行比较,若所有检测点的位置误差均在各自的容差内,则装配变形校正接结束;
否则,重新返回步骤(6)重新执行。
通过验证校正结果,判断校正结果是否正确,提高了该数字化校正方法的可实施性,且有利于提高校正精度。根据检测结果,若校正完成后,存在位置误差超出容差的检测点,则重新进行校正。重新校正时,直接从步骤(6)开始,不需要另外建立数字化校正模型。校正完成后,再次判断校正完成后各个检测点的位置误差是否在各自的容差内,并根据判断结果进一步循环进行,直至所有的检测点的位置误差是否均在各自的容差内。
作为优选,所述的容差为±0.5mm。
本发明中的容差为±0.5mm应理解为所有的检测点的容差的上限为±0.5mm。由于检测点的位置和类型不同,因此实际应用中各个检测点的容差也是相互独立的。对于不重要的检测点其容差可能较大,可能为±0.5mm,对于重要的检测点其容差可能较小,可能为±0.05mm。
与现有技术相比,本发明的优点在于:
(a)采用偏最小二乘回归反演建模方法建立了大型飞机壁板装配变形量与六轴数控定位器各运动参数之间的关系得到飞机壁板装配变形的数字化校正模型,并利用该数字化校正模型,通过六轴数控定位器的协调运动,成功解决了大型飞机壁板装配变形校正和准确定位问题;
(b)可有效降低大型飞机壁板因装配变形引起的装配应力,提升飞机大部件的装配质量。
附图说明
图1为大型飞机壁板整体结构示意图;
图2为大型飞机壁板的工艺接头结构示意图;
图3为六轴数控定位器结构示意图;
图4为大型飞机壁板定位调姿系统的示意图;
图5为基于六轴数控定位器的飞机壁板装配变形的数字化校正方法的流程图;
图中:大型飞机壁板1,蒙皮2,长桁3,隔框4,角片5,工艺接头6,工艺球头7,接头本体8,六轴数控定位器9,X向移动轴10,Y向移动轴11,Z向移动轴12,X向转动轴13,Y向转动轴14,Z向转动轴15,大型飞机壁板定位调姿系统16,控制系统计算机17,测量系统计算机18,检测点19,现场总线20,TCP/IP21,激光跟踪仪22。
具体实施方式
下面将结合附图和具体实施例对本发明进行详细说明。
如图1所示,大型飞机壁板1主要由蒙皮2、长桁3、隔框4、角片5组成。
如图2图所示,大型飞机壁板的工艺接头6主要由工艺球头7和接头本体8组成。
如图3所示,六轴数控定位器9包括X向移动轴10、Y向移动轴11、Z向移动轴12、X向转动轴13、Y向转动轴14、Z向转动轴15,各轴独立运动,并由控制系统计算机17实现操纵控制。其中,X向移动轴10、Y向移动轴11、Z向移动轴12构成的坐标系如图3所示。
工艺球头7球铰于相应的六轴数控定位器9的末端。
如图4图4所示,大型飞机壁板定位调姿系统16包括测量系统计算机18、控制系统计算机17以及激光跟踪仪22通过TCP/IP21实现网络通信与数据传递,控制系统计算机17则通过现场总线20实现对六轴数控定位器9的运动控制,测量系统计算机18指令激光跟踪仪22测得的大型飞机壁板1的检测点19的位置误差数据,并下发至控制系统计算机17,控制系统计算机17根据检测点19位置误差数据计算六轴数控定位器9校形量,并指令六轴数控定位器9移动至指定位置,整个系统形成了“测量-计算-校形”的闭环反馈回路,最终实现大型飞机壁板1的装配变形校正。该定位调姿系统能够直接获取检测点的位置和位置误差。
本实施例的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,如图5所示,包括:
(1)在飞机壁板各个隔框上均匀布置若干个检测点。
本实施例中为20个,各个检测点均匀分布(面分布),各个检测点位于飞机壁板的各个隔框。
(2)根据六轴数控定位器的数量、以及各个六轴数控定位器在在X、Y、Z方向的移动量阈值和转动量阈值确定正交表,正交表的每一行表示六轴数控定位器的平移量和转动量。
(3)将正交表中的每一行数据作为载荷样本,将各个载荷样本加载至飞机壁板的理论有限元模型,获得飞机壁板在各个载荷样本作用下的变形有限元模型。具体包括以下步骤:
(3-1)从理论有限元模型读取各个工艺球头的球心和各个检测点的理论坐标;
(3-2)将正交表中每一行数据作为一个载荷样本,将各个载荷样本施加至工艺球头的球心,并通过有限元模拟得到壁板装配变形的变形有限元模型。
(4)针对每一个载荷样本,利用该载荷样本作用下的变形有限元模型计算当前载荷样本作用下各个检测点的位置误差和各个工艺球头球心的位置及转动误差。具体如下:
(4-1)从变形有限元模型中获取各个检测点和工艺球头的球心的实际坐标;
(4-2)根据各个检测点的理论坐标和实际坐标计算得到姿态变换矩阵和位置平移向量;
(4-3)根据姿态变换矩阵和位置平移向量计算各个检测点的位置误差,以及各个工艺球头的球心的位置及转动误差,其中:
第i个检测点的位置误差为ei,根据公式:
ei=X'Ki-(RXKi+t)
计算得到,其中X'Ki为第i个检测点的实际坐标,XKi为第i个检测点的理论坐标值,R为姿态变换矩阵,t为位置平移向量,i=1,2,……,v,v为检测点的个数;
第j个工艺球头的球心的位置及转动误差fj为:
fj=[fj1,fj2,fj3,fj4,fj5,fj6],
其中,[fj1,fj2,fj3]=X'Bj-(RXBj+t),X'Bj为第j个工艺球头的球心的实际坐标,XBj为第j个工艺球头的球心的理论坐标值,[fj4,fj5,fj6]当前载荷样本对应的第j个六轴数控定位器的转动量,j=1,2,……,w,w为六轴数控定位器的个数。
(5)根据所有载荷样本作用下,各个检测点的位置误差和各个工艺球头球心下的位置及转动误差,采用偏最小二乘回归反演建模方法,建立飞机壁板装配变形的数字化校正模型。具体如下:
(5-1)将所有载荷样本作用下的各个检测点的位置误差组合得到检测点误差矩阵,将所有载荷样本作用下的各个工艺球头的球心的位置及转动误差组合得到球心误差矩阵;
(5-2)对检测点误差矩阵和球心误差矩阵进行标准化处理,并根据标准处理后的检测点误差矩阵和球心误差矩阵,采用偏最小二乘回归建模方法得到飞机壁板装配变形的数字化校正模型:
η=Αε+ηconst
其中,η为六轴数控定位器的校形数据,ε为各个检测点的位置误差,Α为系数矩阵,ηconst为常数项。
(6)获取各个检测点的实际位置误差,并将各个检测点的实际位置误差代入数字化校正模型,计算得到各个六轴数控定位器的校形数据。
各个检测点的实际位置误差采用激光跟踪仪测量系统扫描六轴数控定位器支撑下飞机壁板得到。
(7)根据所述的校形数据,以校形数据为增量,对六轴数控定位器运动进行位置调整,完成大型飞机壁板的装配变形校正。
(8)完成大型飞机壁板的装配变形校正后对校正结果进行检测,具体如下:
获取校正完成后所有检测点的位置误差数据,将各个检测点的位置误差数据与设定的容差进行比较,若所有检测点的位置误差数据均在各自的容差(本实施例中各个检测点的容差为±0.5mm)内,则装配变形校正接结束;
否则,重新返回步骤(6)重新执行。
以上所述仅为本发明的优选实施方式,本发明的保护范围并不仅限于上述实施方式,凡是属于本发明原理的技术方案均属于本发明的保护范围。对于本领域的技术人员而言,在不脱离本发明的原理的前提下进行的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,包括:
(1)在飞机壁板各个隔框上均匀布置若干个检测点;
(2)根据六轴数控定位器的数量,以及各个六轴数控定位器在X、Y、Z方向的移动量阈值和转动量阈值确定正交表,正交表的每一行表示六轴数控定位器的平移量和转动量;
(3)将所述正交表中的每一行数据作为载荷样本,将各个载荷样本加载至飞机壁板的理论有限元模型,获得各个载荷样本作用下飞机壁板变形有限元模型;
(4)针对任意一个载荷样本,利用对应的变形有限元模型计算当前载荷样本作用下各个检测点的位置误差和各个工艺球头球心的位置及转动误差;
(5)根据所有载荷样本作用下,各个检测点的位置误差和各个工艺球头球心的位置及转动误差,采用偏最小二乘回归反演建模方法,建立飞机壁板装配变形的数字化校正模型;
(6)获取各个检测点的实际位置误差,并将各个检测点的实际位置误差代入所述的数字化校正模型,计算得到各个六轴数控定位器的校形数据;
(7)根据所述的校形数据,对六轴数控定位器运动进行位置调整,完成大型飞机壁板的装配变形校正。
2.如权利要求1所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述步骤(3)包括以下步骤:
(3-1)从理论有限元模型获取工艺球头的球心和各个检测点的理论坐标;
(3-2)将正交表中每一行数据作为一个载荷样本,将各个载荷样本施加至工艺球头的球心,并通过有限元模拟得到壁板装配变形的变形有限元模型。
3.如权利要求2所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述步骤(4)通过以下步骤计算各个检测点的位置误差和各个工艺球头球心的位置及转动误差:
(4-1)从变形有限元模型中获取各个检测点和工艺球头的球心的实际坐标;
(4-2)根据各个检测点的理论坐标和实际坐标计算得到姿态变换矩阵和位置平移向量;
(4-3)根据姿态变换矩阵和位置平移向量计算各个检测点的位置误差,以及各个工艺球头的球心的位置及转动误差,其中:
第i个检测点的位置误差为ei,根据公式:
ei=X'Ki-(RXKi+t)
计算得到,其中X'Ki为第i个检测点的实际坐标,XKi为第i个检测点的理论坐标值,R为姿态变换矩阵,t为位置平移向量,i=1,2,……,v,v为检测点的个数;
第j个工艺球头的球心的位置及转动误差fj为:
fj=[fj1,fj2,fj3,fj4,fj5,fj6],
其中,[fj1,fj2,fj3]=X'Bj-(RXBj+t),X'Bj为第j个工艺球头的球心的实际坐标,XBj为第j个工艺球头的球心的理论坐标值,[fj4,fj5,fj6]当前载荷样本对应的第j个六轴数控定位器的转动量,j=1,2,……,w,w为六轴数控定位器的个数。
4.如权利要求3所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述步骤(5)中包括:
(5-1)将所有载荷样本作用下的各个检测点的位置误差组合得到检测点误差矩阵,将所有载荷样本作用下的各个工艺球头的球心的位置及转动误差组合得到球心误差矩阵;
(5-2)对检测点误差矩阵和球心误差矩阵进行标准化处理,并根据标准处理后的检测点误差矩阵和球心误差矩阵,采用偏最小二乘回归建模方法得到飞机壁板装配变形的数字化校正模型:
η=Aε+ηconst
其中,η为六轴数控定位器的校形数据,ε为各个检测点的位置误差,A为系数矩阵,ηconst为常数项。
5.如权利要求4所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述步骤(6)中各个检测点的实际位置误差采用激光跟踪仪测量系统扫描六轴数控定位器支撑下飞机壁板得到。
6.如权利要求5所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述步骤(7)具体如下:
以校形数据为增量,调整六轴数控定位器的位置。
7.如权利要求6所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述的装配变形数字化校正方法在完成大型飞机壁板的装配变形校正后对校正结果进行检测,具体如下:
获取校正完成后所有检测点的位置误差数据,将各个检测点的位置误差数据与设定的容差进行比较,若在所有检测点的位置误差数据均在各自的容差内,则装配变形校正接结束;
否则,重新返回步骤(6)重新执行。
8.如权利要求7所述的基于六轴数控定位器的飞机壁板装配变形的数字化校正方法,其特征在于,所述的容差为±0.5mm。
CN201410171705.XA 2014-04-25 2014-04-25 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法 Active CN103950552B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410171705.XA CN103950552B (zh) 2014-04-25 2014-04-25 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410171705.XA CN103950552B (zh) 2014-04-25 2014-04-25 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法

Publications (2)

Publication Number Publication Date
CN103950552A true CN103950552A (zh) 2014-07-30
CN103950552B CN103950552B (zh) 2015-08-26

Family

ID=51327950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410171705.XA Active CN103950552B (zh) 2014-04-25 2014-04-25 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法

Country Status (1)

Country Link
CN (1) CN103950552B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104816837A (zh) * 2015-05-07 2015-08-05 哈尔滨飞机工业集团有限责任公司 一种铝合金蒙皮数字化装配方法
CN105015800A (zh) * 2015-05-19 2015-11-04 北京星航机电装备有限公司 航天器舱段在地面上的自动化装配系统
CN106626434A (zh) * 2016-11-29 2017-05-10 中国航空工业集团公司沈阳飞机设计研究所 一种飞机复合材料隔框结构制造及设计方法
CN109110103A (zh) * 2018-09-12 2019-01-01 淮北创之社信息科技有限公司 一种飞机机舱壁板对接结构
CN112847174A (zh) * 2019-11-26 2021-05-28 财团法人工业技术研究院 可调式工件支撑系统及方法
CN112977875A (zh) * 2021-04-08 2021-06-18 泉州装备制造研究所 机翼翼盒装配间隙控制装置和多截面约束顺序优化方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5896637A (en) * 1996-09-25 1999-04-27 Mcdonnell Douglas Corporation Assembly tower
CN101362513A (zh) * 2008-09-19 2009-02-11 浙江大学 基于三个定位器的飞机部件位姿调整协同控制方法
CN102745338A (zh) * 2012-07-05 2012-10-24 浙江大学 一种大型飞机中机身数字化装配系统
CN103224036A (zh) * 2013-04-24 2013-07-31 大连理工大学 一种飞行器滚转体模型及制作方法
CN103274055A (zh) * 2013-06-14 2013-09-04 沈阳飞机工业(集团)有限公司 基于室内gps的飞机大部件无应力装配系统及其应用

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5896637A (en) * 1996-09-25 1999-04-27 Mcdonnell Douglas Corporation Assembly tower
CN101362513A (zh) * 2008-09-19 2009-02-11 浙江大学 基于三个定位器的飞机部件位姿调整协同控制方法
CN102745338A (zh) * 2012-07-05 2012-10-24 浙江大学 一种大型飞机中机身数字化装配系统
CN103224036A (zh) * 2013-04-24 2013-07-31 大连理工大学 一种飞行器滚转体模型及制作方法
CN103274055A (zh) * 2013-06-14 2013-09-04 沈阳飞机工业(集团)有限公司 基于室内gps的飞机大部件无应力装配系统及其应用

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104816837A (zh) * 2015-05-07 2015-08-05 哈尔滨飞机工业集团有限责任公司 一种铝合金蒙皮数字化装配方法
CN105015800A (zh) * 2015-05-19 2015-11-04 北京星航机电装备有限公司 航天器舱段在地面上的自动化装配系统
CN106626434A (zh) * 2016-11-29 2017-05-10 中国航空工业集团公司沈阳飞机设计研究所 一种飞机复合材料隔框结构制造及设计方法
CN106626434B (zh) * 2016-11-29 2019-03-08 中国航空工业集团公司沈阳飞机设计研究所 一种飞机复合材料隔框结构制造及设计方法
CN109110103A (zh) * 2018-09-12 2019-01-01 淮北创之社信息科技有限公司 一种飞机机舱壁板对接结构
CN112847174A (zh) * 2019-11-26 2021-05-28 财团法人工业技术研究院 可调式工件支撑系统及方法
CN112977875A (zh) * 2021-04-08 2021-06-18 泉州装备制造研究所 机翼翼盒装配间隙控制装置和多截面约束顺序优化方法
CN112977875B (zh) * 2021-04-08 2022-08-05 泉州装备制造研究所 一种机翼翼盒装配间隙多截面约束顺序优化方法

Also Published As

Publication number Publication date
CN103950552B (zh) 2015-08-26

Similar Documents

Publication Publication Date Title
CN110193829B (zh) 一种耦合运动学与刚度参数辨识的机器人精度控制方法
CN103950552B (zh) 基于六轴数控定位器的飞机壁板装配变形的数字化校正方法
CN103921954B (zh) 基于三轴数控定位器的飞机壁板装配变形的数字化校正方法
CN111367236B (zh) 一种面向加工过程的移动机器人系统标定方法及系统
CN100565406C (zh) 一种基于四个定位器的飞机部件位姿调整系统及方法
Tian et al. Auto-normalization algorithm for robotic precision drilling system in aircraft component assembly
US9435633B2 (en) Quasi-virtual locate/drill/shim process
Peng et al. Total differential methods based universal post processing algorithm considering geometric error for multi-axis NC machine tool
CN103979118B (zh) 一种机翼壁板数字化定位方法以及定位装置
CN106043736A (zh) 一种用于飞机部件调姿的串并联机械结构及调姿方法
Yongguo et al. Trajectory planning algorithm based on quaternion for 6-DOF aircraft wing automatic position and pose adjustment method
CN107414827B (zh) 基于线性反馈控制器的六自由度机械臂自适应检测方法
CN101362514A (zh) 基于三个定位器的飞机部件位姿调整系统及方法
CN110152211A (zh) 一种患者承载医用机械臂误差补偿系统及方法
CN106885676A (zh) 气动载荷产生的非解耦六自由度机构末端位姿误差补偿法
CN111695193A (zh) 一种全局相关三维气动力数学模型的建模方法及系统
Mei et al. Coordination error control for accurate positioning in movable robotic drilling
CN103171772A (zh) 一种飞机机身三段对合及交点孔加工的系统及布局方法
CN105798706A (zh) 一种用于飞机壁板卧式自动钻铆机的末端精度重力补偿方法
Qi et al. Accuracy Improvement Calibrations for the Double‐Position 4‐PPPS Aircraft Docking System
CN116147881B (zh) 一种风洞六分量天平静校系统中复位机构的复位方法
CN112157654A (zh) 一种机器人加工系统定位误差的优化方法
CN112307562B (zh) 综合热变形与重力变形的大型飞机上复杂部件的装配方法
CN107665616B (zh) 一种九自由度运动模拟器相对运动等效方法及系统
CN106840587B (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