Next Article in Journal / Special Issue
Precise In-Process Strain Measurements for the Investigation of Surface Modification Mechanisms
Previous Article in Journal
Manufacture of Profiled Amabilis Fir Deckboards with Reduced Susceptibility to Surface Checking
Previous Article in Special Issue
Error Separation Method for Precision Measurement of the Run-Out of a Microdrill Bit by Using a Laser Scan Micrometer Measurement System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modal Analysis, Metrology, and Error Budgeting of a Precision Motion Stage

Department of Mechanical and Mechatronics Engineering, University of Waterloo, 200 University Avenue West, Waterloo, ON N2L3G1, Canada
*
Author to whom correspondence should be addressed.
J. Manuf. Mater. Process. 2018, 2(1), 8; https://doi.org/10.3390/jmmp2010008
Submission received: 18 December 2017 / Revised: 12 January 2018 / Accepted: 15 January 2018 / Published: 24 January 2018
(This article belongs to the Special Issue Precision Manufacturing)

Abstract

:
In this study, a precision motion stage, whose design utilizes a single shaft supported from the bottom by an air bearing and voice coil actuators in complementary double configuration, is evaluated for its dynamic properties, motion accuracy, and potential machining force response, through modal testing, laser interferometric metrology, and spectral analysis, respectively. Modal testing is carried out using two independent methods, which are both based on impact hammer testing. Results are compared with each other and with the predicted natural frequencies based on design calculations. Laser interferometry has been used with varying optics to measure the geometric errors of motion. Laser interferometry results are merged with measured servo errors, estimated thermal errors, and the predicted dynamic response to machining forces, to compile the error budget. Overall accuracy of the stage is calculated as peak-to-valley 5.7 μm with a 2.3 μm non-repeatable part. The accuracy measured is in line with design calculations which incorporated the accuracy grade of the encoder scale and the dimensional tolerances of structural components. The source of the non-repeatable errors remains mostly equivocal, as they fall in the range of random errors of measurement in laser interferometry like alterations of the laser wavelength due to air turbulence.

1. Introduction

Precision motion stages find extensive application in various industries for carrying out tasks related to manufacturing and inspection, as well as inside commercial products such as optical disk drives [1]. In this study, details of modal testing, laser interferometric metrology, and dynamic response to cutting forces, with subsequent error budgeting for a long-stroke linear nano-positioner is presented. A CAD drawing and photograph of the stage are presented in Figure 1. Some important features of design can be summarized as follows [2]:
  • The nano-positioner utilizes a single shaft engaged to air bushings, instead of the more commonly used double shaft arrangement. The roll resistance is provided by the air bearing at the bottom. This way, the self-aligning property of the air bushings, which are held in the housings using O-rings, is exploited to the greatest extent, making the manufacturing and assembly of the stage easy and low-cost.
  • Actuation is provided by two voice coil actuators (VCA) operating in moving magnet mode in complementary double configuration.
This paper performs a verification study of the nano-positioner, with emphasis on following a systematic approach. The methods can be applied in the metrology and dynamic analysis of similar precision motion systems that are used in precision and ultra-precision manufacturing.
The novel aspects of this work are summarized as follows:
  • A modified ‘peak-picking’ approach has been developed which directly utilizes the accelerance measurement, through its modified formulations for the estimation of modal parameters. This way, the bias at low frequency, observed in the receptance plots obtained from accelerance via double integration, is avoided. Searching the modal analysis literature [3,4,5], the authors were not able to find a peak picking method which directly works with accelerance. With this study, we hope to fill in this gap.
  • Modal testing, laser interferometric metrology, and error budgeting, although being established methods, have been applied to a novel precision motion stage design [2] for the first time. Outcome from these tests has allowed an in-depth evaluation of several design features regarding the overall accuracy and applicability of the motion stage in the micro-machining framework.
  • A hybrid error budget has been compiled which combines the commonly considered quasi-static geometric, thermal, and steady-state servo errors with the dynamic component due to cutting forces. The 3-axis harmonic deflections due to machining forces, at the workpiece level of the stage, could only be predicted using the spatial modal testing results.
In the literature, several works can be found involving the identification of vibratory dynamics of high precision motion stages using finite element analysis (FEA) [6,7,8,9,10,11]. While in [6,7,8,10,11] a precision positioning stage or some of its components were analyzed, in [9] different design alternatives for an ultraprecision micro-milling machine are evaluated. In this paper, instead of FEA, modal analysis is carried out by direct experimentation using impact testing with a hammer and two different accelerometers. Two independent methods for conducting and evaluating the tests are employed, and their results are compared. The main purpose of modal testing has been to verify the dynamic properties of the stage and also facilitate response prediction to multi-axis dynamic disturbances, such as cutting forces, which cannot be directly predicted with model identification that is based only on control system input/output data.
A number of works on the measurement of error motions of precision motion stages have also been published [12,13,14,15,16]. A specially designed laser interferometer for measurement in all six motion axes was employed in [12]. In [13], error motions of a linear stage were measured by comparing results from a laser interferometer, an autocollimator, and capacitance probes. In [14], a two-degree-of-freedom linear encoder capable of simultaneously measuring linear positioning and horizontal straightness errors was used. In [15], out-of-straightness and axis misalignment errors of the carriage slide of a drum roll lathe are measured using two capacitance probes. Reflective-type optical sensors for on-line real-time measurements were proposed in [16]. In this study, error motions of a long-stroke linear nano-positioner are measured using a laser interferometer. Geometric errors obtained this way are combined with servo errors, estimated thermal errors, and the predicted dynamic response due to machining forces to compile the error budget.
The paper is organized as follows. In Section 2 predicted vibratory dynamics, the experimental modal testing and analysis methodology, as well as comparative results from the two independent methods are shown. In Section 3, the laser interferometric measurement procedure is described and the observed error motions are presented. In Section 4, the predicted error budget using the data available at the design phase is shown. Then, the actual error budget is established using the errors due to servo, geometric, thermal, and machining force factors. In Section 5, conclusions are presented.

2. Vibratory Dynamics

In this section, first, the predicted vibratory modes using information available at the design phase are presented. Then, experimental modal testing results are shown.

2.1. Predicted Vibration Modes

A schematic diagram of the motion stage is presented in Figure 2. Catalogue values of stiffness for the air-bushings/bearings [17] and certain dimensions are presented in Table 1. Inertia values of the moving body are calculated using CAD program and presented in Table 2. Flat air bearings are usually not rated for their rotational stiffness. At the time of conducting the design, due to the lack of rotational stiffness data or models concerning flat air bearings, a simple model as shown in Figure 3 was assumed, for estimating the air bearing reaction moment due to the rotational motion. When the stage body is in its rotationally neutral position, the assumed distributed stiffness elements ( w k = k 2 y / L ) are preloaded.
The moment generated due to rotation can be expressed as:
M k = x = L / 2 x = L / 2 ( θ   x ) x w k d x = k 2 y   L 2 θ 12 = k 2 r θ
where k 2 y is the axial stiffness, and k 2 r is the estimated rotational stiffness. The air bearing length ( L ) is different along the X and Z axes ( L 2 x , L 2 z ), which results in different estimations for roll and pitch stiffness. It was observed later in the modal testing results presented in Section 2.6 that this approximation was not very accurate.
The natural frequency predictions for each significant vibration mode are presented in Table 3. For each motion axis, the natural frequency is found using the effective stiffness in that direction due to the bearings and the relevant mass or moment of inertia, using the analogy with a single-degree-of-freedom (SDOF) vibratory system [5]. The lowest predicted natural frequency is associated with the role mode at 538 Hz.

2.2. Summary of Modal Testing Methods

Modal testing of the linear nano-positioner is carried out using two independent methods as summarized in Table 4. The modified peak-picking approach which constitutes Method 1, has been developed as an alternative to the traditional peak-picking approach, with the primary difference being the usage of accelerance FRF directly, rather than converting to receptance first, which can introduce errors due to the measurement noise at low frequency. Details of how each method is employed are discussed in Section 2.3 and Section 2.4, for methods 1 and 2, respectively. Results are compared in Section 2.6.

2.3. Modal Analysis Method 1 (Modified Peak-Picking Approach)

2.3.1. Impact and Measurement Points

For method 1 (modified peak-picking), FRF measurements are taken in 3 different planes: XY, YZ, and XZ (Figure 4). For each measurement plane, accelerometer location and positive direction of acceleration measurement are indicated by A x z , A y z , and, A x y . In each measurement plane, impact locations and directions are indicated by F 1 to F 10 for A x z , F 1 to F 9 for A y z , and F 1 to F 4 for A x y .

2.3.2. Method of Analysis

In the first method, modal parameters, natural frequency ( ω r ) and damping ratio ( ζ r ), are identified using the developed ‘modified peak picking’ approach. ‘Peak-picking’, in the general sense, refers to the usage of graphical features of the real and imaginary parts of the FRF near the natural frequencies to estimate modal parameters [3,4,5]. Different resources may refer to slightly varying formulations of the ‘peak-picking’ method, although they share a similar basic idea. In this paper, a modified approach from traditional ‘peak-picking’ is taken in which the accelerance FRF is employed directly. The derivations of the formulas for this case are presented in the proceeding section side-by-side with formulas that have traditionally been used for receptance. Hence, by avoiding the numerical conversion from accelerance to receptance, the problem of double-integrating the low frequency noise can be largely circumvented through the use of the proceeding modified formulas.
Accelerance FRF between two coordinates (i and k) of a proportionally damped multi-degree of freedom (MDOF) system can be presented as a combination of vibratory modes as [3],
a i k ( ω ) = r = 1 N ω 2 ( 1 / k r )   ψ r i ψ r k 1 ( ω / ω r ) 2 + j 2 ζ r ( ω / ω r )
where, N is the number of modes, ψ r i and ψ r k are the i-th and k-th elements of the r-th mode shape vector, ω r is the natural frequency, ζ r is the damping ratio, and k r is the modal stiffness. In the case of a point FRF ( i = k ), ψ r i ψ r k = 1 can be set without loss of generality, which makes 1 / k r equivalent to the static compliance contribution of each mode. The receptance for the same FRF ( g i k ( ω ) ) can be obtained from the accelerance as,
g i k ( ω ) = a i k ( ω ) ω 2 .
Equation (2) can be separated into real and imaginary parts as,
a i k ( ω ) = r = 1 N ( G r i k ( ω ) + j H r i k ( ω ) ) G r i k ( ω ) = ω r 2 z r 2 ( 1 z r 2 ) ( 1 / k r ) ψ r i ψ r k ( 1 z r 2 ) 2 + ( 2 ζ r z r ) 2 H r i k ( ω ) = 2 ζ r ω r 2 z r 3   ( 1 / k r ) ψ r i ψ r k ( 1 z r 2 ) 2 + ( 2 ζ r z r ) 2
where z r = ω / ω r represents the normalized frequency. If the modes are assumed to be separated from each other (i.e., having sufficiently distant natural frequencies), the real and imaginary plots of the FRF near each eigenfrequency ω = ω r ( z r = 1 ) would resemble the characteristics of a single-degree-of-freedom (SDOF) system as shown in Figure 5 for both accelerance and receptance. Note that the extrema of the real and imaginary parts of the accelerance ( ω r , 1 , ω r , 2 , ω max ) around the natural frequency ( ω r ) attain close but different values than those for receptance ( ω r , 1 , ω r , 2 , ω max ), with the difference depending on the damping ratio ( ζ r ). However, the real part crossing of the abscissa is the same at ω = ω r for both cases.
Important characteristics of the real and imaginary plots of the FRF, and how they can be used to extract modal parameters (i.e., damping ratio, natural frequency, and modal participation factor) are summarized in the following, including the new steps proposed in our modified peak picking method:
  • Setting G r ω = 0 yields two positive roots as, ω r , 1 = ω r / 1 + 2 ζ r and ω r , 2 = ω r / 1 2 ζ r . Hence, in our approach, the frequency values coinciding with the minimum and maximum real components of accelerance around the natural frequency ( ω r ) are used to determine the damping ratio as,
    ζ r = ( ω r , 2 2 ω r , 1 2 ) ( ω r ) 2 4 ω r , 1 2 ω r , 2 2
In traditional peak picking, extrema of the real component of receptance are considered, which occur at the following frequencies, leading to the corresponding damping ratio expression,
ω r , 1 = ω r 1 2 ζ r ,   ω r , 2 = ω r 1 + 2 ζ r ,   ζ r = ( ω r , 2 ) 2 ( ω r , 1 ) 2 4 ω r 2
Figure 6a presents a comparison between the damping ratio values estimated from the real component of accelerance. When the proposed formulation in Equation (5) is used, the damping ratio is estimated correctly. However, if the traditional receptance formulation in Equation (6) is applied together with the approximation that ω r , 1 ω r , 1 and ω r , 2 ω r , 2 , then the damping ratio can only be estimated reasonably well if its actual value is 0.05–0.07 or less. In the case of the actual damping being higher, the use of Equation (6) yields a significant error when peak frequencies of real accelerance are substituted in place of their receptance counterparts. In the comparative modal testing results presented in Section 2.6, damping ratios as high as ζ r = 0.24 were encountered, with the vibration modes originating mainly from air bearing/bushing stiffness and damping properties. The receptance formulation in this case would have given incorrect estimates. Hence, the development of separate formulations for accelerance-based peak-picking was an obvious necessity in this study, and can be applied in other systems as well.
ii.
Setting H r ω = 0 yields only ω max = ω r 2 ζ r 4 ζ r 2 + 1 + 2 ζ r 2 1 as the positive root. For ζ r < < 1 , ω r ω max can be assumed. Hence, the imaginary peak/dip location is used to identify the natural frequency ( ω r ). For receptance, the imaginary peak location is obtained as,
ω max = ω r 2 3 ζ r 4 ζ r 2 + 1 2 3 ζ r 2 + 1 3
In Figure 6b, estimation of the natural frequency using accelerance versus receptance peak values is shown. It is observed that using the imaginary peak yields similar accuracy in both cases. Also, while the accelerance case overestimates, and deteriorates slightly faster for high values of ζ r , the receptance method underestimates the true ω r . Overall, both methods are suitable for peak picking to a certain extent for their type of measurement.
iii.
Limit ω / ω r yields G r ω r 2   ψ r i ψ r k k r and H r 0. Hence, real part of accelerance has residues from the lower frequency modes, and using the horizontal axis crossing of G r for natural frequency estimation would be inaccurate.
iv.
Limit ω / ω r 0 yields G r 0 and H r 0. Hence, in accelerance, higher frequency modes typically do not have an influence on their lower frequency counterparts. In the case of receptance, the situation is reversed in which the higher frequency modes affect the real part only, and lower frequency modes exert very little influence.
For mode shapes to be identified, either the accelerometer location can be fixed and force impacts at different locations can be applied (roving hammer), or the impact location can be fixed while the accelerometer is placed at different points for each measurement (roving accelerometer). Due to the reciprocity rule ( a i k = a k i ), results from the two cases should be equivalent in a linear system. Roving hammer measurements (for the same number of measurement points) can be carried out more quickly, as impacting at a point does not require any significant preparation. On the other hand, in the roving accelerometer case, more time is needed to properly mount the accelerometer at each measurement point, generally using wax. If one wants to determine mode shapes in three dimensions, which allows for a full three dimensional display of the vibratory motions, the response at every measurement point has to be measured in all three orthogonal axes. For the roving hammer case, this requires impacts in three orthogonal directions to be applied at each measurement point. This is very cumbersome; first, due to the difficulty of adjusting the orthogonal impact directions, second, due to the likelihood of some points being impossible to reach from all three directions. In such cases, it is much more advantageous to use a tri-axial accelerometer, which can output accelerations in all three axes at the same time, in roving accelerometer configuration. This way, both the problem of orienting measurement axes is solved, and the possibility of being obstructed by the measured structure is minimized, as the accelerometer is both smaller, and stays in place during measurement. In this paper, as the mode shapes are manually sketched in method 1, roving hammer configuration is used to obtain two dimensional mode shapes using hammer impacts from a single direction for each measurement point. On the other hand, taking advantage of the availability of three dimensional automated calculation and animation of mode shapes, in method 2, roving accelerometer configuration is used with a tri-axial accelerometer.
Denoting the accelerometer location as ‘o’, the accelerance FRF is given by a o k ( ω ) . As the imaginary peak/dip approximately occurs at ω r , the value of the peak/dip ( Q r o k ) can be expressed as:
H r o k ( ω r ) = Q r o k = ω r 2 ψ r o ψ r k 2 k r ζ r
Value of the imaginary peak/dip measured for a number of impact points, k = 1 , , m , can be related to the mode shape ( ψ r ) as,
[ Q r o 1 Q r o 2 Q r o m ] = ω r 2 ψ r o 2 k r ζ r [ ψ r 1 ψ r 2 ψ r m ]
As the mode shapes, which are essentially eigenvectors of the system dynamics, can be scaled by any constant factor, the imaginary peak/dip values can be directly used for visualizing the elements of the mode shape vector. An example case of how the modes are sketched is illustrated in Figure 7, for the YZ measurement plane. The values of Q r o k can be carried on the undeformed sketch of the structure using a graphical scaling factor. The deformed body is sketched using the displaced points, matching the displacements in their respective axes. For the actual analysis, additional points such as F 6 to F 9 are also considered.

2.4. Modal Analysis Method 2 (Software Package)

2.4.1. Impact and Measurement Points

The experimental setup for method 2 (software package) is presented in Figure 8a. In this method, the stage is impacted at the locations F1 and F2, as shown in Figure 8b. FRF’s are measured from the 3-axis accelerometer roved through A 1 A 8 , which totals to 24 FRF’s for each impact point. A related module of Test.Lab® was used to acquire and view the FRF’s.

2.4.2. Method of Analysis

In method 2, modal identification has been carried out using automated and typically complex algorithms that have been documented in the literature, and implemented inside the commercially available vibration analysis software package. Natural frequencies and damping ratios were identified using the ‘PolyMAX’ module within LMS Test.Lab®. The proprietary ‘PolyMAX’ algorithm carries out a similar operation to the commonly used least-squares time domain complex exponential method, in the frequency domain [18]. The resulting stabilization diagram is interpreted for natural frequencies and damping ratios. For constructing the mode shapes, these identified parameters are used in the least-squares frequency domain (LSFD) algorithm [19], which finds the best fit to the modal displacement vector based on the agreement between the measured and fitted FRF’s. The software package allows either complex or real mode shapes to be fit. In this paper, complex mode shapes are enabled to test the proportional damping assumption. Complexity of mode shapes is rated using ‘modal phase collinearity ( M P C )’ and ‘mean phase deviation ( M P D )’ [19]. M P C and M P D rate the complexity of the mode on a scale 0 to 100%, and 0°–90°, respectively. Having obtained a minimum M P C of 96.5%, and a maximum M P D of 12° in the set of identified mode shape vectors, the proportional damping assumption used in method 1 (Section 2.3) is observed to be justifiable. Identified mode shapes can be animated as a 3D video, and screenshots of the animated mode shapes are presented in Section 2.6.

2.5. Measurements from the Encoder

For the identification of axial modes in the X-direction or modes which have significant displacement components along the stage’s direction of sensitivity, measurements that are parallel to the encoder axis are needed. Such modes are critical for the positioning control stability, as they directly enter the control loop through the encoder measurement. In this regard, position readings from the encoder scale, evaluated by the DSpace® DS3002 encoder interface board, were fed to CutPRO®’s MalTF interface, as a position measurement, using the DS2102 digital to analog converter. The boards (DS3002, DS2102) ran at sampling frequency of 20 kHz. The same impact points in three planes mentioned for method 1 (Figure 4) were used, with the accelerometer replaced by the encoder. Receptances acquired this way did not yield any vibratory modes in the 0–2000 Hz range. Eventually, position control bandwidth in the X-axis could be increased up to 650 Hz without experiencing any interactions with vibratory modes, affirming these results. The bandwidth was mainly limited by the phase advance that can be contributed by the control scheme at the desired cross-over frequency, while ensuring that amplification of measurement noise through feedback did not display a significantly deteriorating effect.

2.6. Comparative Results and Discussion

Comparative modal testing results from method 1 (modified peak-picking) and method 2 (LMS Test.Lab®) are presented in Figure 9, along with natural frequency predictions made at the design phase (Section 2.1). It is observed that the identified natural frequencies for methods 1 and 2 are rather close, except for a slightly larger deviation for the first mode. Damping ratios are also observed to be close. On the other hand, large discrepancies between the experimentally identified and predicted natural frequencies can be noticed. This discrepancy is especially critical in the case of the first mode (roll). While methods 1 and 2 have measured 65 Hz and 79 Hz, respectively, the initial theoretical prediction was 538 Hz. The roll motion is only constrained by the flat air bearing at the bottom, and such bearings are usually not rated for rotational stiffness.
For the theoretical calculations, a simple model assuming distributed stiffness was used. Apparently, the actual rotational stiffness of the air bearing is much lower than approximated, which can be attributed to the distortion of the air cushion at the bearing interface and changing air flow conditions dependent on fly height and orientation. In order to make more accurate predictions, more detailed analyses and data on the rotational stiffness of the air bearing is necessary.
Identified modes 2 (horizontal) and 3 (vertical) also imply deviation of the actual vibratory dynamics from the predicted ones. The fact that natural frequencies measured for these two modes are close suggests symmetry in the actual system in the horizontal and vertical directions. On the other hand, such symmetry was not predicted due to the assumed contribution of the flat air bearing stiffness in the vertical direction. The effective normal stiffness of the air bearing appears to be lower than the catalogue value, which may be due to a higher gap in the final assembly. The air bearing still provides some stiffness, as evident from the slightly higher natural frequency identified for the vertical mode (mode 3). The reason for both modes 2 and 3 having lower natural frequency than expected can be due to a lower effective stiffness of the air bushings, likely due to the shaft being manufactured closer to the minimum diameter within the tolerance range.
Contrary to the other modes, the natural frequency identified for mode 4 (pitch) is higher than the theoretical prediction. When the main compliances causing a vibratory mode shape are due to the bearings, mode shapes assume rigid-body motion-like patterns, as assumed for the theoretical predictions. However, motions in each degree of freedom (linear and rotational axes) are not totally decoupled as it was assumed. This can be observed in modes 3 and 4, which have motions in both vertical and pitch directions. This can be the reason mode 4 attained a higher natural frequency than expected. Also, elastic deformations of the structural components couple with bearing compliances to alter vibratory dynamics. Mainly elastic modes were not identified in the 0–2000 Hz range, but a dominantly axial elastic mode was identified at 2765 Hz using method 1.
Overall, these modal measurements help assess the validity of many of the assumptions made during the design phase, and also highlight what additional knowledge is needed for making more realistic predictions. They also allow the response to a complex disturbance input to be evaluated in a fairly detailed manner, as performed in Section 4.2.3.
The procedure outlined in this section can be applied to other similar motion systems as follows:
  • For a quick and straightforward assessment of the most prominent vibration modes (starting from the first mode) the ‘peak-picking’ procedure can be applied, as it is shown to yield sufficient accuracy. This information can be incorporated in the determination of the control bandwidth.
  • For further investigation of the vibratory dynamics, typically with >5 modes, as well as for the representation of mode shapes in three dimensions, a software package similar to the one utilized in this paper can be used.
  • Results from the two methods can be combined in the assessment of design features, like the magnitude and geometry of the compliances of the stage and bearings, as exemplified in this study by the less than ideal roll resistance observed through modal testing.

