Abstract:
A method is disclosed whereby a laser-based spherical coordinate measurement system is dynamically calibrated. A mechanical oscillator, such as, but not limited to, a Foucault pendulum is used to generate periodic motions which can be fitted to Fourier series models. The residuals between the experimental measurements and the model can provide information which can be used to calibrate the instrument. The calibration information is used to augment the ASME B89.4.19-2006 standard to improve sensitivity to cyclic errors and include the servo systems.

Description:
FIELD OF INVENTION 
       [0001]    The invention is drawn to calibration of a laser-based spherical coordinate measurement system, more commonly known as a laser tracker, by using a mechanical oscillator to produce a rich set of well behaved measurement data. 
       BACKGROUND OF THE INVENTION 
       [0002]    Laser-Based Spherical Coordinate Measurement Systems, more commonly known in the art as Laser Trackers, were introduced in the mid 1980s. See U.S. Pat. No. 4,714,339 Three and Five Axis Laser Tracking Systems to Lau and Hocken; U.S. Pat. No. 5,764,360 Electro-Optical Measuring Device for Absolute Distances to Meier; and U.S. Pat. No. 7,352,446 Absolute Distance Meter that Measures a Moving Retroreflector to Bridges and Hoffer; all three of which are incorporated by reference herein. 
         [0003]    Earlier instruments were limited to laser interferometers, which integrated length measurements in the radial direction by fringe counting. Later instruments include absolute distance measurement, which allows the beam to switch between a plurality of targets, or for the beam to be broken. See Large-Scale Metrology An Update Estler, Edmundson, Peggs, and Parker, Annals of the CIRP, Vol. 51/2/2002, which is incorporated by reference herein, for a review of the technology. 
         [0004]    These instruments are calibrated under The American Society of Mechanical Engineers Standard ASME B89.4.19-2006, incorporated by reference herein. The Foreward to the Standard states;
       Performance evaluation of a laser tracker presents challenges different from those associated with conventional Cartesian CMMs. Because of the very large working volume, no full-scale three-dimensional calibrated artifacts exist, and the design of the laser beam steering system is such that individual parametric errors cannot, in general, be isolated and measured individually. For any coordinate measurement system, a fundamental requirement is a test of its ability to realize the SI unit of length, the meter. In a laser tracker, the length scale is often a laser interferometer and usually one does not have a significantly more accurate reference interferometer with which to perform such a test.   For these reasons, the performance evaluation tests in this Standard consist primarily of point-to-point length measurements using calibrated artifacts that can be realized in a number of ways. Measured lengths are compared with manufacturers Maximum Permissible Error (MPE) specifications in order to decide conformance. Realization of the SI meter can be evaluated in a number of ways, including calibration of the laser interferometer, measurement of a series of short calibrated reference lengths, or measurement of a series of long calibrated reference lengths. Procedures are also included for testing the absolute distance measurement (ADM) capability of laser trackers that include this option.       
 
       Sources of Error 
       [0007]    In order to calibrate an instrument, it is instructive to understand the principles by which the instrument performs its functions.  Electronic Distance Measurement , J. M. Rueger, Springer-Verlag, third edition, 1990, is a standard reference on Electronic Distance Measurement (EDM). Chapter 6 of  Electronic Distance Measurement , C. D. Burnside, Blackwell Scientific Publications, Cambridge, Mass., incorporated by reference herein, specifically addresses calibration.  Optical Tooling for Precise Manufacturing and Alignment , Philip Kissam, McGraw-Hill, 1962; and  The Surveying Handbook , Russell C. Brinker and Roy Minnick, Chapman &amp; Hall, second edition, 1995 provide useful information on surveying, and angle measurement in particular. 
         [0008]    Burnside identifies the three main sources of error in EDM as the: Scale Error, Zero Error, and Cyclic Error. The Scale Error is primarily due to drift in the modulation frequency. The Zero Error is primarily due to error in the zero point of the instrument, i.e., a constant offset in the distance. The cyclic error is primarily due to crosstalk between the modulating electronics and the detected optical signal intermediate frequency, IF, as well as multiple reflections. The Cyclic Error is typically the most subtle source of error and one which the ASME Standard is not particularly sensitive to detect. 
         [0009]    In a typical phase measurement instrument, a laser or LED is modulated at f 1 . The reflected beam is detected and mixed with a frequency f 2  offset from f 1  by the IF. For example, in U.S. Pat. No. 5,455,670 Optical Electronic Distance Measuring Apparatus with Movable Mirror, to Payne, Parker, and Bradley, incorporated by reference herein, f 1 =1500.000 MHz, f 2 =1500.001 MHz, which produces a 1 kHz IF, which is digitized by the electronics. Measurements reduce to measuring the phase difference between the IF and a reference phase of the generating electronics. 
         [0010]    If the beam is pointed into free space, with no reflected beam, there will still be a slight constant IF component detected. In a typical measurement, the detected signal will be a combination of the reflected beam in combination with the constant noise signal. If the two are exactly in phase or exactly 180 degrees out of phase, the error in a phase measurement will be zero. If the noise is 90 degrees out of phase, the noise will introduce a slight error in the detected phase. 
         [0011]    The net result is that a plot of measured distance vs true distance will exhibit a small sinusoidal residual that is periodic in distance. For example, in the &#39;670 case, there will be a small sinusoidal residual with a period of approximately 100 mm in distance. The amplitude of the error will be sensitive to the amplitude of the reflected signal, i.e., the error will be more significant for distant measurements due to the fact that the noise is a constant and the return signal is lower amplitude for distant targets. 
         [0012]    It will be recognized by those skilled in the art that the cyclic errors can most easily be extracted by measuring many points over a distance corresponding to a few cycles rather than widely separated points, where for the &#39;670 case one cycle corresponds to approximately 100 mm. 
         [0013]    EDMs can also experience errors due to multiple reflections. These errors produce errors at harmonics of the modulation frequency. For example, in the &#39;670 case, harmonics would be seen at 50 mm, 33.3 mm, etc. 
         [0014]    Sources of error in the angle measurements can be due to: the azimuth and elevation encoders, non orthogonality of the azimuth and elevation axes, bearings, servo tuning and control, latency in the timing of the distance and angle measurements for a moving target, hysteresis, etc. For example, if the rotor and stator are not mounted concentric, there will be a cyclic error in angle. If the servo is not properly tuned, there can be lag or oscillations about the commanded position. 
         [0015]    Calibrations meeting ASME B89.4.19-2006 specifications are typically performed by the National Institute of Standards and Technology (NIST), or the manufacturer under laboratory conditions. The requirements are beyond the means of most laser tracker users. A simple method for checking a laser tracker in the field, or providing enhanced calibration capabilities in the lab would be desirable. 
       BRIEF SUMMARY OF THE INVENTION 
       [0016]    A method is disclosed whereby a laser-based spherical coordinate measurement system is dynamically calibrated. A mechanical oscillator, such as, but not limited to, a Foucault pendulum is used to generate periodic motions which can be fitted to Fourier series models. The residuals between the experimental measurements and the model can provide information which can be used to calibrate the instrument. The calibration information is used to augment the ASME B89.4.19-2006 standard to improve sensitivity to cyclic errors and include the servo systems. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0017]      FIG. 1  shows a typical configuration of the invention. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0018]    While the invention is described in the preferred embodiment as a laser tracker, and a pendulum, it will be understood that it is not limited to a laser or a pendulum. For example LEDs, microwaves, or other sources of electromagnetic radiation, and other mechanical oscillators may be employed. 
         [0019]    In order to detect all of the sources of error, it would be ideal to test a laser tracker in a large volume of known coordinates, and under dynamic tracking conditions, i.e., include the servo system in the calibration. Due to the large volume, it is non trivial to generate such a rich set of known coordinates and trajectories. The ASME standard includes provisions for Long Reference Lengths for Ranging Tests in section 6.4.4. These tests expose cyclic errors, but require an additional laser interferometer, and a precision rail. Angle measurements require precision length artifacts. The standard is silent as to dynamics of the servo system. 
         [0020]    Poirier, and Yaghjian addressed a related problem in U.S. Pat. No. 4,704,614 Apparatus for Scanning and Measuring the Near-Field Radiation of an Antenna, incorporated by reference herein. They used a Foucault pendulum to scan the radiation pattern of a large antenna by using the rotation of the earth to produce a relative rotation over the field of the stationary antenna. 
         [0021]    Mechanical oscillators, such as, but not limited to, pendulums are well known in the art that oscillate at well behaved frequencies. For example, a simple pendulum oscillates, to a first approximation, with a period 
         [0000]    
       
         
           
             
               
                 
                   T 
                   = 
                   
                     2 
                      
                     
                       
                         π 
                          
                         
                           ( 
                           
                             L 
                             g 
                           
                           ) 
                         
                       
                       
                         1 
                         / 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where T is the period, L is the length, and g is the acceleration of gravity. 
         [0022]    For small oscillations, T is independent of the amplitude of the oscillation. It is well known that for larger amplitudes, higher order terms must be included, as taught in U.S. Pat. No. 4,738,032 Inertial Pendulum to Elmer, included by reference herein. Moreover, g depends on the latitude, elevation, earth composition, proximity to mountains, etc. However, mechanical oscillations are repeatable enough to serve as a reference for time measurement, and measurements of g, and as will be shown hereinbelow can be adapted for calibrating a laser tracker. 
         [0023]    The Foucault pendulum, as described in Vol. 5, Encyclopaedic Dictionary of Physics, incorporated by reference herein, typically is constructed by a spherically symmetric mass supported by a music wire pivoting on knife edge bearings. By carefully releasing the mass from an initial position, the pendulum will oscillate in a plane. As Foucault demonstrated, the plane of the pendulum will rotate, with respect to the earth at an angular velocity 
         [0000]      ωsin λ  (2) 
         [0000]    where ω is the angular velocity of the earth, and λ is the latitude. 
         [0024]    As described in the same reference, an impulse normal to the plane, or an imperfect release, will result in an elliptical path which will rotate at an angular speed roughly proportional to the area of the ellipse, i.e., not in a simple vertical plane. Under ideal conditions, the elliptical path produces a circular path. It will be recognized that in addition to oscillating in the polar coordinates, the system will oscillate slightly in the radial direction due to oscillations in the forces on the support wire, the bearing mechanism, and the natural frequency vibration of the music wire. 
         [0025]    Despite all the complexities of the motion of the Foucault pendulum, the motions are deterministic, and in combination with accurate time synchronization, can be used to generate a large-scale 3-D calibration mechanism for laser trackers which provides a rich set of data. For even a modest length L, the pendulum maps out a volume sufficient to map multiple cycles of the EDM. For a factory environment, lengths L of tens of meters are practical. 
         [0026]    For the purpose of this discussion, consider the following gedanken experiment generally depicted as  10  in  FIG. 1  Assume a Foucault pendulum  11  suspended by a wire  12  hanging from a precision bearing  13 , such as a compound knife edge, with a mass  14  with a center of mass  15 . Assume a virtual omnidirectional retroreflector  16 , as will be explained in more detail hereinbelow, is virtually at the center of mass  15 . Assume an absolute, or interferometer, laser tracker  20  mounted on a stable mount  21  follows the retroreflector  16 . Starting at the rest position of the pendulum, acquire distance r, azimuth angle α, elevation angle θ, and time t measurements in the tracker coordinate system, e.g., (r,α,θ,t). 
         [0027]    While tracking, translate the pendulum to a release fixture  30 , such as a string, attached to a stable mount  31  where the pendulum  11  is carefully released, by burning the string (as performed by Foucault in the original experiment), so as to minimize rotation about the L axis  12 . The laser tracker  20  will map a complex path incorporating all of the complexities described in the references cited hereinabove. 
         [0028]    It will be recognized by those skilled in the art that r(t), α(t), and θ(t) will be smooth functions that can be fit to periodic functions, such as a power series or harmonic sinusoidal functions, i.e., a Fourier series. This is relatively easy to do with commercially available software such as Mathematica™, available from Wolfram Research, Champaign, Ill. It will also be recognized that (r,α,θ,t) can be converted into Cartesian coordinates (x,y,z,t), which will also be smooth functions x(t), y(t), and z(t), which can also be fit to harmonic sinusoidal functions. The motions can also be fit advantageously to a generalized coordinate system (as would be used in solving by Lagrange&#39;s equations in a classical mechanics text), (L,β,γ,t), where L is the length from the support bearing to the retroreflector (at the center of mass), and β and γ are spherical polar coordinates, which can also be fit to smooth harmonic sinusoidal functions L(t), β(t), and γ(t). 
         [0029]    For a perfect laser tracker, and an ideal pendulum moving in a vertical place, L(t) will be a constant, and β(t), and γ(t) will be sinusoidal functions of period T where 
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                     = 
                     
                       2 
                        
                       
                         
                           π 
                            
                           
                             ( 
                             
                               L 
                               g 
                             
                             ) 
                           
                         
                         
                           1 
                           / 
                           2 
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   or 
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
             
               
                 
                   
                     L 
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   ≈ 
                   L 
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
             
               
                 
                   
                     β 
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   ≈ 
                   
                     
                       A 
                       1 
                     
                      
                     
                       sin 
                        
                       
                         ( 
                         
                           
                             2 
                              
                             π 
                              
                             
                                 
                             
                              
                             
                               t 
                               / 
                               T 
                             
                           
                           + 
                           
                             φ 
                             1 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
             
               
                 
                   
                     γ 
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   ≈ 
                   
                     
                       A 
                       2 
                     
                      
                     
                       sin 
                        
                       
                         ( 
                         
                           
                             2 
                              
                             π 
                              
                             
                                 
                             
                              
                             
                               t 
                               / 
                               T 
                             
                           
                           + 
                           
                             φ 
                             2 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where A 1  and A 2  are amplitudes and φ 1  and φ 2  are phase angles which depend on the initial conditions at release. In the special case 
         [0000]      A 1 =A 2   (7) 
         [0000]      |φ 2 −φ 1 |=π/2  (8) 
         [0000]    the pendulum will map a circle. In the more general case the pendulum will map an ellipse, and higher order harmonic terms will be necessary to produce a good fit. For larger amplitudes, T will depend on A 1  and A 2 . 
         [0030]    It is not necessary to know a priori the actual values for L, g, T, A 1 , A 2 , φ 1 , φ 2 , or the additional terms needed for a best fit, in order to calibrate some aspects of the laser tracker. The facts that L(t), β(t) and γ(t) are all periodic functions of the same period T, and are all smooth functions are sufficient to identify some errors. 
         [0031]    Symmetry can be used as a constraint in the fitting algorithm. For a well behaved Foucault pendulum, the first principles of the physics dictates that the motions should be symmetrical about the rest position, i.e., β=0 and γ=0. For example, if β(t) passes through the origin at t′ 
         [0000]      β( t ′+τ)=−β( t′ −τ)  (9) 
         [0000]    and if γ(t) passes through the origin at t″ 
         [0000]      γ( t ″+τ)=−γ( t″ -τ).  (10) 
         [0000]    The fact that L(t), β(t) and γ(t) are all periodic functions of the same period T can also be used as a constraint in the fitting algorithm. 
         [0032]    The residuals of the measured data to the best fit functions will be an indication of the instrument errors. Moreover, those skilled in the art will understand how the errors relate to the sources of the errors. 
         [0033]    Any anomalies due to the pendulum should be repeatable and invariant for different locations of the tracker. Since the period is fixed by the length, the maximum velocity can be adjusted by the release position, and used to expose tracker servo problems. Multiple trackers may be used to check against each other, or in cooperation to produce coordinates based on lengths alone (multilateration), or angles alone (triangulation), or combinations. Residuals between the various fits will expose errors. 
         [0034]    Most trackers can not track at the zenith angle. However, since the motions are well behaved, an absolute distance instrument can predict the location at a future time and pick up the target after rotating 180 degrees in azimuth. For example, an instrument located directly under the pendulum could easily detect zero point errors and cyclic encoder errors by plunging the instrument, as is typically performed to test a theodolite. 
         [0035]    In an actual reduction to practice, the problem of mechanically locating a retroreflector at the center of mass is non trival due to the Abbé error. However, Bridges, Brown, and Ackerson in U.S. Pat. No. 5,861,956 Retroreflector for use with Tooling Ball, resolve the problem for a single tracker; and Parker in U.S. Pat. No. 7,110,053 Multidirectional Retroreflectors, both of which are incorporated by reference herein, resolves the problem for multiple trackers or large motions which produce large tracking angles. It will be understood that in either architecture, care must be taken to dynamically balance the pendulum for the additional retroreflector(s). 
         [0036]    It will be noted that the spirit and scope of the invention is not limited to a Foucault pendulum. For example, the same principles can be applied to; plane pendulums, torsional pendulums, compound pendulums, a mass hanging on a spring, tuning forks, vibrating beams, etc. It will also be noted that the invention is not limited to a laser tracker. For example, the method could also be employed to calibrate a motorized total station surveying instrument.