CN1132042C - 多焦光学元件的设计方法 - Google Patents
多焦光学元件的设计方法 Download PDFInfo
- Publication number
- CN1132042C CN1132042C CN008002320A CN00800232A CN1132042C CN 1132042 C CN1132042 C CN 1132042C CN 008002320 A CN008002320 A CN 008002320A CN 00800232 A CN00800232 A CN 00800232A CN 1132042 C CN1132042 C CN 1132042C
- Authority
- CN
- China
- Prior art keywords
- function
- group
- steps
- partiald
- derivative
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02C—SPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
- G02C7/00—Optical parts
- G02C7/02—Lenses; Lens systems ; Methods of designing lenses
- G02C7/06—Lenses; Lens systems ; Methods of designing lenses bifocal; multifocal ; progressive
-
- G—PHYSICS
- G02—OPTICS
- G02C—SPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
- G02C7/00—Optical parts
- G02C7/02—Lenses; Lens systems ; Methods of designing lenses
- G02C7/06—Lenses; Lens systems ; Methods of designing lenses bifocal; multifocal ; progressive
- G02C7/061—Spectacle lenses with progressively varying focal power
-
- G—PHYSICS
- G02—OPTICS
- G02C—SPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
- G02C7/00—Optical parts
- G02C7/02—Lenses; Lens systems ; Methods of designing lenses
- G02C7/024—Methods of designing ophthalmic lenses
- G02C7/028—Special mathematical design techniques
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Mathematical Physics (AREA)
- Eyeglasses (AREA)
Abstract
一种设计多焦光学元件表面的方法,包括把一个区域划分为多个三角形(τi);定义一组函数,每个函数表示一个所述三角形(τi)的表面(li)一部分;和优化所述函数。另一实施例包括把一个区域划分为多个多边形(200);定义一组函数,每个函数表示一个多边形(202)的表面一部分。在多边形边界处为所述函数值和函数的一阶导数值设定连续性限制。对所述函数在多边形边界处的二阶导数值不需要连续性限制。在满足所设连续性限制的情况下对函数进行优化。
Description
技术领域
本发明涉及一种多焦光学元件的设计方法。
背景技术
光学元件是那些具有光的反射或折射表面、基本上透明的材料部件,例如:反射镜、透射镜、分束器和准直器。光学元件用于各个方面,包括:望远镜、显微镜、照相机和眼镜。
光学元件以其光学性能和表面光学特性为特征。当入射光通过该光学元件时,它的光学性能例如像散、光焦度和光谱表明入射光波前在通过光学元件时是如何变形的。表面光学特性例如表面像散、表面光焦度以及表面像散和表面光焦度的梯度表示出该光学元件的表面几何特性。光学性能和表面光学特性是密切相关的,但不是相同的。
多焦光学元件具有多于一个光焦度。例如:双焦透镜具有两个子区,每个子区具有不同的光焦度。双焦眼镜可以用于例如:在一个子区矫正近视(近视),而在另一个子区矫正远视(眼睛改变其晶状体形状的能力丧失)。但遗憾的是,由于从一个子区向另一个子区转换时,这种光焦度的突变使得许多人发现配戴双焦眼镜不舒适。这促使了渐变眼镜的发展,这种渐变眼镜是一种多焦眼镜,其中在该透镜上,从一个点到另一个点的光焦度的变化是平滑的。
光学元件表面上的任意一点的表面光焦度由该表面的平均曲率来定义。一个渐变透镜具有变化的光焦度,所以,它具有可变的曲率,因而被定义为非球面。然而,由于这种渐变透镜的表面,或至少它的一关键部分被定义为非球面,所以,它在许多点都具有两个不同的主曲率k1和k2。该光学元件表面的任意一点的表面像散由这两个主曲率k1和k2之差的绝对值定义。
表面f上的点(x,y)处的主曲率H、高斯曲率G和主曲率k1和k2的定义由方程式1A-1D给出:
由于早期在设计渐变透镜时,主要设计目的在于实现:
a)光焦度缓慢变化;
b)使像散最小;
c)降低各种光学象差例如螺纹畸变、双目不均衡、等等。
已经提出了许多不同的方法以实现这些目的。Maitenaz的美国专利US3687528描述了这样一种技术,其中本底曲线(子午线)从该透镜的上部延伸到它的下部。该透镜的表面由该子午线确定,使得曲率逐渐变化(光焦度也因此而变化)。沿着该子午线本身,主曲率k1和k2满足k1=k2。透镜表面以不同的方法在水平方向上从子午线延伸。对于从该子午线的这种延伸给出了明确的公式。Maitenaz得到了具有相当稳定的光焦度的在上部的一个区域和在下部的另一个区域。而且,在该子午线附近的像散相对很小。
许多渐变透镜的设计把该渐变透镜明确地分为三个区:用于远视的上部区、用于近视的下部区和连接以上两个区的过渡区。上部区和下部区提供所必要的清晰视觉。对于上部区和下部区许多设计都采用球形表面。设计过程中的主要工作是确定一个好的过渡区。
Guilino的美国专利US4315673和Winthrop的美国专利US4861153描述了一种通过应用明确的过渡区公式而得到一个平滑转换区域的方法。1986年出版的“蔡司信息97”(Zeiss information 97),第55-59页,Furter和G.Furter的美国专利US4606622“Zeiss Gradal HS-具有最大配戴舒适度的渐变附加透镜”描述了一种方法,其中该透镜设计者在多个特殊点确定了过渡区中透镜表面的值。然后,整个表面就由这种仿样函数的方法生成。为了提高所生成的表面的性能,设计者对某些特定点的透镜表面值作了修正。
然而,Barkan等人的美国专利US4838675描述了另一种方法。用一个基准表面函数来描述一个具有用于远视的上部区、用于近视的下部区以及过渡区的渐变透镜。利用定义该透镜中的一个子区的优化函数来计算一种改进型的渐变透镜,在该子区中这个优化函数被加入到这个基准表面函数中。
在1998年的“计算机辅助设计30”(Computer Aided Design 30),第595-602页由J.Loos、G.Greiner和H.P.Seidel发表的题为“渐变透镜设计的一种变化方法(A variational approach to progressive lens design)”和在1994年的出版的由P.J.Laurent等人编著的“几何设计中的曲线和表面”(Curves and Surfacesin Geometric Design)一书的第467-474页中由M.Tazeroualti所写的题为“渐变透镜的设计”一文中描述了一种不同的技术。该透镜表面由多个仿样函数方程的组合所确定,因而,该表面必须被认为是一个矩形,而它被分为更小的矩形。这种方法对于要确定那些形状不是矩形的透镜时不能直接适用。确定一个成本函数并考虑到仿样函数因子以使该表面的成本函数降到最小。这种方法不会利用到该表面的边界条件,因而,不能用该方法设计那些在边界需要特定形状的透镜。这种方法采用三次仿样函数提供一个h4的一个准确值,其中h是小矩形对角线与大矩形对角线之比。
Kaga等人的欧洲专利申请EP744646描述了这样一种方法,其中该渐变透镜的表面被划分为多个矩形。在这些矩形的边界,该表面必须是连续的,一阶导数连续而且二阶导数也连续。由于光焦度与曲率相关,而曲率是由该表面的二阶导数所确定的,所以,对于二阶导数连续性的限制看上去是很自然的。实际上,许多设计渐变透镜的方法都包括了该透镜表面的二阶导数的连续性限制。
发明内容
本发明的目的在于提供一种新的设计多焦光学元件的方法。
因而根据本发明的优选实施例提供了一种设计多焦光学元件表面的方法。该方法包括这样的步骤,即把一个区域划分为多个三角形,定义一组函数,其中每个函数表示一个三角形表面的一部分,优化这些函数,从而得到该表面。
而且,根据本发明的一个优选实施例,该光学元件为透镜。
另外,根据本发明的一个优选实施例,该透镜为渐变眼镜。
另外,根据本发明的一个优选实施例,所说的区域是二维区域,而且这些函数是二维的五次多项式。
另外,根据本发明的一个优选实施例,该方法还包括基于所生成的表面而进行的模拟该光学元件光学特性的步骤,如果所生成的表面不满意,将调整这些函数的参数值而重新进行优化和模拟的步骤。
此外,根据本发明的一个优选实施例,该优化步骤还包括确定这组函数的一个函数并优化该函数的步骤。
此外,根据本发明的一个优选实施例,该函数是一个成本函数,优化该函数的步骤是使该函数值最小的步骤。
而且,根据本发明的一个优选实施例,该方法还包括在该函数中所含的加权因子进行选择步骤。
另外,根据本发明的一个优选实施例,这些加权因子是根据该表面上至少一个位置、该组函数和这组函数的一阶导数而变化的。
根据本发明另一个优选实施例还提供了一种设计多焦光学元件表面的方法。该方法包括的步骤有:把一个区域划分为多个多边形,定义一组函数,其中每个函数表示一个多边形表面的一部分。在这些多边形的边界上决定这些函数值和这些函数的一阶导数值的连续性限制。在多边形边界上,这些函数的二阶导数值的连续性限制不需要。优化这些函数,以满足所规定的连续性限制,由此得到该表面。
而且,根据本发明的一个优选实施例,该方法还包括决定多边形顶点处该函数的二阶导数值的连续性限制的步骤。
另外,根据本发明的一个优选实施例,所说的多边形为三角形。
根据本发明又一个优选实施例还提供了一种设计多焦光学元件表面的方法。该方法包括的步骤有:把一个区域划分为多个多边形,定义一组函数,其中每个函数表示一个多边形表面的一部分。由多边形的边界所组成的子集确定一组曲线。至少一组含有该函数、该函数的一阶导数和该函数的二阶导数、或其任何组合的值是由选定的、与该组曲线相应的多边形的顶点所规定的。优化这些函数,以满足所规定的值,由此得到该表面。
而且,根据本发明的一个优选实施例,所说的多边形为三角形。
另外,根据本发明的一个优选实施例,所说的多边形为矩形。
根据本发明的另一个优选实施例还提供了一种具有复表面的多焦光学元件,其中至少一个表面是由以上所述方法设计的。
附图说明
通过下面结合附图进行详细的描述,可以对本发明有全面地理解。其中:
图1是根据本发明优选实施例的一个光学元件表面的设计过程的示意流程框图;
图2是根据本发明优选实施例而得的如图1所示的设计过程中的一部分的示意流程图;
图3是根据本发明优选实施例而得的一个透镜表面的示意图;
图4是为帮助理解本发明所作的划分为多个三角形的一个区域的示意图。
具体实施方式
本发明是一种设计多焦光学元件的新方法。该方法可使得设计者把重点放在该光学元件表面的特殊区域,其中得到满意的光学特性是更重要的。该方法也使得设计者可以控制沿着表面曲线和该曲线的邻域内的光学性能。
在设计一种渐变透镜的特定情况下,本发明不需要把该透镜分为三个区:用于远视的上部区、用于近视的下部区和连接以上所说的两个区的过渡区。相反,从整体上看整个透镜表面最佳。
本发明的方法涉及一种带有反馈的设计过程。这种方法如图1中所示,现在参考该图进行描述。设计者拟定在该设计过程中所需的变量值(步骤100),然后,这些参数用于计算该表面(步骤102)。在步骤104中,设计者根据步骤102的计算结果来评价该表面。设计者检查该表面的光学特性是否满意。表面光学特性包括表面光焦度、表面像散性以及该表面光焦度和表面像散性的梯度,其中表面光焦度由该表面的主曲线确定,表面像散性由两个主曲率之差的绝对值确定。而且,设计者可以用一个模拟器来模拟至少一个通过该光学元件的波前,得出该光学元件的光学特性结果。如果,所计算的表面不满意(步骤106),设计者可以修正这些变量值(步骤108),然后,重新计算该表面(步骤102)。重复该程序直到设计者对所计算的表面满意为止。
图2是该设计过程的一部分的流程示意图,在图2中对参数值的指定(步骤100)和表面的计算(步骤102)的表示更为详细,现在,将参考它进行解释。
图3是透镜的示意图,在以下的描述中还将参考图3。该透镜具有一个已知的表面g和一个未知的表面f。该透镜表面f被认为是在(x,y)平面内的一个区域D。该区域D被划分为多边形的集合(步骤200)。在本发明的一个优选实施例中,该多边形是三角形τi,i=1,...N,其中N是三角形的总数。三角形的数目不受限制,而且这些三角形也不需要相同,这样为任意形状的区域D的划分提供了灵活性。
表面f是由一组“微透镜”li的拼合而成,其中的每个微透镜li均由一个三角形τi确定。在该三角形τi中点(x,y)上方的表面f的高度由方程2中所表示的li(x,y)给出:
f(x,y)=li(x,y), (x,y)∈τi (2)
由于三角形不必互相相同,所以小的微透镜可以用于那些需要更为细微的细部设计的透镜表面区域,而且,大的微透镜可以用于那些粗略的细部设计就足够的透镜表面区域。
每个微透镜li表面由一个五次多项式表示(步骤202),这样遍及它所处的三角形,这就是说,它是一个含有x、y的多项式,其中所包括的xjym形式的组合以j+m来算小于或等于5。在该三角形τi中点(x,y)上方的表面高度由方程3中的li(x,y)给出:
li(x,y)=ai 00+ai 10x+ai 01y+ai 20x2+ai 11xy+ai 02y2
+ai 30x3+ai 21x2y+ai 12xy2+ai 03y3+ai 40x4 (3)
+ai 31x3y+aj 22x2y2+ai 13xy3+ai 04y4+ai 50x5
+ai 41x4y+ai 32x3y2+ai 23x2y3+ai 14xy4+ai 05y5方程3可以重新写为方程4那样的更简洁的形式: 其中,ai jm是第i个微透镜li的因子。表面f因此而由这套局域多项式中的因子ai jm全部确定。
我们将会理解,如果把该区域D划分为例如1000个三角形,那么,确定该表面f时将涉及到确定21000个因子的值。不加上别的限制,将会出现表面f的许多种可能的解。
正如以上所述,最自然的限定就是这些微透镜表面li被拼接起来后使得所得到的表面f沿着该三角形τi的边界是连续的、微分可连续、而且二次微分可连续。由于f的二阶导数是曲率,而该曲率与光焦度相关,而且光焦度又应该是平滑变化的,所以,自然需要该表面f沿各边界二阶导数连续。
本发明的一个重要的特征在于只要微透镜表面li被拼接起来后使得所得到的表面f沿着该三角形τi的边界是连续的且是微分可连续的就足以。这在方程5A--5C中表示:
li(x,y):lj(x,y),(x,y)∈τi∩τj; (5A)
其中,三角形τi和τj的相交处即为这两个三角形的公共边界。
根据本发明的另一个优选实施例,微透镜表面li被拼接起来后使得所得到的表面f在该三角形τi的顶点是二次微分可连续的条件加入到方程5A--5C的连续性限制中,这由方程5D—5F表示:
其中,三角形τi和τj和τk的相交处即为这三个三角形的公共顶点。将会发现,在一定的极端情况下,其中的一个顶点仅为两个三角形的顶点时,方程5D—5F的连续性限定同样适用于这些顶点。将会发现,方程5D—5F的连续性限定与已有技术中的由沿整个边界所给出的二阶导数的连续性限定相比,其约束性要小很多。
由于目的是得到一个光焦度变化的渐变透镜,所以,限定该表面f使它的主曲率尽可能与预先所确定的函数P(x,y)相匹配是合理的,该函数P(x,y)用于确定在每个点(x,y)的适宜的光焦度(步骤204)。另外,由于另一个目的是把该渐变透镜的像散降到最小,所以,限定该表面f以使它的主曲率之差为最小也是合理的。而且,如前所述,设计者可能希望关注该透镜表面的某些特定区域。这由加权因子wi,1和wi,2的选择(步骤204)来实现,加权因子表征每个微透镜li的相对“强度”以及每个微透镜li的光焦度和像散性的相对重要性。所关注的微透镜的加权因子将大于次要的微透镜的加权因子。对于一个特定的微透镜li,其中光焦度比像散性更重要,加权因子wi,2将大于加权因子wi,1。与之类似地,对于一个特定的微透镜li,其中像散性比光焦度更重要时,加权因子wi,1将大于加权因子wi,2。
因此,以方程5A—5F为条件,用于表示该表面f的因子ai jm必须使由方程式6所给出的成本函数E最小:
其中,k1(x,y)和k2(x,y)是在点(x,y)处的表面f的主曲率,dxdy是三角形τi表面区域单元。预先确定的加权因子wi,1和wi,2依赖于x,y,可能还依赖于f和f的一阶导数。这种对f和f的一阶导数的依赖可能来自于多种情形。例如:如果积分是在微透镜li表面进行的,而不是在整个三角形τi中进行的,那么,表面区域单元dxdy将由
代替。因子 可以包括在加权因子中,使该表面区域单元为dxdy。对于表面面积大于别的表面面积的微透镜,该因子
将有效地增加微透镜的权重。依赖于f和f的一阶导数的的加权因子的另一个例子是对于当设计者想在该表面f上作突变时的情况。
使成本函数E最小的表面f不必在点(x,y)处具有光焦度P(x,y)。如果设计者想使该透镜具有特殊的形状和/或在某些点具有特定的光学性能,那么,设计者将确定(步骤208)一沿三角形τi边界子集的有限曲线组Cm,m=1,...,Nc,其中,Nc是这些曲线的总数。设计者指定(步骤208)曲线条件,它们将决定该表面f的任何或一些值以及在曲线Cm顶点Vc处的一阶和二阶导数。这样有效地限定了透镜的形状以及顶点Vc处和曲线Cm邻域的光学性能。将会看到,曲线Cm不必在该表面上划定子区。
根据平板加工技术领域中的公知技术,表面f可以被表示为由整个三角形τi定义的多个五次多项式的拼接,这样表面f将满足方程式5A-5F的连续性限定。在J.H.Argyris,I.Fried和D.W.Scharpf发表在1968年的“航空杂志”(TheAeronautical Joumal)第72卷第701-709页的“The TUBA Family of Plate Elementsfor the Matrix Displacement Method”一文中描述该技术,在此作为参考。建立满足方程式5A—5F连续性的Ns个五次形状多项式Sk(x,y)。根据本发明的一个优选实施例,表面f由方程式7所给出的这些形状多项式的线性组合来表示:
其中dk是该线性组合中的未知因子。然后,通过设计,方程式7的表面f满足方程式5A—5F的连续性限制。如以下所述,由设计者所规定的曲线参数变动范围有效地确定了一些未知因子dk的值。当方程式7的表面f被替代为方程式6的成本函数E时,结果是关于因子dk的成本函数E的表达式(步骤206),其中的一些dk是未知的。现在,为了确定余下的那些未知因子dk,成本函数E的最小化可以不必关心方程式5A—5F的连续性限定和曲线的限制就可进行。
正如Argyris等人的文章中所提到的那样,方程式5A—5F的连续性限制相当于满足f和在三角形τi顶点处它本身它的一阶导数和二阶导数的连续性以及在三角形τi各边中点处相对于其法线f的导数的连续性。所以,对于全套的三角形τi,方程式5A—5F的连续性限定相当于对于每个顶点有6个连续性限定和对每个中点有1个连续性限定。F的6个精确值、它的每个顶点处的一阶和二阶导数以及与每个中点法线相关的f的导数的精确值是自变量,被认为是自由度。
将会发现,由于许多顶点是至少一个三角形的共用顶点,所以,三角形τi组中的顶点数Nv小于三角形数N的3倍。我们还会看到,三角形τi组中的中点数Nm也小于三角形数N的3倍。全套三角形τi的总的自由度数由方程8给出:
NS=6Nv+NM (8)
图4是一个被划分为多个三角形的区域的示意图,在此以它作进一步的参考,正如图4所示,为每一个自由度建立一个形状多项式Sk(x,y)。三角形τ1具有三个顶点V1、V2和V3以及三个中点M1、M2和M3。三角形τ2具有三个顶点V1、V2和V4以及三个中点M1、M4和M5。顶点V1为6个三角形τ1-τ6的共用顶点,中点M1是三角形τ1和τ2的公共中点,中点M2是三角形τ1和τ7的公共中点,中点M3是三角形τ1和τ6的公共中点。
其中的一个自由度Ns与顶点V1处的表面f值有关。建立一个形状多项式S1(x,y)以使它的值在顶点V1处为1而在该套三角形中的其它任何一个顶点处为0,在每一个顶点它的一阶和二阶导数为0,而且,与每个中点处的法线相关的导数为0。这将通过确定每个三角形中的一个五次多项式和为这些五次多项式的以片计算的组合设定值S1(x,y)来完成。三角形τ1相对于它具有21个自由度,它包括f值,它在每个顶点V1、V2和V3的一阶和二阶导数以及与每个中点M1、M2和M3处的法线相关的f的导数值。所以,对于三角形τ1,通过设定在顶点V1处的值为1,在顶点V2和V3处的值为0,在顶点V1、V2和V3处的一阶和二阶导数为0以及在中点M1、M2和M3处与法线相关的导数为0而可以完全确定一个具有21个因子的关于x和y的五次多项式。与此类似,对于三角形τ2,通过设定在顶点V1处的值为1,在顶点V2和V4处的值为0,在顶点V1、V2和V4处的一阶和二阶导数为0以及在中点M1、M4和M5处与法线相关的导数为0而可以完全确定一个关于x和y的五次多项式。将会看到,关于三角形τ1的五次多项式和关于三角形τ2的五次多项式满足方程5A—5F的连续性限定。对于该组中的其余的三角形可以完全确定出类似的关于x和y的五次多项式,形状多项式S1(x,y)被定义为这些五次多项式的计件方式的组合。
相似地,建立第二形状多项式S2(x,y)以使在每个顶点处的值为0,在顶点V1处的一阶导数为1而在别的顶点处为0,它的二阶导数在每个顶点处为0,在每个中点处与法线相关的导数为0。
将会看到,通过设定,每个形状多项式Sk(x,y)都满足方程5A-5F的连续性限制,因此,由方程7所给出的线性组合也满足这些条件。
还将看到,只有在那些在该三角形的顶点或中点处自由度为1的三角形中,每个五次形状多项式Sk(x,y)才是非零的。
设计者所确定的曲线的限制对于f值和在该曲线Cm的顶点Vc处它的一阶和二阶导数的一些可能值是指定的值。每个曲线限制确定f或在该套三角形中某个特定顶点的一阶或二阶导数,因而,将完全确定(步骤208)在方程7的线性组合中所出现的某个特定因子dk的值。其余未知的因子dk的数量由曲线限制的数目Ns来表示。
在平板加工技术领域中公开了建立满足方程5A—5F的连续性限制的五次形状多项式的其它技术。例如:在G.R.Cowper、E.Kosko、G.M.Lindberg和M.D.Olson发表在1969年的“AIAA Joumal”的第7卷第10篇第1957--1965页中题为“Static and Dynamic Applications of a High-Precision Trianglar PlateBending Element”一文中描述了这样的技术,在此作为参考。
由Argyris的形状多项式提供一个h6的精确值,而由Cowper的形状多项式提供一个h5的精确值,其中,h是一个三角形的边界与区域D的一个尺寸的比值。
将会看到,任一个满足方程5A--5F连续性限制的五次多项式都可以用作对于f的一个线性组合的基本成分,而代替以上所述的形状函数。
将会看到,除三角形以外的多边形用于划分区域D,与带多边形相关的自由度的数量大于21,所以,考虑到所有自由度,多项式将不得不高于五次。
如前所述,方程7的线性组合被替代为方程6的成本函数E,其结果是关于因子dk的成本函数E的表达式,其中的一些dk是未知的。现在,为了确定余下的那些未知因子dk,成本函数E的最小化可以不必关心方程式5A—5F的连续性限定和曲线的限制就可进行。从变量的微分计算中我们可以知道:相当于解一个特定的方程,即对于未知的表面f的一个已知的Euler-Lagrange方程。在I.M.Gelfand和S.V.Fomin在Prentice Hall于1963年发表的“变量微分学”(Calculus of Variations)一书中对此进行了解释。
对于方程6所给出的成本函数,Euler-Lagrange方程是非线性的。而且,整个三角形τi的积分包括一个关于未知因子dk的非线性表达式,由此要清楚地计算出来是很困难的。所以,应用一种迭代法来求解(步骤210)这个关于未知因子dk的Euler-Lagrange方程。在本发明的一个优选实施例中,采用的是Newton法,但是,也可以采用其它适宜的迭代方法。
迭代法包括生成表面序列f1,f2,…,使得该序列汇聚成所要的方案f。对于该迭代法,选择初始表面f1的方法有很多,例如:可以选一个在x和y平面上的面作为该初始表面f1。该序列中的每个表面都是线性代数方程组,该方程组以一个矩阵为特征。对于这种出现在迭代法中的矩阵的建立在有限元文献中进行了描述,例如:TJ.R.Hughes在Prentice Hall发表的于1987年所写的“有限元法”(The finite elementsMethod)一书和由W.G.Strang和G.Fix合著的在Wellesley Cambridge发表的于1973年所写的“有限元法的分析”(An analysis tothe finite elementsMethod)一书,在此引用作为参考。
利用有限元法的一个好处是用于迭代过程中的矩阵具有一个特定的结构,其中,除了靠近对角线的窄带外其余所有的矩阵元素都为零。这种矩阵被认为是“带矩阵”(band matrices)。这种特征大大地降低了求解方程组的复杂性,而且因此,当要表示透镜表面时,可以用很多的自由度。
原则上说,迭代过程包括无限多步,因而,需要确定一个标准以使在某个特定点的迭代过程中止,而且以所得到的表面作为结果。这种标准的一个例子是把一个表面fM与前一个表面fM-1进行比较。这种比较的一种方式是找出表面fM和表面fM-1的线性组合因子dk之差。如果,该差的平方和小于所规定的极限值,那么,该序列被认为是收敛于该表面fM。
设计者可能会发现:对于一个解的收敛速度太慢,例如:表面fM和前表面fM-1之间的差太大。根据本发明的一个优选实施例,在转入最初的问题之前对于迭代过程中的一步或多步,表面序列的收敛性可以通过把由方程6所给出的成本函数E用变化成本函数E′来代替而增加。该变化成本函数E′是含有分别为方程1C和1D的主曲率k1和k2的方程6的成本函数,只是其中的k1和k2被方程9C和9D中所给出的k1′和k2′所代替:
将会看到,当β的值为1时,方程9A—9D与方程1A—1D相同,因而,该变化的成本函数E′与方程6所给出的成本函数相同。当β的值为0时,该变化的成本函数E′为方程6所给出的成本函数的线性化。如前所述,在转入初始问题之前,对于迭代过程中的一步或几步,该表面序列的收敛性可以通过利用β的值处于0和1之间,包括0而除去1的变化成本函数E′而增加。也可以看到,当β的值大于1时,k1′和k2′将失去它们的几何意义。
可以看到,以上所述的表面计算可以用于计算任何一种光学元件的单个未知表面,其中该光学元件可能包含多个折射面。
根据本发明的一个优选实施例,在图1的步骤100中指定的参数、步骤108中修正的参数和用于步骤102计算中的参数是把区域D划分为三角形τi、适宜的光焦度函数P(x,y)、加权因子wi,1和wi,2、曲线Cm和曲线限制。
另外,如上所述的有关图1中的步骤104,一旦透镜表面f计算得出,设计者就可以检查所计算得到的该表面的表面特性和模拟整个透镜的光学特性。检查的表面特性包括表面光焦度和表面像散性。例如:所计算得到的表面的表面像散性可能变化太快,会使得基于所计算得到的透镜的配戴者感到不适。
根据本发明的一个优选实施例,这种模拟是模拟当从人眼观看时该透镜的光学特性。一个或多个波前通过该透镜。在本发明的一个优选实施例中,用一组球面波通过该透镜。通常,当波前通过该透镜时它们发生变形。在该透镜的另一边所得到的变形波前可以由包括像散性、光焦度、光谱等多种特性来反映。变形波前的计算可以由公知的光线跟踪技术完成。在D.Malacara和Z.Malacara发表在1992年的“光学工程”(Optical Engineering)的第31卷第7章第1551-1555页,题为“Testing and centering of lenses by means of a Hartmann test withfour holes”一文中描述了一种确定变形波前特性的技术。
本领域的技术人员将看到,本发明并不限制在以上所示的特例和描述中。本发明的范围而是由权利要求书来限定。
Claims (14)
1、一种设计多焦光学元件的表面的方法,该方法包括以下步骤:
把一个区域划分为多个三角形;
定义一组函数,其中每个函数表示一个所述三角形上所述表面的一部分;和
优化所述函数,由此生成所述表面。
2、一种如权利要求1所述的方法,其中所述的光学元件是渐变眼镜。
3、一种如权利要求1所述的方法,还包括以下步骤:
模拟基于所生成的表面的所述光学元件的光学特性;和
如果对所生成的表面不满意,调整所述函数的参数值,重复所述的优化和模拟步骤。
4、一种设计多焦光学元件表面的方法,该方法包括以下步骤:
把一个区域划分为多个多边形;
定义一组函数,其中每个函数表示一个所述多边形上所述表面的一部分;
为所述函数值和所述函数在所述多边形边界处的一阶导数值设定连续性限制;
对于所述函数在所述多边形边界处的二阶导数值不需要连续性限制;和
优化受所设连续性限制的所述函数,由此生成所述表面。
5、一种如权利要求4所述的方法,其中所述的光学元件为渐变眼镜。
6、一种如权利要求4所述的方法,还包括以下步骤:
模拟基于所生成的表面的所述光学元件的光学特性;和
如果对所生成的表面不满意,调整所述函数中的参数值,重复所述的优化和模拟步骤。
7、一种如权利要求4所述的方法,其中所述的优化步骤还包括以下步骤:
定义所述函数组的一个函数;和
优化所述函数。
8、一种设计多焦光学元件表面的方法,该方法包括以下步骤:
把一个区域划分为多个多边形;
定义一组函数,其中每个函数表示一个所述多边形上所述表面的一部分;
确定一组沿所述多边形边界子集的曲线;
在与所述曲线组相应的多边形的选定顶点处,指定至少一组含有所述函数、所述函数的一阶导数和所述函数的二阶导数或它们的任何组合的值;和
优化所述函数以满足所指定的值,由此生成所述表面。
9、一种如权利要求8所述的方法,其中所述的光学元件为渐变眼镜。
10、一种如权利要求8所述的方法,还包括以下步骤:
模拟基于所生成表面的所述光学元件的光学特性;和
如果对所生成的表面不满意,调整所述函数中的参数值,重复所述的优化和模拟步骤。
11、一种如权利要求8所述的方法,其中所述优化步骤还包括以下步骤:
定义所述函数组的一个函数;和
优化所述函数。
12、一种具有多个表面的多焦光学元件,其中,所述表面中至少一个是根据以下的方法设计的,该方法包括以下步骤:
把一个区域划分为多个三角形;
定义一组函数,其中每个函数表示一个所述三角形上所述表面的一部分;和
优化所述函数,由此生成所述表面。
13、一种具有多个表面的多焦光学元件,其中,所述表面中至少一个是根据以下方法设计的,该方法包括以下步骤:
把一个区域划分为多个多边形;
定义一组函数,其中每个函数表示一个所述多边形上所述表面的一部分;
为所述函数值和所述函数在所述多边形边界处的一阶导数值设定连续性限制;
对于所述函数在所述多边形边界处的二阶导数值不需要连续性限制;和
优化受所设连续性限制的所述函数,由此生成所述表面。
14、一种具有多个表面的多焦光学元件,其中,所述表面中至少一个是根据以下的方法设计的,该方法包括以下步骤:
把一个区域划分为多个多边形;
定义一组函数,其中每个函数表示一个所述多边形上所述表面的一部分;
确定一组沿所述多边形边界子集的曲线;
在与所述曲线组相应的多边形的选定顶点处,指定至少一组含有所述函数、所述函数的一阶导数和所述函数的二阶导数或它们的任何组合的值;和
优化所述函数以满足所指定的值,由此生成所述表面。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/262,341 | 1999-03-04 | ||
US09/262,341 US6302540B1 (en) | 1999-03-04 | 1999-03-04 | Method for the design of multifocal optical elements |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1302386A CN1302386A (zh) | 2001-07-04 |
CN1132042C true CN1132042C (zh) | 2003-12-24 |
Family
ID=22997100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN008002320A Expired - Fee Related CN1132042C (zh) | 1999-03-04 | 2000-02-08 | 多焦光学元件的设计方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US6302540B1 (zh) |
EP (1) | EP1076840A4 (zh) |
JP (1) | JP2002538507A (zh) |
KR (1) | KR20010024989A (zh) |
CN (1) | CN1132042C (zh) |
AU (1) | AU763552B2 (zh) |
BR (1) | BR0005226A (zh) |
CA (1) | CA2329714A1 (zh) |
IL (1) | IL139188A (zh) |
WO (1) | WO2000052517A1 (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6366823B1 (en) * | 1998-07-30 | 2002-04-02 | Asahi Kogaku Kogyo Kabushiki Kaisha | Design method for optical curved surface |
AU6575301A (en) * | 2000-04-25 | 2001-11-07 | Rodenstock Optik G | Method for calculating a progressive spectacle lens and method for producing a spectacle lens of this type |
US6655803B1 (en) * | 2000-06-01 | 2003-12-02 | Inray Ltd. | Wavefront method for designing optical elements |
WO2002088828A1 (fr) * | 2001-04-26 | 2002-11-07 | Hoya Corporation | Procede de conception de verres de lunettes et verres de lunettes |
US6554427B1 (en) * | 2001-12-11 | 2003-04-29 | Johnson & Johnson Vision Care, Inc. | Methods for designing contact lenses |
US6755524B2 (en) * | 2001-12-12 | 2004-06-29 | Inray Ltd. | Ophthalmic optical elements and methods for the design thereof |
US6786600B2 (en) | 2002-05-01 | 2004-09-07 | Shamir Optical Industry | Methods for generating a progressive surface and for production of multifocal progressive lenses |
ATE524762T1 (de) * | 2002-05-31 | 2011-09-15 | Crossbows Optical Ltd | Gleitsichtglas |
WO2004034131A1 (ja) * | 2002-10-08 | 2004-04-22 | Hoya Corporation | 眼鏡レンズの光学値の決定方法、眼鏡レンズの製造方法、眼鏡レンズ及びその受発注システム |
US20050052615A1 (en) * | 2003-09-05 | 2005-03-10 | Regents Of The University Of Minnesota | Multifocal optical device design |
US7150527B2 (en) * | 2004-08-17 | 2006-12-19 | Eyal Optical Industry | Progressive power semi-finished spectacle lens blanks for astigmatic presbyopes |
EP1862110A1 (en) * | 2006-05-29 | 2007-12-05 | Essilor International (Compagnie Generale D'optique) | Method for optimizing eyeglass lenses |
EP2075622A3 (fr) | 2007-11-27 | 2010-01-13 | Optofinanz AG | Lentilles ophtalmiques multifocales progressives à distribution constante |
EP2255243A2 (en) * | 2008-02-20 | 2010-12-01 | Swissphonics SA | Progressive multifocal ophthalmic lenses identically optimized over a wide range of refractive indices, base curves and additions |
EP2297661B1 (en) * | 2008-04-23 | 2018-06-13 | The Cleveland Clinic Foundation | Method for modeling biomechanical properties of an eye |
KR101515197B1 (ko) * | 2008-11-10 | 2015-04-24 | 삼성전자주식회사 | 초점 심도 특성으로부터 물리적 렌즈를 설계하는 방법 및 그를 이용해 설계된 연장된 초점 심도를 갖는 축대칭 렌즈 |
EP2207118A1 (en) * | 2008-12-31 | 2010-07-14 | Essilor International (Compagnie Générale D'Optique) | Method for calculating a system, for example an optical system |
US9354455B2 (en) * | 2009-01-23 | 2016-05-31 | Rodenstock Gmbh | Controlling designs using a polygonal design |
EP2270577A1 (en) * | 2009-06-30 | 2011-01-05 | Essilor International (Compagnie Générale D'Optique) | Method of and apparatus for generating a surface of an optical lens |
DE102011009473B4 (de) * | 2010-04-28 | 2022-03-17 | Rodenstock Gmbh | Computerimplementiertes Verfahren zur Berechnung eines Brillenglases mit blickwinkelabhängigen Verordnungsdaten, Vorrichtung zum Berechnen oder Optimieren eines Brillenglases, Computerprogrammerzeugnis, Speichermedium, Verfahren zum Herstellen eines Brillenglases sowie Verwendung eines Brillenglases |
EP2490065A1 (en) * | 2011-02-18 | 2012-08-22 | ESSILOR INTERNATIONAL (Compagnie Générale d'Optique) | A method for determining target optical functions |
US10473952B2 (en) * | 2013-06-13 | 2019-11-12 | Essilor International | Method for optimizing an optical surface |
AU2014280078B2 (en) * | 2013-06-13 | 2017-12-14 | Essilor International | Method for optimizing a measured contour of a spectacle frame |
US10734639B2 (en) | 2013-07-03 | 2020-08-04 | California Institute Of Technology | Carbon nanotubes—graphene hybrid structures for separator free silicon—sulfur batteries |
EP3315948B1 (de) * | 2016-10-26 | 2019-09-04 | Heraeus Quarzglas GmbH & Co. KG | Verfahren zur ermittlung des brechzahlprofils eines zylinderförmigen optischen gegenstandes |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2013548A6 (es) * | 1989-06-12 | 1990-05-01 | Indo Int Sa | Perfeccionamientos en lentes progresivas. |
FR2683642B1 (fr) * | 1991-11-12 | 1994-01-14 | Essilor Internal Cie Gle Optique | Lentille ophtalmique multifocale progressive. |
US5452031A (en) * | 1993-05-05 | 1995-09-19 | Boston Eye Technology, Inc. | Contact lens and a method for manufacturing contact lens |
US5561754A (en) * | 1993-08-17 | 1996-10-01 | Iowa State University Research Foundation, Inc. | Area preserving transformation system for press forming blank development |
US6567576B2 (en) | 2001-02-05 | 2003-05-20 | Jds Uniphase Inc. | Optical switch matrix with failure protection |
US6567323B2 (en) | 2001-02-27 | 2003-05-20 | Texas Instruments Incorporated | Memory circuit redundancy control |
-
1999
- 1999-03-04 US US09/262,341 patent/US6302540B1/en not_active Expired - Lifetime
-
2000
- 2000-02-08 BR BR0005226-4A patent/BR0005226A/pt not_active IP Right Cessation
- 2000-02-08 EP EP00902842A patent/EP1076840A4/en not_active Withdrawn
- 2000-02-08 CA CA002329714A patent/CA2329714A1/en not_active Abandoned
- 2000-02-08 JP JP2000602673A patent/JP2002538507A/ja active Pending
- 2000-02-08 CN CN008002320A patent/CN1132042C/zh not_active Expired - Fee Related
- 2000-02-08 AU AU24580/00A patent/AU763552B2/en not_active Ceased
- 2000-02-08 IL IL13918800A patent/IL139188A/en not_active IP Right Cessation
- 2000-02-08 KR KR1020007012230A patent/KR20010024989A/ko not_active Application Discontinuation
- 2000-02-08 WO PCT/IL2000/000082 patent/WO2000052517A1/en not_active Application Discontinuation
Also Published As
Publication number | Publication date |
---|---|
KR20010024989A (ko) | 2001-03-26 |
AU2458000A (en) | 2000-09-21 |
CA2329714A1 (en) | 2000-09-08 |
JP2002538507A (ja) | 2002-11-12 |
IL139188A (en) | 2004-05-12 |
BR0005226A (pt) | 2001-01-02 |
EP1076840A1 (en) | 2001-02-21 |
EP1076840A4 (en) | 2004-03-03 |
CN1302386A (zh) | 2001-07-04 |
US6302540B1 (en) | 2001-10-16 |
WO2000052517A1 (en) | 2000-09-08 |
IL139188A0 (en) | 2001-11-25 |
AU763552B2 (en) | 2003-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1132042C (zh) | 多焦光学元件的设计方法 | |
US9104041B2 (en) | Optical system determination according to advanced criteria | |
US6382789B1 (en) | Toric ophthalmic lenses | |
EP2384479B1 (en) | Method for calculating an optical system | |
US9547183B2 (en) | Method for determining an ophthalmic lens | |
US9523864B2 (en) | Method for determining a progressive ophthalmic lens | |
US20080013038A1 (en) | Method for defining a supplementary face for spectacles | |
EP2676165A1 (en) | A method for determining target optical functions | |
US8608312B2 (en) | Method for determining an aspherization layer for an ophthalmic lens | |
WO2006064384A2 (en) | Spectacle lenses incorporating atoric surfaces | |
CN1268964C (zh) | 一对多焦渐变眼镜片 | |
Xiang et al. | Comparison of different power distributions for designing freeform progressive addition lenses based on a minimization error function model | |
CN1854787A (zh) | 非球面镜片的设计方法 | |
CN112946922B (zh) | 具有矫正散光的自由曲面渐进式镜片的设计方法及镜片 | |
Xia et al. | Surface astigmatism correction using segmented freeform surfaces for a progressive addition lens | |
US6974214B2 (en) | B-spline optimization of optical glasses especially spectacle glasses with planes of symmetry | |
CN116057456A (zh) | 一种用于提供成品单光眼科镜片的计算机实现的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |