CN109202539A - 一种复合材料弱刚度异形结构在线检测方法 - Google Patents
一种复合材料弱刚度异形结构在线检测方法 Download PDFInfo
- Publication number
- CN109202539A CN109202539A CN201810966708.0A CN201810966708A CN109202539A CN 109202539 A CN109202539 A CN 109202539A CN 201810966708 A CN201810966708 A CN 201810966708A CN 109202539 A CN109202539 A CN 109202539A
- Authority
- CN
- China
- Prior art keywords
- point
- measurement
- interference
- measurement point
- 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.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/20—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring workpiece characteristics, e.g. contour, dimension, hardness
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种复合材料弱刚度异形结构在线检测方法,该检测方法首先对待测产品的测量点坐标进行采样;其次对采样后的测量点坐标进行优化,然后利用C空间理论对测量路径进行干涉碰撞检测,得到无干涉在在位检测路径;最后通过非线性最小二乘法进行误差溯源和分离,消除位姿误差对评定精度的影响。本发明能够完成检测路径的规划和生成测量点云、获得高效无干涉的在线检测路径,消除位姿误差对评定精度的影响。
Description
技术领域
本发明涉及一种面向异形零件的加工精度在线检测方法,尤其是针对复合材料弱刚度结构的精密零件的加工变形进行快速在线检测,属于面向复合材料高精密零件在线检测技术。
背景技术
随着材料科学和精密加工技术的不断发展,C-SiC等复合材料零件具有质量轻、寿命长和良好的热力学性能等优点被广泛应用于航空、航天等领域,其中许多高精度关键零件为了满足各类空间约束和特殊使用功能(高温、高压和剧烈摩擦),通常具有几何结构复杂、测量可达性差、薄壁和弱刚度等特点,在加工和测量过程中此类零件易受环境因素干扰,给精度质量等检测带来较大困难。
目前数控精密加工中应用的检测手段主要是基于三坐标测量机和专用测量仪器的离线检测,其主要局限为:(1)测量完成后,工件在机床进行二次装夹,导致原有定位信息丢失、弱刚度精密零件受力变形,难以保证后续加工精度,且增加工序辅助时间;(2)对于几何特征较多的复杂轮廓精密零件,检测量大,检测效率低;(3)对于空间几何结构复杂的精密零件,难以消除检测过程的***误差,缺少适用的误差评定方法。
发明内容
有鉴于此,本发明提供了一种复合材料弱刚度异形结构在线检测方法,该方法能够完成检测路径的规划和生成测量点云、获得高效无干涉的在线检测路径,消除位姿误差对评定精度的影响。
一种复合材料弱刚度异形结构在线检测方法,该检测方法的实现步骤如下:
步骤一:对待测产品的测量点坐标进行采样;
步骤二:对采样后的测量点坐标进行优化,使坐标更加准确;
步骤三:利用C空间理论对测量路径进行干涉碰撞检测,得到无干涉在位检测路径;
步骤四:通过非线性最小二乘法进行误差溯源和分离,消除位姿误差对评定精度的影响。
进一步地,所述步骤一中采样的过程如下:
第一步:对被测表面曲线进行参数化规范,使相邻两测量点间参数之差为常数;
第二步:在曲线上选取n个测量点,计算各点曲率半径;
第三步:计算相邻测量点之间弦差,其中,R表示相邻两测量点的曲率半径的平均值,L表示圆心O到弦AiAi+1的距离;
第四步:对测量点进行第一次筛选,如果弦差h不小于H,则将该采样点纳入测量点集;反之,则舍去该采样点,并对下一个采样点以同样方法进行判断;
第五步:对测量点进行第二次筛选,如果最大容许步长S大于初始设定值,将采样点纳入到测量点集;反之,则舍去该采样点,并对下一个采样点进行同样方法进行判断。
进一步地,所述步骤二通过齐次坐标变换生成测量点,得到可直接用于数控编程的测量位置数据,包括以下步骤:
第一步:根据设定的测量步长和步骤一中求得的弦差,得到测头触点MC;
第二步:将测头触点沿测量轴半径方向偏置一个位置矢量,得到测量位置点ML,即ML=MC+U;
其中,测量位置点坐标可以表示为:
进一步地,所述步骤三中得到无干涉在在位检测路径的过程如下:
第一步:测头B的任一形位q可用坐标系{B}相对于{A}的位置和姿态来表示,因此B的任一形位由位置矢量和旋转矩阵R两个参数表示,记为B(R);同时,将测头B在运动过程中所有不与其他刚体产生接触的形位称为其C空间模型;取单位球为参考空间,建立待测零件A与测头B在参考空间中的C空间模型,通过布尔求交运算,进行干涉碰撞检验;将待测零件表面分解为n个三角形面块,记为Fi,i∈[1,n]。面块Fi的形心记为Pi,i∈[1,n];
第二步:计算点Pi坐标,i∈[1,n],并计算该点处的法向
第三步:以参考坐标原点为球心建立单位球,对点Pi的法向量进行平移,使其延长方向通过单位球球心O,并与单位球表面交于点Qi,i∈[1,n];
第四步:在单位球表面以Qi为极点,得到n个半球面Si,i∈[1,n];
第五步:对半球面区域求交,得到零件的C空间模型A;
第六步:记测头轴线矢量为对所有形位q重复第三、四、五步,得到测头的C空间模型B;
第七步:通过A和B的布尔求交运算,进行干涉碰撞检验:如果则此测量路径中无干涉碰撞;否则,对重新规划测量路径并重复上述步骤,进而得到无干涉的测量路径。
进一步地,所述步骤四的误差溯源和分离过程如下:
第一步:被测线轮廓度误差与相应的被测曲线名义尺寸相比是微量,评定点N至理论曲线的距离其中fix,fiy为f(x,y)在点N处对x和y的一阶导数;
第二步:建立最小二乘法评定模;n为测量点数目;
第三步:引入位姿误差u=[dx,dy,θ]T,测量基准与理想基准之间的偏差和相应的曲线名义尺寸相比是微量,得出齐次坐标变换表达式,
第四步:根据误差评定模型优化目标函数,
采用非线性最小二乘法表示,
第五步:进行优化求解,在迭代点x(k)取目标函数f(x)二次近δ(k)为f(x)在点x(k)处的一个下降方向,在该搜索方向进行线性搜索,x(k+1)=x(k)+αkδ(k),并循环迭代,位姿误差逐渐趋向于0。
有益效果:
1、本发明针对检测过程中采样点之间的过大弦差,通过几何关系分析,综合考虑测头的尺寸参数和模型尺寸参数,进行检测路径规划,提高测量效率。
2、本发明能够有效避免在线检测过程中测头与零件的干涉碰撞问题,提高路径检测的效率和精度,得到无干涉碰撞的检测路径。
3、本发明在评定精密微小型零件的轮廓误差过程中,通过位姿误差溯源和分离,消除位姿误差对评定精度的影响。
附图说明
图1为本发明的在线检测方法流程图;
图2为对测量点坐标进行采样的示意图;
图3为对测量点坐标进行采样的流程图;
图4为零件干涉检验示意图;
图5为弱刚度零件在线检测结构图;
图6为零件干涉检验过程示意图
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种复合材料弱刚度异形结构在线检测方法,该方法的检测步骤如附图1所示,
步骤一:对待测产品的测量点坐标进行采样,以测量精度为前提,综合考虑测量步长、测头矢量等参数进行采样,包括以下内容:
1)对被测表面曲线进行参数化规范,使相邻两测量点间参数之差为常数;
2)在曲线上选取n个测量点,计算各点曲率半径;
3)计算相邻测量点之间弦差,其中,R表示相邻两测量点的曲率半径的平均值,L表示圆心O到弦AiAi+1的距离,如图2所示;
4)对测量点进行第一次筛选。如果弦差h不小于H,则将该采样点纳入测量点集;反之,则舍去该采样点,并对下一个采样点以同样方法进行判断;
5)对测量点进行第二次筛选。如果最大容许步长S大于初始设定值,将采样点纳入到测量点集;反之,则舍去该采样点,并对下一个采样点进行同样方法进行判断;完整的测量点采样流程如图3所示。
步骤二:通过齐次坐标变换生成测量点,得到可直接用于数控编程的测量位置数据。包括以下内容:
(1)根据设定的测量步长和步骤一中求得的弦差,得到测头触点MC;
(2)将测头触点沿测量轴半径方向偏置一个位置矢量,得到测量位置点ML,即ML=MC+U;
其中,测量位置点坐标可以表示为:
步骤三:基于C空间理论的对测量路径进行干涉碰撞检测,得到无干涉的在位检测路径,包括以下内容:
建立C空间模型C空间理论通常用于工程机构、切削刀具以及测量头等刚***置和姿态的数学表达。在C空间中,刚体结构被映射成一个点,干涉物体映射成C障碍物。通过建立刚体结构C空间模型,可以大大简化其运动规划的计算过程。
测头B的任一形位q可用坐标系{B}相对于{A}的位置和姿态来表示,因此B的任一形位由位置矢量和旋转矩阵R两个参数表示,记为B(R)。同时,将测头B在运动过程中所有不与其他刚体产生接触的形位称为其C空间模型。取单位球为参考空间,建立待测零件A与测头B在参考空间中的C空间模型,通过布尔求交运算,进行干涉碰撞检验。具体方法如下:
1)将待测零件表面分解为n个三角形面块,记为Fi,i∈[1,n]。面块Fi的形心记为Pi,i∈[1,n];
2)计算点Pi坐标,i∈[1,n],并计算该点处的法向
3)以参考坐标原点为球心建立单位球,对点Pi的法向量进行平移,使其延长方向通过单位球球心O,并与单位球表面交于点Qi,i∈[1,n];
4)在单位球表面以Qi为极点,得到n个半球面Si,i∈[1,n];
5)对半球面区域求交,得到零件的C空间模型A;
6)记测头轴线矢量为对所有形位q重复步骤3)、4)、5),得到测头的C空间模型B;
7)通过A和B的布尔求交运算,进行干涉碰撞检验:如果则此测量路径中无干涉碰撞;否则,对重新规划测量路径并重复上述步骤,进而得到无干涉的测量路径,干涉检测实例图如图4所示。
步骤四:通过非线性最小二乘法进行误差溯源和分离,消除位姿误差对评定精度的影响。
1)被测线轮廓度误差与相应的被测曲线名义尺寸相比是微量,评定点N至理论曲线的距离其中fix,fiy为f(x,y)在点N处对x和y的一阶导数;
2)建立最小二乘法评定模型,n为测量点数目;
3)引入位姿误差u=[dx,dy,θ]T,测量基准与理想基准之间的偏差和相应的曲线名义尺寸相比是微量,得出齐次坐标变换表达式,
4)根据误差评定模型优化目标函数,
采用非线性最小二乘法表示,
5)进行优化求解,在迭代点x(k)取目标函数f(x)二次近δ(k)为f(x)在点x(k)处的一个下降方向,在该搜索方向进行线性搜索,x(k+1)=x(k)+αkδ(k),并循环迭代,位姿误差逐渐趋向于0。
如附图5所示,弱刚度环类零件在线检测***包括数控机床、接触式测头、信号接受和处理装置(数控***、电脑及相应软件),在线检测***包括NC程序编写、干涉检验、测头测量、数据分析校验、加工变形误差分析和误差补偿(确定补偿量和建构零件模型)等过程,其流程图如图6所示。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种复合材料弱刚度异形结构在线检测方法,其特征在于,该检测方法的实现步骤如下:
步骤一:对待测产品的测量点坐标进行采样;
步骤二:对采样后的测量点坐标进行优化,使坐标更加准确;
步骤三:利用C空间理论对测量路径进行干涉碰撞检测,得到无干涉在在位检测路径;
步骤四:通过非线性最小二乘法进行误差溯源和分离,消除位姿误差对评定精度的影响。
2.如权利要求1所述的在线检测方法,其特征在于,所述步骤一中采样的过程如下:
第一步:对被测表面曲线进行参数化规范,使相邻两测量点间参数之差为常数;
第二步:在曲线上选取n个测量点,计算各点曲率半径;
第三步:计算相邻测量点之间弦差,其中,R表示相邻两测量点的曲率半径的平均值,L表示圆心O到弦AiAi+1的距离;
第四步:对测量点进行第一次筛选,如果弦差h不小于H,则将该采样点纳入测量点集;反之,则舍去该采样点,并对下一个采样点以同样方法进行判断;
第五步:对测量点进行第二次筛选,如果最大容许步长S大于初始设定值,将采样点纳入到测量点集;反之,则舍去该采样点,并对下一个采样点进行同样方法进行判断。
3.如权利要求2所述的在线检测方法,其特征在于,所述步骤二通过齐次坐标变换生成测量点,得到可直接用于数控编程的测量位置数据,包括以下步骤:
第一步:根据设定的测量步长和步骤一中求得的弦差,得到测头触点MC;
第二步:将测头触点沿测量轴半径方向偏置一个位置矢量,得到测量位置点ML,即ML=MC+U;
其中,测量位置点坐标可以表示为:
4.如权利要求3所述的在线检测方法,其特征在于,所述步骤三中得到无干涉在在位检测路径的过程如下:
第一步:测头B的任一形位q可用坐标系{B}相对于{A}的位置和姿态来表示,因此B的任一形位由位置矢量和旋转矩阵R两个参数表示,记为同时,将测头B在运动过程中所有不与其他刚体产生接触的形位称为其C空间模型;取单位球为参考空间,建立待测零件A与测头B在参考空间中的C空间模型,通过布尔求交运算,进行干涉碰撞检验;将待测零件表面分解为n个三角形面块,记为Fi,i∈[1,n]。面块Fi的形心记为Pi,i∈[1,n];
第二步:计算点Pi坐标,i∈[1,n],并计算该点处的法向
第三步:以参考坐标原点为球心建立单位球,对点Pi的法向量进行平移,使其延长方向通过单位球球心O,并与单位球表面交于点Qi,i∈[1,n];
第四步:在单位球表面以Qi为极点,得到n个半球面Si,i∈[1,n];
第五步:对半球面区域求交,得到零件的C空间模型A;
第六步:记测头轴线矢量为对所有形位q重复第三、四、五步,得到测头的C空间模型B;
第七步:通过A和B的布尔求交运算,进行干涉碰撞检验:如果则此测量路径中无干涉碰撞;否则,对重新规划测量路径并重复上述步骤,进而得到无干涉的测量路径。
5.如权利要求4所述的在线检测方法,其特征在于,所述步骤四的误差溯源和分离过程如下:
第一步:被测线轮廓度误差与相应的被测曲线名义尺寸相比是微量,评定点N至理论曲线的距离其中fix,fiy为f(x,y)在点N处对x和y的一阶导数;
第二步:建立最小二乘法评定模;n为测量点数目;
第三步:引入位姿误差u=[dx,dy,θ]T,测量基准与理想基准之间的偏差和相应的曲线名义尺寸相比是微量,得出齐次坐标变换表达式,
第四步:根据误差评定模型优化目标函数,采用非线性最小二乘法表示,
第五步:进行优化求解,在迭代点x(k)取目标函数f(x)二次近δ(k)为f(x)在点x(k)处的一个下降方向,在该搜索方向进行线性搜索,x(k+1)=x(k)+αkδ(k),并循环迭代,位姿误差逐渐趋向于0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810966708.0A CN109202539B (zh) | 2018-08-23 | 2018-08-23 | 一种复合材料弱刚度异形结构在线检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810966708.0A CN109202539B (zh) | 2018-08-23 | 2018-08-23 | 一种复合材料弱刚度异形结构在线检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109202539A true CN109202539A (zh) | 2019-01-15 |
CN109202539B CN109202539B (zh) | 2020-10-30 |
Family
ID=64989603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810966708.0A Active CN109202539B (zh) | 2018-08-23 | 2018-08-23 | 一种复合材料弱刚度异形结构在线检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109202539B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059879A (zh) * | 2019-04-16 | 2019-07-26 | 上海理工大学 | 车身三坐标测量的自动规划方法 |
CN113884042A (zh) * | 2021-10-26 | 2022-01-04 | 上海理工大学 | 面向电池包的三坐标测量自动规划方法 |
CN114115123A (zh) * | 2021-11-16 | 2022-03-01 | 上海交通大学 | 航空大型薄壁类非刚体零件的参数化数控加工方法和*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0736527A (ja) * | 1993-06-25 | 1995-02-07 | Nikon Corp | 形状測定方法および装置 |
CN101342664A (zh) * | 2008-08-20 | 2009-01-14 | 天津大学 | 一种复杂空间型面在机质量检测*** |
CN101898317A (zh) * | 2010-06-10 | 2010-12-01 | 重庆大学 | 一种虚拟式数控机床在线检测***和方法 |
CN102785129A (zh) * | 2012-07-30 | 2012-11-21 | 广东工业大学 | 复杂零件的曲面加工精度的在线检测方法 |
CN102814512A (zh) * | 2012-08-24 | 2012-12-12 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种针对发动机压气机盘类零件辐板型面的在线测量方法 |
EP2570236A1 (en) * | 2011-09-15 | 2013-03-20 | Ideko, S. Coop | Method of measuring and aligning parts for machining in a machine tool |
-
2018
- 2018-08-23 CN CN201810966708.0A patent/CN109202539B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0736527A (ja) * | 1993-06-25 | 1995-02-07 | Nikon Corp | 形状測定方法および装置 |
CN101342664A (zh) * | 2008-08-20 | 2009-01-14 | 天津大学 | 一种复杂空间型面在机质量检测*** |
CN101898317A (zh) * | 2010-06-10 | 2010-12-01 | 重庆大学 | 一种虚拟式数控机床在线检测***和方法 |
EP2570236A1 (en) * | 2011-09-15 | 2013-03-20 | Ideko, S. Coop | Method of measuring and aligning parts for machining in a machine tool |
CN102785129A (zh) * | 2012-07-30 | 2012-11-21 | 广东工业大学 | 复杂零件的曲面加工精度的在线检测方法 |
CN102814512A (zh) * | 2012-08-24 | 2012-12-12 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种针对发动机压气机盘类零件辐板型面的在线测量方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059879A (zh) * | 2019-04-16 | 2019-07-26 | 上海理工大学 | 车身三坐标测量的自动规划方法 |
CN110059879B (zh) * | 2019-04-16 | 2023-08-18 | 上海理工大学 | 车身三坐标测量的自动规划方法 |
CN113884042A (zh) * | 2021-10-26 | 2022-01-04 | 上海理工大学 | 面向电池包的三坐标测量自动规划方法 |
CN113884042B (zh) * | 2021-10-26 | 2023-07-25 | 上海理工大学 | 面向电池包的三坐标测量自动规划方法 |
CN114115123A (zh) * | 2021-11-16 | 2022-03-01 | 上海交通大学 | 航空大型薄壁类非刚体零件的参数化数控加工方法和*** |
CN114115123B (zh) * | 2021-11-16 | 2024-04-09 | 上海交通大学 | 航空大型薄壁类非刚体零件的参数化数控加工方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN109202539B (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107315391B (zh) | 一种数控机床在线检测的预行程误差补偿方法 | |
CN102785128B (zh) | 面向数控车床的零件加工精度在线检测***及检测方法 | |
US5111590A (en) | Measuring method of machine tool accuracy using a computer aided kinematic transducer link and its apparatus | |
Maisano et al. | Indoor GPS: system functionality and initial performance evaluation | |
CN102785129B (zh) | 复杂零件的曲面加工精度的在线检测方法 | |
CN103777570B (zh) | 基于nurbs曲面的加工误差快速检测补偿方法 | |
Uekita et al. | On-machine dimensional measurement of large parts by compensating for volumetric errors of machine tools | |
CN109202539A (zh) | 一种复合材料弱刚度异形结构在线检测方法 | |
CN102927952A (zh) | 一种航空发动机机匣锥形外壁直径在线检测方法 | |
CN104985482B (zh) | 一种五轴加工中心在机检测复杂型面方法 | |
CN110076631A (zh) | 复杂薄壁结构零件壁厚在机测量方法 | |
Stepien | In situ measurement of cylindricity—Problems and solutions | |
CN104166373A (zh) | 数控机床加工误差在线检测方法及*** | |
Manikandan et al. | Modelling of dimensional and geometric error prediction in turning of thin-walled components | |
CN111797482B (zh) | 用于笛卡尔坐标系轮廓扫描测量的敏感误差权重计算方法 | |
CN202656009U (zh) | 面向数控车床的零件加工精度在线检测*** | |
Li et al. | A high-speed in situ measuring method for inner dimension inspection | |
CN102269581A (zh) | 一种建筑构件加工精度测评方法 | |
CN101498572A (zh) | 凸轮二维曲线测量方法 | |
Zhuang et al. | A state-of-the-art review on the research and application of on-machine measurement with a touch-trigger probe | |
CN112114557B (zh) | 五轴联动数控机床动态精度检测方法、***及储存介质 | |
CN104596461A (zh) | 用于检测三轴金刚石车床定位精度的特征样件及方法 | |
Zhuang et al. | Uncertainty propagation and assessment of five-axis on-machine measurement error | |
Hu et al. | Precision measurement method of “cradle-type” five-axis machining center | |
Papananias et al. | Evaluation of automated flexible gauge performance using experimental designs |
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 |