3. Laser Interferometric Metrology

Laser interferometric metrology has been used to measure the geometric errors along the linear path of the motion stage using varying combinations of the laser source and the measurement optics (retroreflector, Wollaston prism, straightness retroreflector, angular retroreflector). The results from this section are used in compiling the error budget in Section 4.

3.1. Methodology of Measurements

Definitions for the machine tool coordinate system and the error motions are presented in Figure 10 [20]. Measurement target positions are chosen with some randomness in the spacing between each other as recommended in the standard [21], and are presented in Table 5. The stage dwells for 2 s at each position, and the average of measurement between 1.0–1.8 s is used. Each point is crossed 5 times in backward and forward directions.
During the tests, ambient pressure, temperature, and relative humidity were monitored. These values were later used in the Edlen equation [22] to obtain the laser wavelength. The locations of the measurement optics with respect to the center of gravity of the stage are presented in Figure 11. Abbe compensation is applied to the linear measurements using the angular measurements and the moment arms to bring them to the common point at the center of gravity. Linear errors in each axis can be expressed in terms of angular deviations and the 3D rotation matrix simplified for small angular deviations [23] as,
{ δ x δ y δ z } = [ 0 θ c θ b θ c 0 θ a θ b θ a 0 ] [ L x L y L z ]
where δ x , y , z are the linear errors, θ a , b , c are the angular deviations in each rotation axis, and L x , y , z are the moment arms.

3.2. Evaluation of the Results

The deviation at each point is given by d i j where i = 1, …, 18 is the index of measurement position and j = 1, …, 5 is the index of pass. Forward ( d ¯ ), backward ( d ¯ ), and unidirectional ( d ¯ ) mean positional deviations are expressed as:
d ¯ i = 1 n j = 1 n d i j d ¯ i = 1 n j = 1 n d i j d ¯ i = d ¯ i + d ¯ i 2
The estimator of standard uncertainty is found as,
s i = 1 n 1 j = 1 n ( d i j d ¯ i ) 2
where the backward direction ( ) version is found by replacing ( ) with ( ). For each error measurement, a standard final plot with deviations at each point at each pass ( d i j ), mean deviations ( d ¯ i , d ¯ i , d ¯ i ), and mean deviations with uncertainty range ( d ¯ ± 2 s , d ¯ ± 2 s ) are produced. The rest of the parameters related to accuracy and their definitions are presented in Table 6 [21].

3.3. Measurement Results

Linear and angular accuracies are presented in Table 7 and Table 8, respectively, as PV magnitudes. Detailed discussion of each error component is presented in the subsequent sections.

3.3.1. Linear Positioning Error (EXX)

Measurement results for EXX are presented in Figure 12. Abbe compensation is applied by adding the multiplication of the pitch error (ECX) by the 18.5 mm moment arm shown in Figure 11a, to the raw measurement result. The moment arm is measured from the encoder position, as the positioning error of the encoder scale w.r.t. the center of gravity is taken care of by the servo. The uncertainty of measurement due to the misalignment of the laser beam is given in the standard [21] as:
U = 0.3   ε 2 L
where U is the uncertainty in (μm), ε is the misalignment in (mm), and L is the measurement length in (m). As the laser encoder system used [24] provides an LED lamp for the indication of the level of alignment, ε < 1 mm can be assumed [21]. On the other hand, this still implies up to 15 μm possible error. Hence, the approximately 2 μm total slope removed from the Abbe compensated EXX (Figure 12b) can easily be the cosine error due to misalignment, justifying its removal. Alternatively, it can be due to the encoder scale not being mounted at an absolute right angle to the motion stage axis. In the latter case, it is an actual error rather than a measurement error and it should not be removed. However, it would still be a repeatable error that can be compensated by shifting the position reference.
After the removal of the slope, the bidirectional systematic error is given by E = 1.2 μm (Table 7). It is compatible with the expectations as the Heidenhain® LIP501 R encoder scale is graded for ±1 μm accuracy [25], which corresponds to PV 2 μm possible error. The repeatability is evaluated as R = 0.7 μm, which is poorer than expected for the typical optical encoder. Considering the approximately 300 mm dead path, the repeatability corresponds to about ±1 ppm deviation, which can be attributed to environmental disturbances on the laser measurement, such as the turbulence of the ambient air. Altogether, the accuracy is evaluated as A = 1.8 μm.

3.3.2. Straightness Error in Y (EYX)

Measurement results for EYX are presented in Figure 13. Abbe compensation has been applied to the raw results by subtracting the multiplication of the pitch error (ECX) by the 45.5 mm moment arm. As straightness measurements are prone to large errors due to the misalignment of the straightness retroreflector, the constant slope portion has been removed after Abbe compensation (Figure 13c). EYX results indicate a possible curvature on the main shaft, which might have been introduced during the machining of the flat surfaces, used for mounting the top and bottom plates. A least-squares best fit line has been used to shift the zero of the deviations axis. The repeatability of R = 0.5 μm is more likely due to laser measurement errors as in EXX. An overall accuracy of A = 2.5 μm has been measured.

3.3.3. Straightness error in Z (EZX)

Measurement results for EZX are presented in Figure 14. Abbe compensation has been applied by adding the multiplication of the yaw error (EBX) by the 45.5 mm moment arm to the raw result. As in EYX, the slope due to the misalignment of the optics is removed and the reference is shifted to the least-squares best fit line. A repeatability of R = 0.9 μm and overall accuracy of A = 1.6 μm are observed.

3.3.4. Yaw Error (EBX)

Measurement results for EBX are presented in Figure 15. A repeatability of R = 3.4 μm/m and an accuracy of A = 7.2 μm/m have been obtained. As the angular measurements are immune to misalignment errors, no corrections of the raw results have been made.

3.3.5. Pitch Error (ECX)

Measurement results for ECX are presented in Figure 16. The large angle measured can be due to the main shaft having a curvature, and resulting in a semicircular motion of the stage. This is also a potential explanation of the parabola-like EYX measured (Figure 13c).

4. Error Budget

The error budget is made up of components related to geometric, thermal, and servo errors, as well as the dynamic response due to machining forces. In this section, first, a predicted error budget is presented which uses nominal properties of the components. Then, measured/estimated actual errors are used for the actual error budget. Results from the two cases are compared.

4.1. Predicted Error Budget

In the predicted error budget, peak-to-valley (PV) error magnitudes are used likewise the actual one. The items included can be summarized as follows:
  • The PV error due to position sensor resolution is given by E r e s = 0.97 nm, derived from the 4096 times arctangent interpolation of the 4 μm measurement signal period.
  • The linear encoder scale is rated for ±1 μm grating error [25], hence the PV error due to encoder grating defects is given by E g r a t i n g = 2000 nm.
  • The main shaft of the motion stage acts as the guideway for the air bushings. It is manufactured to a cylindricity tolerance of 5 μm as shown in Figure 17. Assuming that roughly half of the possible errors due to the errors in the shaft dimension are cancelled by its realignment in the air gap, PV errors in Y and Z directions can be assumed to be E Y = 2.5 μm, and E Z = 2.5 μm, respectively.
