*Article* **A New Computational Model of Step Gauge Calibration Based on the Synthesis Technology of Multi-Path Laser Interferometers**

## **Guoying Ren 1,2,\*, Xinghua Qu <sup>1</sup> and Xiangjun Chen <sup>1</sup>**


Received: 20 February 2020; Accepted: 17 March 2020; Published: 19 March 2020

**Featured Application: Inspired by the principle of closed-loop vector model, this paper proposes a new computational model to realize the high precision calibration of the step gauge based on the commercial CMM platform and laser interference technology. According to this method, only real-time measurements of the attitude and displacement changes in the movement of the step gauge are needed for the accurate length value of the step gauge to be obtained.**

**Abstract:** A step gauge is a commonly used length standard for international comparison, and its calibration accuracy is often used as a sign to measure a country's length Calibration and Measurement Capability (CMC). Based on this, some developed countries and developing countries all over the world have been carrying out the research of precision calibration technology for step gauge. On the basis of summarizing the current situation of step gauge calibration technology in other countries, this paper presents a new computational model of step gauge calibration based on the Synthesis Technology of Multi-Path Laser Interferometers (SMLI) and an auto-collimator, which can synthesize the three laser light paths into the measured centerline of step gauge. It is very important to obtain a good measurement accuracy for the step gauge, conformed to the Abbe principle, no matter where it is installed on the CMM measurement platform. In this paper, the development of the mathematical model, the data collection algorithms, data analysis techniques, and measurement uncertainty budgets are discussed. Finally, the experimental measurement is carried out and the measurement accuracy is verified to be effective. The results show that this method can effectively avoid the influence of Abbe error in length measurement, and significantly enhance the calibration accuracy of the step gauge.

**Keywords:** metrology; step gauge; length calibration; multi-path laser synthesis technology

## **1. Introduction**

As a universal length standard for calibrating instruments, step gauges are widely used in the precision calibration and error compensation of high-precision instruments, such as the Coordinate Measuring Machine (CMM), numerical control machine tool, optical instruments and so on.

Moreover, the International Bureau of Metrology (BIPM) has also listed the step gauge as a key length comparison, making it an important component in measuring the metrological capacity of a national geometric measurement calibration laboratory [1–4]. Many ambitious countries have established their own step gauge calibration devices [5–7].

On the basis of Leitz CMM, the National Metrology Institute of Japan (NMIJ) realize their step gauge calibration by use of the large displacement measurement of a laser interference measuring system, and the micro displacement of the probe of the CMM, thereby obtaining the dimension values of the measured step gauge. The measurement uncertainty of this device is U= (0.06 + 0.22L) μm [8].

Based on the Moore M48 CMM, the National Institute of Standards and Technology (NIST) has developed a step gauge calibration device using a specially designed induction type probe that provides a short-term repeatability of 7 nm [9].

The above two calibration systems are based on the measurement platform of CMM, but the Abbe error, pitch and yaw errors are considered only as the error sources of measurement uncertainty, rather than being directly measured and compensated.

The National Institute of Metrology of Finland (VTT MIKES) has constructed a step gauge calibration system on a vibration isolated stone table with a 1D movable guide rail, rather than CMM, to eliminate mechanical disturbances [10]. The measurement uncertainty of the system is U = Q[0.090; 0.14L] μm. This system also considers the pitch and yaw errors as the error sources of measurement uncertainty, and so does not directly measure them.

At Physikalisch-Technische Bundesanstalt (PTB), the Nanometer Comparator is upgraded with a tactile sensor system to add the capability to measure step gauges. The 1σ repeatability of the device is better than 7 nm [11]. The system adopts the precise displacement measurement method based on the principle of vacuum interference, which has high precision, but is relatively complex.

Based on the analysis of the advantages and disadvantages of the above-mentioned instruments, this paper proposes a new computational model of step gauge calibration based on the Synthesis Technology of Multi-Path Laser Interferometers (SMLI), in which three of four interferometers are used for distance measurement, and the fourth channel is used for the wavelength compensation of the laser [12–14], thereby establishing our own step gauge calibration device on the platform of Leitz CMM. By adding four laser interferometers in the measurement platform, the spatial position relationship model of three out of four lasers is constructed, and the synthetic optical path of these lasers is translated to the measurement line of the measured step gauge, which can effectively eliminate the influence of the Abbe error of the CMM on the measurement result, and compensate the errors of pitch and yaw. Based on the presented calculation model of the closed-loop vector principle, the error caused by pitch and yaw can be corrected to the measured length in real time by measuring attitude change with an auto-collimator to obtain a better calibration accuracy.

