CN103869395A - 一种对称多边形压弯反射镜 - Google Patents

一种对称多边形压弯反射镜 Download PDF

Info

Publication number
CN103869395A
CN103869395A CN201410096401.1A CN201410096401A CN103869395A CN 103869395 A CN103869395 A CN 103869395A CN 201410096401 A CN201410096401 A CN 201410096401A CN 103869395 A CN103869395 A CN 103869395A
Authority
CN
China
Prior art keywords
catoptron
value
summit
bending
symmetry
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
Application number
CN201410096401.1A
Other languages
English (en)
Other versions
CN103869395B (zh
Inventor
李明
张瑶
盛伟繁
刘鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of High Energy Physics of CAS
Original Assignee
Institute of High Energy Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of High Energy Physics of CAS filed Critical Institute of High Energy Physics of CAS
Priority to CN201410096401.1A priority Critical patent/CN103869395B/zh
Publication of CN103869395A publication Critical patent/CN103869395A/zh
Application granted granted Critical
Publication of CN103869395B publication Critical patent/CN103869395B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mechanical Optical Scanning Systems (AREA)
  • Optical Elements Other Than Lenses (AREA)

Abstract

本发明公开了一种对称多边形压弯反射镜,用于将光源发射的光反射聚焦、发散或准直,所述反射镜为等厚的片体,且所述片体的反射面为边数大于4的轴对称多边形,对所述片体施以平行于反射面并垂直于对称轴的弯矩,使所述反射镜弹性压弯后的镜面面形趋近于一条二次曲线。本发明采用多边形压弯反射镜,压弯面形趋近于一条理想二次曲线,该多边形可根据实际需要调整,从而可在压弯时使得面形精度较高,而且多边形压弯反射镜可通过传统工艺加工成型。

Description

一种对称多边形压弯反射镜
技术领域
本发明涉及反射镜,具体涉及多边形压弯反射镜。
背景技术
各种光学工程尤其是同步辐射X射线多采用反射镜实现聚焦、发散或准直等功能,因此也被称为聚焦镜、准直镜等等。从理论上讲,采用曲边压弯镜作为反射镜可以获得无面形误差的理想面形,但是曲边镜的镜体加工工艺难度很大,成本很高,周期很长。尤其对于较长的镜体(如1m长),更是无法进行加工。因此,由于加工难度、成本等原因,目前市场上还没有曲边镜出现。目前,限于镜体加工的工艺难度,有人采用了梯形压弯镜(或矩形压弯镜)。梯形压弯镜和矩形压弯镜虽然可通过传统工艺完成加工,但是其压弯的面形精度较低,对于要求较高的光学系统(比如低发射度的第三代同步辐射光源的高品质相空间调制)则无法满足聚焦、发散或准直等要求。
发明内容
针对现有技术中存在的问题,本发明的目的为提供一种易于加工且压弯面形精度较高的对称多边形压弯反射镜。
为实现上述目的,本发明的技术方案如下:
一种对称多边形压弯反射镜,用于将光源发射的光反射聚焦、发散或准直,所述反射镜为等厚的片体,且所述片体的反射面为边数大于4的轴对称多边形,对所述片体施以平行于反射面并垂直于对称轴的弯矩,使所述反射镜弹性压弯后的镜面面形趋近于一条二次曲线,该二次曲线的方程为:
com ( x ) ≡ con ( p , q , θ ; x ) = ( p + q ) ( ( p - q ) x cos θ + 2 ( - pq + pq ( pq - x 2 - px cos θ + qx cos θ ) ) ) sin θ - ( p + q ) 2 + ( p - q ) 2 sin 2 θ
其中,p为物距、q为像距、θ为镜面中心处光线掠入射角、x为沿镜面对称轴方向,以对称轴中心为原点的坐标值,当pq>0时为椭圆,pq<0时为双曲线,p或q趋向于±∞时,该式极限为抛物线。
进一步,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure BDA0000477490350000021
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标
Figure BDA0000477490350000022
其余顶点xi,i=2…,n-1,可在
Figure BDA0000477490350000023
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure BDA0000477490350000024
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)采用材料力学梁弯曲的近似理论得到所述反射镜的压弯面形与所述二次曲线的面形之间面形斜率误差均方根值,该值为:
&sigma; SlopeError ( w i , i = 1 . . . n , M 0 f , k Mf ) = ( 1 L &times; &Integral; - L 2 L 2 ( &Integral; 0 x ( M 0 f ( 1 + k Mf x &prime; ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x &prime; ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; x &prime; ) ) d x &prime; ) 2 dx ) 1 / 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(3)根据所述步骤(2)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σStopeError(wi,i=1...n,M0f,kMf'的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
进一步,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure BDA0000477490350000027
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标
Figure BDA0000477490350000028
其余顶点xi,i=2…,n-1,可在
Figure BDA0000477490350000029
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure BDA00004774903500000210
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)采用材料力学梁弯曲的近似理论得到所述反射镜的压弯面形与所述二次曲线的面形之间面形曲率误差均方根值,该值为:
&sigma; CurvatureError ( w i , i = 1 . . . n , M 0 f , k Mf ) = ( 1 L &times; &Integral; - L 2 L 2 ( M 0 f ( 1 + k Mf x ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; ) ) 2 dx ) 1 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(3)根据所述步骤(2)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σCurvature Error(wi,i=1...n,M0f,kMf'的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
进一步,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure BDA0000477490350000033
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标其余顶点xi,i=2…,n-1,可在范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure BDA0000477490350000036
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)计算理想宽度分布函数:
W ( x ) = 12 M 0 f ( 1 + k Mf x ) ET 3 con &prime; &prime; ( p d , q d , &theta; d ; x )
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,M0f和kMf可取使镜体符合梁近似的任意值,由此确定两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值;
(3)计算反射镜宽度分布与上述理想宽度分布偏差的均方根值,该值为:
&sigma; widthError ( w i , i = 1 . . . n ) = ( 1 L &times; &Integral; - L 2 L 2 ( W d ( x i , w i , i = 1 . . . n ; x ) - W ( x ) ) 2 dx ) 1 / 2
其中,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(4)根据所述步骤(3)中的公式,计算σwidthError(wi,i=1...n)的最小值,得出最优的反射镜各顶点位置。
进一步,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure BDA0000477490350000043
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标其余顶点xi,i=2…,n-1,可在
Figure BDA0000477490350000045
范围任意由小到大取为互异的值,也可简单取为等间距值即i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)任意设定具备实际意义的wi、M0f和kMf的初始值,然后采用有限单元法结构分析,以压弯面形y(x)相对于所述二次曲线的面形m阶导数误差均方根值 &sigma; FigureError = ( 1 L &times; &Integral; - L 2 L 2 ( y ( m ) ( x ) - com ( m ) ( p d , q d , &theta; d ; x ) ) 2 dx ) 1 / 2 作为目标函数,以wi、M0f和kMf共n+2个参数中任意n个参数为变量,进行最小化的优化计算,其中上标(m)表示该函数对x的m阶导数,m取为非负整数,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值;
(3)根据所述步骤(1)、(2)中的计算得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
本发明与现有技术相比,本发明采用多边形压弯反射镜,压弯面形趋近于理想二次曲线,该多边形设计可根据实际需要调整,从而可在压弯时使得面形精度较高,而且多边形压弯反射镜可通过传统工艺加工成型,简单方便,因此本发明的多边形压弯反射镜更加便于市场化。
附图说明
下面结合附图对本发明作进一步详细说明:
图1为本发明的对称多边形压弯反射镜的使用状态示意图;
图2为本发明的对称多边形压弯反射镜的结构及压弯示意图;
图3为本发明的对称多边形压弯反射镜具体实施例中采用第一种方法和第二种方法时面形曲率误差示意图;
图4为本发明的对称多边形压弯反射镜具体实施例中采用第一种方法和第二种方法得到的反射镜的结构示意图;
图5为本发明的对称多边形压弯反射镜具体实施例中采用第三种方法时面形斜率误差示意图;
图6为本发明的对称多边形压弯反射镜具体实施例中采用第三种方法得到的反射镜的结构示意图;
图7为本发明的对称多边形压弯反射镜具体实施例中采用第四种方法时面形斜率误差示意图;
图8为本发明的对称多边形压弯反射镜具体实施例中采用第四种方法得到的反射镜的结构示意图。
具体实施方式
体现本发明特征与优点的典型实施例将在以下的说明中详细叙述。应理解的是本发明能够在不同的实施例上具有各种的变化,其皆不脱离本发明的范围,且其中的说明及附图在本质上是当作说明之用,而非用以限制本发明。
本发明的对称多边形压弯反射镜一种多边形压弯反射镜,用于将一光源发射的光反射聚焦,如图1所示,其聚焦公式为:
1 p + 1 q + 1 f - - - ( 1 )
其中,p为物距,q为像距,f为反射镜的焦距。镜面中心处光线掠入射角为θ,镜面上最大掠入射角θmax在镜子一端位置,在掠入射角小于临界掠入射角时可发生全反射,临界掠入射角较小,一般为mrad至百mrad水平,由所反射X射线波长及镜面材料决定,因此X射线反射镜多为长条形。我们以镜子中心为原点,沿镜面对称轴方向为x轴,沿垂直镜面方向为y轴建立xyz右手直角三维坐标系,并在xy平面内得到物点坐标为(-pcosθ,psinθ),像点坐标为(qcosθ,qsinθ)。通过镜面上某点(x,y)的光线光程表达式为:
s = ( x + p cos &theta; ) 2 + ( y - p sin &theta; ) 2 + ( x - q cos &theta; ) 2 + ( y - q sin &theta; ) 2 - - - ( 2 )
依费马原理,光程最短,s对x的全微分为0,得到理想镜面的二次曲线方程:
com ( x ) &equiv; con ( p , q , &theta; ; x ) = ( p + q ) ( ( p - q ) x cos &theta; + 2 ( - pq + pq ( pq - x 2 - px cos &theta; + qx cos &theta; ) ) ) sin &theta; - ( p + q ) 2 + ( p - q ) 2 sin 2 &theta; - - - ( 3 )
其中,p为物距、q为像距、θ为镜面中心处光线掠入射角、x为沿镜面对称轴方向以对称轴中心为原点的坐标值。当pq>0时为椭圆,pq<0时为双曲线,p或q趋向于±∞时,该式极限为抛物线。
本发明的反射镜为等厚的片体,且所述片体的反射面为边数大于4的轴对称多边形,并对片体施以平行于反射面并垂直于对称轴的弯矩(如图2所示)。本发明的反射镜弹性压弯后的镜面面形近似于式(3)所示出二次曲线。
本发明中,设定反射镜多边形的对称两侧顶点数各为n(包含落在对称轴上的顶点),每对对称的顶点间距为wi,位置坐标为
Figure BDA0000477490350000063
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf
上述各对顶点x坐标xi,i=1,2…,n+1如下取值:
两端顶点坐标
Figure BDA0000477490350000064
其余顶点xi,i=2…,n-1,可在
Figure BDA0000477490350000065
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure BDA0000477490350000066
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
可得镜体宽度函数Wd(x)为n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
上述各参数wi、M0f和kMf通过下述过程获得:
第一种方法:
(11)采用材料力学梁弯曲的近似理论得到反射镜的压弯挠度微分方程:
y &prime; &prime; ( x ) = M ( x ) EI ( x ) - - - ( 4 )
其中,M(x)为x位置的总弯矩,仅考虑两端弯矩时M(x)=M0f(1+kMfx);I(x)为惯性矩;E为杨氏模量。镜子x处惯性矩I(x)的物理定义式为
Figure BDA0000477490350000073
其中,W(x)为镜子在x处的宽度,T(x)为镜子在x处的厚度。对于等厚度镜子,T(x)=T为常数,则W(x)与I(x)成正比。
(12)根据步骤(11)得出压弯面形,从而获得该压弯面形与所述二次曲线的面形之间面形斜率误差均方根值,该值为:
&sigma; SlopeError ( w i , i = 1 . . . n , M 0 f , k Mf ) = ( 1 L &times; &Integral; - L 2 L 2 ( &Integral; 0 x ( M 0 f ( 1 + k Mf x &prime; ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x &prime; ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; x &prime; ) ) d x &prime; ) 2 dx ) 1 / 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值;
(13)根据步骤(12)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σStopeError(wi,i=1...n,M0f,kMf)的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
第二种方法:
(21)采用材料力学梁弯曲的近似理论得到反射镜的压弯挠度微分方程:
Figure BDA0000477490350000072
其中,M(x)为x位置的总弯矩,仅考虑两端弯矩时M(x)=M0f(1+kMfx);I(x)为惯性矩;E为杨氏模量。镜子x处惯性矩I(x)的物理定义式为
Figure BDA0000477490350000074
其中,W(x)为镜子在x处的宽度,T(x)为镜子在x处的厚度。对于等厚度镜子,T(x)=T为常数,则W(x)与I(x)成正比。
(22)根据步骤(11)得出压弯面形,从而获得该压弯面形与所述二次曲线的面形之间面形曲率误差均方根值,该值为:
σCUrvature Error(wi,i=1...n,M0f,kMf)
= ( 1 L &times; &Integral; - L 2 L 2 ( M 0 f ( 1 + k Mf x ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; x ) ) 2 dx ) 1 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值;
(23)根据步骤(12)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σCurvature Error(wi,i=i...n,M0f,kMf)的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
第三种方法:
(31)首先计算理想宽度分布函数:
W ( x ) = 12 M 0 f ( 1 + k Mf x ) ET 3 con &prime; &prime; ( p d , q d , &theta; d ; x ) - - - ( 5 )
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,M0f和kMf可取使镜体符合梁近似的任意值;将(5)式代入(4)式可知,该曲边镜压弯面形为理想二次曲线。
(32)计算反射镜宽度分布与上述理想宽度分布偏差的均方根值,该值为:
&sigma; widthError ( w i , i = 1 . . . n ) = ( 1 L &times; &Integral; - L 2 L 2 ( W d ( x i , w i , i = 1 . . . n ; x ) - W ( x ) ) 2 dx ) 1 / 2
(33)根据所述步骤(32)中的公式,计算σwidth Error(wi,i=1...n)的最小值,确定反射镜各顶点位置;将(4)式两端求倒数,并考虑小压弯量及等厚度T,可得压弯面形曲率半径:
&rho; ( x ) &cong; 1 y &prime; &prime; ( x ) = EI ( x ) M ( x ) &Proportional; W ( x )
可见,保持施加力矩M(x)不变,当反射镜宽度与理想宽度偏差均方根最小时,其压弯面形曲率半径误差达到最小值;因此,所确定反射镜各顶点位置即为最优,两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值保持与步骤(31)中理想宽度计算时相同。
第四种方法:
(41)任意设定具备实际意义的wi、M0f和kMf的初始值(可选择矩形镜、梯形镜或曲边镜的力学参数和宽度函数的采样值作为初始值),然后采用有限单元法结构分析,以压弯面形y(x)相对于所述二次曲线的面形m阶导数误差均方根值 &sigma; FigureError = ( 1 L &times; &Integral; - L 2 L 2 ( y ( m ) ( x ) - com ( m ) ( p d , q d , &theta; d ; x ) ) 2 dx ) 1 / 2 作为目标函数,以wi、M0f和kMf共n+2个参数中任意n个参数为变量,进行最小化的优化计算,其中上标(m)表示该函数对x的m阶导数,m取为非负整数,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值;
(42)根据所述步骤(41)中的计算得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
本发明与现有技术相比,本发明采用多边形压弯反射镜,压弯面形近似于理想二次曲线,该多边形可根据实际需要调整,从而可在压弯时使得面形精度较高,而且多边形压弯反射镜可通过传统工艺加工成型,简单方便,因此本发明的多边形压弯反射镜更加便于市场化。
下面以p=40.18m、q=0.12m、θ=0.003rad、L=0.12m、T=0.005m、E=1.301×1011Pa压弯聚焦镜(可得理想二次曲线为一个椭圆)为例:
根据第一种方法或第二种方法获得,n=9时,面形斜率误差如图3所示,其面形斜率误差均方根值为0.0664μrad。此压弯面形精度已高于当前国际上实现的最高精度。
镜子形状如图4所示,反射面为18边形,其9对对称顶点坐标为(-0.06,±0.003598),(-0.045,±0.004797),(-0.03,±0.005577),(-0.015,±0.005975),(0,±0.006026),(0.015,±0.005767),(0.03,±0.005239),(0.045,±0.004489),(0.06,±0.003572)(单位均为m)。镜子中心弯矩M0f=0.2039N·m,弯矩分布的相对斜率kMf=11.27m-1,相应x=-0.06m端弯矩为0.06608N·m,x=0.06m端弯矩为0.3418N·m。
根据第三种方法获得,n=9时,面形斜率误差如图5所示,其面形斜率误差均方根值为0.0648μrad。此压弯面形精度已高于当前国际上实现的最高精度。
镜子形状如图6所示,反射面为18边形,其9对对称顶点坐标为(-0.06,±0.003589),(-0.045,±0.004781),(-0.03,±0.00556),(-0.015,±0.005957),(0,±0.006007),(0.015,±0.005749),(0.03,±0.005223),(0.045,±0.004475),(0.06,±0.003562)(单位均为m)。镜子中心弯矩M0f=0.2033N·m,弯矩分布的相对斜率kMf=11.27m-1,相应x=-0.06m端弯矩为0.06587N·m,x=0.06m端弯矩为0.3408N·m。
根据第四种方法获得,n=9时,面形斜率误差如图7所示,其面形斜率误差均方根值为0.0415μrad。此压弯面形精度已高于当前国际上实现的最高精度。
镜子形状如图8所示,反射面为18边形,其9对对称顶点坐标为(-0.06,±0.003637),(-0.045,±0.004809),(-0.03,±0.005575),(-0.015,±0.005965),(0,±0.006012),(0.015,±0.005753),(0.03,±0.005229),(0.045,±0.004484),(0.06,±0.003572)(单位均为m)。镜子中心弯矩M0f=0.2033N·m,弯矩分布的相对斜率kMf=11.27m-1,相应x=-0.06m端弯矩为0.06587N·m,x=0.06m端弯矩为0.3408N·m。
综上所述,本发明提供的多边形压弯反射镜具有与曲边镜体几乎等同的压弯面形精度,同时无需高难度的镜体加工工艺即可实现不同尺寸镜体的加工,大幅降低了加工难度和成本。
本发明的技术方案已由优选实施例揭示如上。本领域技术人员应当意识到在不脱离本发明所附的权利要求所揭示的本发明的范围和精神的情况下所作的更动与润饰,均属本发明的权利要求的保护范围之内。

