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:
where
A is the average gray scale of the image,
B is the gray scale modulation of the image, and
and
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:
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:
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.
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:
where the four coefficients
can be determined by the following constraints:
- (1)
The function is equal to the known value at each interpolation point, i.e.:
- (2)
The first and second order derivatives of the function at the two end points
and
of the interpolation interval are equal, i.e.:
Computationally, in order to solve , 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,
S2…
Si 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
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.
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
, then the point cloud center of mass
is given by:
(2) Calculate the covariance matrix.
Covariance of the two axes
representing the point cloud data distribution correlation in these two axes, the formula is as follows:
The covariance matrix
C is calculated as follows:
Compute the PCA transformation matrix.
Let the eigenvalues of the covariance matrix
C be
, and
; and let the corresponding eigenvectors be
, respectively. Then, the chi-square coordinate transformation matrix
TPCA is formulated as follows:
The coordinates of the transformed point cloud are .
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:
(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.
(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).