留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高精度等晕角测量关键技术研究

强希文 吴敏 宗飞 翟胜伟 胡月宏 封双连 赵军卫 常金勇

李恒, 李洋, 佟建华, 等. 铜簇修饰微电极小型硝酸根电化学检测系统[J]. 强激光与粒子束, 2015, 27: 024154. doi: 10.11884/HPLPB201527.024154
引用本文: 强希文, 吴敏, 宗飞, 等. 高精度等晕角测量关键技术研究[J]. 强激光与粒子束, 2021, 33: 081008. doi: 10.11884/HPLPB202133.210215
Li Heng, Li Yang, Tong Jianhua, et al. Miniaturized electrochemical measurement system for nitrate detection using renewable copper-clusters modified microelectrode[J]. High Power Laser and Particle Beams, 2015, 27: 024154. doi: 10.11884/HPLPB201527.024154
Citation: Qiang Xiwen, Wu Min, Zong Fei, et al. High-precision measurement technique of isoplanatic angle[J]. High Power Laser and Particle Beams, 2021, 33: 081008. doi: 10.11884/HPLPB202133.210215

高精度等晕角测量关键技术研究

doi: 10.11884/HPLPB202133.210215
详细信息
    作者简介:

    强希文(1967—),男,研究员,主要从事大气光学与自适应光学、激光与物质相互作用等方面研究

  • 中图分类号: P412