Claims (5)

1.一种对称多边形压弯反射镜,用于将光源发射的光反射聚焦、发散或准直,其特征在于,所述反射镜为等厚的片体,且所述片体的反射面为边数大于4的轴对称多边形,对所述片体施以平行于反射面并垂直于对称轴的弯矩,使所述反射镜弹性压弯后的镜面面形趋近于一条二次曲线,该二次曲线的方程为:
com ( x ) &equiv; con ( p , q , &theta; ; x ) = ( p + q ) ( ( p - q ) x cos &theta; + 2 ( - pq + pq ( pq - x 2 - px cos &theta; + qx cos &theta; ) ) ) sin &theta; - ( p + q ) 2 + ( p - q ) 2 sin 2 &theta;
其中,p为物距、q为像距、θ为镜面中心处光线掠入射角、x为沿镜面对称轴方向以对称轴中心为原点的坐标值,当pq>0时为椭圆,pq<0时为双曲线,p或q趋向于±∞时,该式极限为抛物线。
2.如权利要求1所述的对称多边形压弯反射镜,其特征在于,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure FDA0000477490340000012
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标其余顶点xi,i=2…,n-1,可在
Figure FDA0000477490340000013
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure FDA0000477490340000014
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)采用材料力学梁弯曲的近似理论得到所述反射镜的压弯面形与所述二次曲线的面形之间面形斜率误差均方根值,该值为:
&sigma; SlopeError ( w i , i = 1 . . . n , M 0 f , k Mf ) = ( 1 L &times; &Integral; - L 2 L 2 ( &Integral; 0 x ( M 0 f ( 1 + k Mf x &prime; ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x &prime; ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; x &prime; ) ) d x &prime; ) 2 dx ) 1 / 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(3)根据所述步骤(2)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σSlopeBrror(wi,i=1...n,M0f,kMf)的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
3.如权利要求1所述的对称多边形压弯反射镜,其特征在于,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure FDA0000477490340000021
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标
Figure FDA0000477490340000022
其余顶点xi,i=2…,n-1,可在
Figure FDA0000477490340000023
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure FDA0000477490340000024
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)采用材料力学梁弯曲的近似理论得到所述反射镜的压弯面形与所述二次曲线的面形之间面形曲率误差均方根值,该值为:
&sigma; CurvatureError ( w i , i = 1 . . . n , M 0 f , k Mf ) = ( 1 L &times; &Integral; - L 2 L 2 ( M 0 f ( 1 + k Mf x ) E &times; T 3 W d ( x i , w i , i = 1 . . . n ; x ) 12 - con &prime; &prime; ( p d , q d , &theta; d ; ) ) 2 dx ) 1 2
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(3)根据所述步骤(2)中的公式,任意选定wi、M0f和kMf共n+2个参数中的2个赋予具备实际意义的值,计算σCurvature Error(wi,i=1...n,M0f,kMf)的最小值,得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
4.如权利要求1所述的对称多边形压弯反射镜,其特征在于,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure FDA0000477490340000027
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标
Figure FDA0000477490340000031
其余顶点xi,i=2…,n-1,可在
Figure FDA0000477490340000032
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure FDA0000477490340000033
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)计算理想宽度分布函数:
W ( x ) = 12 M 0 f ( 1 + k Mf x ) ET 3 con &prime; &prime; ( p d , q d , &theta; d ; x )
其中,E为镜体材料杨氏模量,T为镜体厚度,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值,M0f和kMf可取使镜体符合梁近似的任意值,由此确定两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值;
(3)计算反射镜宽度分布与上述理想宽度分布偏差的均方根值,该值为:
&sigma; widthError ( w i , i = 1 . . . n ) = ( 1 L &times; &Integral; - L 2 L 2 ( W d ( x i , w i , i = 1 . . . n ; x ) - W ( x ) ) 2 dx ) 1 / 2
其中,Wd为镜体宽度函数即n个点(xi,wi)的线性插值函数:
W d ( x ) &equiv; W d ( x i , w i , i = 1 . . . n ; x ) = w i + x - x i x i + 1 - x i ( w i + 1 - w i ) , x i &le; x &le; x i + 1 , i = 1,2 . . . , n - 1 ;
(4)根据所述步骤(3)中的公式,计算σwidth Error(wi,i=1...n)的最小值,得出最优的反射镜各顶点位置。
5.如权利要求1所述的对称多边形压弯反射镜,其特征在于,所述对称多边形中任一侧顶点和落在对称轴上的顶点的总数为n,每对对称的顶点间距为wi,位置坐标为
Figure FDA0000477490340000037
i=1,2…,n,在对称轴两端施加在所述反射镜中心处的弯矩值为M0f,两端施加在镜面上弯矩分布的相对斜率值为kMf;上述各参数xi、wi、M0f和kMf通过下述过程获得:
(1)上述各对顶点x坐标xi,i=1,2…,n如下取值:
两端顶点坐标
Figure FDA0000477490340000038
其余顶点xi,i=2…,n-1,可在
Figure FDA0000477490340000039
范围任意由小到大取为互异的值,也可简单取为等间距值即
Figure FDA00004774903400000310
i=1,2…,n,其中L为反射镜沿对称轴方向的长度;
(2)任意设定具备实际意义的wi、M0f和kMf的初始值,然后采用有限单元法结构分析,以压弯面形y(x)相对于所述二次曲线的面形m阶导数误差均方根值 &sigma; FigureError = ( 1 L &times; &Integral; - L 2 L 2 ( y ( m ) ( x ) - com ( m ) ( p d , q d , &theta; d ; x ) ) 2 dx ) 1 / 2 作为目标函数,以wi、M0f和kMf共n+2个参数中任意n个参数为变量,进行最小化的优化计算,其中上标(m)表示该函数对x的m阶导数,m取为非负整数,下标d代表设计值,pd、qd、θd依次表示物距、像距和镜面中心处掠入射角的设计值;
(3)根据所述步骤(1)、(2)中的计算得出最优的反射镜各顶点位置、两端施加在所述反射镜中心处的弯矩值以及两端施加在镜面上弯矩分布的相对斜率值。
CN201410096401.1A 2014-03-14 2014-03-14 一种多边形压弯反射镜 Expired - Fee Related CN103869395B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410096401.1A CN103869395B (zh) 2014-03-14 2014-03-14 一种多边形压弯反射镜

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410096401.1A CN103869395B (zh) 2014-03-14 2014-03-14 一种多边形压弯反射镜