This paper will firstly introduce the composition and structure of the developed step gauge measurement system by in detail. Then, the principle and implementation of the measurement method are be addressed systematically. Finally, the experimental results are given and discussed to provide a conclusion.

## **2. The Mathematical Modeling and Simulation**

In order to eliminate the Abbe measurement error, we designed a mathematical model of SMLI which can synthesize the three laser light paths into the one measured center line of the step gauge, as shown in Figure 1. The P1, P2 and P3 are the installation points of three reflector mirrors, and the O is the measurement point of the step gauge. The three reflector mirrors and the step gauge are installed in the same moveable platform. The P1P <sup>1</sup> is the first laser measurement path during the measurement that is marked as l1*u*<sup>1</sup> (l1 is the measurement value of the first interferometer, and *u*<sup>1</sup> is its unit directional vector), the P2P <sup>2</sup> is the second laser measurement path that is marked as l2*u*<sup>2</sup> (l2 is the measurement value of the second interferometer, and *u*<sup>2</sup> is its unit directional vector), the P3P 3 is the third laser measurement path that is marked as l3*u*<sup>3</sup> (l3 is the measurement value of the third interferometer, and *u*<sup>3</sup> is its unit directional vector), and the O O (the O point is the first measurement position and the O' is the nth measurement position) is the synthesized measurement path of the measured step gauge that is marked as r (r is the only unknown quantity). The β<sup>i</sup> is the angle between the direction of OPi and X axis of coordinate system, which is known and definite. The θ<sup>i</sup> is the angle between the direction of PiPi' (i = 1,2,3)and the Z axis of the coordinate system, which represents the change of pitch angle and yaw angle during the measurement and is measured by the auto-collimator.

**Figure 1.** Mathematical model of Synthesis Technology of Multi-Path Laser Interferometers (SMLI).

According to Figure 1, the mathematical model based on the closed-loop vector principle [15–17] is as follows.

$$
\overrightarrow{\dot{\mathbf{r}}} = \overrightarrow{\mathbf{e}}\_i + l\_i \overrightarrow{\mathbf{u}}\_i - \overrightarrow{\mathbf{e}}\_i,\\
\text{i } = 1, 2, 3,
\end{pmatrix} \tag{1}
$$

where <sup>→</sup> *e<sup>i</sup>* = *ei* cos <sup>β</sup>*<sup>i</sup>* sin <sup>β</sup>*<sup>i</sup>* <sup>0</sup> *<sup>T</sup>* (*<sup>i</sup>* <sup>=</sup> 1, 2, 3), → *OPi* = → *ei*,it is known.

where <sup>→</sup> *u<sup>i</sup>* = (sin θ*<sup>i</sup>* cos β*<sup>i</sup>* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* cos θ*i*) *<sup>T</sup>*,it is known. Then, *<sup>l</sup>*<sup>i</sup> is also known. Changing Equation (1) to Equation (2):

$$\left(\left(r - \mathbf{e}\_i - l\_i \mathbf{u}\_i\right)^T \left(r - \mathbf{e}\_i - l\_i \mathbf{u}\_i\right) = \mathbf{e}\_i^2\right. \tag{2}$$

Expanding Equation (2) to Equation (3) is done as follows.

$$\left[\left(\mathbf{r} - \mathbf{e}\_i\right)^T - \left(l\_i \mathbf{u}\_i\right)^T\right] \left[\left(\mathbf{r} - \mathbf{e}\_i\right) - \left(l\_i \mathbf{u}\_i\right)\right] - \mathbf{e}\_i^2 = \mathbf{0}.\tag{3}$$

Then, Equation (3) is expanded as follows.

$$\begin{split} & \implies \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right)^{T} \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right) - \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right)^{T} \left(l\_{i}\boldsymbol{u}\_{i}\right) - \left(l\_{i}\boldsymbol{u}\_{i}\right)^{T} \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right) + l\_{i}^{2} - \boldsymbol{\mathsf{e}}\_{i}^{2} = \boldsymbol{0}, \\ & \implies \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right)^{T} \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right) + l\_{i}^{2} - \boldsymbol{\mathsf{e}}\_{i}^{2} - \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right)^{T} \left(l\_{i}\boldsymbol{u}\_{i}\right) - \left(l\_{i}\boldsymbol{u}\_{i}\right)^{T} \left(\boldsymbol{r} - \boldsymbol{\mathsf{e}}\_{i}\right) = \boldsymbol{0}. \end{split}$$

