气动热环境下共形整流罩热辐射特性研究

薛文慧, 王惠, 包春慧, 范志刚

薛文慧, 王惠, 包春慧, 范志刚. 气动热环境下共形整流罩热辐射特性研究[J]. 应用光学, 2017, 38(6): 999-1005. DOI: 10.5768/JAO201738.0606003
引用本文: 薛文慧, 王惠, 包春慧, 范志刚. 气动热环境下共形整流罩热辐射特性研究[J]. 应用光学, 2017, 38(6): 999-1005. DOI: 10.5768/JAO201738.0606003
Xue Wenhui, Wang Hui, Bao Chunhui, Fan Zhigang. Thermal radiation characteristics of conformal dome in aero-dynamic environment[J]. Journal of Applied Optics, 2017, 38(6): 999-1005. DOI: 10.5768/JAO201738.0606003
Citation: Xue Wenhui, Wang Hui, Bao Chunhui, Fan Zhigang. Thermal radiation characteristics of conformal dome in aero-dynamic environment[J]. Journal of Applied Optics, 2017, 38(6): 999-1005. DOI: 10.5768/JAO201738.0606003

气动热环境下共形整流罩热辐射特性研究

基金项目: 

国家自然科学基金 61275020

航天科技创新基金 JTBL11610

航空科学基金 20130177002

详细信息
    作者简介:

    薛文慧(1993—),女,山西朔州人,硕士研究生,主要从事气动光学方面的研究工作。E-mail:1392181520@qq.com

    范志刚(1966—),男,黑龙江哈尔滨人,教授,博导,主要从事气动光学、光电检测方面的研究工作。E-mail:fzg@hit.edu.cn

  • 中图分类号: TN206;O435.1

