Advanced Stamping Manufacturing Engineering, Auburn Hills, MI

Similar documents
Auburn Hills, MI, USA, 2 Dept. of Mechanical Engineering, Oakland University, Rochester, MI, USA

Enhanced two-frequency phase-shifting method

An Innovative Three-dimensional Profilometer for Surface Profile Measurement Using Digital Fringe Projection and Phase Shifting

SIMULATION AND VISUALIZATION IN THE EDUCATION OF COHERENT OPTICS

Applications of Temporal Phase Shift Shearography

INFLUENCE OF CURVATURE ILLUMINATION WAVEFRONT IN QUANTITATIVE SHEAROGRAPHY NDT MEASUREMENT

Laser sensors. Transmitter. Receiver. Basilio Bona ROBOTICA 03CFIOR

Range Sensors (time of flight) (1)

Study on Gear Chamfering Method based on Vision Measurement

Edge and local feature detection - 2. Importance of edge detection in computer vision

Full surface strain measurement using shearography

Complex Sensors: Cameras, Visual Sensing. The Robotics Primer (Ch. 9) ECE 497: Introduction to Mobile Robotics -Visual Sensors

Chapter 37. Wave Optics

Measurements using three-dimensional product imaging

Draft SPOTS Standard Part III (7)

Chapter 3 Image Registration. Chapter 3 Image Registration

Stereo imaging ideal geometry

Depth. Common Classification Tasks. Example: AlexNet. Another Example: Inception. Another Example: Inception. Depth

Method for improving sinusoidal quality of error diffusion binary encoded fringe used in phase measurement profilometry

A RADIAL WHITE LIGHT INTERFEROMETER FOR MEASUREMENT OF CYLINDRICAL GEOMETRIES

17th World Conference on Nondestructive Testing, Oct 2008, Shanghai, China

Structured Light. Tobias Nöll Thanks to Marc Pollefeys, David Nister and David Lowe

A Comparison of Results Obtained from Alternative Digital Shearography and ESPI Inspection Methods

10/5/09 1. d = 2. Range Sensors (time of flight) (2) Ultrasonic Sensor (time of flight, sound) (1) Ultrasonic Sensor (time of flight, sound) (2) 4.1.

Model-Based Segmentation of Impression Marks

Testing spherical surfaces: a fast, quasi-absolute technique

AP* Optics Free Response Questions

Full-field optical methods for mechanical engineering: essential concepts to find one way

A Method of weld Edge Extraction in the X-ray Linear Diode Arrays. Real-time imaging

A Survey of Light Source Detection Methods

STUDY ON LASER SPECKLE CORRELATION METHOD APPLIED IN TRIANGULATION DISPLACEMENT MEASUREMENT

QUANTITATIVE ANALYSIS OF A CLASS OF SUBSURFACE CRACKS USING

ksa MOS Ultra-Scan Performance Test Data

Minimizing Noise and Bias in 3D DIC. Correlated Solutions, Inc.

Subpixel Corner Detection Using Spatial Moment 1)

specular diffuse reflection.

Coherent digital demodulation of single-camera N-projections for 3D-object shape measurement: Co-phased profilometry

OPTICAL TOOL FOR IMPACT DAMAGE CHARACTERIZATION ON AIRCRAFT FUSELAGE

3D Computer Vision. Structured Light I. Prof. Didier Stricker. Kaiserlautern University.

Distortion Correction for Conical Multiplex Holography Using Direct Object-Image Relationship

Engineered Diffusers Intensity vs Irradiance

Invited Paper. Katherine Creath and James C. Wyant WYKO Corporation 2650 East Elvira Road, Tucson, Arizona ABSTRACT 1. INTRODUCTION.

Chapter 37. Interference of Light Waves

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

IMGS Solution Set #9

Projector Calibration for Pattern Projection Systems

Extracting Sound Information from High-speed Video Using Three-dimensional Shape Measurement Method

Chapter 12 Notes: Optics

Multi-component shearography using optical fibre imaging-bundles

UNIT 102-9: INTERFERENCE AND DIFFRACTION

Digital Shearographic Non-destructive Testing Laboratory

Traditional pipeline for shape capture. Active range scanning. New pipeline for shape capture. Draw contours on model. Digitize lines with touch probe

Dynamic three-dimensional sensing for specular surface with monoscopic fringe reflectometry

A Novel Double Triangulation 3D Camera Design

A three-step system calibration procedure with error compensation for 3D shape measurement