The predicted error budget is presented in Table 9. The PV predicted error is 4.1 μm. The rationale behind using the sum of root-mean-square (RMS) values and the final prediction being the mean of the arithmetic and RMS sums has been explained in Section 4.2.4.

4.2. Actual Error Budget

4.2.1. Geometric Component

Summary of geometric accuracies w.r.t. the center of gravity are presented in Table 10. When it is assumed that compensation of the systematic error in all 5 axes is possible, i.e., if the actuator is used as part of a multi-axis micro machine-tool, the geometric errors reduce to just the repeatability value. In the error budget, both cases are considered separately to demonstrate the extent to which the system can be corrected with compensation.

4.2.2. Thermal Component

Thermal Disturbance Sensitivities

The linear nano-positioner body is made up of aluminum, steel, permanent magnet material, and encoder glass. A diagram showing key variables relevant to thermal disturbance sensitivities is presented in Figure 18. The point of interest is defined at the center of the upper surface of the top plate as shown.
The point of interest is sensitive to the following thermal disturbances:
  • In the case of thermal expansion along the X-axis, combined effect of the expansion of the encoder scale and the top plate needs to be considered. If the encoder scale is thought of as fixed at its center to the top plate, the deviation in X-positioning would be represented by,
    δ x = ( α a l α e n c ) L 2 Δ T
    where α e n c is the coefficient of thermal expansion of the glass encoder scale specified by the manufacturer [25], and α a l is the coefficient of thermal expansion of Aluminum 6061. On the other hand, the encoder scale is held by clamps which do not exert a significant pressure on the scale and the scale can be thought of as decoupled from the top plate. In that case, X-axis positioning error needs to be revised as,
    δ x = α a l L 2 Δ T
    which corresponds to a worse scenario due to α a l > ( α a l α e n c ) .
  • Along the Y-axis, thermal expansion of the stage would push the point of interest upwards by
    δ y = α a l L 1 Δ T 2
    where L 1 is the thickness of the stage, and Δ T is the temperature variation. If the stage body was only constrained by the air bearing at the bottom, the thermal expansion would be δ y = α a l L 1 Δ T . As the air bushings are expected to counteract this, L 1 / 2 is used to approximate the equilibrium position.
  • Thermal expansion along the Z axis does not affect the point of interest.
    The thermal sensitivity in X and Y axes can be defined as,
    γ x = δ x Δ T = α a l L 2 ,   γ y = δ y Δ T = α a l L 1 2
    The total linear thermal sensitivity of positioning can be expressed as,
    γ T = γ x 2 + γ y 2
Values of the thermal sensitivities and parameters used in calculations are summarized in Table 11.

Thermal Disturbances and the Resulting Thermal Error

At the specified peak current, VCA’s dissipate 1.19 W each, which can result in considerable thermal disturbances in a micro/nano-precision setting. The thermal stability of the nano-positioner is tested using a jerk limited cubic acceleration profile trajectory with the maximum acceleration, A = 5000 mm/s2 (which is slightly short of the maximum achievable acceleration of 6480 mm/s2), and the resulting high feedrate, F = 200 mm/s, constrained by the stroke length limitation. With a 0.05 s dwell period prescribed at the far end of the stroke, the period of back and forth motions is determined as 0.35 s. Thermocouples were mounted at the adjustment ring near the VCA core and the top plate on the main shaft. The stage was run for 8 h, for which the collected temperature readings are presented in Figure 19.
From the plot, it can be inferred that after the initial warming-up phase, temperature keeps rising steadily at a rate of approximately 0.1 °C/h. As the adjustment ring is closer to the VCA’s, it heats up in advance. The ambient temperature could possibly be kept steadier with better air-conditioning and ventilation equipment, which would also contribute to the cooling of the stage. However, the proximity of equipment such as the amplifiers, voltage supply, controller, and the PC to the experimental setup makes it a challenge. It should also be noted that the compressed air supply at the air-bushing/shaft interface contributes to the cooling of the stage due to the high heat transfer coefficient associated with the forced convection in the resulting annular duct. However, this effect does not amount to an observable isolation of the top plate from the heat generated by the VCA’s, due to the low flow rate of the air bushings. Assuming that the micro-milling operation for which the nano-positioner is intended can involve a calibration cycle every 2 h., thermal variation of the stage can be budgeted using Δ T = 0.2 °C. Hence, the resulting dimensional error is given by γ T × Δ T = 92 nm.

4.2.3. Machining Force Component

Machining forces contribute to motion errors by causing dynamic deflections (i.e., vibrations). The directions of anticipated machining force components acting on the payload were presented in Figure 1. The machining force input in the X-direction enters the control system as shown in the block diagram of Figure 20 [2]. In the block diagram, K a (A/V) is the current control response, K f (N/V) is the actuator force factor, 1 / m s 2 stands for the rigid body dynamics of the actuator body floating on the air bearings, and G ( s ) stands for the vibratory modes in the axial direction. The overall plant is given by,
G p = K a K f   ( 1 m s 2 + G ( s ) )
which was measured using frequency response, for controller design purposes. The dynamic compliance in the axial direction is given by the control system’s disturbance rejection function ( G F x ) as,
G F x = e F x = ( 1 m s 2 + G ( s ) ) 1 + C G p
where G ( s ) 0 can be set in the 0–2000 Hz range, as the measurements of Section 2.5 revealed no vibration modes.
In the case of F y and F z , dynamic compliances are given by the deformability of the structure as a function of the excitation frequency, which was captured in the FRF measurements in Section 2. Referring to Figure 4, FRF measured between F”1-Axy can be used as G F y , and the FRF measured between F1-Axz can be used as G F z . As the FRF’s are measured in terms of accelerance, conversion to receptance needs to be carried out for the position response. Instead of directly dividing the accelerance magnitudes by the square of the frequency (as suggested by the mathematical definition), receptances are synthesized from modal parameters which were fitted using peak-picking, in order to prevent distortions in the low-frequency region. Resulting G F x , G F y , G F z are presented in Figure 21a–c. The compliance in the Z-direction is about an order of magnitude higher than the others, due to the low effective roll resistance provided by the flat air bearing, as discussed in Section 2.6.
For the quantification of forces, a slotting case was considered with sinusoidal force components as follows:
F x = 1.0 + 1.0 × sin ( ω t )   ( N ) , F y = 0.2 + 0.2 × sin ( ω t )   ( N ) , F z = 1.0 + 1.0 × sin ( ω t )   ( N )
The resulting deflections can be obtained from δ x ,   y ,   z = G F x ,   y ,   z × F x ,   y ,   z . Phase differences in the timing of the forces can be neglected for a worst-case scenario to obtain the total deflection as,
δ = δ x 2 + δ y 2 + δ z 2
Depending on the nature of the cutting operation, the above force model may of course be replaced with a more suitable one as needed [26,27]. The resulting deflection ( δ ) as a function of the excitation frequency is presented in Figure 21d. The reported values comprise both the static response to the constant part and the sinusoidal part which changes as a function of the spindle frequency. In micro-milling, due to the compliance of the tools used, the machining forces are limited to 1–2 N in order to obtain the desired tolerances in the finished part. This is achieved by limiting the chip load using high spindle speeds [28]. Assuming a 2 teeth cutter, and a spindle speed of 50,000 rpm, a tooth passing frequency of ω = 1667 Hz is obtained, for which the total deflection is given by δ = 1.748 μm.

4.2.4. Compilation of the Error Budget

