专利名称：光学元件表面波纹度对其激光损伤阈值影响的评价方法及由此获得元件最佳加工参数的方法
技术领域：
本发明涉及一种光学元件的表面质量评价方法以及一种获得元件最佳加工参数 的方法。
背景技术：
聚变能源清洁、无污染且几乎取之不尽，是未来解决能源问题的理想途径，目前， 各发达国家对其均给予高度重视。利用激光控制核聚变所需的高功率固体激光驱动器在最 后阶段照射核靶丸时需要提供很高的能量，以实现核聚变点火(点火所需能量为3 IOMJ/ cm2，3 5ns)。但是，由于目前激光驱动器中所采用的各类强光光学元件的激光损伤阈值 比较低(如KDP晶体实际阈值为12 20J/Cm2、lnS)，极大地限制了超大功率固体激光器的能 量输出，使核靶丸很难达到核聚变点火时所需能量要求。现阶段，对于光学元件激光损伤机 理的研究重点集中在材料内部稳态和瞬态缺陷(如位错、内部微裂纹等)、内部是否有杂质 及含量大小(如包裹体、有机物及各种杂质离子等)、材料内部发生的雪崩电离及其引致的 热效应等因素上。对这些影响因素进行合理的控制与消除能够在一定程度上提高元件的抗 激光损伤性，但结果依然远远小于元件的理论阈值(如KDP晶体理论阈值为14(T200 J/cm2、 Ins)。如何进一步提高此类光学元件的激光损伤阈值已成为聚变能否成功的关键技术。研 究表明，光学元件机械加工表面质量(如粗糙度、小尺度波纹等)同样会对其激光损伤阈值 有重要影响。美国“国家点火装置”中的大尺寸KDP晶体元件均采用目前代表世界最高加 工水平的“LLNL”实验室提供的超精密加工方法，以保证被加工元件具有良好的表面质量。 因此，深入了解诸如KDP晶体等强光光学元件机械加工表面质量对其激光损伤阈值的影响 机理并提出可靠的评价方法是目前急需解决的一个关键问题，它也是限制光学元件加工精 度的一个重要因素。
发明内容
本发明的目的是解决目前尚无用于评价光学元件表面波纹度对光学元件的激光 损伤阈值的影响程度的问题，以及由此导致用于获得元件最优加工参数的方法中存在结果 不准确的问题，提供了一种光学元件表面波纹度对其激光损伤阈值影响的评价方法及由此 获得元件最佳加工参数的方法。光学元件表面波纹度对其激光损伤阈值影响的评价方法，它的过程如下 步骤一、利用检测仪器，获取光学元件原始加工表面的形貌数据矩阵；
步骤二、根据步骤一获得的形貌数据矩阵，获得光学元件原始加工表面的功率谱密度 曲线，进而获得光学元件原始加工表面的各个特征频率以及每个特征频率的幅值；
步骤三、对步骤二获得的每个特征频率，采用二维连续小波变换法提取并再现各特征 频率的三维形貌，并利用傅立叶模方法计算每个特征频率对应的光学元件内部的光强分 布；步骤四、根据步骤三获得的每个特征频率对应的光学元件内部的光强分布，获得每个 特征频率对应的的光学元件内部的光强最大值，进而获得每个特征频率对应的相对激光损 伤阈值；
步骤五、对步骤四获得的每个特征频率对应的相对激光损伤阈值进行比较筛选，获得 所有相对激光损伤阈值中的最小值，并将所述最小值作为此次对光学元件评价的结果。由上述光学元件表面波纹度对其激光损伤阈值影响的评价方法获得元件最佳加 工参数的方法，它的过程如下
步骤A1、令
表示机床的加工参数组，其中Cg为加工参数的总数；获得所
述加工参数组
中每个加工参数的实际取值范围，其中，参数
c(i)的取值范围为
和1/2(0均为实数；
步骤A2、对每个加工参数,获得一个优选值，具体过程为
在每个加工参数C(Z)的取值范围[W1OX巧(0]内，选取g(￡)个等间距的点，在
其余加工参数均固定的情况下，分别在该g(O个点的条件下，对光学元件进行加工，获
得5(O个光学元件；利用光学元件表面波纹度对其激光损伤阈值影响的评价方法，获得
g(0个光学元件中每个光学元件的评价结果，即获得每个光学元件的相对激光损伤阈 值，然后通过比较筛选，获得相对激光损伤阈值最小的光学元件，并将该光学元件对应的加
工参数作为加工参数的优选值，其中
步骤A4、根据步骤二获得的每个加工参数的优选值,获得优选参数组
，所述优选参数组
即为待加工
元件的最佳加工参数组。本发明的有益效果本发明的评价方法，能够用于评价光学元件表面波纹度对光 学元件的激光损伤阈值的影响程度；本发明的获得元件最佳加工参数的方法，利用上述评 价方法，能够获得高加工质量的光学元件。
图1为小尺度波纹的物理模型图；图2为原始加工表面轮廓的功率谱密度曲线图； 图3为利用白光干涉仪获得的原始加工表面的三维形貌结果图；图4为图3中(SSTym)-1W 特征频率的三维形貌图；图5为ΟθβμπιΓ1的特征频率的三维形貌图；图6为(220 μ Hir1W 特征频率的三维形貌图；图7为(176 μ πι)—1的特征频率的三维形貌图；图8为^Ζ.δμπιΓ1的特征频率的三维形貌图；图9为相对激光损伤阈值随空间频率对应波长的变化曲线图； 图10为ΟΑμπιΓ1的特征频率对应的光强分布图；图11为(δδ μπιΓ1的特征频率对应的 光强分布图；图12为^Ζ.δμπιΓ1的特征频率对应的光强分布图；图13为(ΙΠμπιΓ1的 特征频率对应的光强分布图；图14为(176 μ m)—1的特征频率对应的光强分布图；图15为 (335 μ πι)—1的特征频率对应的光强分布图；图16和图17为KDP晶体的激光损伤阈值实验 损伤点的形貌图；图18为KDP晶体的理论激光损伤阈值与实验得到的相对激光损伤阈值的 对比图；图19为特征频率出现次数的统计图；图20为主导特征频率出现次数的统计图；图 21为
实验测得的波纹度与进给量、背吃刀量的关系图；图22为本发明的评价方法的流程图。
具体实施例方式具体实施方式
一本实施方式的光学元件表面波纹度对其激光损伤阈值影响的 评价方法，它的过程如下
步骤一、利用检测仪器，获取光学元件原始加工表面的形貌数据矩阵； 步骤二、根据步骤一获得的形貌数据矩阵，获得光学元件原始加工表面的功率谱密度 曲线，进而获得光学元件原始加工表面的各个特征频率以及每个特征频率的幅值；
步骤三、对步骤二获得的每个特征频率，采用二维连续小波变换法提取并再现各特征 频率的三维形貌，并利用傅立叶模方法计算每个特征频率对应的光学元件内部的光强分 布；
步骤四、根据步骤三获得的每个特征频率对应的光学元件内部的光强分布，获得每个 特征频率对应的的光学元件内部的光强最大值，进而获得每个特征频率对应的相对激光损 伤阈值；
步骤五、对步骤四获得的每个特征频率对应的相对激光损伤阈值进行比较筛选，获得 所有相对激光损伤阈值中的最小值，并将所述最小值作为此次对光学元件评价的结果。步骤一中所述检测仪器为白光干涉仪与原子力显微镜(AFM)，所述白光干涉仪的 型号为Taylorsurf CCI 2000，所述原子力显微镜采用美国DI公司生产的NanoscopeIII型 原子力显微镜。
具体实施方式
二 本实施方式是对实施方式一的光学元件表面波纹度对其激光损 伤阈值影响的评价方法的进一步说明，步骤二所述内容的具体过程为
令ζ Cr)表示步骤一获得的光学元件原始加工表面的形貌数据矩阵，其中Z Cr)中包含 了#个数据点，且每相邻两个数据点具有相同的采样间隔Δχ，整体采样长度为Δχ；
功率谱密度定义为波前各频率分量傅立叶频谱振幅的平方，是对光学元件 空间域上的表面轮廓函数作傅立叶变换在频率域上的结果，其一维的定义形式为
PSDO) = ㈧]2/Δι,,其中，r为空间频率，Δ κ为频率间隔为畸变波前的
Fourier振巾畐。实际采用如下公式获得光学元件原始加工表面的功率谱密度曲线
,)表示二维基本小波函数的尺度伸缩、坐标旋转及二维位移，
上式中j为波数j=2vi;，/;= /0V Δ χ)为空间频率， 为米样点的序数，且
由功率谱密度曲线，即获得光学元件原始加工表面的各个特征频率/y，再根据下式 计算获得每个特征频率/y的幅值
其中，」/为取样频率。理论上，在整个频谱范围内，功率谱密度曲线的所有极值点都可认为是特征频率， 本实施方式中实际选取峰值较大且较明显的那些极值点作为特征频率，它们在表面组成信 息中所占的比例也较大。
具体实施方式
三本实施方式是对实施方式一或二的光学元件表面波纹度对其激 光损伤阈值影响的评价方法的进一步说明，步骤三所述内容的具体过程为
二维连续小波变换(CWT2D，Continuous Wavelet Transform 2D)的一般形式为
其中，I= [J^xjr是平面直角坐标，f(x)表示二维信号，Wf(a,9yb)表示二维
连续小波变换， =Plf^f是& , Jt2方向上的位移，χ、表达式中的上标Γ表示转置，
为坐标旋转因子，为坐标系逆时针旋转角 a
φ[αΛ ·；1 (χ- )]为(χ- b)的共轭；
特征频率/；是与尺度a —一对应的，尺度与频率之间关系式为
其中，/；为所采用小波基函数的原始中心频率；△为测量仪器的采样周期；
度，
其中，/7=2,3，…，Λτ2;产表示第/7层一个周期内介质与空气分界面的坐标，第/7层的 实际介电常数为e Cr) ε0, f。为真空介电常数，第/7层的实际磁导率为// (χ)//。，//。为真 空磁导率，ε b为基体材料相对介电常数，μ b为基体材料相对磁导率； 将第/7层的相对介电常数和相对磁导率一起表示为傅立叶模形式为
公式Al
其中，Λ为傅里叶级数编号，εη为相对介电常数傅里叶展开后的第n项，μη为相对磁
导率傅里叶展开后的第n项， 由几何关系得
为第-层的 禾口；
对于Mexican2D小波，参考matlab小波工具箱，计算得到原始中心频率fc: 0. 25 ； 将所述原始中心频率/；、采样周期」以及待考察的特征频率/；代入上述关系式，即可 得到特征频率/；对应的尺度A ；再利用YAW小波工具箱，即可完成对所述形貌数据矩阵的 二维连续小波变换；
用正弦波来近似光学元件原始加工表面上的各频率的小尺度波纹，并建立小尺度波纹 的物理模型，所述小尺度波纹的物理模型位于χ-y-Z空间坐标系中，如图1所示，图中给出 了小尺度波纹2个周期的形貌，小尺度波纹的截面在χ -ζ平面，方向为小尺度波纹的波 纹线方向，小尺度波纹的基底平面垂直ζ轴、并沿χ轴方向周期变化，该周期为r;入射光波 以θ角入射到小尺度波纹表面、并通过小尺度波纹；用水平多分层形状对小尺度波纹的形 貌进行拟合，拟合的近似程度与分层数及剖分方法有关；
为便于计算，采用阶梯进行剖分处理，即沿ζ轴将所求空间分成^层，第1层为入射空 气层，第弋层为出射空气层，第Λτ 层为基底层，第2层至第怂-2层为小尺度波纹层，由此， 将整个小尺度波纹近场分布问题分解为求解分层的非均勻介质场的问题；
小尺度波纹层的相对介电常数ε Cr)及相对磁导率#⑴均具有周期性7，即 f Cr) = f {x+T)，μ Cr) =U {χ+Τ)，对于第/7 层有
其中Zp代表第P层上界面的ζ坐标，代表第/7-1层上界面的ζ坐标d代表小尺度 波纹幅值；
11其中，万即表示万(X)，"即表示"(χ), σ=批 ε0 ; ^ = 2双/Α, Λ为单位振幅的入射
光波长；、=O 0+Λ m/τ, a =Sin θ , θ为入射光与ζ轴的夹角， =0，士 1，士 2，…，士 I···， 为傅立叶模编号，I为计算时的截断常数# 、 、^;分别为电场的x、y和ζ向分量;力 、 hym、hzm分别为磁场的χ、y和ζ向分量；K表示第/7层波数的ζ分量，它是待求值； 每一分层中的电磁场满足麦克斯韦方程组
其中，i 是磁感应强度，々是电位移矢量；
考虑到小尺度波纹层在一个周期内、两种介质在χ方向上结合部的不连续性，利用傅 立叶因式分解“逆规则”原理，将公式Al和公式A2代入公式A3，可得TE波的本征方程
其中，j = μ~ι (ε- αμια) ^、只和#分别是根据傅立叶因式分解原理得到的系 数矩阵，上角标-1代表求逆运算；上述TE波的本征方程为广义本征方程，对其求解可以得 到由mi个本征矢组成的矩阵IFg和由正负各mi个本征值组成的对角阵卩+和卩，正
值Γ+代表上行波，负值Γ—代表下行波；
当各分层区域的本征模式场确定后，模式场的通解即为这些本征模式场的线性叠加，
对于第/7层，电场强度_7分量丑f和磁场强度χ分量Jif的最终表达式为
EJ = exp(ik0ax.) W^11' exp ikQr: (ζ-zp)] up + exp ik0F^\z-zv ]] j jcrfi了 = exp(ikQaχ)^-WFk j| exp ikQF^ (z-zp)] up - exp ikQFf (z_z J] dpj 其中if.cF为两个列矢量，由上行波的各本征模式场的振幅系数组成，Ζ由下行波的
由于小尺度波纹带来的介电常数与磁导率的周期性，使得电磁场空间分布也具有周期 性，即i Cr)=Z C^7)，//Cr)=//C^7)，其中i Cr)为电场强度，//⑴为磁场强度，因此，只需在一 个周期内讨论电场与磁场的分布情况；
第/7层电磁场一起表示为傅立叶模形式为 公式A2 该评价方法基于严格的电磁场理论，不具有实验破坏性，不会影响后续元件的正 常使用，评价方法原理简单、速度快，结果准确可靠。
具体实施方式
四本实施方式是对实施方式三的光学元件表面波纹度对其激光损 伤阈值影响的评价方法的进一步说明，二维基本小波函数一‘)采用MeXiCan2D小波基函 数，其表达式为
上式中， 和厶分别表示频域面坐标。
具体实施方式
五本实施方式是对实施方式一至四所述的任意一种光学元件表面 波纹度对其激光损伤阈值影响的评价方法的进一步说明，步骤四所述的获得每个特征频率 对应的相对激光损伤阈值的具体过程为
令ιρβη 表示理想情况下的晶体内部光强评价值，所述理想情况即无小尺度波纹时，
定义光强调制度为
其中，I 表示同等入射条件下、光波经小尺度波纹调制后的晶体内部的光强最大 max
值，从元件可能被破坏的角度看，调制度可用于判断系统运行的安全性，即调制度越大，光 学元件越容易发生诱导损伤破坏；
令/π表示理想情况下晶体的激光损伤阈值，J表示入射光波的光强，由于小尺度波
纹的调制作用，使得晶体内部的光强最大值为
当光强最大值与晶体的激光损伤阈值相等的时候，KDP晶体发生激光诱导损伤，一般为
体损伤，即
各本征模式场的振幅系数组成，可利用反射透射系数阵递推算法(RTCM)对上式求解，进而 获得整个空间的电磁场分布；
根据整个空间的电磁场分布，进而可以得到光学元件内部的光强分布/如下式令/々表示晶体实际的激光损伤阈值，定义相对激光损伤阈值RT (relative
由上述分析可知， RT 二 VM1。
具体实施方式
六由实施方式一所述的光学元件表面波纹度对其激光损伤阈值影 响的评价方法获得元件最佳加工参数的方法，它的过程如下
步骤Al、令C = |c(i),i = 1，2,.",Cq}表示机床的加工参数组，其中^为加工参数 的总数，由实际情况决定；获得所述加工参数组C = |ct/) i = 1,2,cq丨111每个加工参
数的实际取值范围，其中，参数的取值范围为[^ ，巧⑩],U1(I)和rj2(i)均
为实数，它们均由实际条件决定；
步骤A2、对每个加工参数C(Z),获得一个优选值，具体过程为
在每个加工参数的取值范围X^T2CO]内，选取g(z)个等间距的点，在
其余加工参数均固定的情况下，分别在该g(O个点的条件下，对光学元件进行加工，获
得个光学元件；利用光学元件表面波纹度对其激光损伤阈值影响的评价方法，获得
g(l)个光学元件中每个光学元件的评价结果，即获得每个光学元件的相对激光损伤阈 值，然后通过比较筛选，获得相对激光损伤阈值最小的光学元件，并将该光学元件对应的加
工参数LZ1i作为加工参数￡(/)的优选值，其中U1i eft/^XC/W)]；
步骤A4、根据步骤二获得的每个加工参数的优选值,获得优选参数组
U'，所述优选参数组CZ1 =IlriJ=I,.,Cq^ 即为待加工
元件的最佳加工参数组。本发明的获得元件最佳加工参数的方法，利用实施方式一的评价方法，能够获得 高加工质量的光学元件。
具体实施方式
七本实施方式是对实施方式六的获得元件最佳加工参数的方法的 进一步限定，其特征在于所述机床采用KDP晶体超精密加工机床，获得该机床在前角r = -45°时，KDP晶体的最佳加工参数组为
其中， 表示背吃刀量，/表示进给量。应用本实施方式，在KDP晶体超精密加工机床上对KDP样件进行SPDT法加工试 验，得到某一原始加工表面轮廓功率谱密度曲线，如图2所示。图3是利用白光干涉仪获 得的该加工表面的三维形貌结果图，图4至图8是图3中主要特征频率的三维形貌图。图 2中有几个较明显的峰值，其所对应空间频率的小尺度波纹成分在原始表面中占有较大比 例，该原始表面主要由这些小尺度波纹叠加而成。图9为相对阈值随小尺度波纹空间周期 的变化曲线；图10至图15为几个敏感周期所对应的晶体内部的光强分布。将相对阈值的最小值作为该机械加工表面对光学元件激光损伤阈值影响的评价 参数(亦即该机械加工表面的相对阈值)。对比图2和图9可见，空间周期为92. 5μπι的小 尺度波纹成分虽然在该原始加工表面中所占的比例不是最大的，但其对该光学元件损伤阈 值的影响却是最大的，元件该表面的激光损伤阈值主要由该频率信息决定。此外，空间周期 为117μπι和176 μ m的成分也不容忽视。若采取适当的检测手段找到引入这些空间频率成 分小尺度波纹的加工因素(如进给量、主轴跳动及导轨直线度等)，就能有效地采取措施对 光学元件的加工工艺过程进行改进，从而有效地提高光学元件的损伤阈值。在KDP专用超精密加工机床上采用SPDT法对晶体样件待加工表面进行了分区变 参数加工试验，并在成都光学精密研究中心进行了实际激光损伤阈值测定实验，损伤点实 际形貌见图16和图17。利用上述评价方法对各加工区的相对激光损伤阈值，并与实验中实际测得的阈值 结果进行了对比，图18为理论计算与实验结果的对比曲线，其中，“▼”点为理论计算值， “·”点为实验数据，Sl为理论拟合结果，S2为实验拟合结果，由图18知，理论评价结果与 实验结果吻合得很好，从而验证了理论及评价方法的正确性和可行性。由此可知，评价光学 元件已加工表面质量对其抗激光损伤性的影响程度可通过计算其相对激光损伤阈值来间 接说明。实际加工过程中，引入小尺度波纹的因素很多，例如导轨直线度，主轴的轴向跳 动、摆动，工件装卡变形，环境振动等。然而，限于现阶段的技术水平，完全控制上述各种 因素是很难的。因此，找出哪些因素对KDP晶体激光损伤阈值的影响是主要的，哪些是次要 的，对于实际应用具有重要意义。通过分区变参数加工试验，我们共获得了 21种不同的加 工表面并利用上述评价方法计算了每一原始表面的功率谱密度曲线和相对激光损伤阈值 曲线。图19是每一分区原始表面特征周期统计图，图20为决定表面相对阈值的主导周期 统计图。出现的次数从一定程度上能反映出具有该空间周期的小尺度波纹成分在KDP晶 体已加工表面中出现的概率，概率越大，说明该空间周期成分的随机性越小，亦即与加工工 艺过程中一些不变因素的相关性越大(如主轴转速、机床导轨直线度、工件装卡变形等)。由 图19可见，KDP晶体试件的表面轮廓中存在很多出现概率较大的周期性小尺度波纹成分， 其中102. 3μπι、125μπι和194. 1 μ m出现概率最大，且特征周期多集中在约90 μ πΓ350 μ m 范围内。由图20可见，引起KDP晶体损伤的小尺度波纹成分的空间周期出现在92. 1 μ m、 102. 9 μ m、116. 7 μ m和145. 8 μ m的概率较大，且集中分布于约90 μ πΓ 50 μ m范围内。因 此，引入这些空间周期的加工因素是降低KDP晶体激光损伤阈值的关键因素。此外，从图19 和图20的对比中，我们看到，具有主导周期92. 1μπι、102.9μπι和145. 8 μ m的小尺度波纹 成分不仅是引起KDP晶体损伤的主要因素，而且出现的概率也较大，是我们应极力避免引入的敏感对象。另外，计算结果表明，在90 μ πΓ 50 μ m的敏感区间内，元件的相对阈值随小尺度 波纹幅值的增大而减小，且当波纹幅值小降到IOnm以下时，KDP元件的相对激光损伤阈值 将提高到98%以上。因此，我们还可以通过尽量减小波纹幅值来提高其激光损伤阈值。以 加工表面三维轮廓算术平均偏差代表波纹度，图21是实验中测得的波纹度与进给量以及 背吃刀量的关系，其中，图21中的“ ”点为背吃刀量ap=10 μ m时得到的实验数据点，“■” 点为背吃刀量ap=15 μ m时得到的实验数据点，“▼”点为背吃刀量ap=20 μ m时得到的实验 数据点；由图21可见，针对我们试验所使用的KDP晶体超精密加工专用机床，在前角Y = -45°时其最优加工参数组合为
5 μ m ≤ ap ≤ 15 μ m, 3 μ m/r ≤ f ≤8 μ m/r
在此加工组合下，可以稳定地保证KDP元件加工表面粗糙度Ra=3飞nm，表面波纹度 Sa ^ 10 nm。在光学元件运行之前，提前对其激光损伤阈值作出评价，若不满足要求可对光学 元件进行二次加工，以免造成其不可恢复性损伤。本发明在光学元件运行之前，提前对其激光损伤阈值作出评价，若不满足要求可 对光学元件进行二次加工，以免造成其不可恢复性损伤。
权利要求
光学元件表面波纹度对其激光损伤阈值影响的评价方法，其特征在于它的过程如下步骤一、利用检测仪器，获取光学元件原始加工表面的形貌数据矩阵；步骤二、根据步骤一获得的形貌数据矩阵，获得光学元件原始加工表面的功率谱密度曲线，进而获得光学元件原始加工表面的各个特征频率以及每个特征频率的幅值；步骤三、对步骤二获得的每个特征频率，采用二维连续小波变换法提取并再现各特征频率的三维形貌，并利用傅立叶模方法计算每个特征频率对应的光学元件内部的光强分布；步骤四、根据步骤三获得的每个特征频率对应的光学元件内部的光强分布，获得每个特征频率对应的的光学元件内部的光强最大值，进而获得每个特征频率对应的相对激光损伤阈值；步骤五、对步骤四获得的每个特征频率对应的相对激光损伤阈值进行比较筛选，获得所有相对激光损伤阈值中的最小值，并将所述最小值作为此次对光学元件评价的结果。
2.根据权利要求1所述的光学元件表面波纹度对其激光损伤阈值影响的评价方法， 其特征在于步骤二所述内容的具体过程为令ζ Cr)表示步骤一获得的光学元件原始加工表面的形貌数据矩阵，其中ζ Cr)中包含 #个数据点，且每相邻两个数据点具有相同的采样间隔Δχ，整体采样长度为Δχ；采用如下公式获得光学元件原始加工表面的功率谱密度曲线 上式中，k为波数，k=2π fm’fm:m/、Ν Δ χ)为空间频率， 为采样点的序数，且-N/2≤ m ≤Ν/2 ；由功率谱密度曲线，即获得光学元件原始加工表面的各个特征频率/卩，再根据下式 计算获得每个特征频率A的幅值 其中，△f为取样频率。
3.根据权利要求1所述的光学元件表面波纹度对其激光损伤阈值影响的评价方法， 其特征在于步骤三所述内容的具体过程为二维连续小波变换的一般形式为 其中，ι = [J^x2 f是平面直角坐标，f(x)表示二维信号，wf(a^~b)表示二维连续小波变换 广是Jc1，jc2方向上的位移，Λ、&表达式中的上标r表示转置， 为尺度因子， 为坐标旋转因子，为坐标系逆时针旋转角 度， 表示二维基本小波函数的尺度伸缩、坐标旋转及二维位移， 的共轭；特征频率/；是与尺度a —一对应的，尺度与频率之间关系式为 f f 其中，/；为所采用小波基函数的原始中心频率；△为测量仪器的采样周期； 将所述原始中心频率/;、采样周期」以及待考察的特征频率代入上述关系式，即得 到特征频率/；对应的尺度α ；再利用yaw小波工具箱，完成对所述形貌数据矩阵的二维连 续小波变换；用正弦波来近似光学元件原始加工表面上的各频率的小尺度波纹，并建立小尺度波纹 的物理模型，所述小尺度波纹的物理模型位于x-y-z空间坐标系中，小尺度波纹的截面在 χ - ζ平面，方向为小尺度波纹的波纹线方向，小尺度波纹的基底平面垂直ζ轴、并沿χ轴 方向周期变化，该周期为r;入射光波以θ角入射到小尺度波纹表面、并通过小尺度波纹； 用水平多分层形状对小尺度波纹的形貌进行拟合；采用阶梯进行剖分处理，即沿Z轴将所求空间分成层，第1层为入射空气层，第层 为出射空气层，第Λτ 层为基底层，第2层至第Ρ0-2层为小尺度波纹层；小尺度波纹层的相对介电常数ε Cr)及相对磁导率#⑴均具有周期性7，即 ε ij) = e {x+T)，μ Cr) =U {x+T)，对于第/7 层有 其中，/7=2,3，…，表示第/7层一个周期内介质与空气分界面的坐标，第/7层的 实际介电常数为e⑴ε0, f。为真空介电常数，第/7层的实际磁导率为// (Χ) μ ,以0为真 空磁导率，ε b为基体材料相对介电常数，μ b为基体材料相对磁导率； 将第/7层的相对介电常数和相对磁导率一起表示为傅立叶模形式为 其中，Z7为傅里叶级数编号，εη为相对介电常数傅里叶展开后的第η项，μη为相对磁导率傅里叶展开后的第η项，[ε^μηΥ为第P层的4禾； 由几何关系得 其中ζ代表第β层上界面的ζ坐标，ζ—1代表第β-Ι层上界面的ζ坐标d代表小尺度 波纹幅值； 其中，万即表示万(X)，"即表示"(χ), σ=^μ0 ε ；A0 = 2 r/Λ , a为单位振幅的入射光波长=α 0+Λ m/T, a =sin θ , θ 为入射光与 ζ 轴的夹角， =0，士 1，士 2，…，士 I···， 为傅立叶模编号，I为计算时的截断常数# 、 、^;分别为电场的x、y和ζ向分量;力 、 hym、hzm分别为磁场的χ、y和ζ向分量；K表示第/7层波数的ζ分量，它是待求值； 每一分层中的电磁场满足麦克斯韦方程组其中，i 是磁感应强度，々是电位移矢量；利用傅立叶因式分解“逆规则”原理，得TE波的本征方程 其中，乂 =^r1 一- ζ Λζ) H、ρ和#分别是根据傅立叶因式分解原理得到的系 数矩阵，上角标-1代表求逆运算；上述TE波的本征方程为广义本征方程，对其求解得到由 1+1个本征矢组成的矩阵评；和由正负各2#+1个本征值组成的对角阵‘「+和J^，正值r+代表上行波，负值r_代表下行波；对于第层，电场强度 分量和磁场强度χ分量iif的最终表达式为 其中if ,cf为两个列矢量，由上行波的各本征模式场的振幅系数组成，Ζ由下行波的 各本征模式场的振幅系数组成，利用反射透射系数阵递推算法对上式求解，进而获得整个 空间的电磁场分布；根据整个空间的电磁场分布，进而得到光学元件内部的光强分布如下式
4.根据权利要求3所述的光学元件表面波纹度对其激光损伤阈值影响的评价方法， 其特征在于二维基本小波函数采用MeXican2D小波基函数，其表达式为 上式中， 和厶表示频域面坐标。
5.根据权利要求1所述的光学元件表面波纹度对其激光损伤阈值影响的评价方法， 其特征在于步骤四所述的获得每个特征频率对应的相对激光损伤阈值的具体过程为令 表示理想情况下的晶体内部光强评价值，定义光强调制度为 其中，I 表示同等入射条件下、光波经小尺度波纹调制后的晶体内部的光强最大 max值； 令I0表示理想情况下晶体的激光损伤阈值，j表示入射光波的光强，晶体内部的光 强最大值为 当光强最大值与晶体的激光损伤阈值相等的时候，KDP晶体发生激光诱导损伤，一般为 体损伤，即 令表示晶体实际的激光损伤阈值，定义相对激光损伤阈值RT为RT 二 Ia 丨 1Q ；亦有 ^r = IZA/"
6.由权利要求1所述的光学元件表面波纹度对其激光损伤阈值影响的评价方法获得 元件最佳加工参数的方法，其特征在于它的过程如下步骤Al、令C == 1，%...,cq}表示机床的加工参数组，其中为加工参数的总数；获得所述加工参数组C = ^c(I)J = 1，X...，cq |·中每个加工参数的实际取值范围，其中，参数的取值范围为[^! ，巧⑩],CZ1(I)和1/2( )均为实数； 步骤Α2、对每个加工参数C(Z),获得一个优选值，具体过程为 在每个加工参数C7(Z)的取值范围[巧(O，CZ2(0]内，选取g{i)个等间距的点，在 其余加工参数均固定的情况下，分别在该个点的条件下，对光学元件进行加工，获 得个光学元件；利用光学元件表面波纹度对其激光损伤阈值影响的评价方法，获得g(i)个光学元件中每个光学元件的评价结果，即获得每个光学元件的相对激光损伤阈 值，然后通过比较筛选，获得相对激光损伤阈值最小的光学元件，并将该光学元件对应的加工参数CZ1i作为加工参数的优选值，其中CZ1i ej^TjCO，·^ ]；步骤A4、根据步骤二获得的每个加工参数i^f)的优选值 ,获得优选参数组U' -{C/^i^lj，···，^^，所述优选参数组I/ ={t/3.，f =1，2，.··， |·即为待加工元件的最佳加工参数组。
7.根据权利要求6所述的获得元件最佳加工参数的方法，其特征在于所述机床采用 KDP晶体超精密加工机床，获得该机床在前角r= -45°时，KDP晶体的最佳加工参数组为5 μ m ^ ap ^ 15 μ m ； 3 μ m/r ^ / ^ 8 μ m/r ；其中， 表示背吃刀量，/表示进给量。
全文摘要
光学元件表面波纹度对其激光损伤阈值影响的评价方法及由此获得元件最佳加工参数的方法，涉及一种光学元件的表面质量评价方法以及一种获得元件最佳加工参数的方法。所述评价方法为首先获得原始加工表面的形貌数据矩阵，然后利用功率谱密度法、二维连续小波变换法及傅立叶模方法，获得各个特征频率对应的相对激光损伤阈值，然后选择其中的最小值作为评价结果；所述获得元件最佳加工参数的方法，即利用该评价方法，通过比较各种加工参数条件下得到的光学元件的相对激光损伤阈值，进而获得最优加工参数。本发明可用于评价光学元件的质量，并可用于指导光学元件的加工过程。
文档编号G01B21/30GK101887171SQ20101022249
公开日2010年11月17日 申请日期2010年7月9日 优先权日2010年7月9日
发明者姜伟, 李明全, 梁迎春, 王健, 许乔, 陈明君 申请人:哈尔滨工业大学