CN114370826B - 基于热胀系数反求的测量基准偏差修正方法 - Google Patents
基于热胀系数反求的测量基准偏差修正方法 Download PDFInfo
- Publication number
- CN114370826B CN114370826B CN202111569796.9A CN202111569796A CN114370826B CN 114370826 B CN114370826 B CN 114370826B CN 202111569796 A CN202111569796 A CN 202111569796A CN 114370826 B CN114370826 B CN 114370826B
- Authority
- CN
- China
- Prior art keywords
- reference point
- temperature
- thermal expansion
- measurement
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
- G01B21/04—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
- G01B21/042—Calibration or calibration artifacts
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明基于热胀系数反求的测量基准偏差修正方法,属于几何量测量领域,涉及一种基于热胀系数反求的测量基准偏差修正方法。该方法首先在工装局部主动温控条件下同步测量温度和基准点坐标;其次重建温度场函数,并基于温度和基准点坐标偏移间的映射关系反求最优热胀系数。然后,采用比例缩放法修正基准点坐标;最后以最小二乘法建立测量基准与数模基准间的转换关系,求解基准修正前后的转换误差。该方法能够实现对工装热胀系数的现场标定,完成非均匀温度环境中基准偏差修正。可提高测量基准的可靠性和一致性,可应用于现场非均匀温度条件大型构件制造中测量基准与数模基准匹配、激光跟踪仪转站误差修正等工程领域。
Description
技术领域
本发明属于几何量测量领域,涉及一种基于热胀系数反求的测量基准偏差修正方法。
背景技术
随着我国远程宽体客机、载人空间站等航空航天重大战略的推进和实施,大型构件产品的制造要求也越来越高。在大型构件的加工装配过程中,测量基准用于测量数据与构件数模间的匹配,进而指导制造工艺,因此其精度和稳定性是保障制造质量和工艺一致性的前提。此类基准由构件工装上的若干基准点组成,且与数模基准间存在已知的转换关系。然而,受现场非受控、非均匀的时变温度影响,工装在长周期制造中容易发生膨缩变形,导致测量基准发生动态偏移,严重影响测量基准的精度和稳定性。因此,为了保障大型构件产品的加工装配质量,必须考虑现场温度对测量基准的影响,修正基准偏差。
目前,基于线性热胀系数的热变形计算方法已应用于测量偏差的修正中。例如,浙江大学的杨宝旒、俞慈君等,在《航空学报》2015年第36卷第10期的“激光跟踪仪转站热变形误差建模与补偿方法”一文中,公开了一种基于飞机壁板工装单位温度热变形系数矩阵补偿理论坐标的方法。该方法基于空间中温度均匀分布假设,通过实验法推导出了单位热变形系数矩阵,对测量基准点进行了线性修正。然而,装配车间中的环境温度是非均匀的,在垂直地面方向存在较大的温度梯度,因此基于均匀温度假设的修正方法难以适用于现场工况。
针对非均匀温度场中构件热变形的计算,南京航空航天大学的曾捷等人,在专利《薄壁圆管结构轴向热变形计算方法》、专利号CN201610476246.5中,公开了一种通过有限元仿真模拟薄构件热变形的方法。该方法将查表得到的材料热胀系数代入有限元模型进行计算,但由于工装材料成分差异、测量系统误差传递等因素,现场工况下工装材料的真实热胀系数与理论值间存在较大差异,基于查表法的热变形计算方法精度低;而目前材料热胀系数的检定需在恒温计量室进行,步骤繁琐,成本高昂。
发明内容
本发明的目的是克服现有技术的不足,发明一种基于热胀系数反求的测量基准偏差修正方法。该方法首先基于主动温控以构建工装的局部温度场,通过对离散点温度采集和基准点坐标测量获取反求数据,其次基于最小二乘法重建温度场函数并反求热胀系数,然后基于比例缩放法修正基准点坐标,最后与数模坐标系进行基准转换。本发明能够实现对工装热胀系数的现场标定,并修正非均匀温度场中工装热变形引起的基准偏差,可提高测量基准的可靠性,服务于长周期大温差现场环境下大型构件的精密制造,具有广泛的应用前景。
本发明采用的技术方案是一种基于热胀系数反求的测量基准偏差修正方法,其特征是,该方法首先对工装局部结构进行主动温控,并同步测量温度和基准点坐标;其次重建温度场函数,并基于温度和基准点坐标偏移间的映射关系反求最优热胀系数;然后在实际工况下,基于温度测量数据和反求的最优热胀系数,采用比例缩放法修正基准点坐标;最后以最小二乘法建立测量坐标系与数模坐标系间的转换关系,求解基准修正前后的转换误差;该方法的具体步骤如下:
第一步,主动温控下的温度-坐标同步测量
在工装上放置两个激光跟踪仪靶球,分别作为基准点A和基准点B,在两点之间布设m个温度传感器,然后使用加热装置对工装进行加热,通过温度传感器获取tj,j=1,2,…n时刻的温度数据集TEMPj:
TEMPj=[T(x1,tj) T(x2,tj) … T(xi,tj) … T(xm,tj)]T (1)
式中,xi,i=1,2,…m为温度传感器的位置,tj,j=1,2,…n为温度数据的采样时刻,T(xi,tj)为第i个温度传感器在tj时刻测的温度值;
温度传感器测量的同时,使用激光跟踪仪连续测量基准点A和基准点B的空间坐标,按照公式(2)计算加热过程中基准点A和基准点B间距离变化数据集DELT:
第二步,温度场重建与热胀系数反求
在tj时刻,求解式(3)所示的最小二乘方程,对温度数据进行拟合,重建此时刻的温度场函数f(x,tj):
tj时刻,计算在温度场f(x,tj)作用下基准点B相对基准点A的偏移量为:
式中,α为工装材料的热胀系数,为待求参数;
计算n个温度场函数与t1时刻温度场函数差的积分:
由公式(4)可知DELT与INT_TEMP间是线性关系,比例系数即为热胀系数,因此建立线性回归模型以求解最优热胀系数αbest:
DELT=αbest×INT_TEMP+b (6)
第三步,基于比例缩放法的基准点偏差修正
在飞机构件装配的实际工况下,假设温度传感器在tn+1时刻测量的工装温度值为:
TEMPn+1=[T(x1,tn+1) T(x2,tn+1) … T(xm,tn+1)] (7)
式中,αbest为根据公式(6)求得的最优热胀系数,f(x,tn+1)为tn+1时刻的温度场函数;Tnom为假设不发生热变形的参考温度值,一般为20℃;Ln+1为基准点A和基准点B间的距离:
计算基准点B相对基准点A的比例缩放系数:
则修正后的基准点B坐标为:
第四步,测量基准与数模基准转换偏差计算
使用激光跟踪仪测量工装上N个基准点的坐标,获得坐标集S:
采用所述步骤1~步骤3的方法,获得修正后的坐标集Sc:
上述N个基准点在工装数模坐标系中的坐标集为:
建立如公式(15)所示的最小二乘模型,求解测量坐标系和数模坐标系间的基准转换参数(旋转矩阵R和平移向量T),并计算修正前的基准转换偏差E:
E=||D-(R×S+T)|| (16)
式中,S(k)和D(k)分别表示测量坐标集S和D中第k列向量;
同理,以公式(13)中修正后的坐标集Sc替换公式(15)中的坐标集S,计算测量基准修正后的基准转换偏差Ec:
Ec=||D-(Rc×Sc+Tc)|| (18)
由公式(18)计算出测量基准修正后的基准转换偏差Ec。
本发明的有益效果是该方法虑及了工装热胀系数实际值与理论值的差异,提出了基于局部主动温控和现场测量实验的热胀系数反求方法,可实现对工装热胀系数的现场精确标定,有效避免了热胀系数值不准确、误差大的问题;并通过融合离散温度数据与坐标测量数据,实现了测量基准点偏差的有效修正,适用于非均匀温度环境;此外,除温度传感器和加热装置外,本发明没有引入其他额外设备,且对温度测量精度和温度场重建方法要求低,成本低廉,便于现场快速实施;本发明可应用于现场非均匀温度条件大型构件装配中测量基准与数模基准匹配、激光跟踪仪转站误差修正等,具有广阔的应用前景。
附图说明
图1是基于热胀系数反求的测量基准偏差修正方法流程图;
图2是工装局部主动温控下的温度-坐标同步测量示意图;其中,1-工装;2-加热带;2.1-加热带控制器;3-温度传感器;3.1-温度数据采集卡;4-激光跟踪仪;A-基准点A;B-基准点B;
图3是温度场函数重建结果示意图;其中,1-工装局部;A-基准点A;B-基准点B;3-温度传感器;T(x1,t1)-位置坐标为x1的温度传感在t1时刻测得的温度值;f(x,t1)-t1时刻重建的温度场函数;T(x1,t210)-位置坐标为x1的温度传感器在t210时刻测得的温度值;f(x,t210)-t210时刻重建的温度场函数;
图4是基于线性回归的热胀系数反求示意图;
图5是实际工况下激光跟踪仪测量工装基准点示意图;其中,1-工装;3-温度传感器;4-激光跟踪仪;P1~P6-基准点;
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施方式。
本实施例选用的工装主要尺寸为:长5m,宽0.5m,高1.5~2.5m,材料为2024航空铝;选用温度传感器为K型热电偶,测量精度为±1.5℃;加热装置为矩形硅橡胶加热带,加热范围20~100℃,主要尺寸为800×50×2mm;激光跟踪仪型号为Leica AT960MR,其测量不确定度U=±(15μm+6μm/m);用作基准点的靶球尺寸为0.5英寸。
附图1是方法的流程图,该方法首先对工装局部结构进行主动温控,并同步测量温度和基准点坐标;其次重建温度场函数,并基于温度和基准点坐标偏移间的映射关系反求最优热胀系数;然后在实际工况下,基于温度测量数据和反求的最优热胀系数,采用比例缩放法修正基准点坐标;最后以最小二乘法建立测量坐标系与数模坐标系间的转换关系,求解基准修正前后的转换误差;具体步骤如下:
第一步,主动温控下的温度-坐标同步测量
如附图2所示,在工装1两端放置两个激光跟踪仪靶球,分别作为基准点A和基准点B,在两点之间等间距布设6个温度传感器3,然后在工装上粘贴加热带2对工装进行加热,通过温度传感器3获取300组基准点A和基准点B间的温度变化数据集:
TEMP1=[23.0 22.5 22.5 22.4 22.5 22.6]T
TEMP210=[28.5 30.9 31.8 29.2 33.7 31.5]T
TEMP300=[29.8 31.8 32.4 30.5 34.7 34.0]T
温度传感器测量的同时,使用激光跟踪仪连续测量基准点A和基准点B的空间坐标,按照公式(2)计算加热过程中基准点A和基准点B间距离变化数据集:
DELT=[0 0.0039 … 0.1996 … 0.2283]T (300×1)
第二步,温度场重建与热胀系数反求
以t210时刻为例,求解式(3)所示的最小二乘方程,对温度数据进行一次拟合,重建此时刻的温度场函数f(x,t210),如附图3所示:
f(x,t210)=0.0034x+29.3935
按照公式(5)计算500组温度场函数与t1时刻温度场函数差的积分,其中t1、t210、t300时刻的值为:
INT_TEMP=[0 … 7.312×103 … 8.421×103]T
由公式(4)可知DELT与INT_TEMP间是线性关系,比例系数即为热胀系数,如附图4所示,建立线性回归模型以求解最优热胀系数:
DELT=(2.733×10-5)×INT_TEMP-0.0026
因此最优热胀系数αbest=2.733×10-5/℃。
第三步,基于比例缩放法的基准点偏差修正
在飞机构件装配的实际工况下,假设温度传感器在t501时刻测量的工装温度值为:
TEMP406=[32.8 36.4 37.2 34.1 41.0 37.4]
同时,激光跟踪仪测得基准点A坐标为:
测得基准点B的坐标为:
计算基准点B相对基准点A的距离L301=875.968mm,根据公式(8)算得的偏差为:
式中,αbest为最优热胀系数,f(x,t301)为重建的温度场函数;20为假设不发生热变形的参考温度值。
根据公式(10)计算基准点B相对基准点A的比例缩放系数:
根据公式(11)计算修正后的基准点B坐标为:
第四步,测量基准与数模基准转换偏差计算
如附图4所示,使用激光跟踪仪测量工装上个6基准点的坐标,获得坐标集S:
采用所述步骤1~步骤3的方法,获得修正后的坐标集Sc:
上述6个基准点在工装数模坐标系中的坐标集为:
根据修正前测量坐标集S和数模中坐标集D,建立公式(15)所示的最小二乘模型,求解测量坐标系和数模坐标系间的基准转换参数:
根据公式(16)计算修正前的基准转换偏差:
E=0.4608mm
同理,以公式(13)中修正后的坐标集Sc替换公式(15)中的坐标集S,计算测量坐标系和数模坐标系间的基准转换参数:
根据公式(18)计算测量基准修正后的基准转换偏差:
Ec=0.0048mm
本发明虑及了工装热胀系数实际值与理论值的差异,可实现对工装热胀系数的现场标定,进而精确计算工装的热变形偏差;基于此,提出的基于比例缩放法的基准点修正方法,可适用于非均匀温度环境;此外,本发明没有引入其他昂贵设备,且对温度测量精度和温度场重建方法要求低,成本低廉,便于现场快速实施。本发明可应用于现场非均匀温度条件大型构件装配中测量基准与数模基准匹配、激光跟踪仪转站误差修正等,具有广阔的应用前景。
Claims (1)
1.一种基于热胀系数反求的测量基准偏差修正方法,其特征是,该方法首先对工装局部结构进行主动温控,并同步测量温度和基准点坐标;其次重建温度场函数,并基于温度和基准点坐标偏移间的映射关系反求最优热胀系数;然后在实际工况下,基于温度测量数据和反求的最优热胀系数,采用比例缩放法修正基准点坐标;最后以最小二乘法建立测量坐标系与数模坐标系间的转换关系,求解基准修正前后的转换误差;该方法的具体步骤如下:
第一步,主动温控下的温度-坐标同步测量
在工装上放置两个激光跟踪仪靶球,分别作为基准点A和基准点B,在两点之间布设m个温度传感器,然后使用加热装置对工装进行加热,通过温度传感器获取tj,j=1,2,…n时刻的温度数据集TEMPj:
TEMPj=[T(x1,tj) T(x2,tj)…T(xi,tj)…T(xm,tj)]T (1)
式中,xi,i=1,2,…m为温度传感器的位置,tj,j=1,2,…n为温度数据的采样时刻,T(xi,tj)为第i个温度传感器在tj时刻测的温度值;
温度传感器测量的同时,使用激光跟踪仪连续测量基准点A和基准点B的空间坐标,按照公式(2)计算加热过程中基准点A和基准点B间距离变化数据集DELT:
第二步,温度场重建与热胀系数反求
在tj时刻,求解式(3)所示的最小二乘方程,对温度数据进行拟合,重建此时刻的温度场函数f(x,tj):
tj时刻,计算在温度场f(x,tj)作用下基准点B相对基准点A的偏移量为:
式中,α为工装材料的热胀系数,为待求参数;
计算n个温度场函数与t1时刻温度场函数差的积分:
由公式(4)可知DELT与INT_TEMP间是线性关系,比例系数即为热胀系数,因此建立线性回归模型以求解最优热胀系数αbest:
DELT=αbest×INT_TEMP+b (6)
第三步,基于比例缩放法的基准点偏差修正
在飞机构件装配的实际工况下,假设温度传感器在tn+1时刻测量的工装温度值为:
TEMPn+1=[T(x1,tn+1) T(x2,tn+1)…T(xm,tn+1)] (7)
式中,αbest为根据公式(6)求得的最优热胀系数,f(x,tn+1)为tn+1时刻的温度场函数;Tnom为假设不发生热变形的参考温度值,为20℃;
Ln+1为基准点A和基准点B间的距离:
计算基准点B相对基准点A的比例缩放系数:
则修正后的基准点B坐标为:
第四步,测量基准与数模基准转换偏差计算
使用激光跟踪仪测量工装上N个基准点的坐标,获得坐标集S:
采用所述第一步~第三步的方法,获得修正后的坐标集Sc:
上述N个基准点在工装数模坐标系中的坐标集为:
建立如公式(15)所示的最小二乘模型,求解测量坐标系和数模坐标系间的基准转换参数:旋转矩阵R和平移向量T,并计算修正前的基准转换偏差E:
E=||D-(R×S+T)|| (16)
式中,S(k)和D(k)分别表示测量坐标集S和D中第k列向量;
同理,以公式(13)中修正后的坐标集Sc替换公式(15)中的坐标集S,计算测量基准修正后的基准转换偏差Ec:
Ec=||D-(Rc×Sc+Tc)|| (18)
由公式(18)计算出测量基准修正后的基准转换偏差Ec。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111569796.9A CN114370826B (zh) | 2021-12-21 | 2021-12-21 | 基于热胀系数反求的测量基准偏差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111569796.9A CN114370826B (zh) | 2021-12-21 | 2021-12-21 | 基于热胀系数反求的测量基准偏差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114370826A CN114370826A (zh) | 2022-04-19 |
CN114370826B true CN114370826B (zh) | 2022-09-20 |
Family
ID=81139597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111569796.9A Active CN114370826B (zh) | 2021-12-21 | 2021-12-21 | 基于热胀系数反求的测量基准偏差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114370826B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115115653A (zh) * | 2022-06-13 | 2022-09-27 | 广东众志检测仪器有限公司 | 一种冷热冲击试验箱精细化温度校准方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156410B (zh) * | 2016-06-24 | 2019-08-16 | 南京航空航天大学 | 薄壁圆管结构轴向热变形计算方法 |
JP2018128345A (ja) * | 2017-02-08 | 2018-08-16 | 地方独立行政法人東京都立産業技術研究センター | 温度補正方法、温度補正プログラム、及び座標測定機 |
CN110728088B (zh) * | 2019-09-27 | 2021-06-04 | 清华大学 | 工件三维热膨胀变形的跟踪仪转站参数优化方法及装置 |
CN111692976B (zh) * | 2020-06-08 | 2022-09-09 | 中国科学院合肥物质科学研究院 | 一种温度形变自补偿的数显长度基准装置 |
CN112729215B (zh) * | 2020-11-27 | 2022-06-14 | 成都飞机工业(集团)有限责任公司 | 一种基于热膨胀系数修正基准点坐标的测量方法 |
CN112668220B (zh) * | 2020-12-23 | 2022-12-09 | 天津大学 | 一种基于有限元分析的航天装置结构三维热变形测量方法 |
-
2021
- 2021-12-21 CN CN202111569796.9A patent/CN114370826B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114370826A (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105675143B (zh) | 一种真空黑体辐射源 | |
Zhong et al. | Position geometric error modeling, identification and compensation for large 5-axis machining center prototype | |
CN114370826B (zh) | 基于热胀系数反求的测量基准偏差修正方法 | |
CN105388937B (zh) | 一种用于气体传感器的精准恒温控制方法及装置 | |
CN112668220B (zh) | 一种基于有限元分析的航天装置结构三维热变形测量方法 | |
CN106885511A (zh) | 应用于量测的双球杆系统及其误差补偿方法 | |
CN107702758B (zh) | 一种基于热容法大型贮箱推进剂剩余量精确测量方法 | |
CN110631710B (zh) | 一种无挡片红外测温方法及装置 | |
CN115307866B (zh) | 风洞天平体轴系弹性角在线校准装置及方法 | |
CN114131595B (zh) | 一种机器人6d位姿估计系统及方法 | |
CN113483900A (zh) | 基于黑体点在线校准的红外辐射铝合金板温度场测量方法 | |
CN113465595A (zh) | 一种惯性定位定向设备惯组参数温补标定方法 | |
CN111581720B (zh) | 一种评定飞行器全方程控热试验温度数据不确定度的方法 | |
CN115962867A (zh) | 空间科学实验温度测量系统的数据采集标定方法及系统 | |
CN105136304B (zh) | 温度校准方法及系统 | |
CN106500641A (zh) | 关节式坐标测量机的热形变误差补偿方法 | |
CN113513999B (zh) | 一种航天器结构静力试验大视场实时变形测量系统及方法 | |
CN210038533U (zh) | 一种数控机床温升补偿系统 | |
CN114740798B (zh) | 一种数控设备群协同生产线测量场构建方法 | |
CN113405452B (zh) | 一种基于温度补偿的数字化工装飞机坐标系标定方法 | |
CN112629429B (zh) | 一种真空和变温环境下整星热变形测量装置、系统及方法 | |
CN108548839A (zh) | 一种高精度零部件热变形测试装置及测试方法 | |
CN111692976B (zh) | 一种温度形变自补偿的数显长度基准装置 | |
CN113405453A (zh) | 一种基于温度补偿的数字化工装飞机坐标系恢复方法 | |
Zhu et al. | Accuracy improving deformation measurement system for large components in thermal vacuum using close-range photogrammetry |
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 |