CN110132195B - 叶片截面接触式扫描测量的探针测球三维半径补偿方法 - Google Patents
叶片截面接触式扫描测量的探针测球三维半径补偿方法 Download PDFInfo
- Publication number
- CN110132195B CN110132195B CN201910461858.0A CN201910461858A CN110132195B CN 110132195 B CN110132195 B CN 110132195B CN 201910461858 A CN201910461858 A CN 201910461858A CN 110132195 B CN110132195 B CN 110132195B
- Authority
- CN
- China
- Prior art keywords
- point
- vector
- discrete
- plane
- matrix
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005259 measurement Methods 0.000 title claims abstract description 17
- 239000000523 sample Substances 0.000 title claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000001514 detection method Methods 0.000 description 13
- 239000010979 ruby Substances 0.000 description 4
- 229910001750 ruby Inorganic materials 0.000 description 4
- 238000005299 abrasion Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- A Measuring Device Byusing Mechanical Method (AREA)
Abstract
Description
技术领域
本发明涉及三坐标测量机测量技术领域,尤其是一种叶片截面接触式扫描测量的探针测球三维半径补偿方法。
背景技术
叶片在航空发动机、涡轮机、风力发电等领域有重要应用,同时对叶片的加工及检测精度要求日益增高。在当前叶片检测中主要使用接触式和非接触式两种检测方法,前者检测精度较后者有几十微米的误差,只能对一些精度要求等级低或大工件应用。对高精度叶片检测,三坐标测量机是首选工具,在扫描过程中为降低测针磨损,在测针尖端增加了抗摩损的红宝石球作为测球,由于三坐标测量机读取的是红宝石球中心的坐标值,在降低磨损的同时带来了红宝石球半径补偿的问题。
当前测头测针半径补偿方法主要分为曲面补偿和曲线补偿。曲面补偿需要测量叶身数值,再向内偏置红宝石球半径,这就需要使用五轴测头,摆扫测量,这种测量方式不仅数据量大,且测头本身价格昂贵,所以目前叶片检测多数使用截面线扫描检测,并进行曲线补偿。曲线补偿方法是一种二维补偿,没有考虑叶片扭转产生的余弦角带来的余弦误差。针对此问题,专利CN 104330068 A提出一种使用截面扫描方式,测量多个截面线并拟合测量曲面,再按照表面法矢偏置的方法,再截取偏置后截面线。此方法理论上可以减小余弦误差,但需要测量多个截面,测量效率难以保障。
发明内容
本发明的目的在于克服现有技术中存在的不足,提供一种叶片截面接触式扫描测量的探针测球三维半径补偿方法,以解决叶片截面扫描在半径补偿时出现的余弦误差。本发明采用的技术方案是:
一种叶片截面接触式扫描测量的探针测球三维半径补偿方法,包括如下步骤:
步骤S101、以待测叶片坐标系为本方法使用坐标系XYZ,将三坐标测量的每个测球中心点M作为NURBS曲线的型值点P,利用累积弦长法计算NURBS曲线的节点矢量K,并计算NURBS曲线的基函数矩阵Ni,m;根据型值点矩阵反算控制顶点矩阵D,使拟合的NURBS曲线过型值点;
步骤S102、等距离散NURBS曲线,每个离散点表示为点O,计算离散点平行于XY平面的法向矢量Nxy;
步骤S103、离散点O沿Z轴方向在叶片截面所在平面的投影点为O',计算向量A点为测球与叶片实际接触点,令向量在XY平面的方向等于Nxy的方向,向量在Z轴方向的分量为0,则与所成的补偿向量即为补偿三维半径补偿向量。
进一步地,步骤S101具体包括:
以待测叶片坐标系为本方法使用坐标系XYZ,待测叶片截面平行于X/Y轴所在平面即XY平面,待测叶片截面所在平面垂直于Z轴;
将三坐标测量的每个测球中心点M作为型值点P,共有n+1个型值点,NURBS曲线的次数m;
NURBS曲线的方程C表达为:
其中,k为变量,di为控制顶点,wi为权重,K=[k0,…,ki,…,kn+m+1]为NURBS曲线的节点矢量,可表达为:
Pi为型值点坐标;
公式(1)中,Ni,m为NURBS曲线的基函数矩阵,表达为公式(3),或公式(5)的系数矩阵;
控制顶点矩阵D由公式(1)中的控制顶点表达为:
D=[d0 … dn]T (4)
公式(4)由公式(5)求解,并且公式(5)的系数矩阵是公式(3)的一种简便表达;
进一步地,步骤S102具体包括:
等距离散NURBS曲线,共离散l个点,离散距离通过离散点数量确定,离散点记为Ci(k),i∈[1,l],每个离散点表示为点O,所有离散点沿Z向投影到叶片截面所在平面,在叶片截面所在平面内计算每个离散点的切向矢量:
Ti(k)=Ci'(k) (6)
其中,Ci'(k)为Ci(k)的导数,则平行于XY平面的法向矢量Nxy表达为:
法向矢量Nxy表示为Nxy=(nx,ny),nx表示法向矢量在X轴方向的分量,ny表示法向矢量在Y轴方向的分量。
其中,z为离散点的Z轴方向数值,z0为所测截面的Z轴方向数值,r为探针测球的物理半径。
本发明的优点:本发明提出的叶片截面接触式扫描测量的探针测球三维半径补偿方法,利用探针扫描截面时测球球心位置与截面在Z轴方向存在的一个高度差,在二维补偿的基础上增加了这个高度信息,提高了叶片截面扫描中的半径补偿精度。
附图说明
图1为本发明的流程图。
图2为本发明的补偿向量示意图。
具体实施方式
下面结合具体附图和实施例对本发明作进一步说明。
本发明提出一种叶片截面接触式扫描测量的探针测球三维半径补偿方法,包括如下步骤:
步骤S101、以待测叶片坐标系为本方法使用坐标系XYZ,将三坐标测量的每个测球中心点M作为NURBS(非均匀有理B样条,Non-Uniform Rational B-Splines)曲线的型值点P,利用累积弦长法计算NURBS曲线的节点矢量K,并计算NURBS曲线的基函数矩阵Ni,m;根据型值点矩阵反算控制顶点矩阵D,使拟合的NURBS曲线过型值点;
具体地,
以待测叶片坐标系为本方法使用坐标系XYZ,待测叶片截面平行于X/Y轴所在平面即XY平面,待测叶片截面所在平面垂直于Z轴;
将三坐标测量的每个测球中心点M作为型值点P,共有n+1个型值点,本发明使用3次NURBS曲线,即NURBS曲线的次数m取3;
NURBS曲线的方程C表达为:
其中,k为变量,di为控制顶点,wi为权重,本发明权重取1,K=[k0,…,ki,…,kn+m+1]为NURBS曲线的节点矢量,可表达为:
Pi为型值点坐标;
公式(1)中,Ni,m为NURBS曲线的基函数矩阵,可表达为公式(3),其更简便的表达方式为公式(5)的系数矩阵;
控制顶点矩阵D由公式(1)中的控制顶点表达为:
D=[d0 … dn]T (4)
公式(4)由公式(5)求解,并且公式(5)的系数矩阵是公式(3)的一种简便表达;
步骤S102、等距离散NURBS曲线,每个离散点表示为点O,计算离散点平行于XY平面的法向矢量Nxy;
具体包括:
等距离散NURBS曲线,共离散l个点,离散距离通过离散点数量确定,离散点记为Ci(k),i∈[1,l],每个离散点表示为点O,所有离散点沿Z向投影到叶片截面所在平面,在叶片截面所在平面内计算每个离散点的切向矢量:
Ti(k)=Ci'(k) (6)
其中,Ci'(k)为Ci(k)的导数,则平行于XY平面的法向矢量Nxy可表达为:
法向矢量Nxy表示为Nxy=(nx,ny),nx表示法向矢量在X轴方向的分量,ny表示法向矢量在Y轴方向的分量;
步骤S103、点O沿Z轴方向在XY平面的投影点为O',计算向量A点为测球与叶片实际接触点,令向量在XY平面的方向等于Nxy的方向,向量在Z轴方向的分量为0,则与所成的补偿向量即为补偿三维半径补偿向量。如图2所示。
步骤S103、离散点O沿Z轴方向在叶片截面所在平面的投影点为O',计算向量A点为测球与叶片实际接触点,令向量在XY平面的方向等于Nxy的方向,向量在Z轴方向的分量为0,则与所成的补偿向量即为补偿三维半径补偿向量。如图2所示。
具体包括:
其中,z为离散点的Z轴方向数值,z0为所测截面的Z轴方向数值,r为探针测球的物理半径。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种叶片截面接触式扫描测量的探针测球三维半径补偿方法,其特征在于,包括如下步骤:
步骤S101、以待测叶片坐标系为本方法使用坐标系XYZ,将三坐标测量的每个测球中心点M作为NURBS曲线的型值点P,利用累积弦长法计算NURBS曲线的节点矢量K,并计算NURBS曲线的基函数矩阵Ni,m;根据型值点矩阵反算控制顶点矩阵D,使拟合的NURBS曲线过型值点;
步骤S102、等距离散NURBS曲线,每个离散点表示为点O,计算离散点平行于XY平面的法向矢量Nxy;
步骤S103、离散点O沿Z轴方向在叶片截面所在平面的投影点为O',计算向量A点为测球与叶片实际接触点,令向量在XY平面的方向等于Nxy的方向,向量在Z轴方向的分量为0,则与所成的补偿向量即为补偿三维半径补偿向量;
步骤S101具体包括:
以待测叶片坐标系为本方法使用坐标系XYZ,待测叶片截面平行于X/Y轴所在平面即XY平面,待测叶片截面所在平面垂直于Z轴;
将三坐标测量的每个测球中心点M作为型值点P,共有n+1个型值点,NURBS曲线的次数m;
NURBS曲线的方程C表达为:
其中,k为变量,di为控制顶点,wi为权重,K=[k0,…,ki,…,kn+m+1]为NURBS曲线的节点矢量,可表达为:
Pi为型值点坐标;
公式(1)中,Ni,m为NURBS曲线的基函数矩阵,表达为公式(3),或公式(5)的系数矩阵;
控制顶点矩阵D由公式(1)中的控制顶点表达为:
D=[d0…dn]T (4)
公式(4)由公式(5)求解,并且公式(5)的系数矩阵是公式(3)的一种简便表达;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910461858.0A CN110132195B (zh) | 2019-05-30 | 2019-05-30 | 叶片截面接触式扫描测量的探针测球三维半径补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910461858.0A CN110132195B (zh) | 2019-05-30 | 2019-05-30 | 叶片截面接触式扫描测量的探针测球三维半径补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110132195A CN110132195A (zh) | 2019-08-16 |
CN110132195B true CN110132195B (zh) | 2021-04-20 |
Family
ID=67582874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910461858.0A Active CN110132195B (zh) | 2019-05-30 | 2019-05-30 | 叶片截面接触式扫描测量的探针测球三维半径补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110132195B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112033338B (zh) * | 2020-09-14 | 2021-12-28 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种叶片类曲面接触式扫描测量测头半径面补偿方法 |
CN117606327B (zh) * | 2024-01-24 | 2024-04-05 | 天津工业大学 | 一种航空发动机叶片型面检测设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4612622A (en) * | 1983-07-12 | 1986-09-16 | United Technologies Corporation | Probe for coordinate measuring machine |
CN103777570A (zh) * | 2014-01-07 | 2014-05-07 | 浙江大学 | 基于nurbs曲面的加工误差快速检测补偿方法 |
CN106202822A (zh) * | 2016-07-27 | 2016-12-07 | 西北工业大学 | 面向叶片自适应加工的b样条曲面模型重构方法 |
CN106393106A (zh) * | 2016-09-13 | 2017-02-15 | 东南大学 | 参数自适应密化的机器人nurbs曲线运动插补方法 |
CN107238364A (zh) * | 2017-06-30 | 2017-10-10 | 四川大学 | 接触式测量测针球头半径精密补偿方法 |
CN108106522A (zh) * | 2017-11-29 | 2018-06-01 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种非规则曲面的三维测量方法 |
CN108692644A (zh) * | 2018-03-26 | 2018-10-23 | 华中科技大学 | 一种复杂曲面三坐标测量装置及误差补偿方法 |
CN109341632A (zh) * | 2018-11-29 | 2019-02-15 | 株洲中航动力精密铸造有限公司 | 用于涡轮叶片叶型测量的测头及三坐标测量仪 |
-
2019
- 2019-05-30 CN CN201910461858.0A patent/CN110132195B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4612622A (en) * | 1983-07-12 | 1986-09-16 | United Technologies Corporation | Probe for coordinate measuring machine |
CN103777570A (zh) * | 2014-01-07 | 2014-05-07 | 浙江大学 | 基于nurbs曲面的加工误差快速检测补偿方法 |
CN106202822A (zh) * | 2016-07-27 | 2016-12-07 | 西北工业大学 | 面向叶片自适应加工的b样条曲面模型重构方法 |
CN106393106A (zh) * | 2016-09-13 | 2017-02-15 | 东南大学 | 参数自适应密化的机器人nurbs曲线运动插补方法 |
CN107238364A (zh) * | 2017-06-30 | 2017-10-10 | 四川大学 | 接触式测量测针球头半径精密补偿方法 |
CN108106522A (zh) * | 2017-11-29 | 2018-06-01 | 中国航发沈阳黎明航空发动机有限责任公司 | 一种非规则曲面的三维测量方法 |
CN108692644A (zh) * | 2018-03-26 | 2018-10-23 | 华中科技大学 | 一种复杂曲面三坐标测量装置及误差补偿方法 |
CN109341632A (zh) * | 2018-11-29 | 2019-02-15 | 株洲中航动力精密铸造有限公司 | 用于涡轮叶片叶型测量的测头及三坐标测量仪 |
Non-Patent Citations (2)
Title |
---|
复杂型面测量中测头中心的轨迹曲面及测头半径的三维补偿;王建华等;《计量学报》;19940430;第15卷(第2期);正文第1-3节 * |
自由曲面的二维自适应测量及测球半径的三维补偿;刘荷辉等;《机械工程学报》;20040229;第40卷(第2期);正文第1-2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN110132195A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110108208B (zh) | 五轴非接触式测量机的误差补偿方法 | |
JP4660779B2 (ja) | 移動装置の位置誤差評価方法およびその評価結果に基づく移動精度向上方法 | |
CN110108207B (zh) | 基于探针的旋转轴回转中心线几何误差标定方法 | |
CN110500978A (zh) | 点激光传感器的光束方向矢量和零点位置在线标定方法 | |
CN110132195B (zh) | 叶片截面接触式扫描测量的探针测球三维半径补偿方法 | |
CN106406237B (zh) | 一种具有自由曲面金属零件的加工方法 | |
CN102873586B (zh) | 数控加工工件曲率半径快速在线测量装置 | |
CN114012585B (zh) | 一种双摆轴式五轴磁流变机床抛光点位置标定方法 | |
CN106705880B (zh) | 一种大口径反射镜面形轮廓在位检测方法及装置 | |
Li et al. | Monocular-vision-based contouring error detection and compensation for CNC machine tools | |
CN111536877A (zh) | 一种三坐标测量机上的线激光传感器姿态标定方法 | |
CN115682989A (zh) | 一种涡轮叶片基于六点定位的形面测量方法 | |
CN115319727A (zh) | 一种基于位姿约束和力感知的机器人标定方法 | |
Moona et al. | Measurement uncertainty assessment of articulated arm coordinate measuring machine for length measurement errors using Monte Carlo simulation | |
Guo et al. | Continuous measurements with single setup for position-dependent geometric errors of rotary axes on five-axis machine tools by a laser displacement sensor | |
CN113916128A (zh) | 一种基于光笔式视觉测量***的提高精度的方法 | |
CN117053688A (zh) | 一种线激光齿轮测量中心的空间误差建模方法及测量中心 | |
CN111895921A (zh) | 一种提高***对高度差测量精度的补偿模型 | |
CN110186405B (zh) | 叶片截面接触式扫描探针测球三维半径补偿及交叉补偿点纠正方法 | |
Wang et al. | A novel 3D radius compensation method of probe stylus tip in the free-form surface profile curve scanning measurement | |
Guo et al. | Contact on-machine measurement probe error correction method for optical aspheric surface ultraprecision machining | |
Harris et al. | Geometric and quasi-static thermal error compensation for a laser digitizer equipped coordinate measuring machine | |
CN112595280B (zh) | 一种成角度复杂面形测量方法 | |
CN114549521A (zh) | 基于加工g代码引导的碳管阵蜂窝面形精度计算方法 | |
CN113427320A (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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190816 Assignee: Jiangsu Jihui Huake Intelligent Equipment Technology Co., Ltd. Assignor: Wuxi research institute of the Central China University of Science and Technology Contract record no.: X2019980000288 Denomination of invention: The probe of blade profile contact type scanning measurement surveys ball three-dimensional radius compensation method License type: Common License Record date: 20191023 |
|
GR01 | Patent grant | ||
GR01 | Patent grant |