Open Access
Issue
Mechanics & Industry
Volume 22, 2021
Article Number 48
Number of page(s) 12
DOI https://doi.org/10.1051/meca/2021047
Published online 10 December 2021

© Y. Li et al., Published by EDP Sciences 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

High-performance motion systems, such as the wafer stage in the semiconductor industry continuously, call for higher throughput and accuracy, resulting in more challenges for the design in next-generation motion stages [1]. Meanwhile, these increasingly stringent performance requirements result in a desire for increased acceleration and control performance [2]. The increased acceleration implies a more lightweight system design or higher actuation forces, causing that the flexible dynamic behavior is more prominent or even appear at lower frequencies [37]. A higher control bandwidth helps improve the control performance of the system, but meanwhile, it will aggravate the flexibility characteristics. The non-negligible flexibility will lead to the displacement gap between the point of interest and the point of control to become larger. The feedback displacement obtained based on the rigid-body assumption will not meet the system control requirements [8]. The point-of-control (POC) displacement as the input of the controller can be the displacement of the measurement point or the displacement calculated by measurement points. The point of interest is the system performance point. Its control performance determines the system performance. So, the POI displacement should be essentially used as the input of the controller. However, due to the various constraints of the system structure, the POI displacement cannot be used as the input of the system feedback. In the lithograph system [911], as shown in Figure 1, the sensors are located at the edge of the wafer stage due to the structural design. The points of interest are located in the exposure area of the wafer stage [9,12], and the POC displacements are calculated by measurement points. Therefore, it is necessary to propose a method to improve the POI displacement estimation accuracy in ultra-precise motion systems with flexibility.

Common methods of POI displacement estimation mainly include model-based and data-based methods. The model-based methods mainly estimate POI displacement by constructing observers [1315]. Observers usually include the undisturbed observers [13,14] and the disturbance observers [1517]. The undisturbed observer only consists of the controlled object itself, while the disturbance observer is composed of the controlled object and the disturbance model. In reference [13], the importance of the POI displacement estimation using the observer was evaluated from the positioning accuracy of the system through simulation analysis; the paper [14] analyzes the importance of the POI displacement estimation using the observer from the comparison of the magnitude of the flexibility of the closed-loop transfer function by building a single-degree-of-freedom experimental setup and introducing active control; the paper [15] illustrates the effectiveness of the POI displacement estimation using observer by comparing the amount of output information without the observer and with the observer. The advantage of this method is that the internal state change mechanism of the system can be understood through the observer construction process. Simultaneously, because it is highly dependent on the model, the accuracy of the model determines the accuracy of the POI displacement estimation. Usually, the disturbance model is unpredictable and stochastic. It is mainly constructed by the equivalent method based on engineering experience, such as the airflow disturbance is equivalent to Brownian motion. However, the deviation caused by equivalence, unpredictability, and randomness can lead to the inaccuracy of the disturbance observer model and thus affect the accuracy of the POI displacement estimation. The modeling methods of the controlled objects mainly include the finite element modeling method [16] and the system identification method [14,18]. For the ultra-precise motion system with flexible characteristics, the complexity of the model increases, and the difficulty of identification increases accordingly. The identification model acquired based on the existing sensors and the corresponding layout information can no longer meet the actual control requirements. For example, in the literature [18], 16 additional local sensors are added to obtain more details of the model, improving the identification accuracy.

Compared with the model-based approach, the data-based approach starts directly from the available data of the system. This data itself contains information about the structural changes, unmodeled dynamics and unknown disturbances of the system. This approach gets rid of the dependence on the controlled object model. Thus, the robustness of the estimation accuracy can be improved. In the literatures [10,19,20], the layout of the existing limited number of sensors is optimized, the linear least-squares method is used to obtain the mapping matrix between the measured and desired points, and the cloud map of estimation errors is reduced from tens of nanometers to the hundred picometer level by simulation analysis. The goal of the investigated literature is to propose an enlightening mathematical method, but there is a certain distance from practical engineering applications. Similar to the model-based method, there will be incomplete sample data for high-precision systems with prominent flexibility in practice. In the literature [10], the amount of sensor output information increases by optimizing the sensor layout, but for example, the measurement sensors in the wafer system mainly include laser interferometers or plane gratings [21]. Limited by the installation space and the volume of the sensor itself, these sensors can only be located on the edge of the wafer stage, and the number of the sensors cannot be added arbitrarily. When the flexible characteristics of the structure are continuously excited, the dynamic deformation of the structure becomes more complicated, and the amount of information transmitted by the high-precision displacement sensor is not enough to reflect the POI displacement.