Suppose the coordinate of point O is (0,0,0), and the coordinate of point O' is (x,y,z), then:

$$\mathbf{e}\_i \mathbf{G}\_i - l\_i \mathbf{r}^T \mathbf{u}\_i + l\_i \mathbf{e}\_i^T \mathbf{u}\_i - l\_i \mathbf{u}\_i^T \mathbf{r} + l\_i \mathbf{u}\_i^T \mathbf{e}\_i = \mathbf{0},\tag{4}$$

where:

*Gi* = (*r* − *ei*) *<sup>T</sup>*(*<sup>r</sup>* <sup>−</sup> *<sup>e</sup>i*) <sup>+</sup> *<sup>l</sup>* 2 *<sup>i</sup>* <sup>−</sup> *<sup>e</sup>*<sup>2</sup> *<sup>i</sup>* <sup>=</sup> *<sup>r</sup>T<sup>r</sup>* <sup>−</sup> *<sup>r</sup>Te<sup>i</sup>* <sup>−</sup> *<sup>e</sup><sup>T</sup> <sup>i</sup> <sup>r</sup>* <sup>+</sup> *<sup>e</sup><sup>T</sup> <sup>i</sup> e<sup>i</sup>* + *l* 2 *<sup>i</sup>* <sup>−</sup> *<sup>e</sup>*<sup>2</sup> *i* = *<sup>x</sup>*<sup>2</sup> + *<sup>y</sup>*<sup>2</sup> + *<sup>z</sup>*<sup>2</sup> <sup>−</sup> *xyz* ⎛ ⎜⎜⎜⎜⎜⎜⎜⎜⎝ *ei* cos β*<sup>i</sup> ei* sin β*<sup>i</sup>* 0 ⎞ ⎟⎟⎟⎟⎟⎟⎟⎟⎠ − *ei* cos <sup>β</sup>*<sup>i</sup> ei* sin <sup>β</sup>*<sup>i</sup>* <sup>0</sup> ⎛ ⎜⎜⎜⎜⎜⎜⎜⎜⎝ *x y z* ⎞ ⎟⎟⎟⎟⎟⎟⎟⎟⎠ + *e*<sup>2</sup> *<sup>i</sup>* + *l* 2 *<sup>i</sup>* <sup>−</sup> *<sup>e</sup>*<sup>2</sup> *i* = *<sup>x</sup>*<sup>2</sup> + *<sup>y</sup>*<sup>2</sup> + *<sup>z</sup>*<sup>2</sup> <sup>−</sup> <sup>2</sup>*xei* cos <sup>β</sup>*<sup>i</sup>* <sup>−</sup> <sup>2</sup>*yei* sin <sup>β</sup>*<sup>i</sup>* + *<sup>l</sup>* 2 *i lirTu<sup>i</sup>* = *li*( *xyz* )( sin θ*<sup>i</sup>* cos β sin θ*<sup>i</sup>* sin β*<sup>i</sup>* cos θ *i*) *T* = *li*(*x* sin θ*<sup>i</sup>* cos β*<sup>i</sup>* + *y* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* + *z* cos θ*i*) *lie<sup>T</sup> <sup>i</sup> u<sup>i</sup>* = 1*i*( *ei* cos β*<sup>i</sup> ei* sin β*<sup>i</sup>* 0 )( sin θ*<sup>i</sup>* cos β*<sup>i</sup>* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* cos θ*<sup>i</sup>* ) *T* = *li*(*ei* sin θ*<sup>i</sup>* cos2 β*<sup>i</sup>* + *ei* sin θ*<sup>i</sup>* sin<sup>2</sup> β*i*) = *eili* sin θ*<sup>i</sup> liu<sup>T</sup> <sup>i</sup> r* = *li*( sin θ*<sup>i</sup>* cos β*<sup>i</sup>* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* cos θ*<sup>i</sup>* )( *xyz* ) *T* = *li*(*x* sin θ*<sup>i</sup>* cos β*<sup>i</sup>* + *y* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* + *z* cos θ*i*) *liu<sup>T</sup> <sup>i</sup> e<sup>i</sup>* = *li*( sin θ*<sup>i</sup>* cos β*<sup>i</sup>* sin θ*<sup>i</sup>* sin β*<sup>i</sup>* cos θ*<sup>i</sup>* )( *ei* cos β*<sup>i</sup> ei* sin β*<sup>i</sup>* 0 ) *T* = *li*(*ei* sin θ*<sup>i</sup>* cos2 β*<sup>i</sup>* + *ei* sin θ*<sup>i</sup>* sin2 β*i*) = *eili* sin θ*<sup>i</sup>*