Publications (2)

Publication Number Publication Date
CN103869395A true CN103869395A (zh) 2014-06-18
CN103869395B CN103869395B (zh) 2017-01-18

Family

ID=50908126

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410096401.1A Expired - Fee Related CN103869395B (zh) 2014-03-14 2014-03-14 一种多边形压弯反射镜

Country Status (1)

Country Link
CN (1) CN103869395B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62502779A (ja) * 1985-04-30 1987-10-22 インスチツ−ト オブスチエイ フイジキ アカデミ− ナウク エスエスエステル 単色放射集束法およびこの方法を実施する光学位相素子
CN102053352A (zh) * 2010-12-09 2011-05-11 中国科学院高能物理研究所 具有最大适用范围的同步辐射曲边聚焦镜设计方法
CN102103257A (zh) * 2011-01-28 2011-06-22 中国科学院高能物理研究所 同步辐射曲边聚焦镜重力协弯联合最大适用范围设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62502779A (ja) * 1985-04-30 1987-10-22 インスチツ−ト オブスチエイ フイジキ アカデミ− ナウク エスエスエステル 単色放射集束法およびこの方法を実施する光学位相素子
CN102053352A (zh) * 2010-12-09 2011-05-11 中国科学院高能物理研究所 具有最大适用范围的同步辐射曲边聚焦镜设计方法
CN102103257A (zh) * 2011-01-28 2011-06-22 中国科学院高能物理研究所 同步辐射曲边聚焦镜重力协弯联合最大适用范围设计方法