Therefore, whether based on model-based [18] or data-based approach [10,19], for ultra-precise motion systems with prominent flexibility characteristics, it is necessary to increase the output of the model information, ensuring the accuracy of model identification and the completeness of model output data. The local compact vibration sensors attached to the surface of the structure can be added to obtain structural flexibility information. Acceleration and resistance strain gauges are often used for vibration. For example, in references [22,23], the specific points in the measurement bridge where the sensors cannot be installed are used to obtain the vibration information of the measured object by fusing acceleration sensors and resistive strain gauges, and the laser displacement sensors are used as reference sensors to achieve the estimation of the unmeasurable points. In reference [24], to improve the control accuracy, the motion displacement of a flexible beam is obtained through the resistance strain gauge, acceleration sensor and gyroscope as the system feedback inputs. However, the local sensors in such applications exhibit low measurement accuracy and are not suitable for the ultra-precise motion systems. The optical sensor as a local sensor is with high-resolution characteristics in addition to being small and basically free from structural constraints. The functional optical fiber sensor combined with the phase modulation principle has a resolution of 0.01 µε [25,26].

According to the above review, this paper proposes a POI displacement estimation method for ultra-precision motion systems based on the functional optical fiber sensors and the phase modulation principle. The main contents are as follows: Section 2 presents the specific measurement principle based on this measurement scheme and the layout optimization of the functional optical fiber sensor. Then, a POI displacement estimation model is proposed using a time-series dynamic neural network based on the data approach. In Section 3, the optimal optical fiber sensor positions and lengths corresponding to the critical modes are determined using the simulation data of the controlled object by the proposed sensor layout optimization method. In Section 4, based on the above contents, an experimental setup is designed, and a POI displacement estimation model is constructed based on the optical fiber sensor output data using a data-based approach. The validity of the POI displacement estimation model based on the optical fiber sensor is proved through comparative analysis.

thumbnail Fig. 1

POC and POI in the schematic of the flexible wafer-stage system.

2 Principle

2.1 Optical fiber sensor selection and measurement principle

The design of the measurement scheme based on the optical fiber sensor revolves around two aspects: the selection of the sensor and the measurement principle.

2.1.1 Sensor selection

In the ultra-precision motion system, the positioning accuracy is at the nanometer level. Therefore, in terms of sensor selection, a functional optic fiber sensor is selected from a high-resolution perspective [25,26]. The functional sensor uses the characteristics of the optical fiber itself as a sensitive element to modulate the light transmitted within the optical fiber, making the intensity, phase, frequency or polarization state of the transmitted light change, and then the modulated signal is demodulated to obtain the measured signal [27,28].

2.1.2 Measurement principle

The phase modulation method is widely used in optical fiber sensing fields such as fiber-optic sensors and fiber-optic hydrophones for its simple hardware structure and easy multiplexing [29]. The basic sensing mechanism of the phase-modulated optic fiber sensor is to change the phase of the light wave propagating in the optical fiber by the action of the measured energy field, and then use interferometric techniques to convert the phase transformation into a change in light intensity to detect the physical quantity to be measured.

