1、第45卷 第8期 红 外 技 术 Vol.45 No.8 2023年8月 Infrared Technology Aug.2023 828 红外热成像折转光学系统的光轴静态敏感度分析 杨 丹1，金 宁1，杨开宇1，李 晶1，董树林1，胡健钏1，李晓君2（1.昆明物理研究所，云南 昆明 650223；2.空装成都局驻昆明地区军代表室，云南 昆明 650223）摘要：红外热成像折转光学系统在复杂环境条件下，光轴容易因为光学元件的偏心或倾斜而发生漂移，影响系统对目标的指示精度。在红外热成像系统设计之初对光学系统开展光轴静态敏感度分析，能够识别出光学系统的敏感点，为满足光轴稳定性的结构优化设计提供约束
2、条件。通过基于旋转矩阵的坐标变换，建立了光学元件旋转过程量和倾斜状态量的转换关系，从而实现了光学元件在任意方向倾斜的空间姿态模拟，确保了光轴敏感度公差分析中的蒙特卡罗采样与结构设计的约束条件相对应，并在此基础上搭建了对红外折转光学系统光轴静态敏感度分析的流程，编制了程序。用所编程序对某典型红外热成像折转光学系统进行实例分析，根据光轴稳定性的指标要求，依次对光学系统中各光学件的偏心量和倾斜量进行了光轴的灵敏度和反灵敏度分析，得出了初始公差限，再针对初始公差限数据进行了任意方向采样的蒙特卡罗分析，最终得出了各光学元件能够满足光轴稳定性指标的偏心和倾斜公差限数据，通过建立多重坐标系的方法验证了所得数
3、据的准确性，为指导光机热优化设计奠定了基础。关键词：红外折转光学系统；光轴漂移；静态敏感度；公差分析；蒙特卡罗采样 中图分类号：O435 文献标识码：A 文章编号：1001-8891(2023)08-0828-09 Optical Axis Static Sensitivity Analysis for Infrared Thermal Imaging Folding Optical System YANG Dan1，JIN Ning1，YANG Kaiyu1，LI Jing1，DONG Shulin1，HU Jianchuan1，LI Xiaojun2(1.Kunming Institute
4、 of Physics,Kunming 650223,China;2.Military Representative Office of Chengdu Bureau of AIR Force Stationed in Kunming Area,Kunming 650223,China)Abstract:The optical axis of an infrared thermal imaging folding optical system is prone to shift owing to decentering of the tilt of optical components und
5、er complex environmental conditions,which affects the indication accuracy of the system for the target.Static sensitivity analysis of the optical axis for the optical system at the beginning of the design of the infrared thermal imaging system is useful for identifying the sensitive points of the op
6、tical system and provides constraints for the structural optimization design to meet the stability of the optical axis.The conversion relationship between the rotation process and spatial state quantities of the optical components was established by coordinate transformation based on the rotation ma
7、trix to simulate the spatial attitude of the optical component tilted in any direction and ensure that the Monte Carlo sampling in the optical axis sensitivity analysis corresponds to the constraint conditions of the structural design.On this basis,the flow of the static sensitivity analysis of the
8、optical axis of the infrared folding optical system was established,and a program was compiled.A typical infrared thermal imaging folding optical system was analyzed using this program.According to the index requirements of the optical axis stability,the optical axis sensitivity and inverse sensitiv
9、ity of the decenter and tilt of each optical component in the optical system were analyzed,and the initial tolerance limit was obtained.Then,Monte Carlo analysis sampling could be performed in any direction according to the initial tolerance limit data;thus,the decenter and tilt tolerance limit data
10、 that meet the optical axis stability index could be obtained,and the accuracy of the obtained data was verified by establishing a multi-coordinate system.Static sensitivity 收稿日期：2022-12-27；修订日期：2023-02-15.作者简介：杨丹（1991-），女，云南昆明人，工程师，硕士，主要从事光学系统设计和仿真方面的研究。E-mail：。通信作者：杨开宇（1984-），女，云南昆明人，正高级工程师，博士，主要从
11、事光学系统设计和仿真方面的研究。E-mail：。基金项目：国家重点研发计划（2017YFA0701200）。第45卷 第8期 Vol.45 No.8 2023 年 8 月 杨 丹等：红外热成像折转光学系统的光轴静态敏感度分析 Aug.2023 829 analysis provides a foundation for guiding the design of opticalmechanical thermal optimization.Key words:infrared folding optical system,optical axis shift,static sensitivit
12、y,tolerance analysis,Monte Carlo sampling 0 引言 红外热成像系统工作在较为严酷的使用环境中时，受安装力、高低温和冲击振动等因素影响，内部的机械支撑结构产生微小的变形，使光学元件出现空间状态的偏移（包含偏心或倾斜），导致系统的光轴发生漂移，影响系统对目标的指示精度1-6。红外折转光学系统具有体积小、布局紧凑的优势，但在工程实践中发现，相较于直线形光学系统，光轴在环境影响下的漂移量相对更大。为了满足使用要求，光轴的最大漂移角度即光轴稳定性必须限定在允许范围之内。因此，在系统设计阶段，需严格分析各光学元件的偏移误差对系统光轴带来的影响，识别出光学系统在静态
13、下的敏感点，同时基于光轴稳定性指标分析出各个镜片偏心和倾斜的合理公差限，从而能够得到光学系统所容许各支撑结构件最大的形变范围，以作为结构优化的约束条件7-14，为光轴稳定性的设计提供指导。对于折转光学系统，由于其非旋转对称性，光学元件的偏移状态需要考虑偏移量和偏移方向两方面信息。也就是说，在分析某个元件允许的最大偏心量或倾斜量时，需要对所有可能的偏心方向或倾斜方向进行采样，选取最严苛的情况作为结构优化的约束。国内外有关光学系统公差灵敏度的研究，在考虑元件的偏心量或倾斜量时大多仅针对 X-轴方向或 Y-轴方向进行采样，而极少考虑其他倾斜方向，对折转光学系统的适用性较低15-17。目前主流的光学软
14、件具有完备分析公差的流程，但其对光学元件倾斜状态的建模是先后绕 X、Y、Z 各个坐标轴的旋转形成，其旋转值为过程量，不能代表旋转后形成的空间状态量18-19。利用现有软件进行公差分析时，只能完全按照建模的输入方式，对光学元件旋转过程量采样，无法对其倾斜的空间状态量采样；而对结构优化来说，所需的约束条件为状态量，这就导致采样结果与约束条件不能对应，对结构设计的指导存在局限性。本文基于公差分析的思想，立足于红外折转光学系统元件的偏心、倾斜静态公差，建立光学元件旋转过程量与旋转形成的空间状态量（即倾斜量和倾斜方向）之间的转换关系，实现对光学元件在任意方向倾斜状态的准确采样；在此基础上，建立光学元件偏
15、移对光轴稳定性的影响的灵敏度分析、反灵敏度分析及蒙特卡罗分析整个公差计算流程，以获得折转光学系统的静态敏感点、以及结构优化设计所需的约束条件。1 光学元件旋转过程量和倾斜状态量的转换关系建立 蒙特卡罗公差分析是光轴静态敏感度分析中最为重要的一步，是最终获得光学系统中各镜片偏心和倾斜合理公差限的必要过程。在红外折转光学系统光轴敏感度公差的蒙特卡罗分析中，为了和结构优化设计能够对应，需在每个元件垂直于光轴的参考平面的360所有方向上，对倾斜量进行蒙特卡罗采样，然后根据采样的状态构建出对应的倾斜模型，从而能够对该模型进行非顺序光线追迹，计算出光轴的最终漂移情况。对倾斜模型的构建，需要通过旋转元件来完
16、成。在几何空间中，每个元件自带一个右手局部坐标系O-XYZ，其中 Z-轴沿光轴方向。元件的旋转运动姿态采用欧拉角来表述，旋转后的最终姿态除了与角度相关，还与旋转顺序相关20。在对元件的倾斜进行建模时，假设旋转按图 1 所规定的顺序，即元件先绕自身坐标轴的 X-轴左旋角度，再绕旋转后形成的坐标轴的 Y-轴左旋角度，再绕旋转后形成的坐标轴的 Z-轴右旋角度，得到元件的新坐标系 O-X1Y1Z1，则空间中某点 W 在原坐标系 O-XYZ 和新坐标系 O-X1Y1Z1中的坐标描述对应关系为：111100cos0sin0cossin0100sincossin0coscossin0sincos0001XY
17、ZXYZ-(1)式中：X,Y,Z是点 W 在原坐标系 O-XYZ 中的坐标，X1,Y1,Z1是点 W 在旋转后坐标系 O-X1Y1Z1中的坐标。第45卷 第8期 红 外 技 术 Vol.45 No.8 2023年8月 Infrared Technology Aug.2023 830 图 1 光学元件的旋转规则 Fig.1 Rotation rules of optical component 旋转后，新 Z-轴的方向向量为1n，将O-X1Y1Z1中的Z-轴的单位向量0,0,1代入上式中X1,Y1,Z1，可求出经过旋转后，1n在原始坐标系O-XYZ下的向量坐标X,Y,Z：sincossincos
18、cosXYZ (2)由于Z-轴为元件的光轴方向，因此绕Z-轴的旋转并不影响方向向量最后的朝向，所以由式(2)可以看到，方程组是一个与无关的等式。此时，若要描述元件旋转后的倾斜姿态，可使用空间状态量，即描述倾斜方向的u角度，以及描述倾斜量值的v角度，如图2所示，其中，倾斜角v代表光学元件倾斜前后的夹角，方向角u代表元件倾斜的参考轴线与原坐标系Y-轴的夹角。经过旋转后的1n在原始坐标系O-XYZ下的向量坐标X,Y,Z也可以由状态量u、v表示：cossinsinsincosXuvYuvZv (3)图 2 旋转后的状态量 Fig.2 Status parameters after rotation 如
19、果要实现任意方向倾斜采样的蒙特卡罗分析，需要建立光学元件旋转过程，和空间状态量的u，v之间的转换关系。因此联立式(2)与式(3)，得到了u，v与，之间的对应关系：arcsin(cossin)sinsinarcsin()cosuvuvb (4)这样，通过旋转矩阵的推导，解出了光学元件倾斜状态量（方向角u、倾斜角v）到倾斜建模过程量（旋转角，）的转换对应公式。在任意方向蒙特卡罗分析中，只需对方向角u和倾斜角v进行采样，就可利用上述公式将采样数据转换建模需要的旋转角，从而构建出元件倾斜的模型，进行非顺序光线追迹。2 光轴静态敏感度分析流程构建 2.1 基本思路 对红外折转光学系统光轴静态敏感度分析，
20、本质上是公差分析，需将光学元件的偏心、倾斜公差引入光学系统中，分析由元件偏移引起的光轴漂移与光轴稳定性指标的匹配情况。构建红外折转光学系统的光轴敏感度公差分析流程的基本思路如图3所示，输入为光学系统的仿真模型以及系统的光轴稳定性指标，经过光轴敏感度分析以后，最终输出的是光学系统中各个元件合理的公差限数据。该思路对光轴敏感度分析分解为三步：灵敏度分析-反灵敏度分析-蒙特卡罗分析21，每一步的分析需要基于上一步分析的结果数据作为输入，而输出的数据将作为下一步分析的输入。流程中首先开展灵敏度公差分析，根据光轴稳定性指标进行初步分解，将光学系统的公差大致设定在某个范围，用非顺序建模实现对光学元件空间姿
21、态的精确模拟，并通过非顺序光线追迹，仿真计算这些公差所造成的光轴漂移量数据（如图4所示），此数据即为各元件偏心和倾斜的光轴灵敏度数据，以此识别敏感元件。第45卷 第8期 Vol.45 No.8 2023 年 8 月 杨 丹等：红外热成像折转光学系统的光轴静态敏感度分析 Aug.2023 831 图 3 构建红外光学系统光轴敏感度公差分析流程的基本思路 Fig.3 Basic idea of constructing tolerance analysis flow of optical axis sensitivity for infrared optical system (a)(b)图 4
22、透镜倾斜(a)和透镜偏心(b)引起的光轴漂移 Fig.4 Optical axis shift caused by tilted lens(a)and decentered lens(b)反灵敏度分析则假设光学系统的评价标准（即光轴漂移量）变化特定的数值，评估造成相应变化的各项公差范围，称之为光轴反灵敏度数据，可根据上一步光轴灵敏度数据的多次计算迭代求得。光轴的反灵敏度数据，可作为蒙特卡罗公差分析中的输入，即各光学件偏心和倾斜的光差限初值，以公差限初值作为反射镜各运动分量的蒙特卡洛采样区间，在区间内进行各元件偏心和倾斜随机取值并叠加成组合误差，进行了蒙特卡洛分析，得出光轴漂移量的初步统计评估结
23、果，与系统的光轴稳定性的指标要求对比后，放宽或收紧采样区间，进行多轮分析，最终得到合理的公差限。2.2 程序设计 根据上述思路，对光学系统进行光轴敏感度分析的程序流程如图5所示，基于Matlab编制计算程序，通过调用LightTools软件实现光学元件偏移及非顺序光线追迹，从而完成光轴漂移仿真数据计算。该流程以光学系统的光轴稳定性指标Iaxis作为流程输入，以光学系统中各镜片的偏心和倾斜合理公差限为输出。程序中，调取各透镜的原始偏心值D0和原始倾斜值T0，分别叠加标准偏心值Dstan和标准倾斜值Tstan，通过对轴上主光线追迹计算光轴漂移，并进行公差分析，得到光学元件的灵敏度数据。基于光学元件
24、的灵敏度数据和系统的光轴稳定性指标Iaxis，计算出光学元件的反灵敏度数据。基于光学元件的光轴反灵敏度数据，得到公差限初值0，即光学元件各运动分量的蒙特卡洛采样区间。给定循环次数Kloop，在公差限初值范围内对每一光学件的偏心方向u，偏心量v，旋转方向m，旋转量n四个值分别产生Kloop组随机数。偏心方向u和旋转方向v在0360的范围内采样，偏心量v和旋转量n则在公差限初值的范围内采样。由偏心方向u，偏心量v，旋转方向m，旋转量n四个值，通过上节中的数据变换，计算得出光学元件建模所需的过程量输入，即沿X-轴方向的平移量x，沿Y方向的平移量y，绕X-轴的旋转量，绕X-轴的旋转量。m，n与x，y的
25、对应关系可由简单的三角函数求得：cossinxnmynm (5)式中：u，v与，的对应关系由上节中推导得出的光学元件旋转过程量和倾斜状态量的转换关系式(4)求得。第45卷 第8期 红 外 技 术 Vol.45 No.8 2023年8月 Infrared Technology Aug.2023 832 图 5 光轴稳定性公差分析流程 Fig.5 Flow diagram for tolerance analysis of stability of optical axis 基于计算得出的位置与旋转角度输入，进行光学元件偏移状态的建模和光轴漂移的仿真分析，得到各光学元件偏心与倾斜随机采样数据搭配后
26、的光线输出数据，对输出数据进行统计分析，得出光轴漂移量的蒙特卡洛分析结果。将输出的蒙特卡罗分析结果与光轴漂移量的指标要求进行对比，得出新一轮蒙特卡罗分析的公差限调整系数Pn+1，如下式计算：Pn+1Iaxis/In (6)式中：In为本轮光轴漂移量蒙特卡洛分析结果。若Pn+1已接近1，则说明输入的公差限已经合理，无需再调整。否则，则需按照下式放宽或收紧公差限：n+1Pn+1n (7)式中：n为本轮的公差限；n+1为新一轮公差限。获取新一轮公差限后，再次展开蒙特卡罗分析，重复迭代此过程，直至公差限调整系数P接近1，即光轴漂移量的蒙特卡罗统计结果与系统光轴稳定性指标要求接近，则可将该轮的公差限作为
27、光学元件偏移的约束值。3 光轴敏感度分析实例 根据上述搭建的光轴敏感度分析流程和计算程序，选取了光轴较为敏感的某型红外热成像光学系统作为分析对象，光学系统的仿真模型如图6所示。该系统为制冷型中波热像仪，为减少体积，采用了两个折转反射镜形成U型折转光路。本例以系统的光轴稳定性指标Iaxis0.25 mrad为例，对系统中8个光学透镜和2个折转反射镜的倾斜和偏心的公差限进行分析计算。图 6 光学系统的静态仿真模型 Fig.6 Static simulation model of the optical system 在倾斜和偏心的灵敏度分析中，本例的标准偏心值Dstan设置为0.01 mm，标准倾
28、斜值Tstan设置为0.01。在倾斜和偏心的反灵敏度分析中，由于系统的光轴稳定性指标Iaxis0.25 mrad，而光学系统中包括反射镜在内，共有10个光学件，且每个光学件包含两个公差（偏心和倾斜），这样一来，共由20个公差对光轴漂移产生贡献，因此，将0.25 mrad的指标均分成20份，以光轴漂移量为0.0125 mrad作为反灵敏度分析的标准。图7(a)为倾斜灵敏度分析结果，横坐标是光学元件的序号，纵坐标是光轴漂移量；图7(b)为倾斜反灵敏度分析结果，横坐标是光学元件的序号，纵坐标是光学件的倾斜量。从倾斜灵敏度折线图中可以较为直观地看出，当每个透镜同时倾斜0.01，对光轴漂移量产生的贡献存
29、在很大差异。对比倾斜灵敏度折线图与倾斜反灵敏度折线图中可以发现，倾斜反灵敏度的数据与倾斜灵敏度数据正好相反，第2反射镜的倾斜灵敏度最高，倾斜反灵敏度最低，第5透镜的倾斜灵敏度最低，倾斜反灵敏度最高。从图中能够识别出，第1反射镜和第2反射镜是对倾斜最敏感的关键元件。第45卷 第8期 Vol.45 No.8 2023 年 8 月 杨 丹等：红外热成像折转光学系统的光轴静态敏感度分析 Aug.2023 833 图8(a)为偏心灵敏度分析结果，横坐标是光学件的序号，纵坐标是光轴漂移量；图8(b)为偏心反灵敏度分析结果，横坐标是光学件的序号，纵坐标是光学件的偏心量。对比偏心灵敏度折线图与偏心反灵敏度折线
30、图中可以发现，偏心反灵敏度的数据与偏心灵敏度数据同样正好相反，第2透镜偏心灵敏度最高，偏心反灵敏度最低，第5透镜的偏心灵敏度最低，偏心反灵敏度最高。从图中能够识别出，第2透镜和第3透镜是对偏心最敏感的关键元件。(a)(b)图 7 倾斜灵敏度(a)和倾斜反灵敏度(b)折线图 Fig.7 Line chart of tilt sensitivity(a)and tilt anti-sensitivity(b)(a)(b)图 8 偏心灵敏度(a)和偏心反灵敏度(b)折线图 Fig.8 Line chart of decenter sensitivity(a)and decenter anti-sen
31、sitivity(b)将光轴反灵敏度分析数据设置为各光学元件偏心和倾斜的初始公差限，在8个透镜与2个反射镜的偏心与倾斜初始公差限范围内进行蒙特卡罗分析，本例中循环次数Kloop设置为10000次。图9(a)为蒙特卡罗分析后输出的10000组光学系统光轴漂移量数据的分布区间直方图，横坐标为光轴漂移量，纵坐标为分布在相应光轴漂移量区间的蒙特卡罗实验次数，由图可知，分布在光轴漂移量为0.021 mrad附近区间的实验次数最多。图9(b)为数据的累积分布函数图，横坐标为光轴漂移量，纵坐标为累计分布概率，可以看到，50%的概率光轴移动量在0.025 mrad范围内，70%的概率光轴移动量在0.033 m
32、rad范围内，90%的概率光轴移动量在0.0456 mrad范围内。本文中以90%的概率为例，因此基于初始公差限，光学系统的光轴移动量评估结果大概率分布在0.0456 mrad范围内，而对比本例设置的指标，系统光轴稳定性Iaxis0.25 mrad。因此基于蒙特卡罗分析评估的光轴漂移量结果还远低于光轴一致性指标，说明设置的初始公差限过于保守，因此按照式(6)计算，公差调整系数等于5.48，并按照式(7)放宽各光学件的偏心与倾斜公差限。根据调整后的公差限，进行蒙特卡罗随机采样。新一轮蒙特卡罗分析的输出数据分布区间直方图及累积分布函数图如图10(a)和图10(b)所示。本轮输出的光轴漂移量数据以9
33、0%的概率分布在00.256 mrad范围内，基本与光学系统的光轴一致性指标相匹配。因此将本轮输入的公差限作为最终各光学件的倾斜和偏心合理公差。表1为采用本文程序对目标光学系统依次进行了灵敏度分析、反灵敏度分析、蒙特卡罗分析之后，得出的各个光学元件的合理公差限。此公差限可作为结构优化设计的约束条件，只要由结构件引起的光学元第45卷 第8期 红 外 技 术 Vol.45 No.8 2023年8月 Infrared Technology Aug.2023 834 件偏心和倾斜不超过表1数据的范围，即可保证光学系统的光轴漂移量不超出系统光轴稳定性0.25 mrad的指标要求。从表1数据和前面灵敏度、
34、反灵敏度分析结果可以看出，结构优化设计需要着重考虑高低温或外力的作用下，第1、2反射镜的倾斜、以及第2、3透镜的偏心。因此，该分析识别出了光学系统的敏感点，同时能够指导结构优化设计，以满足光轴稳定性指标。(a)(b)图 9 光轴漂移量分布直方图(a)和累积分布图(b)Fig.9 Distribution histogram(a)and cumulative distribution(b)of optical axis shift (a)(b)图 10 调整公差限后的光轴漂移量分布直方图(a)和光轴漂移量累积分布图(b)Fig.10 Distribution histogram(a)and cu
35、mulative distribution(b)of optical axis shift after adjusting tolerance 表 1 各光学件倾斜与偏心最终公差限 Table 1 Final tolerances of tilt and decenter of each optical component LENS Tilt tolerance/Decenter tolerance/mm Lens 1 0.17811389 0.01035527 Lens 2 0.19549745 0.00539823 Lens 3 0.15172749 0.00718520 Lens 4 1
36、.50604475 0.01545493 Lens 5 1.55067910 0.02021238 Mirror 1 0.02377837-Mirror 2 0.01895890-Lens 6 0.23078320 0.00937822 Lens 7 0.11583835 0.00953447 Lens 8 1.05775137 0.01556741 4 数据准确度验证 本文为实现任意角度倾斜的蒙特卡罗分析，采用数据转换的思路，利用坐标变换推导出了方向角和倾斜角到旋转建模角度的转换关系式。为了对文中角度转换公式的准确性进行验证，利用上一章实例中输出的公差限，并通过上文中提到的通过手动建立多重临
37、时坐标系的替代方法，对实例中光学系统再次进行蒙特卡洛分析，观察是否能得出一致的光轴移动量分布数据。在仿真软件中，如图11所示，根据每一个光学元件的原始坐标系分别建立一个相同的临时坐标系，并使临时坐标系绕Y-轴左旋90旋转，使临时坐标系第45卷 第8期 Vol.45 No.8 2023 年 8 月 杨 丹等：红外热成像折转光学系统的光轴静态敏感度分析 Aug.2023 835 的X-轴朝向原内层坐标系的Z-轴，临时坐标系的Y-轴朝向原内层坐标系的Y-轴，临时坐标系的Z-轴朝向原内层坐标系的X-轴负方向，将原始坐标系嵌套在临时坐标系的内层并绑定，让外层坐标系带动光学元件一同旋转，实现让光学件先绕光
38、轴旋转，在绕与光轴垂直方向的坐标轴旋转。以实例中基于最后一轮公差限产生随机采样数据，输入到临时坐标系中，进行蒙特卡洛分析。得到的数据分布区间直方图如图12(a)，与上一章中得出的数据分布区间直方图如图12(b)进行对比，可以看出，两种方法得出的结果分布完全一致，证明了数据转换方法的准确性。此外，相较于手动建立多重临时坐标系的方法，本文建立的方法适用性更强、效率更高，省去了很多繁琐过程，不需要为每一个新的光学系统的所有光学元件手动建立与其坐标相对应的两重临时坐标系，避免了数据输入或匹配错误。图 11 临时坐标系示意图 Fig.11 Sketch of temporary coordinate s
39、ystem (a)(b)图 12 数据转换方法得出(a)和临时坐标系方法得出(b)的分布直方图 Fig.12 Histogram of distribution obtained by data conversion method(a)and temporary coordinate system method(b)5 结论 在对红外热成像折转光学系统进行光轴稳定性的结构优化过程中，非常有必要针对光学系统开展光轴的静态敏感度分析，基于系统的光轴稳定性指标分析出各个光学元件偏心和倾斜的合理公差限，从而指导结构优化。为解决蒙特卡罗分析中对折转光学系统元件的任意方向倾斜采样的难点，通过基于旋转矩阵的
40、坐标变换，建立了光学元件旋转过程量，和空间状态量的u，v相互转换的方法，使蒙特卡罗分析能够模拟真实情况下的光学元件运动，计算效率高且具有通用性。在建立的数据转换方法的基础上，搭建了对红外光学系统光轴静态敏感度分析的流程，并运用Matlab对LightTools进行二次开发，实现了程序化。与此同时，以一个光轴较为敏感的实际光学系统为例，运用建立的方法对其开展了光轴静态敏感度分析，得出了光学系统的敏感点和各光学元件的合理公差限，并结合建立临时外重坐标系的方法，对数据准确性进行了验证，证实了公式推导以及程序的可靠性和准确性。文中所建立的光轴静态敏感度分析流程为实现完整的光机热优化设计奠定了基础，它虽
41、然是针对光学系统的光轴分析所搭建的，但与此同时，文中提出的任意方向倾斜的蒙特卡罗采样方法也为折转光学系统像质的公差分析提供了解决思路。后续将对如何优化文中方法分析出的敏感元件开展研究，以期进一步提高光学系统的光轴稳定性。参考文献：1 郭夏锐,廖志波,王春雨,等.光轴一致性误差对空间透射式系统像差和质量的影响J.红外与激光工程,2012,41(2):437-441.第45卷 第8期 红 外 技 术 Vol.45 No.8 2023年8月 Infrared Technology Aug.2023 836 GUO Xiarui,LIAO Zhibo,WANG Chunyu,et al.Aberrat
42、ion and alignment quality investigation of space transmission optical system with centering error of optical axesJ.Infrared and Laser Engineering,2012,41(2):437-441.2 马宏川,范宏波,林宇,等.热像仪光机热集成分析综述J.红外技术,2019,41(2):134-141.MA Hongchuan,FAN Hongbo,LIN Yu,et al.Review of thermal-structural-optical integrat
43、ion analysis of thermal imagerJ.Infrared Technology,2019,41(2):134-141.3 王阳,黄煜,李占峰,等.利用恒星对天文观测系统光轴平行性检校J.红外与激光工程,2017,46(5):0517003.WANG Yang,HUANG Yu,LI Zhanfeng,et al.Calibration of optical axis parallelism by using star for astronomical observation systemJ.Infrared and Laser Engineering,2017,46(5
44、):0517003.4 陈志斌,肖文健,马东玺,等.大间距多光轴一致性野外在线检测方法J.光学学报,2017,37(1):0112006.CHEN Zhibin,XIAO Wenjian,MA Dongxi,et al.A method for large distance multi-optical axis parallelism online detectionJ.Acta Optica Sinica,2017,37(1):0112006.5 杨耀山,赵泉,王建龙,等.红外光学系统光轴稳定性建模与分析J.电光与控制,2017,24(12):82-84.YANG Yaoshan,ZHAO
45、Quan,WANG Jianlong,et al.Optical axis stability modeling and analysis of infrared optical systemJ.Electronics Optics Control,2017,24(12):82-84.6 隋杰,程会艳,余成武,等.星敏感器光轴热稳定性仿真分析方法J.空间控制技术与应用,2017,43(4):38-41.SUI Jie,CHENG Huiyan,YU Chengwu,et al.A thermal stability analysis and simulation method for bore
46、sight axis of star sensorJ.Aerospace Control and Application,2017,43(4):38-41.7 董雪岩,许方宇,陈骥,等.1 m 红外太阳望远镜镜面位置误差致光轴偏移分析J.红外技术,2016,38(10):870-876.DONG Xueyan,XU Fangyu,CHEN Ji,et al.Analysis on optic axis deviation induced by lens position disorder of 1 m new vacuum solar telescopeJ.Infrared Technolog
47、y,2016,38(10):870-876.8 Kalikivayi V,Kumar V C P,Kannan K,et al.Tolerance analysis of misalignment in an optical system using ShackHartmann wavefront sensor:experimental studyJ.Optical Engineering,2015,54(7):075104.9 高磊.基于蒙特卡洛法的公差分析及优化设计方法研究D.哈尔滨:哈尔滨理工大学,2015.GAO Lei.Research on the Tolerance Analys
48、is Based on Monte Carlo Method and Optimal Design MethodD.Harbin:Harbin University of Science and Technology,2015.10 刘琳,张兴德,贺谊亮.基于蒙特卡洛模拟法的红外光学系统公差分析J.激光与红外,2020,40(5):497-499.LIU Lin,ZHANG Xingde,HE Yiliang.Monte Carlo simulation and its application in the IR optical systemJ.Laser Infrared,2020,40(5
49、):497-499.11 袁贵星,王平.蒙特卡洛模拟及其在公差设计中的应用J.天津科技大学学报,2008,23(2):61-64.YUAN Guixing,WANG Ping.Monte Carlo simulation and its application in tolerance designJ.Journal of Tianjin University of Science&Technology,2008,23(2):61-64.12 YANG K Y,JIN N,YANG D,et al.Monte Carlo calculation on stray radiation scatt
50、ered by the baffle with a heterotypic surface for uncooled infrared systemC/Proceedings of SPIE,2018,10815:108150R.13 苏朋.365 nm光刻照明系统中变焦系统的设计及公差分析J.红外与激光工程,2022,51(7):20210524.SU Peng.Design and tolerance analysis of the zoom system in 365 nm UV lithography illumination systemJ.Infrared and Laser En