Research on high deterministic removal characteristics of fused quartz using RF focused ion beam source
-
摘要: 在离子束抛光工艺过程中,材料确定性去除特性对预测光学元件的各工位材料去除量和驻留时间具有极其重要的作用。采用射频离子源对熔石英光学元件的离子束刻蚀特性进行了研究,利用ZYGO激光干涉仪获得准确的去除函数,系统分析了气体流量、屏栅电压、离子束入射角和工作距离等因素对熔石英去除函数的影响,并分析了各单一工艺因素微小扰动时,材料峰值去除率、半高宽和体积去除率的相对变化率。实验结果表明,相同工作真空条件下,工作气体质量流量的微小变化对去除函数影响极小,在典型的工艺条件下,屏栅电压在±5 V、离子束入射角±1°、工作距离在±0.5 mm范围内变化时,熔石英峰值去除率、体积去除率和峰值半高宽的相对变化均小于5%,去除函数具有较好的确定性和稳定性。Abstract: The characterization of the deterministic material removal plays an important role in the ion-beam figuring (IBF) polishing process, which is used to predict the material removal and help to determine the dwell time needed at different positions of the workpiece. The ion beam etching characteristics of fused quartz optical components were studied by radio frequency(RF) ion source, the precise removal function was calculated and obtained by using ZYGO interferometer, and the influences of the plasma source operational parameters such as work gas mass flow rate, screen voltage, incidence angle of ion beam and work distance on the material removal function of fused quarts were investigated. Furthermore, the relative change rates of the material peak removal rate, the full-width-at-half-maximum (FWHM) of peak and the material volume removal rate were measured and calculated when individual work parameter was perturbed. The experiment results show that the variation of removal functions varies slightly with the gas mass flow rate under the same work vacuum pressure. It is found that the relative errors of the removal function's peak removal efficiency, the FWHM and the volume removal efficiency can all controlled in 5% when the screen voltage fluctuation is ±5 V, the incidence angle deviation of ion beam is ±1° and the work distance deviation is ±0.5 mm, respectively. Hence the removal function of IBF for fused quartz by using RF focusing ion source has better certainty and stability.
-
Keywords:
- radio frequency ion source /
- removal function /
- ion beam figuring /
- fused quartz
-
引言
飞行目标的偏航角、俯仰角等姿态特征参数是火箭发射装备、发射火药及国防工程设计、定型、验收和故障诊断等过程的主要参数,准确地测量这些参数对提高武器效能和国防工程的防御能力意义重大[1-2]。目前,靶场姿态测量以两站以上交会测量为主,但在设备故障及捕获目标失误的情况下,单站测姿必须发挥重要作用,故研究单站测姿非常有必要。
综合现有文献的单站测姿方法,基本可划分为需要建立精确目标模型前提下的单站模型匹配法; 不需要相机内参数的DLT方法[3],但该方法至少需要5对以上的物像点对应; 仅适用于投影轴对称目标的类似比例方法包括长宽比、椭圆长短轴比例法[4]等,该方法将透视投影做了平行投影的近似; 为突破该方法仅适合投影轴对称目标的局限,文献[5]利用了主轴以外的侧翼等信息,提出利用弱透视投影理论求解飞机等非投影轴对称目标的姿态角,但对模型缺乏系统、完整性分析; 文献[6]中使用的像点坐标差值比例法,该方法保持了透视模型的准确度,且利用了经纬仪较为准确的测角信息,但至少需要4个特征点,并对特征点的选取有一定要求,另外用到了站点的位置信息; 需要利用相机内外参数的PNP[7-9]以及光雷一体化方法[10],其中PNP也需要3对以上的控制点对应,且大多基于迭代的方法。
本文测姿方法立足靶场实际,不需要建立目标模型,在只有单站图像序列,仅仅依靠某些先验特征长度、摄像系统相机内外参数以及成像信息的情况下,不需要利用测站位置信息。在保持透视投影模型准确度的前提下,对于投影轴对称目标,可以精确定位弹首、弹尾2个特征点作为输入条件,较准确地获取目标的姿态角。
1 基本原理
本文以中轴线为透视投影对象,将已知长度的目标向已知内、外相机参数的姿态测量系统像面投影,获得模拟投影像长; 同时提取实际图像序列中的目标实际像长,两者匹配相等,即可获取包括参考物距、目标深度等关键中间参数,然后利用关键参数,可计算出目标姿态角。大致流程图如图 1所示。
2 算法基础—透视投影2种等效方法
透视投影2种等效方式是本文获取关键中间参数的重要环节,可等效为2种方式。
2.1 目标跨度内的物面等效为垂直于光轴的若干连续物面组
经过目标垂直于主光轴的平面作为近似目标平面组,这样的平面连续且有无限个,故本来连续的目标物可以离散化为无数个连续的近似目标平面,透视投影则为无数个离散目标微元在像面上投影的累加。
2.2 以任意的其中一个等效物面为基准
可以近似目标平面组的任意一个平面为等效目标平面,按照透视投影的关系,连接光心及目标点的直线与等效目标平面相交,交点则为该目标点在等效目标平面的等效物点,如此获得等效物面内的等效物长。
3 目标姿态求解
3.1 方法分析及关键过程设置
本文方法涉及的主要坐标系,除了基准坐标系、摄像机坐标系外,还涉及目标体坐标系,与测量坐标系平行的过渡坐标系(零姿态坐标系),O0-XpYpZp所有坐标系均符合右手坐标系。其中目标体坐标系O0-X0Y0Z0以体现目标姿态状态的中轴线为主轴建立; 设置初始姿态下的体坐标系与摄像机坐标系平行。
3.2 算法推导
已知条件:
物方 LAB=L;
像方 AB在像面的成像情况:可利用的包括A点成像AA(x1, y1),B点成像BB(x2, y2);
相机内参数 像主点(x0, y0)以及焦距f已知,畸变暂不考虑;
相机外参数 光轴指向α、β角,即为摄像机坐标系相对于基准坐标系的旋转角度。
求解:目标中轴线相对于基准坐标系的姿态ψ和φ。
透视投影如图 2所示,始终设B所在物面为参考物面,A向B所在的参考物面投影点为D,∠ABD=θ,由上面的直线透视投影原理,计算过程为
1) 计算中间变量u及sinθ
$$ \left\{ \begin{align} & \sqrt{{{({{x}_{1}}-{{x}_{2}})}^{2}}+{{({{y}_{1}}-{{y}_{2}})}^{2}}}\times \mu \times {{10}^{-6}}=1 \\ & \frac{u}{{{10}^{3}}f}=\frac{\sqrt{{{(\frac{u}{\text{cos}{{\alpha }_{1}}})}^{2}}+{{(\frac{u}{\text{cos}{{\alpha }_{2}}})}^{2}}-\frac{u}{u+L~\text{sin}\theta }\left[ \frac{{{(u+L\text{sin}\theta )}^{2}}}{{{(\text{cos}{{\alpha }_{1}})}^{2}}}+\frac{{{u}^{2}}}{(\text{cos}{{\alpha }_{2}})}-{{L}^{2}} \right]}}{\mu \sqrt{{{({{x}_{1}}-{{x}_{2}})}^{2}}+{{({{y}_{1}}-{{y}_{2}})}^{2}}}} \\ \end{align} \right. $$ (1) 2) 计算关键特征点在摄像机坐标系下的坐标;
3) 计算目标相对于摄像机坐标系的姿态角;
4) 计算目标中轴线矢量在过渡坐标系下的坐标,即计算过渡坐标系下的目标中轴线矢量;
5) 经进一步坐标系旋转,计算单位中轴线矢量在基准坐标系下的坐标(X, Y,Z), 姿态角为
$$ \left\{ \begin{array}{l} \psi = {\tan ^{ - 1}}\frac{{{Z^\prime }}}{{{X^\prime }}} + k\pi \\ \phi = {\tan ^{ - 1}}\frac{{{Y^\prime }}}{{\sqrt {{X^\prime }^2 + {Z^\prime }^2} }} \end{array} \right.(k以实际象限取值) $$ (2) 4 试验验证与误差分析
4.1 试验验证
为了验证姿态模型的精度,我们在靶场距离姿态测量系统测试点1.7 km左右处放置了一门120火炮,以炮管为目标进行测量和验证,炮管的的方位和俯仰角度值由全站仪测得,姿态测量系统对炮管进行静态或动态扫描拍摄,获得了系列不同角度值的测量图像,以其中2组为例,如图 3所示。
其中姿态测量系统内外参数为:拍摄频率为100 fps,可见光高速电视分辨率1 536×1024,焦距值为1 523.1 mm,光轴指向如表 1所示,被研究的关键段长度约为2.655 m,用本文的算法进行姿态处理结果如表 1所示。
表 1 算法结果单位(°) Table 1. Results of algorithmα β 偏航绝对误差 俯仰角误差 1 323.063 3 -0.913 1 -1.975 6 1.5 2 323.053 2 -0.940 2 1.575 6 -0.8 3 323.102 1 -0.878 9 1.224 4 0.8 4 323.047 3 -0.861 6 -1.075 6 0.6 5 323.018 3 -0.870 3 1.775 6 0.3 6 323.019 9 -0.824 8 1.424 4 0.8 7 323.029 -0.867 1.924 4 1.2 8 323.052 -0.844 2.024 4 1.1 RMSE 1.7 1.0 从表 1中可得,姿态角度误差值偏航角约为1.7°,俯仰角误差约1.0°,主要原因在于偏航角主要靠对目标深度的鉴别力,而实际物距相对目标的比值大,故系统对偏航角的分辨力差。下面按照输入参数对模型进行误差分析与仿真。
4.2 误差分析
从模型(2)分析
$$ \left\{ \begin{align} & \psi =f(l, {{x}_{1}}, {{y}_{1}}, {{x}_{2}}, {{y}_{2}}, f, \alpha , \beta ) \\ & \phi =g(L, {{x}_{1}}, {{y}_{1}}, {{x}_{2}}, {{y}_{2}}, f, \alpha , \beta ) \\ \end{align} \right. $$ (3) 下面逐一分析关键模型参数对姿态角的影响。
1) 分析经纬仪指向误差对姿态角误差的影响。以表 1中经纬仪的测角数据为基准,按照目前设备的指向精度为指标参考,左右放缩0.3°为极限标准,分别以光轴指向的方位角和俯仰角2种情况进行分析,经纬仪指向误差对姿态角的测量误差趋势如图 4所示。
从图 4可知,测角精度对整个姿态结果有一定的影响,以目前的指向精度指标不超过0.3°为依据,对姿态角的影响不超过0.3°,光轴指向的方位角测角误差对偏航角影响大于俯仰角; 光轴指向的俯仰角测角误差对俯仰角的影响则大于偏航角。
2) 分析作用距离对姿态角误差的影响。从模型上分析,对作用距离敏感性较强的是偏航角,故作用距离对姿态角的影响主要以偏航角展开,作用距离对姿态角误差的影响如图 5所示。
从图 5中可以看出,作用距离对姿态角误差影响比较大,当横轴比值为4 800时,则姿态角绝对误差可达1°,当比值为7 200时,绝对误差可达2°。
3) 分析像点判读误差对姿态角的影响。以文中投影轴对称目标试验为基础,物距大约为1 700 m左右,在此物距及其它系统参数条件下,以表 1中的某一组数据为基础,2个端点像点反向移动下图中横轴的像素数,像素判读误差对姿态角误差影响如图 6所示。
从图 6可知,像素判读误差对姿态角误差影响占绝对比重,故在测量中,尽量利用目前比较成熟的亚像素定位精度技术,提高特征像点定位精度。
5 结论
文章提出了直线透视投影的2种形式,以长度已知的特征直线为研究对象,通过直线透视投影,可直接获取像长,或在某等效物面透视投影得到等效物长。以此为基础,通过像机内参数和外参数建立了投影轴对称的姿态测量模型; 并以炮管试验为依托,获得序列图像和相对真值,对透视算法进行了数据验证,结果为在试验设置的条件下投影轴对称姿态精度偏航角约为1.7°,俯仰角精度约为1.0°; 最后对算法进行了误差分析及实验仿真,分析了关键因素对姿态角误差的影响。
对于靶场姿态测量,文章的算法简单、方便,不需要3D建模及测试站点信息,测试精度满足靶场单站测姿需求,且该算法同样适用于非投影轴对称目标,为靶场姿态测量提供了一种新的思路和方法。
-
-
[1] ARNOLD T, BÖHM G, FECHNER R, et al. Ultra-precision surface finishing by ion beam and plasma jet techniques: status and outlook[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2010, 616(2/3): 147-156. http://cn.bing.com/academic/profile?id=096d17779b3399df844c1ba271a3a04c&encoded=0&v=paper_preview&mkt=zh-cn
[2] ALLEN L N. Progress in ion figuring large optics[J]. SPIE, 1994, 2428: 237-247. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC026084881
[3] 李圣怡, 戴一凡, 解旭辉, 等.中大口径非球面镜制造与测量新技术[M].北京:国防工业出版社, 2011. LI Shengyi, DAI Yifan, XIE Xuhui, et al. New technology for manufacturing and measurement of large and middle-scale aspheric surfaces[M]. Beijing: National Defense Industry Press, 2011.
[4] FAWCETT S C. Ion beam figuring of small optical components[J]. Optical Engineering, 1995, 34(12): 3565. doi: 10.1117/12.215648
[5] ALLEN L N, HANNON J H, WAMBACH R W. Final surface error correction of an off-axis aspheric petal by ionFiguring[J]. SPIE, 1991, 1543:190-200. https://www.spiedigitallibrary.org/redirect/proceedings/proceeding?doi=10.1117/12.51181
[6] 戴一帆, 周林, 解旭辉, 等.离子束修形技术[J].应用光学, 2011, 32(4):753-760. doi: 10.3969/j.issn.1002-2082.2011.04.032 DAI Yifan, ZHOU Lin, XIE Xuhui, et al. Ion beam figuring technology[J]. Journal of Applied Optics, 2011, 32(4): 753-760. doi: 10.3969/j.issn.1002-2082.2011.04.032
[7] 王东方.大气等离子体加工熔石英技术研究[J].陕西理工学院学报:自然科学版, 2013, 29(2):1-5. http://d.old.wanfangdata.com.cn/Periodical/sxgxyxb201302001 WANG Dongfang. Research on atmospheric pressure plasma processing technology of fused silicon[J]. Journal of Shaanxi University of Technology: Natural Science Edition, 2013, 29(2):1-5. http://d.old.wanfangdata.com.cn/Periodical/sxgxyxb201302001
[8] PHILIP S H.Optica needs for futrue space telescopes[J]. SPIE, 2003, 5180:1-5.
[9] ALLEN D M, SHORE P, EVANS R W, et al. Ion beam, focused ion beam, and plasma discharge machining[J]. CIRP Annals, 2009, 58(2): 647-662. doi: 10.1016/j.cirp.2009.09.007
[10] 谌贵阳, 解旭辉, 周林, 等.光学加工射频离子源的去除函数与稳定性研究[J].航空精密制造技术, 2016, 52(5): 19-22, 40. doi: 10.3969/j.issn.1003-5451.2016.05.005 CHEN Guiyang, XIE Xuhui, ZHOU Lin, et al. Removal function and stability study of RF ion source for optics figuring[J]. Aviation Precision Manufacturing Technology, 2016, 52(5): 19-22, 40. doi: 10.3969/j.issn.1003-5451.2016.05.005
[11] 刘金声.离子束技术及应用[M].北京:国防工业出版社, 1995. LIU Jinsheng.Ion beam technology and application[M]. Beijing:National Defense Industry Press, 1995.
[12] 张振鹏, 解旭辉, 周林, 等.光学加工射频中和器的结构设计和性能研究[J].航空精密制造技术, 2017, 53(3): 8-11, 15. doi: 10.3969/j.issn.1003-5451.2017.03.003 ZHANG Zhenpeng, XIE Xuhui, ZHOU Lin, et al. Structure design and performance research of RF neutralizer for optics figuring[J]. Aviation Precision Manufacturing Technology, 2017, 53(3): 8-11, 15. doi: 10.3969/j.issn.1003-5451.2017.03.003
[13] 唐瓦, 邓伟杰, 郑立功, 等.离子束抛光去除函数计算与抛光实验[J].光学精密工程, 2015, 23(1): 31-39. http://d.old.wanfangdata.com.cn/Periodical/gxjmgc201501005 TANG Wa, DENG Weijie, ZHENG Ligong, et al. Calculation of removal function of ion beam Figuring and polishing experiment[J]. Optics and Precision Engineering, 2015, 23(1): 31-39. http://d.old.wanfangdata.com.cn/Periodical/gxjmgc201501005
[14] ALLEN L N, KEIM R E. An ion figuring system for large optic fabrication[J]. SPIE, 1989, 1168: 33-50. http://cn.bing.com/academic/profile?id=af7c27e0cecf6abbfab242e0f04a4977&encoded=0&v=paper_preview&mkt=zh-cn
[15] 孔兵, 王昭, 谭玉山.激光光斑的高斯拟合[J].激光技术, 2002, 26(4): 277-278. doi: 10.3969/j.issn.1001-3806.2002.04.008 KONG Bing, WANG Zhao, TAN Yushan. Gaussian fitting technique of laser spot[J]. Laser Technology, 2002, 26(4): 277-278. doi: 10.3969/j.issn.1001-3806.2002.04.008
[16] 段沽坪.离子束抛光去除特性研究[D].成都: 中国科学院研究生院(光电技术研究所), 2013. http://cdmd.cnki.com.cn/Article/CDMD-80151-1013026467.htm DUAN Guping.Study of removal characteristics of ion beam polishing[D]. Chengdu: Graduate School of Chinese Academy Sciences(Institute of Optics and electronics), 2013. http://cdmd.cnki.com.cn/Article/CDMD-80151-1013026467.htm
[17] KATARDJIEV I V. A kinematic model of surface evolution during growth and erosion: numerical analysis[J]. Journal of Vacuum Science & Technology A: Vacuum, Surfaces, and Films, 1989, 7(6): 3222-3232. http://cn.bing.com/academic/profile?id=d4b3a5eebf8d3e32cb4590a78cf0f7c9&encoded=0&v=paper_preview&mkt=zh-cn
[18] STEINBR CHEL C. Universal energy dependence of physical and ion-enhanced chemical etch yields at low ion energy[J]. Applied Physics Letters, 1989, 55(19): 1960-1962. doi: 10.1063/1.102336
[19] SIGMUND P. Sputtering by particle bombardment I[M]. Berlin: Springer-Verag, 1981.
[20] 廖文林, 戴一帆, 周林, 等.离子束作用下的光学表面粗糙度演变研究[J].应用光学, 2010, 31(6): 1041-1045. doi: 10.3969/j.issn.1002-2082.2010.06.035 LIAO Wenlin, DAI Yifan, ZHOU Lin, et al. Optical surface roughness in ion beam process[J]. Journal of Applied Optics, 2010, 31(6): 1041-1045. doi: 10.3969/j.issn.1002-2082.2010.06.035
-
期刊类型引用(4)
1. 钟志,王琛,刘磊,孙岩,刘彬,单明广. 数字全息声场成像测量技术综述(特邀). 激光与光电子学进展. 2024(24): 11-24 . 百度学术
2. 席锋,秦华锋,胡莉,林睿. 可视化聚焦纹影系统实现流场中的三维扰动成像. 光学学报. 2021(09): 90-95 . 百度学术
3. 陈水桥,程诗卓,郭鹍,陆子毅,许熠辉. 非均匀介质驻波声场纹影法成像及其声速测量的改善. 实验室研究与探索. 2020(12): 46-50 . 百度学术
4. 罗劲明,薛凌峰,郑景辉,曾锦富,李子文,纪桂彬. 空气折射率梯度分布的可视化实验. 物理实验. 2018(01): 19-21+28 . 百度学术
其他类型引用(2)