VARIABLE-STEP-DISTANCE MICRO-MILLING REPAIR CUTTER PATH GENERATING METHOD FOR DAMAGE POINTS ON SURFACE OF OPTICAL CRYSTAL

Abstract

A variable-step-distance micro-milling repair cutter path generating method for damage points on a surface of an optical crystal related to a field of optical material and optical element surface repair and includes steps of establishing a mathematical model of a repair profile; determining discrete contact points between a cutter and the repair profile to obtain a cutter contact control point set by a GPR path generating method to control a movement trend of a pseudo-random path; interpolating the cutter position control point set into a spatial curve by a NURBS modeling method; creating a UG curve in a UG software according to the mathematical model, and using the UG curve as the repair path to perform a machining process simulation. The method has good elimination effects on cutter marks with constant period and improves the ability of the KDP crystal to resist strong laser damage.

Claims

1. A variable-step-distance micro-milling repair cutter path generating method for damage points on a surface of an optical crystal, comprising: step 1: designing a damage repair profile with suitable depth and width according to a laser damage degree of a surface of a potassium dihydrogen phosphate (KDP) crystal and morphological characteristics of damage points; establishing a mathematical model of the damage repair profile; determining a conical repair profile as a laser-friendly repair profile; establishing a mathematical model of the conical repair profile, and a control equation is as follows: $\begin{array}{cc}z=f\left(x,y\right)=-\frac{D}{2h}\sqrt{x^2+y^2}+h& \left(1\right)\end{array}$ wherein x, y, and z respectively represent a Cartesian coordinate value of any point on the conical repair profile; D is a width of the conical repair profile; h is a depth of the conical repair profile; since the mathematical model of the conical repair profile is a revolving model, if ρ.sup.2=x.sup.2+y.sup.2, a coordinates of any point P on the conical repair profile is P=(x,y,z); the point P is expressed in polar coordinates as P=(ρ, θ) and satisfies: $\hspace{1em}\begin{array}{cc}\{\begin{array}{c}x=\rho \cos \theta \\ y=\rho \sin \theta \\ z=f\left(x,y\right)\end{array}& \left(2\right)\end{array}$ step 2: determining Gaussian pseudo-random (GPR) path generation parameters according to a processing requirement; determining discrete contact points between a cutter and the conical repair profile when the cutter is milling the conical repair profile to obtain a cutter contact control point set by a GPR path generating method; wherein the cutter contact control point set is configured to control a movement trend of a pseudo-random path; a specific process is as follow: adopting a spiral repair path with a small cutting depth and a small step distance; generating pseudo-random cutter contact control points by adjusting polar angles and polar diameters under a GPR strategy; when establishing a GPR path, a first cutter contact control point is initialized, and the first cutter contact control point is initialized to a point on an edge of the conical repair profile where a polar angle is 0, wherein
P.sub.CC.sup.0=(D/2,0)  (3) other cutter contact control points are generated by a polar diameter and polar angle control equation as follow: $\begin{array}{cc}\{\begin{array}{c}{\theta }_i={\theta }_{i-1}+\frac{K}{\lg \left({\rho }_{i-1}\right)+b}\\ {\rho }_i=\frac{D}{2}-\frac{{\theta }_i}{2\pi }\mu -{\delta }_i\end{array},i=1,2,3,\mathrm{.Math.}& \left(4\right)\end{array}$ above formula gives a distribution law of an i.sup.th cutter contact control point Pi CC=(ρ.sub.i,θ.sub.i); in a circumferential direction of the conical repair profile, a polar angle θi gives a stepped size $\frac{K}{\lg \left({\rho }_{i-1}\right)+b}$ related to ρ.sub.i−1 on a basis of θ.sub.i−1; K and b are polar angle adjustment parameters; by adjusting values of K and b, a sparseness of the cutter contact control points in the conical repair profile under different polar diameters is adjusted; a pseudo-random processing is performed on the polar diameters in a radial direction of the conical repair profile; where $\frac{D}{2}-\frac{{\theta }_i}{2\pi }\mu$ is a pole diameter of the cutter contact control points on the spiral repair path, and δ.sub.i is an offset of the cutter contact control points to the spiral repair path; when θ.sub.i∈[0, 2π), δ.sub.i=0; when θ.sub.i∈[2π, +∞), δ.sub.i is a pseudo-random number that satisfies a Gaussian distribution δ.sub.i˜N(μ,σ.sup.2); σ represents a standard deviation and is a pseudo-random adjustment parameter configured to adjust a change degree of the step distance; bandwidth constraint and adjacent constraint are performed on a generation range of the pseudo-random number; in (5), ε is a half-bandwidth of the pseudo-random number, and η is an adjacent constraint parameter of the cutter contact control points, $\hspace{1em}\begin{array}{cc}\{\begin{array}{c}\left|{\delta }_i\right|<E\\ |{\delta }_i-{\delta }_{i-1}|<\eta \end{array}& \left(5\right)\end{array}$ a spiral-like distribution of cutter contact control points is continuously introduced by performing formulas (4)-(5) until the cutter contact control points tangent to a bottom of the revolving model to form the cutter contact control point set; step 3: calculating central positions of a ball-end cutter one-to-one corresponding to the cutter contact control points to form a cutter position control point set by the established mathematical model of the repair profile and a size of a selected KDP crystal micro-milling cutter; wherein the central positions are defined as the cutter position control points; wherein the KDP crystal micro-milling cutter is a ball-end micro-milling cutter made of cubic boron nitride (CBN); the cutter position control points P.sub.CL are obtained to form the cutter position control point set according to a direction of a curved surface at each of the cutter position control points and a radius Rc of the ball-end micro-milling cutter, where a calculation equation is as follows: $\begin{array}{cc}\{\begin{array}{c}{P}_{CL}^{\left(j\right)}=P_{CC}^{\left(j\right)}+R_c\frac{\stackrel{.fwdarw.}{v^{\left(j\right)}}}{\left|\stackrel{.fwdarw.}{v^{\left(j\right)}}\right|}\\ \stackrel{.fwdarw.}{v^{\left(j\right)}}=\left(\frac{\partial f}{\partial x},\frac{\partial f}{\partial y},-1\right){|}_{x=x_i,y=y_i}\end{array}& \left(6\right)\end{array}$ step 4: interpolating the cutter position control point set into a spatial curve by a non-uniform rational B-splines (NURBS) modeling method; a mathematical model of the spatial curve is multiple cubic equation sets controlled by a unique parameter; wherein in the GPR strategy, the cutter position control points are the discrete contact points on the spiral repair path and are treated as value points; there are n+1 points in the cutter position control point set, the n+1 points are discrete value points of a k.sup.th NURBS curve; the k.sup.th NURBS curve is interpolated as the spiral repair path; a rational polynomial vector function of the NURBS curve is as follows: $\begin{array}{cc}p\left(u\right)=\frac{{.Math.}_{i=0}^n{\omega }_i{d}_i{N}_{i,k}\left(u\right)}{{.Math.}_{i=0}^n{\omega }_i{N}_{i,k}\left(u\right)}& \left(7\right)\end{array}$ wherein d.sub.i i is a curve control vertex and is obtained by an inverse calculation of a cutter position control points P.sub.CL; ω.sub.i is defined as a weight factor and is configured to assign weights to the curve control vertex; in a KDP crystal repair model, each curve control vertex di has a same weight when generating the NURBS curve, so each weight factor ω.sub.i=1 (i=0, 1, . . . , n); Ni,k(u) is a k.sup.th standard B-spline basis function defined by a node vector U=[u.sub.0, u.sub.1, . . . , u.sub.n+k+1] according to a recurrence formula Cox-De Boor; wherein $\hspace{1em}\begin{array}{cc}\{\begin{array}{c}{N}_{i,0}\left(u\right)=\{\begin{array}{c}1,u_i\le u\le u_{i+1}\\ 0,\mathrm{other}\end{array}\\ N_{i,k}\left(u\right)=\frac{u-u_i}{u_{i+k}-u_i}{N}_{i,k-1}\left(u\right)+\frac{u_{i+k+1}-u}{u_{i+k+1}-u_{i+1}}{N}_{i+1,k-1}\left(u\right)\\ \mathrm{define}\frac{0}{0}=0\end{array}& \left(8\right)\end{array}$ in the GPR strategy, the n+1 cutter position control points P.sub.CL have a node vector U=[u0, u1, . . . , u.sub.n+6], in order to determine a u.sub.i+3 corresponding to Pi CL (i=0, 1, . . . , n), the value points are parameterized; considering a large change in distances between the cutter position control points, a cumulative chord length parameterization method of a formula (9) is applied to reflect a distribution of data points according to a chord length: $\begin{array}{cc}\{\begin{array}{c}{u}_0=u_1=u_2=u_3=0\\ \begin{array}{c}{u}_{i+3}=u_{i+2}+\sqrt{\left|P_{CL}^i-P_{CL}^{i-1}\right|}/{.Math.}_{i=1}^n\sqrt{\left|P_{CL}^i-P_{CL}^{i-1}\right|}\\ i=1,2,\mathrm{.Math.},n-1\end{array},\\ u_{n+3}=u_{n+4}=u_{n+5}=u_{n+6}=1\end{array}& \left(9\right)\end{array}$ as long as there are control vertices di, a NURBS repair path of the ball-end micro-milling cutters is calculated; then weighted control vertices of the NURBS curve are inversely calculated; for a cubic NURBS curve generated by n+1 P.sub.CL, there are n+3 control vertices, and the cubic NURBS curve is defined as a formula (10), where u∈[u.sub.i, u.sub.i+1]⊂[u.sub.3, u.sub.n+3]; $\begin{array}{cc}{p}_i\left(u\right)=\underset{j=i-3}{\overset{i}{.Math.}}{d}_j{N}_{j,k}\left(u\right)& \left(10\right)\\ p\left(u_{i+3}\right)=\underset{j=i}{\overset{i+3}{.Math.}}{d}_j{N}_{j,3}\left(u_{i+3}\right)=P_{CL}^i,i=0,1,\mathrm{.Math.},n& \left(11\right)\end{array}$ substituting nodal values in a definition domain of the cubic NURBS curve into the equations in turn, and the cubic NURBS curve of different segments are continuous at the value points, which is defined as a formula (11); n+1 equations in the above formulas are not enough to determine n+3 unknown weighted control vertices, so two additional equations (12) are given by boundary conditions: $\begin{array}{cc}{p}_0^{\prime }\left(0\right)=\frac{3\frac{{\omega }_1}{{\omega }_0}\left(d_1-d_0\right)}{u_4-u_3}.p_n^{\prime }\left(1\right)=\frac{3\frac{{\omega }_{n+1}}{{\omega }_{n+2}}\left(d_{n+2}-d_{n+1}\right)}{u_{n+3}-u_{n+2}}& \left(12\right)\\ \left[\begin{array}{ccccc}1& & & & \\ a_2& b_2& c_2& & \\ & \ddots & \ddots & \ddots & \\ & & a_n& b_n& c_n\\ & & & & 1\end{array}\right].Math.\left[\begin{array}{c}{d}_1\\ d_2\\ \mathrm{.Math.}\\ d_n\\ d_{n+1}\end{array}\right]=\left[\begin{array}{c}{e}_1\\ e_2\\ \mathrm{.Math.}\\ e_n\\ e_{n+1}\end{array}\right]& \left(13\right)\end{array}$ linear equation sets in the equation (13) are obtained by the two additional equations (12); and all weighted control vertices are obtained; then the weighted control vertices are substituted into the equation (10) to obtain a NURBS curve equations set after interpolation of P.sub.CL; there are n sets of the linear equation sets, and each of the linear equation sets has three equations; the three equations are the equations of u on x, y, and z components of a curve between adjacent P.sub.CL; the three equations are second-order continuous at the P.sub.C1; parameter equation is used as the repair path of the KDP crystal to realize digitization of the repair path, and a only related parameter is u; step 5: creating a UG curve in a UG software according to the mathematical model of the spatial curve, and using the UG curve as the repair path to performed a machining process simulation; if a simulation result meets safety and manufacturability requirements of a machining cutter path, a step 6 is performed; if not, go back to the step 2 to check the processing requirement and adjust the GPR path generation parameters; wherein in the step 5, a calculated NURBS curve equation is imported into a CAM module of the UG software; the machining process simulation is performed in a curve-driven manner, then whether the simulation result meets the safety and manufacturability requirements of the machining cutter path is checked; if the simulation result meets the safety and manufacturability requirements of the machining cutter path, the step 6 is performed; if not, return to step 2 to check the process and modify the GPR path generation parameters σ; step 6, converting the machining process simulation of the step 5 into general NC codes, and performing precision micro-milling repair on the KDP crystal on a KDP crystal repair machine; wherein in the step 6, a conventional repair process simulation and a GPR variable-step-distance repair process simulation are converted into the NC codes through a post-processor of the UG software, and the precision micro-milling repair of the KDP crystal is performed on the KDP crystal repair machine.

2. The variable-step-distance micro-milling repair cutter path generating method for the damage points on the surface of the optical crystal according to claim 1, wherein the damage repair profile in the step 1 is the conical repair profile; a depth and a width of the conical repair profile are selected according to morphology, size and distribution of the damage points on the surface of the KDP crystal; and the mathematical model of the conical repair profile is a quadric surface equation; wherein the processing requirement in step 2 is a cutting step distance corresponding to a conventional spiral milling path; the GPR path generation parameters comprise polar angle adjustment parameters and polar diameter offset; according to generated parameters, a series of coordinates of the discrete contact points of a surface of the conical repair profile is calculated by a GPR algorithm; the discrete contact points are the cutter contact control points; when the cutter is milling the conical repair profile, a contact point set to a theoretical surface is the spatial curve, and the cutter contact control points are configured to control a trend of the spatial curve; the micro-milling cutter in the step 3 is the ball-end micro-milling cutter made of CBN, the ball-end micro-milling cutter is configured to realize a plastic domain processing of a repair profile on the surface of the KDP crystal; when the ball-end micro-milling cutter performs milling motion on a surface of the conical repair profile of the KDP crystal, a sphere center motion path of the ball-end micro-milling cutter is the spatial curve that defined as a cutter motion path or the cutter path; a point at a distance of the cutter radius outward from a normal direction of a cutter contact point is defined as a cutter position control point; the cutter position control points are configured to control the trend of the cutter path.

3. The variable-step-distance micro-milling repair cutter path generating method for the damage points on the surface of the optical crystal according to claim 2, wherein in the step 4, the cutter position control points are acted as curve value points, by the NURBS modeling method, to interpolate the k.sup.th NURBS spline curve.

4. The variable-step-distance micro-milling repair cutter path generating method for the damage points on the surface of the optical crystal according to claim 3, wherein in step 5 further comprises creating the UG curve in the UG software according to the mathematical model of the spatial curve, saving obtained curve equations as an .exp format file, and importing the .exp format file into the UG software in an expression form to make a regular curve; and performing an intuitive machining process simulation to confirm safety and reliability of the machining process by driving the regular curve; wherein in step 6, a process of converting the machining process simulation of step 5 into the NC codes depends on the post-processor of the UG software.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0045] FIG. 1 is a flow chart of a GPR micro-milling repair cutter path generating method for damage points on a surface of a large-diameter KDP crystal.

[0046] FIG. 2 is a schematic diagram of a variable-step-distance cutter repair path.

[0047] FIG. 3A is a schematic diagram of an obtained repair morphology of the KDP crystal after the conventional repair,

[0048] FIG. 3B is a schematic diagram of an obtained a repair morphology of the KDP crystal after the variable-step-distance repair of the present disclosure.

[0049] FIG. 4 a schematic diagram showing PSD curves of different repair methods.

DETAILED DESCRIPTION

[0050] In an example analysis of GPR micro-milling repair cutter path generating method for damage points on a surface of a KDP crystal, an experimental verification of a variable-step-distance repair path generating method is verified according to a process shown in FIG. 1.

[0051] A process of the variable-step-distance micro-milling repair cutter path generating method for the damage points on the surface of the optical crystal is as follow:

[0052] 1) designing a damage repair profile with a suitable depth and width according to a laser damage degree of the surface of the KDP crystal and the morphological characteristics of damage points; establishing a mathematical model of the damage repair profile;

[0053] A laser damage experiment shows that laser damage due to ablation often occurs on the surface if the KDP crystal, which is characterized by large damage area and shallow depth. The conical repair profile is able to remove the damaged area well. Verification by optical simulation shows that the conical repair profile is a laser-friendly repair profile, which effectively improves ability of the KDP crystal to resist laser damage after repair.

[0054] Establishing a mathematical model of the conical repair profile by computer technologies, and its control equation is as follows:

[00013]$\begin{array}{cc}z=f\left(x,y\right)=-\frac{D}{2h}\sqrt{x^2+y^2}+h& \left(1\right)\end{array}$

[0055] x, y, and z respectively represent a Cartesian coordinate value of any point on the conical repair profile. D is a width of the conical repair profile. h is a depth of the conical repair profile. Since the mathematical model of the conical repair profile is a revolving model, if ρ2=x2+y2, a coordinates of any point P on the conical repair profile is P=(x,y,z). The point P is expressed in polar coordinates as P=(ρ, θ) and satisfies:

[00014]$\begin{array}{cc}\{\begin{array}{c}x=\mathrm{\rho cos\theta }\\ y=\mathrm{\rho sin\theta }\\ z=f\left(x,y\right)\end{array}& \left(2\right)\end{array}$

[0056] 2) In a conventional repair method, UG software is adopted to generate two processes of cavity milling (layer milling) and area profile milling (spiral milling). A purpose of the cavity milling is to quickly remove ablated portions of the damaged points and initially mill out the conical repair profile. Therefore, a large step distance and large depth of cut repair process is used in a molding area. A purpose of area profile milling is to further improve a surface quality of the repair profile and process a “laser-friendly” smooth repair profile, so a spiral repair path with a small cutting depth and a small step distance is adopted.

[0057] Pseudo-random cutter contact control points are generated by adjusting polar angles and polar diameters under a GPR strategy. When establishing a GPR path, a first cutter contact control point is initialized, and the first cutter contact control point is initialized to a point on an edge of the conical repair profile; where a polar angle is 0, where

P.sub.CC.sup.0=(D/2,0)  (3)

[0058] Other cutter contact control points are generated by a polar diameter and polar angle control equation as follow:

[00015]$\begin{array}{cc}\{\begin{array}{c}{\theta }_i={\theta }_{i-1}+\frac{K}{\lg \left({\rho }_{i-1}\right)+b}\\ {\rho }_i=\frac{D}{2}-\frac{\theta }{2\pi }\mu -{\delta }_i\end{array},i=1,2,3,\mathrm{.Math.}& \left(4\right)\end{array}$

[0059] The above formula gives a distribution law of an i.sup.th cutter contact control point Pi CC=(ρ.sub.i,θ.sub.i). In a circumferential direction of the conical repair profile, a polar angle δi gives a stepped size

[00016]$\frac{K}{\lg \left({\rho }_{i-1}\right)+b}$

related to ρ.sub.i−1 on a basis of θ.sub.i−1. K and b are polar angle adjustment parameters. By adjusting values of K and b, a sparseness of the cutter contact control points in the conical repair profile under different polar diameters is adjusted. A pseudo-random processing is performed on the polar diameters in a radial direction of the conical repair profile.

[00017]$\frac{D}{2}-\frac{{\theta }_i}{2\pi }\mu$

is a pole diameter of the cutter contact control points on the spiral repair path, and δ.sub.i is an offset of the cutter contact control points to the spiral repair path. When θ.sub.i∈[0, 2π), δ.sub.i=0. When θ.sub.i∈[2π, +∞), δ.sub.i is a pseudo-random number that satisfies a Gaussian distribution δ.sub.i˜N(μ,σ.sup.2). σ represents a standard deviation and is a pseudo-random adjustment parameter configured to adjust a change degree of the step distance;

[0060] In order to prevent the pseudo-random cutter contact control points from appearing in uncontrollable positions with a small probability, bandwidth constraint and adjacent constraint as shown in formula (5) are performed on a generation range of the pseudo-random number. In (5), ε is a half-bandwidth of the pseudo-random number, and η is an adjacent constraint parameter of the cutter contact control points.

[00018]$\begin{array}{cc}\{\begin{array}{c}.Math.{\delta }_i.Math.<\mathrm{.Math.}\\ .Math.{\delta }_i-{\delta }_{i-1}.Math.<\eta \end{array}& \left(5\right)\end{array}$

[0061] A spiral-like distribution of cutter contact control points is continuously introduced by performing formulas (4)-(6) until the cutter contact control points tangent to a bottom of the revolving model to form the cutter contact control point set;

[0062] 3) The KDP crystal micro-milling cutter is a ball-end micro-milling cutter made of CBN. Theoretically, when doing milling motion, the ball-end micro-milling cutter and the conical repair profile are only tangent at one point, The cutter position control points P.sub.CL are obtained to form the cutter position control point set according to a direction of a curved surface at each of the cutter position control points and a radius Rc of the ball-end micro-milling cutter, where a calculation equation is as follows:

[00019]$\begin{array}{cc}\{\begin{array}{c}{P}_{\mathrm{CL}}^{\left(i\right)}=P_{\mathrm{CC}}^{\left(i\right)}+R_c\frac{\stackrel{.fwdarw.}{v^{\left(i\right)}}}{.Math.\stackrel{.fwdarw.}{v^{\left(i\right)}}.Math.}\\ \stackrel{.fwdarw.}{v^{\left(i\right)}}=\left(\frac{\partial f}{\partial x},\frac{\partial f}{\partial y},-1\right){|}_{x=x_i,y=y_i}\end{array}& \left(6\right)\end{array}$

[0063] 4) There are two forms of NURBS curve in the actual use. One form is that the control vertex is given to find a curve expression or any point on the NURBS curve, which is defined as a forward calculation problem. The other form is to find the curve expression or any point on the NURBS curve when several value points on the NUBRS curve are given. In the process, it is necessary to inversely calculate the curve control vertex, so it is defined as an inverse calculation problem. In the GPR strategy, the cutter position control points are the discrete contact points on the spiral repair path and are treated as value points. Therefore, the interpolation of the repair path in the GPR strategy is the inverse calculation problem.

[0064] There are n+1 points in the cutter position control point set. The n+1 points are discrete value points of a NURBS curve. The NURBS curve is a k.sup.th NURBS curve and is interpolated as the spiral repair path. A rational polynomial vector function of the NURBS curve is as follows:

[00020]$\begin{array}{cc}p\left(u\right)=\frac{{.Math.}_{i=0}^n{\omega }_i{d}_i{N}_{i,k}\left(u\right)}{{.Math.}_{i=0}^n{\omega }_i{N}_{i,k}\left(u\right)}& \left(7\right)\end{array}$

[0065] d.sub.i i is a curve control vertex and is obtained by an inverse calculation of a cutter position control points P.sub.CL. ω.sub.i is defined as a weight factor and is configured to assign weights to the curve control vertex; in a KDP crystal repair model. Each curve control vertex di has a same weight when generating the NURBS curve, so each weight factor ω.sub.i=1 (i=0, 1, . . . , n); Ni,k(u) is a k.sup.th standard B-spline basis function defined by a node vector U=[u.sub.0, u.sub.1, . . . , u.sub.n+k+1] according to a recurrence formula Cox-De Boor; wherein

[00021]$\begin{array}{cc}\{\begin{array}{c}{N}_{i,0}\left(u\right)=\{\begin{array}{c}1,u_i\le u\le u_{i+1}\\ 0,\mathrm{other}\end{array}\\ N_{i,k}\left(u\right)=\frac{u-u_i}{u_{i+k}-u_i}{N}_{i,k-1}\left(u\right)+\frac{u_{i+k+1}-u}{u_{i+k+1}-u_{i+1}}{N}_{i+1,k-1}\left(u\right)\\ \mathrm{define}\frac{0}{0}=0\end{array}& \left(8\right)\end{array}$

[0066] According to the above formulas, if k=3, the NURBS curve is a cubic NURBS curve. In the GPR strategy, the n+1 cutter position control points P.sub.CL have a node vector U=[u0, u1, . . . , un+6]. In order to determine a u.sub.i+3 corresponding to Pi CL (i=0, 1, . . . , n), the value points are parameterized. Considering a large change in distances between the cutter position control points, a cumulative chord length parameterization method of a formula (9) is applied to reflect a distribution of data points according to a chord length.

[00022]$\begin{array}{cc}\{\begin{array}{c}{u}_0=u_1=u_2=u_3=0\\ u_{i+3}=u_{i+2}+\sqrt{.Math.P_{\mathrm{CL}}^i-P_{\mathrm{CL}}^{i-1}.Math.}/\\ {.Math.}_{i=1}^n\sqrt{.Math.P_{\mathrm{CL}}^i-P_{\mathrm{CL}}^{i-1}.Math.},i=1,2,\mathrm{.Math.},n-1\\ u_{n+3}=u_{n+4}=u_{n+5}=u_{n+6}=1\end{array}& \left(9\right)\end{array}$

[0067] As long as there are control vertices di, a NURBS repair path of the ball-end micro-milling cutters is calculated. Then weighted control vertices of the NURBS curve are inversely calculated;

[0068] For a cubic NURBS curve generated by n+1 P.sub.CL, there are n+3 control vertices, and the cubic NURBS curve is defined as a formula (10), where u∈[u.sub.i, u.sub.i+1]⊂[u.sub.3, u.sub.n+3].

[00023]$\begin{array}{cc}{p}_i\left(u\right)=\underset{j=i-3}{\overset{i}{.Math.}}{d}_j{N}_{j,k}\left(u\right)& \left(10\right)\\ p\left(u_{i+3}\right)=\underset{j=i}{\overset{i+3}{.Math.}}{d}_j{N}_{j,3}\left(u_{i+3}\right)=P_{CL}^i,i=0,1,\mathrm{.Math.},n.& \left(11\right)\end{array}$

[0069] Substituting nodal values in a definition domain of the cubic NURBS curve into the equations in turn, and curves of different segments are continuous at the value points P.sub.CL, which is defined as a formula (11); n+1 equations in the above formulas are not enough to determine n+3 unknown weighted control vertices, so two additional equations (12) given by boundary conditions:

[00024]$\begin{array}{cc}{p}_0^{\prime }\left(0\right)=\frac{3\frac{{\omega }_1}{{\omega }_0}\left(d_1-d_0\right)}{u_4-u_3},p_n^{\prime }\left(1\right)=\frac{3\frac{{\omega }_{n+1}}{{\omega }_{n+2}}\left(d_{n+2}-d_{n+1}\right)}{u_{n+3}-u_{n+2}}& \left(12\right)\\ \left[\begin{array}{ccccc}1& & & & \\ a_2& b_2& c_2& & \\ & \ddots & \ddots & \ddots & \\ & & a_n& b_n& c_n\\ & & & & 1\end{array}\right].Math.\left[\begin{array}{c}{d}_1\\ d_2\\ \mathrm{.Math.}\\ d_n\\ d_{n+1}\end{array}\right]=\left[\begin{array}{c}{e}_1\\ e_2\\ \mathrm{.Math.}\\ e_n\\ e_{n+1}\end{array}\right]& \left(13\right)\end{array}$

[0070] Linear equation sets in the equation (13) are obtained by additional conditions; and all weighted control vertices are obtained. Then the weighted control vertices are substituted into the equation (10) to obtain a NURBS curve equations set after interpolation of P.sub.CL. There are n sets of the linear equation sets, and each of the linear equation sets has three equations. The three equations are the equations of u on x, y, and z components of a curve between adjacent P.sub.CL. The three equations are second-order continuous at the P.sub.C1. The parameter equation is used as the repair path of the KDP crystal to realize digitization of the repair path, and an only related parameter is u.

[0071] 5) The conventional repair process is to establish the damage points and a cutter model in the UG software. After setting process parameters, the post-processor generates the NC codes that is able to be recognized by the KCP crystal repair machine. The GPR repair method is to calculate a repair path according to a geometric model, and any point on the repair path has a value u corresponding to it. Therefore, in theory, the coordinates of any point is able to be output to the NC codes.

[0072] The present disclosure combines numerical calculation and computer-aided manufacturing technology and in the present disclosure, a calculated NURBS curve equation is imported into a CAM module of the UG software; as shown in FIG. 2, the machining process simulation is performed in a curve-driven manner, then whether the simulation result meets the safety and manufacturability requirements of the machining cutter path is checked; if the simulation result meets the safety and manufacturability requirements of the machining cutter path, the step 6 is performed; if not, return to step 2 to check the process and modify the GPR path generation parameters a.

[0073] 6) A conventional repair process simulation and a GPR variable-step-distance repair process simulation are converted into the NC codes through the post-processor of the UG software, and the precision micro-milling repair of the KDP crystal is performed on the KDP crystal repair machine.

[0074] The width of the conical repair profile is 800 μm, the depth of the width of the conical repair profile is 40 μm, and an average step distance μ of the area profile milling and GPR repair methods are both 20 μm, taking σ=5 μm, ε=2σ, and η as a large number. After the processing was completed, a white light interferometer is configured to measure morphology information of the repaired surface of the KDP crystal. The morphology comparison is shown in FIG. 3A and FIG. 3B.

[0075] From the morphology comparison of FIG. 3A and FIG. 3B, it can be seen that a period of cutter marks is significantly reduced in the repair profile processed by the GPR repair method. In order to analyze periodicity of the cutter marks, the present disclosure adopts the power spectral density (14) to characterize it. The power spectral density is a way to extract frequency information in wave signals through Fourier transform, which reflects different frequency components in complex wave signals.

[00025]$\begin{array}{cc}PSD\left(f\right)=\frac{1}{2\pi N}|{.Math.}_{n=0}^N\left(S\left(x_i\right)-P\left(x_i\right)\right)e^{-j2\pi \mathrm{fn}}{|}^2& \left(14\right)\end{array}$

[0076] N is the number of radial sampling points, f is the spatial frequency (the reciprocal of the distance). S (xi) is a profile of a radial cross-section. P (xi) is the minimum square line of S (xi). The calculated PSD curves of the two repair profile are shown in FIG. 4.

[0077] When the step distance of the repair path is 20 μm, an inherent knife frequency of the repair path is 0.05 μm-1. An existence of the frequency information greatly enhances interference diffraction effect of the KDP crystal in the laser environment. By calculating the PSD curve of the entire repair profile, the PSD value of a sensitive frequency band is reduced from 1.746 μm*μm2 to 0.297 μm*μm2 after processed by the GPR repair method, that is, a drop is 83.0%. It should be noted that the PSD value of middle and low frequency bands has increased to a certain extent. This is because that a residual height of the GPR repair profile is high in some positions and is low in some other positions, which eventually causes straightness of the sampling points to decrease and the low frequency components to increase. However, these residual cutter marks of low frequency and nanometer amplitude do not negatively affect the laser damage resistance of the repaired KDP crystal.

[0078] The above steps adopts the process flow of the present disclosure to realize generation of variable-step-distance repair path for the KDP crystal, which proves that the present disclosure has a good elimination effect on cutter marks with constant period and improves the ability of the KDP crystal to resist strong laser damage.