Image Formation. Antonino Furnari. Image Processing Lab Dipartimento di Matematica e Informatica Università degli Studi di Catania

REMOTE SENSING OF SURFACE STRUCTURES

THE EXPERIMENTAL STUDY OF MATERIAL SEPARATION IN PERSONAL TYRES

To see how a sharp edge or an aperture affect light. To analyze single-slit diffraction and calculate the intensity of the light

Orthogonal Projection Matrices. Angel and Shreiner: Interactive Computer Graphics 7E Addison-Wesley 2015

Addressing High Precision Automated Optical Inspection Challenges with Unique 3D Technology Solution

Three-Dimensional Measurement of Objects in Liquid with an Unknown Refractive Index Using Fisheye Stereo Camera

Graphics and Interaction Transformation geometry and homogeneous coordinates

Eric Lindmark, Ph.D.

D&S Technical Note 09-2 D&S A Proposed Correction to Reflectance Measurements of Profiled Surfaces. Introduction

Multi-view Surface Inspection Using a Rotating Table

Camera Model and Calibration

ENGR142 PHYS 115 Geometrical Optics and Lenses

A Fringe Center Detection Technique Based on a Sub-Pixel Resolution, and Its Applications Using Sinusoidal Gratings

Optic Flow and Basics Towards Horn-Schunck 1

A 100Hz Real-time Sensing System of Textured Range Images

Chapter 36. Image Formation

Comparative study on passive and active projector nonlinear gamma calibration

Chapter 24. Wave Optics

Shearography as an Industrial Application Including 3D Result Mapping

THE VIEWING TRANSFORMATION

Optik. Analia Sicardi-Segade a,, Amalia Martínez-García a, Noel-Ivan Toto-Arellano b, J.A. Rayas a. a r t i c l e i n f o.

Advanced Vision Guided Robotics. David Bruce Engineering Manager FANUC America Corporation

Akrometrix Testing Applications

Reflection and Refraction of Light

Inspection of Laser Generated Lamb Waves using Shearographic Interferometry

1. (25pts) Answer the following questions. Justify your answers. (Use the space provided below and the next page)

ACTA TECHNICA CORVINIENSIS Bulletin of Engineering Tome X [2017] Fascicule 2 [April June] ISSN:

1 (5 max) 2 (10 max) 3 (20 max) 4 (30 max) 5 (10 max) 6 (15 extra max) total (75 max + 15 extra)

A New Technique of Extraction of Edge Detection Using Digital Image Processing

Comparative measurement of in plane strain by shearography and electronic speckle pattern interferometry

Optical 3D Sensors for Real Applications Potentials and Limits Applications pratiques de capteurs optiques tridimensionnels: potentiel et limites

New Opportunities for 3D SPI

3D data merging using Holoimage

16. Holography. Dennis Gabor (1947) Nobel Prize in Physics (1971)

ACTIVITY TWO CONSTANT VELOCITY IN TWO DIRECTIONS

Polarization of Light

Flatness Measurement of a Moving Object Using Shadow Moiré Technique with Phase Shifting

Coherent Gradient Sensing Microscopy: Microinterferometric Technique. for Quantitative Cell Detection

Seam tracking for fillet welds with scanner optics

Stevens High School AP Physics II Work for Not-school

Contour LS-K Optical Surface Profiler

Introduction to Computer Vision. Introduction CMPSCI 591A/691A CMPSCI 570/670. Image Formation

Massachusetts Institute of Technology. Department of Computer Science and Electrical Engineering /6.866 Machine Vision Quiz I

Digital Volume Correlation for Materials Characterization

Transcription:

