CN102620917A - 透射式光学元件光致热变形像质分析方法 - Google Patents
透射式光学元件光致热变形像质分析方法 Download PDFInfo
- Publication number
- CN102620917A CN102620917A CN2012101034824A CN201210103482A CN102620917A CN 102620917 A CN102620917 A CN 102620917A CN 2012101034824 A CN2012101034824 A CN 2012101034824A CN 201210103482 A CN201210103482 A CN 201210103482A CN 102620917 A CN102620917 A CN 102620917A
- Authority
- CN
- China
- Prior art keywords
- optical element
- optical
- image quality
- thermal deformation
- deformation image
- 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
Abstract
本发明属于光学仪器热光学领域,涉及一种透射式光学元件光致热变形像质分析方法。该方法把大气压、功率、波长、折射率温度系数等因素考虑在内,通过有限元软件(I-DEAS)建立热分析模型,考虑各个面的辐射换热和传导换热,求解出光学元件温度场分布和各个节点的坐标位置与变形量,再通过Zernike多项式实现光学表面面型拟合,并由光学系统本身热源特点出发,提出一种新的圆柱坐标拟合方法用于实现温度场的精确拟合,由此计算出不同通光孔径位置光线通过的几何路程及材质的折射率函数,从而得出波前畸变及其PV、RMS值。本方法实现简单、精确度高,程序通用性好,适用于任何透射式光学元件波像差求取。
Description
技术领域
本发明属于光学仪器热光学领域,涉及一种透射式光学元件光致热变形像质分析方法。
背景技术
近年来,随着科学技术的发展,一些学者越来越重视温度对光学系统的影响,不仅考虑了温度变化产生的面形畸变,而且也把温度引起的折射率梯度因素考虑在内。常用光学设计软件(Code-V,Zemax等)一般认为光学元件温度为均匀分布,不能任意设定温度梯度分布。美国Sigmadyne公司的商业软件SigFit和RMR Design Group公司的PC Fringe,这两个软件使用了相同的数据传输技术,即利用了干涉检验接口,并通过干涉图文件导入光学设计软件。此种方法是在随机的光线追迹过程中计算光程变化,因此具有局限性。而对于含有吸收型光强调制光学元件的光学系统来说,元件内温度梯度变化很大,将产生明显的折射率梯度,从而影响系统光束质量。
因此,本发明目的在于提供一种更加精确快捷计算波前畸变的热光学分析方法。能够更加精确的仿真温度变化后系统波前畸变,为光机系统设计提供准确的反馈信息,因此具有重要的意义。
发明内容
当光学系统受光照射时,光学元件内部将产生不均匀的温度场,从而引起光学镜面热变形及光学元件折射率的变化。因此需要对光学系统进行热光学分析,得出光学系统受热后产生的光波前畸变。本发明目的在于提供一种更加精确计算波前畸变的热光学分析方法。
当光束射入到透射式光学元件各个面时,由于介质对光能的吸收,光束的强度将随着射入介质深度的增加而逐渐减弱,光能转化为热能,将产生温度梯度变化从而导致光学元件折射率的变化。考虑上述条件,本发明将通过有限元软件(I-DEAS)建立热分析模型,不仅考虑各个面的辐射换热和传导换热,而且也考虑了光强随射入介质深度变化等因素,计算每个单元的热流密度得出光学元件温度场分布并通过热弹性计算得到各个节点的坐标位置与变形量,再通过Zernike多项式实现光学表面面型拟合,并由光学系统本身热源特性出发,提出一种新方法用于实现温度场的精确拟合,即建立以光轴方向为z轴的圆柱坐标下的温度基底函数系。该表达式是将Zernike多项式的系数变换为沿光轴方向位置的函数,如式1:
使用该式对节点温度进行最小二乘法拟合,即能够得到某一光轴位置截面温度分布解析表达式,也能够得到某一孔径位置沿轴方向温度分布解析表达式。本方法在不同系数个数情况下拟合精度均优于现在用于温度场拟合的方法,尤其是拟合系数个数较多(>100)时,本文方法能够达到极高的拟合精度。由此计算出不同通光孔径位置光线通过的几何路程及材质的折射率函数,从而得出波前畸变及其PV、RMS值。本方法能够更加精确地衡量光学系统出射光束质量,且适用于任何透射式光学元件。
本发明能够解决的技术问题和积极效果是:
1、该方法是从热源本身出发,将Zernike多项式的系数变换为沿光轴方向位置的函数,也能够得到某一孔径位置沿轴方向温度分布解析表达式。用此方法拟合温度场更加精确,且方法简单,易于实现。
2、适用于任何透射式光学元件波前畸变的求取,适用性强。
具体实施方式
1.根据光学系统结构设计方案,建立有限元模型。
2.由光源功率及材质吸收率计算有限元各单元热流密度作为热分析软件(I-DEAS)的输入边界条件,考虑光学元件与周围环境的辐射换热,生成热分析模型求解系统稳态温度分布,并通过热弹性计算得到各节点位置变形量。
3.利用节点变形后坐标,拟合各光学元件变形后表面面型。
4.利用节点温度,由光学系统本身热源特点出发,使用式1圆柱坐标温度基底函数系实现温度场的精确拟合。
5.在孔径内均匀采样,在各孔径位置拟合沿轴折射率曲线,求取数值积分,得到该孔径位置的光程。最终合成系统出射波面波前畸变,求取PV、RMS值。
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101034824A CN102620917A (zh) | 2012-04-11 | 2012-04-11 | 透射式光学元件光致热变形像质分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101034824A CN102620917A (zh) | 2012-04-11 | 2012-04-11 | 透射式光学元件光致热变形像质分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102620917A true CN102620917A (zh) | 2012-08-01 |
Family
ID=46560967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012101034824A Pending CN102620917A (zh) | 2012-04-11 | 2012-04-11 | 透射式光学元件光致热变形像质分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102620917A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902857A (zh) * | 2012-10-08 | 2013-01-30 | 华中科技大学 | 一种光学元件支撑参数的蒙特卡洛分析方法及系统 |
CN108919486A (zh) * | 2018-08-20 | 2018-11-30 | 南京理工大学 | 低温真空环境下光学窗口的热光学特性表征与优化方法 |
CN109883993A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN109883994A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN109883996A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN110095858A (zh) * | 2018-12-12 | 2019-08-06 | 中国科学院紫金山天文台 | 自适应光学变形镜弹性模态像差表征方法 |
CN111158140A (zh) * | 2019-12-26 | 2020-05-15 | 哈尔滨新光光电科技股份有限公司 | 一种用于高能激光光学系统的杂散光分析抑制方法及高能激光光学系统 |
CN115452326A (zh) * | 2022-09-02 | 2022-12-09 | 上海索辰信息科技股份有限公司 | 一种基于无序折射率的热光效应分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6104537A (en) * | 1998-01-07 | 2000-08-15 | Olympus Optical Co., Ltd. | Image-forming optical system having a decentered optical element in a first lens group |
EP1413870A1 (en) * | 2001-07-05 | 2004-04-28 | Nikon Corporation | Optical member for optical lithography and evaluation method therefor |
CN102200432A (zh) * | 2010-03-23 | 2011-09-28 | 富士胶片株式会社 | 非球面物体测量方法和设备 |
TW201140019A (en) * | 2010-02-25 | 2011-11-16 | Nikon Corp | Measuring method and measuring apparatus of pupil transmittance distribution, exposure method and exposure apparatus, and device manufacturing method |
-
2012
- 2012-04-11 CN CN2012101034824A patent/CN102620917A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6104537A (en) * | 1998-01-07 | 2000-08-15 | Olympus Optical Co., Ltd. | Image-forming optical system having a decentered optical element in a first lens group |
EP1413870A1 (en) * | 2001-07-05 | 2004-04-28 | Nikon Corporation | Optical member for optical lithography and evaluation method therefor |
TW201140019A (en) * | 2010-02-25 | 2011-11-16 | Nikon Corp | Measuring method and measuring apparatus of pupil transmittance distribution, exposure method and exposure apparatus, and device manufacturing method |
CN102200432A (zh) * | 2010-03-23 | 2011-09-28 | 富士胶片株式会社 | 非球面物体测量方法和设备 |
Non-Patent Citations (4)
Title |
---|
丁延卫 等: "空间光学窗口的热光学灵敏度分析", 《光电工程》, vol. 29, no. 5, 31 October 2002 (2002-10-31), pages 15 - 18 * |
巩盾 等: "离轴三反射系统的热光学分析和温控指标的制定", 《光学精密工程》, vol. 19, no. 6, 30 June 2011 (2011-06-30), pages 1213 - 1220 * |
罗传伟 等: "光学系统折射率温度效应的模拟计算", 《应用光学》, vol. 29, no. 2, 31 March 2008 (2008-03-31), pages 234 - 239 * |
陈恩涛 等: "空间光学遥感器光-机-热集成分析方法研究", 《宇航学报》, vol. 26, no. 1, 31 January 2005 (2005-01-31), pages 66 - 70 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902857A (zh) * | 2012-10-08 | 2013-01-30 | 华中科技大学 | 一种光学元件支撑参数的蒙特卡洛分析方法及系统 |
CN108919486A (zh) * | 2018-08-20 | 2018-11-30 | 南京理工大学 | 低温真空环境下光学窗口的热光学特性表征与优化方法 |
CN108919486B (zh) * | 2018-08-20 | 2021-03-26 | 南京理工大学 | 低温真空环境下光学窗口的热光学特性表征与优化方法 |
CN110095858A (zh) * | 2018-12-12 | 2019-08-06 | 中国科学院紫金山天文台 | 自适应光学变形镜弹性模态像差表征方法 |
CN110095858B (zh) * | 2018-12-12 | 2021-06-08 | 中国科学院紫金山天文台 | 自适应光学变形镜弹性模态像差表征方法 |
CN109883993A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN109883994A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN109883996A (zh) * | 2019-03-05 | 2019-06-14 | 中国计量大学 | 基于哈特曼光线追迹的非均匀介质场的测量方法 |
CN111158140A (zh) * | 2019-12-26 | 2020-05-15 | 哈尔滨新光光电科技股份有限公司 | 一种用于高能激光光学系统的杂散光分析抑制方法及高能激光光学系统 |
CN111158140B (zh) * | 2019-12-26 | 2022-03-08 | 哈尔滨新光光电科技股份有限公司 | 一种用于高能激光光学系统的杂散光分析抑制方法及高能激光光学系统 |
CN115452326A (zh) * | 2022-09-02 | 2022-12-09 | 上海索辰信息科技股份有限公司 | 一种基于无序折射率的热光效应分析方法 |
CN115452326B (zh) * | 2022-09-02 | 2023-08-15 | 上海索辰信息科技股份有限公司 | 一种基于无序折射率的热光效应分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102620917A (zh) | 透射式光学元件光致热变形像质分析方法 | |
Xu et al. | Influence of altitude on aero-optic imaging deviation | |
Karpeev et al. | Fibre sensors based on transverse mode selection | |
van Hinsberg et al. | Density measurements using near-field background-oriented schlieren | |
CN107421436B (zh) | 基于空间光调制器参考面的非球面干涉测量系统及方法 | |
Zhu et al. | A new optical evaluation approach for parabolic trough collectors: first-principle optical intercept calculation | |
CN101290218B (zh) | 用于校正非球面非零位检测时的原理误差的方法 | |
CN101980056A (zh) | 基于m2因子和光闪烁指数确定大气湍流参数方法及装置 | |
CN115421298B (zh) | 复杂激光内通道热效应影响下的光链路传输计算方法 | |
Tao et al. | Annealing of compression molded aspherical glass lenses | |
Xiao et al. | Joint influences of aerodynamic flow field and aerodynamic heating of the dome on imaging quality degradation of airborne optical systems | |
Camillo et al. | Combined Experimental and Numerical Investigation of a Hypersonic Turbulent Boundary Layer by Means of FLDI and Large-Eddy Simulations | |
Niewiem et al. | Variation of structured laser beam pattern and optimization for an alignment reference line creation | |
Wang et al. | On the accuracy of Malley probe measurements of aero-optical effects: a numerical investigation | |
Gavril'eva et al. | Approaches to cross-talk noise reduction in modal holographic wavefront sensors | |
CN207351658U (zh) | 一种楔形透镜透射波前的测量装置 | |
CN105301035A (zh) | 光刻投影物镜热效应评估方法 | |
CN103984094A (zh) | 一种光学系统热性能仿真方法 | |
Dixit et al. | Systematic characterization of near-index-matched optics based atmospheric turbulence simulator | |
Yu | Image processing methods based on physical models | |
Hu et al. | Circular Subaperture Stitching Interferometry Based on Polarization Grating and Virtual–Real Combination Interferometer | |
Gao et al. | An analytic expression for the beam wander of a Gaussian wave propagating through scale-dependent anisotropic turbulence | |
CN104977156A (zh) | 一种用于CCD远场法光束质量β因子测量的校准系统 | |
CN105021292A (zh) | 一种获取同步辐射部分相干光的传播特性的方法 | |
CN204924616U (zh) | 一种用于CCD远场法光束质量β因子测量的校准系统 |
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: 20120801 |