技术领域
[0001] 本
发明涉及反射镜或者透镜镜面面形
精度分析技术领域，特别涉及一种能够直接表征镜面的刚体位移与镜面的柔性
变形的光机集成分析面形误差处理方法。
背景技术
[0002] 为了优化镜体结构，保证镜面的PV值，RMS值满足设计要求，需要采用
有限元分析进行分析。有限元分析仅仅能够得到镜面表面
节点的位移与转
角，不能直接分离出镜体的刚体位移、球差以及面形误差。为了进行光机集成分析，需要采用一定的计算流程来保证计算结果的合理性。
[0003] 为了计算镜面的刚体位移，以往的刚体位移计算方法通常采用均匀压强场计算节点的加权因子或者采用不加加权因子进行计算，或者在球
坐标系中计算刚体位移。之前的刚体位移计算方法仅仅计算一次，最终得到的刚体位移量与面形误差量在表示面形误差时有遗漏。该
专利将镜面的刚体位移、球差以及面形误差分别分离出来，在镜面面形优化设计过程中应用非常方便。
[0004]
现有技术通常不采用加权因子或者加权因子计算方法采用均匀压
力工况下计算的加权因子，该加权因子对在均匀压强或者
真空工作的镜体不适用。通过采用镜体
顶点坐标系下平面投影，采用
重心与各边中点分割的方法将面积分割给三角形或者四边形的顶点节点。该计算方法得到的面积加权因子在刚体位移计算、Zernike系数拟合过程中均可以采用。
[0005] 现有技术在计算刚体位移时均采用单次计算，会忽略x y方向的柔性位移。造成计算的刚体位移不准确。通过该计算方法计算的刚体位移量值更大，而且最终去除刚体位移之后x y方向无柔性位移，表述的面形误差更准确。
[0006] 除此之外还论述了采用减去Zernike系数球差来替代减去镜面
曲率变化的技术。镜面曲率变化计算过程中通常需要利用到镜面公式，该镜面公式通常不是线性的，而拟合Zernike系数的球差项可以通过线性最小二乘法计算得到。通过该方法可以在有限元分析
软件镜面结构优化中使用。
发明内容
[0007] 本发明旨在克服现有技术的
缺陷，提供一种能够直接表征镜面的刚体位移与镜面的柔性变形的光机集成分析面形误差处理方法。
[0008] 为实现上述目的，本发明采用以下技术方案：提供一种光机集成分析面误差处理方法，包括如下步骤：
[0009] (1)
修改光学软件中的镜面的有限元模型，通过添加虚拟面，使分析镜面无倾斜与偏轴，方便导入刚体位移，依据光学软件中镜面的各个顶点坐标系，在有限元模型中同样建立顶点坐标系；
[0010] (2)在有限元模型分析结果中分离出镜面各个节点的刚体位移；
[0011] (3)根据镜面节点在坐标系上的投影，计算出各个节点的加权系数；
[0012] (4)计算镜面的全部刚体位移，镜面的六个刚体位移表示为：Tx,Ty,Tz,Rx,Ry,Rz，这些刚体位移导致单个节点的位移量为：dx,dy,dz：
[0013]
[0014]
[0015]
[0016] (5)去除镜体轴线偏移与倾斜之后的各个节点
位置为dx′i，dy′i，dz′i，由以下公式计算得到各个节点去除刚体轴线偏移与倾斜的节点位置：
[0017] dx′i＝dxi-Tx-ziRy+yiRz
[0018] dy′i＝dyi-Ty+ziRx-xiRz；
[0019] (6)将各个节点剩余柔性变形转换为z向变形，通过泰勒斯公式展开，将上一步计算得到的dx′i，dy′i，转换为各个节点单个z方向的位移，设镜面的公式为：
[0020] 则：
[0021]
[0022] 其中，c表示：圆锥面的曲率，也即是
曲率半径的倒数；
[0023] 其中k表示：圆锥面的圆锥系数。
[0024] (7)通过dx′、dy′与dz′，重新计算镜面的刚体位移，得到镜面的刚体位移T′x,T′y,T′z,R′x,R′y,R′z；
[0025] (8)重新计算镜面的轴线偏移[T′x T′y 0]与倾斜为[ziR′y-yiR′z -ziR′x+xiR′z 0]，减去偏移与倾斜之后节点的x与y值分别为：
[0026]
[0027]
[0028] (9)通过泰勒斯公式展开，将上一步得到的柔性变形 转化到z方向的位移：
[0029]
[0030] (10)利用T′x,T′y,T′z,R′x,R′y,R′z计算z方向减去刚体位移之后的z方向位移：
[0031]
[0032] (11)利用最优化
算法计算镜面曲率变化，设镜面公式为：
[0033]
[0034] 其中，Sj、rj分别表示：曲面上第j个节点z坐标值与该节点在xy平面上[0035] 的投影与原点的距离，
[0036] 拟合球面曲率，球心不在轴线上时，球心位置为 利用如下公式：
[0037]
[0038] 其中E表示理论z坐标值与实际计算z值偏差的加权平方和；
[0039] wj表示面积加权系数；
[0040] c0表示设计的圆锥面曲率；
[0042] 通过最优化寻优算法，求E得最小值，可得到最优的
[0043] 通过拟合球差项来减去曲率变化导致的变形量，Zernike多项式前四项为：
[0044]
[0045] 其中A00为偏移，A11与B11为倾斜项，A20为球差项；ρi与θi为第i个节点极坐标表示的半径与角度，为拟合出Zernike系数的前四项，需使下式最小：
[0046]
[0047] 由：
[0048]
[0049] 得到最优解 与
[0050] (12)计算减去曲率变化或者球差之后的镜面剩余误差值，计算公式如下式所示：
[0051]
[0052] 其中，z方向的刚体位移同时更新为：
[0053] 或：
[0054]
[0055] 采用如下公式计算镜面的RMS值与PV值：
[0056]
[0057] PV＝max(dzerror)-min(dzerror)。
[0058] 在有限元模型中根据镜面单元格属性或者节点位置将镜面的节点分离出来，将各个节点的坐标转换到顶点坐标系中，根据节点编号，在有限元模型分析结果中分离出镜面各个节点的刚体位移。
[0059] 将镜面单元格投影到垂直于顶点坐标系Z轴的平面上，利用三角形或者四边形重心位置与各个边中点分割成小
块面积，将各个小块面积划分给相近的节点，计算出各个节点的加权系数；
[0060] 计算镜面的全部刚体位移，镜面的六个刚体位移表示为：Tx,Ty,Tz,Rx,Ry,Rz，这些刚体位移导致单个节点的位移量为：dx,dy,dz：
[0061]
[0062]
[0063]
[0064] 为拟合得到最佳的刚体位移，需要使下式最小
[0065]
[0066] 由
[0067]
[0068] 可以写为：
[0069]
[0070] 由上式可知刚体位移计算是线性计算，通过逆矩阵就可以计算得到，在计算得到刚体位移之后可以确定该镜面的轴线偏移[Tx Ty 0]与倾斜为[ziRy-yiRz -ziRx+xiRz 0]。
[0071] 通过dx′、dy′与dz′，利用步骤4中的公式重新计算镜面的刚体位移，得到镜面的刚体位移T′x,T′y,T′z,R′x,R′y,R′z。
[0072] 本发明的有益效果在于：
[0073] 光学软件中镜体表面面形精度仅仅依靠光轴方向变化来表示，以往面形计算软件不能明确的区分出轴线偏移与倾斜。
[0074] 本发明面形处理步骤中严格区分出刚体位移
中轴线偏移与倾斜，然后将剩余误差转换到z轴单个方向，然后在该
基础上计算刚体位移，x、y向减去刚体位移误差之后，剩余残差值为零，仅仅剩余z向误差。
[0075] 本发明的面形刚体计算方法中进行了两次刚体位移计算，第二次的刚体位移量与第一次相比，刚体位移量值上更大，能够将更多的柔性变形转换为刚体位移，并不遗漏任何的柔性变形。
[0076] 在得到减去刚体位移之后的面形值之后，拟合镜面曲率变化或者球差，其中采用拟合球差的方法全部为线性运算，计算过程更简单，不需要采用最优化算法寻找最优值。能够在有限元软件中采用MPC实现。
附图说明
[0077] 图1所示为本发明光机集成分析面误差处理方法的流程
框图。
[0078] 图2所示为镜面的有限元模型及顶点坐标系示意图。
[0079] 图3所示为有限元分析得到的各个节点z向位移示意图。
[0080] 图4所示为各个节点的面积加权系数示意图。
[0081] 图5a所示为第一次转换之后的镜面节点z向位移示意图。
[0082] 图5b所示为第二次转换之后的镜面节点z向位移示意图。
[0083] 图6所示为减去刚体位移之后z向各个节点位移示意图。
[0084] 图7所示为减去曲率变化之后z向各个节点位移示意图。
具体实施方式
[0085] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体
实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，而不构成对本发明的限制。
[0086] 首先对本发明中出现的相关技术进行解释：
[0087] 有限元分析：有限元分析是将连续体分成有限部分，用有限的单元格的
应力与应变来反映连续体的应力与应变。
[0088] 面形精度：面型精度是采用精密
抛光表面与理想镜面的偏差量来表征的数据。表征面型精度有两个参数：PV和RMS值。
[0089] 顶点坐标系：以镜面与光轴交点为坐标原点，光轴方向为z向。在光学软件中为各个镜面的局部坐标系。
[0090] 刚体位移：刚性位移就是镜面的形状不发生变化产生的位移
[0091] 球差：不同孔径角的光线交光轴于不同点上，相对于理想象点的位置有不同的偏离。
[0092] 镜面柔性误差：镜面受力后发生的不规则变形，导致镜面与理想镜面发生的偏移，不包含刚体位移。
[0093] 偏轴镜面：使用的镜面区域与光轴z向不对称。
[0094] 如图1所示，本发明实施例提供的光机集成分析面误差处理方法，包括如下步骤：
[0095] 1、如图2所示，修改CODE V、Zemax等光学软件中的光学模型，通过添加各个虚拟面，使分析镜面尽量无倾斜与偏轴，方便导入刚体位移。依据光学软件中的顶点坐标系，在有限元模型中同样建立顶点坐标系。
[0096] 2、在有限元模型中根据单元格属性或者节点位置将镜面表面的节点分离出来。将各个节点的坐标转换到顶点坐标系中。根据节点编号，在有限元分析结果中分离出镜面各个节点的位移。如图3所示，表示各个节点z向位移，该位移是
温度升高导致的。
[0097] 3、如图4所示，将镜体表面单元格投影到垂直于顶点坐标系Z轴的平面上，利用三角形或者四边形重心位置与各个边中点分割面积，将各个小块面积划分给相近的节点。采用这种方法计算出各个节点的加权系数。
[0098] 4、计算镜面的全部刚体位移，镜体的六个刚体位移表示为：Tx,Ty,Tz,Rx,Ry,Rz。这些刚体位移导致单个节点的位移量为：dx,dy,dz
[0099]
[0100]
[0101]
[0102] 为拟合得到最佳的刚体位移，需要使下式最小
[0103]
[0104] 由：
[0105]
[0106] 可以写为：
[0107]
[0108] 由上式可知刚体位移计算是线性计算，通过逆矩阵就可以计算得到。
[0109] 在计算得到刚体位移之后可以确定该镜面的轴线偏移[Tx Ty 0]与倾斜为[ziRy-yiRz -ziRx+xiRz 0]。
[0110] 5、去除镜体轴线偏移与倾斜之后的各个节点位置为dx′、dy′与dz′。由以下公式计算得到各个节点去除刚体轴线偏移与倾斜的节点位置：
[0111] dx′i＝dxi-Tx-ziRy+yiRz
[0112] dy′i＝dyi-Ty+ziRx-xiRz
[0113] 6、如图5a所示为第一次转换之后的镜面节点z向位移示意图。将各个节点剩余柔性变形转换为z向变形。将各个节点剩余柔性变形转换为z向变形，通过泰勒斯公式展开，将上一步计算得到的dx′i，dy′i，转换为各个节点单个z方向的位移，设镜面的公式为：
[0114] 则：
[0115]
[0116] 其中，c表示：圆锥面的曲率，也即是曲率半径的倒数；
[0117] 其中k表示：圆锥面的圆锥系数。
[0118] 7、通过dx′、dy′与dz′，重新计算镜面的刚体位移，得到镜面的刚体位移T′x,T′y,T′z,R′x,R′y,R′z。
[0119] 8、重新计算镜面的轴线偏移[T′x T′y 0]与倾斜为[ziR′y-yiR′z -ziR′x+xiR′z 0]，减去偏移与倾斜之后节点的x与y值分别为：
[0120]
[0121]
[0122] 9、如图5b所示为第二次转换之后的镜面节点z向位移示意图，通过泰勒斯公式展开，将上一步得到的柔性变形 转化到z方向的位移。公式如图同步骤6所示：
[0123]
[0124] 10、如图6所示，利用T′x,T′y,T′z,R′x,R′y,R′z计算z方向减去刚体位移之后的z方向位移：
[0125]
[0126] 11、利用最优化算法计算镜面曲率变化。设镜面公式为 其中，Sj、rj分别表示：曲面上第j个节点z坐标值与该节点在xy平面上的投影与原点的距离，[0127] 拟合球面曲率，球心不在轴线上时，球心位置为 利用如下公式：
[0128]
[0129] 其中E表示理论z坐标值与实际计算z值偏差的加权平方和；
[0130] wj表示面积加权系数；
[0131] c0表示设计的圆锥面曲率；
[0132] 表示最优化算法求得的圆锥面曲率。
[0133] 通过最优化寻优算法，求E得最小值，可得到最优的
[0134] 通过拟合球差项来减去曲率变化导致的变形量，Zernike多项式前四项为：
[0135]
[0136] 其中A00为偏移，A11与B11为倾斜项，A20为球差项；ρi与θi为第i个节点极坐标表示的半径与角度，为拟合出Zernike系数的前四项，需使下式最小：
[0137]
[0138] 由：
[0139]
[0140] 得到最优解 与
[0141] (12)计算减去曲率变化或者球差之后的镜面剩余误差值，如图7所示。计算公式如下式所示：
[0142]
[0143] 其中，z方向的刚体位移同时更新为：
[0144] 或：
[0145]
[0146] 采用如下公式计算镜面的RMS值与PV值：
[0147]
[0148] PV＝max(dzerror)-min(dzerror)。
[0149] 以上所述本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所作出的各种其他相应的改变与变形，均应包含在本发明
权利要求的保护范围内。