RECENT DEVELOPMENT FOR SURFACE DISTORTION MEASUREMENT L.X. Yang 1, C.Q. Du 2 and F. L. Cheng 2 1 Dep. of Mechanical Engineering, Oakland University, Rochester, MI 2 DaimlerChrysler Corporation, Advanced Stamping Manufacturing Engineering, Auburn Hills, MI Abstract: Small distortions on the surface of sheet metal stamping parts may have significant impact on the optical appearance of the final product after finish treatment. Recently, there are considerable research activities in the development of optical methods for surface distortion measurement because of the fullfield and nondestructive nature of these methods. Optical methods such as Lights and Mirrors, Defracto, 3D-Digitazition for contour measurement, 3D-Digitazition for surface distortion measurement etc. are emerging as strong candidates for industrial inspections of surface distortion. Because each technique has always its unique advantages and disadvantage for different applications, the selection of one technique depends on the particular objective as well as limitations of each technique. This paper will give a through review of the above optical methods for surface distortion measurement. In particular, the recent development 3D-Digitazition for surface distortion measurement will be discussed in details. The fundamental and methodology are demonstrated by example of surface distortion measurement on door panel of a car body. Introduction: The demands for greater quality on the surface sheet metal stamping parts has created a need for better techniques of nondestructive and not-contacting surface inspection. Nowadays, different technologies have been used for body surface inspection. They can be summarized as follows: Lights and mirrors, D Sight from Defracto and 3D digitization for measuring surface shape. This method of Lights and Mirror has been widely used to inspect surface distortion by observing (through human eyes) the reflected light from mirror-like surface. If the surface to be inspected is not smooth enough, a thin oil film should be spread on the surface so that the applied light can be reflected from the surface. The measuring sensitive of this technique greatly depends on the observation direction, the bigger the angle is, the higher the sensitivity. The advantage of this method is its simplification. It is able to find a small surface distortion, such as a dent or a surface low. However, it is a qualitative method and impossible to indicate the depth of the lows, a success entirely depends on observing and illuminating direction and the experiences of the operator. Moreover, the method is suited only for a mirror-like surface, for the surfaces like aluminum and sheet-metal material, a surface pre-processing (e.g. applying a oil on the surface) are required. The technique of Defracto technique has also been used for detecting a surface distortion by automotive industries. Although the measuring sensitivity of this technique is very high (it has a capability to find a surface distortion less than 25 microns in depth), the results obtained depend greatly on the experiences of the operator. Similar to Light and mirror method, the measuring sensitive greatly depends on the observing and illuminating direction, and it is a qualitative method and suited only for a mirror-like surface Nowadays, the demands for greater quality on the surface sheet metal stamping parts has created a need for a technique of quantitative measurement of surface distortion. The mission of the measurement is not only to find the surface distortion on sheet metal stamping parts, but also to improve the manufacturing quality, e.g. to improve the quality of the die-sets. 3D digitization for measuring surface shape is a relatively new technique to deliver surface geometry quantitatively and has been applied for contour and 3D-shape measurement. In this paper, we have developed this technique for surface distortions measurement so that the surface low, highs and dents can be evaluated quantitatively. Fundamental of 3D Digitization: 3D digitization is a technique for 3D shape (contour) measurement based on the optical triangulation. Fig. 1 shows the fundamental of spot of light triangulation. A relationship between a displacement x (usually called distortion) and the surface depth change z can be derived as following:

x = M a (1) where x is the distortion, M is the magnification of the camera used and a is a distance from observed point on the object surface to optical axis of the digital camera. Because a is equal to z sin θ from geometry, the relationship between the distortion x and the surface depth change z is given by: x = M sin θ z = k z (2) z Optical light θ a x Lens of a camera Detector where k is called the sensitivity factor. The fundamental of 3D digitization is to get the contour information z by measuring the fringe distortion x. Fig. 1 shows only the spot of light Fig. 1 Principle of Spot of light triangulation triangulation. In order to measure the contour information on a whole surface, a fringe pattern created by a projector should be projected onto the object surface. Fig. 2 shows an example of 3D digitization for whole field 3D shape measurement. While a fringe pattern is projected onto a flat surface, the fringes are a pattern that looks like a group of parallel lines (cf. Fig. 2a). While the fringe pattern is projected onto a 3D-object surface, the fringe pattern is distorted and it looks like Fig. 2b. By measuring the magnitude and direction of the distortions at each position, the 3D-shape of object to be tested can be determined in full field. Fig. 2 a fringe pattern projected on a flat surface, the fringe pattern projected on 3D-object A key issue to determine the 3D-shape is how to determine the distortion x as showed in Fig. (2b). In the following part of the paper, we will discussed how to quantitatively determine the distortion x. An equation of the intensity expression for the fringe pattern on the flat object surface can be explained by 1-2 : I = a(x,y) + b(x,y) cos ( P 2π x ) or I = a(x,y) + b(x,y) cos f (3) where f = ( P 2π x) is called a phase distribution of a intensity equation, P is the pitch of the fringe pattern, i.e. the distance between two adjacent fringes, a(x,y) is the background of the fringe pattern and b(x,y) is

