-
根据马吕斯定律,光学系统成完善像的条件是入射波面和出射波面对应点间的光程均为定值[12]。如图1所示,在头戴显示系统中,对于任意一个视场,入射波面为平面波,从平面波发出的平行光经过反射镜反射后到达目标点(理想像点),其中每一条光线都具有相同的光程。可以根据每个视场中由主光线确定的光程来计算该视场自由曲面上的特征数据点,拟合成自由曲面。
图 1 基于马吕斯定律的三维直接设计方法示意图
Figure 1. Diagram of three-dimensional direct design method based on Marius's law
基于马吕斯定律的三维直接设计方法流程如图2所示。首先定义特征光线和特征光线的起点,确定第i+1个视场下未知自由曲面的参考点,即参考光线与前
$ i$ 个视场对应的自由曲面的交点。然后计算第i+1个视场参考光线的光程,并根据马吕斯定律求解该视场未知自由曲面上的特征数据点。将前i+1个视场的所有特征数据点用多项式拟合成自由曲面,作为寻找下一个参考点的参考面。待所有视场的特征数据点都计算完毕,得到含有自由曲面的初始结构。图 2 基于马吕斯定律的三维直接设计方法流程图
Figure 2. Flow chart of three-dimensional direct design method based on Marius's law
(1)确定特征光线和特征光线的起点
如图3(a)所示,选取M个视场,在极坐标下对每个视场的光瞳分别均匀采样K个数据点,对应该视场下K条特征光线;在采样过程中,为了确保采样均匀,在光瞳内采样点个数随着半径的增大而增加。定义特征光线为入射至自由曲面反射镜Ω的平行光
${r}_{ij}\left(i=1,\mathrm{ }\mathrm{ }\cdots ,M;j=0,\mathrm{ }1,\mathrm{ }2,\mathrm{ }\cdots ,K-1\right)$ 以及经自由曲面反射镜的反射光线${r}_{ij}^{\mathrm{\text{'}}}\left(i=1,\mathrm{ }2,\mathrm{ }\cdots ,M;j=0,\mathrm{ }1,\mathrm{ }\mathrm{ }\cdots , K-1\right)$ ,定义特征数据点为特征光线与自由曲面反射镜Ω的交点${P}_{ij}\left(i=1,\mathrm{ }2,\mathrm{ }\cdots ,M;j=0,\mathrm{ }1,\mathrm{ }\mathrm{ }\cdots ,K-1\right)$ 。图 3 (a) 定义光瞳上的采样点;(b) 定义特征光线和特征光线的起点
Figure 3. (a) Define the data points on pupil; (b) Define the feature rays and the starting points of feature rays
如图3(b)所示,对于每个采样视场i,波前平面经过光瞳中心
$ {S}_{i0} $ ,并与主光线$ {r}_{i0} $ 垂直,因此该视场下的所有特征光线与波前平面相交且垂直,交点${S}_{ij}(i=1,\mathrm{ }2,\mathrm{ }\cdots ,M;j=0,\mathrm{ }1,\mathrm{ }\mathrm{ }\cdots ,K-1)$ 定义为特征光线的起点。定义
${T}_{i}(i=1,\mathrm{ }2,\mathrm{ }\cdots ,M)$ 为各视场特征光线在像面上的理想像点,根据物像关系和像面位置计算理想像点的位置。对于来自第i个视场$ \left({\omega }_{x},{\omega }_{y}\right) $ 的特征光线,$ {T}_{i,ideal} $ 的全局坐标为:$$ \begin{split} {T_{i,ideal}} =& \left[ \begin{gathered} {T_{i,ideal,x}} \\ {T_{i,ideal,y}} \\ {T_{i,ideal,z}} \\ \end{gathered} \right] =\\ & \left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&{\cos \alpha } \\ 0&{ - \sin \alpha } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {f\tan {\omega _x}} \\ {f\tan {\omega _y}} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{x_0}} \\ {{y_0}} \\ {{z_0}} \end{array}} \right] \\ \end{split} $$ (1) 式中:
$ f $ 为系统的焦距;$ \left({x}_{0},{y}_{0}{,z}_{0}\right) $ 是像面中心的坐标,即零视场$ \left(0°,\mathrm{ }0°\right) $ 的理想像点坐标,根据无遮拦条件来确定其坐标;$ \mathrm{\alpha } $ 为像面关于x轴的倾斜角度。(2)求解未知自由曲面上的特征数据点
根据特征光线的起点
$ {S}_{ij} $ 、自由曲面反射镜Ω的已知特征数据点$ {P}_{10} $ 和理想目标点$ {T}_{1} $ ,建立相邻两个采样视场之间的迭代关系,计算自由曲面反射镜Ω的所有特征数据点。如图4(a)、图5所示,在第1个采样视场中,定义从出瞳中心
$ {S}_{10} $ 发出的主光线$ {r}_{10} $ 为参考光线,定义第1个采样视场下的特征数据点$ {P}_{1 j} $ 构成的曲面为自由曲面${\varOmega }_{1}$ ;参考光线$ {r}_{10} $ 与自由曲面${\varOmega }_{1}$ 的交点作为参考点$ {P}_{10} $ ,自由曲面${\varOmega }_{1}$ 的位置根据系统尺寸确定;参考光线$ {r}_{10} $ 在参考点$ {P}_{10} $ 处反射后到达理想目标点$ {T}_{1} $ ;计算第1个采样视场下经过自由曲面${\varOmega }_{1}$ 反射的光程:图 4 (a) 自由曲面单反射面光学系统YOZ平面示意图;(b) 全视场下自由曲面构建示意图
Figure 4. (a) Layout of freeform surface single-mirror optical system in YOZ plane; (b) Diagram of freeform surface construction in full FOV
$$ OP = n\left| {{S_{10}}{P_{10}}} \right| + n\left| {{P_{10}}{T_1}} \right| $$ (2) 式中:
$ n $ 为介质的折射率,系统在空气中,则$ n=1 $ ;$ {T}_{1} $ 为理想目标点,通过物像关系求解。根据马吕斯定律,该采样视场中的所有特征光线具有相同的光程,在该视场中求得K个特征数据点${P}_{1 j}(j=\mathrm{0,1},\mathrm{ }\cdots ,K-1)$ ,它满足以下条件:$$ n\left|{S}_{1j}{P}_{1j}\right|+n\left|{P}_{1j}{T}_{1}\right|=OP $$ (3) 将第1个视场下的所有特征数据点
$ {P}_{1 j} $ 拟合成用多项式表征的自由曲面${\varOmega }_{1}$ ;在得到前
$ i\left(1\leqslant i \leqslant M-1\right) $ 个视场下的自由曲面$\sum _{1}^{i}\varOmega$ 后,在第$ i+1 $ 个采样视场中,将从出瞳中心$ {S}_{i+\mathrm{1,0}} $ 发出的主光线$ {r}_{i+\mathrm{1,0}} $ 作为参考光线,定义第$ i+1 $ 个采样视场下的特征数据点$ {P}_{i+1,j} $ 构成的曲面为自由曲面${\varOmega }_{1+1}$ ;参考光线$ {r}_{i+\mathrm{1,0}} $ 与已知自由曲面$\sum _{1}^{i}\varOmega$ 的交点作为参考点$ {P}_{i+\mathrm{1,0}} $ ,参考光线$ {r}_{i+\mathrm{1,0}} $ 在参考点$ {P}_{i+\mathrm{1,0}} $ 处反射后到达理想目标点$ {T}_{i+1} $ ;计算第 i + 1 个视场参考光线的光程,根据马吕斯定律,求解第$ i+1 $ 个采样视场下自由曲面${\varOmega }_{1+1}$ 上特征数据点$ {P}_{i+1,j} $ 。将前
$ i+1\left(1\leqslant i \leqslant M-1\right) $ 个视场下的所有特征数据点${P}_{1 j}\mathrm{、}{P}_{2 j}\cdots {P}_{i+1,j}$ 拟合成用多项式表征的自由曲面$\sum _{1}^{i+1}\mathrm{\varOmega }$ ;重复上述步骤构建前$ i $ 个视场和第$ i+1 $ 个视场下的自由曲面的迭代关系,分别得到M个视场下自由曲面Ω的特征数据点坐标${P}_{1 j}\mathrm{、}{P}_{2 j}\cdots {P}_{i,j}\cdots {P}_{M,j}$ 。(3)将特征数据点拟合成自由曲面
用最小二乘法将自由曲面反射镜Ω的所有特征数据点拟合成用泽尼克多项式(公式(4))表征的自由曲面,得到离轴头戴显示光学系统初始结构;利用评价函数对拟合的多项式系数进行优化,确定所述多项式的最优拟合系数,生成所需的自由曲面反射镜,得到离轴头戴显示光学系统最终结构。
$$ \begin{split} \varphi (x,y) =& \sum\limits_{j = 1}^n {{k_j}{z_j}(x,y)} = {k_1}{z_1}(x,y) + \\ & {k_2}{z_2}(x,y) + \cdots + {k_n}{z_n}(x,y) \\ \end{split} $$ (4) 上述计算方法可用于多个自由曲面反射镜的计算。当系统中存在N个待求自由曲面反射镜时,首先将自由曲面反射镜F1作为第一个待求自由曲面,将其他N−1个反射镜设为已知球面;根据费马原理和马吕斯定律,以及特征光线的起点、相邻两个采样视场下自由曲面反射镜F1的特征数据点和理想目标点,建立相邻两个采样视场之间的迭代关系,在全视场范围内计算第1个自由曲面反射镜F1的特征数据点,将该自由曲面反射镜F1的所有特征数据点进行多项式拟合,得到多项式表征的所述自由曲面反射镜F1。 将自由曲面反射镜F2作为第2个待求自由曲面,保持其他N−2个初始球面反射镜及自由曲面反射镜F1不变,计算得到自由曲面反射镜F2。重复上述步骤,完成N个待求自由曲面反射镜的计算。
以双反射面自由曲面系统为例,如图5所示,P和Q分别在自由曲面反射镜F1和F2上。将自由曲面反射镜F1作为第一个待求自由曲面,F2设为已知球面。根据费马原理,对于任意一个视场,
$ {Q}_{i0} $ 在球面F2上,并使$ {P}_{i0} $ 至理想目标点$ {T}_{i} $ 的距离最近:图 5 自由曲面双反射面光学系统YOZ平面示意图
Figure 5. Layout of freeform surface two-mirror optical system in YOZ plane
$$ n\left| {{P_{i0}}{Q_{i0}}} \right| + n\left| {{Q_{i0}}{T_i}} \right| = \min $$ (5) 计算参考光线的光程
$ OP $ :$$ OP = n\left| {{S_{i0}}{P_{i0}}} \right| + n\left| {{P_{i0}}{Q_{i0}}} \right| + n\left| {{Q_{i0}}{T_i}} \right| $$ (6) 根据马吕斯定律,该采样视场中的所有特征光线具有相同的光程,在该视场中找到其他K−1个特征数据点
$ {P}_{ij}(j=1,\mathrm{ }2,\mathrm{ }\cdots ,K-1) $ ,它满足以下条件:$$ n\left| {{S_{ij}}{P_{ij}}} \right| + n\left| {{P_{ij}}{Q_{ij}}} \right| + n\left| {{Q_{ij}}{T_i}} \right| = OP $$ (7) -
对于大多数头戴显示系统,全视场角通常约为20°~40°[13-15];而太大的视场角可能会引起视觉疲劳,因为眼睛需要扫描较大角度。另外对于单反射面头戴显示系统,当视场角过大时,图像源必须设置为曲面才能满足要求,文中设计采用平面微显示器作为图像源,因此选择中等要求的24°对角视场,用上述方法进行设计,设计指标如表1所示。X方向视场角为19.12°,Y方向视场角为14.4°。
表 1 自由曲面离轴单反头戴显示光学系统设计指标
Table 1. Specifications of freeform surface off-axis single-mirror HMD optical system design
Parameter Specifications Full field of view/(°) 24(19.12×14.4) Exit pupil diameter/mm 3 Effective focal length/mm 14.25 Pixel size/μm 14.1 Distance of exit pupil/mm >15 Spectral range Visible light band Image quality >30%@35 lp/mm Distortion <5% Size of the micro display 0.25 in(6.35 mm) 对165个视场进行采样,分别为:(0°, 0°),(0°, ±1°),···,(0°, ±6°),(0°, ±7.2°); (1°, 0°),(1°, ±1°),···,(1°, ±6°),(1°, ±7.2°); ···(9°, 0°),(9°, ±1°),···,(9°, ±6°),(9°, ±7.2°); (9.56°, 0°),(9.56°, ±1°),···,(9.56°, ±6°),(9.56°, ±7.2°)。除边缘视场外,其他视场按1°间隔进行采样。每个视场选取63条特征光线,共选取10395条不同视场不同孔径位置的特征光线。按前述的设计方法求解全视场的自由曲面特征数据点,并将其拟合成用泽尼克多项式表征的自由曲面,自由曲面矢高如图6(a)所示,拟合系数如表2所示。在光学软件中构建自由曲面头戴显示光学系统初始结构,如图6(b)所示。图6(c)为初始结构的调制传递函数(MTF)图。
图 6 (a) 拟合得到的自由曲面矢高图;(b) 光学系统光路图;(c) 光学系统的调制传递函数曲线
Figure 6. (a) Sag map of freeform surface fitted; (b) Layout of the optical system; (c) MTF curves of the optical system
表 2 自由曲面反射镜的泽尼克多项式系数
Table 2. Zernike polynomial coefficients of freeform surface mirror
Polynomial $ {Z}_{i}(x,y) $ Polynomial coefficients $ {C}_{j} $ $ {Z}_{1} $ 1 $ {C}_{1} $ 18.472 $ {Z}_{2} $ 2$x$ $ {C}_{2} $ 0 ${{Z} }_{3}$ 2$ {y} $ $ {{C}}_{3} $ $ 9.011\times {10}^{-2} $ ${{Z} }_{4}$ $\sqrt{3}\left[2\left({{x} }^{2}+{{y} }^{2}\right)-1\right]$ ${{C} }_{4}$ $ -5.269\times {10}^{-3} $ ${{Z} }_{5}$ $2\sqrt{6}{x}{y}$ ${{C} }_{5}$ 0 ${{Z} }_{6}$ $\sqrt{6}\left({{x} }^{2}-{{y} }^{2}\right)$ ${{C} }_{6}$ $ 1.499\times {10}^{-4} $ ${{Z} }_{7}$ $ \sqrt{8}\left[3\left({{x}}^{2}+{{y}}^{2}\right)-2\right]{y} $ $ {{C}}_{7} $ $ 1.149\times {10}^{-5} $ ${{Z} }_{8}$ $ \sqrt{8}\left[3\left({{x}}^{2}+{{y}}^{2}\right)-2\right]{x} $ $ {{C}}_{8} $ 0 $ {{Z}}_{9} $ $ \sqrt{8}\left(3{{x}}^{2}-{{y}}^{2}\right){y} $ $ {{C}}_{9} $ $ 1.748\times {10}^{-6} $ $ {{Z}}_{10} $ $ \sqrt{8}\left({{x}}^{2}-{3{y}}^{2}\right){x} $ $ {{C}}_{10} $ 0 $ {{Z}}_{11} $ $ \sqrt{5}\left[6{\left({{x}}^{2}+{{y}}^{2}\right)}^{2}-6\left({{x}}^{2}+{{y}}^{2}\right)+1\right] $ $ {{C}}_{11} $ $ -6.282\times {10}^{-7} $ $ {{Z}}_{12} $ $ \sqrt{10}\left({{x}}^{2}-{{y}}^{2}\right)\left[4\left({{x}}^{2}+{{y}}^{2}\right)-3\right] $ $ {{C}}_{12} $ $ 3.439\times {10}^{-8} $ $ {{Z}}_{13} $ $ 2\sqrt{10}\left[4\left({{x}}^{2}+{{y}}^{2}\right)-3\right]{x}{y} $ $ {{C}}_{13} $ 0 $ {{Z}}_{14} $ $ \sqrt{10}\left({{x}}^{4}-{6{x}}^{2}{{y}}^{2}+{{y}}^{4}\right) $ $ {{C}}_{14} $ $ -9.838\times {10}^{-7} $ $ {{Z}}_{15} $ $ 4\sqrt{10}\left({{x}}^{2}-{{y}}^{2}\right){x}{y} $ $ {{C}}_{15} $ 0 根据设计指标设置优化评价函数,用光学设计软件进一步优化,得到如图7(a)所示的最终结构。图7(b)为用三维直接设计方法设计得到的最终系统网格畸变图,最大畸变约为−4.69%。图7(c)为系统的MTF曲线,MTF在35 lp/mm时约为0.3。可以推断出,用基于马吕斯定律的自由曲面三维直接设计方法能够直接设计得到自由曲面单反射面头戴显示光学系统的初始结构,且该结构优化后光学性能满足要求。
-
用上述方法设计了一个双反射面头戴显示系统,设计指标如表3所示。放松状态下的眼睛的瞳孔直径约为2~3 mm,确定出瞳直径为8 mm,全视场角约为28°,相应X方向的视场角为23°,Y方向的视场角为16°。
表 3 自由曲面双反射头戴显示光学系统设计指标
Table 3. Specifications of freeform surface two-mirror HMD optical system design
Parameters Specifications Micro display Micro display OLED Size 0.61 in (12.36 mm×9.36 mm) Resolution/pixel 800×600 Pixel size/μm 15 Head-mounted display system Structure Off-axis two-reflective system Exit pupil diameter/mm 8 F/# 4 Distance of exit pupil/mm >30 Full field of view/(°) 28(23×16) Performance MTF >0.4@30 lp/mm,with a 4 mm exit pupil Maximum distortion <5% 按照第1节的设计流程,设初始结构的两个反射镜都为球面,且曲率半径相等。根据焦距
$ {f}{{{'}}}=32\; \mathrm{m}\mathrm{m} $ ,球面之间的空气间隔$ d=-8\; \mathrm{m}\mathrm{m} $ ,以及公式(8)可计算两个球面的曲率半径为$ {R}_{1}=-120\; \mathrm{m}\mathrm{m} $ ,$ {R}_{2}=120\; \mathrm{m}\mathrm{m} $ 。$$ {f'} = - \frac{{f_1'f_2'}}{D} = \frac{{n{R_1}{R_2}}}{{(n - 1)[n({R_2} - {R_1}) + (n - 1)d]}} $$ (8) 设计结果如图8所示。图8(a)中,两个反射面均为球面的系统作为初始结构;将球面F2固定,计算自由曲面F1,如图8(b)所示;再将已知的自由曲面F1固定,计算自由曲面F2,得到F1和F2均为自由曲面的光学系统,如图8(c)所示。三种系统各视场平均 RMS 光斑半径的对比如图9所示。F1, F2均为球面时,平均 RMS 光斑半径为171.6 μm,此时系统存在较大的场曲和像散;F1为自由曲面时,平均 RMS 光斑半径减小至49.3 μm,这是由于自由曲面F1校正了部分像散;F1, F2均为自由曲面时,成像质量进一步改善,平均 RMS 光斑半径减小至21.5 μm。F1, F2均为自由曲面的系统,可以作为初始结构进一步优化。
图 8 自由曲面双反射面头戴显示系统光路图。 (a) F1, F2均为球面;(b) F1为自由曲面;(c) F1, F2均为自由曲面
Figure 8. Layouts of freeform surface two-mirror HMD system. (a) Both F1 and F2 are spheres; (b) F1 is a freeform surface; (c) Both F1 and F2 are freeform surfaces
图 9 三种光学系统平均 RMS 光斑半径对比
Figure 9. Comparison of the average RMS spot diameter between three optical systems
将图8(c)的自由曲面初始结构进行优化,得到如图10(a)所示的最终光学系统。图10(b)和10(c)显示,双反射面头戴显示系统的最大畸变约为−4.63%,MTF在33 lp/mm时大于0.4。由此可见,用三维直接设计方法设计的自由曲面双反射面头戴显示光学系统作为一个良好的初始结构,优化后其光学性能满足要求。
-
对单反射面头戴显示系统和双反射面头戴显示系统分别进行公差敏感度分析和1000次蒙特卡洛分析,公差项包括镜面加工公差和系统装调公差。单反射面头戴显示系统中,对成像质量影响最大的公差项为自由曲面M1的X方向偏心;双反射面头戴显示系统中,对成像质量影响最大的公差项为自由曲面F1的Y方向偏心和自由曲面F2的Y方向偏心。
通过收紧最敏感公差项,得到如表4和表5所示的公差分配以及表6和表7所示的公差分析结果。其中,当自由曲面M1、F1和F2系数公差为±0.05时,引起面型误差PV变化在1/4 λ范围内。公差分析表明,引入公差后单反射面头戴显示系统全视场MTF均大于0.3,相对于理论值的变化量均小于5.1%;引入公差后双反射面头戴显示系统边缘视场MTF大于0.35,相对于理论值的变化量为12.4%;其他采样视场的MTF均大于0.35,变化量均在10%以下。
表 4 自由曲面单反射面头戴显示光学系统公差分配
Table 4. Tolerance data for freeform surface single-mirror HMD optical system
Manufacture tolerance Alignment tolerance Shape error PV/λ Air thickness
/mmDecenter X
/mmDecenter Y
/mmTilt X/(°) Tilt Y/(°) Tilt Z/(°) M1 1/4 - ±0.2 ±0.2 ±0.2 ±0.2 ±0.2 表 5 自由曲面双反射面头戴显示光学系统公差分配
Table 5. Tolerance data for freeform surface two-mirror HMD optical system
Manufacture tolerance Alignment tolerance Shape error PV/λ Air thickness
/mmDecenter X
/mmDecenter Y
/mmTilt X/(°) Tilt Y/(°) Tilt Z/(°) F1 1/4 ±0.2 ±0.2 ±0.15 ±0.2 ±0.2 ±0.2 F2 1/4 ±0.2 ±0.2 ±0.1 ±0.2 ±0.2 ±0.25 表 6 自由曲面单反射面头戴显示光学系统公差分析结果
Table 6. Tolerance analysis of freeform surface single-mirror HMD optical system
Field X field
/(°)Y field
/(°)Theoretical MTF avg Toleranced MTF avg MTF change 1 0 0 0.378 0.368 2.7% 2 0 7.2 0.416 0.405 2.6% 3 0 −7.2 0.495 0.482 2.6% 4 6.76 0 0.407 0.401 1.5% 5 6.76 7.2 0.356 0.338 5.1% 6 6.76 −7.2 0.535 0.530 0.9% 7 9.56 0 0.455 0.438 3.7% 8 9.56 7.2 0.328 0.315 3.9% 9 9.56 −7.2 0.532 0.528 0.7% 表 7 自由曲面双反射面头戴显示光学系统公差分析结果
Table 7. Tolerance analysis for freeform surface two-mirror HMD optical system
Field X field
/(°)Y field
/(°)Theoretical MTF avg Toleranced MTF avg MTF change 1 0 0 0.516 0.478 7.4% 2 0 8 0.560 0.506 9.6% 3 0 −8 0.574 0.533 7.1% 4 11.5 0 0.453 0.397 12.4% 5 11.5 8 0.579 0.539 6.9% 6 11.5 −8 0.407 0.370 9.1% 7 5.75 0 0.569 0.526 7.6% 8 5.75 8 0.738 0.687 6.9% 9 5.75 −8 0.505 0.462 8.5% -
对自由曲面的面型加工性进行分析,图11为单反射面头戴显示系统自由曲面矢高图和有效通光口径内斜率图,即矢高相对于Y轴的变化率。图12(a)、(b)分别为双反射面头戴显示系统自由曲面F1、F2矢高图。图12(c)、(d)分别为F1、F2在有效通光口径内斜率图。由图可知,矢高在有效通光口径内均不存在突变和拐点,整体变化趋势较为平缓,有利于自由曲面的加工。
Design method of freeform surface for off-axis reflective head-mounted display optical system
-
摘要: 为克服现有离轴头戴显示光学系统设计方法不能直接设计全视场全孔径自由曲面反射镜的问题,提出了基于马吕斯定律的自由曲面三维直接设计方法。首先根据马吕斯定律求解自由曲面在全视场和全孔径范围内的所有特征数据点,然后将特征数据点拟合成用多项式表征的自由曲面,直接获得成像质量良好的自由曲面离轴系统初始结构,最后利用评价函数对拟合的多项式系数进行优化,确定最佳拟合系数,生成所需的自由曲面反射镜,得到离轴头戴显示光学系统最终结构。该方法简化了设计流程,计算效率高。基于提出的方法分别设计了单反射面和双反射面的头戴显示光学系统,单反射面系统的出瞳直径为3 mm,视场角19.12°×14.4°;双反射面系统的出瞳直径为8 mm,视场角23°×16°。设计结果表明,用该方法设计的单/双反射面头戴显示系统,其成像质量良好,系统结构紧凑。公差分析表明,引入公差后单/双反射面头戴显示系统最终可实现全视场调制传递函数(MTF)大于0.3 lp/mm和0.35 lp/mm。Abstract: In order to overcome the shortcomings of the existing direct design method of freeform surface and the large amount of calculation, according to the structural characteristics of the off-axis display system, a three-dimensional direct design method of freeform surface based on Marius's law is proposed. First, according to Marius's law, all the characteristic data points of the freeform surface are solved in the full field of view and the full aperture. Then the characteristic data points are fitted into a freeform surface characterized by polynomials, and the initial structure of freeform off-axis system with good imaging quality is obtained directly. Finally, the polynomial coefficient is optimized with the merit function and the best fitting coefficient is determined. The required freeform surface mirror is generated and the final structure of the off-axis head-mounted display optical system is obtained. This method simplifies the design process and has high calculation efficiency. Based on the proposed method, a single-mirror head-mounted display optical system and a two-mirror head-mounted display optical system are designed with an entrance pupil diameter of 3 mm, a field angle of 19.12°×14.4° and an entrance pupil diameter of 8 mm, a field angle of 23°×16°, respectively. The design results show that single/two-mirror head-mounted display system designed by proposed method has good imaging quality and compact system structure. Tolerance analysis shows that the single/two head-mounted display system can finally achieve MTF greater than 0.3 lp/mm and 0.35 lp/mm in full field of view after introducing tolerances.
-
Key words:
- optical design /
- freeform surface /
- Marius's law /
- off-axis display system
-
表 1 自由曲面离轴单反头戴显示光学系统设计指标
Table 1. Specifications of freeform surface off-axis single-mirror HMD optical system design
Parameter Specifications Full field of view/(°) 24(19.12×14.4) Exit pupil diameter/mm 3 Effective focal length/mm 14.25 Pixel size/μm 14.1 Distance of exit pupil/mm >15 Spectral range Visible light band Image quality >30%@35 lp/mm Distortion <5% Size of the micro display 0.25 in(6.35 mm) 表 2 自由曲面反射镜的泽尼克多项式系数
Table 2. Zernike polynomial coefficients of freeform surface mirror
Polynomial $ {Z}_{i}(x,y) $ Polynomial coefficients $ {C}_{j} $ $ {Z}_{1} $ 1 $ {C}_{1} $ 18.472 $ {Z}_{2} $ 2 $x$ $ {C}_{2} $ 0 ${{Z} }_{3}$ 2 $ {y} $ $ {{C}}_{3} $ $ 9.011\times {10}^{-2} $ ${{Z} }_{4}$ $\sqrt{3}\left[2\left({{x} }^{2}+{{y} }^{2}\right)-1\right]$ ${{C} }_{4}$ $ -5.269\times {10}^{-3} $ ${{Z} }_{5}$ $2\sqrt{6}{x}{y}$ ${{C} }_{5}$ 0 ${{Z} }_{6}$ $\sqrt{6}\left({{x} }^{2}-{{y} }^{2}\right)$ ${{C} }_{6}$ $ 1.499\times {10}^{-4} $ ${{Z} }_{7}$ $ \sqrt{8}\left[3\left({{x}}^{2}+{{y}}^{2}\right)-2\right]{y} $ $ {{C}}_{7} $ $ 1.149\times {10}^{-5} $ ${{Z} }_{8}$ $ \sqrt{8}\left[3\left({{x}}^{2}+{{y}}^{2}\right)-2\right]{x} $ $ {{C}}_{8} $ 0 $ {{Z}}_{9} $ $ \sqrt{8}\left(3{{x}}^{2}-{{y}}^{2}\right){y} $ $ {{C}}_{9} $ $ 1.748\times {10}^{-6} $ $ {{Z}}_{10} $ $ \sqrt{8}\left({{x}}^{2}-{3{y}}^{2}\right){x} $ $ {{C}}_{10} $ 0 $ {{Z}}_{11} $ $ \sqrt{5}\left[6{\left({{x}}^{2}+{{y}}^{2}\right)}^{2}-6\left({{x}}^{2}+{{y}}^{2}\right)+1\right] $ $ {{C}}_{11} $ $ -6.282\times {10}^{-7} $ $ {{Z}}_{12} $ $ \sqrt{10}\left({{x}}^{2}-{{y}}^{2}\right)\left[4\left({{x}}^{2}+{{y}}^{2}\right)-3\right] $ $ {{C}}_{12} $ $ 3.439\times {10}^{-8} $ $ {{Z}}_{13} $ $ 2\sqrt{10}\left[4\left({{x}}^{2}+{{y}}^{2}\right)-3\right]{x}{y} $ $ {{C}}_{13} $ 0 $ {{Z}}_{14} $ $ \sqrt{10}\left({{x}}^{4}-{6{x}}^{2}{{y}}^{2}+{{y}}^{4}\right) $ $ {{C}}_{14} $ $ -9.838\times {10}^{-7} $ $ {{Z}}_{15} $ $ 4\sqrt{10}\left({{x}}^{2}-{{y}}^{2}\right){x}{y} $ $ {{C}}_{15} $ 0 表 3 自由曲面双反射头戴显示光学系统设计指标
Table 3. Specifications of freeform surface two-mirror HMD optical system design
Parameters Specifications Micro display Micro display OLED Size 0.61 in (12.36 mm×9.36 mm) Resolution/pixel 800×600 Pixel size/μm 15 Head-mounted display system Structure Off-axis two-reflective system Exit pupil diameter/mm 8 F/# 4 Distance of exit pupil/mm >30 Full field of view/(°) 28(23×16) Performance MTF >0.4@30 lp/mm,with a 4 mm exit pupil Maximum distortion <5% 表 4 自由曲面单反射面头戴显示光学系统公差分配
Table 4. Tolerance data for freeform surface single-mirror HMD optical system
Manufacture tolerance Alignment tolerance Shape error PV/λ Air thickness
/mmDecenter X
/mmDecenter Y
/mmTilt X/(°) Tilt Y/(°) Tilt Z/(°) M1 1/4 - ±0.2 ±0.2 ±0.2 ±0.2 ±0.2 表 5 自由曲面双反射面头戴显示光学系统公差分配
Table 5. Tolerance data for freeform surface two-mirror HMD optical system
Manufacture tolerance Alignment tolerance Shape error PV/λ Air thickness
/mmDecenter X
/mmDecenter Y
/mmTilt X/(°) Tilt Y/(°) Tilt Z/(°) F1 1/4 ±0.2 ±0.2 ±0.15 ±0.2 ±0.2 ±0.2 F2 1/4 ±0.2 ±0.2 ±0.1 ±0.2 ±0.2 ±0.25 表 6 自由曲面单反射面头戴显示光学系统公差分析结果
Table 6. Tolerance analysis of freeform surface single-mirror HMD optical system
Field X field
/(°)Y field
/(°)Theoretical MTF avg Toleranced MTF avg MTF change 1 0 0 0.378 0.368 2.7% 2 0 7.2 0.416 0.405 2.6% 3 0 −7.2 0.495 0.482 2.6% 4 6.76 0 0.407 0.401 1.5% 5 6.76 7.2 0.356 0.338 5.1% 6 6.76 −7.2 0.535 0.530 0.9% 7 9.56 0 0.455 0.438 3.7% 8 9.56 7.2 0.328 0.315 3.9% 9 9.56 −7.2 0.532 0.528 0.7% 表 7 自由曲面双反射面头戴显示光学系统公差分析结果
Table 7. Tolerance analysis for freeform surface two-mirror HMD optical system
Field X field
/(°)Y field
/(°)Theoretical MTF avg Toleranced MTF avg MTF change 1 0 0 0.516 0.478 7.4% 2 0 8 0.560 0.506 9.6% 3 0 −8 0.574 0.533 7.1% 4 11.5 0 0.453 0.397 12.4% 5 11.5 8 0.579 0.539 6.9% 6 11.5 −8 0.407 0.370 9.1% 7 5.75 0 0.569 0.526 7.6% 8 5.75 8 0.738 0.687 6.9% 9 5.75 −8 0.505 0.462 8.5% -
[1] Zhang Nan, Chang Junlei, Li Qinglin, et al. High-stability opto-mechanical structure design of dual channel off-axis remote sensing camera [J]. Infrared and Laser Engineering, 2019, 48(4): 0418005. (in Chinese) doi: 10.3788/IRLA201948.0418005 [2] Wang J, Zhou Q, Chen J, et al. Design of a see-through off-axis head-mounted-display optical system with an ellipsoidal surface [J]. Current Optics and Photonics, 2018, 2(3): 280-285. [3] Nie Y, Mohedano R, Benítez P, et al. Multifield direct design method for ultrashort throw ratio projection optics with two tailored mirrors [J]. Applied Optics, 2016, 55(14): 3794-3800. doi: 10.1364/AO.55.003794 [4] Cheng D, Wang Y, Hua H, et al. Design of an optical see-through head-mounted display with a low f-number and large field of view using a freeform prism [J]. Applied Optics, 2009, 48(14): 2655-2668. doi: 10.1364/AO.48.002655 [5] Yang T, Cheng D, Wang Y, et al. Design method of non-symmetric imaging systems consisting of multiple flat phase elements [J]. Optics Express, 2018, 26(19): 25347-25363. doi: 10.1364/OE.26.025347 [6] Chen Bingxu, Liao Zhiyuan, Cao Chao, et al. Design of the freeform imaging system with large field of view and large relative aperture [J]. Infrared and Laser Engineering, 2020, 49(8): 20200005. (in Chinese) doi: 10.3788/IRLA20200005 [7] Pan J W, Che-Wen C, Huang K D, et al. Demonstration of a broad band spectral head-mounted display with freeform mirrors [J]. Optics Express, 2014, 22(11): 12785-12798. doi: 10.1364/OE.22.012785 [8] Liu J, Miñano J C, Benítez P, et al. Single optical surface imaging designs with unconstrained object to image mapping[C]//SPIE, 2012, 8550: 855011. [9] González-Acua R G, Chaparro-Romo H A, Gutiérrez-Vega J C. General formula to design a freeform singlet free of spherical aberration and astigmatism [J]. Applied Optics, 2019, 58(4): 1010-1015. doi: 10.1364/AO.58.001010 [10] 苏宙平. 一种离轴多反自由曲面光学系统初始结构的设计方法: 中国, CN110927964 A[P]. 2020-03-27. [11] Wu W, Wang H, Jin G, et al. Fast automatic design method for freeform imaging systems through system construction and correction [J]. Optics Letters, 2020, 45(18): 5140-5143. doi: 10.1364/OL.398924 [12] 张以谟. 应用光学[M]. 北京: 电子工业出版社, 2008. Zhang Y M. Applied Optics[M]. Beijing: Publishing House of Electronics Industry, 2008. (in Chinese) [13] Cakmakci O, Rolland J P. Head-worn displays: A review [J]. Journal of Display Technology, 2006, 2(3): 199-216. doi: 10.1109/JDT.2006.879846 [14] Cakmakci O, Rolland J P. Examples of HWD Architectures: Low-, Mid-, and Wide-field of View Designs[M]. Berlin: Springer-Verlag, 2016. [15] Kress B C. Optical Architectures For Augmented-, Virtual-, and Mixed-reality Headsets [M]. Bellingham: SPIE Press, 2020.