Also Published As

Publication number Publication date
CN103869395B (zh) 2017-01-18

Similar Documents

Publication Publication Date Title
CN105445918B (zh) 自由曲面离轴三反光学系统
CN107290845A (zh) 自由曲面离轴三反光学系统
JP2007528486A (ja) ビーム調整システム
Herkommer Phase space optics: an alternate approach to freeform optical systems
Oliker et al. Ray mapping and illumination control
JP3902796B2 (ja) 回折測距システム用可変ピッチ格子
JP2006235620A (ja) コリメーティングレンズ構造
JP6543825B2 (ja) マイクロレンズアレイ
JP2012194543A (ja) 回折光学素子及び計測装置
JP5942190B2 (ja) 二重反射型x線ミラーを用いた斜入射x線結像光学装置
Konyakhin et al. Optic-electronic sensor for measuring the deformations of the axle at the radio-telescope
CN103869395A (zh) 一种对称多边形压弯反射镜
CN108292026B (zh) 无热光学组件
CN103869394A (zh) 一种多边形压弯反射镜
CN104570305B (zh) 自由曲面透镜及含有该自由曲面透镜的成像系统
CN111029906A (zh) 一种激光器的矫正系统、光源系统及投影装置
Kallmann et al. Simulative and experimental characterization of an adaptive astigmatic membrane mirror
Vannoni et al. Adaptive x-ray mirror tuning simulation through influence functions’ modeling and error function minimization
CN104483725B (zh) 采用非球面单透镜作为准直结构的光栅曝光系统
CN209640610U (zh) 一种用于产生远距离稳定传输贝塞尔光束的装置
CN104898284A (zh) 一种非球面扩束镜
Völl et al. Free-form optics for non-idealized light sources in 3D: a phase-space approach
CN111316154A (zh) 光学元件以及激光照射装置
CN110058416B (zh) 一种用于产生远距离稳定传输贝塞尔光束的装置及方法
RU152500U1 (ru) Голографический коллиматорный прицел

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170118

Termination date: 20200314

CF01 Termination of patent right due to non-payment of annual fee