The error budget is presented in Table 12. Original errors are presented at the left, and non-repeatable errors after the assumed ideal geometric compensation are presented at the right. Servo error is obtained from trajectory following tests conducted with the stage. Thermal errors and machining force deflections are considered as non-repeatable. The error budget is compiled for the point of interest at the center of the top plate, therefore EXX is modified from the values obtained for the center of gravity (Table 10) using the Abbe moment arm and ECX.
When combining errors, the arithmetic sum gives the worst case, while the RMS sum provides a more realistic figure, assuming that the different error sources are probabilistic and not correlated. In practice, average of the two is used as a suitable estimate [29]. For the RMS sum to be applicable, variance of the presented error components has to be found, instead of PV magnitude. Assuming a uniform probability distribution centered at zero, RMS sum of the errors can be expressed as [30]:
E R M S = 1 2 3   [ i = 1 N ( E i ) 2 ] 1 / 2
where E i are the individual PV errors.
The error budget indicates a PV 5.7 μm deviation with the repeatable errors conserved, and 2.3 μm when 5 DOF compensation is assumed to be available. With these results, the nano-positioner can be concluded to be marginally suitable for micro-milling applications. Comparing measured values to the estimations from Section 4.1, linear and straightness errors are observed to be of similar magnitude. The non-repeatable parts observed in measurement were not anticipated in predictions. The source of non-repeatable errors is mostly equivocal, as they fall in the range of random errors of measurement in the applied laser interferometric method, for example due to the turbulence of ambient air. The systematic part of the errors in Y and Z axes are smaller than predicted, as the predictions employed the whole tolerance range of the shaft. In reality, within the tolerance range, the shaft may be manufactured to a better uniformity. In EXX, the systematic error measured is larger than estimated error due to the observed curvature of the shaft mostly evident in the EYX and ECX measurements. This affects the point of interest through the moment arm. Servo and thermal errors are observed to contribute a relatively small portion of the overall error. Dynamic deflections due to machining forces are anticipated as a major component of the error, mostly due to the compliance in the roll direction.
The error budget represents data that goes into the core of the performance specifications for a precision motion stage, in terms of accuracy. As such stages also constitute motion axes of micro-/nano-machine tools, they directly affect the achievable manufacturing accuracy. In this study, modal testing has been incorporated in the error budget to determine the dynamic deflections due to machining forces, in this regard. The new procedures proposed in this paper (i.e., modified peak picking method and merging of dynamic response predictions into an error budget) are applicable to a wide set of precision motion stages and machine tools.

5. Conclusions

The long-stroke linear nano-positioner design has been verified through modal testing, and laser interferometric measurements with subsequent error budgeting. The lowest resonance is observed to be a roll mode at 65 Hz, which was predicted to occur at 538 Hz before the stage was built. It is inferred that the roll/pitch resistance of the flat air bearing is much lower than the prediction, which highlights the need for additional data and modeling tools to aid designers of precision equipment. Nevertheless, this roll mode did not have any negative influence on the servo performance of the stage, which could achieve 450–650 Hz crossover frequency [2].
The error budget concludes PV 5.7 μm accuracy, which has a non-repeatable part of PV 2.3 μm. Systematic part of the errors could be reasonably predicted at the design phase, and they came more or less from the manufacturing tolerance and the assembly process. The semi-circular trend in the motion most likely originates from a curvature on the main shaft. There has also been a degree of uncertainty regarding the origin of the non-repeatable geometric errors. In the future, investment into better environmental control in the laboratory may help obviate some of the random sources of laser interferometric measurement error.
Overall, the modal analyses, laser interferometry, spectral analysis, and error budgeting has enabled the validation of many of the assumptions and predictions made during the design of the precision stage. The measurements and analyses have also clearly indicated the limitations of, and possible improvements to, this type of linear motion stage, such as enhancing the stiffness in the roll direction. The results further demonstrate the need for additional modelling capability and data regarding the multi-directional stiffness behavior of porous flat air bearings.

Acknowledgments

The authors would like to thank B. Tryggvason for his inputs into the development of the precision stage, and also University of Waterloo technicians R. Wagner and N. Griffett for their meticulous assistance in building and instrumenting it. This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) through Discovery grant RGPIN-03879 and Engage grant EGP #436910-12.

Author Contributions

This research has been carried out by Ahmet Okyay, as part of his PhD studies under the supervision of Kaan Erkorkmaz and Mir Behrad Khamesee. The supervisory roles of Kaan Erkorkmaz and Mir Behrad Khamesee included, but were not limited to, the structuring of the research problem, determination of the methods, advice on the relevant literature, direct participation in certain experiments, and the evaluation of results. The bulk of the research results have been obtained, and the write-up of the paper has been carried out by Ahmet Okyay with substantial editorial suggestions from the other authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tan, K.K.; Lee, T.H.; Huang, S. Precision Motion Control: Design and Implementation, 2nd ed.; Springer: Berlin, Germany, 2008. [Google Scholar]
  2. Okyay, A. Mechatronic Design, Dynamics, Controls, and Metrology of a Long-Stroke Linear Nano-Positioner. Ph.D. Thesis, University of Waterloo, Waterloo, ON, Canada, 2016. [Google Scholar]
  3. Ewins, D.J. Modal Testing: Theory and Practice; RSP: Taunton, UK, 1986. [Google Scholar]
  4. Schmitz, T.L.; Smith, K.S. Machining Dynamics: Frequency Response to Improved Productivity; Springer: Berlin, Germany, 2009. [Google Scholar]
  5. De Silva, C.W. Vibration: Fundamentals and Practice; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  6. Chang, S.H.; Tseng, C.K.; Chien, H.C. An ultra-precision XYθZ piezo-micropositioner, I. design and analysis. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1999, 46, 897–905. [Google Scholar] [CrossRef] [PubMed]
  7. Tenzer, P.E.; Mrad, R.B. A systematic procedure for the design of piezoelectric inchworm precision positioners. IEEE-ASME Trans. Mechatron. 2004, 9, 427–435. [Google Scholar] [CrossRef]
  8. Dong, W.; Tang, J.; ElDeeb, Y. Design of a linear-motion dual-stage actuation system for precision control. Smart Mater. Struct. 2009, 18, 095035. [Google Scholar] [CrossRef]
  9. Huo, D.; Cheng, K.; Wardle, F. A holistic integrated dynamic design and modelling approach applied to the development of ultraprecision micro-milling machines. Int. J. Mach. Tool Manuf. 2010, 50, 335–343. [Google Scholar] [CrossRef]
  10. Moon, J.H.; Pahk, H.J.; Lee, B.G. Design, modeling, and testing of a novel 6-dof micropositioning stage with low profile and low parasitic motion. Int. J. Adv. Manuf. Technol. 2011, 55, 163–176. [Google Scholar] [CrossRef]
  11. Li, Y.; Huang, J.; Tang, H. A compliant parallel xy micromotion stage with complete kinematic decoupling. IEEE Trans. Autom. Sci. Eng. 2012, 9, 538–553. [Google Scholar] [CrossRef]
  12. Liu, C.H.; Jywe, W.Y.; Hsu, C.C.; Hsu, T.H. Development of a laser-based high-precision six-degrees-of-freedom motion errors measuring system for linear stage. Rev. Sci. Instrum. 2005, 76, 055110. [Google Scholar] [CrossRef]
  13. Gao, W.; Arai, Y.; Shibuya, A.; Kiyono, S.; Park, C.H. Measurement of multi-degree-of-freedom error motions of a precision linear air-bearing stage. Precis. Eng. 2006, 30, 96–103. [Google Scholar] [CrossRef]
  14. Kimura, A.; Gao, W.; Lijiang, Z. Position and out-of-straightness measurement of a precision linear air-bearing stage by using a two-degree-of-freedom linear encoder. Meas. Sci. Technol. 2010, 21, 054005. [Google Scholar] [CrossRef]
  15. Lee, J.C.; Gao, W.; Shimizu, Y.; Hwang, J.; Oh, J.S.; Park, C.H. Precision measurement of carriage slide motion error of a drum roll lathe. Precis. Eng. 2012, 36, 244–251. [Google Scholar] [CrossRef]
  16. Lee, H.W.; Liu, C.H. High precision optical sensors for real-time on-line measurement of straightness and angular errors for smart manufacturing. Smart Sci. 2016, 4, 134–141. [Google Scholar] [CrossRef]
  17. New Way Air Bearings. Catalogue; New Way Air Bearings: Aston, PA, USA, 2017; Available online: http://www. newwayairbearings.com/catalog/components (accessed on 17 January 2018).
  18. Peeters, B.; Van der Auweraer, H.; Guillaume, P.; Leuridan, J. The polymax frequency-domain method: A new standard for modal parameter estimation? Shock Vib. 2004, 11, 395–409. [Google Scholar] [CrossRef]
  19. Heylen, W.; Lammens, S.; Sas, P. Modal Analysis Theory and Testing; Department of Mechanical Engineering, Katholieke Universiteit Leuven: Leuven, Belgium, 1995. [Google Scholar]
  20. International Organization for Standardization. ISO/DIS 230-1 Test Code for Machine Tools, Part 1: Geometric Accuracy of Machines Operating Under No-Load or Quasi-Static Conditions; ISO: Geneva, Switzerland, 2009. [Google Scholar]
  21. International Organization for Standardization. ISO 230-2 Test Code for Machine Tools, Part 2: Determination of Accuracy and Repeatability of Positioning Numerically Controlled Axes; ISO: Geneva, Switzerland, 2006. [Google Scholar]
  22. National Institute of Standards and Technology (NIST). Engineering Metrology Toolbox; NIST: Gaithersburg, MD, USA, 2017. Available online: http://emtoolbox.nist.gov/Wavelength/Edlen.asp (accessed on 17 January 2018).
  23. Kuipers, J.B. Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality; Princeton University Press: Princeton, NJ, USA, 1999. [Google Scholar]
  24. Renishaw. RLE Fiber Optic Laser Encoder; Renishaw: Hong Kong, China, 2015. [Google Scholar]
  25. Dr. Johannes Heidenhain GmbH. Exposed Linear Encoders; Dr. Johannes Heidenhain GmbH: Traunreut, Germany, 2000. [Google Scholar]
  26. Altintas, Y. Manufacturing Automation: Metal Cutting Mechanics, Machine Tool Vibrations, and CNC Design, 2nd ed.; Cambridge University Press: New York, NY, USA, 2012. [Google Scholar]
  27. Chae, J.; Park, S.S.; Freiheit, T. Investigation of micro-cutting operations. Int. J. Mach. Tool Manuf. 2006, 46, 313–332. [Google Scholar] [CrossRef]
  28. Cheng, K.; Huo, D. (Eds.) Micro-Cutting: Fundamentals and Application; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  29. Thompson, D.C. The design of an ultra-precision CNC measuring machine. CIRP Ann. Manuf. Technol. 1989, 38, 501–504. [Google Scholar] [CrossRef]
  30. Shen, Y.L. Comparison of combinatorial rules for machine error budgets. CIRP Ann. Manuf. Technol. 1993, 42, 619–622. [Google Scholar] [CrossRef]