Thermal radiation characteristics of conformal dome in aero-dynamic environment

  • 摘要: 高温整流罩产生强烈的红外辐射,在探测器接收面上形成背景噪声,严重影响成像质量。为了评估气动热环境下高速飞行器共形整流罩热辐射对探测器性能的影响,建立了气动热环境下整流罩热流固耦合计算模型、整流罩热辐射发射以及传输模型,计算了整流罩非均匀温度场、形变场、应力以及应变场,仿真得到了共形整流罩热辐射对高速飞行器光学系统成像质量的影响。研究结果表明:随着时间的不断增长,探测器接收面的最大辐照度也逐渐增大。在第10 s时,0°攻角下共形整流罩热干扰辐照度最小值为0.094 W/m2,最大值为0.108 W/m2。和相同工况下的共形整流罩相比,球形整流罩10 s时的最大辐照度分别是椭球面和抛物面形整流罩最大辐照度的近12倍和7倍,即共形整流罩产生的干扰对探测系统的影响较小,不会对探测器的探测性能产生致命性的影响。
    Abstract: The high-temperature dome produces intense infrared radiation, which forms the background noise on the detector receiving surface and seriously affects the image quality. To assess the performance degradation as a result of aerodynamic heating of the conformal dome flying at high speed, this paper established the multi-physics coupling model of conformal domes under aerodynamic heating environment, the dome thermal radiation emission and transmission model. Based on these models, the temperature field, the stress field, the strain field and the deformation field of conformal domes were obtained, and the imaging quality of high-speed aircraft affected by heat radiation of conformal dome was analyzed. Results show that as time goes on, the maximum irradiance of the detector is increasing. At 10 s, the minimum value of the thermal interference irradiance of the conformal shroud at 0° attack angle is 0.094 W/m2, and the maximum value is 0.108 W/m2. Compared with conformal domes under the same working condition, the maximum irradiance of the spherical dome at 10 s is about 12 times and 7 times of the maximum irradiance of ellipsoid and parabolic dome, respectively. Conditions thermal radiation interference of the conformal dome has less effect on the detection system comparing with spherical dome, and does not have a fatal effect on the detection performance of the detector.
  • 光干涉测量方法具有非接触、高效率、高灵敏度等特点,也就必然成为新世纪精密测量技术的重要发展方向[1]。干涉图相位提取是干涉测量中的关键一步,相位提取精度的优劣直接影响最终检测的精度[2]。快速傅里叶变换法(fast fourier transform,FFT)具有处理瞬变波面相位的能力,可有效克服机械振动和气流的影响[3-5],因此成为大口径光学元件测量的常用方法之一。

    快速傅里叶变换法(FFT)于20世纪80年代被提出并应用于条纹图的相位提取中[6],对干涉条纹图进行FFT相位提取法的目的是提取出频谱中包含被测波面信息的正一级频谱,也称为滤波过程[7-8],滤波过程的好坏直接影响相位提取精度的高低。目前国内外文献针对FFT相位提取技术主要围绕算法中几个关键步骤进行研究,如文献[9]对干涉图延拓进行分析;文献[10]提出了一种基于FFT时移特性的叠栅条纹细分方法;文献[11]较为综合地分析了影响FFT相位提取法的边缘误差、窗函数、滤波器设计、干涉图延拓及载波条纹数等因素;文献[12]将傅里叶变换相位提取法应用于仿真及实验采集电子散斑干涉(electronic speckle interference,ESPI)信息处理中。为得到更高的相位提取精度,对二维FFT相位提取法的进一步研究仍是国内外研究热点。

    针对FFT相位提取技术的滤波过程进行研究,通过选取不同滤波窗对计算机仿真干涉图进行相位提取,对比不同类型的滤波窗滤出的正一级频谱分布形状及其对相位提取精度的影响,最后利用实验采集干涉图对算法的可靠性进行验证。

    对于一般的实验干涉仪,引入空间载频${f_x}$${f_y}$后的干涉条纹的强度分布可以表示为

    $$\begin{split} f(x,y) = & {I_d}(x,y) + c(x,y) \cdot {{\rm e}^{{\rm i}[{\rm{2}}{\text{π}} {f_x}x + {\rm{2}}{\text{π}}{f_y}y]}} + {c^ * }(x,y) \cdot \\ & {{\rm e}^{ - {\rm i}[{\rm{2}}{\text{π}} {f_x}x + {\rm{2}}{\text{π}} {f_y}y]}} \\ \end{split} $$ (1)

    其中

    $$c(x,y) = {1 / 2} \cdot b(x,y){{\rm e}^{{\rm i}[\varphi (x,y)]}} $$ (2)

    式中:${I_d}(x,y)$是干涉条纹的背景光强分布;$b(x,y)$是干涉条纹的调制度分布;$\varphi (x,y)$为含有待测波面相位信息的相位分布函数;$ * $表示复共轭[13]

    对(1)式进行傅里叶变换,将空域中的干涉条纹信息转移到频域内,得到的频谱分布函数表示为

    $$\begin{split} {\rm{F(}}{f_{\rm{1}}}{\rm{,}}{f_{\rm{2}}}{\rm{)}} = & {\rm{A(}}{f_{\rm{1}}}{\rm{,}}{f_{\rm{2}}}{\rm{)}} + C({f_{\rm{1}}} - {f_x},{f_{\rm{2}}} - {f_y}) + \\ & C*({f_{\rm{1}}} + {f_x},{f_{\rm{2}}} + {f_y}) \end{split} $$ (3)

    式中$A({f_1},{f_2})$$C({f_1} - {f_x},{f_2} - {f_y})$$C*({f_1} + {f_x},{f_2} + {f_y})$分别表示零级、正一级、负一级频谱分布函数。选取一个中心频率为$({f_x},{f_y})$的滤波器将正一级频谱$C({f_1} - {f_x},{f_2} - {f_y})$提取出来并平移至原点,得到$C({f_x},{f_y})$,对其进行二维傅里叶逆变换,即可得到$C(x,y)$。由(2)式可得所求相位

    $$\varphi (x,y) = \arctan \frac{{{\rm{lm}}[c(x,y)]}}{{{\rm Re} [c(x,y)]}} $$ (4)

    式中${\rm Re} [c(x,y)]$${\rm{lm}}[c(x,y)]$分别为$C(x,y)$的实部和虚部。

    图1(a)为仿真干涉图的频谱分布,要准确无误地从该频谱图中提取正一级频谱,就要选取适当的滤波窗,但若选取滤波窗不当,就会造成有效频谱信息的缺失[14]。如图1(b)所示为选取滤波窗为圆域形窗口时滤出的正一级频谱,其频谱信息出现明显截断现象,为了减少截断误差,有必要对几种通用的滤波器进行分析与比较。目前二维FFT进行干涉测试数据处理时,使用较多的滤波器有海明(Hamming)窗、汉宁(Hanning)窗、高斯(Gaussian)窗、布莱克曼(Blackman)窗。

    图  1  FFT法滤波阶段频谱分布
    Figure  1.  Filtering stage spectrum distribution of FFT method

    根据干涉图频谱分布的这一特点,选取与其分布形状相似的几种滤波窗函数,表1为几种滤波窗函数的一维表达式。

    表  1  几种滤波窗函数
    Table  1.  Several filter window functions
    Window typeFunctionParameter
    Blackman$\omega (n) = \left\{ \begin{array}{l} 0.42 - 0.5\cos (\frac{{2{\text{π}} n}}{{N - 1}}) + \\ 0.08\cos (\frac{{4{\text{π}} n}}{{N - 1}})\quad \quad \quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $main lobe width ${\rm{12}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 56 dB
    Gaussian$\begin{array}{l} \omega (n) = \;\exp \left[ { - \frac{1}{2}{{\left( {\frac{{n - N/2}}{{\sigma N/2}}} \right)}^2}} \right]\quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \sigma \leqslant 0.5 \\ \end{array} $main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 44 dB
    Hanning$\omega (n) = \left\{ \begin{array}{l} 0.5\left[ {{\rm{1}} - \cos (\frac{{2{\text{π}} n}}{{N - 1}})} \right]\quad \quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 31 dB
    Hamming$\omega (n) = \left\{ \begin{array}{l} 0.5{\rm{4 - 0}}{\rm{.46}}\cos (\frac{{2{\text{π}} n}}{{N - 1}})\quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 41 dB
    下载: 导出CSV 
    | 显示表格

    根据表1中各类窗函数的分布函数绘制时频振幅响应分布,选取N=30作为采样点,为直观起见,绘制其一维空/频分布如图2,其对应二维空/频分布,可由一维空/频分布经旋转得到。

    图  2  几种滤波窗函数振幅响应
    Figure  2.  Spectral response of several filter window functions

    根据图2多种滤波窗函数的时域/频域响应及表1中参数可知,不同窗函数主要差别在于主瓣宽度及旁瓣衰减速率。其中主瓣宽度主要影响信号的能量分布,旁瓣峰值衰减速率影响能量的泄露程度[15],旁瓣越高,能量泄露越严重,衰减速率就越慢[16]。对于FFT滤波过程,要保证滤出正一级频谱的同时滤掉其他频谱成分,且要最大程度地保留正一级频谱沿x轴和y轴延伸区域内的频谱成分,需综合考虑滤波窗的频谱分布形状及对应参数。

    针对第一节提出的频谱截断现象,选取Hamming、Hanning、Gaussian、Blackman四种与正一级频谱分布形状相近的滤波窗对仿真干涉图的频谱图(图1(a))进行滤波处理,其中仿真干涉图尺寸为300×300像素,零频频谱中心频率位于(151,151),正一级频谱中心频率位于(164,164),即滤波窗中心频率确定。根据距离公式求得二者像素间距为18.384 7像素,则滤出正一级频谱对应窗口函数的滤波半径可以此为参考。如图3分别为其对应频谱响应(滤波半径取值为18)及选取该窗口滤出的正一级频谱分布图。

    分析图3(a)(d)可得:经Hamming窗和Hanning窗提取的正一级频谱较完整地保留了原始频谱成分,没有其他频谱成分的混入。分析其造成的原因主要是由于主瓣宽度覆盖了正一级频谱带宽,且第一旁瓣高度较低,未出现频谱泄露现象,即未混入零频及负一级频谱;而经Gaussian、Blackman窗提取的正一级频谱分布虽较完整地保留了原始频谱成分,但有其他频谱成分的混入。分析对应窗口函数及频谱响应分布得到,其主瓣宽度覆盖正一级频谱宽度,但由于其第一旁瓣有一定幅值,造成了在此窗函数xy方向上引入较多频谱泄露,即混入的零频及正一级频谱。

    图  3  不同滤波窗及其提取模拟干涉图的正一级频谱分布
    Figure  3.  Positive first-order spectral distribution of simulated interferogram extracted by different filtering windows

    将不同滤波窗函数提取出的正一级频谱进行移中并做逆傅里叶变换,利用反正切函数求得包裹相位,再利用离散余弦(discrete-cosin-transform,DCT)解包裹法提取连续相位,最后进行波面拟合,得到波面分布。计算不同滤波窗口选取下的FFT相位提取法得到波面的波面峰谷值(peak-to-valley,PV)和波面均方根值(root-mean-square,RMS),分析相位提取精度并与原始波面进行残差计算。其中原始波面分布评价参数PVa为0.246 4λ,RMSa为0.056 4λ,经FFT法得到的波面分布评价参数为PV1、RMS1,记ΔPV=PV1−PVa,ΔRMS=RMS1−RMSa。其值如表2所示。

    表  2  不同滤波窗下FFT算法结果比较 (λ=632.8 nm)
    Table  2.  Comparison of FFT algorithm results in different filtering windows (λ=632.8 nm)
    Window TypePV1/λ∆PV/λRMS1/λ∆RMS/λ
    Hamming0.257 90.008 50.056 3−0.000 1
    Hanning0.258 60.012 20.056 3−0.000 1
    Gaussian0.262 00.015 60.056 3−0.000 1
    Blackman0.248 30.001 90.056 2−0.000 2
    下载: 导出CSV 
    | 显示表格

    对比上述ΔPV及ΔRMS可知,选取Blackman窗及Hamming窗处理得到波面的评价参数在数值上较其他窗口函数更逼近原始波面的评价参数,结果优于Gaussian窗及Hanning窗。其中最优窗口类型为Hamming窗,其波面峰谷值残差为0.008 5λ,波面均方根值残差为0.000 1λ。

    为验证仿真结果的可行性和精度,分别采用二维傅里叶变换法(2D-FFT)和移相干涉法(PSI)对实验采集的干涉图进行处理并比较其测试结果。实验阶段所用的菲索干涉仪是美国Zygo公司生产的Verifire PE激光干涉仪,该干涉仪是共光路面形计量干涉仪,光源为低功率的632.8 nm的氦氖激光,采用PZT移相干涉原理可以对多种光学元件的面形进行检测,其干涉仪装置如图4所示。

    图  4  Zygo激光干涉仪装置图
    Figure  4.  Diagram of Zygo laser interferometer device

    使用Zygo干涉仪采集实验干涉图如图5(a)所示,其尺寸为300×300像素。经消除白噪声后的条纹图如图5(b),对处理后的条纹图进行二维傅里叶变换,得到频谱分布如图5(c)。选取Hamming窗将正一级频谱取出并滤掉其他频谱,并将正一级频谱平移到整个频谱中心,如图5(d)(e);对移中后的正一级频谱做二维傅里叶逆变换,利用反正切函数求取包裹相位如图5(f),采用离散余弦解包裹法进行处理,其结果如图5(g)

    图  5  FFT相位提取法各阶段处理结果图
    Figure  5.  Diagram of interferogram processing result by using FFT

    选取其他3种滤波窗(Hanning、Blackman、Gaussian)对实验采集条纹图进行上述相位提取的算法处理,图6(a)(d)所示为上述4种滤波窗提取出的实验干涉图正一级频谱分布。

    图6进行分析可得,实验干涉图在不同滤波窗处理下的正一级频谱分布与仿真干涉图的结果呈现相同的趋势。即Hamming窗和Hanning窗较完整地保留正一级频谱且无其他频谱成分混入,而Blackman窗及Gaussian窗仍有其他频谱成分混入。

    图  6  不同滤波窗提取实验干涉图的正一级频谱分布
    Figure  6.  Positive first-order spectral distribution of experimental interferogram extracted by different filtering windows

    表3为用FFT相位提取法(不同滤波窗口)得到的波面及PSI法测得波面的评价参数汇总。

    表  3  FFT(不同滤波窗)&PSI算法的结果比较 (λ=632.8 nm)
    Table  3.  Comparsion of results in FFT(different filtering windows) and PSI (λ=632.8 nm)
    MethodPV/λRMS/λ
    PSI0.445 40.002 7
    FFT (Hamming Window)0.436 10.002 2
    FFT (Hanning Window)0.417 70.002 8
    FFT (Gaussian Window)0.326 20.002 2
    FFT (Blackman Window)0.527 30.002 9
    下载: 导出CSV 
    | 显示表格

    表3可知,选取Hanning窗及Hamming窗处理得到波面的评价参数在数值上较其他窗口函数更逼近PSI法所得波面评价参数,结果优于Blackman窗及Gaussian窗。其中最优窗口类型为Hamming窗,其波面峰谷值残差为0.009 3λ,波面均方根值残差为0.000 5λ。

    针对二维傅里叶变换相位提取算法滤波过程中滤波窗的选取进行研究。首先分析干涉条纹图的频谱分布及各类滤波窗口函数的参数及其频谱响应;其次进行仿真实验,针对干涉图的频谱分布特点选取合适的窗口进行滤波处理并进行相位提取,分析不同窗口选取下对相位提取精度的影响;最后利用立式菲索干涉仪进行实验验证并与移相干涉测量法所测的结果进行对比。实验及仿真结果表明:选取Hamming窗提取的正一级频谱较完整地保留了原始频谱成分,相位提取精度优于0.01λ,可进一步应用于大口径光学元件的测量中。

  • 图  1   利用APDL建立椭球形整流罩

    Figure  1.   Flow chart of establishing ellipsoid fairing by using APDL

    图  2   椭球形整流罩第10 s温度场分布

    Figure  2.   Temperature field distribution of ellipsoid fairing 10 s

    图  3   椭球形整流罩第10 s等效应力场分布

    Figure  3.   Equivalent stress field distribution of the ellipsoid fairing 10 s

    图  4   椭球形整流罩第10 s等效应变场分布

    Figure  4.   Equivalent strain field distribution of ellipsoid fairing 10 s

    图  5   椭球形整流罩第10 s等效形变场分布

    Figure  5.   Equivalent deformation field distribution of ellipsoid fairing 10 s

    图  6   0°攻角2 s时椭球面整流罩热辐射在探测器接收面上3维及其对应的2维投影分布

    Figure  6.   On detector, 3-d and 2-d projection distribution of ellipsoid fairing thermal radiation at attack angle 0°and at time of 2 s

    图  7   0°攻角10 s时椭球形整流罩热辐射在探测器接收面上3维及其对应的2维投影分布

    Figure  7.   On detector, 3-d and 2-d projection distribution of ellipsoid fairing thermal radiation at attack angle 0°and at time of 10 s

    图  8   0°攻角1 s~10 s探测器接收面上最大辐照度变化趋势示意图

    Figure  8.   From 1~10 s, variation tendency of the biggest irradiance of detector surface at attack angle 0°

    图  9   10°攻角2 s时椭球形整流罩热辐射在探测器接收面上3维及其对应的2维投影分布

    Figure  9.   On detector, 3-d and 2-d projection distribution of ellipsoid fairing thermal radiation at attack angle 10°and at time of 2 s

    图  10   10°攻角10 s时椭球形整流罩热辐射在探测器接收面上3维及其对应的2维投影分布

    Figure  10.   On detector, 3-d and 2-d projection distribution of ellipsoid fairing thermal radiation at attack angle 10°and at time of 10 s

    图  11   10°攻角1 s~10 s探测器接收面上最大辐照度变化趋势示意图

    Figure  11.   From 1~10 s, variation tendency of the biggest irradiance of detector surface at attack angle 10°

    图  12   0°攻角3种典型面形共形整流罩热辐射随时间变化趋势图

    Figure  12.   Three typical surface shape conformal fairing heat radiation varies with time trend at attack angle 0°

    表  1   300 K蓝宝石晶体的物理特性

    Table  1   Physical properties of sapphire crystal at 300 K

    密度
    /(kg/m3)
    熔点/K 膨胀系数
    /(10-6K-1)
    比热容
    /(J/(kg·K))
    弹性模量
    /GPa
    导热系数
    /(W/(m·K))
    屈服强度
    /(MPa)
    泊松比
    3 980 2 050 5.3 750 344 36 300 0.27
    下载: 导出CSV
  • [1] 殷兴良.现代光学新分支学科——气动光学[J].中国工程科学, 2006, 7(12):1-6. doi: 10.3969/j.issn.1009-1742.2006.12.001

    Yin Xinliang.A new subdiscipline of contem porary optics-aero-optics[J].Enfineering Science, 2006, 7(12):1-6. doi: 10.3969/j.issn.1009-1742.2006.12.001

    [2] 马毅飞, 赵文平.窗口辐射对红外成像探测影响的研究[J].系统工程与电子技术, 2005, 27(3):427-430. doi: 10.3321/j.issn:1001-506X.2005.03.013

    Ma Yifei, Zhao Wenping. Effects of window radiation on infrared imaging detection[J].Systems Engineeting and Electronics, 2005, 27(3):427-430. doi: 10.3321/j.issn:1001-506X.2005.03.013

    [3]

    Burns S P, Howell J R, Klein D E. Assessment of the swartz-wendroff finite element formulation for combined mode heat transfer, radiative heat transfer current research[J]. ASME HTD, 2002, 276:1-7.

    [4]

    Brandon S, Derby J J. A finite element method for internal radiative heat transfer and its application to analysis of the growth of semitransparent crystals[J].ASME HTD, Fundamentals of Radiation Heat Transfer, 1992, 106:1-6.

    [5] 范志刚, 于翠萍, 肖昊苏, 等.气动光学头罩热辐射效应数值仿真研究[J].应用光学, 2012, 33(4):649-653. http://d.old.wanfangdata.com.cn/Periodical/yygx201204001

    Fan Zhigamg, Yu Cuiping, Xiao Haosu, et al.Numerical simulation of termal radiation effect about optical dome in aerodynamic thermal environment[J].Journal of Applied Optics, 2012, 33(4):649-653. http://d.old.wanfangdata.com.cn/Periodical/yygx201204001

    [6] 刘珂, 陈宝国, 李丽娟.空空导弹红外导引头技术发展趋势及关键技术[J].激光与红外, 2011, 41(10):1117-1121. doi: 10.3969/j.issn.1001-5078.2011.10.012

    Liu Ke, Chen Baoguo, Li Lijuan.Development tendency and key technology of IR seeker for air-to-air missile[J].Laser & Infrared, 2011, 41(10):1117-1121. doi: 10.3969/j.issn.1001-5078.2011.10.012

    [7]

    Frase B S, Hemingway A. High performance faceted domes for tactical and strategic missiles[J]. SPIE, 1994, 2286:486-492. http://cn.bing.com/academic/profile?id=ed22b60a788f39ad134eb870b4c73d05&encoded=0&v=paper_preview&mkt=zh-cn

    [8]

    Knapp D J, Mills J P, Hegg R G, et al. Conformal optics risk reduction demonstration[J]. SPIE, 2001, 4375:147. http://cn.bing.com/academic/profile?id=7de2efb40080c09dba9797a53e90a2b6&encoded=0&v=paper_preview&mkt=zh-cn

    [9]

    Shorey A B, Kordonski W, Tracy J, et al. Developments in the finishing of domes and conformal optics[J]. SPIE, 2007, 6545:65450Q1-65450Q9. http://cn.bing.com/academic/profile?id=abbf2cdb9a1e169e8d03e44586dc4971&encoded=0&v=paper_preview&mkt=zh-cn

    [10]

    Shannon R R. Overview of conformal optics[J]. SPIE, 1999, 3705:180-188. doi: 10.1117-12.354622/

    [11] 王勖成.有限单元法[M].北京:清华大学出版社, 2003:441-464.

    Wang Maocheng.Finite element method[M].Beijing:Tsinghua University Press, 2003:441-464.

    [12]

    Wakaki M, Kudo K, Shibuya T.Physical properties and data of optical materials[C]. Boca Raton: CRC Press, 2007: 383-535.

    [13]

    Visknata R, Anderson E E. Heat transfer semitransparent solids[J].Advances in Heat Transfer, 1975, 11:317-441. doi: 10.1016/S0065-2717(08)70077-7

    [14] 张琳.梯度折射率介质内辐射传递方程数值模拟的有限元法[D].哈尔滨: 哈尔滨工业大学学位论文, 2009.

    Zhang Lin.Finite element method for solving radiative transfer equation in graned index mediium[D].Haerbin: Harbin Institute of Technology, 2009: 32-63.

  • 期刊类型引用(7)

    1. 张瑜,张峰,郭芮,苏瑛,张云龙,许增奇,王阜超. 磁流变抛光二次非球面工件位置对刀校正. 应用光学. 2022(03): 518-523 . 本站查看
    2. 刘春阳,李洛克. 基于轮廓测量的集料表面粗糙度表征方法试验研究. 工业安全与环保. 2021(03): 15-19+25 . 百度学术
    3. 夏良静. 工艺参数对光学元件加工精度影响的数学模型研究. 激光杂志. 2020(08): 201-204 . 百度学术
    4. 张晋烨,张海云,赵玉刚,张军,王显达. 磨料振动对曲轴抛光的影响. 组合机床与自动化加工技术. 2019(01): 46-48 . 百度学术
    5. 高丰,李琦,向阳,宁旭. 高清折反式视频显微镜. 光学技术. 2019(02): 136-141 . 百度学术
    6. 刘红丽. 散粒磨料振动光饰技术专利分析. 科技风. 2019(14): 141 . 百度学术
    7. 马旭德,张东梅,尚春民. 固着磨料球面高速研磨加工技术研究. 新型工业化. 2018(04): 34-42 . 百度学术

    其他类型引用(4)

图(12)  /  表(1)
计量
  • 文章访问数:  943
  • HTML全文浏览量:  202
  • PDF下载量:  172
  • 被引次数: 11
出版历程
  • 收稿日期:  2016-09-19
  • 修回日期:  2017-05-23
  • 刊出日期:  2017-10-31

目录

/

返回文章
返回