Then, Equation (4) can be extended to Equation (5).

$$\cos^2\theta + y^2 + z^2 - 2\mathbf{x}\boldsymbol{\varepsilon}\_i\cos\beta\_i - 2y\boldsymbol{\varepsilon}\_i\sin\beta\_i + l\_i^2 - 2l\_i(\mathbf{x}\cos\beta\_i + y\sin\beta\_i - \boldsymbol{\varepsilon}\_i)\sin\theta\_i - 2l\_i\boldsymbol{\varepsilon}\cos\theta\_i = 0 \tag{5}$$

Considering that Equation (5) is the function of (x,y,z), we can get Equation (6):

$$f\_i(\mathbf{x}, y, z) = \mathbf{x}^2 + y^2 + z^2 + l\_i^2 - 2(l\_i \sin \theta\_l + e\_l)\cos \beta\_l \mathbf{x} - 2(l\_i \sin \theta\_l + e\_l)\sin \beta\_l y + 2l\_i e\_i \sin \theta\_l - 2l\_i x \cos \theta\_l = 0. \tag{6}$$

Let t=(x y z), t (k) is the k-th approximation of the solution t, *f* i(t) is expanded in vector form by Taylor series to Equation (7).

$$f\_i(t) \approx f\_i(t^{(k)}) + \nabla f\_i(t^{(k)})^T (t - t^{(k)}),\tag{7}$$

.

that is:

$$f'(t) = \begin{pmatrix} \nabla f\_1(t)^T \\ \nabla f\_2(t)^T \\ \vdots \\ \nabla f\_n(t)^T \end{pmatrix} = \begin{pmatrix} \frac{\partial f\_1(t)}{\partial t\_1} & \frac{\partial f\_1(t)}{\partial t\_2} & \cdots & \frac{\partial f\_1(t)}{\partial t\_n} \\ \frac{\partial f\_2(t)}{\partial t\_1} & \frac{\partial f\_2(t)}{\partial t\_2} & \cdots & \frac{\partial f\_2(t)}{\partial t\_n} \\ \vdots & \vdots & \ddots & \vdots \\ \frac{\partial f\_n(t)}{\partial t\_1} & \frac{\partial f\_n(t)}{\partial t\_2} & \cdots & \frac{\partial f\_n(t)}{\partial t\_n} \end{pmatrix}$$

By using the Newton iterative method to solve the formula, the solution t can easily be calculated:*t* (*k*+1) = *t* (*k*) + Δ*t* (*k*).

The procedure of the MATLAB simulation is as follows.

(1) The initial values of coordinate point *<sup>t</sup>*<sup>0</sup> <sup>=</sup> *x*<sup>0</sup> *y*<sup>0</sup> *z*<sup>0</sup> *<sup>T</sup>* are assigned for the Newton iterative method;

(2) Calculate the initial value of the function according to the initial value *t*0;

(3) Calculating the initial Hesse matrix *f* <sup>0</sup>(*t*);

(4) According to the formula *f* (*t* (*k*))Δ*t* (*k*) = <sup>−</sup>*f*(*<sup>t</sup>* (*k*)), the value Δ*t* (0) can be calculated;

(5) According to the formula *t* (*k*+1) = *t* (*k*) + Δ*t* (*k*), the value *t* (1) can be obtained;

(6) Judge whether the end condition of the iteration is met. If not, return to step (2) to continue iteration.