Figure 1. CAD drawing and photograph of the long-stroke linear nano-positioner.
Figure 1. CAD drawing and photograph of the long-stroke linear nano-positioner.
Jmmp 02 00008 g001
Figure 2. Air-bushing/bearing arrangement of the motion stage.
Figure 2. Air-bushing/bearing arrangement of the motion stage.
Jmmp 02 00008 g002
Figure 3. Coarse approximation for air bearing rotational stiffness.
Figure 3. Coarse approximation for air bearing rotational stiffness.
Jmmp 02 00008 g003
Figure 4. Impact and measurement locations for method 1 (modified peak-picking).
Figure 4. Impact and measurement locations for method 1 (modified peak-picking).
Jmmp 02 00008 g004
Figure 5. Real and imaginary part plots of accelerance and receptance of a sample FRF.
Figure 5. Real and imaginary part plots of accelerance and receptance of a sample FRF.
Jmmp 02 00008 g005
Figure 6. Comparison of using different peak-picking formulations: (a) Damping ratio estimation; (b) Natural frequency estimation.
Figure 6. Comparison of using different peak-picking formulations: (a) Damping ratio estimation; (b) Natural frequency estimation.
Jmmp 02 00008 g006
Figure 7. Using the imaginary peak/dip values to determine the modal displacements.
Figure 7. Using the imaginary peak/dip values to determine the modal displacements.
Jmmp 02 00008 g007
Figure 8. Experimental procedure for method 2: (a) Photograph of the impact hammer and the tri-axis accelerometer at the A4 position; (b) Impact and measurement points.
Figure 8. Experimental procedure for method 2: (a) Photograph of the impact hammer and the tri-axis accelerometer at the A4 position; (b) Impact and measurement points.
Jmmp 02 00008 g008
Figure 9. Comparative modal testing results from the two methods and predictions at the design phase. Dimensions in mm.
Figure 9. Comparative modal testing results from the two methods and predictions at the design phase. Dimensions in mm.
Jmmp 02 00008 g009
Figure 10. Machine tool coordinate system and error motions (copied by Waterloo University with the permission of the Standards Council of Canada (SCC) on behalf of ISO) [20].
Figure 10. Machine tool coordinate system and error motions (copied by Waterloo University with the permission of the Standards Council of Canada (SCC) on behalf of ISO) [20].
Jmmp 02 00008 g010
Figure 11. Locations of the measurement optics with respect to the center of gravity: (a) Retroreflector; (b) Wollaston prism. Dimensions in mm.
Figure 11. Locations of the measurement optics with respect to the center of gravity: (a) Retroreflector; (b) Wollaston prism. Dimensions in mm.
Jmmp 02 00008 g011
Figure 12. EXX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Figure 12. EXX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Jmmp 02 00008 g012
Figure 13. EYX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Figure 13. EYX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Jmmp 02 00008 g013
Figure 14. EZX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Figure 14. EZX results: (a) Raw measurement; (b) Abbe compensated with straight line fit; (c) ISO presentation after slope removal.
Jmmp 02 00008 g014
Figure 15. EBX results in ISO presentation.
Figure 15. EBX results in ISO presentation.
Jmmp 02 00008 g015
Figure 16. ECX results in ISO presentation.
Figure 16. ECX results in ISO presentation.
Jmmp 02 00008 g016
Figure 17. Tolerances on the main shaft.
Figure 17. Tolerances on the main shaft.
Jmmp 02 00008 g017
Figure 18. Diagram of the stage with key variables relevant to thermal disturbance sensitivities. Dimensions in mm.
Figure 18. Diagram of the stage with key variables relevant to thermal disturbance sensitivities. Dimensions in mm.
Jmmp 02 00008 g018
Figure 19. Temperature readings over 8 hours of continuous operation.
Figure 19. Temperature readings over 8 hours of continuous operation.
Jmmp 02 00008 g019
Figure 20. Control block diagram with the machining force disturbance.
Figure 20. Control block diagram with the machining force disturbance.
Jmmp 02 00008 g020
Figure 21. Frequency domain functions: (ac) Transfer function between forces and deflections in the X-, Y-, Z-directions; (d) Total deflection as a function of frequency.
Figure 21. Frequency domain functions: (ac) Transfer function between forces and deflections in the X-, Y-, Z-directions; (d) Total deflection as a function of frequency.
Jmmp 02 00008 g021
Table 1. Air-bushing/bearing stiffness properties and dimensions.
Table 1. Air-bushing/bearing stiffness properties and dimensions.
PropertySymbolValue
Air bushing axial stiffness k 1 y z 23 N/μm
Air bushing rotational stiffness k 1 b c 2.8 Nm/mrad
Air bearing axial stiffness k 2 y 35 N/μm
Air bearing roll stiffness (estimated) k 2 a 4.7 Nm/mrad
Air bearing pitch stiffness (estimated) k 2 c 7.3 Nm/mrad
Air bearing length along the X-axis L 2 x 50 mm
Air bearing length along the Z-axis L 2 z 40 mm
Distance between the middle of the air bushings along the X-axis L x 180 mm
Table 2. Inertia properties of the moving body.
Table 2. Inertia properties of the moving body.
PropertySymbolValue
Mass m 1.318 kg
Moment of inertia in A (Roll) I a 409 kg mm2
Moment of inertia in B (Yaw) I b 17,439 kg mm2
Moment of inertia in C (Pitch) I c 17,454 kg mm2
Table 3. Prediction of the natural frequencies.
Table 3. Prediction of the natural frequencies.
DirectionExpressionNatural Frequency (Hz)
Y (Vertical) ω y = 2 k 1 y z + k 2 y m 1248
Z (Horizontal) ω z = 2 k 1 y z m 940
A (Roll) ω a = k 2 a I a 538
B (Yaw) ω b = 2 k 1 b c + 2 ( L x 2 ) 2 k 1 y z I b 741
C (Pitch) ω c = 2 k 1 b c + 2 ( L x 2 ) 2 k 1 y z + k 2 c I c 748
Table 4. Comparison of the two independent methods used in modal testing.
Table 4. Comparison of the two independent methods used in modal testing.
FeatureMethod 1 (Modified Peak-Picking)Method 2 (Software Package)
Frequency response function (FRF) acquisition systemCutPRO® MalTF module
by Manufacturing Automation Laboratory (MAL), Inc.
LMS Test.Lab®
by Siemens-PLM Software
Testing procedureRoving hammerRoving accelerometer
Accelerometer typeDytran® 3035AG
(1-channel)
PCB Electronics® 356A02
(3-channel)
Impact hammer typeDytran® 5800SLDytran® 5800SL
Identification of the natural frequency and damping ratioModified peak-picking methodPolyMAX [18]
Identification of the mode shape vectorsModified peak-picking methodLeast-Squares Frequency Domain (LSFD) [19]
Presentation of the mode shapesManual 2D drawingsAutomated 3D animations
Table 5. Target positions used in error measurements.
Table 5. Target positions used in error measurements.
i P i (mm)i P i (mm)i P i (mm)
11.50977.5481313.545
22.50888.5011414.509
33.50799.5491515.533
44.5101010.5151616.547
55.5481111.5261717.549
66.5011212.5431818.505
Table 6. Error motion parameters [21].
Table 6. Error motion parameters [21].
ParameterDefinitionFormula
B Reversal value B i = d ¯ i d ¯ i , B = max { B i }
B ¯ Mean reversal value B ¯ = 1 m i = 1 m B i
M Range mean bidirectional positional deviation M = max { d ¯ i } min { d ¯ i }
E Systematic positional deviation E = max { d ¯ i ;   d ¯ i } min { d ¯ i ;   d ¯ i }
R Repeatability of positioning R i = 4 s i , R i = 4 s i ,
R i = max { 2 s i + 2 s i + | B i | ;   R i ;   R i } , R = max { R i }
A Accuracy A = max { d ¯ i + 2 s i ;   d ¯ i + 2 s i }
min { d ¯ i 2 s i ;   d ¯ i 2 s i }
Table 7. Linear accuracies.
Table 7. Linear accuracies.
EXX (Linear)EYX (Vertical)EZX (Horizontal)
(μm) bi bi bi
B N/AN/A0.1N/AN/A0.2N/AN/A0.1
B ¯ N/AN/A0.0N/AN/A0.1N/AN/A0.0
M N/AN/A1.2N/AN/A2.1N/AN/A0.8
E 1.21.21.22.12.22.30.80.80.9
R 0.70.70.70.40.30.50.80.90.9
A 1.61.71.82.32.42.51.61.51.6
Table 8. Angular accuracies.
Table 8. Angular accuracies.
EBX (Yaw)ECX (Pitch)
(μm/m) bi bi
B N/AN/A0.7NANA2.0
B ¯ N/AN/A0.0NANA−1.0
M N/AN/A4.5NANA163.4
E 4.44.54.5163.7163.2163.8
R 3.23.43.46.64.56.6
A 6.67.27.2167.0167.2167.5
Table 9. Predicted error budget.
Table 9. Predicted error budget.
Error ComponentsPV Magnitude (nm)
Position sensor resolution ( E r e s )0.97
Position sensor grating error ( E g r a t i n g )2000
Y straightness ( E Y )2500
Z straightness ( E Z )2500
Arithmetic sum7001
RMS sum1173
Mean4087
Table 10. Summary of geometric accuracies.
Table 10. Summary of geometric accuracies.
ComponentAccuracy (A)Repeatability (R)Units
EXX1.80.7(µm)
EYX2.50.5(µm)
EZX1.60.9(µm)
EBX7.23.4(µm/m)
ECX167.56.6(µm/m)
Table 11. Thermal sensitivities and parameters used in calculations.
Table 11. Thermal sensitivities and parameters used in calculations.
QuantitySymbolValue
Thermal coefficient of expansion of Aluminum 6061 α a l 23.5 ppm/K
Thermal coefficient of expansion of the glass encoder scale α e n c 8 ppm/K
Thickness of the moving body L 1 33.7 mm
Distance between the center of the top plate and the encoder scale L 2 10 mm
Thermal sensitivity along the X-axis γ x 235 nm/K
Thermal sensitivity along the Y-axis γ y 396 nm/K
Total thermal sensitivity γ T 460 nm/K
Table 12. The error budget.
Table 12. The error budget.
PV Magnitude (nm)
Repeatable Errors ConservedRepeatable Errors SubtractedEstimated at Design Phase
Linear (EXX)39646752001
Straightness
 Vertical (EYX)25034862500
 Horizontal (EZX)16209362500