related to the contrast of the pattern. While the fringe pattern is projected onto a 3D-object surface, the fringes are distorted and the intensity equation on the 3D-object surface becomes: I = a(x,y) + b(x,y) cos [ P 2π (x x)] = a(x,y) + b(x,y) cos { P 2π [x k(x,y) z(x,y)]} or I = a(x,y) + b(x,y) cos c (4) where c = P 2π [x k(x,y) z(x,y)] is called the phase distribution on 3D object, k(x,y) is the sensitivity factor at each point which has been explained in equation (2), z(x,y) is the contour information which we want to measure. Assuming c and f is determinable, the subtraction between c and f produces: f - c = { P 2π x - P 2π [x k(x,y) z(x,y)]} = P 2π k(x,y) z(x,y) (5) Equation 5 can be rewritten by: z(x,y) = ( f - c )/ P 2π k(x,y) = ( f - c )/K(x,y) (6) where z(x,y) is the contour data, f and c are the phase distribution of the intensity equation shown in (3) and (4), respectively. The capital letter K(x,y) is the sensitivity factor for a whole-field measuring system, it is direct proportional to k(x,y) - the sensitivity factor for point measurement, but indirect proportional to the pitch P. Now, the key for determine contour z(x,y) becomes: (1) To quantitatively measure f and c (2) To determine K(x,y), also called calibration. In order to determine the phase distribution f or c from an intensity equation, the phase shift technique should be used. The Phase shift technique is method to quantitatively determine phase values of f or c from three or four intensity images (equations) 3. The following shows briefly the fundamental for calculating phase distribution f by recording four intensity images: I 1 (x,y) = a(x,y) + b(x,y) cos ( f ) I 2 (x,y) = a(x,y) + b(x,y) cos ( f +90 ) = a(x,y) - b(x,y) sin f I 3 (x,y) = a(x,y) + b(x,y) cos ( f +180 ) = a(x,y) - b(x,y) cos f (7) I 4 (x,y) = a(x,y) + b(x,y) cos ( f +270 ) = a(x,y) + b(x,y) sin f The phase shifting 90, 180 and 270 can be obtained by moving the fringe pattern to small displacements P/4, P/2 and 3P/4, respectively. The phase distribution f can be determined now from the above equations as follows: I 4( x, y) - I 2( x, y) f = arc tan [ ] (8) I ( x, y) - I ( x, y) 1 3 In same way, the phase distribution c in the intensity equation I can be determined too.

To determine (calibrate) the system factor K(x,y), a plane surface is measured twice and between the two measurements, the plane is moved by a constant distance, e.g. 10 mm. Instead of phase measurements of a flat and a curved surface, e.g. f and c, the phase measurements on two flat surfaces on different position with a constant distance, i.e. f1 and f2, are performed. Because of a constant movement, z(x,y) is known, e.g. 10 mm. Obviously, K(x,y) can be determined from equation (6): K(x,y) = ( f1 - f2 )/ z(x,y) = ( f1 - f2 )/10 (9) According to the definition of K(x,y) [cf. equation (6) and (2)], K(x,y) is defined as K(x,y) = (2π/P) k(x,y) = (2π/P) M sin θ(x,y). M is the magnification of the camera used. The M value will keep constant if the camera zoom setting does not change. It is also easy to keep a constant pitch P. Therefore, the calibration factor K(x,y) depends mainly on the angle θ(x,y). If the camera and the projector position keep constant and the distance from object surface to the measuring system doesn t change, the sensitivity factor K(x,y) will keep constant. A new calibration is required only if the position relationship among the camera, the projector and distance is changed. Of course, a new calibration is also required if any one of zoom and pitch P has been modified. (c) (d) Fig. 3 Principle of 3D-digitization for whole field 3D-shape measurement f of a flat surface (reference plane), c on 3D-object, (c) the phase difference f - c and (d) display of 3D-shape [including an operation of ( f - c )/K(x,y)] Let us look at equation (6) again. So far, we have shown how to quantitatively measure f and c as well as how to determine the sensitivity factor K(x,y). Obviously, the contour data z(x,y) can now be determined quantitatively. Fig. 3 shows an example of 3D-digitization for whole field 3D-shape measurement. Although a phase distribution f is displayed in Fig. 3, the phase distribution f of a flat surface doesn t change with different object. Therefore, this image can be saved in computer and utilized