In order to verify the accuracy of the mathematical model, according to the actual operation state and the actual situation of the measuring device, we design two sets of the angles θ, θ ∈ [−2◦, 0◦] and θ ∈ [0◦, 2◦] respectively, and generate the deflection angle according to the interval of 0.5◦. The simulation results are shown in Figure 2.

**Figure 2.** Simulation distribution of measurement length deviation of different distances of step gauge with different attitude angle θ; (**a**). θ is a negative angle, which means the normal direction of the measured surface and the starting zero plane of the step gauge is opposite; (**b**) when θ is a positive angle, the normal direction of the measured surface and the starting zero plane of the step gauge are the same.

It can be seen from Figure 2 that the maximum length error of the simulation results is less than 3 nm.

## **3. Device and Experiment**

According to the above-mentioned mathematical model, the step gauge calibration system is developed, as shown in Figure 3.

**Figure 3.** The step gauge measurement system; (**a**) actual calibration scenario of step gauge measurement system; (**b**) the critical optical system. Among this, 1 is the refractive index tracker to compensate the laser wavelength change, 2 is the plane interferometer, 3, 4 and 5 are the three-way length measuring interferometers from a single-frequency laser.

The trigger signal of the probe of CMM is used to determine the position of the measuring surface of the step gauge.

Before the experiment, it is necessary to calibrate the trigger accuracy of CMM probe. In this system, we adopt the trigger probe with a parallel leaf spring structure, and the displacement change of the probe can be measured by a Linear Variable Differential Transformer (LVDT).

The measuring principle of the probe is shown in Figure 4. Given the Uprim signal at the input end of LVDT, when the probe moves to the left side, it drives the core to move to the left, which leads to the greater voltage Uleft in the left than the voltage Uright in the right, so the output voltage Usec is displayed as a sinusoidal positive signal; otherwise, it is a sinusoidal negative signal. After shaping, subdividing and amplifying the output voltage Usec, it is converted into digital voltage signal via an 18-bit analog-to-digital conversion(AD) card first, and then converted into the moved displacement value.

**Figure 4.** Measuring principle diagram of precision displacement probe for Coordinate Measuring Machine (CMM). (**a**) Mechanical schematic diagram of the CMM probe; (**b**) circuit schematic diagram of the CMM probe; (**c**) schematic diagram of the voltage of the probe signal changing with time during the contact measurement of the probe.

In this system, the trigger sampling interval of the probe is in the linear region of the falling edge of the probe voltage signal shown in Figure 4c.

The key technology of the step gauge measurement system in realizing the high measurement accuracy is to realize the synchronous sampling and signal synthesis of the probe voltage signal of the CMM and the displacement measurement data of the laser interferometer. The probe of CMM produces the voltage signal with linear change in the measurement process, as shown in Figure 5a. Through synchronous dynamic sampling, the corresponding laser interferometer displacement measurement signal is obtained, as shown in Figure 5b. After the two signals are synthesized, the constant coordinates of the measuring points at the time of triggering are obtained, as shown in Figure 5c.

**Figure 5.** Signal synthetic diagram of probe voltage and displacement of laser interferometer: (**a**) the sampled voltage value via AD card after triggering during measurement; (**b**) the sampled laser value after triggering; (**c**) the combined value of laser value and voltage value obtained by synchronous sampling after triggering.

Before calibration, we first use the known high-level standard ball (its sphericity is 0.03 μm, its diameter is 29.98604 mm) to calibrate the diameter of the CMM probe. The two-side measurement points of the standard ball are measured 10 times, and the measurement results are shown in Table 1 below.


**Table 1.** The measurement results for calibrate the diameter of the CMM probe.

Then, the calibrated diameter of the probe can be calculated by Equation (8).

*D*probe = *L*leftLaser − *R*rightLaser − *D*sphereNom. = 713.81026 − 678.82776 − 29.98604 = 4.99646 mm. (8)

In order to verify the correctness of the mathematic model and this device, a step gauge made by the KOBA company of Germany is calibrated. The maximum deviation of measurement results between this device and the Deutscher Kalibrierdienst (DKD) device is less than 0.5 μm, as shown in Figure 6.

**Figure 6.** The variation trend of differences between Deutscher Kalibrierdienst (DKD) data and our device with measurement length.

In order to verify the measurement accuracy of this device, a second-grade standard gauge block of 500 mm was repeatedly calibrated 10 times. The deviation between the measurement results of this device and the nominal value of the gauge block is shown in Figure 7.

