CN102853780A - 抛光工件表面轮廓各频段误差的分离方法 - Google Patents
抛光工件表面轮廓各频段误差的分离方法 Download PDFInfo
- Publication number
- CN102853780A CN102853780A CN2012103315733A CN201210331573A CN102853780A CN 102853780 A CN102853780 A CN 102853780A CN 2012103315733 A CN2012103315733 A CN 2012103315733A CN 201210331573 A CN201210331573 A CN 201210331573A CN 102853780 A CN102853780 A CN 102853780A
- Authority
- CN
- China
- Prior art keywords
- error
- frequency
- curve
- rank
- imf
- 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.)
- Pending
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
抛光工件表面轮廓各频段误差的分离方法,涉及光学表面加工方法。1)对采用光学检测装置检测得到的被测光学元件表面的检测曲线进行拟合,得到拟合曲线和残差;2)将残差进行经验模态分解,得到各阶固有模态函数和剩余信号;3)计算各阶固有模态函数的瞬时频率曲线及其平均波长;4)根据各阶固有模态函数的平均波长及其特征将所有固有模态函数分为高、中、低频组,将高频和中频组进行组内叠加得到光学表面轮廓的高频和中频误差;将低频组、剩余信号及拟合曲线叠加得到包含低频误差的面形。
Description
技术领域
本发明涉及光学表面加工方法,尤其是涉及抛光工件表面轮廓各频段误差的分离方法。
背景技术
小磨头确定性抛光是目前加工光学非球面的可行方法,目的是通过采用比光学元件尺寸小很多的加工工具进行定区域修研,可以大幅度提高加工精度,但同时会导致光学元件表面形成各频段误差分布不均匀的现象普遍存在,在小工具停留时间长的区域中频误差含量较大,停留时间短或不停留的区域,中频误差含量相对较小。光学元件的中频误差会导致的光束高频调制和非线性增益,将造成光学元件的丝状破坏和降低光束的可聚焦功率,对极紫外或新兴的X射线光学而言,中频段误差产生的小角度散射会极大地降低光学系统的分辨率。
为确保高能量激光系统和高分辨率系统的性能,需要对光学元件表面轮廓各频段误差进行分离,识别中频误差频率和方位,指导补偿加工。
光学元件波面数据的功率谱密度(Power Spectral Density,PSD)是目前国际上广泛采用的衡量光学元件中频误差的方法。但PSD基于傅里叶变换,会将局部瞬时波动平均到全局范围内,是一种综合评价手段,无法实现误差频率和方位的准确识别。
申请号为200810030817.8的中国专利申请公开一种确定性光学加工条件下频段误差分布特性的分析方法,它首先基于面形误差数据的PSD选择敏感频段,然后选择基本小波对加工误差敏感频段进行分析,确定误差分布区域。该方法存在不足之处:
其一:对面形误差数据计算PSD之后,将局部瞬时波动平均到全局范围内,会严重降低信噪比,忽略重要信息;
其二:敏感频段选择和小波基函数的选择对分析结果的影响较大。
发明内容
本发明的目的在于克服现有误差评价手段和识别方法的不足,提供抛光工件表面轮廓各频段误差的分离方法。
该方法可以自适应地分离光学元件表面轮廓的高、中、低频误差,识别中频误差的频率和方位,为确定性抛光提供指导。
本发明所述抛光工件表面轮廓各频段误差的分离方法,包括以下步骤:
1)对采用光学检测装置检测得到的被测光学元件表面的检测曲线进行拟合,得到拟合曲线和残差;
2)将残差进行经验模态分解,得到各阶固有模态函数和剩余信号;
3)计算各阶固有模态函数的瞬时频率曲线及其平均波长;
4)根据各阶固有模态函数的平均波长及其特征将所有固有模态函数分为高、中、低频组,将高频和中频组进行组内叠加得到光学表面轮廓的高频和中频误差;将低频组、剩余信号及拟合曲线叠加得到包含低频误差的面形。
本发明是一种抛光工件表面轮廓各频段误差的分离方法,它对面形拟合残差进行经验模态分解(Empirical mode decomposition,EMD),得到一系列固有模态函数(Intrinsic modefunctions,IMF),根据各阶IMF特征,识别不同空间位置存在的误差及其波动频率,并将其合成得到光学表面轮廓的高、中、低频误差。与现有技术相比,本发明具有如下突出优点:
1、本发明将光学元件表面的各频段误差进行分离,便于更加准确的评价高端光学元件表面质量;
2、利用EMD技术将面形拟合残差分解为一系列单分量的IMF,根据各阶IMF能够准确判断中频误差分布方位及频率,从而指导补偿加工和识别误差来源;
3、对光学表面轮廓的高、中、低频误差的分离采取基于空间信号自身局部特征的自适应滤波方式,排除普通滤波中频段选择对结果的影响。
附图说明
图1是本发明实施例的第9阶IMF。横坐标为扫描轨迹X(mm),纵坐标为第9阶IMF的幅值imf(mm)。
图2是本发明实施例的第9阶IMF的瞬时频率曲线。横坐标为扫描轨迹X(mm),纵坐标为第9阶IMF的瞬时频率infreq(1/mm)。
图3是本发明实施例的第10阶IMF。横坐标为扫描轨迹X(mm),纵坐标为第10阶IMF的幅值imf(mm)。
图4是本发明实施例的10阶第IMF的瞬时频率曲线。横坐标为扫描轨迹X(mm),纵坐标为第10阶IMF的瞬时频率infreq(1/mm)。
图5是本发明实施例的第11阶IMF。横坐标为扫描轨迹X(mm),纵坐标为第11阶IMF的幅值imf(mm)。
图6是本发明实施例的第11阶IMF的瞬时频率曲线。横坐标为扫描轨迹X(mm),纵坐标为第11阶IMF的瞬时频率infreq(1/mm)。
图7是本发明实施例的第12阶IMF。横坐标为扫描轨迹X(mm),纵坐标为第12阶IMF的幅值imf(mm)。
图8是本发明实施例的第12阶IMF的瞬时频率曲线。横坐标为扫描轨迹X(mm),纵坐标为第12阶IMF的瞬时频率infreq(1/mm)。
图9是本发明实施例的第13阶IMF。横坐标为扫描轨迹X(mm),纵坐标为第13阶IMF的幅值imf(mm)。
图10是本发明实施例的第13阶IMF的瞬时频率曲线。横坐标为扫描轨迹X(mm),纵坐标为第12阶IMF的瞬时频率infreq(1/mm)。
图11是本发明实施例的光学元件表面轮廓低频误差分离结果。横坐标为扫描轨迹X(mm),纵坐标为包含低频误差Z(mm)。
图12是本发明实施例的光学元件表面轮廓中频误差分离结果。横坐标为扫描轨迹X(mm),纵坐标为包含中频误差Z(mm)。
图13是本发明实施例的光学元件表面轮廓高频误差分离结果。横坐标为扫描轨迹X(mm),纵坐标为包含高频误差Z(mm)。
具体实施方式
本发明实施例是对轴对称非球面光学元件表面轮廓各频段误差进行分离:包括以下步骤:
1)对采用光学检测装置检测得到的被测光学元件表面的检测曲线进行拟合,得到拟合曲线和残差;
所述拟合是采用下述光学元件子午线通用方程进行拟合;
式中,r为顶点曲率半径,e为偏心率,ai(i=1~6)为高次项系数,x,z分别为子午线的横、纵坐标。通过逐步拟合的方式求取各系数,具体拟合步骤是首先将上述公式中的高次项置零,根据实际曲线数据最小二乘拟合得到顶点曲率半径和偏心率,再逐个添加高次项改变拟合公式,并将上次拟合求得的参数作为初始赋值进行再次拟合;
将拟合得到的各系数及原始曲线数据代入拟合公式得到拟合曲线c0,在原始曲线数据中减去c0得到残差s;
2)将残差s进行经验模态分解(EMD分解),具体步骤如下:
(1)找到s所有极大值点和极小值点,利用三次样条拟合得到s的上、下包络线smax和smin,并计算它们的平均值曲线m11;
在s中减去m11得到h11;将h11代替s重复上面的步骤,重复迭代:
h1k(x)=h1(k-1)(x)-m1k(x)k≥2
直到满足:
此时h1k是第一阶IMF:c1;
(2)从残差s中减去c1得到第一阶剩余信号r1,将r1代替s重复步骤(1),依次得到二阶IMF,直到n阶IMF和剩余信号rn,则残差s可以表示为各阶IMF和剩余信号之和的形式:
根据上述步骤完成残差s的经验模态分解,共得到16阶IMF,其中有代表性的第9~13阶IMF分别如图1、图3、图5、图7、图9所示。第9阶及以下IMF为无规律的高频波动。第10阶IMF在70mm~80mm存在某些局部波动,除此之外无显著或有规律的波动形式。第11阶IMF在70mm~100mm段有显著波动,说明在该范围内存在显著的表面误差。第12阶IMF在35mm~70mm段有小幅波动,说明该区域也包含某种频率的表面误差,但是幅度相比70mm~100mm段较小。第13阶IMF基本属于无规律小幅波动,13阶以上IMF为低频无规律小幅波动。
3)计算IMF的瞬时频率曲线,首先求各阶IMF的解析信号:
根据上述公式计算各阶IMF的瞬时频率曲线,第9~13阶IMF的瞬时频率曲线分别如图2、图4、图6、图8、图10所示。求各阶IMF瞬时频率曲线的平均值:
其中符号<·>x表示对所有x求均值。则第i阶IMF的平均波长为:
根据上述过程计算得到所有IMF的平均波长如表1所示,各阶IMF的平均波长数值递增。第9、10、13阶IMF基本为无规律波动,平均波长分别为0.319mm、0.649mm和7.39mm。对于无规律IMF,用平均波长评价其整体的波动形式,对于有特征的IMF,观察特征区域的瞬时频率。第11阶IMF在70mm~100mm区域内,瞬时频率曲线较为平坦,说明该段内为固定频率的单一波动形式,该区域瞬时频率均值为0.276[1/mm],波长约为3.62mm。第12阶IMF在35mm~70mm区域内存在一定波动,其对应的瞬时频率数值有一定变化,说明该区域存在一定幅度的变频波动形式,该区域瞬时频率均值为0.232[1/mm],波长约为4.31mm。
4)高、中、低频误差的分解。
美国劳伦斯·利弗莫尔实验室在研制NIF过程中,将空间波长l划分为:
高频带:l<0.12mm,对应为f=l/l>8.33mm-1;
中频带:0.12mm≤l≤33mm,对应为0.03mm-1≤f<8.33mm-1;
低频带:l>33mm,对应为f<0.03mm-1。
基于这一准则,首先根据表1中列出的各阶IMF的平均波长将IMF初步划分为高频组(第1~8阶)、中频组(第9~15阶)和低频组(第16阶),然后根据各阶IMF的特点对分组进行调整,由于9、10阶仍为无规律波动,因此,将其归为高频组。对高频组(第1~10阶)内的IMF叠加得到光学表面高频误差;第9~15阶IMF叠加得到中频误差;低频组内叠加之后与剩余信号rn和拟合面形c0合成得到包含低频误差的面形曲线,图11~13为合成之后的光学元件表面轮廓低、中、高频误差分离结果,其中(a)、(b)、(c)分别表示包含低频误差的面形、中频误差和高频误差。
本发明将光学元件表面的各频段误差进行分离,可更加准确的评价高端光学元件表面质量。相比于其它频段误差分布特性分析方法,本发明根据各阶IMF能够准确判断中频误差分布方位及频率,指导补偿加工和识别误差来源;对光学表面轮廓的高、中、低频误差的分离采取基于空间信号自身局部特征的自适应滤波方式,排除普通滤波中频段选择对结果的影响。
表1
本发明实施例中的检测曲线拟合残差的各阶IMF平均波长参见表1。
Claims (1)
1.抛光工件表面轮廓各频段误差的分离方法,其特征在于包括以下步骤:
1)采用光学检测装置得到被测光学元件表面的检测曲线,对检测曲线进行拟合,得到拟合曲线和残差;
2)将残差进行经验模态分解得到各阶固有模态函数和剩余信号;
3)计算各阶固有模态函数的瞬时频率曲线及其平均波长;
4)根据各阶固有模态函数的平均波长及其特征将所有固有模态函数分为高、中、低频组,将高频和中频组进行组内叠加得到光学表面轮廓的高频和中频误差;将低频组、剩余信号及拟合曲线叠加得到包含低频误差的面形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012103315733A CN102853780A (zh) | 2012-09-07 | 2012-09-07 | 抛光工件表面轮廓各频段误差的分离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012103315733A CN102853780A (zh) | 2012-09-07 | 2012-09-07 | 抛光工件表面轮廓各频段误差的分离方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102853780A true CN102853780A (zh) | 2013-01-02 |
Family
ID=47400585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012103315733A Pending CN102853780A (zh) | 2012-09-07 | 2012-09-07 | 抛光工件表面轮廓各频段误差的分离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102853780A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104034299A (zh) * | 2014-05-27 | 2014-09-10 | 杭州电子科技大学 | 基于经验模态分解的圆度误差评定方法 |
CN104279983A (zh) * | 2014-11-05 | 2015-01-14 | 哈尔滨工业大学 | 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 |
CN105823432A (zh) * | 2016-04-26 | 2016-08-03 | 北京理工大学 | 非连续加工表面系统误差和随机误差分离方法 |
CN109238198A (zh) * | 2018-08-28 | 2019-01-18 | 广西科技大学 | 获取批量加工的曲面零件的最大加工误差的方法 |
CN110333285A (zh) * | 2019-07-04 | 2019-10-15 | 大连海洋大学 | 基于变分模态分解的超声兰姆波缺陷信号识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101240998A (zh) * | 2008-03-14 | 2008-08-13 | 中国人民解放军国防科学技术大学 | 确定性光学加工条件下频段误差分布特性的分析方法 |
CN101986097A (zh) * | 2010-07-09 | 2011-03-16 | 浙江大学 | 在球面面形干涉检测中高精度消除离焦误差及倾斜误差的方法 |
JP2012132682A (ja) * | 2010-12-17 | 2012-07-12 | Canon Inc | 計測方法及び計測装置 |
-
2012
- 2012-09-07 CN CN2012103315733A patent/CN102853780A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101240998A (zh) * | 2008-03-14 | 2008-08-13 | 中国人民解放军国防科学技术大学 | 确定性光学加工条件下频段误差分布特性的分析方法 |
CN101986097A (zh) * | 2010-07-09 | 2011-03-16 | 浙江大学 | 在球面面形干涉检测中高精度消除离焦误差及倾斜误差的方法 |
JP2012132682A (ja) * | 2010-12-17 | 2012-07-12 | Canon Inc | 計測方法及び計測装置 |
Non-Patent Citations (2)
Title |
---|
D. M. AIKENS ET AL.: "The use of Power Spectral Density (PSD) functions in specifying optics for the National Ignition Facility", 《PROCEEDING OF SPIE》 * |
NORDEN E HUANG ET AL.: "The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis", 《PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104034299A (zh) * | 2014-05-27 | 2014-09-10 | 杭州电子科技大学 | 基于经验模态分解的圆度误差评定方法 |
CN104034299B (zh) * | 2014-05-27 | 2016-08-17 | 杭州电子科技大学 | 基于经验模态分解的圆度误差评定方法 |
CN104279983A (zh) * | 2014-11-05 | 2015-01-14 | 哈尔滨工业大学 | 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 |
CN104279983B (zh) * | 2014-11-05 | 2016-10-05 | 哈尔滨工业大学 | 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 |
CN105823432A (zh) * | 2016-04-26 | 2016-08-03 | 北京理工大学 | 非连续加工表面系统误差和随机误差分离方法 |
CN105823432B (zh) * | 2016-04-26 | 2018-06-19 | 北京理工大学 | 非连续加工表面系统误差和随机误差分离方法 |
CN109238198A (zh) * | 2018-08-28 | 2019-01-18 | 广西科技大学 | 获取批量加工的曲面零件的最大加工误差的方法 |
CN110333285A (zh) * | 2019-07-04 | 2019-10-15 | 大连海洋大学 | 基于变分模态分解的超声兰姆波缺陷信号识别方法 |
CN110333285B (zh) * | 2019-07-04 | 2021-07-27 | 大连海洋大学 | 基于变分模态分解的超声兰姆波缺陷信号识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102853780A (zh) | 抛光工件表面轮廓各频段误差的分离方法 | |
US6717655B2 (en) | Micropulse lidar system | |
CN103198322A (zh) | 基于机器视觉的磁瓦表面缺陷特征提取及缺陷分类方法 | |
CN102837228B (zh) | 非球面光学零件短流程加工方法 | |
CN108804814B (zh) | 一种提取路面谱参数的方法 | |
CN105113375A (zh) | 一种基于线结构光的路面裂缝检测系统及其检测方法 | |
CN101995437B (zh) | 一种拉深件裂纹声发射信号的特征提取方法 | |
CN106447086A (zh) | 一种基于风电场数据预处理的风电功率组合预测方法 | |
CN103969162A (zh) | 一种基于数据融合的矿井煤尘浓度测量的方法及装置 | |
CN101240998A (zh) | 确定性光学加工条件下频段误差分布特性的分析方法 | |
CN103464906A (zh) | 激光焊接质量在线检测方法 | |
CN107607195B (zh) | 一种基于复振幅实时获取的光束质量测量方法 | |
CN111809463B (zh) | 一种基于ai方法的钢轨智能打磨系统及相应的打磨方法 | |
EP0837453A3 (en) | Speech analysis method and speech encoding method and apparatus | |
CN101524784A (zh) | 一种基于聚氨酯橡胶模的激光冲击成形方法和装置 | |
CN103839106A (zh) | 一种基于遗传算法优化bp神经网络的球磨机负荷检测方法 | |
CN107238500A (zh) | 汽车操纵稳定性试验快速评价系统建立方法 | |
CN103868884A (zh) | 一种基于调制系数的气体吸收率在线测量方法 | |
CN103440657A (zh) | 一种在线道路裂缝筛查方法 | |
CN102982339A (zh) | 一种高光谱特征变量选取的方法 | |
CN105891787B (zh) | 基于最小二乘逼近的一阶海杂波检测方法 | |
CN105975669A (zh) | 基于cae碰撞仿真的汽车零部件损伤评定方法及装置 | |
CN108032145B (zh) | 射频离子源工作参数优化方法 | |
CN103203670B (zh) | 基于最大信息熵与方向散度的磨削颤振预测方法 | |
CN102519577A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130102 |