UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Enhancing the robustness of ESPI measurements using digital image correlation Bingleman, Luke 2010

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata

Download

Media
24-ubc_2010_fall_bingleman_luke.pdf [ 8.46MB ]
Metadata
JSON: 24-1.0071157.json
JSON-LD: 24-1.0071157-ld.json
RDF/XML (Pretty): 24-1.0071157-rdf.xml
RDF/JSON: 24-1.0071157-rdf.json
Turtle: 24-1.0071157-turtle.txt
N-Triples: 24-1.0071157-rdf-ntriples.txt
Original Record: 24-1.0071157-source.json
Full Text
24-1.0071157-fulltext.txt
Citation
24-1.0071157.ris

Full Text

      ENHANCING THE ROBUSTNESS OF ESPI MEASUREMENTS USING DIGITAL IMAGE  CORRELATION     by  Luke Bingleman  B.E.Sc., The University of Western Ontario, 2008  A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF  THE REQUIREMENTS FOR THE DEGREE OF  MASTER OF APPLIED SCIENCE  in  THE FACULTY OF GRADUATE STUDIES  (Mechanical Engineering)             THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver)     August 2010  © Luke Bingleman, 2010 ii    Abstract Electronic Speckle Pattern Interferometry (ESPI) provides a sensitive technique for measuring  surface deformations.  The technique involves comparison of the speckle phase angles within  surface images measured before and after material deformation.  This phase angle comparison  requires that the speckle positions be consistent in all images.  A lateral shift between images  by just one pixel substantially degrades ESPI measurements, while a shift of two or more pixels  typically causes complete speckle decorrelation and compromises the measurement entirely.   To prevent such lateral motions, the specimen and the optical system must be rigidly fixed.   This requirement typically prevents use of the ESPI method in applications outside laboratories  or where it is necessary to remove the specimen from the optical setup between ESPI  measurements.  Here, Digital Image Correlation (DIC) is used to track speckle motion caused by  specimen displacement between ESPI measurements.  The measured images can then be  mathematically shifted to restore the original speckle locations, thereby recorrelating the ESPI  measurements.  Examples are presented where ESPI measurements are successfully made with  specimen shifts in excess of 60 pixels.  A novel ESPI measurement technique where the specimen  is removed in between ESPI measurements is also developed and validated.          iii    Table of Contents Abstract ............................................................................................................................................... ii  List of Tables ........................................................................................................................................ v  List of Figures ...................................................................................................................................... vi  Nomenclature ..................................................................................................................................... ix  Acknowledgements ............................................................................................................................ xi  Dedication ......................................................................................................................................... xii  1  Introduction ................................................................................................................................. 1  1.1  Electronic Speckle Pattern Interferometry ................................................................................... 1  1.2  Challenges of ESPI ......................................................................................................................... 3  1.3  Proposed Method:  Hybrid ESPI ‐ Digital Image Correlation ........................................................ 4  1.4  Residual Stresses and the Hole Drilling Method ........................................................................... 5  1.5  Literature Review .......................................................................................................................... 7  1.5.1  ESPI Robustness ........................................................................................................................ 7  1.5.2  Digital Image Correlation Using Laser Speckles ........................................................................ 9  1.6  Research Goals ............................................................................................................................ 10  1.7  Overview of Experiments ............................................................................................................ 11  2  Electronic Speckle Pattern Interferometry (ESPI) ........................................................................ 12  2.1  The Nature of Light and the Laser Speckle Phenomenon ........................................................... 12  2.2  Speckle Size ................................................................................................................................. 14  2.3  A Typical ESPI Optical Setup ........................................................................................................ 16  2.4  Phase Stepping ............................................................................................................................ 18  2.5  Fringe Patterns and Phase Unwrapping ...................................................................................... 22  2.6  Speckle Decorrelation and Its Effects ......................................................................................... 24  3  Digital Image Correlation & ESPI Fringe Correction ..................................................................... 27  3.1  Digital Image Correlation ............................................................................................................ 27  3.1.1  Digital Image Correlation Using Painted Speckles .................................................................. 31  3.1.2  Digital Image Correlation Using Laser Speckles ...................................................................... 32  3.2  ESPI Fringe Correction Procedure ............................................................................................... 33  4  Experimental Validation ............................................................................................................. 40  4.1  Hybrid ESPI – DIC Apparatus ....................................................................................................... 40  iv    4.2  ESPI System Validation ................................................................................................................ 42  4.3  DIC Measurements Using ESPI Images ....................................................................................... 45  4.4  Measuring Surface Motions Using EDIC & the Effect of Speckle Size ......................................... 50  4.5  Null Fringe Correction ................................................................................................................. 60  4.5.1  Phase Gradients – Cause and Correction ................................................................................ 62  4.6  Validation of ESPI Fringe Correction Using Hole Drilling ............................................................ 69  4.7  Specimen Removal and Replacement ......................................................................................... 74  5  Conclusions and Future Work ..................................................................................................... 80  5.1  Conclusions ................................................................................................................................. 80  5.2  Future Work ................................................................................................................................ 84  References ......................................................................................................................................... 86        v    List of Tables Table 4.1 ‐ Apparatus components ............................................................................................... 42  Table 4.2 ‐ Calibrated vs. measured stresses ................................................................................ 44  Table 4.3 ‐ Camera parameters and speckle sizes ........................................................................ 51  Table 4.4 ‐ EDIC measurement standard deviation for painted and laser speckles ..................... 56  Table 4.5 ‐ Range and deviation of fringe correction process ...................................................... 73  Table 4.6 ‐ Removal and replacement vs. baseline numerical comparison ................................. 77     vi    List of Figures Figure 1.1 ‐ A typical ESPI setup (Adapted from [7]) ...................................................................... 2  Figure 1.2 ‐ A laser speckle image ................................................................................................... 3  Figure 1.3 ‐ ESPI fringe data from hole drilling ............................................................................... 3  Figure 2.1 – The speckle phenomenon (Adapted from [13]) ....................................................... 13  Figure 2.2 – Specimen under coherent and incoherent illumination ........................................... 14  Figure 2.3 ‐ Speckle image optical setup (Adapted from [14]) ..................................................... 15  Figure 2.4 – Specimen image with various speckle sizes .............................................................. 15  Figure 2.5 ‐ A typical ESPI optical setup (Adapted from [15]) ...................................................... 17  Figure 2.6 ‐ ESPI setup including propagation vectors (Adapted from [15]) ................................ 22  Figure 2.7 ‐ Phase stepping schematic (Adapted from [18]) ........................................................ 19  Figure 2.8 – Fringe unwrapping .................................................................................................... 23  Figure 2.9 ‐ Speckle decorrelation ................................................................................................ 24  Figure 2.10 ‐ The effect of speckle decorrelation ......................................................................... 25  Figure 3.1 ‐ Digital image storage of a painted speckle image ..................................................... 27  Figure 3.2 ‐ The correlation coefficient (Adapted from [21]) ....................................................... 28  Figure 3.3 – Template matching using the correlation coefficient .............................................. 29  Figure 3.4 ‐ Determining template displacement ......................................................................... 30  Figure 3.5 ‐ The correspondence problem ................................................................................... 31  Figure 3.6 ‐ The ESPI fringe correction procedure ........................................................................ 34  Figure 3.7 ‐ DIC subset selection ................................................................................................... 35  Figure 3.8 ‐ NCC distribution and maximum ................................................................................. 35  vii    Figure 3.9 ‐ NCC maximum and surrounding pixels ...................................................................... 36  Figure 3.10 ‐ Information loss due to image shifting .................................................................... 38  Figure 4.1 ‐ Apparatus schematic ................................................................................................. 40  Figure 4.2 – ESPI beam splitter assembly ..................................................................................... 41  Figure 4.3 – Adjustable mirror ...................................................................................................... 41  Figure 4.4 ‐ Drill assembly ............................................................................................................. 41  Figure 4.5 ‐ Fringe result using calibrated specimen .................................................................... 43  Figure 4.6 ‐ Sections of a phase stepped set of images ................................................................ 45  Figure 4.7 ‐ Standard deviation of DIC displacement estimates vs. amount of image shift ........ 47  Figure 4.8 ‐ Average correlation coefficient R2 vs. amount of image shift................................... 48  Figure 4.9 ‐ Coordinate system used to describe specimen motion. ........................................... 50  Figure 4.10 ‐ X Direction EDIC Measurements.............................................................................. 52  Figure 4.11 ‐ Y Direction EDIC Measurements .............................................................................. 55  Figure 4.12 ‐ Measurement range vs. speckle size ....................................................................... 56  Figure 4.13 – Null fringe correction .............................................................................................. 61  Figure 4.14 ‐ Conical asymmetrical illumination of a flat specimen undergoing a shift d ........... 62  Figure 4.15 ‐ Path length change of beam 1 due to specimen shift ............................................. 62  Figure 4.16 ‐ Conical beam illumination ....................................................................................... 63  Figure 4.17 ‐ Conical beam phase gradient validation results ...................................................... 65  Figure 4.18 ‐ Phase gradients for various specimen shifts ........................................................... 66  Figure 4.19 – Phase gradient subtraction ..................................................................................... 68  Figure 4.20 ‐ Recovered stress solution vs. X direction specimen shift ....................................... 71  viii    Figure 4.21 ‐ Recovered stress solution vs. Y direction specimen shift ........................................ 71  Figure 4.22 ‐ Recovered stress solution vs. +Z direction specimen shift ...................................... 72  Figure 4.23 – Recovered stress solution vs. ‐Z direction specimen shift ...................................... 72  Figure 4.24 ‐ Removal & replacement baseline tests with the specimen fixed ........................... 75  Figure 4.25 ‐ Specimen fixture ...................................................................................................... 76  Figure 4.26 ‐ Removal and replacement residual stress results ................................................... 77  Figure 4.27 – Removal and replacement fringe results ................................................................ 78   ix    Nomenclature mA   Complex amplitude of light wave  ma   Amplitude of light wave  m   Phase of wave  I   Light intensity  S   Speckle size     Source wavelength  F   Focal length of lens  a   Aperture diameter  M   Magnification  od   Distance from object to lens  id   Distance from lens to imaging plane  f   f‐number of lens  1 2,m mA A   Complex amplitudes of reference and illumination beams  1 2,I I   Intensities of reference and illumination beams     Phase of interference     Phase change of interference pattern due to surface displacement  , 'I I   Intensity distribution before and after surface displacement  '   Phase of interference after surface displacement  nI   Intensity distribution of nth phase stepped image  A   Mean of phase stepped set  B   Modulation (amplitude) of phase stepped set  n   Phase step of nth phase stepped image (0°,90°,180°, 270°)  sK    Sensitivity vector, bisector of illumination and reference beams  sd    Surface displacement at specimen  1 2,k k     Propagation vectors of reference and illumination beams     Fringe pattern, expressed on range of [0,1]  ( , )c s t   Correlation coefficient  ( , )f x y   Potential match area from the target image  ( , )w x y   Template being matched from the base image  ( , )NCC s t   Normalized cross correlation coefficient   z   Normalized cross correlation coefficient distribution  x    , ,x y xy     In plane stress components  yS   Material yield strength  2R   Pearson product moment correlation coefficient (equivalent to NCC)  pd   Physical pixel size     Optical path length change  1 2,    Angles of incidence of reference and illumination beams (beam 1 and 2)  1 2,    Cone angles of reference and illumination beams  ,L R    Change in optical path length at right and left side of image  ,L R    Change in phase at right and left side of image     Phase change across image due to specimen shift  D   Image width at specimen  d   Specimen shift  L   Average beam path length (illumination and reference)  , ,P Q R   2D linear equation coefficients describing phase gradient      xi    Acknowledgements I’d like to thank everyone at UBC who has helped me over the course of this research – faculty,  staff and my fellow students.  Your kindness and support was much appreciated and you have  my sincere gratitude.  In particular, many thanks to the mechanical engineering shop staff who  provided much needed instruction and advice during the apparatus design and fabrication  stages of my research.   Special thanks to Dr. Schajer, who has helped me expand my engineering knowledge immensely  and exposed me to the many facets of research in an academic setting.  I am especially grateful  for his patience, compassion and insight; without which this research would never have  reached fruition.  Finally thanks to my parents who have supported me unquestionably throughout my  endeavours.  Without their encouragement and support this would not have been possible and  I am always grateful.  xii    Dedication   To My Father              1    1 Introduction 1.1 Electronic Speckle Pattern Interferometry Interferometry is a powerful optical tool that has many applications within engineering and  physics.  When two at least partially coherent light waves are superimposed, a spatially  modulated intensity pattern is observed in the region of superposition.  A common example of  this phenomenon is the coloured fringes that occur at the surface of an oil slick.  Electronic speckle pattern interferometry (ESPI) is an interferometric technique that utilizes a  highly coherent light source, generally a laser.  The laser light is split and recombined at an  object’s surface.  The resulting interference, or speckle, pattern can then be used to infer  information about the surface itself.  Most commonly in‐plane or out‐of‐plane displacements  and or vibrations are measured.   Figure 1.1 shows a typical ESPI setup.  In this case  displacements occurring at the object’s surface in the direction of the sensitivity vector can be  measured.   Some typical applications of ESPI techniques include; full field measurement of surface  deformations [1], material defect detection [2,3] and vibration mode shape measurement.  [4]   Although these are typical uses of ESPI, the range of applications is extensive.  More unique  applications of ESPI include the detection of buried landmines [5] and the condition monitoring  of historical mosaics [6] for example.  Within this thesis the particular ESPI application studied is  the measurement of residual stress.    2      Figure 1.1 ‐ A typical ESPI setup (Adapted from [7])    The speckle images recorded within ESPI measurements are inherently random in nature.  The  phase of one speckle relative to its neighbours’ is completely random.  The speckle images are  interference patterns created by laser light and modulated by the laser light and surface  roughness of the object under investigation.  An enlarged view of a typical speckle image is  shown in Figure 1.2.  Although the speckle pattern itself is random, small surface deformations  cause each speckle to change phase consistently.  This phase shift is then visible by subtraction,  and the randomness of the images is effectively removed.  Problems arise when rigid body   motions  changing one set o occurs, th and the m   1.2 C One of th environm and in pa rigid‐bod measure characte Figure cause the ob  the pixel lo f images to  e phase ch easuremen hallenge e major cha ental distur rticular area y motions c ment can ea ristic speckle  1.3 ‐ ESPI fr dr ject to shift cation of th the other.  W ange inform t is compro s of ES llenges inhe bances.  In   under inve an comprom sily be lost.  pattern cre inge data fr illing    significantl e speckles fr hen this  ation is lost mised.     PI rent to ESP order to obt stigation, m ise the mea   This is due  ated by the om hole  y,  om    I measurem ain accurate ust not shif surement e to the way s  surface fea investiga about the comparin after som changes  phase ac surface in of a fringe F ents is that   results fro t during the ntirely.  If s urface info tures of the tion is comp  surface de g surface sp e mechanic cause the sp cordingly wi formation c  pattern as  igure 1.2 ‐ A they are ver m ESPI data  measureme hifts do occu rmation is re  specimen u letely rando formations  eckle patte al change.  T eckle patte th the surfa an be extrac shown in Fig  laser speck y sensitive t , the specim nt.  Even sm r the  corded.  Th nder  m.  Informa is obtained  rns before a hese minut rn to change ce, so that t ted in the fo ure 1.3.  Ea le image 3  o  en,  all  e  tion  by  nd  e    he  rm  ch  4    light and dark fringe represents surface displacements of half the wavelength of the light source  used.  Problems can arise however when the surface, and therefore the speckle pattern, is shifted  significantly.  If a speckle physically moves from one pixel location to another from image to  image, the critical phase information is lost and the resulting fringe pattern is severely  damaged.  Often, due to the extremely small (micrometer) scale of ESPI measurements,  environmental disturbances such as air currents or vibrations are enough to induce rigid body  motions that result in speckle decorrelation and corrupt the ESPI measurement.  All of the preceding examples of ESPI applications are currently performed in a laboratory  setting, where the environment can be precisely controlled.  However, in an industrial or  outdoor environment, there exists disturbances that will cause rigid body motions during  measurements and this will corrupt the measurement data.  This problem seriously limits the  use of ESPI measurements in field applications.    1.3 Proposed Method:  Hybrid ESPI - Digital Image Correlation Digital Image Correlation DIC is a well established, non contact optical technique.  It involves storing digital images of an  object in various states, and then mathematically comparing them to obtain information about  object displacements or deformations.  Frequently, DIC methods are used to compare many  sub regions of images to obtain full field displacement or deformation data. [8] In order to  compare and match images accurately, the surface of the object under investigation must  5    contain a pattern or optical features.  DIC can be performed with a wide variety of surface  patterns, including grids and dots, however it is most effective when using a completely  random pattern.  Recently, [9‐11] work has been done involving utilizing painted and laser  speckle patterns as the surface pattern for DIC measurements, leading to the emergence of  digital speckle correlation (DSC) measurements.     The method used within this thesis is a novel hybrid ESPI‐digital image correlation (EDIC)  technique.  An ESPI measurement is performed normally, and the specimen planar rigid body  motions that occur during measurement are determined from the existing ESPI speckle images  using DIC.  Any rotations that may have occurred are not addressed using this procedure.  Using  this DIC data, the ESPI images can be mathematically shifted to compensate for any rigid body  motions that may have occurred during measurement.  From these new corrected images,  accurate displacement data can be obtained.      1.4 Residual Stresses and the Hole Drilling Method Residual stresses are self‐equilibrating stresses that are locked into a material and may exist  without any external loading or thermal gradients.  The majority of manufacturing operations,  including turning, grinding, heat‐treating, surface hardening and welding, introduce residual  stresses into the workpiece.  Residual stresses can be quite large and can significantly influence  the behaviour of a material, especially with respect to fatigue life and dimensional stability.   Additionally residual stresses can impair material strength, and in extreme cases cause major  structural failure.  Conversely, residual stresses can have beneficial effects, for example, fatigue  life improvement due to shot peening.    6    The hole drilling method for measuring residual stresses is a very popular, well‐established  method.  One of its main advantages is that it inflicts minimal damage to the specimen.  A small  hole, typically a few mm in diameter, is drilled into the specimen to a depth approximately  equal to its diameter.  The surface deformations around the hole that result due to residual  stress relaxation are then measured.   Traditionally strain gauges are used to measure these  surface displacements.  From these surface displacements, the stress that existed within the  removed material can then be determined via an inverse calculation.    With the advent of low cost, high resolution CCD cameras, ESPI has recently been applied to  measure surface deformations within the hole drilling method. [12] Using an optical method  such as ESPI presents many advantages when compared with traditional strain gauge methods.   The main advantages include the relatively high speed at which measurements can be  performed.  Since all the surface information is contained within images captured via the  camera, there is no need to attach anything to the surface of the specimen itself.  The per‐ measurement cost using an ESPI setup is also very low.  The initial cost is significantly higher  when comparing ESPI to traditional strain gauge methods; however in an industrial setting in  which many measurements need to be performed, the additional cost is easily mitigated.    The challenges of ESPI persist when performing hole drilling measurements; if the specimen  moves even slightly before and after the hole is drilled, the entire measurement is  compromised.  This necessitates elaborate, extremely rigid fixtures to hold the specimen  accurately, as well as careful control of the measurement environment.  As well, since the  specimen must be so rigidly fixed, the mechanism to drill the hole must be integrated into the  7    ESPI system.  The action of drilling the hole may itself cause the specimen to shift,  compromising the measurement.  More practically, the chips from the drilling can also damage  the specimen’s surface and interfere with the ESPI measurement.     If it were possible to relax the constraints on specimen shifting before and after hole drilling, it  would allow ESPI to be performed in less isolated environments.  It would also be possible to  perform much more rapid ESPI measurements by eliminating the need for a very rigid fixture to  hold the specimen.  By removing the specimen to drill the hole and then replacing it, the drilling  procedure could be separated from the ESPI system and performed in an area where the  vibrations, cutting forces and chips from the process would not disturb the measurement.   Currently due to the inherent sensitivity of the ESPI method this is almost impossible.  It  therefore would be very valuable to develop a method to mitigate the negative effects of these  rigid body motions so that accurate ESPI measurements can be performed even under non‐ ideal conditions, and so that specimen removal and replacement would be possible.   1.5 Literature Review 1.5.1 ESPI Robustness The sensitivity to outside disturbances and delicate nature of ESPI equipment and  measurements are well known.  Currently the predominant method of circumventing these  limitations is to perform measurements in a very controlled, isolated environment.  Some work  however has been done with the goal of increasing ESPI robustness to provide a more direct  solution to the delicate nature of the measurement.    8    Ritter et al. [13] first investigated using a small Michelson interferometer in conjunction with a  conventional out of plane ESPI system to remove the effect of out of plane vibrations between  the ESPI system and the specimen.  This was accomplished by using the information from the  Michelson interferometer as an additional signal fed into the phase stepper of the larger  system.  By synchronizing the phase stepping mechanism with the environmental vibrations,  detrimental effects of these vibrations could be reduced for environmental vibrations with  frequencies up to 100 Hz.   Chau [14] proposed an alternate method of correcting for environmental vibrations within  digital shearing speckle interferometry (DSSI), a measurement technique similar to ESPI.  This  process continually monitored a DIC subset relatively unaffected by the measured  displacements within the images acquired.  By continuously computing the correlation  coefficient, an algorithm was developed that only acquired images at moments when the  subsets happened to be well correlated.  Using this result, DSSI fringes were able to be obtained  for vibrations with an amplitude of up to 1 mm and frequencies of 25 Hz.  Findeis et al. [15] performed experiments using ESPI to measure thermal stresses in the  presence of mechanical environmental vibrations induced by a compressor.  By increasing the  shutter speed of the camera significantly and utilizing a variety of mechanical and opto‐ electrical shutters, the effect of the environmental vibrations was mitigated significantly and  reasonable measurements were possible for the environmental disturbances tested.   Mechanical shutters were constructed in the form of discs with various slots to mechanically  simulate pulsed laser illumination.  By synchronizing the frequency of the rotating slots with the  9    exposure time of the camera, the effects of the environmental vibrations were able to be  somewhat mitigated.  For the best performance however this type of solution was found to be  unwieldy and too large for more portable systems.  Optical shutters were also investigated in  the form of a liquid crystal modulator (LCM) and lithium niobate crystal.  The LCM proved too  slow to be effective, however the lithium niobate crystal was able to accommodate the much  shorter exposure times required to diminish the effect of the environmental vibrations.     Reu and Hansche [16] recently suggested that speckle correlation may be recovered for ESPI  measurements with displacements on the same order as the speckle size at the specimen.   Similar to the method developed in this thesis, by utilizing DIC data in conjunction with ESPI  images, they showed speckle correlation could be reasonably recovered. Unfortunately details  of the results and methods were not discussed and the idea appears not to have been pursued  any further.  1.5.2 Digital Image Correlation Using Laser Speckles The use of laser speckles within DIC measurements has experienced rapid growth in recent  years, mainly due to the increased popularity and growth of digital image correlation methods  in general.  Initial work was performed by Peters, Ranson, Chu and Sutton [17‐19]  who first  examined using DIC methods in combination with laser speckles to measure full field  displacements and rigid body motions for specimens with optically rough surfaces.    Takai and Asakura [20]  examined the effectiveness of using laser speckle based DIC  measurements to measure full field displacements of a beam in bending and a specimen  undergoing thermal expansion.  They determined that the upper limit of measurable  10    displacement was controlled by the DIC subset size used for matching, and the lower limit of  displacement was somewhat controlled by the speckle size and pixel resolution of the imaging  system.  Chen and Chiang [21] as well as Sjödahl and Benckert [22] examined the optimum laser speckle  size for displacement, and work by both authors suggested speckle size selection according to  the Nyquist, or Shannon sampling theorem produced reasonably accurate results.  Chen and  Chiang examined in plane displacements of 200µm of an aluminum plate and performed  various measurements using a constant speckle size and varying the camera used in order to  vary the sampling frequency.  They determined using a camera with a physical pixel size of  approximately half (0.528) the speckle size produced the best results, where the correlation  coefficient peak was the best defined.  Another important conclusion from this work involved  identifying the spatial dependence of laser speckles; the pattern will not translate exactly with  the surface of a specimen. Taking a slightly different approach, Sjödahl and Benckert performed  a similar investigation.  They examined the effectiveness of DIC measurements on computer  generated speckle images with speckles of various sizes and using various subset sizes.  They  determined that using speckles of approximately twice the pixel dimensions of the camera  produced good results. Also, by keeping the ratio of subset to speckle size above 10, the DIC  measurement was successful about 80% of the time.  1.6 Research Goals Currently ESPI is largely limited to laboratory use within very controlled environments.  In many  applications however, measurements must be performed in non‐ideal environments.  These  11    environments can introduce rigid body motions during measurements that completely  compromise the measurement itself.   Although previous work has been done to minimize the  effect of environmental displacements, specimen shifting still remains a significant problem for  ESPI measurements.  Therefore, the main goal of this research is to develop a method to  counteract the effect of rigid body motions on ESPI measurements, with the ultimate goal of  helping to move ESPI from a fragile laboratory technique to a more robust field measurement  method.    In the case of ESPI applied to hole‐drilling, this new method would provide better resistance to  environment‐induced rigid‐body motions.  By counteracting specimen shifting, specimen  removal and replacement would also be possible in simple fixtures.  This would allow drilling  and its negative effects to be separated from the ESPI system as well as increasing  measurement speed and accuracy.  1.7 Overview of Experiments First, the performance and accuracy of DIC using ESPI speckle images (EDIC) will be evaluated  using the experimental apparatus.  The effectiveness of using DIC with both laser and painted  speckles will also be investigated.  The performance of the EDIC method for residual stress hole  drilling measurements will then be examined by performing measurements involving various  rigid body motions and disturbances.  Finally the possibility of the removal and replacement of  the specimen after drilling will be investigated.    12    2 Electronic Speckle Pattern Interferometry (ESPI) 2.1 The Nature of Light and the Laser Speckle Phenomenon Light transmission can be modeled as the propagation of a harmonic plane wave. Its amplitude  and phase can be mathematically represented by a complex phasor:      expሺ ሻm m mA a i   (2.1)  where:   = Amplitude of the wave      = Phase of wave  CCD cameras and similar sensors measure light intensity, which can be expressed as:    2mI A   (2.2)  A random ‘speckle’ pattern is created when illuminating a surface with a coherent source, i.e., a  source where all parts of the light have fixed relative phase.  The pattern forms due to the  constructive and destructive interference of light waves reflecting from the surface features.  A  typical imaging setup with a detector array and lens is shown in Figure 2.1.  If light waves  interfere at the imaging plane constructively, then a bright speckle will be formed, if they  interfere destructively then a dark speckle will be formed.    ma m                      The cohe relations construct coherent dark and specimen speckle p the surfa Since the recorded pattern it the surfa Fig rence of the hip between ive or destr  source and  light speckl  illuminated attern is no ce itself.     phase or am  as values o self can be  ce moves, t ure 2.1 – Th  source is im  the rays of uctive inter  the speckle es caused by  using conv t visible, an plitude of  f intensity, w considered  he speckles  e speckle p portant be  light shown ference. Fig  pattern is c  the cohere entional wh d the only fe the light can hich can be as a set of s move with i henomeno cause witho  in Figure 2 ure 2.2 (a) s learly visible nt interfere ite light.  Si atures in th not be mea  expressed  urface mark t.  The speck n (Adapted  ut it there i .1, and thus hows a spec .  It is chara nce.  Figure nce the sour e image are sured direct as shown in ers virtually les also cha from [26])  s no fixed p  no fixed po imen illumin cterized by   2.2 (b) show ce is not co  the physica ly, the spec  equation (2  attached to nge intensit hase  sitions of  ated with a the random s the same herent, a  l markings o kle pattern  .2).  The sp  the surface y with moti 13        n  is  eckle  .  As  on  14    (a)  (b) due to their spatial dependence; as the surface moves the incident light will interfere  differently causing the speckles to change phase from light to dark.  This is the governing  principle behind several speckle pattern techniques including ESPI.     2.2 Speckle Size The size of the speckles within an image can also be controlled by varying the optical setup  used. Considering the optical setup shown in Figure 2.3, the size of the speckles contained  within an image can be estimated as: [27]    (1 ) (1 )F MS f M a        (2.3)  where:   = Speckle size      = Source wavelength  = focal length of the lens  = aperture size  = magnification, defined as    = f number of lens defined as    Figure 2.2 – (a) Specimen illuminated with coherent light       (b) Specimen illuminated with white light  S F i o d d a M f F a  (a   For an op size of an on the ty advantag speckle s ) Fig tical system  image can  pe of measu eous ‐ the r ize.  Figure 2 Fig ure 2.3 ‐ Spe  with a set  be controlle rement, co esolution of .4 shows tw ure 2.4 ‐ Spe ckle image source and w d by varying ntrolling the  several typ o images ta ckle image (b)  optical setu ith compo  the f numb  size of the  es of speckl ken of the s  taken with  p (Adapted nents fixed a er of the im speckles wit e measurem ame specim (a) f = 8 (b)   from [27])  t set distan aging syste hin an imag ents is cont en with on f = 32    ces, the spe m.  Dependi e can be  rolled by th ly the f num 15  ckle  ng  e  ber  16    of the lens and the camera exposure time varied to let in equivalent amounts of light.   Comparing the two images it is quite apparent how varying the speckle size can affect the  resolution of the image – the hole clearly visible in Figure 2.4 (a) is nearly completely obscured  in Figure 2.4 (b).  2.3 A Typical ESPI Optical Setup Figure 2.5 shows a typical optical setup used in ESPI measurements.  There are several main  components: a coherent light source, a beam splitter, a phase stepper and a camera to record  the speckle images.  ESPI is typically used to evaluate full‐field surface displacements or  vibrations.  The general procedure for an ESPI measurement is to record the speckle pattern  created by splitting and recombining light from a coherent source.   Speckle patterns are  recorded before and after some mechanical displacement.  These speckle pattern recorded  second is modulated by surface displacements; by examining the phase change of the speckles  before and after displacement occurs, information about the surface displacement can be  inferred.       The inten illuminat   Applying   Fig sity distribu ion and refe  equation (2 ure 2.5 ‐ A t tion measu rence beam .2),  ypical ESPI  red on the C s as:  I  1I I I  optical setu CD can be e 2 1 2m mA A 2 1 22 cI I p (Adapted  xpressed as os   from [28])   the superposition of th 17  e  (2.4)  (2.5)  18    where   is the phase of the interference pattern created by the reference and illumination  beams and  1I and 2I  are the intensities of the illumination and reference beams respectively.   For a surface displacement  d , the intensity distribution becomes:     1 2 1 2' 2 cosI I I I I        (2.6)  where    is the phase change caused by to the motion of the surface.  A challenge of ESPI  measurements is how to determine this phase change   .  Since the only measurable quantity  is intensity, an additional technique is needed to evaluate the phase change of the speckles due  to the movement of the measured surface.  2.4 Phase Stepping The observed intensity distribution corresponding to the speckle pattern created by splitting  and recombining the coherent light source is described by equation (2.5). After surface  displacements occur, the phase difference between the two beams becomes:    1 2 1 2' 2 cos 'I I I I I      (2.7)  where:      ,      = the intensities of the illumination and reference beam respectively    '      = phase of interference pattern after surface displacement  By combining equations (2.5)  and (2.7), it can be seen that:    '       (2.8)  2I1I   Thus, th phase ch expresse   To solve  This is do Sets of im images b   With  n   be rewrit e phase can  ange due to d as:  for the phas ne by phase Fig ages, I and efore displa being a kno ten using th be measure  the displac e difference  stepping [2 ure 2.6 ‐ Ph  J are record cement, the wn set of ph e additive c d before ( ements   I  , a system o 9] schemat ase steppin ed before a  set of imag nI A ase steps, u osine identi ) and after   can be dete cosA B  f equations ically shown g schemati nd after the es can then cos(B    sually in 90 ty as:  ( ' ) the surf rmined.  Sim  in the form  in Figure 2 c (Adapted f  surface dis  be written  )n ° increment ace displace plifying, (2  of (2.9) mu .6:  rom [30])  placement.  as:  s.  Equation ment then  .5) can be  st be forme    Considering (  (2.10) can a 19  the  (2.9)  d.    the  2.10)  lso  20      cos sinn n nI A u v      (2.11)  where:  cos( )u B      sin( )v B     The system of equations can be re‐written as:    1 11 1 cos sin 1 cos sinn n n I A u I v d Gm                                  (2.12)  where:     Forming a least squares solution by minimizing the error yields:    min 0 T T error d Gm dd Gm d Gm dm G Gm G d           (2.13)  Simplifying equation (2.13) yields:    2 2 cos sin cos cos sin cos cos sinsin sin cos sin n n n n n n n n n n nn n n n n IA u I v I                                                    (2.14)  1 11 1 cos sin , , 1 cos sinn n n I A d G m u I v                                 21    Equation (2.14) is the general form of a least squares solution for a phase stepping scheme  involving n  steps.  In this work a common 5 step scheme [29] is used with  n =  [0°,90°,180°,270°,360°].  Using these values, equation (2.14) can be solved to yield the desired  phase difference between the illumination and reference beams as: [31]    4 2 1 2 3 4 5 7( )tan 4 6 4 I I I I I I I         (2.15)  By considering both phase stepped sets of images,   and  ' can be determined, yielding the  phase change due to the change in optical path length,        .  Finally, this change can be related  to actual surface displacement by using the sensitivity vector  sK  :     2 ssK d        (2.16)  where:     2k   and  1k  are defined as the propagation vectors of the illumination and reference beams  respectively, as illustrated in Figure 2.7.         2 nk    2 1s K k k    22                  Figure 2.7 ‐ ESPI setup including propagation vectors (Adapted from [28])    This phase difference can be plotted on the range of [0,1] forming a fringe pattern as:    1 (1 cos ) 2       (2.17)    Figure 2.8 (a) shows an example synthetic fringe pattern. The pattern can be visualized as a  topographic map of the surface displacement in the direction of the sensitivity vector, where  each pair of light and dark fringes represents a displacement of a wavelength of the light  source.    2.5 Fringe Patterns and Phase Unwrapping The phase difference evaluated in equation (2.15) lies on the range of [ , ]  .  Angles outside  this range are “wrapped” into the range by addition or subtraction of multiples of 2π.  To  determine the absolute phase difference due to change in optical path length   over the  entire image this data must be unwrapped.  Figure 2.8 (b) shows the phase difference along the  23    line shown in Figure 2.8 (a).  The phase experiences several jumps as the unwrapped phase  reaches the limits of [ , ]  .  An “unwrapping” process is required to add  2 where these  jumps occur to restore the unwrapped “true” value of the phase difference that is not  constrained on the range [ , ]  .  This unwrapped phase is then used to determine the surface  displacements.  Here, a common phase unwrapping algorithm [32] is used.  (a)   (I made these figures myself using synthetic fringes and MATLAB)  (a)  (b) (c) Figure 2.8 ‐ (a) Fringe pattern (b) Wrapped phase difference (c) Unwrapped phase difference  2.6 S  ESPI invo phase ch during al that the s           Figure 2. occupyin could mo case, the relative t surface d Speckle d disturban decorrela peckle D lves measu ange can be l measurem peckle mov 9 (a) illustra g the pixels  ve from the  pixels at th o the origin isplacemen ecorrelatio ces.  If the  ted, the ph ecorrel ring the pha  measured p ents.  Speck es to a diffe tes a typical marked 1‐4  upper left i e upper left  al speckle.  T t informatio n is the reas specimen sh ase informa F ation an se change o roviding th le decorrela rent pixel lo  speckle pat .  If the spec n Figure 2.9 then measu hus, the sp n is lost.    on that ESP ifts, even by tion is lost a igure 2.9 ‐ S (a) (a) d Its Eff f speckles c e speckles r tion occurs  cation.  tern.  One s imen displa  (a) to the u re a differe eckle phase I measurem  very small  nd the fring peckle deco ects aused by su emain at the when the su peckle in pa ces to the ri pper right  i nt speckle w  becomes de ents are so s amounts, th e pattern is rrelation rface movem  same pixe rface move rticular is m ght, the spe n Figure 2.9 hose phase correlated  ensitive to  e speckles b  destroyed.  ( ent.  This  l locations  s far enough arked ‘X’  ckle marked  (b).  In that  is random  and local  environmen ecome    b)  24   so   “X”    tal  25    Figure 2.10 ‐ The effect of speckle decorrelation   Figure 2.10 shows the effect of speckle decorrelation on an ESPI fringe pattern.  The white  circles and dimensions representing the hole drilled were added to the figure separately. Figure  2.10 (a) shows an undisturbed ESPI hole drilling fringe pattern where no rigid body motions  have occurred during the measurement. Figure 2.10 (b) shows the same fringe pattern if a  motion of one pixel (15 m ) to the right had occurred after drilling.  The pixel decorrelation has  introduced substantial noise and the quality of the fringe pattern is greatly decreased.  Figure  2.10 (c) shows the fringe pattern where the surface has moved 2 pixels (30 m ) to the right  after drilling.  The thin line in Figure 2.10 (c) is placed for scale and is 2 pixels wide.  In this case  the fringe pattern is completely destroyed with only noise remaining.  These figures clearly  illustrate how severely speckle decorrelation can affect ESPI measurements.  Even very small  motions of the surface can completely destroy the measurements.  Speckle motion can be caused by a number of factors.  These include air turbulence, mechanical  vibrations, specimen motion within its fixture or, in the case of hole drilling ESPI measurements,  the drilling process.  Simply touching the workpiece can even cause sufficient motion to  decorrelate the measurement.  Speckle decorrelation is a serious problem for ESPI  Undisturbed ESPI fringe  pattern   Shifted +1 pixel right after  drilling   Shifted +2 pixels right after  drilling   (a)   (b)  (c)  26    measurements that is currently preventing more widespread use of the method.  If specimen  motion could be counteracted, then the robustness of the ESPI method could be significantly  improved.  Being able to recover accurate measurements even in the presence of rigid body  motions would represent a solid step towards moving ESPI from a sensitive lab technique to  more robust field use.                             3 Dig 3.1 D Digital Im between from the displacem Figure 3. the came values re             ital Ima igital Im age Correla  a pair of im  first image  ent has oc 1.  Each pixe ra used.  0 r present sha   Figu ge Corr age Cor tion (DIC) is ages.  It is b is selected a curred.  Digi l holds a va epresents b des of grey. re 3.1 ‐ Dig elation relation  a popular o ased on the nd a match tal images a lue from 0 t lack, and th   ital image s & ESPI  ptical meth  principle of  is found wit re recorded o a maximu e maximum torage of a  Fringe C od used to d  image matc hin a secon  as intensity m value det  value repre painted spe orrecti etermine re hing; a tem d image afte  distribution ermined by  sents white ckle image  on lative moti plate or sub r some  s, as shown the bit dept .  Intermedi   27  on  set   in  h of  ate  28    A correlation coefficient which relates the template to its potential match can be defined as  [33]:          ( , ) ( , ) ( , ) x y c s t f x y w x s y t     (3.1)    Figure 3.2 ‐ The correlation coefficient (Adapted from [33])  where the terms are defined as shown in Figure 3.2 and the sums are computed over the pixel  area currently covered by the template,  ( , )w x s y t  .   ( , )f x y  represents the second image  where the template is being searched for, with size  M x N .   ( , )w x y  is the template, or feature  of interest, selected from the first image.  By moving a template sized area,  ( , )w x s y t  , over  29    some range of  sand  t  in , the distribution of the correlation coefficient can be built up.   For each position  ( , )s t  of the template, the correlation coefficient,  ( , )c s t , can be computed.   This coefficient will reach a maximum where the best match of  ( , )w x y  and   is found.   This process is illustrated in Figure 3.3.    Figure 3.3 (c) shows the correlation coefficient distribution, with higher values represented as  brighter areas.  The correlation coefficient is shown to reach a maximum at the location where  the template is matched correctly; however areas where poorer matches were found, around  the other characters, appear as well.  One of the difficulties with equation (4.1) is its dependence on the amplitudes of  and  .  If these functions change amplitudes the values of the coefficient will change  accordingly.  In order to remove this dependence, equation (4.1) can be normalized to form the  normalized correlation coefficient (NCC) [33]:    ( , )f x y ( , )f x y ( , )f x y ( , )w x y (c)(b)(a)  Figure 3.3 – Template matching using the correlation coefficient   (a) Base image  ( , )f x y (b) Template  ( , )w x y (c) Correlation coefficient distribution  ( , )c s t    (Adapted from [34])  Figure   Equation template falls on t template By repea displacem equation partially  padded w Another  features   3.4 ‐ Deter NCC  (3.2) remov  and area u he range of   from one im ting this pro ent data fo  (3.2), the m be placed ov ith blank va important c within the im mining tem ( , ) x s t    es the depe nder review [‐1,1].  By fin age to ano cess for ma r the entire atching pro er areas wh lues; howe onsideration ages.  Alth plate displa ( , ) x y y f f x y     ndence on   and norma ding the m ther can be ny small tem  image.  It s cess suffers ere there a ver accuracy  when perf ough image cement 2 ( , ) x y x y f f        amplitude b lizing to a un aximum of t  determined plates, it is hould be no  near the ed re no image  still suffers orming DIC  s with a var ( , ( , w x s y t w x s y      y removing  it magnitud he NCC, the , as shown  In this ca enclosin template using the displace image 1   possible to  ted howeve ge of the im  data.  In th  greatly.  measureme iety of featu 1/ 2 ) ) w t w      the mean in e.  This new  displaceme in Figure 3.4 se, the bou g the letter  .  By determ  NCC within ment   d  of  to 2 can be  build up a s r that due to ages since t ese cases th nts is the ch res such as  2   tensities of   coefficient nt of the  .    nded area  ‘A’ is the  ining its mo  image 2, th the feature  determined et of   the nature he template e images ar aracter of t grids, lines o 30  (3.2)  the   now  tion  e  from  .     of   will  e  he  r   dots can  informat several id is known If the tem to any on a perfect case, the determin In order t random s that no fa template 3.1.1 D A popula the surfa exhibit th informat simple sp speckles  not chan be used in D ion content entical max  as the corre plate, indic e of the pos  match will   motion of t ed as a scal o avoid this peckle patt lse positive  displaceme igital Ima r choice for  ce under inv e desired q ion content  ray can is su are physical ge significan IC, the best .  If a period ima as the t spondence ated by the itions indic be found.  T he template ar multiple o  problem, h erns are oft s are report nt solution  ge Correl image corre estigation.  ualities for D so long as t fficient.  Th ; they are d tly as the su  results are  ic pattern is emplate is p  problem in   dark box, is ated by the  herefore, in  can only be f the patte igh informa en used to e ed and that is unique.  ation Usin lation is a ra  An example IC measure he size of th e advantage urable and h rface unde obtained us  used, it is p eriodically  image matc  moved  vectors,   this    rn size.    tion  nsure   the  g Painted ndomly spr  is given in  ments – the e paint spec  of a painte ave no spat r investigatio F p ing random ossible for t matched ac hing and illu  Speckles ayed painte Figure 3.1.   y are rando kles is contr d pattern li ial depende n moves, u igure 3.5 ‐ T roblem   patterns wi he coefficie cording to th strated in F  d speckle pa Painted spe m, and have olled.   Ofte es with the f nce.  Their p nless very la he corresp th high  nt to registe e pattern.   igure 3.5.  ttern applie ckle pattern  high  n paint from act that the ositions als rge  ondence  31  r  This  d to  s   a    o do  32    deformations occur.  This durability and consistency aids in producing accurate DIC results.   Because of the desirable characteristics, recently the use of painted speckle DIC has been  applied within various areas of experimental mechanics including tensile testing [9] , strain field  measurement [11] and residual stress measurement. [35,36]  3.1.2 Digital Image Correlation Using Laser Speckles The characteristic speckle pattern created using a laser light source also contains many of the  same characteristics as a painted speckle pattern – the pattern is highly random, unique and  has dense information content.  All of these qualities make laser speckle images an attractive  choice as the basis for DIC measurements.  Laser speckles also have the added advantage of not  needing physical application; all that is required is a coherent monochromatic light source and  the object to have a suitable surface finish with roughness on the scale of the source  wavelength. The main disadvantage of using laser speckles is the impermanence of the pattern.   As the surface moves, the speckles move along with the surface.  However, the local  illumination gradually changes from place to place, and so the pattern of the speckles  correspondingly changes.  In addition, the phase of the speckle pattern can change, causing  light speckles to become dark, and vice‐versa.  These behaviours present obvious problems for  DIC measurements as the template matching process suffers due to both the change in phase  and speckle shape.  Despite these challenges, laser speckle DIC measurements have been  applied in similar areas to that of painted speckle DIC, with the majority of the applications  aimed at strain and displacement measurements [10,24]. Due to the influence of the surface  features on the laser speckle pattern, additional applications involving surface condition  33    monitoring [25] have also recently been explored.  These measurements involve monitoring the  correlation coefficient for a certain area as the surface condition changes.  3.2 ESPI Fringe Correction Procedure Currently due to their inherent sensitivity to relative motions and environmental disturbances,  ESPI measurements are limited largely to laboratory use in very controlled environments.   Specimen rigid‐body motions cause speckle decorrelation, which corrupts the critical phase  information necessary to determine the surface displacement.  The ESPI fringe correction  procedure proposed here represents a method to correct rigid‐body motions by re‐correlating  speckle images and therefore restoring the ESPI fringe map.  A similar idea was first suggested  by Reu and Hansche [16] but was never pursued.  The procedure makes use of DIC methods to first determine the rigid‐body motions that have  occurred during the ESPI measurement.  This is done using the captured ESPI image sets as the  basis of DIC displacement measurements.  By using these laser speckle images, the rigid body  motions that occurred between the image sets can be determined.  Once these motions are  known, the affected ESPI image set can be mathematically shifted in order to ensure the  speckle phase information is re‐correlated. Figure 3.6 shows an outline of the fringe correction  procedure.    In order t are selec used.  Th in the eq the searc represen search ar the maxi for one s   o determin ted circumfe e search ran uivalent for h area.  Figu ts a separat ea by shiftin mum value  ubset.  Figure 3. e the rigid b rentially ar ge can be s m of the Pea re 3.7 show e subset.  Fo g the temp of the NCC i 6 ‐ The ESPI ody motion ound the ce et according rson produ s an examp r each subs late pixel by s found.  Fig  fringe corre s of the spec ntre of the m ly, and the  ct‐moment  le of how th et, a NCC di  pixel over t ure 3.8 show ction proce imen, corre easureme normalized  correlation  e subsets a stribution is he search a s a NCC dis dure  lation subse nt; a 60x60  cross correl coefficient,  re selected ‐  built up ove rea.  From t tribution an ts, or patch pixel subset ation coeffic is formed ov  each box  r the select hat distribut d its maxim 34    es,   is  ient,  er  ed  ion,  um               Figure Figure 3.7 ‐  3.8 ‐ NCC d DIC subset s istribution  election  and maximum 35   Forming  achieve s surround constant   The NCC  pixel whe         Since the ordinates the NCC dist ub‐pixel res ing 8 pixels. s, a‐f:  maximum a re the max    pixels are u  of each of  Fig ribution usi olution, a 2   Starting w nd surround imum value  nit spaced,  the 9 pixels  ure 3.9 ‐ NC ng equation D quadratic  ith the gene z a bx   ing pixels a of the NCC  a matrix eq as:  C maximum  (3.2) allow function is f ral form of a 2cy dx ex  re defined a is located.  uation can b  and surrou s matching t it to the ma  2D quadra 2y fy s shown in  e formed by nding pixel o integer pi ximum valu tic, there ar Figure 3.9, w  considerin s xel values.  T e and the  e 6 unknow ith 5 being  g the co‐ 36  o  n  (3.3)  the  37      1 2 3 4 5 6 7 8 9 1 1 1 1 1 1 1 0 1 0 0 1 1 1 1 1 1 1 1 1 0 1 0 0 1 0 0 0 0 0 1 1 0 1 0 0 1 1 1 1 1 1 1 0 1 0 0 1 1 1 1 1 1 1 z z a z b z c z d z e z f z z                                                                  (3.4)  This system is overdetermined since there are 9 equations and only 6 coefficients.  Here, it is  chosen to fit the coefficients a b c d f exactly to the central five points 2 4 5 6 8 and coefficient e  to the average of the corner points 1 3 7 9.  This procedure focuses the fit on the central points.   The coefficients can then be determined by addition and subtraction as:    5 6 4 8 2 4 6 5 1 3 7 9 2 8 5 ( ) / 2 ( ) / 2 ( ) / 2 ( ) / 4 ( ) / 2 a z b z z c z z d z z z e z z z z f z z z                  (3.5)  Once the quadratic is formed, directional derivatives are calculated in order to find its  maximum.  Continuing from equation (3.3):    2 0 2 0 z b dx ey x z c ex fy y             (3.6)   Equation relative t   These res Once the according Mathem speckle c they orig  To summ develope process d Figure s (3.6) can f o the maxim ults are add  relevant rig ly.  Bi‐linea atical shiftin orrelation.   inally did wi arize, digita d.  This algo iscussed.  T  3.10 ‐ Info image inally be sol um pixel as ed to the m id body mo r interpolati g of the ima The speckle thin the firs l image cor rithm will b he aim of th rmation los  shifting  ved to yield  :  x y   aximum int tions are de on is used t ge sets to a s in the seco t set.  This p relation has e used as th e process is s due to  the co‐ordin 2 2 2 4 2 4 bf ce df e be cd df e      eger pixel v termined, th o account fo ccount for t nd set of im ixel re‐regis be recove solution to edge, ther shifted are possible th the loss of can becom  been invest e main com  to mitigate ates of the alue to allow e affected  r the sub‐p he rigid‐bod ages now o tration allow red and the   be largely  e is no infor a, so image ere.  Figure  pixels near  e quite larg igated and  ponent of t  the negativ  maximum o  sub‐pixel i set of ESPI im ixel image sh y motions r ccupy the s s accurate  surface disp restored.  A mation to m  recorrelatio  3.9 shows a the image e e for large i a basic algo he ESPI fring e effects of  f the NCC  nterpolation ages is shi ifts.    e‐establishe ame pixels a phase data  lacement  t the image  ove into th n is not  n example  dge   This lo mage shifts. rithm  e correction specimen ri 38  (3.7)  .   fted  s  s  to  e  of  ss      gid  39    body motions within ESPI measurements.  If this fringe correction process can be proven  effective, the robustness of ESPI measurements will be significantly improved, potentially  allowing measurements to be performed in more adverse environments as well as other  related benefits.      40    4 Experimental Validation 4.1 Hybrid ESPI – DIC Apparatus An experimental apparatus capable of making hole drilling ESPI measurements was constructed  to test the feasibility and effectiveness of the proposed DIC fringe correction procedure.  Figure  4.1 illustrates the layout and major components of the experimental apparatus.    Figure 4.1 ‐ Apparatus schematic    Table 4.1 lists the major components of the apparatus. The ESPI portion of the  apparatus – encompassing the beam splitter mirror and piezoelectric actuator, was custom  designed and built.  Chip suction was also implemented within the drilling apparatus to   minimize impact o illustrate           F  the spread  n optical co  the apparat igure 4.2 –  of the meta mponents a us.  ESPI beam s llic particles nd measure plitter asse Figure 4.  cut from th d surfaces. T mbly 4 ‐ Drill asse e hole and  he photos i Fig mbly to reduce th n Figure 4.2 ure 4.3 – A eir adverse  ‐ Figure 4.4 djustable m 41  irror 42    Table 4.1 ‐ Apparatus components  Component  Make/Model  Specifications  Laser  JDS Uniphase  CDPS532M  532nm wavelength, 50mW  Camera  Prosilica EC 750  640 x 480 resolution, 8 bit grayscale color depth  with telecentric lens  Piezoelectric  Actuator  Piezomechanik PSt  150/7x7/7  9µm maximum extension  3‐Axis Stage  Daedal  Drill and specimen mounting possible  High Speed Drill  NSK EMS‐3041  25k rpm max speed  Drill Stage  Newport CMA‐25CCCL  Spring return, lead screw driven    The objective of the apparatus design was flexibility.  Each beam can be adjusted  independently, as can the vertical mirror placed near the specimen.  This allows the area of the  specimen under investigation to be uniformly illuminated by coherent beams oriented with the  proper incidence angles.  By repositioning and adjusting the components, a range of specimens  can easily be accommodated.  Attaching the selected specimen directly to the 3 axis table  allows known specimen rigid‐body motions to be accurately created. This allows the  effectiveness of EDIC as well as the DIC fringe correction procedure to be explored  systematically.  4.2 ESPI System Validation Upon completion of the ESPI apparatus, it was necessary to first validate that it provided  accurate measurements.  Since residual stress hole drilling measurements are used as the test  case for DIC fringe correction within this work, measurement of a specimen with a known  residual stress field was examined.     The spec National  By measu comparin Figure 4. calibrate Although drilling re work dea provide E described imen used w Laboratory, ring the spe g the result 5 shows the d specimen.  frequently  sidual stres ls with ESPI SPI fringes.  in [12].    as a ring an  New Mexic cimen usin s to those o  ESPI fringe    Figure 4.5 used to com s solution w  measureme   All residua d plug spec o.  This spec g the newly  btained in L result from   ‐ Fringe res pare the ac ill not be di nts; residua l stresses ca imen provid imen had p constructed os Alamos t a hole drillin ult using ca curacy of ES scussed in d l stress hole lculated are ed by Dr. M reviously be  experimen he new app g measurem librated spe PI fringes w etail within   drilling is o  done so exa . Steinzig of en calibrate tal apparatu aratus could ent perfor cimen  ithin this re this thesis.   nly used as  ctly in the m  the Los Ala d at Los Ala s and   be verified med on the    search, the  The focus o a test case t anner  43  mos  mos.   .   hole  f this  o  44    The induced stress field within the specimen according to the measurements and calculations  performed in Los Alamos is compared with the solution measured using the experimental  apparatus in Table 4.2.     Table 4.2 ‐ Calibrated vs. measured stresses             The measured stresses agree with the calibrated values within about 4MPa.  This is an  encouraging result and gives credibility to the experimental system.  The stress difference is  modest compared with practical residual stresses, which can be in the range of some hundreds  of MPa.   Expressing the difference as a percentage of the yield strength of the material gives a  more relative indication of the scale of the differences in the measurement.  These errors are  all less than 2%, and this again indicates the measurements are in good agreement.  Comparing the two results they are shown to be in good agreement, differing by only a few  MPa when comparing each of the three in plane stress components.  Since the capacity of the  constructed system to measure residual stresses accurately has now been proven, the  investigation of ESPI fringe correction using this apparatus can be explored with greater  confidence.    x y  xy Calibrated Stress (MPa) 35.6 57.5 ‐0.3  Measured Stress (MPa) 31.4 53.2 ‐2.9  Difference (MPa)  4.2 4.3 2.6  % Difference w.r.t.  yS (%)  1.5 1.5 0.9  45    4.3 DIC Measurements Using ESPI Images Examining a typical ESPI measurement, a large amount of data is available for analysis.  When  using a 4‐step phase stepping algorithm, 8 total images are therefore available for each ESPI  measurement; 4 captured initially ‐ the reference set, and 4 captured after some mechanical  change ‐ the fringe set.  Each of the images within each set is phase stepped according to the  scheme described in section 2.4.  Due to the phase dependent nature of laser speckles, this  causes the images within each set to vary significantly.  Figure 4.6 shows sections from a set of  4 phase‐stepped images.  These sections are taken from the same location of each image.  The  phase‐ stepping among these images causes them to appear very different, even though they  share the same underlying speckle pattern.  Speckles are seen to change phase from light to  dark as well as to change shape from image to image.        This image variation becomes the first challenge when attempting to utilize ESPI images within  DIC measurements.  When examining an arbitrary pair of images, one from the reference set  and one from the fringe set, it is possible that the speckles could be in phase, out of phase, or  somewhere in between.  For successful DIC measurements it is imperative that the pattern  being matched remains consistent from image to image.  Because phase variation has a direct  effect on intensity, using these phase‐stepped images directly is not suitable for DIC image  Figure 4.6 ‐ Sections of a phase stepped set of images  46    analysis.  Therefore, it is essential to develop some means of mitigating this pattern fluctuation  due to phase variation.  Returning to equation (2.10), the intensity distribution created by the illumination and  reference beams within a phase stepped image can be expressed as:        cos( )n nI A B       (4.1)      In order to construct images suitable for DIC analysis, the quantities  A  and  B  which are seen  to be phase independent, are of interest.  By solving equation (2.18) in the case of a 4 step  algorithm, these quantities can be expressed in terms of the four phase‐stepped image  intensities at a given pixel as:    2 2 4 2 3 11 2 3 4 (     )     (     ) 4 2 I - I I II I I IA B        (4.2)   Therefore,  A and  B ’images’ can be constructed from the 4 images within a phase stepped set.    To test the feasibility and effectiveness of utilizing these new phase independent images within  DIC measurements, statistical tests of DIC capabilities were conducted on an example specimen  mechanically displaced by eighteen 0.1mm increments up to 1.8mm using the three‐axis table.   In order to minimize the backlash errors within the table, the specimen was moved backwards  and forwards each time before recording images to take up slack in the table drive system.  Sets  of four phase‐stepped ESPI images were measured after each displacement.  For the  47    0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0 20 40 60 80 100 120 St an da rd  De vi at io n ( pi xe ls ) Image Shift (pixels)  "A" Image Min. Image Max Image "B" Image magnification used, a 0.1mm displacement corresponded to 6.1 pixels.  DIC evaluations of the  image displacement were done at one hundred 50x50 pixel patches in a 10x10 grid spanning  the images.  The standard deviations of these DIC evaluations and the average correlation  coefficients  2R  within the patches gave indications of the statistical quality of the results.   Figure 4.7 shows the standard deviation of the DIC displacement estimates vs. amount of image  shift.  A lower standard deviation indicates less statistical uncertainty and hence greater  expected computational precision.  In general, the standard deviation increases with image  shift, indicating that the speckle pattern does not simply translate, as would happen with  images of displaced physical surface features.  ESPI speckle patterns from displaced surfaces  slowly change due to differences in the illumination at different points in space.      Figure 4.7 ‐ Standard deviation of DIC displacement estimates vs. amount of image  shift 48      The lowest of the four lines in Figure 4.7 shows the DIC results using the constant “A” from  equation (4.2).  The adjacent bold line shows the results from “B” image, which are comparable  to the “A” results but deteriorate at higher image shifts.  The “Min. Image” and “Max. Image”  lines respectively show the minimum and maximum standard deviations observed among the  16 possible combinations of four initial stepped images with four final stepped images.  The  “Min. Image” results derive from the ESPI images that happen to be in phase, and give results  similar to those from the “A” images.  However, it is not known in advance which particular ESPI  images are in phase.  Thus, it would be necessary to test all image combinations to find them,  Figure 4.8 ‐ Average correlation coefficient R2 vs. amount of image shift  0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 20 40 60 80 100 120 Co rr el at io n C oe ffi ci en t R 2 Image Shift (pixels)  "A" Image Min. Image Max Image "B" Image 49    which is not a very practical procedure.  The “Max. Image” curve represents case when the  chosen images are out of phase, and thus are poorly correlated.  The large “Max. Image”  standard deviations at small image shifts occur because the speckle patterns remain mostly  intact in this range, and so adverse relative phase has a large effect.  At larger image shifts, the  speckle patterns become distorted, and so the additional distortion caused by adverse phase is  no longer so influential.  Figure 4.8 shows the average correlation coefficient  2R  for the same correlation patches as in  Figure 4.7.  Higher  2R  is desirable, and thus the trends in Figure 4.8 support the conclusions  from Figure 4.7.  The correlations decrease with amount of image shift, showing the changes in  speckle patterns with surface displacement.  The  2R  values are much lower than the range 0.95  to 0.99 typically achieved with images of physical surface features [8].  This lower correlation  occurs because speckle patterns respond to local changes in illumination as well as to surface  displacements.  Image “A” in Figure 4.8 achieve the best (highest) correlation coefficients;  image “B” and the raw ESPI images are consistently lower.    Thus, based on both Figure 4.7 and  Figure 4.8, the general use of “A” images for DIC analysis is indicated. Based on this result  whenever an ESPI based DIC (EDIC) calculation is performed, an “A” image will first be formed  for each phase stepped set of images and then the DIC method will be applied.         4.4 M S Having n measure is to inve these qu involving that the s effect of  the table inherent  selected  used to a edges of  an averag converte Figure 4. to des easurin peckle S ow establish ments, the f stigate the r alities comp  convention ize of the s laser speckl  was run ba to the table as before on void measu the image.   e recorded d into pixel  9 ‐ Coordina cribe specim g Surfac ize ed an effec ringe correc ange and ac are to those al painted [ peckles used e size on ED ckwards and ’s drive syst  a circular g rements los Using these .   The DIC r units using a te system u en motion e Motion tive method tion proces curacy of th  of physical 37] and lase  affects me IC measurem To thi a flat s Consid specim move both t  forwards f em.  For eac rid as show ing correlat  subsets 18  esult was th  calibration sed  .  s Using  of using ph s is one step ese measur  feature bas r [22] speck asurement  ents will a s end, a seri pecimen m ering a co‐o en as show d from an in he x and y d or each set i h phase ste n in Figure 3 ion due to th individual D en compare  done by re  EDIC & ase stepped  closer to re ements.  Of ed DIC.  In a le based DIC accuracy.  B lso be invest es of experi ounted on t rdinate sys n in Figure  itial position irections in n order to m pped set, 18 .6.   The cen e search ar IC measurem d with the a cording ima  the Effe  ESPI image ality.  The n  particular i ddition, pre  measurem earing this i igated.    ments was c he 3‐axis m tem oriente 4.9, the spe  over a rang  0.1mm incr inimize the  60x60 pixe tral area of ea being co ents were  pplied displ ges of a prec ct of s within DIC ext logical s nterest is ho vious work  ents sugges n mind, the  onducted u otorized tab d on the  cimen was  e of 3mm i ements.  Ag  slack errors l subsets we  the image w nstrained by performed  acement,  ision scale  50    tep  w  ts  sing  le.   n  ain,    re  as   the  and  51    placed on the surface of the specimen.  The entire process was then repeated for a range of  laser speckle sizes.  Recalling equation (2.3), the variation in speckle size was achieved by  changing the f# of the lens of the CCD camera.   In order to compensate for the change in light  incident to the CCD due to this change, the exposure time of the camera was adjusted  accordingly to achieve a consistent average intensity over the entire image.  The camera  parameters for the range of speckle sizes examined are listed in Table 4.3.  The speckle size  listed represents an estimate calculated using equation (2.3).   Table 4.3 ‐ Camera parameters and speckle sizes    A higher f# represents a smaller aperture and therefore larger speckles.  By closing the aperture  down smaller speckles are unable to form and all that remains are the larger ones [38].  In  addition to laser speckles, a painted speckle pattern was also examined.  This pattern was  created directly on the surface of the specimen using a can of conventional black spray paint.   When recording images of the painted pattern, only a single image was required as no phase  stepping occurred. Specimen illumination was accomplished using diffuse white light from a  fluorescent lamp.  The experimental results for specimen shifts in the x direction are shown in  Figure 4.10.  The ideal results obtained from the table itself are shown as a line with unity slope.   The accuracy of the various EDIC measurements can then be measured by how closely the  f   5.6  8  11  16  22  27  32  Exposure time (ms)  1  3  7  18  40  70  125  Speckle Size,  S , (px/s)  0.61  0.87  1.19  1.73  2.38  2.92  3.46  Mean Intensity,  I   101  108  104.7  101.3  100.1  101.4  101.2  52    corresponding data matches the table data.  Examining Figure 4.10, it becomes very evident  that speckle size has a significant effect on the effectiveness of the EDIC measurements.  Both  the largest and smallest speckles, 3.46px/s and 0.61px/s respectively, perform poorly.  The  large speckles maintain a reasonably accurate result until a displacement of 82.6 px, or 1.4mm.   Beyond that point the solution deteriorates rapidly and deviates from the true displacement  more than any speckle size tested.  On the other end of the size spectrum, the smallest speckles  lose accuracy at even smaller displacements, deviating at displacements of 70.9 px or 1.1mm.   When compared to the largest speckles, the smallest speckles tested did provide a closer  Figure 4.10 ‐ X Direction EDIC Measurements  ‐100 ‐50 0 50 100 150 200 0 50 100 150 200 M ea su re d D is pa lc em en t (p x) Table Displacement (px) Table 0.61 px/s 0.87 px/s 1.19 px/s 1.73 px/s 2.38 px/s 2.92 px/s 3.46 px/s Painted 53    agreement with the applied displacements however.  Although neither the smallest or largest  speckles performed well, the results improved within the middle of the speckle size range.   Speckles 0.87px/s in size, the next smallest speckles, began deviating at displacements of  94.5px followed by the second largest speckles of 2.92px/s at 129px displacement.  It is evident examining the data that the trend is towards an optimum speckle size that  maximizes the EDIC effectiveness.   Consistent with this trend, the datasets representing the  midsize speckles, 1.19px/s and 1.73px/s are shown to perform the best, maintaining good  agreement with the table displacement right up to the maximum displacement measured of  3mm or 177.2px.   Over this range the calculation error remains sub pixel except in a few  outlying cases.   The painted speckles are shown to give accurate results, matching the applied displacements  very closely.  This increased accuracy is attributed to the fact that painted speckles are physical  features that move with the specimen surface.   They are subject to a different limitation  caused by the possible presence of geometrical distortions due to aberrations in the imaging  optics.  However, this is a minor effect with good quality optics.  This quality is in stark contrast  with laser speckles, which as previously discussed in section 2, do have a spatial dependence  relating to the phase of the interference.  This causes the speckle pattern to change as the  specimen moves, even though the phase dependence has been mitigated by using the  technique described in section 4.3.  Since DIC measurements rely on the consistency of the  pattern being matched this difference indicates that painted speckles will always yield more  accurate results than laser speckles.   This is confirmed in the experimental results.   54    When examining EDIC measurements in the y direction, much the same trend is observed.   These experimental results are shown in Figure 4.11.  Again, the largest and smallest speckles  lose accuracy the fastest but in this case 2.92px/s is the next speckle size to fail as opposed to  0.87px/s as in the x direction.  The general trend however remains the same, with 1.19px/s and  1.73px/s again performing the best of those speckle sizes tested.  It should be noted that the  last data points shown in Figure 4.11 are especially inaccurate because the search area of  several DIC subsets was curtailed by the bounds of the image causing the average  measurement to be severely affected.  This is due to the 4:3 aspect ratio of the image itself; the  increased data available in the x direction allows displacements to be measured at greater  distances than in the y direction.  This trend is clearly visible when examining the EDIC results.   Although painted DIC remains accurate to 2.8mm this is most likely due to the greater strength  of correlation when using painted physical features when compared to the spatially dependent  laser speckles.  Even with the subset being only partially matched, the correlation using painted  speckles is strong enough to maintain a reasonably accurate result.    55         Figure 4.12 shows the measurement range of each speckle size in both the x and y directions.   The trend towards an optimal speckle size around 1.19px/s and 1.73px/s, where the  measurement range is greatest, is clearly shown for both directions.      ‐150 ‐100 ‐50 0 50 100 150 200 250 0 50 100 150 200 M ea su re d D is pl ac em en t (p x) Table Displacement (px) Table 0.61 px/s 0.87 px/s 1.19 px/s 1.73 px/s 2.38 px/s 2.92 px/s 3.46 px/s Painted Figure 4.11 ‐ Y Direction EDIC Measurements 56    To examine more clearly the effectiveness of the most accurate speckle sizes F11 and F16 as  Examining the data for 1.19px/s, 1.73px/s as well as the painted speckles, the absolute average  error and total error were calculated as shown in Table 4.4, omitting the last measurements in  the y direction because of the aspect ratio issue.    Table 4.4 ‐ EDIC measurement standard deviation for painted and laser speckles    1.19px/s 1.73px/s Painted Average Error (px) 0.60  0.58  0.46  Total Error (px)  32.01  30.56  25.10    0 20 40 60 80 100 120 140 160 180 200 0 0.5 1 1.5 2 2.5 3 3.5 4 M ea su re m en t R an ge  (p x) Speckle Size (px/s) X Direction Y Direction Figure 4.12 ‐ Measurement range vs. Speckle Size 57    These data yield much clearer results showing that the laser speckle size of 1.73px/s, where  16f   yields speckles that are the most accurate when used within EDIC measurements.  Using  this size the average error of 0.58px is slightly lower when compared with the average error of  0.60px when using speckles 1.19px in size.  This trend is also echoed when considering the total  error across all measurements.  At 1.73px/s, the average speckle has a size that is larger than  the conventional size of 1px.  Also, at 1.2px/s the EDIC error is very similar.  This suggests  speckles 1‐2px in size are the most effective for EDIC measurements.  Painted speckles again, as  expected, are shown to be more effective than laser speckles, with both an average error and  total error lower than either laser speckle size.  The average error of 0.46px for painted  speckles is relatively high when compared with typical error associated with the DIC method,  however errors in the range of 0.001 to 0.5px have been reported within previous work [8,23],  which is consistent with this result.    The deterioration in DIC correlation accuracy at the extremes of speckle size conforms to  theoretical expectations.  For a given patch size, DIC correlation accuracy varies approximately  inversely with the average linear size of the features within the analyzed images.  Thus, the  larger speckle sizes associated with larger  f  give lower correlation accuracy.  The DIC accuracy  increases as the speckle size reduces with smaller f .  However, with the smallest values of  f ,  the average speckle size reduces below the CCD pixel size and cannot be adequately resolved  by the camera, thereby also reducing correlation accuracy.    This concept was investigated for single beam DSC measurements by Feiel and Wilksch [24]  who, based on the work of Sjödahl and Benckert [22] proposed an optimum f number be  58    pd  M chosen according to the following criteria, in order to minimize aliasing errors while maintaining  the highest resolution possible:    2 (1 ) pdf M     (4.3)  Where:  = Physical pixel size  = Magnification      = Source wavelength  Equation (4.3) is based on the Nyquist criterion, which states that in order to accurately  measure a signal containing a maximum frequency of Y , the signal must be sampled at a  frequency of  2Y .  Physically, the speckle pattern contains speckles at a frequency of the  inverse of the speckle size as defined in equation (2.3).  The sampling frequency in this case is  defined as the inverse of the pixel size.  By this reasoning, equation (4.3) was developed to yield  the corresponding f number in line with the Nyquist criterion and speckle size as defined by  equation (2.3).  Using equation (4.3) for the case of the experimental apparatus, the optimum  18.49f  shows good agreement with the experimental result of  16f  . This indicates that  EDIC measurements share similarities to single beam DSC measurements.    The particular DIC algorithm used here has been applied in recent work [39] with good results,  so the errors in displacement measurement are most certainly due to physical imperfections  within the measurement setup.  Determining the field of view using a scale introduces certain  inaccuracies inherent to the scale itself.  Due to physical constraints the camera itself had to be  placed at a slight angle to the specimen, decreasing the CCD surface sensitivity to planar  displacements.  Finally, the 3‐axis table itself is not perfectly accurate and the lead screw drive  59    system contains inherent backlash errors. DIC accuracy could be improved by using a higher  resolution camera with a smaller physical pixel size and greater bit depth with a higher  magnification lens [8]. In the fringe correction application discussed here, DIC is used only to  evaluate bulk specimen motions but not to measure specimen surface displacements directly.   Therefore, for this application, the lower accuracy achieved using ESPI images is acceptable.   Although more complex DIC algorithms exist [23], this particular one is attractive for its ease of  implementation as well as speed.  A typical DIC measurement encompassing 18 individual  subsets was completed in a few seconds on a 2.14 Ghz processor once the images were  recorded.    Through these experiments three key results are obtained.  First, using the constructed  experimental apparatus, surface displacements can be measured accurately with sub pixel  resolution using EDIC in excess of 3mm or 177.2px in the x direction and up to 2.1mm or  137.0px in the y direction.  This difference in measurement range can be attributed to the  aspect ratio of the images used as well as the optical setup.  Second, the size of the laser  speckles used has a very significant affect on the accuracy and range of the EDIC measurement,  with the best results being obtained for this setup using an f number of 16 which corresponds  to an average speckle size of approximately 1.73px.  Third, using painted speckles consistently  produces more accurate DIC results than the proposed EDIC method.  By comparison, using this  experimental apparatus, both painted DIC and EDIC measurements produce errors on the same  scale.  Considering the accuracy when using painted speckles, it is on the lower end of the  reported range.  This is mainly due to various inaccuracies inherent to the apparatus and  method itself.  The slightly lower accuracy of EDIC should not be a significant concern within the  60    ESPI fringe correction method since the hole drilling surface displacements are measured using  ESPI and DIC is used only to re‐correlate the ESPI images.      4.5 Null Fringe Correction Previous experiments have shown that EDIC measurements are capable of measuring surface  displacements of significant distances with sub‐pixel accuracy.  This is an encouraging result and  bodes well for the fringe correction process.  A simplified fringe correction is now investigated –  one in which no measurement is actually performed and the data is blank.  Ideally if the fringe  correction is successful, the recorrelation process will yield a completely null fringe.  These  experiments will provide a good starting point for examining the viability of the fringe  correction process.    The procedure used for these tests is similar to that used in section 4.4.  Using a nominally flat  specimen, a set of phase‐stepped images is first recorded to serve as the reference set.  The  specimen is then displaced a known amount using the table and a second set of fringe images is  recorded.  The table displacement is calculated using an EDIC solution Based on 18  circumferentially spaced 60x60 pixel patches.  Finally the fringe set of images are shifted to  reverse the table motion and re‐correlate the laser speckles.  Figure 4.13 shows an example set  of measurements.  61    (a)  (b) (c) Figure 4.13 ‐ (a) Baseline zero fringe (b) Corrupted fringe after 0.5mm displacement (c) Corrected fringe after image shifting  Figure 4.13(a) shows the baseline zero fringe obtained by recording consecutive sets of phase  stepped images without moving the table.  In the ideal case this fringe would be completely  blank; however slight texture created by measurement noise is visible.  This noise is caused by  the superposition of a random speckle pattern on a regular CCD pixel grid. Figure 4.13 (b) shows  the fringe obtained after the specimen has been shifted 0.5mm in the positive x direction.  As  expected, the speckle correlation is completely lost and the resulting fringe pattern is  destroyed.  Finally, after performing the EDIC measurement, the fringe images are shifted  accordingly to correct for the 0.5mm displacement. The resulting fringe is shown in Figure 4.13  (c), which shows that the speckle correlation has been restored, but with a large phase gradient  superimposed across the image.  This initial result shows promise – speckle correlation can  indeed be recovered, however the cause and correction of the resulting phase gradient needs  to be determined.      4.5.1  P Upon inv visible in of severa beam ex source, w reference one of th corrected perfectly conical b 4.14.  If t d , the p shown in beam 2 d therefore    The bea being con in Figure  hase Gra estigation, t  Figure 4.13  l factors.  W pander is pla hich makes  beams con e factors th  ESPI fringe  flat specim eams at ang he specimen ath length o  Figure 4.15 ecreases by  changes by ms can be id ical in shap 4.16. The co dients – C he cause of (c) was foun ithin the ex ced directly  both the ill ical in shap at causes ph s.  Consider en imaged a les  1  and   is then shi f beam 1 inc  and similar   2sind  .  T :   ealized as o e due to the ne angles o ause and  the phase g d to be a co perimental   in front of  umination a e.  This coni ase gradien  the ideal ca symmetrica 2as shown fted right by reases by d ly the path l he relative  (d  riginating fr  beam expa f the beams  Correctio radient  mbination  apparatus, a the laser  nd  cal shape is  ts within  se of a  lly by two   in Figure   a distance  1sin as  ength of  path length 1sin sin  om point so nder and ill  can be des n     2) urces at a d uminating a cribed as  Figure 4. illuminat undergo Figure 4 beam 1  istance L fro n image  D   14 ‐ Conical ion of a flat ing a shift d .15 ‐ Path le due to spec m the surfa wide as sho  asymmetri  specimen    ngth chang imen shift 62  (4.4)  ce,  wn  cal  e of  63      1 1 2 2cos , cosD DL L       (4.5)  Because of the conical shape of the beam, the  beam angles, or angles of incidence, vary across  the image.  Considering only a one‐dimensional  image, the beam angles on the left side will be  1 2 1 2,2 2     and on the right they will be  1 2 1 2,2 2     respectively.  Considering the  change in path length at the left and right side of the image yields:        1 2 1 2 1 2 1 2 sin 2 2 sin 2 2 L R d d                                           (4.6)  Finally calculating the total path length based on equation (4.6) and simplifying using  trigonometric identities yields:    2 12 12 cos sin cos sin2 2L R d             (4.7)  For small cone angles, and substituting equation (4.5) this yields:     2 22 1cos cosL R dDL        (4.8)  Utilizing trigonometric identities, equation (4.8) simplifies to:  Figure 4.16 ‐ Conical  beam illumination  64         1 2 1 2sin sinL R dDL           (4.9)  Phase difference can then be related to the change in path length using the wavelength of the  source as:      2L R L R               (4.10)  Simplifying equation (4.10) yields:       1 2 1 22 sin sindDL          (4.11)  Where:   = Specimen shift distance  = Image width at specimen  = Source wavelength  = Beam path length  = Beam 1 angle of incidence      = Beam 2 angle of incidence  The phase gradient due to the conical illumination beams and angles of incidence can be  calculated using equation (4.11).  To validate this equation, an experiment was performed  where the flat mirror that reflects the illumination beam was mounted on a precision turntable.   By rotating the mirror set angular increments, the angle of incidence of the reference beam was  varied.  Sets of images were recorded before and after a constant specimen shift of 0.25mm in  the x direction.  The fringe set of images were then shifted back into alignment according to the  EDIC result and the number fringes measured numerically using a least squares solution in the  same manner described later in this section.  Figure 4.17 shows the results of the mirror turning  experiment.  d D λ L 1 2 65    As shown in the figure, the experimental data matches the results from equation (4.11) quite  well.  The difference between the two results may partly be due to the difficulty of accurately  identifying the illumination distance “L” because the position of the source focal point could not  be exactly identified.        The remaining phase gradients are caused by the change in the relative path length of the  illumination and reference beams from one side of the image to the other.  These changes in  path length can be attributed to a combination of factors.  Any specimen curvature in the z  direction will cause phase gradients as the specimen is moved since the relative path lengths of  ‐5 ‐4 ‐3 ‐2 ‐1 0 1 2 3 4 5 ‐15 ‐10 ‐5 0 5 10 15 Fr in ge s (2 πr ad ) Δθ, (°) Ideal Data Experimental data Parallel Illumination Figure 4.17 ‐ Conical beam phase gradient validation results  66    the illumination and reference beams will change according to the curvature.  Specimen motion  that is not completely parallel to the x‐y plane of the specimen will also change the relative  path lengths across an image and introduce a phase gradient, as will specimen rotation.    To summarize; if the angles of incidence of the two beams is not constant across the specimen  surface and specimen motion occurs, the path lengths of the two beams will not change equally  and a phase gradient will occur.  Specimen rotation will also cause phase gradients for the same  reason.  For nominally flat specimens, adding lenses to create parallel illuminating beams and  adjusting the optical equipment to create constant angles of incidence will allow planar  motions to produce very little relative phase change.  After adding the lenses and adjusting the  vertical mirror near the specimen to create nominally equal angles of incidence, the phase  gradients of various planar motions in the x and y direction after fringe correction are shown in  Figure 4.18.     0.2mm (12.9 pixels)  0.4mm (25.5 pixels) 0.6mm (37.6 pixels) 0.8mm (50.4 pixels)   Figure 4.18 ‐ Phase gradients for various specimen shifts  The top row of figures shows specimen motions in the x direction, and the bottom row shifts in  the y direction.  The relative size of the shifts can be seen by observing the band of information  67    loss visible at the edge of the fringe patterns.   The noise content clearly increases with  specimen shift due to the speckle pattern gradually changing shape.  As the speckles change  shape, some of the speckle correlation is lost upon shifting the images back and the noise is  increased.  Therefore it is expected that fringe correction will only be feasible up to a certain  range which is to be investigated in future experiments.  Comparison of Figure 4.13 with Figure 4.18 shows that the residual phase gradient has been  largely corrected.  A slight phase gradient is still visible with larger specimen motions however,  even with the addition of parallel lenses and beam adjustment.  It should also be noted that  using physical methods to correct phase gradients requires the equipment to be adjusted  independently for each specimen measured.  The fringes in Figure 4.18 were obtained after  careful adjustment to illustrate that this type of gradient correction is possible.  In a more  practical application however it is very likely that phase gradients will still occur; even for the  case shown after careful adjustment slight phase gradients exist.   Therefore an additional  method to systematically correct this artifact would be very beneficial.    To this end, phase gradient correction by mathematical methods was investigated.  A least‐ squares solution, similar to that used to solve for pixel phase (equations (2.17)‐(2.20)) is used to  estimate these phase gradients.  The gradient is considered as a bi‐linear function in the form  of:    P Qx Ry      (4.12)  Considering an area within the image bounded by two concentric circles containing n  pixels:  68      11 1 22 2 1 1 1 n n n x y P x y Q R x y                                       (4.13)  where   ,n nx y are the co‐ordinates of the pixel with the origin equal to that of the circles  bounding the area and  n being the phase of the pixel.    This equation is over determined and can be solved using the least‐square method.  The  polynomial coefficients can be determined as:    2 21 x y P Q R x y            (4.14)  Once the phase gradient is calculated using equation (4.14) it can be subtracted from the fringe  pattern, yielding a more accurate result of the true measurement fringe.  Figure 4.19 (a) shows  a corrected zero fringe with a specimen motion of 0.5mm; the phase gradient is clearly visible.   Figure 4.19 (b) shows the same fringe after the phase gradient is subtracted using the least  squares solution.    (a)  (b) Figure 4.19 ‐ (a) Corrected fringe with phase gradient (b) Corrected  fringe after phase gradient subtraction  69    The fringe after phase subtraction is clearly much cleaner; the gradient is removed leaving a  fringe pattern much closer to the ideal case of a completely white image.  Correcting for the phase gradients in the case of residual stress hole drilling measurements  serves mainly to help avoid unwrapping errors.  The stress calculation itself has provisions for  removing the effect of the gradients; however unwrapping errors can still damage the stress  solution.  By removing the phase gradients, the possibility of extraneous unwrapping errors is  significantly reduced.  For ESPI measurements of other kinds however, mitigating phase  gradients may prove to be more significant.     To summarize, although phase gradients are created using the proposed ESPI fringe correction  process, they can be corrected using a combination of physical and mathematical methods.   These experiments show that after such treatment, a restored zero fringe is shown to have  good agreement with the ideal case of a completely null image.  Noise is clearly visible and will  increase with larger specimen shifts as the speckles change shape, but within a certain range,  the ESPI fringe correction process is shown to be feasible for the case of blank fringes.  4.6 Validation of ESPI Fringe Correction Using Hole Drilling The fringe correction process has been shown to be feasible for correcting a null measurement  fringe which is an important result.  Ideally however this process will be used to correct rigid  body motions affecting actual ESPI measurements.  Therefore, in order to test the effectiveness  of the process in an actual application, the process was applied to actual ESPI hole drilling  measurements.  These series of experiments involved comparing a baseline measurement  performed under controlled conditions with a corrupted measurement containing some  70    specimen rigid‐body motion.  The fringe correction procedure was performed on the corrupted  measurement, and the recovered fringe pattern was used as the basis of the residual stress  solution.  This recovered solution was then compared with the solution obtained using the  baseline measurement to establish the effectiveness of the correction process.  By repeating  this experiment for various specimen motions, a range over which fringe correction is possible  was determined.  The specimen used for these tests was the same calibrated specimen used in section 4.2.  Since  this specimen contained a well‐known residual stress distribution, a measurement which  matches this stress was known to be accurate and serves as the baseline.  So long as the  recovered measurement matches the baseline, its accuracy can be safely assumed.  The initial tests performed examined the maximum range of specimen motion that could be  corrected.  The results for the various specimen directions are shown in Figures 4.20‐4.23.  The  dark black lines represent the baseline measurement performed in section 4.2.  If the residual  stress solution remains close to the line, then the process is successful.   Both positive and  negative Z direction specimen shifts were considered since it was thought the depth of field of  the lens may have an effect on the range of correction.  Examining the results, it’s quite  apparent that the procedure is effective over a specific range, after which the solution  deteriorates rapidly.  If a deviation of ± 3 MPa is considered, the correction is effective over the  ranges shown in Table 4.5.  The standard deviation is also tabulated for each stress component  within this table.     71    ‐20 0 20 40 60 80 100 120 0 0.5 1 1.5 2 2.5 St re ss  (M pa ) Y Direction Specimen Shift (mm) σX σY τxy   Figure 4.20 ‐ Recovered stress solution vs. X direction specimen shift Figure 4.21 ‐ Recovered stress solution vs. Y direction specimen shift  ‐20 0 20 40 60 80 100 120 0 0.5 1 1.5 2 2.5 St re ss  (M pa ) X Direction Specimen Shift (mm) σX σY τxy y xy x y xy x 72      ‐20 0 20 40 60 80 100 120 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 St re ss  (M pa ) +Z Direction Specimen Shift (mm) σX σY τxy ‐20 0 20 40 60 80 100 120 0 0.5 1 1.5 2 St re ss  (M pa ) ‐Z Direction Specimen Shift (mm) σX σY τxy Figure 4.22 – Recovered stress solution vs. +Z direction specimen shift Figure 4.23 ‐ Recovered stress solution vs. ‐Z direction specimen shift y x x y x x 73    Table 4.5 ‐ Range and deviation of fringe correction process  Specimen Shift Direction Range (mm) Standard Deviation (Mpa)  x   y   xy   X  1.4mm  0.283 0.642 0.243  Y  1.2mm  1.086 0.928 0.356  +Z  0.8mm  0.516 1.065 0.186  ‐Z  0.9mm  0.445 1.230 0.339    The data in the table show that specimen motions in the out‐of‐plane (Z) direction have the  most detrimental effect on the measurement recovery, as the range of correction in these  directions were both significantly lower than the in plane (X,Y) directions.  This follows from the  fact that specimen motions in the z direction change the path length of both ESPI beams much  more than in‐plane motions.  As the pattern changes the re‐correlation process suffers and the  measurement is lost.  Another trend is observed examining the standard deviations of the  various in plane stress components.  The y component of stress exhibits the highest deviation  across nearly all specimen motions.  This may be attributed to the fact that the optical setup of  the ESPI system is aligned to measure displacements in the x direction only.  For y‐stresses,  such displacements occur only through the action of Poisson’s ratio, and are only about one  third the size of the corresponding displacements for x‐stresses.  Thus, measurement of x‐ displacements gives sensitivity to y‐stresses of about one third the sensitivity to x‐stresses.  This  trend is evident when considering the graphs.  The y component of stress consistently deviates  from the true solution first while the x component of stress, which is more directly measured,  maintains agreement some distance longer.  If the system was adapted to measure all three  components of stress evenly, by changing the sensitivity vector such that x and y displacements  74    are measured equally, the range of correction for this application could most likely be  extended.  4.7 Specimen Removal and Replacement As well as increasing the resilience of ESPI measurements to specimen rigid body motions, the  ESPI fringe correction process developed here has now made possible a novel measurement  method.  Throughout previous experiments the drilling process was integrated within the  apparatus – the drill was moved into position after the reference set of images was recorded.   This can cause problems however as chips created during drilling can contaminate the optics  and the cutting forces themselves can disturb the optical equipment.  If drilling could be  performed at a separate location these problems could be avoided, and the apparatus itself  could be simplified.  This provided the motivation for the following experiments.   The specimen used for these tests was developed for previous work [40] for hole drilling  residual stress measurements.  To establish a baseline residual stress result, a series of 5  measurements was performed with the specimen fixed rigidly to the optical table.  The results  are shown in Figure 4.24, with the average for each in plane stress component plotted as a  straight line.     75      Figure 4.24 ‐ Removal & replacement baseline tests with the specimen fixed    Examining these results, they are quite closely grouped around the averages indicating  reasonable consistency.  This average stress solution therefore represents a reasonable  baseline residual stress solution for this particular specimen.  If a measurement produced using  specimen removal and replacement can produce a similar solution it can be successfully  validated.  ‐20 0 20 40 60 80 100 120 140 160 1 2 3 4 5 St re ss  (M Pa ) Measurement # σX σY τxy y xy x  To accom hold the  process o existing E with the  replaced repeatab additiona the expe straight l does sho with all t between experime Figure  plish specim specimen in f these exp SPI system. same drillin  in the fixtur ility of the r l 10 times,  riment are s ine.  Consid w promise.  he measure  the baselin nts is show 4.25 ‐ Speci en remova  the path of eriments co   The fixture g system us e and the fr eplacement with a new f hown in Fig ering the res  The remova ments space e tests wher n in Table 4 men fixture l and replac  the ESPI be nsisted of a  was then b ed in previo inge set of i  procedure  ringe set of ure 4.26, wi ults, there  l and replac d closely ar e the specim .6.  ement, an a ams.  This f ligning the s olted to the us experime mages was  the specime  images bein th the avera is definite va ement proc ound the av en was fixe dditional fix ixture is sho pecimen in   optical tabl next m the se record remov bolted differen nts.  After d recorded.   I n was then  g recorded ge of all 10  riation; how ess is show erage.  A nu d, and the  ture was co wn in Figure the fixture w e.  The hole arked on th t of referenc ed.  The spe ed from the  to a separa t location a rilling, the s n order to t removed an  each time.   tests again s ever the m n to have go merical com removal and nstructed to  4.25.  The  ithin the   location wa e specimen e images w cimen was   fixture and te fixture at nd then dril pecimen wa est the  d replaced  The results  hown as a  easurement od repeatab parison   replaceme 76    s   and  as  then     a  led  s  an  of    ility  nt  77    Figure 4.26 ‐ Removal and replacement residual stress results     Table 4.6 ‐ Removal and replacement vs. baseline numerical comparison                Comparing both tests numerically shows that they are in good agreement, with the stress  component in the y direction varying more than the other two.  Again this can be attributed to  the low sensitivity to displacements in the y direction.  This difference may also be attributed  due to the need to reposition the specimen and fixture independently of the previous tests.   ‐20 0 20 40 60 80 100 120 140 1 2 3 4 5 6 7 8 9 10 St re ss  (M Pa ) Measurement # σX σY τxy   Average (MPa)   x   y xy Baseline  129.4 92.4 0.4 Removal & Replacement 123.5 78.2 ‐6.4 Absolute Difference  5.9 14.2 6.9 y xy x 78    Although care was taken to try and return the specimen to the same position used in the fixed  tests, small deviations could cause a slightly different measurement due to the beam angles of  the illumination and reference beam changing.  Since the removal and replacement results are  consistently lower than those of the fixed tests, this could certainly be the case.  A graphical  comparison of the fringes obtained during both the baseline and removal and replacement  tests are shown in Figure 4.27  This figure clearly shows the effectiveness of the fringe correction process as well as good  agreement between the baseline fixed specimen measurements and the removal and  replacement measurements. The effectiveness of the correction process is clearly evident when  considering Figure 4.27 (b) and (c). (b) shows a completely corrupted fringe pattern, while after  correction (c) shows a fringe pattern with good contrast and well defined fringes.  Comparing  Figure 4.27 (a) and (c) the similarity is evident; the fringe patterns show very close agreement,  indicating that a comparable solution can indeed be obtained using specimen removal and  replacement, even using the simple fixture that was constructed.  The noise content is slightly  higher in the corrected fringe however this is typical of the correction process as previously  discussed and does not affect the measurement significantly.    Figure 4.27 ‐ (a) Baseline fringe (b) Removal and replacement fringe before correction         (c) Removal and replacement fringe after correction  79    The experimental results presented in this section work to investigate and validate the various  aspects of the ESPI fringe correction process.  In order generate accurate ESPI data, an  adjustable ESPI hole drilling apparatus was constructed and validated using a calibrated  specimen.  In order to make use of phase stepped ESPI images for DIC measurements (EDIC), a  phase independent method was successfully developed.  The size effect of the speckles within  these ESPI images was investigated and an optimum speckle size of 1.73px, which yielded the  largest EDIC measurement range, was determined for this optical setup.  Null fringe correction,  representing the most basic application of the fringe correction process, was investigated and  shown to be successful; however a phase gradient artifact was encountered.  The cause of  these phase gradients was determined analytically, validated experimentally and corrected  using a combination of both physical and mathematical methods.  Once the fringe correction  process was established, it was applied to ESPI hole drilling measurements.  Using the 3 axis  table, specimen motions were systematically created within hole drilling measurements.  By  comparing recovered solutions after fringe correction with a baseline measurement performed  in controlled conditions, the accuracy of the process could be assessed.   This allowed the  effective range of the correction process to be determined.  Finally, a new ESPI measurement  method involving removal and replacement of the specimen was investigated.  This type of  measurement was made possible due to the flexibility afforded by the fringe correction  process.  It was shown that this type of new measurement is accurate and repeatable when  compared with a baseline measurement performed in the conventional manner.    80    5 Conclusions and Future Work 5.1 Conclusions A novel adaptable ESPI system was designed and constructed  The ESPI apparatus developed here was designed and constructed for use when measuring  residual stresses using the hole drilling method.  The apparatus was validated using a specimen  with a known residual stress field; results were obtained and shown to be in good agreement  with the stresses contained in the specimen.  Phase independent DIC using ESPI phase stepped image sets (EDIC) was developed  By constructing a phase independent image from the available images within a stepped set,  accurate DIC results could be obtained using the existing images within an ESPI measurement.   Statistical tests determined that using an average (“A”) image yielded the best results of any  possible combination of phase‐stepped images or other phase independent quantities.    Using this new EDIC method, the effect of laser speckle size was investigated systematically by  varying the optical configuration of the apparatus.  It was determined that speckle size did have  a significant effect on EDIC effectiveness, and an optimum lens aperture F16 was determined to  yield the most accurate results for this particular apparatus.  This agrees well with previous  work done on single beam DSC measurements.  When compared with conventional painted  speckles, EDIC was found to be less accurate, but the errors within both methods were  comparable.  The difference in accuracy can be attributed to the spatial dependency of laser  speckles.  Physical surface features do not suffer from this dependency and therefore generally  81    yield more accurate DIC results.  Using EDIC, displacements in excess of 150 pixels or 2.5mm  were measured successfully with sub pixel error.  An ESPI fringe correction process was developed to correct for specimen motion  Utilizing EDIC measurements, an ESPI fringe correction process was developed and successfully  applied to residual stress hole drilling measurements.  Significant specimen rigid body motions  occurring between the recordings of phase stepped sets were able to be corrected using this  method.    Phase gradients resulting from this fringe correction were investigated and were found to be  mostly caused by the use of conical illumination beams.  This mechanism was shown  analytically and experimentally.  The phase gradients were largely mitigated by the addition of  lenses to convert the conical illumination into parallel illumination. In addition, the optical  equipment was carefully adjusted to ensure equal angles of incidence for both the reference  and illumination beams.  Even with careful optical adjustments, some gradients remained due  to specimen rotations and curvature.  To address these remaining gradients, mathematical  correction by a least squares solution and subtraction method was developed.  Using a precision 3‐axis table, specimen motions were inserted into actual ESPI residual stress  measurements after the hole drilling had occurred.  Using the fringe correction process,  measurements could be recovered accurately for specimen motions of up to 1.4mm in the x  direction, 1.1mm in the y direction and 0.8mm and 0.9mm in the positive and negative z  directions respectively.  These results were validated by matching them to conventional ESPI  measurements performed on the same specimen.  Without the fringe correction process,  82    motions of approximately 0.030mm in any direction would compromise the measurement  entirely.  The fringe correction process is therefore shown to improve the ability of ESPI to  measure stress accurately, independent of these motions, by 2 orders of magnitude.  This result is encouraging with respect to helping ESPI move to more practical use in the field,  especially in terms of a portable ESPI system.  Within such an application, it would be very  possible for the specimen to shift after drilling, due to either environmental disturbances or the  drilling itself.  Using the fringe correction process developed, no additional equipment or data  are required, and the ability of the system to handle specimen motions is significantly  improved.  Specimen removal and replacement was investigated  A novel residual stress hole drilling technique was also made possible by the fringe correction  process.  Previously the drilling process had been integrated within the experimental apparatus;  however measurements where the specimen was removed for drilling and replaced within the  system were investigated.  A simple fixture was developed for this purpose, and a new  specimen was fabricated.  This specimen was then measured 5 times using the conventional  measurement method with the specimen fixed, and a baseline residual stress solution was  established by taking the average of these tests.  The specimen and fixture were then aligned  within the system and the new process was tested.  In this case the specimen was removed  after the reference set of images was recorded, drilled at a separate location with the same  drilling system used previously, and replaced within the fixture.  The fringe set of images was  then recorded and the fringe correction procedure was applied.  Using this method, results  83    were obtained that were comparable to those done in the baseline tests with the specimen  fixed.  In order to test the repeatability of the method, the specimen was removed and  replaced an additional 9 times and a new set of fringe images was recorded.  The results for all  10 tests were shown to be in good agreement indicating that the process possesses good  repeatability.  This is an exciting result as it creates several new possibilities.  The drilling system can now  easily be separated from the ESPI apparatus.  By doing this, the optical system itself can be  simplified, and components can be placed in better alignment.  In addition, the chips created  during the cutting process can be kept separate from the delicate optical equipment.  Even with  the chip suction system implemented, specimen chips were still a problem within the apparatus  and they would often obstruct the vertical mirror in particular.  This problem can now be  avoided.  In addition, the data corruption caused by the cutting forces induced during the  drilling process can also be avoided using specimen removal and replacement.  Specimen  removal and replacement also opens the possibility of measuring a single specimen at multiple  attitudes within a single system.  For instance, consider a hole drilling measurement performed  on a thin plate.  Two measurements can now be performed simultaneously on the front and  back faces of the plate by simply taking a reference set of images with the specimen in the  fixture, flipping the specimen and recording reference images of the back face.  By repeating  this process after the hole is drilled, displacement data for the front and back face can be  obtained simultaneously.  Previously it would be necessary to devise an extremely complicated  ESPI system to achieve this type of measurement however through this work it is now possible  using the conventional ESPI system constructed.  84    5.2 Future Work Investigate correction of individual stepped images for active correction  This work has been shown to correct specimen motions occurring in between recording sets of  phase‐stepped images.  Vibrations may cause the specimen to shift during the acquisition of a  phase‐stepped set, causing the steps to deviate from the desired 90‐degree increments.  Some  work was done examining the feasibility of correcting individual images within a set, however  the tests done were mainly inconclusive because the induced vibrations were too fast to be  accurately captured by the camera; significant blurring occurred which destroyed the speckle  correlation.  The use of stroboscopic laser illumination in combination with a higher speed  camera could reduce this blurring.  If this is possible, then individual images within a set could  possibly be corrected in the same manner as full sets within this work, so long as the imaging  frequency was high enough to minimize the displacement between each image in a set.  This  could allow another significant increase in ESPI system robustness and could mitigate the effect  of vibrations and other disturbances common to field environments.  Extend fringe correction to in plane rotations  Within this work, a fringe correction process is developed that can account for planar motions  and out of plane rotations.  It is possible however that in plane (xy) specimen rotations could  occur as well.  Therefore, it follows logically that adapting the process to account for these  rotations as well as investigating the performance of rotation correction should be investigated  in the future.    85    Develop a portable ESPI system  The ESPI fringe correction process represents a solid step towards constructing a viable ESPI  system for portable measurements.  Using this method as well as the adjustable design  developed for the experimental apparatus, it become practical to develop a feasible system  that could be used for making mobile residual stress measurements.  If work was done to  develop single image correction for more active specimen motion correction the robustness of  a portable system could be even further increased.    86    References [1] J. De Lemos, H. Mischo, T. Pfeifer, Comparison of Strain/Stress Measurements on Free Form  Surfaces Using ESPI and Strain Gauge Technique, Proceedings of SPIE ‐ The International Society  for Optical Engineering, v 4420, Laser Metrology for Precision Measurement and Inspection in  Industry (2001) 123‐131.   [2] M.V. Rao, R. Samuel, A. Ananthan, S. Dasgupta, P.S. Nair, Electronics Speckle Interferometry  Applications for NDE of Spacecraft Structural Components, Proceedings of SPIE,  v 7155, Ninth  International Symposium on Laser Metrology (2008) 563‐571.  [3] F. Jin, F.P. Chiang, ESPI and Digital Speckle Correlation Applied to Inspection of Crevice  Corrosion on Aging Aircraft, Research in Nondestructive Evaluation v 10 (1998) 63‐73.  [4] G.S. Schajer, M. Steinzig, Sawblade Vibration Mode Shape Measurement Using ESPI, Journal  of Testing and Evaluation, v 36 (2008) 259‐263.  [5] J.M. Sabatier, V. Aranchuk, W.C.K. Alberts II, Rapid High‐Spatial‐Resolution Imaging of Buried  Landmines Using ESPI, Proceedings of SPIE, v 5415, Detection and Remediation Technologies  for Mines and Minelike Targets IX, (2004) 14‐20.   [6] G.S. Spagnolo, D. Ambrosini, D. Paoletti, ESPI for Mosaics Diagnostics, Proceedings of SPIE. v  4778, Interferometry XI: Applications (2002) 377‐384.  [7] M. Steinzig, E. Ponslet, Residual Stress Measurement Using the Hole Drilling Method and  Laser Speckle Interferometry: Part 1, Experimental Techniques, v 27, (2003) 43‐46.   [8] M.A. Sutton, Digital Image Correlation for Shape and Deformation Measurements, Springer  Handbook of Experimental Solid Mechanics, (2008).    [9] D. Amodio, G. Broggiato, F. Campana, G. Newaz, Digital Speckle Correlation for Strain  Measurement by Image Analysis, Experimental Mechanics, v 43, (2003) 396‐402.   [10] I. Yamaguchi, Measurement and Testing by Digital Speckle Correlation, Proceedings of SPIE  v 7129, Seventh International Symposium on Instrumentation and Control Technology:  Optoelectronic Technology and Instruments, Control Theory and Automation, and Space  Exploration, (2008) 71290Z‐9.   [11] L. Muravs'kyi, M. Hvozdyuk, T. Polovynko, Evaluation of The Surface Strains in the  Composites by the Methods of Digital Speckle Correlation, Materials Science, v 43 (2007) 568‐ 573.   87    [12] G. Schajer, M. Steinzig, Full‐Field Calculation of Hole Drilling Residual Stresses From  Electronic Speckle Pattern Interferometry Data, Experimental Mechanics, v 45, (2005) 526‐532.   [13] R. Ritter, K. Galanulis, D. Winter, E. Mueller, B. Breuckmann, Notes on the Application Of  Electronic Speckle Pattern Interferometry, Optics and Lasers in Engineering. v 26, (1997) 283‐ 299.   [14] F.S. Chau, J. Zhou, Combined Digital Speckle Shearing Interferometry and Digital Image  Correlation for Analysing Vibrating Objects, Proceedings of SPIE ‐ The International Society for  Optical Engineering, 2nd International Conference on Experimental Mechanics, v 4317, (2001)  380‐385.  [15] D. Findeis, J. Gryzagoridis, D.R. Rowland, Vibration Isolation Techniques Suitable for  Portable Electronic Speckle Pattern Interferometry, Proceedings of SPIE ‐ The International  Society for Optical Engineering, Nondestructive Evaluation and Health Monitoring of Aerospace  Materials and Civil Infrastructures, v 4704, (2002) 159‐167.   [16] P.L. Reu, B.D. Hansche, Digital Image Correlation Combined With Electronic Speckle Pattern  Interferometery for 3D Deformation Measurement In Small Samples, SEM Annual Conference  and Exposition on Experimental and Applied Mechanics 2006, v 3, (2006) 1396‐1402.  [17] W.H. Peters, W.F. Ranson, Digital Imaging Techniques in Experimental Stress Analysis,  Optical Engineering, v 21 (1982) 427‐431.   [18] M. Sutton, W. Wolters, W. Peters, W. Ranson, S. McNeill, Determination of Displacements  Using an Improved Digital Correlation Method, Image Vision Computing, v 1 (1983) 133‐139.   [19] T. Chu, W. Ranson, M. Sutton, Applications of Digital‐Image‐Correlation Techniques to  Experimental Mechanics, Experimental Mechanics, v 25, (1985) 232‐244.   [20] N. Takai, T. Asakura, Vectorial Measurements of Speckle Displacement by the 2‐D  Electronic Correlation Method, Applied Optics, v 24, (1985) 660‐665.   [21] D.J. Chen, F.P. Chiang, Optimal Sampling And Range of Measurement in Displacement‐Only  Laser‐Speckle Correlation, Experimental Mechanics, v 32, (1992) 145‐153.   [22] M. Sjodahl, L.R. Benckert, Systematic and Random Errors in Electronic Speckle  Photography, Applied Optics, v 33, (1994) 7461‐7471.   [23] P. Bing, Performance of Sub‐Pixel Registration Algorithms in Digital Image Correlation,  Measurement Science and Technology, v 17, (2006) 1615.   [24] R. Feiel, P. Wilksch, High‐Resolution Laser Speckle Correlation for Displacement and Strain  Measurement, Applied Optics, v 39, (2000) 54‐60.   88    [25] M. Sjodahl, L. Larsson, Speckle Correlation Used for Measuring Micro‐Structural Changes in  Paper, Proceedings of SPIE, Speckle Metrology 2003, v 4933, (2003) 33‐38.  [26] G. Cloud, Optical Methods in Experimental Mechanics Part 26: Subjective Speckle,  Experimental Techniques, v 31, (2007) 17‐19.   [27] G. Cloud, Optical Methods in Experimental Mechanics Part 27: Speckle Size Estimates,  Experimental Techniques, v 31, (2007) 19‐22.  [28] D.M. Kennedy, Z. Schauperl, S. Greene, Application of ESPI‐Method for Strain Analysis in  Thin Wall Cylinders, Optics and Lasers in Engineering, v 41, (2004) 585‐594.   [29] J.C. Wyant, Phase‐Shifting Interferometry, Optics 513 Course Notes, (1998), pp. 1‐38.  Available at:  http://www.optics.arizona.edu/jcwyant/Optics513/ChapterNotes/Chapter05/3.PrintedVersion PhaseShiftingInterferometry.pdf .  Retrieved on: July 25th, 2010  [30] Y. An, Residual Stress Measurement Using Cross‐Slitting and ESPI, University of British  Columbia, MASc Thesis, (2008).   [31] R.S. Sirohi, Optical Methods of Measurement : Wholefield Techniques, New York : Marcel  Dekker, (1999).  [32] H. Kadono, H. Takei, S. Toyooka, A Noise‐Immune Method of Phase Unwrapping in Speckle  Interferometry, Optics and Lasers in Engineering, v 26 (1997) 151‐164.    [33] R.C. Gonzalez, Digital Image Processing, Reading, Mass. : Addison‐Wesley, (1992).   [34] E.L. Hall, A Survey of Preprocessing and Feature Extraction Techniques for Radiographic  Images, IEEE Transactions on Computers, v 20, (1971) 1032‐1044.   [35] D.V. Nelson, A. Makino, T. Schmidt, Residual Stress Determination Using Hole Drilling and  3D Image Correlation, Experimental Mechanics, v 46, (2006) 31‐38.   [36] J.D. Lord, D. Penn, P. Whitehead, The Application of Digital Image Correlation for  Measuring Residual Stress by Incremental Hole Drilling, Applied Mechanics and Materials, v 13‐ 14 (2008) 65‐73.    [37] D. Lecompte, A. Smits, S. Bossuyt, H. Sol, J. Vantomme, D. Van Hemelrijck, A.M. Habraken,  Quality Assessment of Speckle Patterns for Digital Image Correlation, Optics and Lasers in  Engineering, v 44, (2006) 1132‐1145.   89    [38] G. Cloud, Optical Methods in Experimental Mechanics Part 24: Demonstrations of Laser  Speckle Phenomena, Experimental Techniques, v 30 (2006) 27‐30.   [39] B. Winiarski, P.J. Withers, Micron‐scale Residual Stress Measurement using Micro‐hole  Drilling and Digital Image Correlation, SEM Annual Conference and Exposition on Experimental  and Applied Mechanics 2010, (2010).  [40] M.C. Cavusoglu, Improvements In Electronic Speckle Pattern Interferometry for Residual  Stress Measurements, MASc Thesis, University of British Columbia, (2007).   

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.24.1-0071157/manifest

Comment

Related Items