When the measured object undergoes flexible vibration, the length and longitudinal direction (axial direction) of the measured object are subjected to mechanical stress. In the modal coordinate, it is characterized by the following relation: { Δ L y ( t ) = i = 1 m ϕ y i q y i ( t ) Δ L z ( t ) = i = 1 m ϕ z i q z i ( t ) Δ L x ( t ) = i = 1 m ϕ x i q x i ( t ) (1)

where ΔL(x,y,z)(t) the change in three directions; x is the longitudinal direction, y and z is the transverse direction. Moreover, the corresponding ith eigenvector and modal displacement are ϕ(xi,yi,zi), q(xi,yi,zi); m denotes the measured object has m modes.

When an optical fiber receives mechanical stresses in the longitudinal (axial) direction, the length, core diameter, and core refractive index of the optical fiber will change, and these transformations will lead to the change in the phase of the optical wave in the optical fiber. The specific expression is [28]: Δ φ = β Δ L + Δ β L = β L Δ L L + L β a Δ a + L β n Δ n , (2)where β = 2π/λ is the propagation constant of light waves in optical fiber, λ = λ0/n is the propagation wavelength of light waves in optical fiber, λ0 is the propagation wavelength of light waves in vacuum, n is the refractive index of light waves, and a is the radius of the optical fiber core. In this formula, the first term is the strain effect, which indicates the phase delay caused by the change of the optical fiber; the second term is the Poisson effect, which is the phase delay caused by the change of the optical fiber radius; the third term is the elasto-optical effect, which represents the phase delay caused by the change of the induced refractive index.

For the single-mode fiber, the strain optical matrix and Hooke's law are introduced into equation (2), and the general form of the single-mode fiber can be obtained [28]: Δ φ = n 2 π λ 0 L ( 1 1 2 n 2 ( ( 1-v ) P 12 v P 11 ) v ( λ 4 a ) 2 1 ( 1 [ λ 4 a ] 2 ) 0.5 1 2 a ) ϵ 3 , (3)where P11 and P12 are the strain optical constants; ε3=ΔL/L is the longitudinal strain of the optical fiber; v=|ε13| is the Poisson's ratio, and ε1 is the transverse strain of the optical fiber.

Combining equations (1) and (3), the relationship between the phase change and the eigenvalue of the object and the modal displacement can be obtained: Δ φ = n k 0 ( 1 1 2 n 2 ( ( 1 - v ) P 12 v P 11 ) v ( λ 4 a ) 2 1 ( 1 [ λ 4 a ] 2 ) 0.5 1 2 a ) i = 1 m ϕ x i q x i ( t ) . (4)

Interferometric techniques, including homodyne interferometry and heterodyne interferometry, are often used to obtain the phase change. The homodyne interferometry technique is susceptible to low-frequency electronic noise, while the heterodyne interferometry technology can avoid the interference of low-frequency electronic noise through modulation and demodulation. From the perspective of anti-interference, the heterodyne interference is selected to obtain the phase change of the optical fiber.

The two coherent light beams are with the frequencies f1 and f2, the light vector amplitudes E1 and E2, and the initial phases θ1 and θ2. When the object vibration causes the optical fiber to be stretched or compressed in the measurement path, the optical path difference is generated in the optical fiber. It means that the phase change of the light vector Δφ (t) is introduced into the measurement light. Then the light vectors of the measured light and the reference light can be expressed as: { E r e f = E 1 e i ( θ 1 + 2 π f 1 t ) E m e a = E 2 e i ( θ 2 + 2 π f 2 t + Δ φ ( t ) ) . (5)

Two light waves interfere in the optical fiber, and the light intensity at a point in the interference field can be expressed as: I = ( E r e f + E m e a ) ( E r e f + E m e a ) = E 1 2 + E 2 2 + 2 E 1 E 2 cos ( θ 1 θ 2 Δ φ ( t ) + 2 π Δ f t ) , (6)where Δf=f1-f2.

Combining equations (4) and (6), the relations among the POI eigenvectors, the phase change and the optical intensity of the fiber light wave can be seen: ϕ x i Δ φ I , (7)

As mentioned above, the phase change is proportional to the strain eigenvector and the modal displacement, meaning that the optical fiber change can characterize the flexible properties of the object.

In summary of the above measurement principles, a frequency synthesizer and two acousto-optic modulators are utilized to complete the acquisition of signals from the optical fiber sensors based on the principle of heterodyne interferometry [30]. First, the light wave from the laser is divided into two by an optical splitter, and the two signals with a certain frequency difference are formed through the frequency synthesizer and the acousto-optic modulator. One is the reference path, and the other is the measurement path. After the reference signal interferes with the measurement signal, the phase information is obtained through the phase board. The specific measurement principle is shown in Figure 2.

thumbnail Fig. 2

Schematic of the functional optical fiber sensor measurement principle.

2.2 Layout optimization of optical fiber sensor

Unlike the point measurement of common sensors, the functional optical fiber sensor measurement is a line measurement. In modal coordinates, the flexible displacement of the analyzed object can be described as the m-order modal superposition. It is necessary to find several critical modes contributing the most to the flexible displacement. Therefore, a sensor optimal layout method with this measurement feature, including optimizing the location and length of the sensor based on critical modes of the measured object in the interest area, is proposed. The details are as follows:

(1) Judgment of critical modes in the interest area

When the structural flexibility is excited, every POI exhibits different dynamic characteristics. The critical modes are determined in the interest area with the orientation of displaying the POI displacement.

The proportion of the motor-driven input corresponding to each mode is set under certain excitation condition: ( p f ) i = f i i = 1 m f i , (8)where, for each mode, fi represents the ith modal force, and (pf)i is the percentage of the ith modal force.

It is assumed that the longitudinal deformation of the POI in the interest area is ΔL. The interest area and frequency bands are pre-defined. The deformation corresponding to the critical mode is determined under the conditions of the specific excitation conditions in equations (9) and (10). These shows that the modal order can be determined by ΔLi, further clarifying the critical modes. Δ L i = l ϕ P O I ( x , y ) i { b m i } { u } 2 ζ i ω i 2 = l ϕ P O I ( x , y ) i [ Φ i ] T [ B ] { u } 2 ζ i ω i 2 , (9) Δ L i ( p f ) i * l ϕ P O I ( x , y ) i 2 ζ i ω i 2 , (10)where [Φi] T[B] {u} is the ith modal force, [Φi] T is the ith transposition of the eigenvalue of the applied force point, [B] is the force distribution matrix, {u} is the input force amplitude; ζi and ωi are the ith damping ratio and eigenvalue, respectively; the sum of the ϕ.POI is determined by the position (x, y) and the length l of the optical fiber sensor.

(2) Determination of optimization goal

There are nf alternative optimization locations; each has mm nodes in the discrete structure. The square of the strain eigenvalue is taken as the optimization target from a strain energy perspective. And considering the effect of the force, the optimal position and length ordering of the critical modes are obtained: E i = max ( ( p f ) i * E i ( x , y , l ) ) = max [ ( ( p f ) i * j j = 1 m m ( - 1 ) i * ( ϕ ϕ s ( x j j , y j j ) ) 2 ) ] , (11)where Ei is the ith optimization target; l is the length of the optical fiber sensor.

In the measurement, the measurement result of each sensor depends not only on the sensor layout but also on the measurement properties of the sensor itself. The measurement variance σ j s of each sensor is taken as a correction factor to be added to the optimization target. If there are nc sensors, the correction coefficient of the jth sensor is as follows: α j = ( 1 σ s ) / j = 1 n c ( 1 σ s ) , (12) γ j = α j max { α 1 , . , .. , α j , . , . , α n } , (13)where γj is the correction coefficient of the jth sensor.

Setting the proportion of the position and length of the sensor as ŋ and the proportion of the property of the sensor itself as (1-ŋ), the final correction coefficient of the sensor where the ith critical mode is located is: κ i = ŋ + ( 1 ŋ ) γ j , (14)

The corrected optimization goal is: E ̂ i ( x , y , l ) = κ i E i ( x , y , l ) . (15)

In order to reduce the redundancy of the output information at the location of each sensor, the correlation is used as a judgment index. The weak correlation value of less than 0.4 is chosen as the threshold value. Assuming that the eigenvalue corresponding to the first-ranked critical mode is Q1, the eigenvalue Q 2 corresponding to the location and length of the second-ranked critical mode is limited by Q 1. Then the location of the second critical mode that has been optimized is screened for the second optimization, as following: ρ 12 = C o v ( Q 1 , Q 2 ) Q 1 Q 2 < 0.4 . (16)

And so on, for the subsequent ith order critical modes, the following conditions should be to be satisfied: ρ 1 i & ρ 2 i & ρ 3 i & ... < 0.4 . (17)

The optimal location and length of the optical fiber sensor corresponding to the critical modes can be obtained in turn through the above optimization.

The signal of the optimal sensor corresponding to each order of the flexible mode inevitably contains the information of the other orders, so the information transmitted by the sensors has a certain redundancy. On the other hand, the multiple sensors can reduce the uncertainty of measurement. Differential entropy can quantify the total uncertainty in the probability distribution of continuous random variables. The smaller the value is, the smaller the uncertainty is. Therefore, the following is proved from the perspective of differential entropy.

Assuming that ΔL conforms to a Gaussian distribution, the possibility that nc(nc≥2) sensors will occur simultaneously is to multiply the probability of nc events.

Let Ni(ΔL) be Gaussian PDFs with arbitrary mean µi and standard deviation σi, { N 1 ( Δ L , μ 1 , σ 1 ) = 1 σ 1 2 π e ( Δ L μ 1 ) 2 2 σ 1 2 . . . N n c ( Δ L , μ n , σ n ) = 1 σ n 2 π e ( Δ L μ n ) 2 2 σ n 2 . (18)

Their product is: N(ΔL,(μ1,μ1),(μ2,μ2),.,.,.,(μn,μn))=i=1n1μi2πμi2e(ΔLμi)22μi2=1(2ππ)n2i=1nμσi2eΔL22εj=1nc(ΔLji=1nc(jnc)μσi2)j=1nc(i=1nc(jnc)μσi2)+j=1nc(ΔLj2i=1nc(jnc)μσi2)j=1nc(i=1nc(jnc)μσi2)2i=1ncμσi2j=1nc(i=1nc(jnc)μσi2),(19)

Then the fusion standard deviation can be obtained as follows: σ f u s e d 2 = 2 i = 1 n c σ i 2 j = 1 n c ( i = 1 n c ( j n ) σ i 2 ) = σ j 2 - σ j 4 ( i = 1 n c ( j n c ) σ i 2 ) j = 1 n c ( i = 1 n c ( j n c ) σ i 2 ) . (20)

From the above derivation, it is clear that the variance of the fusion is smaller than that of any single sensor: σ f u s e d 2 < σ j 2 . (21)

In terms of the differential entropy, it is derived as follows: H ( x ) = N ( Δ L ) lnN ( Δ L ) d Δ L = ln ( 2 π σ ) σ 2 π 2 σ e ( Δ L μ ) 2 2 σ 2 d ( Δ L μ 2 σ ) + 1 σ 2 π 2 σ e ( Δ L μ ) 2 2 σ 2 ( Δ L μ ) 2 2 σ 2 d ( Δ L μ 2 σ ) , (22)

Let z = Δ L μ 2 σ : - e z 2 d z = π . (23)

Then the differential entropy becomes: H ( x ) = ln ( 2 π σ ) σ 2 π 2 σ - e z 2 d z + 1 σ 2 π 2 σ - e z 2 z 2 d z = 1 2 ( ln ( 2 π σ 2 ) + 1 ) . (24)

Combining equations (21) and (24), it can be seen that multiple sensors can effectively reduce the uncertainty of measurement.

2.3 POI displacement estimation model

The displacement induced by flexibility is not instantaneous but a process of constant decay. Therefore, the current moment displacement is related to the displacement of the previous moments. Therefore, the constructed model needs to have dynamic properties. ANN is a non-linear regression model and has been widely used for regression [30,31]. A neural network structure is constructed using the input layer with the input delays, a Gate Recurrent Unit (GRU) layer [32] and the output layer, as shown in Figure 3. The input of the input layer includes not only the current inputs but also the past inputs, while the GRU layer includes not only the output of the input layer but also the output of the hidden layer at the previous moments. These two layers involve the previous inputs, making the model dynamic. The POI model can be expressed: y ̂ t = f ( x t , x t 1 , ... , x t d e l a y , w 1 , w 2 , w 3 , B 1 , B 2 , B 3 ) , (25)where w1, w2, w3 and B1, B2, B3 represent the network weight vectors and the biases of the input-hidden layer, the hidden layer, and the hidden-output layer, respectively; y ̂ t stand for the estimated POI displacement. The relative error (EH2) is applied to estimate the accuracy of the estimated displacements intuitively. E H 2 = y r e f y ̂ 2 y r e f 2 × 100 % , (26)where yref is the reference POI displacement.

thumbnail Fig. 3

Block diagram of the neural network structure.

3 Simulation analysis of functional optical fiber sensor optimization

The simulation object, shown in Figure 4, is taken as the analysis object. In this figure, the central area is the interest area, and the flexible frequencies of interest within the 2500 Hz are the analysis focus. Through simulation analysis, the frequencies involved are 313 Hz, 851 Hz, 1065 Hz, 1469 Hz, 1641 Hz, and the corresponding order ranges from 7th to 11th.

Excluding the drive area of the analyzed object and other measurement areas, there are a total of three areas where the optical fiber sensors can be placed, each area being discretized into six locations, as shown in Figure 5. According to the constraints of the lengths and locations, there are 1275 choices to be selected for each location and a total of 19125 choices for the three regions.

According to the optimization method in Section 2, the critical modes of the desired region are first determined. Assuming a constant force input of 1 N, the critical modes of the interest area are determined and ordered as 7th, 8th, 11th, 9th, 10th. The critical modes become the 7th mode and the 11th mode by adding the force constraint.

Since the measured sensor output information is not available, the correction coefficients of the sensors are set as: { κ 1 = 1 κ 2 = 1 . (27)

Then, the optimization target and correlation constraints are set based on the correlation. The optimization results and corresponding strain eigenvalue of the optical fiber sensors can be obtained for the 7th mode and the 11th mode, as shown in Table 1 and Figure 6, respectively.

thumbnail Fig. 4

Analyzed object with the area of interest.

Table 1

Optical fiber sensor optimization result.

thumbnail Fig. 5

Optical fiber sensor optimization − alternative areas.

thumbnail Fig. 6

Critical mode ordering: (a) under the force-free constraints; (b) under the force constraints.

4 Experimental analysis

To verify the role of the optical fiber sensor in estimating the POI displacement, an experimental device, as shown in Figure 7, was built according to the optical fiber measurement principle diagram in Section 2 and the sensor optimization results in Section 3. This experimental device mainly includes a set of capacitive sensors, optical fiber sensors and drive motors. The position and number of three sets of the capacitive sensors are assumed to be unchangeable. Sensors 1# and 3# are used for closed-loop control, and their solution results based on the rigid-body assumption are equivalent to the POC displacement; Sensor 2# is a reference sensor, and its measurement result is equivalent to the POI displacement. The two sets of optical fiber sensors correspond to the 7th and 11th critical modes, respectively. Their measuring surface is on the same side as the measuring surface of the capacitive sensors. There are three sets of drive motors in total. Motors 1# and 3# are used for closed-loop control, and motor 2# is an excitation motor, which is used to excite the analysis object when the system is under closed-loop control. The data acquisition for the capacitance sensors and the optical fiber sensors is performed using an AD board and a Zygo4104c board with a sampling period of 5000 Hz, respectively; the drive motor is driven by a linear driver VAREDAN-LA525.

thumbnail Fig. 7

Experimental setup.

thumbnail Fig. 8

Analysis of optical fiber sensor characteristics: (a) sensor 1#; (b) sensor 2#.

4.1 Correction of sensor optimization results

According to the measurement principle diagram in Section 2, the two optical fiber sensor signals to be arranged were recorded without the active disturbance, as shown in Figure 8. The optical fiber sensors have high sensitivity and are susceptible to external disturbances, causing that the signal of the sensors have low-frequency drift. This low-frequency drift over time cannot characterize the measurement capability of the sensor. Therefore, the low-frequency drift needs to be eliminated. The signal was removed baseline drift based on the wavelet packet method [33,34] with the wavelet basis function ‘Db44’ and the six-layer decomposition. After removing the low-frequency drift, it can be observed in the histogram that the signal is normally distributed. So, { σ 1 = 12.78 σ 2 = 8.21 σ f u s e d = 7.19 { H ( σ 1 ) = 3.97 H ( σ 2 ) = 3.52 H ( σ f u s e d ) = 3.39 . (28)

Based on the variances of the signals, it is known that { γ 1 = 0.6424 γ 2 = 1 . (29)

The sensor properties are set to 20%, then { κ 1 = 0.9285 κ 2 = 1 .(30)

Optical fiber sensor #1 corresponds to the 7th critical mode and Optical fiber sensor #2 corresponds to the 11th critical mode. Then the corresponding eigenvalues after the correction are:

Comparing Figures 9 and 10, there is a certain effect on the magnitude of sensor #1 under the constraints of the sensor properties. However, it does not change the performance of the most critical modes on this sensor. If the sensor properties affect the information acquisition of the mode, it is necessary to adjust the corresponding relationship with the sensors according to the order of the modes or to replace the sensor with the awful properties.

thumbnail Fig. 9

Optimal result of optical fiber sensor: the optimal strain vector.

thumbnail Fig. 10

Fiber optimization result: optimal strain vector.

4.2 Signal analysis

The signals of the optical fiber sensors and the capacitive sensors are acquired under closed-loop control conditions from the host computer under a multi-sine sweep signal active disturbance outputting from motor 2#. The spectrogram of the signals is illustrated in Figure 11.

Based on the extracted signals, the following conclusions are drawn.

  • The actual analyzed object and the simulated object cannot be the same, so there is a certain gap between the obtained frequency. But this gap is not enough to affect the modal ordering for this structure in this paper.

  • In the amplitude spectral density plot, the capacitive sensor output information contains flexible information; therefore, the POC displacement obtained by capacitive sensor #1 and capacitive sensor #2 based on the rigid-body assumption is not equal to the POI displacement of capacitive sensor #2.

  • Compared to optical fiber sensor #2, optical fiber sensor #1 has a higher amplitude in the low-frequency band and lower amplitude in the high-frequency band, being consistent with the optimized results of the optical fiber sensors in Section 3.

  • The signals from the optical fiber sensors include all flexible information outputted by capacitive sensor #2, while the flexible information output from capacitive sensor #2 is not reflected in the capacitive sensors #1 and #3. In addition to the flexibility information of the measured object, the optical fiber sensors also contain additional flexible information belonging to the support base, being analyzed by the accelerometers, as shown in Figure 12. Compared to the more random signals in capacitive sensors 1# and 2#, these deterministic signals can be easily separated from the original signals.

thumbnail Fig. 11

Spectrogram of the signals: (a) capacitive sensors; (b) optical fiber sensors.

thumbnail Fig. 12

Analysis of controlled objects by accelerometers: (a) experiment setup; (b) spectrogram of the signals.

4.3 POI displacement estimation model

A POI displacement estimation model was constructed using the neural network structure mentioned in Section 2. This network structure includes five inputs with 60 delays and 512 GRU neurons. The five inputs are capacitive sensors #1 and #3, optical fiber sensors #1 and #2, and the exciting force. The excitation force is a complicated multi-sine excitation force with an excitation frequency of 10 Hz∼2000 Hz. There are 10000 data sets collected in the experiment, being randomly separated into training and validation sets as 70:30 ratios. Moreover, five other group data sets with 1000 data respectively are taken as the final test sample sets. The data are normalized into a range [-1, 1] to accelerate the convergence rate and improve the accuracy [35]. The training error is 0.9%, while the five test errors are between 3.84% and 4.69%; based on the rigid body assumption estimation method, the results of the five test groups are between 10.22% and 10.86%. The average error of the five tests is shown in Figure 13. Therefore, the average estimation error of the proposed optical fiber sensor-based POI displacement estimation method is reduced by more than 61% compared with the rigid-body assumption-based estimation method.

thumbnail Fig. 13

Average error of the five group tests.

5 Conclusion

This paper proposes a POI displacement estimation method for ultra-precision motion systems based on the functional optical fiber sensor and the phase modulation principle. The core part of the paper is the measurement scheme and the layout optimization analysis based on the functional optical fiber sensor. To verify the effectiveness of the proposed estimation method, an experimental setup with two sets of measurement sensors, including capacitive sensors and optical fiber sensors, is designed. The following conclusions are drawn from the experimental analysis and the results of the POI displacement estimation based on the time-series dynamic model: (1) The analysis of the frequency domain plot of the test data shows that the output information of the optical fiber sensor contains more information about the flexible characteristics of the analyzed object than the output information of the capacitive sensors. (2) Through the estimation of the test data, the five groups of the test errors of the POI displacement estimation method using the optical fiber sensors proposed in this paper is between 3.84% and 4.69%, while the errors are between 10.22% and 10.86% using the rigid-body assumption estimation method. Based on the latter error, the estimation error is reduced by about 61% compared with the rigid-body hypothesis-based estimation method.

This paper does not separate out the signals belonging to the controlled object in the optical fiber sensors when constructing the POI displacement estimation model. However, the valid and deterministic data can make the POI displacement estimation model built by the data-based method have better generalization performance and estimation capability. In addition, there is not any detailed analysis of the model construction, such as the relationship between model complexity and model accuracy, and model generalization ability. Therefore, the specific model structure needs to be further optimized. In the next research work, we will focus on extracting the effective signal and constructing the best POI displacement estimation model to lay the foundation for subsequent applications.

Disclosures

The authors declare that there is no conflict of interest.

References

  1. T. Oomen et al., Next-generation wafer stage motion control: Connecting system identification and robust control, in: 2012 American Control Conference (ACC), Montreal, QC 2012, pp. 2455–2460 [CrossRef] [Google Scholar]
  2. M.J.C. Ronde, Feedforward control for lightweight motion systems, PhD thesis, Eindhoven University of Technology, 2014 [Google Scholar]
  3. T. Oomen et al., Connecting system identification and robust control for next-generation motion control of a wafer stage, IEEE Trans. Control Syst. Technol. 22, 102–118 (2013) [Google Scholar]
  4. T. Oomen et al., Inferential motion control: Identification and robust control framework for positioning an unmeasurable point of interest, IEEE Trans. Control Syst. Technol. 23, 1602–1610 (2014) [Google Scholar]
  5. T. Oomen et al., Inferential motion control: Identification and robust control with unmeasured performance variables, in: 2011 50th IEEE Conference on Decision and Control and European Control Conference, Orlando, FL 2011, pp. 964–969 [CrossRef] [Google Scholar]
  6. M. Dorosti et al., Finite element model reduction and model updating of structures for control, IFAC Proc. 47, 4517–4522 (2014) [CrossRef] [Google Scholar]
  7. G. Van Baars, Getting high tech systems in shape and fit for the future, Mikronie7k 52, 30–36 (2011) [Google Scholar]
  8. Y. Jiang et al., Coupled dynamic modeling and analysis of ultra-precision 6-DOF fine stage, Proc. Chin. Soc. Electr. Eng. 30, 5451–5457 (2014) [Google Scholar]
  9. T. Oomen, E. Grassens, F. Hendriks, Inferential motion control: Identification and robust control framework for positioning an unmeasurable point of interest, IEEE Trans. Control Syst. Technol. 23, 1602–1610 (2014) [Google Scholar]
  10. J.G. Vogel et al., Estimation of the deformation of a plate using a limited number of sensors, in: ASPE Spring Topical Meeting MIT Laboratory for Manufacturing and Productivity Annual Summit, Cambridge, Massachusetts, 2013, pp. 21–23 [Google Scholar]
  11. R.A.C.M. Beerens et al., Positioning system, lithographic apparatus and device manufacturing method, U.S. Patent, 9, 383, 659[P], 2016 [Google Scholar]
  12. T. Castenmiller et al., Towards ultimate optical lithography with NXT: 1950i dual stage immersion platform, Proc. SPIE 7640, Optical Microlithography XXIII, 76401N (2010) [Google Scholar]
  13. K.W. Verkerk, H. Butler, P.P.J. Van den Bosch, Improved disturbance rejection for high precision systems through estimation of the flexible modes[C], 2015 IEEE Conference on Control Applications (CCA), 1191–1196 (2015) [CrossRef] [Google Scholar]
  14. R. Voorhoeve et al., Estimating structural deformations for inferential control: a disturbance observer approach, IFACPapersOnLine 49, 642–648 (2016) [Google Scholar]
  15. Mooren, et al., Compensating quasi-static disturbances for inferential control: an observer-based approach applied to a wafer stage, in: Proceedings of the 2018 IEEJ International Workshop on Sensing, Actuation, Motion Control, and Optimization Tokyo, Japan, 2018 [Google Scholar]
  16. R. de Rozario et al., Identification of position-dependent mechanical systems, in: 35th Benelux Meeting on Systems and Control, Soesterberg, Netherlands, 2016 pp. 22–24 [Google Scholar]
  17. M. Dorosti, Reduced-other model updating for prediction of performance variables in mechanical structures, PhD thesis, Eindhoven University of Technology, 2017 [Google Scholar]
  18. R. Voorhoeve et al., Identifying position-dependent mechanical systems: A modal approach with applications to wafer stage control, arXiv preprint arXiv:1807.06942 (2018) [Google Scholar]
  19. S.H.J. Van der Kleij, Shape estimation of a compliant wafer chuck in lithography systems, Master thesis, Delft University of Technology, 2014 [Google Scholar]
  20. A. Koevoets et al., Thermal-elastic compensation models for position control, in: Proceedings of the 24th ASPE Annual Meeting, United States, 2009 [Google Scholar]
  21. T. Castenmiller, M.F.V.D., Towards ultimate optical lithography with NXT: 1950i dual stage immersion platform, Proc. SPIE 7640, 511–519 (2010) [Google Scholar]
  22. S. Cho, C. Yun, S. Sim, Displacement estimation of bridge structures using data fusion of acceleration and strain measurement incorporating finite element model, Smart Struct. Syst. 15, 645–663 (2015) [CrossRef] [Google Scholar]
  23. S. Cho et al., Reference-free displacement estimation of bridges using kalman filter-based multimetric data fusion, J. Sens, 1–9 (2016) [Google Scholar]
  24. T. Örtel, J.F. Wagner, F. Saupe, Integrated motion measurement illustrated by a cantilever beam, Mech. Syst. Signal Process. 34, 131–145 (2013) [CrossRef] [Google Scholar]
  25. D.A. Krohn, T. MacDougall, A. Mendez, Fiber Optic Sensors: Fundamentals and Applications, SPIE Press, Bellingham, WA, 2014 [CrossRef] [Google Scholar]
  26. R. Di Sante, Fiber optic sensors for structural health monitoring of aircraft composite structures: recent advances and applications, Sensors 15, 18666–18713 (2015) [CrossRef] [PubMed] [Google Scholar]
  27. K.T.V. Grattan, B.T. Meggitt, L.S. Grattan, Optical Fiber Sensor Technology, Springer, 2000 [Google Scholar]
  28. Y. Zhao, Principles and Application Technology of Optical Fiber Sensors, Tsinghua University Press, Beijing, 2007 [Google Scholar]
  29. A. Dandridge, A.B. Tveten, T.G. Giallorenzi, Homodyne demodulation scheme for fiber optic sensors using phase generated carrier, IEEE J. Quant. Electr. 18, 1647–1653 (1982) [CrossRef] [Google Scholar]
  30. Y. Yu, Research on Zygo Interferometer and Its Application, Master thesis, Tianjin University, 2005 [Google Scholar]
  31. A. Prieto, B. Prieto, E.M. Ortigosa et al., Neural networks: An overview of early research, current frameworks and new challenges, Neurocomputing 214, 242–268 (2016) [CrossRef] [Google Scholar]
  32. K. Cho et al., Learning phrase representations using RNN encoder-decoder for statistical machine translation, arXiv preprint arXiv:1406.1078 (2014) [Google Scholar]
  33. I. Daubechies, Ten Lectures on Wavelets, Society for Industrial and Applied Mat, 1992 [CrossRef] [Google Scholar]
  34. R. Yan, Base wavelet selection criteria for non-stationary vibration analysis in bearing health diagnosis, PhD thesis, University of Massachusetts Amherst, 2014 [Google Scholar]
  35. J. Sola, J. Sevilla, Importance of input data normalization for the application of neural networks to complex industrial problems, IEEE Trans. Nucl. Sci. 44, 1464–1468 (1997) [CrossRef] [Google Scholar]

Cite this article as: Y. Li, M. Zhang, Y. Zhu, Research on the estimation method of the point-of-interest (POI) displacement for ultra-precision flexible motion system based on functional optical fiber sensor, Mechanics & Industry 22, 48 (2021)

All Tables

Table 1

Optical fiber sensor optimization result.

All Figures

thumbnail Fig. 1

POC and POI in the schematic of the flexible wafer-stage system.

In the text
thumbnail Fig. 2

Schematic of the functional optical fiber sensor measurement principle.

In the text
thumbnail Fig. 3

Block diagram of the neural network structure.

In the text
thumbnail Fig. 4

Analyzed object with the area of interest.

In the text
thumbnail Fig. 5

Optical fiber sensor optimization − alternative areas.

In the text
thumbnail Fig. 6

Critical mode ordering: (a) under the force-free constraints; (b) under the force constraints.

In the text
thumbnail Fig. 7

Experimental setup.

In the text
thumbnail Fig. 8

Analysis of optical fiber sensor characteristics: (a) sensor 1#; (b) sensor 2#.

In the text
thumbnail Fig. 9

Optimal result of optical fiber sensor: the optimal strain vector.

In the text
thumbnail Fig. 10

Fiber optimization result: optimal strain vector.

In the text
thumbnail Fig. 11

Spectrogram of the signals: (a) capacitive sensors; (b) optical fiber sensors.

In the text
thumbnail Fig. 12

Analysis of controlled objects by accelerometers: (a) experiment setup; (b) spectrogram of the signals.

In the text
thumbnail Fig. 13

Average error of the five group tests.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.