CN113031201A - 一种在轨补偿主镜曲率半径误差影响的调整方法 - Google Patents
一种在轨补偿主镜曲率半径误差影响的调整方法 Download PDFInfo
- Publication number
- CN113031201A CN113031201A CN202110274445.9A CN202110274445A CN113031201A CN 113031201 A CN113031201 A CN 113031201A CN 202110274445 A CN202110274445 A CN 202110274445A CN 113031201 A CN113031201 A CN 113031201A
- Authority
- CN
- China
- Prior art keywords
- mirror
- wavefront aberration
- curvature radius
- orbit
- primary mirror
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000004075 alteration Effects 0.000 claims abstract description 83
- 201000009310 astigmatism Diseases 0.000 claims abstract description 31
- 239000013598 vector Substances 0.000 claims description 15
- 210000001747 pupil Anatomy 0.000 claims description 9
- 206010073261 Ovarian theca cell tumour Diseases 0.000 claims description 8
- 208000001644 thecoma Diseases 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 2
- 206010010071 Coma Diseases 0.000 abstract description 15
- 230000008859 change Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 14
- 230000035945 sensitivity Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 description 1
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B7/00—Mountings, adjusting means, or light-tight connections, for optical elements
- G02B7/18—Mountings, adjusting means, or light-tight connections, for optical elements for prisms; for mirrors
- G02B7/182—Mountings, adjusting means, or light-tight connections, for optical elements for prisms; for mirrors for mirrors
- G02B7/185—Mountings, adjusting means, or light-tight connections, for optical elements for prisms; for mirrors for mirrors with means for adjusting the shape of the mirror surface
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Telescopes (AREA)
- Lenses (AREA)
Abstract
一种在轨补偿主镜曲率半径误差影响的调整方法,涉及空间望远镜主动光学技术领域,解决现有大口径离轴三反空间望远镜发射入轨后,主镜曲率半径改变量超出公差范围时导致波前像差不满足天文观测要求的问题,空间望远镜在轨对宇宙中恒星进行成像获得星点图像进而得到在轨状态下系统波前像差;判断系统波前像差是否满足天文观测要求如果否,将当前系统波前像差值与地面装调状态下系统的波前像差值做差,获得主镜曲率半径误差引入的像散增量和彗差增量;计算主镜曲率半径误差;根据波前像差RMS的最小值,确定次镜调整量,获得系统波前像差RMS的综合贡献量,求取系统波前像差RMS最小值,在轨调整次镜及摆镜,调整完毕后,空间望远镜进入工作状态。
Description
技术领域
本发明涉及空间望远镜主动光学技术领域,具体为通过分析离轴三反空间望远镜在轨波前像差相对于地面装调状态波前像差改变特征来确定像差场的补偿策略。
背景技术
空间大口径离轴三反消像散天文望远镜是天文望远镜发展的重要方向之一。口径增大及主镜轻量化设计使主镜支撑面对诸多挑战:温度场等外界条件发生改变时,主镜曲率半径将发生改变;当曲率半径改变量超出公差范围时,曲率半径误差对离轴三反波前像差的影响将不能被忽略,进而对天文观测产生不利的影响。
当前,地面大口径天文望远镜主镜采用主动支撑的结构形式,通过主动调整主镜背部应力的分布实现对主镜曲率半径及面形调整的目的。然而,空间望远镜受到负载量、稳定性、安全性等诸多因素的限制,导致空间望远镜主镜不宜采用主动支撑的结构形式。因此,有必要提出一种新的方法,以最小的代价实现对主镜曲率半径误差引起的波前像差的变化进行补偿,同时,满足在轨像差补偿高可靠性的技术方案要求。
发明内容
本发明为解决现有大口径离轴三反空间望远镜发射入轨后,主镜曲率半径改变量超出公差范围时导致波前像差不满足天文观测要求的问题,提供一种在轨补偿主镜曲率半径误差影响的调整方法。
一种在轨补偿主镜曲率半径误差影响的调整方法,该方法由以下步骤实现:
步骤一、空间望远镜在轨对宇宙中恒星进行成像获得星点图像,通过星点图像灰度值,获得在轨状态下系统波前像差;
步骤二、判断所述在轨状态下系统波前像差是否满足天文观测要求,如果是,则进行天文观测;如果否,则执行步骤三;
步骤四、根据主镜曲率半径误差与像散增量和彗差增量的解析关系,计算主镜曲率半径误差;解析关系式如下:
式中,Δr为主镜曲率半径误差,为归一化视场矢量,分别为归一化视场里的x轴坐标和y轴坐标,为归一化光瞳矢量,分别为x轴和y轴光瞳偏心矢量,K040r、K222r和K131r分别为波前像差系数W040、W222、W131关于曲率半径的敏感度;
步骤五、根据波前像差RMS的最小值,确定次镜调整量;具体为:
将次镜调整量作为变量,建立系统像散增量、彗差增量与主镜曲率半径误差和次镜轴向调整距离的解析关系,分别由下式表示为:
求取系统波前像差RMS最小值,即:
步骤六、在轨调整次镜及摆镜,调整完毕后,空间望远镜进入工作状态,返回步骤一。
本发明的有益效果:
针对波前像差方程而言,主镜曲率误差主要改变波前像差系数Wklm,而轴向调整次镜位置同样改变波前像差系数Wklm。同时,主镜曲率半径误差和调整次镜轴向位置所引入的离焦像差可忽略,因其可通过摆镜进行校正。故根据主镜曲率半径误差和调整次镜轴向位置所引入的与视场相关的像散改变量和与视场相关的彗差改变量对波前像差RMS的综合贡献量的极小值来确定次镜调整量。
本发明通过改变次镜轴向位置(主动引入波前像差)来补偿主镜曲率半径误差影响的在轨调整方法,主要具备以下优点:
(1)基于解析理论及解析算法,有利于保证在轨实施的可靠性;
(2)根据像差理论,分析主镜曲率半径误差对波前像差影响的内在机理,有针对性的提出补偿方法。
(3)补偿方法易于实现。计算简单,调整方便,提高了在轨补偿的可靠性。
附图说明
图1为本发明所述的一种在轨补偿主镜曲率半径误差影响的调整方法的流程图;
图2中(a)为系统设计状态下的像散分布图,(b)为主镜存在曲率半径误差状态下的像散分布图;
图3中(a)为系统设计状态下的彗差分布图,(b)为主镜存在曲率半径误差状态下的彗差分布图;
图4中(a)为存在曲率半径误差(-0.2mm)的像散分布图,(b)为调整次镜补偿后(-0.11mm)的像散分布图。
图5中(a)为存在曲率半径误差的彗差分布图,(b)为调整次镜补偿后的彗差分布图;
图6中(a)为系统设计状态下波前像差RMS的分布图,(b)为补偿主镜曲率半径误差后,系统波前像差RMS的分布图。
具体实施方式
具体实施方式一、结合图1说明本实施方式,一种在轨补偿主镜曲率半径误差影响的调整方法,该方法由以下步骤实现:
步骤a:空间望远镜指向事先规划好的满足主动光学技术要求的特定天区(该天区可以保证特定视场点存在亮度足够的恒星),即满足主动光学波前检测需求。望远镜在轨对宇宙中恒星进行成像获得星点图像,通过星点图像灰度值(如相位差算法、相位恢复算法等)计算在轨状态下系统波前像差;
步骤b:判断在轨状态下系统波前像差(含有主镜曲率半径误差引入的波前像差)是否满足天文观测要求,满足则进行天文观测;不满足,则将当前系统波前像差值(Zernike系数值)作为在轨补偿主镜曲率半径误差影响的调整方法的输入;
步骤c:将输入的波前像差(Zernike系数值)与不含主镜曲率半径误差状态下的(地面装调状态下的)波前像差(Zernike系数值)做差,获得主镜曲率半径误差引入的像散增量、彗差增量(Hx、Hy分别为x轴和y轴的归一化视场坐标;
步骤d:根据主镜曲率半径误差与像散、彗差增量的解析关系,计算主镜曲率半径误差值,解析关系式如式(1)所示:
式中,Δr为主镜曲率半径误差,为归一化视场矢量,分别为归一化视场里的x轴坐标和y轴坐标,为归一化光瞳矢量,分别为x轴和y轴光瞳偏心矢量,K040r、K222r和K131r分别为波前像差系数W040、W222、W131关于曲率半径的敏感度,其可通过光学设计软件进行建模计算。
步骤e:将次镜调整量作为变量,建立系统像散增量、彗差增量与主镜曲率半径误差值、次镜轴向调整距离的解析关系,如下式:
步骤g:在轨调整次镜及摆镜(根据离焦量调节,采用现有技术实现),调整完毕后,空间望远镜进入工作状态,返回步骤a。
具体实施方式二、结合图1至图6说明本实施方式,本实施方式为具体实施方式一所述的一种在轨补偿主镜曲率半径误差影响的调整方法的实施例:
1、解算主镜曲率半径误差:
根据矢量像差理论,由曲率半径误差引起的像散增量可表示为:
其中,Δr为主镜曲率半径误差,为归一化视场矢量,为归一化光瞳矢量,为光瞳偏心矢量,分别为x轴和y轴光瞳偏心矢量,K040r、K222r、K131r分别为波像差系数W040、W222、W131关于曲率半径的敏感度,可通过光学设计软件进行建模计算。
以类WFIRST空间天文望远镜离轴方案为例,对主镜曲率半径误差对像散分布的影响进行说明,图2中(a)为系统设计状态下的像散分布,图2(b)为主镜存在曲率半径误差状态下的像散分布。
根据矢量像差理论,由曲率半径误差引起的彗差增量可表示为:
本实施方式中结合图3对主镜曲率半径误差对彗差分布的影响进行说明,图3中(a)为系统设计状态下的像散分布,(b)为主镜存在曲率半径误差状态下的彗差分布。
空间望远镜在轨运行时,空间望远镜对星点进行成像,在望远镜系统焦面上将获得星点的像,主动光学技术中的波前解算技术,通过相位差算法、相位恢复算法等,将得到系统当前状态下特定视场的波前像差(条纹Zernike系数值),进而可以确定当前状态下不同视场点像散、彗差与理想状态(地面最佳装调状态)下像散、彗差差异,即像散增量和彗差增量联合(7)式和(8)式,将相对准确地计算出主镜曲率半径误差。
2、解算次镜轴向位姿在轨调整量
大口径离轴三反空间望远镜在轨调整不比于地面装调,通常主镜、次镜、三镜以及焦面不全具有位姿调整能力。一方面,各个元件都配备调整模块,由于发射、运输及在轨重力卸载等的存在,将增大元件相对位姿失调的风险。另一方面,增加调整模块,将相应增加电子学模块,这也将给系统的功耗、热控以及安全性带来巨大的挑战。因此,很多正在或计划在轨应用的空间望远镜只将次镜和摆镜作为调整元件(其中次镜主要用来校正或补偿像差,摆镜主要用来调整后工作距),而主镜、三镜及像面在轨时处于相对固定的状态。
调整次镜对波前像差的影响机理与主镜曲率半径误差对波前像差的影响机理具有相似性,当曲率半径误差与次镜轴向位移偏差(主动调整次镜轴向位置)同时存在时,像散增量、彗差增量(差异量)可表示为:像散、彗差与主镜曲率半径误差值、次镜轴向调整距离的解析关系式为:
主镜曲率半径误差和次镜轴向位姿改变主要影响像散、彗差和离焦,大口径空间望远镜通常配有摆镜,并用其校正离焦像差(不是本实施例关注的点),故在轨补偿曲率误差时,只需要考虑像散和彗差对系统波前像差的贡献。系统波前像差RMS与各类型波前像差具有如下关系:
根据数学极值思想,为了获得极小值,就需要建立ΔRMSsystem关于Δd1d2的导数,当导数值为0时,系统波前像差最小,即可求取次镜的调整距离:
将通过第1步解算出来的主镜曲率半径误差带入(6)时,并通过求解上式,便可计算出次镜调整距离。图4中(a)为存在曲率半径误差(-0.2mm)的像散分布图,(b)为调整次镜补偿后(-0.11mm)的像散分布图。图5中(a)为存在曲率半径误差的彗差分布图,b)为调整次镜补偿后的彗差分布图。通过图6可以看出,补偿后,系统波前像差RMS分布与设计状态基本一致,成像质量满足天文观测需求。
Claims (2)
1.一种在轨补偿主镜曲率半径误差影响的调整方法,其特征是:该方法由以下步骤实现:
步骤一、空间望远镜在轨对宇宙中恒星进行成像获得星点图像,通过星点图像灰度值,获得在轨状态下系统波前像差;
步骤二、判断所述在轨状态下系统波前像差是否满足天文观测要求,如果是,则进行天文观测;如果否,则执行步骤三;
步骤四、根据主镜曲率半径误差与像散增量和彗差增量的解析关系,计算主镜曲率误差;解析关系式如下:
式中,Δr为主镜曲率半径误差,为归一化视场矢量,分别为归一化视场里的x轴坐标和y轴坐标,为归一化光瞳矢量,分别为x轴和y轴光瞳偏心矢量,K040r、K222r和K131r分别为波前像差系数W040、W222、W131关于曲率半径的敏感度;
步骤五、根据波前像差RMS的最小值,确定次镜调整量;具体为:
将次镜调整量作为变量,建立系统波前像差RMS与主镜曲率半径误差和次镜轴向调整距离的解析关系,分别由下式表示为:
求取系统波前像差RMS最小值,即:
步骤六、在轨调整次镜及摆镜,调整完毕后,空间望远镜进入工作状态,返回步骤一。
2.根据权利要求1所述的一种在轨补偿主镜曲率半径误差影响的调整方法,其特征在于:步骤一中,空间望远镜在轨成像前,需要满足主动光学波前检测的要求,即:空间望远镜指向事先规划好的满足主动光学技术要求的天区,该天区保证相应视场点存在亮度足够的恒星。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110274445.9A CN113031201B (zh) | 2021-03-15 | 2021-03-15 | 一种在轨补偿主镜曲率半径误差影响的调整方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110274445.9A CN113031201B (zh) | 2021-03-15 | 2021-03-15 | 一种在轨补偿主镜曲率半径误差影响的调整方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113031201A true CN113031201A (zh) | 2021-06-25 |
CN113031201B CN113031201B (zh) | 2022-02-11 |
Family
ID=76468790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110274445.9A Active CN113031201B (zh) | 2021-03-15 | 2021-03-15 | 一种在轨补偿主镜曲率半径误差影响的调整方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113031201B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114994859A (zh) * | 2022-05-06 | 2022-09-02 | 北京理工大学 | 基于cnn的tma望远镜面形、位姿误差在轨校正方法 |
CN115586640A (zh) * | 2022-10-21 | 2023-01-10 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
CN115641376A (zh) * | 2022-10-17 | 2023-01-24 | 中国科学院长春光学精密机械与物理研究所 | 望远镜在轨位姿失调检测方法、装置、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0909971A1 (fr) * | 1996-07-26 | 1999-04-21 | Giat Industries | Procédé pour corriger la surface d'onde d'un faisceau optique au moyen d'un miroir déformable et mirroir pour la mise en oeuvre du procédé |
US6426834B1 (en) * | 1999-02-13 | 2002-07-30 | Leica Geosystems Ag | Method for correcting optical wavefront errors and optical system, such as a telescope, produced accordingly |
CN101078808A (zh) * | 2007-06-28 | 2007-11-28 | 中国科学院光电技术研究所 | 基于像清晰化原理的自适应光学星体目标成像系统 |
JP2015004644A (ja) * | 2013-06-24 | 2015-01-08 | キヤノン株式会社 | 光学系の偏芯量算出方法及びそれを用いた光学系の調整方法 |
US20170168399A1 (en) * | 2014-09-15 | 2017-06-15 | Carl Zeiss Smt Gmbh | Projection lens, projection exposure apparatus and projection exposure method for euv microlithography |
-
2021
- 2021-03-15 CN CN202110274445.9A patent/CN113031201B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0909971A1 (fr) * | 1996-07-26 | 1999-04-21 | Giat Industries | Procédé pour corriger la surface d'onde d'un faisceau optique au moyen d'un miroir déformable et mirroir pour la mise en oeuvre du procédé |
US6426834B1 (en) * | 1999-02-13 | 2002-07-30 | Leica Geosystems Ag | Method for correcting optical wavefront errors and optical system, such as a telescope, produced accordingly |
CN101078808A (zh) * | 2007-06-28 | 2007-11-28 | 中国科学院光电技术研究所 | 基于像清晰化原理的自适应光学星体目标成像系统 |
JP2015004644A (ja) * | 2013-06-24 | 2015-01-08 | キヤノン株式会社 | 光学系の偏芯量算出方法及びそれを用いた光学系の調整方法 |
US20170168399A1 (en) * | 2014-09-15 | 2017-06-15 | Carl Zeiss Smt Gmbh | Projection lens, projection exposure apparatus and projection exposure method for euv microlithography |
Non-Patent Citations (2)
Title |
---|
HAMID HEMMATI 等: "Telescope Wavefront Aberration Compensation with a Deformable Mirror in an Adaptive Optics System", 《PROCEEDINGS OF SPIE》 * |
赵海波 等: "大型空间望远镜次镜调整机构精度分析与测试", 《光学精密工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114994859A (zh) * | 2022-05-06 | 2022-09-02 | 北京理工大学 | 基于cnn的tma望远镜面形、位姿误差在轨校正方法 |
CN115641376A (zh) * | 2022-10-17 | 2023-01-24 | 中国科学院长春光学精密机械与物理研究所 | 望远镜在轨位姿失调检测方法、装置、设备及介质 |
CN115586640A (zh) * | 2022-10-21 | 2023-01-10 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
CN115586640B (zh) * | 2022-10-21 | 2024-07-09 | 中国人民解放军国防科技大学 | 一种同轴两反系统分步计算的装调方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113031201B (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113031201B (zh) | 一种在轨补偿主镜曲率半径误差影响的调整方法 | |
EP2257847B1 (en) | Low order adaptive optics by translating secondary mirror of off-aperture telescope | |
US4471448A (en) | Method and apparatus for aligning an optical system | |
CN113066021B (zh) | 基于图像特征匹配的空间望远镜在轨像差补偿方法 | |
Diolaiti et al. | Conceptual design and performance of the multiconjugate adaptive optics module for the European Extremely Large Telescope | |
Pariani et al. | MORFEO optical design and performances: status at preliminary design review | |
US10241326B2 (en) | Correcting undesired distortions or aberrations and generating desired wavefronts in optical imaging, sensing, signaling and other applications based on bi-valued walsh functions | |
CN109522573A (zh) | 一种光学遥感相机主动光学系统的仿真方法 | |
Rodeghiero et al. | The impact of ELT distortions and instabilities on future astrometric observations | |
WO2010110934A1 (en) | Real-time optical compensation of orbit-induced distortion effects in long integration time imagers | |
CN109343212B (zh) | 大视场巡天望远镜主动光学系统实现方法 | |
CN117268259A (zh) | 一种基于cnn的大动态范围次镜位置误差检测方法 | |
EP3644107B1 (en) | Adaptive optical system and adaptive optical method | |
CN115682991A (zh) | 在轨检测拼接式望远镜各子镜曲率半径非一致性的方法 | |
CN114994859B (zh) | 基于cnn的tma望远镜面形、位姿误差在轨校正方法 | |
Upton et al. | Active reconstruction and alignment strategies for the Advanced Technology Solar Telescope | |
CN113849765B (zh) | 以全视场椭率最小为校正目标的空间望远镜在轨调校方法 | |
Sahoo et al. | Segment-level thermal sensitivity analysis for exo-Earth imaging | |
Yang et al. | An adaptive laser beam shaping technique based on a genetic algorithm | |
CN117130172B (zh) | 一种拼接式空间望远镜全局装调方法、装置、设备及介质 | |
Diolaiti et al. | The E-ELT multi conjugate adaptive optics module | |
Dai et al. | Active correction experiment on a 1.2 m thin primary mirror | |
Li et al. | Experimental Study on Wavefront Distortion Correction of 40-Elements Adaptive Optical System Based on SPGD Algorithm | |
Li et al. | Testing of the stochastic parallel radient descent algorithm for the alignment of a two-mirror telescope | |
CN116974063A (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 |