**Figure 7.** The trend chart of deviation of measurement result of this device from nominal value.

## **4. Measurement Uncertainty**

From this device, we can evaluate the uncertainty in measurement in accordance with the GUM (ISO, "Guide to the Expression of Uncertainty in Measurement," Switzerland, 1993, corrected and reprinted 1995). The calibration value is obtained by Equation (9) [8,18].

$$L\_{\text{stppage}} = L\_{\text{massaral}} - L\_{\text{massaral}} x\_{\text{stepage}} \Delta t\_{\text{stepage}} + \delta\_{\text{Astopyage}} + \delta\_{\text{Alaser}} + \delta\_{\text{precbining}} + \delta\_{\text{sphre}} \tag{9}$$

where:

*Lstepgage*: the measured distance of the step gauge at the reference temperature of 20 ◦C.

*Lmeasured*: the distance of the step gauge from the calculated measurement value of SMLI.

α*stepgage*: the thermal expansion coefficient of the step gauge.

Δ*tstepgage*: the difference between the step gauge temperature and the reference temperature of 20 ◦C.

δ*Astepgage*: the uncertainty of the step gauge alignment.

δ*Alaser*: the uncertainty of the laser alignment.

δPr*obing*: the reproducibility of the probing.

δ*sphere*: the uncertainty of the standard sphere calibration.

The uncertainty of the step gauge calibration is calculated by Equation (10).

$$\begin{array}{ll} \mathfrak{u}\left(\mathrm{L}\_{\mathrm{stropage}}\right)^{2} = & \mathfrak{u}\left(\mathrm{L}\_{\mathrm{massurrad}}\right)^{2} + \mathrm{L}^{2}\_{\mathrm{Nom.}} \times \mathfrak{u}\left(\mathrm{a}\_{\mathrm{stropage}}\right)^{2} \times \Delta t^{2}\_{\mathrm{stropage}} \\ & + \mathrm{L}^{2}\_{\mathrm{Nom.}} \times \mathrm{a}^{2}\_{\mathrm{stropage}} \times \mathfrak{u}\left(\Delta t\_{\mathrm{stropage}}\right)^{2} + \mathfrak{u}\left(\delta\_{\mathrm{Astropage}}\right)^{2} + \mathfrak{u}\left(\delta\_{\mathrm{Alaser}}\right)^{2} \\ & + \mathfrak{u}\left(\delta\_{\mathrm{Prabing}}\right)^{2} + \mathfrak{u}\left(\delta\_{\mathrm{sphere}}\right)^{2}. \end{array} \tag{10}$$

Table 2 shows the uncertainty components in the step gauge calibration.


**Table 2.** Uncertainty budget in step gauge calibration.

In Table 2, 1–3 are the independent sources of uncertainty, and 4–13 are dependent on thermal expansion uncertainties and have been converted to the length L.

Then, the uncertainty in measurement of the step gauge calibration using our instrument is as follows.

$$\mathcal{U} = (k \times \sqrt{0.054^2 + (0.50L)^2})\,\upmu\text{m; }L\text{-m; }k\text{-coverage factor.}$$

#### **5. Conclusions**

In order to achieve the high calibration accuracy of the step gauge, we proposed acomputational model of step gauge calibration based on SMLI technology, completed the corresponding mathematical modeling and simulation, and established a set of calibration device. We also carried out experimental verification and uncertainty analysis.

From the measurement results, we find that, although the measurement repeatability of this calibration device is stable over a short distance, it has a divergence trend over a long distance. Although the deviation of the measurement results is in accordance with the uncertainty analysis conclusion claimed by the device, the specific causes of this phenomenon will be further analyzed and studied.

In the future, we will further reduce the measurement uncertainty and improve the calibration accuracy according to the error source. It is hoped that this paper will provide a new measurement method reference for related researchers.

**Author Contributions:** Conceptualization, G.R. and X.Q.; data curation, X.C.; methodology, G.R.; software, G.R.; writing – original draft, G.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Key Research and Development Program of China, grant number: 2018YFF0212702.

**Acknowledgments:** This study was supported by the National Key Research and Development Program of China (Grant No.: 2018YFF0212702). The authors would like to thank the other members of the research team for their contributions to this study.

**Conflicts of Interest:** The authors declare no conflict of interest.

## **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
