CN106826400B - 一种复杂曲面组合加工方法 - Google Patents
一种复杂曲面组合加工方法 Download PDFInfo
- Publication number
- CN106826400B CN106826400B CN201610590629.5A CN201610590629A CN106826400B CN 106826400 B CN106826400 B CN 106826400B CN 201610590629 A CN201610590629 A CN 201610590629A CN 106826400 B CN106826400 B CN 106826400B
- Authority
- CN
- China
- Prior art keywords
- face shape
- shape error
- low order
- zernike
- residence time
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B1/00—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B13/00—Machines or devices designed for grinding or polishing optical surfaces on lenses or surfaces of similar shape on other work; Accessories therefor
- B24B13/0031—Machines having several working posts; Feeding and manipulating devices
- B24B13/0037—Machines having several working posts; Feeding and manipulating devices the lenses being worked by different tools, e.g. for rough-grinding, fine-grinding, polishing
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Grinding And Polishing Of Tertiary Curved Surfaces And Surfaces With Complex Shapes (AREA)
Abstract
本发明公开一种复杂曲面组合加工方法,一、将待加工复杂曲面面形误差zini表示为n*3的矩阵形式[xi,yi,zi];二、通过zernike多项式拟合得到各数据点的低阶面形误差zzernike_i,将面形误差zini去除低阶面形误差zzernike_i,剩下中高频误差;三、采用反卷积的方法求解利用应力盘或大磨头去除低阶面形误差的加工驻留时间T1;四、根据加工驻留时间T1与应力盘或大磨头去除函数矩阵计算得到理论材料去除zremoval1,并得到去除低阶面形误差后残留的中高频误差;五、根据中高频误差zmid与小磨头去除函数计算小磨头驻留时间T2,通过小磨头去除函数与其驻留时间T2的计算得到小磨头材料去除zremoval2,分别经过低阶面形去除,中高频材料去除后加工完成,本发明能够实现复杂曲面高效高精度的加工。
Description
技术领域
本发明属于光学系统加工的技术领域,具体涉及一种复杂曲面组合加工方法。
背景技术
非球面光学元件可以有效降低光学系统复杂程度,提升光学系统性能,因此在现代光学系统中得到了广泛的应用。近年来,随着现代光学技术的不断发展,对光学系统的性能要求越来越高,尤其在空间遥感与深空探测方面,对光学系统的分辨率等性能参数提出了更高的要求。
为满足现代光学系统的性能要求,光学系统中核心器件—主反射镜的口径不断增大,由数百毫米增至数米甚至达到数十米。这就给大口径复杂曲面的加工带来了更大的挑战。
复杂曲面的加工主要面临两方面困难:
首先是加工效率问题,光学曲面面积与其半径是平方关系,随着光学曲面半径的增大,需要加工的光学表面积迅速增大,因此需要进一步提升大口径复杂曲面的加工效率。
其次,复杂曲面的加工是一个复杂,漫长的过程,整个加工过程中任何微小的失误都可能带来严重的后果,因此,需要对每一次的加工过程进行准确的仿真,避免实际加工结果与理论计算结果不符。
发明内容
有鉴于此,本发明提供了一种复杂曲面组合加工方法,能够实现复杂曲面高效高精度的加工。
实现本发明的技术方案如下:
一种复杂曲面组合加工方法,包括以下步骤:
步骤一、将待加工复杂曲面面形误差zini表示为n*3的矩阵形式[xi,yi,zi],其中,xi和yi为第i个数据点的坐标,zi为第i个数据点矢高,即第i个点处面形误差,i=1,2,……,n;
步骤二、通过zernike多项式拟合得到各数据点的低阶面形误差zzernike_i,将面形误差zini去除低阶面形误差zzernike_i,剩下中高频误差;
步骤三、采用反卷积的方法求解利用应力盘或大磨头去除低阶面形误差的加工驻留时间T1;
步骤四、根据加工驻留时间T1与应力盘或大磨头去除函数矩阵计算得到理论材料去除zremoval1,并根据式(3)得到去除低阶面形误差后残留的中高频误差;
zmid=zini-zremoval1 (3)
步骤五、根据中高频误差zmid与小磨头去除函数计算小磨头驻留时间T2,通过小磨头去除函数与其驻留时间T2的计算得到小磨头材料去除zremoval2,分别经过低阶面形去除,中高频材料去除后可得到大小磨头组合加工结果。
进一步地,步骤二中获得低阶面形误差zzernike_i具体为:
步骤2.1、将各数据点坐标由直角坐标系(xi,yi)转变为极坐标系(ρi,θi),并将半径ρ归一化;
步骤2.2、使用zernike多项式拟合复杂曲面面形误差zini,至少拟合zernike前9项,各数据点极坐标(ρi,θi)经zernike拟合获得各数据点低阶面形误差zzernike_i。
进一步地,步骤三具体为:
步骤3.1、应力盘或大磨头在工件表面的材料去除可以用应力盘或大磨头的去除函数沿驻留点卷积方程表示为:
E(x,y)=R(x,y)**D(x,y) (1)
其中,E(x,y)为材料去除,R(x,y)为去除函数,D(x,y)驻留时间分布,**表示卷积符号;
步骤3.2、将卷积方程转换为矩阵方程:
[ei]=[rij][tj] (2)
其中,ei表示第i个数据点的面形误差,rij表示应力盘或大磨头在第j个驻留点处对第i个数据点单位时间内的材料去除,tj表示应力盘或大磨头在第j个驻留点处的驻留时间,其中,j=1,2,3…,m;
步骤3.3、将所述低阶面形误差zzernike_i代入[ei],通过去除函数实验获得去除函数的材料去除矩阵[rij],求解[tj],[tj]即为所述加工驻留时间T1。
有益效果:
本发明是将复杂曲面的面形误差离散化为矩阵形式后分离为低阶面形误差与中高频面形误差。并且采用多种加工方式组合加工:使用应力盘或大磨头去除低阶面形误差,用小磨头、磁流变与离子束去除中高频面形误差。使用矩阵反卷积算法求解驻留时间,对加工过程进行理论计算,以获得高效率高精度的加工策略指导加工。
附图说明
图1为本发明方法流程示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
如图1所示,本发明提供了一种复杂曲面组合加工方法,包括以下步骤:
步骤一、将待加工复杂曲面面形误差zini表示为n*3的矩阵形式[xi,yi,zi],其中,xi和yi为第i个数据点的坐标,zi为第i个数据点矢高,即第i个点处面形误差,i=1,2,……,n;
步骤二、将面形误差分离为中高频面形误差和低阶面形误差:通过zernike多项式拟合得到各数据点的低阶面形误差zzernike_i;将面形误差zini去除低阶面形误差zzernike_i,剩下中高频误差;步骤二中获得低阶面形误差zzernike_i具体为:
步骤2.1、将各数据点坐标由直角坐标系(xi,yi)转变为极坐标系(ρi,θi),并将半径ρ归一化;
步骤2.2、使用zernike多项式拟合复杂曲面面形误差zini,至少拟合zernike前9项,根据需要选择拟合精度,最多可拟合至zernike前35项,各数据点极坐标(ρi,θi)经zernike拟合获得各数据点低阶面形误差zzernike_i。
步骤三、采用反卷积的方法求解利用应力盘或大磨头去除低阶面形误差的加工驻留时间T1;步骤三具体为:
步骤3.1、应力盘或大磨头在工件表面的材料去除可以用应力盘或大磨头的去除函数沿驻留点卷积表示:
E(x,y)=R(x,y)**D(x,y) (1)
其中,E(x,y)为材料去除,R(x,y)为去除函数,D(x,y)驻留时间分布,**表示卷积符号;
步骤3.2、将卷积方程转换为矩阵方程:
[ei]=[rij][tj] (2)
其中,ei表示第i个数据点的面形误差,rij表示应力盘或大磨头在第j个驻留点处对第i个数据点单位时间内的材料去除,tj表示应力盘或大磨头在第j个驻留点处的驻留时间,其中,j=1,2,3…,m;即应力盘或大磨头有m个驻留点。
式(2)也可写为:
步骤3.3、将所述低阶面形误差zzernike_i代入[ei],通过去除函数实验(或根据磨头尺寸、运动方式等参数的计算)获得去除函数的材料去除矩阵[rij],求解[tj],[tj]即为所述加工驻留时间T1。因为驻留时间非负,因此可以通过正则化方法或非负最小二乘法求得大磨头最优解T_big为大磨头去除低阶面形误差的驻留时间分布。
步骤四、根据加工驻留时间T1与应力盘或大磨头去除函数矩阵计算得到理论材料去除Zremoval1,并根据式(3)得到去除低阶面形误差后残留的中高频误差;
zmid=zini-zremoval1 (3)
步骤五、根据中高频误差zmid与小磨头去除函数计算小磨头驻留时间T2,通过小磨头去除函数与其驻留时间T2的计算得到小磨头材料去除zremoval2,分别经过低阶面形去除,中高频材料去除后可得到大小磨头组合加工结果。
为实现上述流程,编写了仿真计算软件,其可实现误差分离、单次加工方法驻留时间求解仿真、组合加工方法驻留时间求解仿真。根据过去加工实验建立了应力盘下、小磨头、磁流变加工方法的去除函数数据库,根据实际需要设定去除函数参数并进行仿真计算。
为验证该算法及仿真模型,使用实际面形进行仿真计算,求解驻留时间后,得到虚拟加工结果,加工前面形误差RMS为0.127λ,仿真加工结果面形误差为0.025λ。收敛效率达到80%,验证了组合加工算法的有效性。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种复杂曲面组合加工方法,其特征在于,包括以下步骤:
步骤一、将待加工复杂曲面面形误差zini表示为n*3的矩阵形式[xi,yi,zi],其中,xi和yi为第i个数据点的坐标,zi为第i个数据点矢高,即第i个点处面形误差,i=1,2,……,n;
步骤二、通过zernike多项式拟合得到各数据点的低阶面形误差zzernike_i,将面形误差zini去除低阶面形误差zzernike_i,剩下中高频误差;
步骤三、采用反卷积的方法求解利用应力盘或大磨头去除低阶面形误差的加工驻留时间T1;
步骤四、根据加工驻留时间T1与应力盘或大磨头去除函数矩阵计算得到理论材料去除zremoval1,并根据式(3)得到去除低阶面形误差后残留的中高频误差;
zmid=zini-zremoval1 (3)
步骤五、根据中高频误差zmid与小磨头去除函数计算小磨头驻留时间T2,通过小磨头去除函数与其驻留时间T2的计算得到小磨头材料去除zremoval2,分别经过低阶面形去除,中高频材料去除后可得到大小磨头组合加工结果。
2.如权利要求1所述的一种复杂曲面组合加工方法,其特征在于,
步骤二中获得低阶面形误差zzernike_i具体为:
步骤2.1、将各数据点坐标由直角坐标系(xi,yi)转变为极坐标系(ρi,θi),并将半径ρ归一化;
步骤2.2、使用zernike多项式拟合复杂曲面面形误差zini,至少拟合zernike前9项,各数据点极坐标(ρi,θi)经zernike拟合获得各数据点低阶面形误差zzernike_i。
3.如权利要求1所述的一种复杂曲面组合加工方法,其特征在于,步骤三具体为:
步骤3.1、应力盘或大磨头在工件表面的材料去除可以用应力盘或大磨头的去除函数沿驻留点卷积方程表示为:
E(x,y)=R(x,y)**D(x,y) (1)
其中,E(x,y)为材料去除,R(x,y)为去除函数,D(x,y)驻留时间分布,**为卷积符号;
步骤3.2、将卷积方程转换为矩阵方程:
[ei]=[rij][tj] (2)
其中,ei表示第i个数据点的面形误差,rij表示应力盘或大磨头在第j个驻留点处对第i个数据点单位时间内的材料去除,tj表示应力盘或大磨头在第j个驻留点处的驻留时间,其中,j=1,2,3…,m;
步骤3.3、将所述低阶面形误差zzernike_i代入[ei],通过去除函数实验获得去除函数的材料去除矩阵[rij],求解[tj],[tj]即为所述加工驻留时间T1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590629.5A CN106826400B (zh) | 2016-07-25 | 2016-07-25 | 一种复杂曲面组合加工方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590629.5A CN106826400B (zh) | 2016-07-25 | 2016-07-25 | 一种复杂曲面组合加工方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106826400A CN106826400A (zh) | 2017-06-13 |
CN106826400B true CN106826400B (zh) | 2018-05-01 |
Family
ID=59145852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610590629.5A Active CN106826400B (zh) | 2016-07-25 | 2016-07-25 | 一种复杂曲面组合加工方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106826400B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107520683B (zh) * | 2017-08-31 | 2019-03-19 | 中国工程物理研究院激光聚变研究中心 | 子口径抛光的等高线路径规划方法 |
CN107918178A (zh) * | 2017-12-14 | 2018-04-17 | 中国航空工业集团公司洛阳电光设备研究所 | 一种提高100mm~150mm口径非球面胶合精度的方法 |
CN108555729B (zh) * | 2017-12-20 | 2020-01-10 | 复旦大学 | 光学镜片小磨头光学加工中的边缘误差控制方法 |
PL240456B1 (pl) * | 2017-12-21 | 2022-04-04 | Politechnika Rzeszowska Im Ignacego Lukasiewicza | Sposób nadzorowania procesu szlifowania wgłębnego |
CN109227226B (zh) * | 2018-11-12 | 2020-11-03 | 中国科学院光电技术研究所 | 一种光学元件加工过程中驻留时间的匀滑方法 |
CN109437599A (zh) * | 2018-12-11 | 2019-03-08 | 北京控制工程研究所 | 一种交会对接敏感器用超高精度反射镜一体式加工方法 |
CN110134915B (zh) * | 2019-05-16 | 2022-02-18 | 中国工程物理研究院激光聚变研究中心 | 一种磁流变抛光驻留时间的处理方法及装置 |
CN113275976B (zh) * | 2020-02-20 | 2022-07-05 | 中国科学院长春光学精密机械与物理研究所 | 计算机控制光学面形误差收敛加工方法、装置及介质 |
CN111843754B (zh) * | 2020-07-31 | 2021-10-26 | 中国人民解放军国防科技大学 | 轴类零件外圆圆度确定性修形方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7234810B2 (en) * | 2003-11-14 | 2007-06-26 | Ophthonix, Inc. | System for manufacturing an optical lens |
CN101898324B (zh) * | 2010-07-28 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 高陡度镜面的离子束抛光方法 |
CN102848287B (zh) * | 2012-09-14 | 2014-09-17 | 中国人民解放军国防科学技术大学 | 光学零件中高频误差的组合式加工方法 |
CN103921176B (zh) * | 2014-03-27 | 2017-06-09 | 中国科学院长春光学精密机械与物理研究所 | 适用于超大口径光学加工的磁流变抛光装置 |
CN105643395A (zh) * | 2016-01-04 | 2016-06-08 | 中国科学院长春光学精密机械与物理研究所 | 一种光学自由曲面磨削成型方法 |
-
2016
- 2016-07-25 CN CN201610590629.5A patent/CN106826400B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106826400A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106826400B (zh) | 一种复杂曲面组合加工方法 | |
Krefl et al. | Machine learning of Calabi-Yau volumes | |
CN103198322A (zh) | 基于机器视觉的磁瓦表面缺陷特征提取及缺陷分类方法 | |
CN103886569A (zh) | 连续帧多特征点无人机侦察图像并行与匹配精度约束的拼接方法 | |
CN103400399B (zh) | 一种基于空间矩的线结构光中心提取方法 | |
CN103294899A (zh) | 一种计算小型实验反应堆堆芯中子通量分布的方法 | |
CN109909815A (zh) | 光学复杂曲面元件的磁流变抛光补偿加工方法、系统及介质 | |
Yan et al. | Generation of consistent skin model shape based on FEA method | |
CN103353907A (zh) | 一种法兰连接螺栓极限强度校核的计算方法 | |
CN106294936B (zh) | 一种化学机械研磨模拟方法 | |
CN103927731A (zh) | 无需pos辅助的低空遥感影像快速自动拼接方法 | |
CN103247041A (zh) | 一种基于局部采样的多几何特征点云数据的分割方法 | |
CN101644915A (zh) | 基于去除函数预测模型的计算机控制抛光方法 | |
CN104616349A (zh) | 基于局部曲面变化因子的散乱点云数据精简处理方法 | |
CN103778433B (zh) | 基于点到直线距离的广义点集匹配方法 | |
CN101504770B (zh) | 一种结构光光条中心的提取方法 | |
CN104392130B (zh) | 多轴疲劳的最损伤载荷方向的确定方法及其应用 | |
CN104318022A (zh) | 一种预测工件表面粗糙度并提高切削效率的方法 | |
CN106855895B (zh) | 一种基于离子束技术的复杂曲面去除函数计算方法 | |
Yu et al. | Method for discriminating geometric feasibility in assembly planning based on extended and turning interference matrix | |
Yang et al. | Aerodynamic shape optimisation of benchmark problems using SU2 | |
CN201514207U (zh) | 一种复杂曲面的数字化检测系统 | |
CN102411652B (zh) | 带耳片矩形梁截面尺寸确定方法 | |
CN106250458B (zh) | 一种单位格点几何信息提取方法 | |
CN102542550A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |