Next Article in Journal
Analytic Functions in a Complete Reinhardt Domain Having Bounded L-Index in Joint Variables
Previous Article in Journal
On Variance and Average Moduli of Zeros and Critical Points of Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fringe Projection Profilometry for Three-Dimensional Measurement of Aerospace Blades

1
Center of Ultra-Precision Optoelectronic Instrument Engineering, Harbin Institute of Technology, Harbin 150080, China
2
Key Laboratory of Ultra-Precision Intelligent Instrumentation Engineering, Harbin Institute of Technology, Ministry of Industry and Information Technology, Harbin 150080, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 25 January 2024 / Revised: 8 March 2024 / Accepted: 11 March 2024 / Published: 14 March 2024
(This article belongs to the Section Engineering and Materials)

Abstract

:
The aero-engine serves as the “heart” of an aircraft and is a primary factor determining the aircraft’s performance. Among the crucial components in the core of aero-engines, aero-engine compressor blades stand out as extremely important. They are not only numerous but also characterized by a multitude of parameters, making them the most complex parts in an aero-engine. This paper aims to address the trade-off between accuracy and efficiency in the existing measurement methods for asymmetric blades. Non-contact measurements were conducted using a structured light system composed of a stereo camera and a DLC projector. The point cloud data of the blades are processed using methods such as the PCA (Principal Component Analysis) algorithm, binary search, and least squares fitting. This paper established a fringe-projection profilometry light sensor system for the multi-view measurement of the blades. High-precision rotary tables are utilized to rotate and extract complete spatial point cloud data of aviation blades. Finally, measurements and comparative experiments on the blade body are conducted. The obtained blade point cloud data undergo sorting and denoising processes, resulting in improved measurement accuracy. The measurement error of the blade chord length is 0.001%, the measurement error of blade maximum thickness is 0.895%, compared to CMM (Coordinate Measuring Machine), where the measurement error of chord is 0.06%.

1. Introduction

Aero-engine blades, as one of the most crucial components of the engine, directly impact the efficiency and reliability of the entire engine through their design and performance [1,2]. With the aviation industry’s increasing emphasis on environmental impact and fuel efficiency, higher requirements are placed on the materials, design, and manufacturing technology of aero-engine blades. Blades are among the most precise and essential components of the engine, enduring high temperatures, pressures, and significant centrifugal forces. Their processing quality directly influences the overall efficiency, safety, and stability of the entire aircraft. In particular, the profiles design of the blades directly affects the engine’s efficiency, stability, and durability [3]. Therefore, the development of aero-engine blades has become a significant milestone in the progress of aviation technology. It centrally embodies cutting-edge technologies in various fields, including materials science, fluid mechanics, thermodynamics [4], and precision manufacturing.
In the aviation industry, the precise measurement of engine blades is a crucial step in ensuring their performance and reliability. With the continuous advancement of aviation technology, blade measurement methods are also constantly evolving and innovating to meet increasingly stringent accuracy requirements and challenges posed by complex shapes. Aero-engine blades are characterized by their diverse types, large quantities, complex structures, high geometric precision requirements, and challenging manufacturing processes [5]. The processing precision and manufacturing level of engine blades are crucial factors influencing engine performance, safety, and lifespan. In order to ensure the normal performance of the engine, strict requirements are placed on the manufacturing precision of the blades, which must have precise dimensions, accurate shapes, and stringent surface integrity. Blade measurement primarily involves the assessment of characteristic parameters such as chord length, maximum thickness, and chord line. Additionally, with the development of various high-performance engines, there is a significant increase in demand for blades with different shapes, sizes, and surface qualities. This imposes higher requirements on blade measurement. Therefore, high precision and rapid measurement of blade surfaces are of great significance for improving blade processing quality and ensuring the performance of aero-engines [6,7]. For the measurement methods of the profiles of aero-engine blades, they can be broadly categorized into two types: contact measurement and non-contact measurement. Contact measurement typically refers to methods that use physical contact means to obtain surface data from the blades, with the most typical example being Coordinate Measuring Machines (CMMs) [8,9]. This type of method’s main advantage lies in its high precision and reliability. Through the contact of a probing system with the blade surface, CMM can measure the geometric dimensions and shapes of the blades with extremely high accuracy. However, the drawback of contact measurement is its relatively slow speed, especially when dealing with complex or large blades [10,11]. Modern CMMs have the capability to move at high speeds, reaching several hundred mm/s. However, if accuracy is required, scanning typically needs to be conducted at slower speeds. To achieve acceptable accuracy on parts with tight tolerance, conventional scanning systems perform measurements at low speeds, usually less than 20 mm/s (0.8 in/s). Therefore, CMMs require approximately two hours to measure these blades. Non-contact measurement technologies include laser scanning, structured light scanning, and optical digitization. These methods involve using light or other radiation sources (such as X-rays) to measure the shape of blades without the need for direct contact with the blade surface. Song, X. developed an optical inspection device based on a flexible robotic arm for small compressor blades. Small compressor blades, typically blades with a chord length less than 70mm, are referred to as small blades. This device is utilized for detecting defects on the surface of AEB (Aero-engine Blade) during the production process [12]. This method enhances the measurement performance of small defects, but the measurement apparatus is relatively complex. Fringe projection profilometry is an advanced non-contact optical measurement technology commonly uses for accurately measuring objects with complex geometric shapes, such as aero-engine blades. This technology combines structured light systems with binocular vision principles, enabling the precise reconstruction of the three-dimensional surface profiles of the measured object. Ma Kang Sheng [13] used a multi-frequency heterodyne-phase-shift method with an automatic positioning device for in situ measurement simulation, obtaining the original point cloud data of the blades. Ma Qian Li [14] proposed an algorithm for optimizing the selection of the optimal measurement angles in binocular-structured light, based on visual cones and vector optimization using measurement points. Using an industrial robotic arm to position the structured light at the selected measurement points, multi-view scanning measurements were conducted on aero-engine blades.
Liu X. developed a new technique combining composite stripe coding and stepped phase coding to solve high-quality absolute phase and realize deep learning-based binocular 3D profiles design reconstruction. The method only needs to project the composite stripe coding pattern to obtain the wrapped phase and the stepped phase coding pattern to obtain the stripe order, which solves the effects of defocus and noise on the deviation of wrapped phase and stripe order [15]. For the issue of fewer blade samples, a metric learning for the multi-output Gaussian process regression method (ML_MOGPR) for aero-dynamic performance prediction of the plane cascade is proposed. It shares parameters between multiple output Gaussian distributions during training and measures the similarity between input samples in a new embedding space to reduce bias and improve overall prediction accuracy [16]. For actively controlling non-linear blade vibrations, a rotating blade’s non-linear oscillations are reduced via a time-delayed, non-linear saturation controller (NSC) [17].
After measuring and obtaining blade data, it is necessary to obtain blade characteristic parameters defined under standards for evaluation. Li, X. proposed an improved YOLOv5 method called DDSC-YOLO [18], which integrates several aspects, such as shape feature extraction, computational effort, and measurement performance. Ching Hin Lydia Chan [19] investigated the blade checkpoints’ optimal distribution that will produce the most accurate B-spline fit to the surface of the CAD model, thus shortening the inspection and analysis process without compromising the accuracy. Jing-Jun Li proposed an ultrasonic defect in aero-engine blade point cloud [20] for quantitative evaluation, in which the regionally discretized defective point cloud is clustered into independent defective regions by the DBSCAN clustering algorithm. An averaged differential imaging method is presented and applied to the ultrasonic nondestructive evaluation of wind turbine blades [21]. In the image, we can successfully find line images with damaged parts, and we can estimate the depth of the damaged parts with high accuracy.
Generally, the methods for measuring blade profiles are more perfect at present, but the efficiency of the measurement will inevitably be reduced while ensuring the accuracy of the measurement. In this paper, we aim to solve the problem that the existing blade measurement scheme cannot balance accuracy and efficiency and will adopt a measurement system combining a fringe projection profilometry and a high-precision rotary table, which realizes the function of measuring a blade with high efficiency, precision, and stability. For the aero-engine blade parameters, a synchronized noise reduction algorithm based on the moving least squares method for sorting the blade’s cross-section point cloud and a bifurcated search algorithm for mid-arc are proposed to fit the blade parameters, and the blade parameters are evaluated without the use of complex formulas to fit the blade profiles. Finally, several sets of comparison experiments with the CMM coordinate machine are conducted to verify that the accuracy meets the measurement requirements.
In this paper, we first obtained and unwrapped the phase containing the blade profile information based on the principle of the fringe projection profilometry method and utilized the spline difference for stereo matching. In order to integrate the measured data into the same coordinate system, the rotary table calibration method was proposed. In order to obtain the characteristic parameters of the blade, the cross-section of the blade was extracted using PCA analysis, the point cloud data were sorted, and the noise was reduced in order to improve the accuracy. Finally, the blade body was measured and evaluated against the standard values. In order to verify that the accuracy requirements were met, several sets of comparison experiments with the CMM were conducted, and the measurement requirements were met.

2. Measurement Principle

2.1. Measurement System Overview

The structure of the proposed measurement system based on fringe projection profilometry and an air-bearing rotary table is given in Figure 1. The sample of the asymmetric blade to be measured is mounted on a lifting table in the center of the rotary table with two degrees of freedom guide rails to adjust the horizontal as well as the vertical measurement position of the binocular-structured light sensor. The relative spatial relationship between the rotary table and the binocular-structured light sensor remains constant during the measurement process. A single measurement at a preset angular position and repeated measurements at equal intervals in conjunction with the rotation of the rotary table realize the multi-view rotational measurement of the blade.
The asymmetric blade measurement system is mainly composed of three parts: measurement device part, control and data processing part, and data interaction part. The measuring device mainly consists of a binocular-structured light sensor and a high-precision rotary table. The control and data processing section is an industrial-control computer system. The blade profile measurement needs to obtain the complete cross-section of the blade profile, so it is necessary to use a rotary table to extend the structured light measurement to multiple viewing angles; so, this paper designs a blade profile measurement system with a binocular-structured light and air-bearing rotary table.

2.2. Phase-Shift Method–Principle of Multi-Wavelength Phase Unwrapping

The phase-shifting method projects multiple pairs of phase-shifting sine fringe images within a cycle with a DLP projection camera, and then calculates the phase information of each pixel of the phase-shifting image. This phase information essentially encodes the height information, meaning that the height of the object’s surface modulates the phase value. However, these phase details are wrapped within multiple cycles, hence referred to as the wrapped phase. By calculating the wrapped phase information, the height information of the object can be calculated. According to the mathematical model, the light intensity formula of the phase-shifting fringe image obeys the standard sine distribution as follows:
I x , y = A + B cos φ x , y + φ k
where A is the average gray scale of the image, B is the gray scale modulation of the image, and φ and φ k are the relative phase code value and phase shift value, respectively. The overdetermined solution is computed by calculating the phase-shift equation in more than four steps. After calculating the final phase shift method, the parcel phase is calculated as follows:
φ x = a r c tan B 2 x , y B 1 x , y = a r c tan k = 0 N 1 ( I k sin ( 2 k π N ) ) k = 0 N 1 ( I k cos ( 2 k π N ) )
The wrapped phase is the result obtained after the phase information is modulo 2 π operation. In the phase-shift method, different phase-shift images are obtained by multiple optical projections of the measured object, and these phase-shift images are processed by the above formula to finally obtain the wrapped-phase image of the measured object, and the phase-difference image can be obtained by subtracting two wrapped-phase images, thus realizing the calculation of the surface profile of the object.
For a higher measurement accuracy, this paper chooses the number of phase-shift steps N to be 12 and uses three modulation frequencies for phase expansion so as to minimize the effect of nonlinear distortion. By using a set of phase-shift sequences containing 12 discrete phase-shift values, the original signal is phase shifted to obtain 12 phase-shifted signals. By performing amplitude ratio operations on these 12 signals, a wrapped-phase signal containing multiple phase jumps is obtained. Then, by unwrapping this wrapped-phase signal, the continuous phase information of the original signal can be obtained.
The basic Idea of the method is to obtain an alternate signal with a lower frequency by outperforming two sinusoidal signals with a small difference in frequency, thus realizing the phase unwrapping of the original signal. This method can effectively deal with signals whose phases contain a large number of phase cycles and improves the stability and accuracy of phase unwrapping. The mathematical principle of the dual-frequency outlier is shown in Figure 2.
The schematic diagram of the dual-frequency outlier principle is shown in Figure 2, where φ 1 and φ 2 are the phase functions of the two small cycles, while φ 12 is the phase function of the larger cycle calculated by the outlier method of φ 1 and φ 2, and its value is the least common multiple of φ 1 and φ 2. Assuming that φ 1 < φ 2, the formula is as follows:
φ 12 = φ k φ 2 1 T 12 = 1 T 1 1 T 2 T 12 = T 1 T 2 T 1
Similarly, based on the results of the dual-frequency outlier, the multi-frequency heterodyne can be calculated to achieve the effect of unfolding the global phase.
T 23 = T 2 T 3 T 3 T 1
T 123 = T 12 T 23 T 23 T 12
Therefore, the multi-frequency heterodyne method requires three different cycles of the parcel phase.

2.3. Cubic Spline Interpolation Subpixel Phase Value Stereo Matching Algorithm

The pixel values in the left and right views are interpolated by the cubic spline interpolation method, and the sub-pixel level matching results are obtained. The phase values of the pixels in the views obtained by the left and right cameras are calculated and interpolated by the sub-pixel matching results to obtain the phase values at the sub-pixel level. Based on the binocular camera parameters and the phase values, the 3D coordinates of the corresponding pixels are calculated. To obtain the 3D-point cloud data of the object surface, the 3D coordinates for all pixels are calculated.
The main advantage of the stereo matching algorithm for sub-pixel phase values based on cubic spline interpolation is that it can make full use of the sub-pixel level information between pixels, thus improving the measurement accuracy. Meanwhile, the use of cubic spline interpolation can make the interpolation results smoother and avoid unnecessary noise caused by interpolation errors. The mathematical principle of cubic spline interpolation is as follows: let the pixel horizontal coordinates of the interpolated point be x0, x1…, xn, and the corresponding phase values be y0, y1…, yn. A cubic polynomial is used to approximate the interpolated point in the following form:
s i x = a i + b i x x i + c i x x i 2 + d i x x i 3 , x i x x i + 1
where the four coefficients a i , b i , c i , d i can be determined by the following constraints:
(1)
The function is equal to the known value at each interpolation point, i.e.:
s i x i = y i
(2)
The first and second order derivatives of the function at the two end points x i and x i + 1 of the interpolation interval are equal, i.e.:
s i k x i = s i 1 k x i , s i k x i + 1 = s i + 1 k x i + 1 , k = 0,1 , 2
Computationally, in order to solve a i , b i , c i , d i , a system of linear equations needs to be constructed, and in order to simplify this problem, iterative solution methods are often used to approximate the solution of the system of equations.

2.4. Calibration Method of Relative Position of Rotary Table

In order to integrate all measurements of the blade into the same coordinate system, a standard sphere is used to determine the spatial attitude of the turntable. A standard sphere is configured on the turntable for rotational scanning. A binocular-structured light sensor scans the sphere defined on the turntable to obtain the center position of the sphere at each measurement position. The positions are determined by illuminating the standard sphere using data from the measurement portion of the binocular-structured light beam. As shown in Figure 3, S1, S2Si are the measurement positions and i is the number of measurements. A schematic diagram of the rotary table calibration device is shown in Figure 3.
As shown in Figure 4, the fitted center O′ (x0, y0, z0) and normal vectors E (a, b, c) of the sequence c 1 ~ c N j are obtained using the least squares fitting method. As shown in Figure 5, O′ and E are the center and direction of the rotation axis, respectively. Here, E (a, b, c) is reduced to two parameters, α and β. The angle α starts from the positive direction of the y-axis and rotates counterclockwise.
α = a r c c o s b a 2 + b 2 ,                     a 0 2 π a r c c o s b a 2 + b 2 ,   a < 0 β = a r c c o s c a 2 + b 2 + c 2 .

3. Measurement Method

3.1. Definition of Blade Body Characteristic Coefficients

According to the industrial standard, the parameters to be extracted from the blade body are shown in Figure 6:
In this paper, the blades we measured are made of high-temperature nickel-based alloy. Several parameters that have the greatest impact on the performance of aero-engine blades are selected for evaluation: chord, chord length, mid-arc, and maximum thickness of the blade. Their definitions are shown in Table 1.

3.2. Extraction of Blade Cross Sections by PCA Method

In this paper, a Principal Component Analysis (PCA algorithm) is used to determine the benchmarks for the evaluation of the measurement data. The PCA algorithm is a set of linearly uncorrelated variables converted into a set of variables that may be correlated by orthogonal transformation, so that multiple variables in the data are independent of each other, and the converted variables are called principal components. For point cloud data, the principal components represent the degree of dispersion between the three coordinates, and the PCA algorithm establishes the characteristic coordinate system by the direction with the largest variance of the point cloud coordinates. The direction with the largest variance of point cloud coordinates is the Z-axis of the blade; thus, the PCA algorithm can be used to determine the base coordinate system of the measurement data. The main process of the PCA algorithm is as follows:
(1) Calculate the center of mass position of the point cloud.
Let the blade-point cloud data be C = P 1 , , P Ν , P n = x n , y n , z n , then the point cloud center of mass m = ( x ¯ , y ¯ , z ¯ ) is given by:
m = 1 N n = 1 N P n = 1 N n = 1 N x n , 1 N n = 1 N y n , 1 N n = 1 N z n ,
(2) Calculate the covariance matrix.
Covariance of the two axes cov ( x , y ) representing the point cloud data distribution correlation in these two axes, the formula is as follows:
cov ( x , y ) = cov ( y , x ) = 1 n 1 i = 1 n 1 x i x ¯ y i y ¯
The covariance matrix C is calculated as follows:
C = cov ( x , x ) cov ( x , y ) cov ( x , z ) cov ( y , x ) cov ( y , y ) cov ( x , z ) cov ( z , x ) cov ( z , y ) cov ( z , z )
Compute the PCA transformation matrix.
Let the eigenvalues of the covariance matrix C be λ 1 , λ 2 , λ 3 , and λ 1 , λ 2 , λ 3 ; and let the corresponding eigenvectors be ε 1 , ε 2 , ε 3 , respectively. Then, the chi-square coordinate transformation matrix TPCA is formulated as follows:
T P C A = x ¯ e 2 e 3 e 1 y ¯ z ¯ 1
The coordinates of the transformed point cloud are C P C A = T P C A C .

3.3. Point Cloud Processing and Noise Reduction Sorting

In order to evaluate the cross-section features subsequently, the point cloud needs to be partitioned, and the ordered point cloud can effectively improve the efficiency of the subsequent noise reduction and partitioning algorithm. During the measurement process, due to the angle between the structured light and the main direction of the blade, it is not possible to measure the complete and organized cross-section line data within the same measurement. When the fringe projection profilometry is used to measure the point cloud of the blade, the number of noise points is large, and the result of processing the point cloud will have a great impact; so, it is necessary to perform noise reduction and smoothing filtering on the original point cloud data collected.
The basic steps of the aero-engine blade sequencing synchronous noise reduction method based on the moving least squares method proposed in this paper are as follows:
(1) Weight selection: Select the appropriate threshold H. Starting from the first point of the original data, which is also the first point after sorting, set this point as the current point with serial number i, transverse the remaining point pj, calculate its distance r from this point, and the weight of the remaining point is:
w j = e r 2 H 2
(2) Coordinate transformation: Fit a local least squares line with the weights and compute a rotation matrix to rigidly transform the entire cross-section point cloud to a new coordinate system with the line as the x-axis and the current point as the origin.
(3) Noise reduction: Recalculate the weights of the remaining points after the transformation is completed, fit the parabola with the new weights and replace the current point with the intersection of the parabola and the y-axis of the current coordinate system.
w j = 2 r 3 H 3 3 r 2 H 2 + 1       r H     0           r > H
(4) Sorting: Using the current coordinate system as a reference, the x-axis in a positive direction from the origin of the nearest point for the i + 1 point, as shown in the blue line, represents the current order of connection. The left side of the red circle points to complete the sorting, and the right side is to be sorted. Repeat steps (1) to (4).

4. Experiments

When using fringe projection profilometry to obtain the point cloud of the blade, the number of noise points is large, and it has a significant impact on processing the point cloud. Therefore, it is necessary to perform noise reduction and smoothing filtering on the collected raw point cloud data. First of all, the measurement system is constructed by using a binocular-structured light sensor and high-precision rotary table, and the relative position between the sensor measurement coordinate system and the relative position of the rotary table in the measurement coordinate system are calibrated; then, the aero-engine blade is measured, and the blade profile parameters are obtained from the specified profile after the data splicing, processing, and blade profile characteristic parameter extraction. The blade profile parameters are measured by scanning the measurement system are compared with the three coordinate measurement results to verify the effectiveness of the measurement method proposed in this paper. Finally, the blade parameters measured by the scanning measurement system are compared with the CMM results to verify the effectiveness of the measurement method proposed in this paper. The physical diagram of the blade measurement system is shown in Figure 7.
The measurement system utilizes two cameras with dimensions of 36 mm × 31 mm × 38.8 mm, and the specific parameters of the cameras are listed in Table 2 below.

4.1. Calibration Experiment

An air-float rotary table and a binocular-structured light sensor are used as the main equipment of the proposed measurement system. As shown in Figure 8, a standard sphere is used for the calibration experiment. The standard sphere is placed on the descending table and the air-bearing turntable is rotated by 30° to obtain a set of point clouds of the sphere. As shown in Figure 9, the left and right panels are the pictures obtained by photographing the standard sphere during the calibration of the binocular camera, respectively. The least squares method was used to determine the fitting center O′ (x0, y0, z0) and the fitting plane, which are equal to Zo′ (−21.5382, 0.0423, 1.1557) mm and −0.4428x + 0.0206y − z = 8.3823, i.e., α = 272.6578° and β = 156.0973, and the method described in Section II-D is used to compute the central sequence of the standard spheres C1-CJ in conjunction with the radius of the reference sphere.
We utilized a standard sphere with a known center-to-center distance of 100.125 mm for the turntable calibration experiment. The turntable calibration results were validated by moving the turntable and measuring the standard sphere multiple times. The point cloud data obtained from multiple measurements were stitched together using the turntable calibration results. The fitted sphere center distance was then compared with the known standard sphere center distance. The results, as shown in Table 3, indicate a final measurement error of 0.0329 mm for the sphere center distance, meeting the measurement requirements. The table depicting the tolerance of the blade contour is presented in Table 4.

4.2. Blade Body Measurement and Evaluation

Measurement blades are shown in Figure 10.
After the calibration of the measurement system is completed, the blade scanning measurement is carried out. The blade is adjusted so that it lies within the measurement range of the binocular camera, the appropriate sampling interval for the structured light sensor is set, and measuring from the front side of the blade begins until the entire blade contour is measured. After the blade profile measurement is completed, the splicing of the structured light sensor measurement profile is completed according to the results of the rotary table movement direction calibration, and the whole splicing process is based on the initial measurement position of the blade.
The phase-shift method–principle of multi-wavelength phase unwrapping method described in Section 2.2 is used to calculate the parcel phase, and the results are shown in Figure 11. As shown in Figure 11, the transition from red to blue represents phase values from high to low. According to the discussion in Section 2.2, for this experiment, we chose T1 = 28, T2 = 26, T3 = 24, in units of the pixels of the projection light engine. It can be observed from the figure that accurate phase unwrapping has been achieved, resulting in a well-preserved global phase.
The PCA principal component analysis method described in Section 3.2 is used to determine the benchmark for measurement data evaluation. After repeated measurements through the rotary turntable to obtain the original point cloud after pre-processing and multi-view scanning of the aero-engine blade, in order to extract the final parameters, the cross-section needs to be extracted. The blade cross-section point cloud extraction algorithm is used, and the PCA algorithm is used to align the original point cloud of the aero-engine blade to the datum. The scanned point cloud is aligned to the original point cloud by the global matching algorithm, as shown in Figure 12, and the white point cloud is the scanned point cloud. The original point cloud, after the alignment of the cross-section projection, takes two intervals of 0.3 mm perpendicular to the PCA alignment of the coordinate system after the Z-axis of the plane interception figure in the green area of the cross-section of the point cloud, as shown in Figure 12.
In order to accurately extract the parameters, the cross-section is filtered, noise is reduced, and the 2D-point cloud of the height section is sorted at the same time, as shown in Figure 13.
The parabolic segmentation method is used for the blade basin and blade dorsum to segment the blade basin and blade dorsum, after extracting the mid-arc by the bisecting search method, as shown in the orange curve in Figure 14, and the inner tangent circle at the maximum thickness is shown in the green circle in Figure 14.
The results of this are evaluated on two parameters: chord length and maximum thickness. Random cross-sections of different heights were evaluated, and the error of the maximum thickness was 0.096 mm. The experimental data are shown in Table 5 and Table 6 below.

4.3. Comparative Blade Profiles Experiment

In order to verify whether the accuracy of the measurement system can meet the requirements, a comparative experiment on the same set of blades was carried out with a CMM, which is currently known to have a high measurement accuracy. The blades were fixed on a rotary table, and the CMM probe was contacted to obtain 3D data of the blades.
A least squares fit was used to fit the cross-section of the blade at a certain height and the comparison data are shown in Table 7 below. The precision of the coordinate measuring machine (CMM) we used is within 1 micron.

5. Conclusions

This paper proposes a blade profile measurement method based on fringe projection profilometry and a rotary table. The measurement system adopts a binocular camera to ensure that the measurement accuracy meets the requirements while effectively improving the measurement efficiency of the blade. The method makes full use of the advantages of fringe projection profilometry to reconstruct the overall macroscopic shape information of the blade surface and carries out a multi-view profile measurement of the blade through the rotary table. Such a combination can overcome the problem of insufficient accuracy of fringe projection profilometry at the microscopic scale. The results of our study can be summarized as follows:
(1)
Completed the establishment of the fringe projection profilometry system and the air flotation rotary table model, according to the model to achieve the rotary motion control of the blade; the use of fringe projection profilometry with its high precision, high resolution, and high efficiency maximize the degree of simplification of the measurement process, and at the same time to meet the requirements of precision.
(2)
Aiming at the nonlinear displacement problem of the projector, in order to reduce the nonlinear error of the projector, a wrapped-phase calculation method based on the twelve-step phase-shift method is proposed, which improves the accuracy of the wrapped-phase solution, and a phase-unwrapping algorithm based on the multi-frequency outlier method is adopted, which reduces the problem of the phase cycle hopping error from the principle, and improves the accuracy of phase unwrapping.
(3)
A simultaneous noise-reduction–sorting algorithm is proposed. Noise reduction is accomplished while sorting, and uncomplicated formulas can be used to fit the surface of the blade and the point cloud can be rated directly. The study of sorting and noise reduction algorithms on the blade cross-section is accomplished, and the experimentally obtained point cloud data of the blade cross-section are filtered out of the noise.
(4)
Measurement and comparison experiments are carried out on the blade, and the accuracy of the binocular system measurement is verified using the CMM, and the maximum error of the blade maximum thickness measurement is verified to be 1.626% through experiments. A chord search algorithm was proposed to extract the chord and the chord length, the maximum error of the chord length extraction is verified to be 0.142% through experiments, and the error of the CMM is 0.06%.
Prospect: In this paper, the aero-engine blade profile has been investigated, but there are still some aspects that need to be investigated in more depth, for example, the sensitivity to ambient light and reflections on specific surface materials. Aero blades are usually made of metal or are coated with reflective materials, which may lead to high reflections of structured light on the blade surface, thus interfering with the measurement accuracy. Subsequent polarization filters may be used, or special surface treatment techniques may be applied to reduce reflections. Ambient light in laboratory or factory environments may interfere with the projection and capture of structured light, and we can only control the ambient light conditions by using filters or making measurements in darker environments.

Author Contributions

Writing—original draft, methodology, Z.C.; writing—original draft, visualization, M.Z.; data curation, validation, M.Z.; writing—original draft, experiments, C.S.; writing—review and editing, funding acquisition, Y.L.; data Curation, validation, M.Z.; writing—review and editing, funding acquisition, J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Heilongjiang Province of China (grant number LH2021E060), National Key R&D Program of China (grant number 2021YFF0603200), National Natural Science Foundation Major Research Projects of China (grant number 91960109), National Natural Science Foundation of China (grant number 52275525, 52175498, 52205560, 51975158), and the Fundamental Research Funds for the Central Universities (grant number 2022FRFK060025, XNAUEA5750302321).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

We declare that we have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Amoo, L.M. The design and structural analysis of jet engine fan blade structures. Prog. Aerosp. Sci. 2013, 60, 1–11. [Google Scholar] [CrossRef]
  2. Naghmehsanj, M.; Rahmani, B. Optimal control of supersonic pre-twisted rotating functionally graded thin-walled blades. Struct. Control Health Monit. 2017, 24, e1957. [Google Scholar] [CrossRef]
  3. Li, Z.; Wang, Z.; Huang, J.; Duan, Q.; Qi, M.; Gao, J.; Wang, W.; Dong, Q.; Li, Q.; Ai, S. Blade profile extraction and edge completion method based on structured light measurement point cloud. Precis. Eng. 2023, 86, 225–238. [Google Scholar] [CrossRef]
  4. Mishnaevsky, L., Jr.; Branner, K.; Petersen, H.N.; Beauson, J.; McGugan, M.; Sørensen, B.F. Materials for wind turbine blades: An overview. Materials 2017, 10, 1285. [Google Scholar] [CrossRef] [PubMed]
  5. Kharyton, V. Faults Detection in Blades of an Aviation Engine in Operation. Doctoral Dissertation, École Centrale de Lyon, Écully, France, 2009; pp. 1–4. [Google Scholar]
  6. Feng, G.; Ziyue, P.; Xutao, Z.; Yan, L.; Jihao, D. An adaptive sampling method for accurate measurement of aeroengine blades. Measurement 2021, 173, 108531. [Google Scholar] [CrossRef]
  7. Pan, T.; Yan, Z.; Lu, H.; Li, Q. Numerical investigation on the forced vibration induced by the low engine order under boundary layer ingestion condition. Aerosp. Sci. Technol. 2021, 115, 106802. [Google Scholar] [CrossRef]
  8. Chang, H.C.; Lin, A. Automatic inspection of turbine blades using a 3-axis CMM together with a 2-axis dividing head. Int. J. Adv. Manuf. Technol. 2005, 26, 789–796. [Google Scholar] [CrossRef]
  9. Makem, J.; Ou, H.; Armstrong, C.; Rennie, A.; Nikov, S. A Virtual Inspection Technique for Assessing the Dimensional Accuracy of Forged Compressor Blades Using FE Modeling and CMM Inspection. Int. J. Mater. Form. 2008, 1 (Suppl. S1), 375–378. [Google Scholar] [CrossRef]
  10. Lin, X.; Jiang, S.; Liu, X.; Yang, K. The CMM Measurement Path Planning for Blade Surface Based on the Contour Measurement. In Proceedings of the 2011 Second International Conference on Digital Manufacturing & Automation, Zhangjiajie, China, 5–7 August 2011; pp. 1228–1232. [Google Scholar] [CrossRef]
  11. Wei, X.; Li, B.; Yan, X.; Chen, L.; Su, W.; Xin, M.; Shang, Z. Calibration of Sensor Pose Error in Aero-Engine Blade Measurement. In Proceedings of the 2019 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Auckland, New Zealand, 20–23 May 2019; pp. 1–6. [Google Scholar] [CrossRef]
  12. Song, K.; Sun, X.; Ma, S.; Yan, Y. Surface Defect Detection of Aeroengine Blades Based on Cross-Layer Semantic Guidance. IEEE Trans. Instrum. Meas. 2023, 72, 2514411. [Google Scholar] [CrossRef]
  13. Li, M.; Zhao, W.; Ma, K.; Wang, F. Automatic Adjustment Method for Point Cloud Simplification and Reconstruction Parameters in Structured Light Detection. Chin. J. Sci. Instrum. 2022, 43, 122–130. [Google Scholar] [CrossRef]
  14. Liu, T.; Yu, Z.; Ma, Q. Target distance calculation based on parallel binocular vision. J. Qingdao Univ. Nat. Sci. Ed. 2009, 22, 59–62. [Google Scholar]
  15. Liu, X.; Yang, L.; Chu, X.; Zhou, L. A novel phase unwrap method for binocular structured light 3D reconstruction based on deep learning. Optik 2023, 279, 170727. [Google Scholar] [CrossRef]
  16. Liu, L.; Yang, C.; Xiang, H.; Lin, J. Plane Cascade Aerodynamic Performance Prediction Based on Metric Learning for Multi-Output Gaussian Process Regression. Symmetry 2023, 15, 1692. [Google Scholar] [CrossRef]
  17. Hamed, Y.S.; Kandil, A. Influence of Time Delay on Controlling the Non-Linear Oscillations of a Rotating Blade. Symmetry 2021, 13, 85. [Google Scholar] [CrossRef]
  18. Li, X.; Wang, W.; Sun, L.; Hu, B.; Zhu, L.; Zhang, J. Deep learning-based defects detection of certain aero-engine blades and vanes with DDSC-YOLOv5s. Sci. Rep. 2022, 12, 13067. [Google Scholar] [CrossRef] [PubMed]
  19. Chan, C.H.L.; Wang, Q.; Holden, R.; Huang, S.; Zhao, W. Optimal Number of Control Points for Fitting B-Splines in Wind Turbine Blade Measurement. Int. J. Precis. Eng. Manuf. 2019, 20, 1507–1517. [Google Scholar] [CrossRef]
  20. Li, J.-J.; Yan, C.-F.; Rui, Z.-Y.; Zhang, L.-D.; Wang, Y.-T. A Quantitative Evaluation Method of Aero-engine Blade Defects Based on Ultrasonic C-Scan. In Proceedings of the 2020 IEEE Far East NDT New Technology & Application Forum (FENDT), Kunming, China, 20–22 November 2020; pp. 91–95. [Google Scholar] [CrossRef]
  21. Tsukuda, K.; Egawa, T.; Taniguchi, K.; Hata, Y. Average difference imaging and its application to ultrasonic nondestructive evaluation of wind turbine blade. In Proceedings of the 2012 IEEE International Conference on Systems, Man, and Cybernetics (SMC), Seoul, Republic of Korea, 14–17 October 2012; pp. 2601–2604. [Google Scholar] [CrossRef]
Figure 1. Components and interconnection of the proposed measurement system.
Figure 1. Components and interconnection of the proposed measurement system.
Symmetry 16 00350 g001
Figure 2. Principle diagram of phase unwrapping using dual frequency heterodyne method.
Figure 2. Principle diagram of phase unwrapping using dual frequency heterodyne method.
Symmetry 16 00350 g002
Figure 3. Diagram of experimental setup for turntable calibration.
Figure 3. Diagram of experimental setup for turntable calibration.
Symmetry 16 00350 g003
Figure 4. Determination of the axis center O′ (x0, y0, z0) and normal vector E (a, b, c).
Figure 4. Determination of the axis center O′ (x0, y0, z0) and normal vector E (a, b, c).
Symmetry 16 00350 g004
Figure 5. Spatial direction E of the axis of the rotary table.
Figure 5. Spatial direction E of the axis of the rotary table.
Symmetry 16 00350 g005
Figure 6. Blade profile parameters.
Figure 6. Blade profile parameters.
Symmetry 16 00350 g006
Figure 7. Physical diagram of blade measurement system.
Figure 7. Physical diagram of blade measurement system.
Symmetry 16 00350 g007
Figure 8. Conduct calibration experiments using a pair of standard spheres.
Figure 8. Conduct calibration experiments using a pair of standard spheres.
Symmetry 16 00350 g008
Figure 9. Phase map captured by left and right cameras.
Figure 9. Phase map captured by left and right cameras.
Symmetry 16 00350 g009
Figure 10. Measuring blade.
Figure 10. Measuring blade.
Symmetry 16 00350 g010
Figure 11. (a) Multi-frequency heterodyne unpacking. (b) Adjustment chart.
Figure 11. (a) Multi-frequency heterodyne unpacking. (b) Adjustment chart.
Symmetry 16 00350 g011
Figure 12. Extraction results of aero-engine blade cross-section.
Figure 12. Extraction results of aero-engine blade cross-section.
Symmetry 16 00350 g012
Figure 13. Synchronous noise reduction results of aero-engine blade section sorting.
Figure 13. Synchronous noise reduction results of aero-engine blade section sorting.
Symmetry 16 00350 g013
Figure 14. Search results for midline and maximum blade thickness.
Figure 14. Search results for midline and maximum blade thickness.
Symmetry 16 00350 g014
Table 1. Names and definitions of blade body characteristic parameters.
Table 1. Names and definitions of blade body characteristic parameters.
ParameterDefinition
cordA straight-line tangent to both the leading and trailing edges of the blade
chord lengthThe maximum length of the blade cross-section in the chord direction
camber lineThe line connecting the maximum inscribed circle center of the blade section
Maximum thickness of blade profileThe diameter of the tangential circle within the blade section
Table 2. Digital camera parameters.
Table 2. Digital camera parameters.
Frame RatePixel CountPixel SizePixel DepthExposure Time
15.1 fps26,200 thousand2.5 µm × 2.5 µm8 bit, 12 bit11 μs~1 s
Table 3. Turntable calibration experiment.
Table 3. Turntable calibration experiment.
No.Distance (mm)Error (mm)
1100.10770.0172
2100.13870.0137
3100.20.075
4100.1690.0439
5100.11020.0147
Avg100.145120.0329
Table 4. Tolerance of blade contour.
Table 4. Tolerance of blade contour.
Nominal Size of Chord Length (mm)Accuracy GradeCompressor BladesTurbine Blades
EdgeMidsectionEdgeMidsection
≤1610.050.050.100.13
20.060.060.130.16
30.080.080.160.20
40.100.100.200.25
>16~2510.050.060.130.16
20.060.080.160.20
30.080.100.200.25
40.100.130.250.31
>25~4010.060.080.160.20
20.080.100.200.25
30.100.130.250.31
40.130.160.310.40
>40~6410.060.100.160.25
20.080.130.200.31
30.100.160.250.40
40.130.200.310.50
Table 5. Chord length measurement results.
Table 5. Chord length measurement results.
Chord Length ReferenceMean Measurements ValueDifference Value
5453.92320.142%
53.7553.71760.060%
53.553.49920.001%
53.2553.24590.008%
5352.99940.001%
Table 6. Maximum thickness measurement results.
Table 6. Maximum thickness measurement results.
Maximum Thickness Reference ValueMean Measurements ValueDifference Value
2.82.821780.778%
2.92.947161.626%
33.030241.008%
3.13.114840.479%
3.23.18130.584%
Table 7. Comparison results.
Table 7. Comparison results.
No.Chord Measurement
CMMCameraError
152.79652.7700.05%
252.81552.7640.10%
352.71552.7710.11%
452.74352.7620.04%
552.81252.7630.09%
652.75452.7590.01%
752.72252.7740.10%
852.78252.7550.05%
952.71652.7730.11%
1052.76752.7620.01%
Avg52.76252.7650.01%
Std0.0390.006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Z.; Zhu, M.; Sun, C.; Liu, Y.; Tan, J. Fringe Projection Profilometry for Three-Dimensional Measurement of Aerospace Blades. Symmetry 2024, 16, 350. https://doi.org/10.3390/sym16030350

AMA Style

Chen Z, Zhu M, Sun C, Liu Y, Tan J. Fringe Projection Profilometry for Three-Dimensional Measurement of Aerospace Blades. Symmetry. 2024; 16(3):350. https://doi.org/10.3390/sym16030350

Chicago/Turabian Style

Chen, Ze, Mulun Zhu, Chuanzhi Sun, Yongmeng Liu, and Jiubin Tan. 2024. "Fringe Projection Profilometry for Three-Dimensional Measurement of Aerospace Blades" Symmetry 16, no. 3: 350. https://doi.org/10.3390/sym16030350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop