Investigation of X and Y Conﬁguration Modal and Dynamic Response to Velocity Excitation of the Nanometer Resolution Linear Servo Motor Stage with Quasi-Industrial Guiding System in Quasi-Stable State

: Research institutions and industrial enterprises demand high accuracy and precision positioning systems to fulﬁl cutting edge requirements of up-to-date technological processes in the ﬁeld of metrology and optical fabrication. Linear motor system design with high performance mechanical guiding system and optical encoder ensures nanometer scale precision and constant static error, which can be calibrated by optical instruments. Mechanical guiding systems has its beneﬁts in case of control theory and its stability; unfortunately, on the other hand, there exists high inﬂuence of structure geometry and tribological effects such as friction and modal response. The aforementioned effect cannot be straightforwardly identiﬁed during the assembly process. Degradation of dynamic units can be detected only after certain operating time. Single degree of freedom systems are well investigated and the effect of degradation can be predicted, but there exists a gap in the analysis of nanometer scale multi degree of freedom dynamic systems; therefore, novel diagnostic tools need to be proposed. In this particular paper, dual axes dynamic system analysis will be presented. The main idea is to decouple standard stacked XY stage and analyse X and Y conﬁguration as two different conﬁgurations of the same object, while imitators of corresponding axes are absolutely solid and stationary. As storage and analysis of time domain data is not efﬁcient, main attention will be concentrated on frequency domain data, while, of course, statistical and graphical representation of dynamic response will be presented. Transfer function, dynamic response, spectral analysis of dynamic response, and modal analysis will be presented and discussed. Based on the collected data and its analysis, comparison of X and Y responses to different velocity excitation will be presented. Finally, conclusions and recommendations of novel diagnostic way will be presented.


Introduction
The precision of the manufacturing process, increasing efficiency and machining accuracy, determines the application of digitized precision computer control (CNC) in machine building [1,2]. It is impractical, non-economical and probably impossible to produce a perfect mechanical system. Air-bearing based systems can be a solution but, unfortunately, sharp drawbacks and disadvantages in particular designs occur; frictionless systems are more sensitive to changes of inertia, it is very hard to obtain stability of its dynamic response in a wide range of velocities, it is not compact, pneumatic feedthroughs are required; therefore, friction-based systems can be a very good solution for some of the application.
cally [29]. The characteristic of the vibrational response is that it is harmonic, but stochastic during the time. Vibrations amplify modal response in unexpected frequencies, influence the preload in the guiding system, dynamic response of the stage and finally its static error. Generally speaking, increased level of vibration and its non-linear behaviour, reduces overall positioning accuracy of the translation system. To reduce these effects, Kalman filters (KF), field oriented control (FOC) and sliding mode controller (SMC) are broadly used, while eddy current damping (ECD) can be additionally evaluated [28,30,31].
Due to vibrations, while the system is in motion and especially in start-stop operation, the dynamic error of the system increases in time. Due to vibration magnitude irregularity and frequency irregularity, it is hard to identify degradation of the robotic system by using standard control algorithms; therefore, novel control and diagnostic algorithms are required [32]. The traction harmonics of a linear motor, even at small oscillation amplitudes, cause oscillations that directly affect the mechanical system itself [16,33,34]. They have a cumulative effect that causes already large displacement fluctuations. In precision manufacturing, these substations in smooth motion are not taken into account [35], as the dynamic characteristics of the mechanical system is modeled and controlled by process controller. This is relevant given that the oscillations affect the entire mechanical system, and in particular the traction of the motor through the feedback system [36], so a detailed study of the electromechanical connection caused by the traction harmonic is necessary.
Separate studies have been performed to analyze mechanical oscillation, current circuit permeability [37], current upper limit, and control parameters in a linear motor supply system [15] but have only been defined by feedback errors without detailed analysis of traction harmonics and displacement oscillations. Such studies were performed by Yang et al. [16], where the dependences between air gap fluctuations and the harmonic error is displacement fluctuation were defined. This paper describes the modal and dynamic response of X and Y configuration linear motor stage with quasi-industrial linear guiding system and aims to find out the margins and actual response of the single degree of freedom stage excited in two different configurations: X (while fixed constrain of stage housing is organized to have full contact area with the surface) and Y (while fixed constrain of the stage housing is organized to have limited contact area with the base surface). PID with different configurations of notch filter, low pass filter and second order bi-quad filter is used to ensure system stability, without adjusting the system mechanical construction. To obtain the frequency response functions LMS is excited by current in the limited bandwidth and fixed amplitude (which might be considered as force excitation). The bandwidth is compared with the results obtained by the accelerometers. Spectral and dynamic response is investigated to analyse important X and Y configuration information. In addition, the Fourier spectral analysis of the amplitudes of the displacement dynamic errors is performed.
The obtained results of the investigation are represented by the derived transfer functions of the plant, controllers, the whole open loop LMS. The experimentally identified frequency response functions (FRF) are compared with FRF of the theoretically derived transfer functions of the LMS.

The Object of Investigation
The one degree of freedom linear motor stages with quasi-industrial guiding system, hereafter denoted as LMS.
LMS is assembled on heavy granite base (pos. 2 of Figure 1) along with the load (pos. 4 of Figure 1) and imitators (pos. 5 of Figure 1) to ensure real-world application configuration and resulting conditions. Position and velocity control loops of the platform are both closed by reflective optical linear encoder. The measuring head of the encoder is attached directly to the moving platform (pos. 1 of Figure 1) to ensure the precise displacement feedback. In LMS under consideration, the guiding system is quasi-industrial linear bearings with recirculating balls. Bandwidth and frequency response of the LMS is estimated by accelerometers (pos. 3 of Figure 1). Additionally, frequency response function is estimated by high performance servo controller and power amplifier by exciting moving platform by constant magnitude sine sweep (pos. 6 of Figure 1). The information was acquired by computer, MMI studio diagnostics software and spectrum analyser (pos. 7 of Figure 1). ration and resulting conditions. Position and velocity control loops of the platform are both closed by reflective optical linear encoder. The measuring head of the encoder is attached directly to the moving platform (pos. 1 of Figure 1) to ensure the precise displacement feedback. In LMS under consideration, the guiding system is quasi-industrial linear bearings with recirculating balls. Bandwidth and frequency response of the LMS is estimated by accelerometers (pos. 3 of Figure 1). Additionally, frequency response function is estimated by high performance servo controller and power amplifier by exciting moving platform by constant magnitude sine sweep (pos. 6 of Figure 1). The information was acquired by computer, MMI studio diagnostics software and spectrum analyser (pos.7 of Figure 1).

The Test Equipment
The test stand consists of: granite base by "JFA" (pos. 1 of Figure 2), LMS 8MTL165-300-Len1 by "Standa" (pos. 2 of Figure Figure 2). The investigation was done in "Standa" mechatronic laboratory premises. The general view of the test equipment is shown in Figure 2.

The Test Equipment
The test stand consists of: granite base by "JFA" (pos. 1 of Figure Figure 2). The investigation was done in "Standa" mechatronic laboratory premises. The general view of the test equipment is shown in Figure 2.
Motor windings are directly connected to current amplifier CMHP3A34N0N800WNA NNN, which is part of the motion control setup (pos. 5 of Figure 2). Position and velocity control loops of the platform both are closed by reflective type optical encoder. The encoder interface is analog Sin/Cos which is interpolated directly through CMHP3A34N0N800WNA NNN analog encoder interface. Analog signal is interpolated x4096 times by internal multiplier to ensure encoder resolution equal to 4.8828125 nm. Sin/Cos encoder errors caused by assembly procedures is prevented by adjusting of Lisajou graphs after final assembly of robotic system and adding work load. Motion controller (single board computer or master ACS SPii+, which is built into CMHP3A34N0N800WNANNN) is used to establish the communication with PC (pos. 4 of Figure 2) using TCP/IP interface. Further, by using PC (pos. 4 of Figure 2) and combination based on motion controller ACS SPii+ and current amplifier CMHP3A34N0N800WNANNN (pos. 5 of Figure 2) real time control and diagnostics of LMS (pos. 3 of Figure 2) will be investigated. To obtain the frequency response functions LMS is excited by current in the limited bandwidth and fixed amplitude (which might be considered as force excitation as force have direct dependence on current trough motor force constant K f ). The spectral response is compared with spectrum of dynamic response and results obtained by the accelerometers.
Mathematics 2021, 9, 951 5 of 26 Figure 2. The general view of the test equipment: 1 is granite base by "JFA"; 2 is LMS 8MTL165-300-Len1 by "Standa"; 3 is diagnostics system and motion controller ACS Spii+ by "ACS Motion Control Inc."; 4 is PC with control and diagnostic MMI Studio software by "ACS Motion Control Inc."; 5 is motion controller and amplifier module CMHP3A34N0N800WNANNN by "ACS Motion Control Inc."; 6 is accelerometer system by "Brüel & Kjaer".
Motor windings are directly connected to current amplifier CMHP3A34N0N800WNANNN, which is part of the motion control setup (pos. 5 of Figure 2). Position and velocity control loops of the platform both are closed by reflective type optical encoder. The encoder interface is analog Sin/Cos which is interpolated directly through CMHP3A34N0N800WNANNN analog encoder interface. Analog signal is interpolated x4096 times by internal multiplier to ensure encoder resolution equal to 4.8828125 nm. Sin/Cos encoder errors caused by assembly procedures is prevented by adjusting of Lisajou graphs after final assembly of robotic system and adding work load. Motion controller (single board computer or master ACS SPii+, which is built into CMHP3A34N0N800WNANNN) is used to establish the communication with PC (pos. 4 of Figure 2) using TCP/IP interface. Further, by using PC (pos. 4 of Figure 2) and combination based on motion controller ACS SPii+ and current amplifier CMHP3A34N0N800WNANNN (pos. 5 of Figure 2) real time control and diagnostics of LMS (pos. 3 of Figure 2) will be investigated. To obtain the frequency response functions LMS is excited by current in the limited bandwidth and fixed amplitude (which might be considered as force excitation as force have direct dependence on current trough motor force constant ). The spectral response is compared with spectrum of dynamic response and results obtained by the accelerometers.
In the next stage, spectral and dynamic response is investigated to analyse important X and Y configuration frequency domain information. Additionally, the Fourier spectral analysis of the amplitudes of the displacement dynamic errors is performed.

Mathematical Modelling of the Linear Motor Stage with Quasi-Industrial Guiding System
The block diagram of LMS and control system is shown in Figure 3. LMS with current amplifier (including current control loop which consists of "proportional-integral" controller) is hereafter called the plant-see Figure 3. After tuning of the current loop, the plant is considered as constant, while the influence to the modal and dynamic response In the next stage, spectral and dynamic response is investigated to analyse important X and Y configuration frequency domain information. Additionally, the Fourier spectral analysis of the amplitudes of the displacement dynamic errors is performed.

Mathematical Modelling of the Linear Motor Stage with Quasi-Industrial Guiding System
The block diagram of LMS and control system is shown in Figure 3. LMS with current amplifier (including current control loop which consists of "proportional-integral" controller) is hereafter called the plant-see Figure 3. After tuning of the current loop, the plant is considered as constant, while the influence to the modal and dynamic response of the different velocity excitation and LMS configuration is investigated. As it was already mentioned, in the present article the position and velocity loops together with the velocity loop shaping filters are called controller PID1, when configuration is X; controller PID2, when configuration is Y.
The position and velocity loops together comprise the control system. The position loop corresponds to "proportional" controller, velocity loop which consists of "proportionalintegral" controller and "loop shaping filters". Further, influence on dynamic error of certain configuration and different velocity excitation will be investigated.
The mathematical modelling of LMS is based on the following assumptions: (1) Constant friction forces appear in the moving parts of LMS; therefore, friction might be considered as minor component in mathematical model and identified frequency response function; (2) All parts of LMS are absolutely rigid. That is, only the rigid body dynamics is considered in the present investigation; of the different velocity excitation and LMS configuration is investigated. As it was already mentioned, in the present article the position and velocity loops together with the velocity loop shaping filters are called controller PID1, when configuration is X; controller PID2, when configuration is Y. The position and velocity loops together comprise the control system. The position loop corresponds to "proportional" controller, velocity loop which consists of "proportional-integral" controller and "loop shaping filters". Further, influence on dynamic error of certain configuration and different velocity excitation will be investigated.
The mathematical modelling of LMS is based on the following assumptions: (1) Constant friction forces appear in the moving parts of LMS; therefore, friction might be considered as minor component in mathematical model and identified frequency response function; (2) All parts of LMS are absolutely rigid. That is, only the rigid body dynamics is considered in the present investigation; (3) There is no backlash in the connections of the moving LMS parts; (4) The motion of the all moving LMS parts is straight-linear; (5) The displacements of the all moving LMS parts are identical; (6) The response of current amplifier can be considered as linear after tuning of PI controller in current control loop. Response of current amplifier to certain current excitation have rise time and overshoot which might be neglected. Dead time in PWM amplifier is ignored; (7) The moving mass is constant and amplifier bus voltage is constant; (8) Environment conditions are constant and well controlled; (9) Ironless BLDC motor have three symmetric phases. All phases of the motor are per-

Governing Equations
In the present section, the transfer functions of the open loop system consisting of the plant and the controller are derived analytically. To assess the accuracy of the analytically obtained transfer functions the frequency response functions obtained from the transfer functions were compared with the experimentally identified frequency response functions of the open loop systems of LMS.
The system of the governing differential equations of the LMS consists of three equations that relate the displacement of the LMS platform, the electromechanical parameters of the LMS devices, the current and voltage: From (1)- (3), it follows that: x (m/s) and ..
x (m/s −2 ) are the velocity and acceleration of the LMS platform or the first and the second order time derivatives of the displacement x of the LMS platform, respectively; m totx = 23.5 (kg), m toty = 25 (kg) is the total mass of the parts of the LMS platform that move in the straight-line motion in X and Y configuration; I is current (A); F motor is the motor generated force (N); K f = 12.9 (N/A) is motor force constant; L = 1.24 (mH) is inductance of the motor winding; R = 6.4 (Ohm) is resistance of the motor winding; U = 320.0 (VDC) is nominal bus voltage; K v = 15.0 ((Vs)/m) is motor Back EMF constant. In Equations (1)-(3), current I is time dependent given function, while quantities m tot , K f , K v , L, R and U are time independent constant parameters. The displacement x(t), the velocity dx(t)/dt . x(t) and acceleration d 2 x(t)/dt 2 a(t) of the LMS platform are unknown time functions.
As illustrated in Figure 3, the plant consists of motor, quasi-indsutrial guiding system, current amplifier and current controller, which is considered linear after the current loop is tuned. Current controller might be represented as transfer function given in Equation (5), which represents the PI controller: where, in Equation (5) The current PI controller is tuned empirically by applying periodic square wave into the windings of the motor. K p is considered as maximal after reaching the level of the reference current avoiding overshoot in whole range of current operation. The K i is considered as maximal after reaching minimal rise time of the current signal without destabilization of the actual current.

Theoretical Transfer Functions, Their Frequency Response Functions and Comparison with the Experimentally Identified Frequency Response Functions
The transfer functions obtained in this section will be written as ratios of polynomials and will be illustrated graphically and compared with the experimental results. One degree of freedom mass-spring mechanical system will not be analysed as a set of differential equation. The Laplace transform of differential equation of mechanical system will be considered as obvious. Additionally, ration between input force and output position (transfer function) will not be expressed to not overload the paper by standard analysis. Factorized style transfer function expressions will not be presented here in order to have more compact equations, to simplify the modelling by mathematical software tools. As mechatronic system is operated by controlling the current while current is directly proportional to the force generated by motor winding, it is mandatory to mention that transfer function is ration between force and acceleration. The structure of the controller is presented in Figure 3 while PID coefficients and filter bandwidths are represented in Table 1. The transfer function of the plants, denoted by H PLANT,X and H PLANT,Y that includes the motor current amplifier, see Equation (4), and the current controller see Equation (5), in the frequency domain, can be represented as the following polynomial ratio: 7.884·10 7 ·s + 6.622·10 10 6.126·10 −1 ·s 4 + 1.976·10 4 ·s 3 + 1.037·10 7 ·s 2 (6) H PLANT,Y (s) = 2.586·10 7 ·s + 2.722·10 8 5·580 −1 ·s 4 + 1.800·10 4 ·s 3 + 9.450·10 6 ·s 2 Before starting analysis of the dynamic errors of LMS, the experimental identification of PLANT was conducted by defining the empirical FRFsL PLANT,X andL PLANT,Y . Additionally, these functionsL PLANT,X andL PLANT,Y were compared with the FRF of the theoretical transfer functions H PLANT,X and H PLANT,Y denoted by L PLANT,X , L PLANT,Y . These frequency response functions are depicted in Figure 4.  Before starting analysis of the dynamic errors of LMS, the experimental identification of PLANT was conducted by defining the empirical FRFs ^, and ^, . Additionally, these functions ^, and ^, were compared with the FRF of the theoretical transfer functions  given in Equations (6) and (7).
From Figure 4, we can observe that the experimentally determined behaviour of the real PLANT has a set of the eigenfrequencies that can have an influence on the whole stability, dynamic displacement error, settling time and settling window of LMS. Usually, From Figure 4, we can observe that the experimentally determined behaviour of the real PLANT has a set of the eigenfrequencies that can have an influence on the whole stability, dynamic displacement error, settling time and settling window of LMS. Usually, single degree of freedom systems are considered as modular units for multi axis configuration; therefore, it is important to investigate both configuration while other element in system is stationary. As it was mentioned above, LMS has no transmission and the only friction source is quasi-industrial linear guiding system. Therefore, eigenfrequencies can be caused by the structure of LMS, its spatial configuration, friction, cable management, etc. The more detailed mechanical analysis and design revisions can be considered to damp the parasitic resonances of the system. However, in the present investigation, the influence of the spatial configuration to PLANT modal and dynamic response will be investigated without applying any mechanical or control revisions; therefore, brief recommendations for control theory scientists and engineers will be presented.
Further, high order PID controllers was added for both X (6th order) and Y (8th order) cases of PLANT was added in order to stabilize the PLANT and control its dynamic response see Equations (8) and (9): where coefficients a i , and b j , are following: a 1 = 3.482·10 10 , a 2 = 6.260·10 10 a 3 = 1.245·10 9 , where coefficients a i , and b j , are following: a 1 = 8.085·10 9 , a 2 = 1.032·10 10 a 3 = 2.631·10 8 , In both cases, a relatively similar control strategy was chosen see Figure 5. The main difference between regulation strategies is central frequency of second order NF: 280 Hz for X configuration and 165 Hz for Y configuration. It should be stated that X configuration has higher parasitic frequency than Y configuration LMS system. From Figure 4, it is obvious that X configuration of the stage has less aggressive eigenfrequencies in medium frequency range (50-500 Hz); therefore, higher bandwidth and higher dynamic range can be reached by tunning the proportion and integral regulator coefficients. vestigated without applying any mechanical or control revisions; therefore, mendations for control theory scientists and engineers will be presented.
where coefficients , and , are following: 1 = 8.085 ⋅ 10 9 , 2 = 1.032 2.631 ⋅ 10 8 , 4 = 2.750 ⋅ 10 5 , 5 = 5.108 ⋅ 10 2 , 6 = 2.441 ⋅ 10 1 , 7 = 1.954 1.1977 ⋅ 10 10 , 3 = 6.944 ⋅ 10 7 , 4 = 1.163 ⋅ 10 5 , 5 = 1.177 ⋅ 10 5 , 6 = 7 7 = 2.553 ⋅ 10 −5 , 8 = 4.303 ⋅ 10 −9 In both cases, a relatively similar control strategy was chosen see Figure difference between regulation strategies is central frequency of second order for X configuration and 165 Hz for Y configuration. It should be stated that X tion has higher parasitic frequency than Y configuration LMS system. From F obvious that X configuration of the stage has less aggressive eigenfrequencies frequency range (50-500 Hz); therefore, higher bandwidth and higher dynam be reached by tunning the proportion and integral regulator coefficients. Finally, to generalize, the previously discussed PLANT functions and controllers can be integrated, while transfer functions of open loop systems can be represented by Equations (10) and (11): where coefficients a i , and b j , are: where coefficients a i , and b j , are: a 1 = 8.085·10 9 , a 2 = 1.032·10 10 a 3 = 2.631·10 8 , To sum up all that was modelled and said above, a brief illustrative summary can be represented by Table 1. Table 1 represents bandwidth and dynamic capabilities of the PLANT itself, while indicating its stability margins and controller structure. 20% higher bandwidth and <10% lower load of X configuration stage corresponds to the wider capabilities of the X stage dynamics, while practically such a load and bandwidth differences can be neglected and can be considered in the tolerance.

Results and Discussion
In the present section, spectral and dynamic response of the LMS stages in X and Y configurations is investigated. From Section 3 before, it might be considered that the only difference between investigated LMS systems is configuration (spatial representation) of the stage and its mounting: X configuration is mounted to the granite via full surface area of housing, Y configuration system is mounted only by limited surface area flange imitating real XY configuration. LMS stages in X and Y configurations was examined by measuring: FRF via excitation of constant harmonic force signals, dynamic error of the displacement of the stage moving platform, see detail 1 in Figure 1. The FRF was identified while moving part was moving in order to decrease the influence of stick slip friction to the measurement result. The displacement was entailed by exciting LMS by different velocities (1,5,10,20,100,200, 400) mm/s, the mass of moving platform, the acceleration and the jerk were constant, while after the accusation of the dynamic error data and its statistical and spectral analysis, modal response of the moving platform was estimated and compared with FRF data, dynamic displacement error data transformed to frequency domain by applying Fourier transform. After converting dynamic data of dynamic process to frequency domain, 100 of the most powerful resonant frequencies were chosen and were additionally limited by 30 nm magnitude, which might be considered as non-measurable, therefore it will neglected in further analysis. After accusation of FRF and modal response, equation of resonance may be designed and actual position dynamics might be replicated by summation of harmonic signals. It is considered to be minor part; therefore, most of the analysis would be based on frequency domain data analysis, investigation of the resonance source and its modal shape.
Two configurations of the LMS were analysed: X configuration and Y configuration. Both setups were assembled from exactly same components in the same workshop. In each phase of the experiment the following sequence of activities were arranged: Spectral analysis of the dynamic response by applying Fourier transform of dynamic response data in order to identify most significant modes and its influence to dynamic response; -Simultaneously, during excitation by velocity set, modal response was collected by accelerometers to identify resonances and modal response in "feed" and "pitch" direction; -Spectral response of the three different methods was compared.
First of all, FRF response was analysed in order to identify most significant resonance information limited by bandwidth of 1000 Hz. From Figures 6 and 7 it can be identified that total number of "feed" direction resonances in X configuration is 9, comparing to the Y configuration system, which is 12. It can be also stated that in X configuration, most significant resonances are most likely located at the higher bandwidth 108-923 Hz, while in Y configuration, it is clearly obvious that there exists a narrow sequence of low frequency resonances in the 11.7-62.1 Hz range. Low bandwidth resonances are clearly caused by the Y configuration and its mounting technique. It is more complicated to damp low frequency resonances and those resonances are directly involved in the dynamics of the moving platform.     As we can see from Figures 8 and 9, the dependency of the displacement dynamic error ∆ e on the time is harmonic; therefore, dynamic response can be easily replicated by applying superposition of harmonic functions. The displacement error ∆ e is smaller and seems to be comparable in X and Y configurations at the low excitation velocities, i.e., when .
x ∈ {1,5,10,20,100} mm/s in comparison to the dynamic error ∆ e at the higher excitation velocities, when .
x ∈ {200,400} mm/s. The influence of LMS configurations and mounting technique on the dynamic errors ∆ e is bigger at the high velocity and not so visible at low velocity range.
In Table 2, there are summary estimations of the dynamic displacement error ∆ e of the LMS platform in different configurations X and Y at different velocities   As we can see from Figures 8 and 9, the dependency of the displacement dynamic error on the time is harmonic; therefore, dynamic response can be easily replicated by   As we can see from Figures 8 and 9, the dependency of the displacement dynamic error on the time is harmonic; therefore, dynamic response can be easily replicated by applying superposition of harmonic functions. The displacement error is smaller and  ard deviation and mean values of dynamic displacement error on earlier defined set of velocities.    x = 100 m/s, standard deviation s ∆e of the LMS in X configuration seem to be converging to near zero value 3.1000 µm, while standard deviation of the LMS in Y configuration seems to be diverging from zero with systematic increasement in its value and become peak 9.9000 µm at maximal . x = 400 m/s velocity; it can be stated that X configuration system is less predictable as standard deviation saturates at 100 mm/s velocity and does not significantly increases by adding additional power to the moving platform; while standard deviation in Y configuration system is more functionally predictable and is constantly increasing. Magnitude of error in Y configuration caused by excitation, diverge from stable near zero value and is >300% higher comparing to the same LMS in X configuration at maximal 400 mm/s velocity. • When the velocity . x = 200 m/s, mean value m ∆e of the LMS in X configuration is slowly converging to a positive 0.1244 µm value, while LMS in Y configuration is rapidly converging to a negative −0.3412 µm; the mean value is, significantly, >270% higher in Y configuration comparing to the same system in X configuration. In general, X configuration, which is fixed by using the whole housing surface, settles faster and dynamic error values increase is slower, while standard deviation of the dynamic error is not linear and saturates when velocity profile is systematically increased. In comparison with Y configuration behavior, system which has fixed constraints only through the flange of the limited area, dynamic error behavior increases faster, but is more predictable and regular. In connecting with these conclusions, further dynamic behavior in frequency domain will be investigated to understand the influence of X and Y configuration to dynamic error from frequency domain point of view.
Fourier transform of dynamic response data was made. After that, 100 of the harmonics with the highest amplitude was estimated, to analyse only part of the spectrum and initialize clusters of frequencies which might have highest influence to the dynamic response of the stages. In Figure 11, spectral representation of the dynamic response is presented. From the presented graphical information, it is only clear that most of the parasitic modal frequencies are located at the 10-250 Hz bandwidth in both configurations; the magnitude tendency of Y configuration stage is higher for <400% in peak response and nominally higher for >50%.
Further, more detailed spectral information was presented in Figures 12 and 13. In order to understand the frequency of each mode and its amplitude spectral response was separated for each investigated velocity velocities . x ∈ {1,5,10,20,100,200,400}. presented. From the presented graphical information, it is only clear that most of the parasitic modal frequencies are located at the 10-250 Hz bandwidth in both configurations; the magnitude tendency of Y configuration stage is higher for <400% in peak response and nominally higher for >50%. Further, more detailed spectral information was presented in Figures 12 and 13. In order to understand the frequency of each mode and its amplitude spectral response was separated for each investigated velocity velocities ˙∈ { , , , , , , }. x = 5000 mm/s. asitic modal frequencies are located at the 10-250 Hz bandwidth in both configurations; the magnitude tendency of Y configuration stage is higher for <400% in peak response and nominally higher for >50%. Further, more detailed spectral information was presented in Figures 12 and 13. In order to understand the frequency of each mode and its amplitude spectral response was separated for each investigated velocity velocities ˙∈ { , , , , , , }.   It is important to state that vibrations of magnitude less than 6 counts of encoder (single count encoder resolution is <5 nm) will be neglected in further analysis; therefore, bandwidth of >30 nm magnitude vibration is additional introduced in Table 3. Limited amplitude bandwidth will further help to sort the vibrations modes accrued by accelerometers.
From the investigation of dynamic response earlier, it is clear that, generally speaking, in X configuration, the level of vibration magnitude is significantly lower comparing to the Y configuration. It must be emphasized that spectral response in X configuration at velocity ˙∈ { , , , } mm/s is clearly influenced by higher frequency modes with bandwidth of up to 912 Hz, which is relatively easy to damp; on the other hand, at higher excitation velocity ˙∈ { , , } mm/s, significant influence is exerted by low frequency vibrations, which is obviously limited by bandwidth up to the 0-153 Hz, which is more complicated to damp; therefore, stability margins and overall bandwidth must be sacrificed. It must be added that if bandwidth cut at >30 nm magnitude of vibrations, X configuration spectral response bandwidth would be significantly narrowed and can be described by 0-39 Hz frequency range at all velocities. It is represented in Table 4.  It is important to state that vibrations of magnitude less than 6 counts of encoder (single count encoder resolution is <5 nm) will be neglected in further analysis; therefore, bandwidth of >30 nm magnitude vibration is additional introduced in Table 3. Limited amplitude bandwidth will further help to sort the vibrations modes accrued by accelerometers.
From the investigation of dynamic response earlier, it is clear that, generally speaking, in X configuration, the level of vibration magnitude is significantly lower comparing to the Y configuration. It must be emphasized that spectral response in X configuration at velocity . x ∈ {1,5,10,20} mm/s is clearly influenced by higher frequency modes with bandwidth of up to 912 Hz, which is relatively easy to damp; on the other hand, at higher excitation velocity .
x ∈ {100,200,400} mm/s, significant influence is exerted by low frequency vibrations, which is obviously limited by bandwidth up to the 0-153 Hz, which is more complicated to damp; therefore, stability margins and overall bandwidth must be sacrificed. It must be added that if bandwidth cut at >30 nm magnitude of vibrations, X configuration spectral response bandwidth would be significantly narrowed and can be described by 0-39 Hz frequency range at all velocities. It is represented in Table 4.
The situation in Y configuration is different. Spectral response at velocities .
x ∈ {200,400} mm/s bandwidth is limited by 0-154 Hz. If bandwidth cut at >30 nm magnitude of vibrations, Y configuration spectral response bandwidth as well as X configuration bandwidth will be significantly narrowed and will remain: in 0-39 Hz in X configuration and 0-153.7 Hz frequency at all velocities in Y configuration.
In conclusion, it is obvious that dynamics in X configuration system are significantly influenced by lower bandwidth frequency excitation, which is most likely coming from the granite base and vibration isolation system (mechanical dampers), while Y configuration spectral response bandwidth is wider, so harmonics which are coming from above 39 Hz can be considered as structural resonances caused by mounting, geometry, stage configuration, preload, and friction; therefore, further, it is important to analyse the range 39-154 Hz.  Additionally, modal analysis will be further introduced to estimate the modes which would correspond to the bandwidth accrued by FRF and dynamic analysis of the plant. The modal analysis will be represented by summary in Table 5 and Figures 14 and 15. (e)-mode 5 or "pitch" mode. Figure 14. Illustration of five Modal responses of X configuration stage; (a)-mode 1 or "feed" mode; (b)-mode 2 or "yaw" mode; (c)-mode 3 or "vertical" mode; (d)-mode 4 or "roll" mode; (e)-mode 5 or "pitch" mode.
In Table 5 is denoted first the five most significant modes in both X and Y configura tions. To illustrate structural deformations deflection direction is added to the addition columns and illustrated by Figure 14. The feed corresponds to the deflection to the dire tion of the motion; the yaw corresponds to the deflection around the vertical axis perpen dicular to the motion direction; the pitch corresponds to the deflection around horizont axis perpendicular to the motion direction; the roll corresponds to the deflection aroun horizontal axis parallel to the motion direction; the vertical corresponds to the deflectio to the direction perpendicular to motion direction in vertical plane.
It is clear that yaw mode 126.4 Hz and pitch 136.6 Hz modes in Y configuration ar significant and, most likely, it is important that they be prevented and damped by con troller, mechanical structure and interfaces. According to the literature, it can be state that modes between 100 and 140 Hz are most likely caused by preload in guiding system while preload is well stabilized in X configuration systems; unfortunately, because of th structural dynamics and micro vibrations, preload can variety depending on the motio profile and stage regime. Variation of preload can be a serious cause of degradation o guiding system; therefore, resonances in the range of 100-140 Hz in Y configuration mu be observed and controlled during the whole life cycle of stage in order to prevent degr dation of its kinematic performance.
Finally, modal response variation as a function of excitation velocity set will be pr sented in Figures 16 and 17 and analysed in Tables 6 and 7.  Figure 15. Illustration of five Modal responses of Y configuration stage; (a)-mode 1 or "vertical" mode; (b)-mode 2 or "yaw" mode; (c)-mode 3 or "feed" mode; (d)-mode 4 or "roll" mode; (e)-mode 5 or "pitch" mode. In Table 5 is denoted first the five most significant modes in both X and Y configurations. To illustrate structural deformations deflection direction is added to the additional columns and illustrated by Figure 14. The feed corresponds to the deflection to the direction of the motion; the yaw corresponds to the deflection around the vertical axis perpendicular to the motion direction; the pitch corresponds to the deflection around horizontal axis perpendicular to the motion direction; the roll corresponds to the deflection around horizontal axis parallel to the motion direction; the vertical corresponds to the deflection to the direction perpendicular to motion direction in vertical plane.
It is clear that yaw mode 126.4 Hz and pitch 136.6 Hz modes in Y configuration are significant and, most likely, it is important that they be prevented and damped by controller, mechanical structure and interfaces. According to the literature, it can be stated that modes between 100 and 140 Hz are most likely caused by preload in guiding system, while preload is well stabilized in X configuration systems; unfortunately, because of the structural dynamics and micro vibrations, preload can variety depending on the motion profile and stage regime. Variation of preload can be a serious cause of degradation of guiding system; therefore, resonances in the range of 100-140 Hz in Y configuration must be observed and controlled during the whole life cycle of stage in order to prevent degradation of its kinematic performance.
Finally, modal response variation as a function of excitation velocity set will be presented in Figures 16 and 17 and analysed in Tables 6 and 7 Table 6. Summarized spectral information of the modal response at the certain velocity excitation in bandwidth (0-400) Hz in X configuration.   Table 6. Summarized spectral information of the modal response at the certain velocity excitation in bandwidth (0-400) Hz in X configuration. It is obvious that in Y configuration there is 300% more significant harmonics than in X configuration-33 versus 11. On the other hand, it is important to take into account the bandwidth of open loop system, which is 10 Hz in Y configuration and 12.5 in X configuration. The margin frequency is the frequency which is strongly dependent on control system and its design, while frequencies above of the stability margin might be more interesting for analysis. From the summary presented in Figure 18, we could see that mounting configuration has significant influence on spectral response and, in X configuration, it is possible to estimate stability margin frequency (10 Hz), structural resonances (50 Hz) and guiding system/preload frequency (150 Hz) which have strong dependence on excitation velocity, while in Y configuration there is a number of frequencies which can be considered as a function of velocity, especially: 10 Hz, 90 Hz, 110 Hz, 120 Hz,  140 Hz, 150 Hz, 160 Hz, 260 Hz, 290 Hz, 310 Hz, 380 Hz. 390 Hz. Therefore, it must be stated that it is important to handle Y axes configuration with wider capabilities of bandwidth control and the "health" can be identified by observing modal response in the range of 10-150 Hz in X configuration, while in Y configuration wider bandwidth 10-390 Hz of observation must be considered. Additionally, analysing higher bandwidth in Figures 16 and 17, it can be observed that Y configuration system is less stable in frequency 2000 Hz, 2900 Hz, 3000 Hz.

Conclusions
From the analysis and investigation, the following conclusions can be presented:  After analysis of "feed" direction resonances accrued by FRF response, it can be concluded that X configuration can be characterized by wider bandwidth 108-923 Hz, while Y configuration bandwidth is limited by 11.7-62.1 Hz range in "feed" direction.  Further analysis of dynamic response by exciting LMS with different velocity showed that when the velocity ˙= { , , , , } m/s, LMS in X and Y configurations statistical values of dynamic response are comparable, while when velocity excitation was increased to ˙= { , , } m/s, dynamic response in X and Y configuration become less comparable, predictable and generally non-linear.  It is obvious that there exists exception for velocity ˙= m/s in the investigated LMS and when the velocity ˙= m/s, standard deviation of the LMS in X configuration is converging to zero and settles at 3.1000 μm magnitude, while standard deviation of the LMS in Y configuration diverging from zero with systematic increasement in its value and become peak 9.9000 μm at maximal ˙= m/s velocity. In X configuration LMS, standard deviation saturates at 100 mm/s velocity and does not significantly increases by adding additional power to the moving platform. Standard deviation in Y configuration LMS is linear and is constantly increas-  Table 6. Summarized spectral information of the modal response at the certain velocity excitation in bandwidth (0-400) Hz in X configuration.  Table 7. Summarized spectral information of the modal response at the certain velocity excitation in bandwidth (0-400) Hz in Y configuration.

Conclusions
From the analysis and investigation, the following conclusions can be presented: • After analysis of "feed" direction resonances accrued by FRF response, it can be concluded that X configuration can be characterized by wider bandwidth 108-923 Hz, while Y configuration bandwidth is limited by 11.7-62.1 Hz range in "feed" direction. x = 100 m/s, standard deviation s ∆e of the LMS in X configuration is converging to zero and settles at 3.1000 µm magnitude, while standard deviation of the LMS in Y configuration diverging from zero with systematic increasement in its value and become peak 9.9000 µm at maximal . x = 400 m/s velocity. In X configuration LMS, standard deviation s ∆e saturates at 100 mm/s velocity and does not significantly increases by adding additional power to the moving platform. Standard deviation in Y configuration LMS is linear and is constantly increasing.
Magnitude of error in Y configuration caused by excitation, diverge from stable near zero value and is >300% higher comparing to the same LMS in X configuration at maximal 400 mm/s velocity as in X configuration standard error deviation settles. • Talking about mean error value, it must be noticed that when the velocity reaches .
x = 200 m/s, mean value m ∆e of the LMS in X configuration is slowly converging to positive 0.1244 µm value, while LMS in Y configuration is rapidly converging to negative −0.3412 µm and the mean value is significantly >270% higher in Y configuration comparing to the same system in X configuration. • X configuration, which is fixed by using whole housing surface, settles faster and dynamic error values increase is slower, but standard deviation of the dynamic error is not linear and saturates when velocity in the moving profile is systematically increased and is non-linear. In comparison with Y configuration behaviour, which has fixed constraints only through flange of the limited area, dynamic error behaviour increases faster, but is more predictable and linear. • Spectral response in X configuration at velocity .
x ∈ {1,5,10,20} mm/s is clearly influenced by higher frequency modes with bandwidth of up to 912 Hz which is technically easier to damp; on the other hand, at higher excitation velocity .
x ∈ {100,200,400} mm/s significant influence to a dynamic response is done by low frequency vibrations which is obviously limited by bandwidth 0-153 Hz, which is technically more complex task to damp; therefore, stability margins and overall bandwidth must be sacrificed. It must be added that if bandwidth cut at >30 nm magnitude of vibrations, X configuration spectral response bandwidth would be significantly narrowed and can be described by 0-39 Hz frequency range at all velocities. In Y configuration spectral response at velocities . x ∈ {1,5,10,20,100} mm/s Is limited by 0-1000 Hz bandwidth, while at velocities .
x ∈ {200,400} mm/s bandwidth is limited by 0-154 Hz. If bandwidth cut at >30 nm magnitude of vibrations, Y configuration spectral response bandwidth as well as X configuration bandwidth will be significantly narrowed and will remain 0-153.7 Hz frequency at all velocities in Y configuration.