High-precision measurement technique of isoplanatic angle

  • 摘要: 等晕角是重要大气光学参数之一,对天文光学站址选择、大气湍流效应自适应光学校正效果评价等有着重要的意义。分析了国内外等晕角测量技术,介绍了等晕角测量的基本原理,设计了三环变迹镜并开展了数值验证,分析并提出了高精度等晕角测量仪的主要技术指标,研制了一套高精度等晕角测量仪,结果表明,该仪器能够满足全天时的等晕角测量。
  • 光学晶体加工中工艺参数的选择会直接影响工件表面质量和使用性能。如惯性约束聚变过程中对光学元件可靠性要求极为严格[1]。晶体的高硬脆性及强各向异性等一系列不利因素使加工高质量表面变得困难[2-4]。目前,对晶体材料精密加工而言,单点金刚石超精加工[5-6]相比于磁流变光整加工[7-9]和磨料切割工艺[10]两种常用加工方式,能有效缩短加工周期,表面一致性较好。其过程参数如切削用量、主轴振动、导轨误差、刀具形状、刀具磨损以及切削力的波动将直接影响表面质量。已加工表面微观形貌可看作是加工过程中各参数综合作用的集合,切削用量对微观形貌影响最为明显[11]。本文根据切削用量影响晶体表面微观形貌的频率特征进行分析,获得了切削参数(进给量、切削深度和主轴转速)与实际频率特征关系。

    表面形貌极大程度影响着材料使用功能,而切削用量是影响表面形貌的主要因素。国内外学者已开展机加工过程中切削用量对晶体表面质量及其功能基础性研究。Chen N用自行研制的微型PCD球头铣刀对KDP材料进行加工,获得铣削过程中最优参数并验证了微型PCD球头铣刀对表面微缺陷具有修复能力[12-13]。Zheng W研究SiCp/AI复合材料磨削过程中工艺参数对表面形貌影响,综合分析了进给量、主轴转速、主轴振动及切削深度对表面粗糙度不同影响程度[14]。Wang S对KDP晶体材料进行切削实验,获得进给量对KDP单晶表面质量和切削力的关系并分析了材料韧脆性变化[15]。Tian F发现金刚石刀具车削过程中工件振动对表面形貌有很大影响,建立振动对加工表面影响的理论模型并分析工件和刀具在切削方向上的振动关系[16]。Wang X考虑到加工参数及振动对研磨材料的影响并建立研磨运动模型,发现了磨削表面形貌主要受研磨过程中研磨速度影响[17]。Krolczyk G M对表面形貌参数进行差异化描述,包括表面振幅和粗糙度等参数进行评估[18]。Zhang Q用AISI H13钢球头刀具进行铣削实验,分析了进给量和切削深度对已加工表面形貌的影响关系[19]。Wei W研究了加工参数在高速铣削过程中对铣削力的影响规律[20-21]。但上述研究均未考虑表面形貌是实际频率特征无序叠加的组合体,对微观形貌分析仍停留在表面轮廓纹理分布上。为获得加工过程中切削用量对微观形貌影响规律,本文进行单点金刚石车削KDP单晶实验。PSD结合连续小波方法,进行已加工表面轮廓频率特征的提取,研究切削参数对频率特征的影响。

    功率谱密度(PSD)通常用来分析输入信号的空间频率信息。PSD不仅包含表面垂直方向的高度信息,同时包含重要的横向空间频率分布[22]。表面微观轮廓频率信息由计算一维PSD得到。PSD定义式为

    PSD1(f)=limL|z(f,L)|2L
    (1)

    式中:L是采样长度,f是空间频率,zfL)是表面微观轮廓傅里叶变换的结果。zfL)的计算公式为

    z(f,L)=L/2L/2z(x)exp(j2πfx)dx
    (2)

    式中:zx)是微观轮廓高度信息。

    空间域中,小波变换(WT)定义为分析原始信号fx)的一系列函数,即小波基函数ψx)。ψabx)是小波基ψx)根据尺度因子a和平移因子b构成的一组函数。连续小波变换(CWT)是对空间频率进行多尺度变换[23],能有效地从微观形貌中提取实际频率特征。CWT定义式为

    ψa,b(x)=|a|1/2ψ(xba)
    (3)
    Wψ(f)(a,b)=1af(t)ψ(tba)dt
    (4)

    式中:Wψ是连续小波系数,由ψx)(本文取Mexican小波基)、尺度因子a、平移因子b(本文中b=0)共同决定。

    尺度因子a为轮廓频率特征信息与小波基中心频率的变换关系。频率特征与a对应关系为

    a=fcfsΔ
    (5)

    式中:fc是Mexican小波基中心频率(fc=0.25);Δ是白光干涉仪采样周期;fs是实际频率特征的中心频率。其中fs可表示为

    fs=fLN/2
    (6)

    式中:f是PSD提取表面轮廓的实际频率,N是采样长度上采样点数目。

    式(4)的小波系数Wψ表示原信号与小波基相似程度,与轮廓频率特征无关。需对小波系数逆变换,定义式为

    f(t)=1Cψab[Wψ(f)(a,b)]ψa,b(t)dbdaa2
    (7)

    式中:ft)是小波系数重构的实际频率特征;Cψ是小波容许常数。

    微观形貌处于nm量级时,仍能用描述表面微观轮廓大小的传统统计参数,如轮廓均方根Rq表征轮廓垂直方向的高度信息,可表示为

    Rq=L0[z(x)]2dxL
    (8)

    式中:zx)为轮廓高度。

    本文中KDP单晶表面(001)由单点金刚石车削获得。刀具参数见表1,切削方式为端面车削,选取切削参数见表2

    表  1  圆弧刃天然单晶金刚石车刀参数
    Table  1.  Parameters of diamond turning tool
    nose radius r/mmedge radius rn/nmrake angle γ0/(°)clearance angle α0/(°)
    3.215009
    下载: 导出CSV 
    | 显示表格
    表  2  切削参数及其轮廓均方根比较
    Table  2.  Cutting parameters and RMS of contour
    No.(group)spindle speed n/(r/min)feed rate f /(μm/r)depth of cut ap/μmRMS of contour Rq/nm
    140010317
    635
    937
    23008315
    1224
    1827
    350010525
    40024
    30026
    下载: 导出CSV 
    | 显示表格

    采样样本通过白光干涉仪测量获得,采样面积360 μm×360 μm;采样点数256×256;最小采样周期1.412 μm/pixel(360 μm/255 pixels),如图13所示。根据采样大小确定公式(5),(6)中参数取值Δ=1.412,L=360 μm,N=256。图1(a)左为例,能看到相对运动的走刀痕迹。轮廓纹理需提取垂直走刀方向上高度信息,如图1(a)右所示,轮廓均方根Rq表2

    图  1  白光干涉仪测量的变切削深度时晶体表面和轮廓
    Figure  1.  Crystal surface and contour at different cutting depth measured by a white light interferometer
    图  2  白光干涉仪测量的变进给量时晶体表面和轮廓
    Figure  2.  Crystal surface and contour at different feed rate measured by a white light interferometer
    图  3  白光干涉仪测量的主轴转速改变时晶体表面和轮廓
    Figure  3.  Crystal surface and contour at different spindle speed measured by a white light interferometer

    通过PSD计算获得取样域内频率值及其功率量。根据白光干涉仪测量参数,确定取样周期范围2.824~360 μm,频率取值范围0.002 8~0.354 μm−1,最小频率间隔为0.002 8 μm−1。因频率取值为最小间隔的整数倍,只得到取值范围内的固定频率及其功率量。为弥补轮廓上部分频率信息遗漏,连续小波变换时需对频率段进行提取及重构。式(1),(2)计算表面微观轮廓的功率谱密度大小,如图46所示。

    图  4  变切削深度时PSD计算结果
    Figure  4.  Results of PSD at different depth of cut
    图  5  变进给量时PSD计算结果
    Figure  5.  Results of PSD at different feed rate
    图  6  变主轴转速时PSD计算结果
    Figure  6.  Results of PSD at different spindle speed

    图4在转速n=400 r/min,进给量f=10 μm/r,切削深度ap分别为3、6、9 μm时PSD计算结果。可以看出,切削深度变化时,PSD值在频率0.011 2 μm−1处最大,说明该值是影响表面微观轮廓的主要频率,且随切深变大,功率量也相应增加。频率为0.005 6、0.027 7 μm−1存有少量功率量,是影响表面轮廓的次要频率。

    图5在转速n=300 r/min,切深ap=3 μm,进给量f分别为8、12、18 μm/r时PSD计算结果。可以看到,变进给量时主频值发生变化,进给量f为8、12 μm/r时主频率为0.008 4 μm−1f为18 μm/r时主频率为0.005 6 μm−1。这是因为轮廓细密程度随进给量减少变得细密,采样范围内纹理增多,反应在PSD上为主频值右移。图5(a)(b)也可看出,进给量减小时中频特征(0.011 2、0.013 8 μm−1)有向低频过渡的趋势,低频成为进给量变化时主要频率,其重构波形对表面轮廓影响程度慢慢变小。

    图6在切削深度ap=5 μm,进给量f =10 μm/r,主轴转速n分别为500、400、300 r/min时PSD计算结果。频率0.011 2 μm−1是转速变化时影响微观轮廓的主要频率。图6(a)(b)(c)具有相同的频率分布,表明转速不会影响频率值变化,只影响功率量大小。说明这三个轮廓拥有近似频率特征,其重构波形有相似性。

    超精加工表面形貌可看作是多尺度波纹无序叠加组成。图46选取较为明显的频率值进行频率特征提取及重构。为保证分析过程及结果准确一致,变切削深度和变转速条件下选择低、中、高频率值分别为0.005 6,0.011 2,0.027 7 μm−1;变进给量情况下提取及重构频率值分别为0.005 6,0.008 4,0.011 2,0.013 8,0.027 7 μm−1。式(5),(6)中小波参数及尺度计算结果见表3

    表  3  尺度因子及其计算过程参数
    Table  3.  Scale a and its parameters of calculations
    frequency f/μm−1size L/μmnumber Nwaviness central frequency fsmexihat central frequency fcsampling period Δscale a
    low 0.005 6 360 256 0.015 8 0.25 1.412 11.209 3
    0.008 4 0.023 7 7.472 9
    mid 0.011 2 0.031 2 5.674 8
    0.013 8 0.038 8 4.563 2
    high 0.027 7 0.077 9 2.272 8
    下载: 导出CSV 
    | 显示表格

    式(3),(4)得到的小波系数与表面轮廓并无实质联系。为获得轮廓频率特征信息需对小波系数逆变换,由式(7)重构微观轮廓频率特征。

    切削参数对微观轮廓影响主要体现在频率特征的波长及幅值上。图7为变切削深度时轮廓与实际频率特征进行比较。频率为0.011 2 μm−1重构的波形与原始轮廓最为接近,频率值为0.005 6、0.027 7 μm−1含有少量轮廓信息,与PSD分析结果一致。由图7的原始轮廓比较发现:其周期及振幅随切削深度改变而变化,具体表现为切削越深,轮廓周期明显,幅值变大。这与重构中频0.011 2 μm−1的波形特征保持一致。切削深度对高、低频轮廓影响较小,在幅值和波长上无明显变化。均方根Rq分别为17,35,37 nm,得到ap= 3 μm是本文中较为合适的切削深度。

    图  7  变切削深度时原始轮廓与实际频率特征比较
    Figure  7.  Comparison of surface profile and actual frequency characteristics at different depth of cut

    图8为变进给量时轮廓与实际频率特征进行比较。低频特征(0.008 4、0.005 6 μm−1)的大小及走势与原始轮廓最为接近,视为主频率。对比低频特征发现,在进给量增加时,频率值变小,其特征表现为波长变长,幅值增加,即进给量变化会影响主频率值变化,随进给量增大趋向低频,振幅变大。进给量变化对中、高频无规律变化。均方根Rq分别为15、24、7 nm,得到f=8 μm/r是本文较为适宜的进给量。

    图  8  变进给量时原始轮廓与实际频率特征比较
    Figure  8.  Comparison of surface profile and actual frequency characteristics at different feed rate

    图9为变转速时原始轮廓与实际频率特征进行比较。在变转速情况下,中频特征(0.011 2 μm−1)与原始轮廓走势一致。高频特征(0.027 7 μm−1)反映了轮廓的细节变化,与加工过程中的振动和材料自身性质有关。转速变化对高、低频无明显影响。均方根Rq分别为25,24,26 nm,得到n=400 r/min是本文的理想转速。

    图  9  变主轴转速时原始轮廓与实际频率特征比较
    Figure  9.  Comparison of surface profile and actual frequency characteristics at different spindle speed

    通过分析切削参数变化对KDP晶体已加工表面实际频率特征的影响,可得到以下结论:(1)功率谱密度结合连续小波的方法能够获得采样域内的频率特征,将表面微观轮廓信息按实际频率进行提取及重构。获得影响微观形貌的主频率,建立加工参数与实际频率特征的关系;(2)切削参数对表面微观形貌的影响归结于实际频率特征上。中频特征反映了切深与主轴转速信息;进给量变化会影响低频值发生改变,频率特征随进给量增加时表现为波长变大、幅值增加。通过实际频率特征分析对评价加工工艺和光学晶体表面质量有着重要指导意义;(3)PSD仅能获得取样长度范围内取样步长整数倍的频率,但表面形貌上空间频率是连续分布的。通过连续小波重构可以在某种程度上弥补PSD方法在分析频率特征连续性的不足,提取出微观轮廓上的实际频率特征。

  • 图  1  三环变迹镜示意图

    Figure  1.  Diagrammatic layout of tri-annuli apodizing mask (unit: mm)

    图  2  Wz)的计算值及其拟合曲线比较

    Figure  2.  Calculated value of W(z) compared with its fitting curve

    图  3  三种变迹镜加权函数与Wz)=cz5/3的一致性对比

    Figure  3.  Weighting function of three kinds of apodizing masks compared with W(z)=cz5/3

    图  4  三种变迹镜的等晕角相对误差对比

    Figure  4.  Comparison of relative error of three kinds of apodizing masks

    图  5  不同模型下等晕角随高度变化图

    Figure  5.  Isoplanatic angle vary with altitude in different models

    图  6  高精度等晕角测量仪组成及布局

    Figure  6.  Composition and layout of high-precision measurement technique of isoplanatic angle

    图  7  典型条件下的全天等晕角测量数据

    Figure  7.  Measured data through the day in typical conditions

    表  1  变迹镜结构尺寸

    Table  1.   Size of tri-annuli apodizing mask

    type of apodizing mask radius /mm
    center block35.560
    telescope101.60
    inside/outside radius of inner transmittance ring37.3888/43.8404
    inside/outside radius of middle transmittance ring62.8904/69.2404
    inside/outside radius of outer transmittance ring81.9404/101.6000
    下载: 导出CSV

    表  2  变迹镜结构尺寸

    Table  2.   The size of apodizing mask

    design procedures of
    apodizing mask
    fit coefficient c of
    weighting function
    transmittance area A
    of the ring /m2
    inside/outside radius of middle
    transmittance ring/mm
    single aperture 1.209×10−16 0.0095 radius 55.00
    dual ring 1.551×10−16 0.0205
    inside/outside radius of inner transmittance ring 37.50/50.00
    inside/outside radius of outer transmittance ring 70.00/101.6
    下载: 导出CSV

    表  3  高精度等晕角测量仪的主要技术指标

    Table  3.   Key specifications of the high-precision isoplanometer

    key specificationsdesign value
    dynamic range3~30 μrad(λ=500 nm)
    wave length600~1000 nm
    sample frequency≥300 Hz
    single sample size≥600
    sampling precision12 bit
    variation range of zenith angle5°~45°
    variation range of azimuth0°~360°
    下载: 导出CSV
  • [1] Fried D L. Limiting resolution looking down through the atmosphere[J]. Journal of the Optical Society of America, 1966, 56(10): 1380-1384. doi: 10.1364/JOSA.56.001380
    [2] Loos G C, Hogge C B. Turbulence of the upper atmosphere and isoplanatism[J]. Applied Optics, 1979, 18(15): 2654-2661. doi: 10.1364/AO.18.002654
    [3] Walters D L. Propagation through atmospheric turbulence[M]//High Energy Laser Propagation Handbook, Chapter 5. Ann Arbor: Optimetrics Inc., 1983.
    [4] Stevens K B. Remote measurement of the atmospheric isoplanatic angle and determination of refractive turbulence profiles by direct inversion of the scintillation amplitude covariance function with Tikhonov regularization[D]. Monterey: Monterey Naval Postgraduate School, 1985.
    [5] Eaton F D, Peterson W A, Hines J R, et al. Isoplanatic angle direct measurements and associated atmospheric conditions[J]. Applied Optics, 1985, 24(19): 3264-3273. doi: 10.1364/AO.24.003264
    [6] Walters D L, Bradford L W. Measurements of r0 and θ0: two decades and 18 sites[J]. Applied Optics, 1997, 36(30): 7876-7886. doi: 10.1364/AO.36.007876
    [7] 杨高潮, 刘晓春, 范承玉, 等. 大气相干长度与等晕角的测量[J]. 强激光与粒子束, 1994, 6(2):215-220. (Yang Gaochao, Liu Xiaochun, Fan Chengyu, et al. Measurement of atmospheric coherence length and isoplanatic angle[J]. High Power Laser and Particle Beams, 1994, 6(2): 215-220
    [8] 中国人民解放军63655部队. 一种双环变迹镜及其制备方法和其测量等晕角的方法: 201310277603.1[P]. 2016-03-23.

    No. 63655 Troops of the Chinese People S Liberation Army. Dual-ring apodization lens, manufacturing method thereof and method for measuring isoplanatic angle through dual-ring apodization lens: 201310277603.1[P]. 2016-03-23).
    [9] 中国人民解放军63655部队. 一种基于三环变迹镜的等晕角测量系统及方法: 201310277560.7[P]. 2016-01-06.

    No. 63655 Troops of the Chinese People's Liberation Army. Three-ring apodization filter based system and method used for measuring isoplanatic angles: 201310277560.7[P]. 2016-01-06).
    [10] Tyson R K. Principles of adaptive optics[M]. 4th ed. Boca Raton: CRC Press, 2016.
    [11] Tyson R K. Introduction to adaptive optics[M]. Bellingham: SPIE Press, 2000.
  • 期刊类型引用(0)

    其他类型引用(6)

  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  1002
  • HTML全文浏览量:  287
  • PDF下载量:  67
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-06-01
  • 修回日期:  2021-08-02
  • 网络出版日期:  2021-08-14
  • 刊出日期:  2021-08-15

目录

/

返回文章
返回