Angular(included in Linear)-
Servo3030-
Thermal9292-
Machining force17481748
Total Error
 Arithmetic Sum995739677001
 RMS Sum15186211173
Mean573822944087

Share and Cite

MDPI and ACS Style

Okyay, A.; Erkorkmaz, K.; Khamesee, M.B. Modal Analysis, Metrology, and Error Budgeting of a Precision Motion Stage. J. Manuf. Mater. Process. 2018, 2, 8. https://doi.org/10.3390/jmmp2010008

AMA Style

Okyay A, Erkorkmaz K, Khamesee MB. Modal Analysis, Metrology, and Error Budgeting of a Precision Motion Stage. Journal of Manufacturing and Materials Processing. 2018; 2(1):8. https://doi.org/10.3390/jmmp2010008

Chicago/Turabian Style

Okyay, Ahmet, Kaan Erkorkmaz, and Mir Behrad Khamesee. 2018. "Modal Analysis, Metrology, and Error Budgeting of a Precision Motion Stage" Journal of Manufacturing and Materials Processing 2, no. 1: 8. https://doi.org/10.3390/jmmp2010008

APA Style

Okyay, A., Erkorkmaz, K., & Khamesee, M. B. (2018). Modal Analysis, Metrology, and Error Budgeting of a Precision Motion Stage. Journal of Manufacturing and Materials Processing, 2(1), 8. https://doi.org/10.3390/jmmp2010008

Article Metrics

Back to TopTop