•
Dynamics in X configuration system is significantly influenced by lower bandwidth frequency excitation which is resulted by granite base and vibration isolation system (mechanical dampers), while Y configuration spectral response bandwidth is wider, so vibrations which are caused by 39-154 Hz frequencies can be considered as structural resonances caused by mounting, geometry, stage configuration and must be investigated further. • "Yaw" mode 126.4 Hz and "Pitch" 136.6 Hz modes in Y configuration are significant and must be prevented and damped by: controller or mechanical structure. Variation of preload can be serious cause of degradation of guiding system; therefore, resonances in the range of 100-140 Hz in Y configuration must be observed and controlled during the whole life cycle of stage in order to prevent degradation of its kinematic performance. Analysing Y configuration modal response, 300% more significant harmonics can be observed comparing it to X configuration: 33 compared to 11 in X configuration. • X configuration can be characterized by: control loop bandwidth (10 Hz), structural resonances (50 Hz) and guiding system/preload frequency (150 Hz) which have strong dependence on excitation velocity. Y configuration is represented by: 10 Hz, 90 Hz, 110 Hz, 120 Hz, 140 Hz, 150 Hz, 160 Hz, 260 Hz, 290 Hz, 310 Hz, 380 Hz. 390 Hz frequencies which have strong dependence on excitation velocities. "Health" can be identified by observing modal response in the range of 10-150 Hz in X configuration, while in Y configuration wider bandwidth 10-390 Hz of observation must be considered.

•
Higher bandwidth observations showed that Y configuration system is less stable in frequency 2000 Hz, 2900 Hz, 3000 Hz; therefore, narrowed bandwidth can be updated by individual frequency modes to observe LMS degradation.
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.