for other measurements. That means, K(x,y) and f are determinable before measurement, only c is required for the practical measurement. For determining c, four images need to be recorded due to utilization of the phase shift technique. The time for measurement of c is about 2-3 seconds. Principle of 3D Digitization for surface distortion measurement: Surface distortion described here means a surface low, high and dent. Usually a quantitative value for such surface distortions lies between 25 to 250 microns in depth and a few millimeter in size for dents and from several decade to two hundred millimeters in size for surface lows and highs. If a surface distortion lies on a flat surface, such a distortion is relative easy to be localized by the current 3D-Digitization technique. However, most sheet metal stamping parts in automotive industries, such as door panel etc. is not a flat surface, they are usually a curved surface, a distortion on a curved surface becomes difficult to be identified by the current 3D- Digitization technique. This results from the fundamental of 3D-Digitization technique. As described above, 3D-Digitization technique measures contour information. If the object to be tested has a curved surface, both contour and surface distortion information will be measured. Usually, the contour value is much bigger that the surface distortion value, for instance, the shape value of an object is 10 mm in depth and a distortion is 40 microns. This distortion value occupied only 0.4% of the contour value. For a hardware resolution with 8 bit, there are 256 gray levels lie between the maximal and minimal value, 0.4% of the 256 gray levels is only 1 gray level. Obviously, such a surface distortion is impossible to be identified. If the contour data can be removed, the information of the surface distortion becomes visible. This is just the basic idea of 3D Digitization for surface distortion measurement. Fig. 4 Fundamental of 3D Digitization for surface distortion measurement data containing both contour and distortion, data after removing the contour information. Fig. 4 demonstrates the fundamental of this idea. For simplification only a section profile is discussed. In addition, for better understanding, the surface low has been amplified. Fig. 4a shows a measured curve contains both the contour information and a surface distortion. Our objective is to find a regression curve that represents the ideal curve of surface to be inspected. Because most surfaces of a car body can be expressed by a polynomial with order 2 as shown in equation 10, a regression surface can be created according to regression method of polynomial with order 2, f(x, y) = a x 2 + b y 2 + c xy + d x + e y + f (10) Equation 10 represents a mathematic expression of the regression surface. By subtracting the data of the measured surface from that of the regression surface, the contour data are removed and the invisible information of distortion become visible.

In practical applications, the measured data contain a lot of high-frequency noises. Therefore, before doing regression method, the original measured data should be smoothed by a specially developed smoothing algorithm 4-5. Applications: A commercial measuring system based on the above theory, called ABIS, has been developed by the company of Steinbichler GmbH, Germany. The following measuring results were obtained by utilization of ABIS. Fig. 5a shows a live image of a door panel. The area marked by red line shows the measuring area. Fig. 5c presents the measuring result that clearly shows a surface low. A section profile is demonstrated in Fig. 5c. The section profile indicated the surface low is about 40 microns in depth. The size of the surface low is about 150 mm in horizontal direction and 100 mm in vertical direction. After the measurement of the panel, it is stoned at the measuring area. The result obtained by stoning shows that a surface low is located at the exactly same position as 3D-digitazation indicated. (c) (c) Fig. 5 An example for surface distortion measurement on a door panel, live image, the area indicated by red line is the measuring area, measured result, (c) a section display, and (d) the stoned area after measurement of 3D-digitazation. Conclusions: A new technique for surface distortion has been explained and demonstrated. Because the contour data are removed, the very small information of surface low or high can be displayed very clearly. The measuring sensitivity of 3D-digitazation for surface distortion measurement is almost same as for 3D contour measurement. It can be down to 25 microns for a measuring are of 300x300 mm 2. The measuring result is well repeatable. It is expected that a wide range of applications of the technique will be seen in the near future. Reference 1. Y.Y. Hung, L. Lin, H.M. Shang, and B.G. Park, Prastical three-dimensional computer vision techniques for full-field surface measurement, Optical Engineering, Vol. 39,No.1, 143-149 (2000).

2. F. Chen, G.M. Brown, and M. Song, Overview of three-dimensional shape measurement using optical methods, Optical Engineering, 39 (1) 10 22 (2000). 3. K. Creath, Phase-Shifting Speckle Interferometry, Applied Optics, Vol. 24 (18), 3053 3058 (1985). 4. W. Steinchen, L.X. Yang, Digital Shearography; Theory and Application of Digital Speckle Pattern Shearing Interferometry, SPIE Presss, Bellingham, Wanshington USA (2003). 5. L.X. Yang, Grundlagen und Anwendungen der Phaserschieb-Shearografie zur zerstoerungsfreien Werkstoffpruefung, Dehnungsmessung und Schwingungsanalyse, VDI Fortschritt-Berichte, Reihe 8: Mess-, Steuerungs- und Reqelungstechnik, Nr. 682, VDI Verlag GmbH, Geremany (1997).