Compressor impeller efficiency consistency detection method

Abstract

The invention discloses a method for detecting the efficiency consistency of a gas compressor impeller. The method comprises the steps that 1, sample preprocessing and point cloud scanning are carried out; step 2, color space conversion: converting an RGB color space into an HSV color space and utilizing the characteristic of chromaticity H dimension; 3, deviation parameter calculation: calculating the overall level of the average deviation quantization leaf shape deviation, calculating the extreme degree of the maximum deviation quantization leaf shape deviation, and calculating the sensitive area distribution density of the large deviation ratio quantization leaf shape deviation; step 4, predicting efficiency deviation, and converting geometric deviation data into a visual performance deviation result by establishing a quantitative correlation model of leaf shape deviation and efficiency deviation; and step 5, judging performance consistency according to an efficiency deviation result. By establishing a quantitative correlation model of leaf shape deviation and efficiency deviation and combining color space conversion, precise deviation extraction is achieved, and finally efficient and high-precision performance consistency detection is achieved.

Description

Technical Field

[0001] This invention belongs to the field of testing technology for internal combustion engine turbocharger components, and particularly relates to a method for testing the performance consistency of compressor impellers, which is suitable for fast-paced and high-precision testing of compressor impeller performance consistency in mass production scenarios. Background Technology

[0002] The compressor impeller is the core moving component of an internal combustion engine turbocharger. Its blade profile accuracy directly determines key performance indicators such as post-compression pressure and efficiency, thereby affecting engine power, economy, and emissions compliance. Currently, customers have stringent requirements for turbocharger performance consistency, such as controlling post-compression pressure variation within ΔP=±3kPa and strictly limiting compressor efficiency deviation. Simultaneously, emissions regulations are placing increasing demands on the sensitivity of turbocharger performance, requiring that performance fluctuations of each impeller be within a controllable range.

[0003] Current methods for testing the performance consistency of compressor impellers mainly rely on three types of methods, but all of them have significant drawbacks: 1. Experimental testing method: Pressure ratio and efficiency deviation are obtained through bench testing of the booster. However, the testing accuracy is limited by the bench equipment (e.g., the test deviation of the Kaze test bench is ≥ ±0.55%), which cannot meet the theoretical efficiency deviation requirement of ±0.5%. Moreover, bench testing is costly and time-consuming (testing a single unit takes ≥ 2 hours), making it difficult to adapt to the needs of mass production.

[0004] 2. Simulation calculation method: The performance is predicted by using the process of "point cloud scanning-reverse modeling-CFD calculation". However, the reverse modeling process will generate deviation amplification. The accumulation of reverse deviation will lead to the deviation between the CFD calculation results and the actual performance, resulting in low simulation credibility.

[0005] 3. Blade shape detection method: Based on geometric deviation standards (such as the maximum blade shape deviation of cast impeller ±0.2mm), but this standard only focuses on the maximum blade shape deviation value and does not consider the influence of deviation location (such as deviations at the leading edge and trailing edge, which will increase boundary layer separation loss) and deviation distribution on performance. There are cases of "geometric compliance but performance exceeding the tolerance" and "geometric non-compliance but performance not exceeding the tolerance", so the standard is not reasonable enough.

[0006] In summary, existing testing methods cannot simultaneously meet the testing requirements of "high precision, high efficiency, and low cost." There is an urgent need to develop a new testing technology to establish a direct correlation between blade shape deviation and performance deviation, thereby achieving accurate and rapid testing of compressor impeller performance consistency. Summary of the Invention

[0007] The technical problem to be solved by the present invention is to provide a method for detecting the efficiency consistency of compressor impellers in order to address the above-mentioned shortcomings. By establishing a quantitative correlation model between blade shape deviation and efficiency deviation, and combining color space conversion, the deviation can be accurately extracted, and finally, efficient and high-precision performance consistency detection can be achieved.

[0008] To solve the above technical problems, the present invention adopts the following technical solution: A method for detecting the consistency of compressor impeller efficiency includes the following steps: Step 1, Sample preprocessing and point cloud scanning; Step 2: Color space conversion. Convert the RGB color space to the HSV color space and utilize the characteristics of the chromaticity H dimension to achieve accurate quantitative extraction of leaf shape deviation, while removing background interference. Step 3: Deviation parameter calculation. The average deviation is calculated to quantify the overall level of the leaf shape deviation, the maximum deviation is calculated to quantify the extreme degree of the leaf shape deviation, and the large deviation ratio is calculated to quantify the distribution density of the sensitive area of ​​the leaf shape deviation, thus constructing a comprehensive characterization system for the leaf shape deviation. Step 4, Efficiency Deviation Prediction: By establishing a quantitative correlation model between leaf shape deviation and efficiency deviation, geometric deviation data is transformed into intuitive performance deviation results. Step 5, performance consistency determination, based on efficiency deviation results.

[0009] Furthermore, the specific implementation process of step 1 is as follows: Clean the surface of the compressor impeller to be tested by wiping it with anhydrous ethanol to remove oil stains and sanding the blade roots and edges to remove burrs. Fix the impeller on the testing fixture and ensure that the parallelism between the impeller reference surface and the fixture surface is ≤0.01mm to avoid installation deviations affecting scanning accuracy. A high-precision 3D scanner was used to perform a full-surface scan on the pre-processed impeller: the scanning accuracy was set to be higher than 0.005 mm, the scanning point density was not less than 100 points / mm², and the scanning range covered the main blades and the shunt blades of the impeller; after the scan was completed, the 3D scanner automatically generated 3D point cloud data of the blade shape deviation; the color image representing the blade shape deviation was adjusted so that the leading edge tip and trailing edge tip of the blade were placed as far left and far right as possible in the color area, and then the RGB image was exported.

[0010] Furthermore, the specific implementation process of step 2 includes the following steps: Step 21, data preprocessing, normalization of true color data; ; ; ; Step 22, Calculation of key parameters; Calculate the maximum, minimum, and difference values ​​in the RGB values; Maximum value: ; Minimum value: ; Difference: ; Calculate brightness: ; Calculate saturation: ; Calculate chromaticity: ; Color correction value: , ; Step 23, Background recognition and labeling: Separate the background and remove background interference; The white background of the image is adjusted in chroma. H All dimensions are marked as " NaN By defining parameters after color space conversion, background interference is eliminated, and only leaf shape deviation data of the effective area of ​​the leaf is retained: .

[0011] Furthermore, the specific implementation process of step 2 also includes the following steps: Step 24, Mapping of chromaticity correction values ​​to leaf shape deviation: ; in S d This represents the maximum downward deviation of the blade. S m To determine the maximum upper deviation of the leaf, the if branch in the formula removes invalid background data and checks if the chroma of the current pixel is "". NaN If it is NaN This indicates that the pixel is a white background area, which has no valid leaf shape deviation data and will not participate in subsequent calculations. The else branch converts the color of the valid leaf area into a specific leaf shape deviation. The current pixel is not... NaN This pixel belongs to the effective area of ​​the leaf; When inspecting the main blade, the main blade shape deviation ,in S Md The maximum downward deviation of the main blade. S Mu The maximum upper deviation of the main blade; When inspecting the splitter blades, the blade shape deviation... ,in S Sd The maximum downward deviation of the main blade.S Su The maximum upper deviation of the main blade; Therefore, the matrix for the blade shape deviation of the main blade and the spur blade is: Main leaf shape deviation and color mapping The chromaticity values ​​of the main blade region (0°-360°) are proportionally mapped to the deviation range of the main blade. S Md arrive S Mu The chromaticity values ​​of the leaf area are converted into specific physical deviation data, among which... M Corrected chromaticity values ​​for the main leaf area; Split vane shape deviation and color mapping The chromaticity values ​​of the splitter blade region (0°-360°) are proportionally mapped to the deviation range of the splitter blades. S Sd arrive S Su The chromaticity values ​​of the leaf area are converted into specific physical deviation data, among which... S This is the corrected chromaticity value for the splitter blade region.

[0012] Furthermore, the specific implementation process of step 3 also includes the following steps: Step 31: After the leaf shape deviation mapping is completed, determine the following parameters: Impeller diameter: Enter the impeller diameter according to the impeller model to be tested. d ; The blade consists of a leading edge region and a trailing edge region. The leading edge region is the area 10% of the chord length from the blade inlet end; the trailing edge region is the area 10% of the chord length from the blade outlet end. Large Deviation Identification Criteria: ; Step 32, eliminate the influence of impeller size on average blade profile deviation: count the deviation values ​​of all effective scanning points of the main blade and the splitter blade, and calculate the arithmetic mean; Step 33, Maximum Deviation Calculation: Extract the maximum deviation value from the effective scanning points of the main blade and the splitter blade, eliminate the influence of the impeller size, and quantify the degree of the most serious blade shape deviation in the compressor impeller; ; Step 34: In order to quantify the distribution density of severe blade shape deviation in the sensitive area in the leading edge region and trailing edge region, and thus reflect the potential impact of blade shape deviation on aerodynamic performance, it is necessary to calculate the maximum blade deviation ratio, that is, to calculate the ratio of the number of large deviation pixels in the region that exceed the set threshold to the total number of effective pixels in the region. Step 35, Output Results: Output the key parameters that characterize the consistency deviation of the blade shape performance obtained from the calculation.

[0013] Furthermore, the specific process of implementing step 32 is as follows: Count the number of valid pixels in the corresponding leaf-shaped deviation image: All leaf shape deviations in the image S ij Not for NaN The number of pixels is the total number of effective pixels. Num Where i and j are the row and column indices of the image pixels, i being the row number and j the column number, used to locate a specific pixel in the image. m , n These are the total number of rows and columns of the image, representing the pixel scale of the entire scanned image. S ij : Leaf shape deviation corresponding to the pixel in row i and column j; Calculate the sum of leaf shape deviations in the image: The sum of leaf shape deviations of all pixels in the effective area of ​​the main leaf in the image is calculated. The average blade profile deviation, after eliminating the influence of impeller size by combining the main blades and the splitter blades, is as follows: , where Δ S sum_M Total leaf shape deviation of the main leaf, Δ S sum_S This represents the total blade profile deviation of the splitter blades. Num m The effective number of pixels for the main blade. Num s This represents the effective number of pixels for the splitter blades.

[0014] Furthermore, the specific process of implementing step 34 is as follows: Obtain the column indices of the elements corresponding to the leaf shape deviation in the matrix, and construct a vector matrix of the effective pixels of the leaf: ,in J k It is the first in the matrix k One element; statistics J Number of elements in the vector: ; Identify the row coordinate vector matrix corresponding to the pixels at the leading and trailing edges of the blade: The row coordinate vector matrix of the leading edge region, starting with the effective region matrix of the blade. J In the middle, the pixels in the leading edge region with the largest absolute value of deviation are selected, and a matrix is ​​used. J top10% Mark the locations of these extreme deviations: ,inJ r It is a matrix J top10% The first in r One element, k It is a matrix J Element index, l It is the number of pixels with the maximum deviation in the leading edge region. r It is a matrix J top10% The element index, ranging from 1 to... l ; The row coordinate vector matrix of the trailing edge region, starting with the effective area matrix of the blade. J In the process, the pixels in the trailing edge region with the smallest absolute value of deviation are selected; a matrix is ​​then used. J bottom10% Mark the locations of these optimal regions, focusing on the areas where the leaf shape most closely resembles the standard model: ,in J s It is a matrix J bottom10% The first in s One element, s It is a matrix J bottom10% The element index, the range of values ​​is N - l +1 to N ; Column indexes for the leading and trailing regions: ; Modulus of column index: ; A matrix is ​​constructed by extracting all deviation values ​​that satisfy the large deviation threshold q from the pixel deviations of the leading and trailing edge regions of the blade. S locat_large : ,in S ij_large It is the leaf shape deviation value corresponding to the pixel in the i-th row and j-th column within the leaf's sensitive area, and this deviation value is ≥ q. The first one, the second one, the third one... t Leaf shape deviation value corresponding to each pixel in the sensitive area; Recognition Matrix S locat_large Elements that satisfy the large deviation condition for leaf shape are marked as 1, and elements that do not satisfy the large deviation condition are marked as 0. ; Statistical matrix S locat_large The number of pixels that satisfy the large deviation condition: ; Calculate the percentage of pixels in the leading and trailing edge regions of the blade that satisfy the large deviation condition relative to the total number of pixels in the corresponding leaf shape deviation image. .

[0015] Furthermore, the specific implementation process of step 4 is as follows: Leaf shape efficiency deviation is the result of the superposition of average leaf shape deviation and large leaf shape deviation; therefore, the formula for calculating deviation efficiency can be written as follows: ; in Represents the average leaf shape deviation Contribution to efficiency impact This represents the maximum leaf shape deviation. and the proportion of large deviations The overall contribution to efficiency.

[0016] Furthermore, the specific implementation process of step 4 is as follows: The fitted correlation for the efficiency deviation is: ; Among them, the fitting coefficient A =0.2537, basic deviation B =-0.0018; by When fitting to the form of independent variables, during adjustment a and b After the gain coefficient, it can be compared with... To maintain a good linear relationship, therefore The fitting correlation is: ; Where the fitting coefficient C =0.0112, additional loss base value D =-0.00002, the gain coefficient of the maximum deviation minus the average deviation. a =0.25, the gain coefficient for the large deviation proportion b =0.15.

[0017] The overall blade efficiency prediction bias is: ; .

[0018] Furthermore, the specific implementation process of step 5 is as follows: If: the impeller performance is determined to be consistent and meets the requirements; If the impeller performance is determined to be inconsistent, it is considered a defective product and must be returned to the production stage for reprocessing or scrapped.

[0019] The present invention adopts the above technical solution and has the following technical effects compared with the prior art: 1. High detection accuracy: The leaf shape deviation is accurately quantified through RGB-HSV color space conversion. Combined with the deviation correlation model, the performance deviation prediction accuracy is high, meeting the efficiency deviation detection requirement of ±0.5%. 2. High testing efficiency: The total testing time for a single impeller is ≤15 minutes (including scanning, calculation, and judgment), which is 8 times more efficient than bench testing (≥2 hours), making it suitable for mass production needs; 3. Low testing cost: No need to build an expensive turbocharger test bench (test bench cost ≥ 8 million yuan), only a conventional 3D scanner and computer are needed, reducing equipment investment costs by more than 90%; 4. Strong standard rationality: Establish a direct correlation between blade shape deviation (average deviation, maximum deviation, deviation ratio) and performance deviation, avoiding the problems of "geometric compliance but performance exceeding tolerance" and "geometric non-compliance but performance within tolerance", making the testing standard more in line with actual performance requirements. Attached Figure Description

[0020] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. In all the drawings, similar elements or parts are generally identified by similar reference numerals. In the drawings, the elements or parts are not necessarily drawn to scale.

[0021] Figure 1 This is a flowchart of the detection method in this invention; Figure 2 This is a schematic diagram illustrating the conversion relationship between the RGB color space and the HSV color space in this invention; Figure 3 This is a fitting formula for the efficiency deviation caused by the average leaf shape deviation in this invention. Figure 4 This is a fitting formula for the efficiency deviation caused by the maximum leaf shape deviation and the proportion of large deviations in this invention. Detailed Implementation

[0022] Examples, such as Figure 1 As shown, a method for detecting the consistency of compressor impeller efficiency includes the following steps: Step 1, Sample preprocessing and point cloud scanning; Clean the surface of the compressor impeller to be tested by wiping it with anhydrous ethanol to remove oil stains, and polish it with sandpaper (grit ≥ 800 mesh) to remove burrs from the blade root and blade edge. Fix the impeller on the testing fixture and ensure that the parallelism between the impeller reference surface (such as the hub end face) and the plane of the fixture is ≤ 0.01mm to avoid installation deviations affecting scanning accuracy.

[0023] A high-precision 3D scanner was used to perform a full-surface scan of the pre-processed impeller: the scanning accuracy was set to be higher than 0.005 mm, the scanning point density to be no less than 100 points / mm², and the scanning range covered the main blades and branch blades of the impeller; after scanning, the 3D scanner automatically generated 3D point cloud data of blade shape deviation; the color images representing blade shape deviation were adjusted so that the leading edge and trailing edge tips of the blades were placed as far left and right as possible in the colored area, and then RGB images were exported (image resolution). Different colors in the image represent blade shape deviation in different areas, and the maximum lower deviation of the main blade is recorded as . S Md Maximum upper deviation of the main blade S Mu Maximum downward deviation of the splitter blade S Sd Maximum upper deviation of the splitter blade S Su .

[0024] Step 2: Color space conversion. The RGB color space is converted to the HSV color space, and the characteristics of the chromaticity H dimension are used to achieve accurate quantitative extraction of leaf shape deviation. At the same time, background interference is removed to provide clean and effective data support for subsequent calculation of average deviation, maximum deviation and deviation ratio. like Figure 2 As shown, since the point cloud data obtained by the 3D scanner is a true-color (RGB) image that reflects leaf shape deviation through color, the RGB color space is a three-dimensional matrix (red, green, and blue channels), meaning that each pixel has three corresponding values. Its mathematical expression is: ; This means that it is impossible to use the RGB color space to statistically analyze one-dimensional leaf shape deviation data, and therefore impossible to obtain statistical data on average deviation, maximum deviation, and deviation ratio from true color images.

[0025] This invention converts the RGB color space to the HSV color space (chromaticity). H saturation S ,brightness V ), due to chromaticity H Dimensions possess full color attributes, therefore they can be utilized in the transformation color space. H Dimensional data is statistically analyzed to correspond to the statistics of leaf shape deviation.

[0026] The HSV color space expression is as follows: ; Among them, chromaticity H : Indicates the type of color (such as red, yellow, blue, green, cyan, etc.). It is represented by an angle (0°-360°), forming a color wheel. Within the 360° range, each degree represents a color.

[0027] brightness V : Indicates the brightness or darkness of a color. That is, the degree to which a color is diluted by white light. The lower the brightness, the closer the color is to black (0% is black, 100% is the brightest state of the color).

[0028] Saturation S : Indicates the purity or vividness of a color.

[0029] The color space conversion process is as follows: Step 21, data preprocessing, normalization of true color data; ; ; .

[0030] Step 22, Calculation of key parameters: Calculate the maximum, minimum, and difference values ​​in the RGB values; Maximum value: ; Minimum value: ; Difference: ; Calculate brightness: ; Calculate saturation: ; Calculate chromaticity: ; Color correction value: , .

[0031] Step 23, Background recognition and labeling: Separate the background and remove background interference; The white background of the image is adjusted in chroma. H All dimensions are marked as " NaN By defining parameters after color space conversion, background interference is eliminated, and only leaf shape deviation data of the effective area of ​​the leaf is retained: ; The saturation in the HSV space is set using the above formula. S ,brightness VThe threshold is used to filter out all pixels that meet the white characteristic, and then the chromaticity of these pixels is... H Marked as NaN This achieves separation between the background and the effective area of ​​the blade.

[0032] Step 24, Mapping of chromaticity correction values ​​to leaf shape deviation: ; Where S d This represents the maximum downward deviation of the blade. S m To determine the maximum upper deviation of the leaf, the if branch in the formula removes invalid background data and checks if the chroma of the current pixel is "". NaN If it is NaN This indicates that the pixel is a white background area, which has no valid leaf shape deviation data and will not participate in subsequent calculations. The else branch converts the color of the valid leaf area into a specific leaf shape deviation. The current pixel is not... NaN This pixel belongs to the effective area of ​​the leaf.

[0033] When inspecting the main blade, the main blade shape deviation S Md The maximum downward deviation of the main blade. S Mu The maximum upper deviation of the main blade; When inspecting the splitter blades, the blade shape deviation... ,in S Sd The maximum downward deviation of the main blade. S Su The maximum upper deviation of the main blade; Therefore, the matrix for the blade shape deviation of the main blade and the spur blade is: Main leaf shape deviation and color mapping The chromaticity values ​​of the main blade region (0°-360°) are proportionally mapped to the deviation range of the main blade. S Md arrive S Mu The chromaticity values ​​of the leaf area are converted into specific physical deviation data, among which... M Corrected chromaticity values ​​for the main leaf area; Split vane shape deviation and color mapping The chromaticity values ​​of the splitter blade region (0°-360°) are proportionally mapped to the deviation range of the splitter blades. S Sd arrive S Su The chromaticity values ​​of the leaf area are converted into specific physical deviation data, among which... S This is the corrected chromaticity value for the splitter blade region.

[0034] Step 3, Deviation Parameter Calculation: The average deviation is calculated to quantify the overall level of blade shape deviation, the maximum deviation is calculated to quantify the extreme degree of blade shape deviation, and the large deviation ratio is calculated to quantify the distribution density of the sensitive area of ​​blade shape deviation. A comprehensive characterization system of blade shape deviation is constructed to provide key data support for the subsequent establishment of a correlation model between blade shape deviation and efficiency deviation and accurate prediction of impeller efficiency, while avoiding the limitation of focusing only on the value of a single geometric deviation. Step 31: After the leaf shape deviation mapping is completed, determine the following parameters: Impeller diameter: Enter the impeller diameter according to the impeller model to be tested. d ; Leading edge / trailing edge definition range: The blade includes a leading edge region and a trailing edge region. The leading edge region is the area 10% of the chord length from the blade inlet end. It is the initial action area where the airflow enters the blade. Blade shape deviations in this area (such as bulges and depressions) will directly change the initial flow direction of the airflow, easily causing airflow separation and turbulence loss. It is a highly sensitive area affecting compressor efficiency. The trailing edge region is the area 10% of the chord length from the blade outlet end. It is the final adjustment area where the airflow leaves the blade. Blade shape deviations in this area will affect the uniformity of the airflow, and thus affect the inlet conditions of the subsequent impeller stage. Large Deviation Identification Criteria: .

[0035] Step 32, eliminate the influence of impeller size on average blade profile deviation: count the deviation values ​​of all effective scanning points of the main blade and the splitter blade, and calculate the arithmetic mean; Count the number of valid pixels in the corresponding leaf-shaped deviation image: All leaf shape deviations in the image S ij Not for NaN The number of pixels is the total number of effective pixels. Num ,in i , j It is the row and column index of the image pixels. i It's a line number. j These are column numbers, used to locate a specific pixel in the image. m , n These are the total number of rows and columns of the image, representing the pixel scale of the entire scanned image. S ij : No. i line, number j Leaf shape deviation corresponding to column pixels.

[0036] Calculate the sum of leaf shape deviations in the image: The sum of leaf shape deviations of all pixels in the effective area of ​​the main leaf in the image is calculated.

[0037] The average blade profile deviation, after eliminating the influence of impeller size by combining the main blades and the splitter blades, is as follows: , where Δ S sum_M Total leaf shape deviation of the main leaf, Δ S sum_S This represents the total blade profile deviation of the splitter blades. Num m The effective number of pixels for the main blade. Num s This represents the effective number of pixels for the splitter blades.

[0038] Step 33, Maximum Deviation Calculation: Extract the maximum deviation value from the effective scanning points of the main blade and the splitter blade, eliminate the influence of the impeller size, and quantify the degree of the most serious blade shape deviation in the compressor impeller; .

[0039] Step 34: In order to quantify the distribution density of severe blade shape deviations in the sensitive areas within the leading and trailing edge regions, and thus reflect the potential impact of blade shape deviations on aerodynamic performance, it is necessary to calculate the maximum blade deviation ratio. This is the ratio of the number of pixels with large deviations exceeding a set threshold in the region to the total number of effective pixels in the region. The specific calculation process is as follows: Obtain the column indices of the elements corresponding to the leaf shape deviation in the matrix, and construct a vector matrix of the effective pixels of the leaf: This process clearly defines the pixel location range of the effective area of ​​the blade, allowing subsequent deviation statistical calculations to accurately focus on the blade area and eliminate background interference. J k It is the first in the matrix k Each element.

[0040] statistics J Number of elements in the vector: .

[0041] Identify the row coordinate vector matrix corresponding to the pixels at the leading and trailing edges of the blade: The row coordinate vector matrix of the leading edge region, starting with the effective region matrix of the blade. J In the process, the leading edge region pixels with the largest absolute deviation are selected; a matrix is ​​then used. J top10% By marking the locations of these extreme deviations and focusing on the areas with the most severe blade shape deviations, a precise pixel range is provided for analyzing the impact of local extreme deviations on impeller efficiency. ,in J r It is a matrix J top10%The first in r One element, k It is a matrix J Element index, l It is the number of pixels with the maximum deviation in the leading edge region. r It is a matrix J top10% The element index, ranging from 1 to... l .

[0042] The row coordinate vector matrix of the trailing edge region, starting with the effective area matrix of the blade. J In the process, the pixels in the trailing edge region with the smallest absolute value of deviation are selected; a matrix is ​​then used. J bottom10% Mark the locations of these optimal regions, focusing on the areas where the leaf shape most closely resembles the standard model, to provide a precise pixel range for comparing the performance differences between optimal and extreme deviation regions: ,in J s It is a matrix J bottom10% The first in s One element, s It is a matrix J bottom10% The element index, the range of values ​​is N - l +1 to N .

[0043] Column indexes for the leading and trailing regions: ; Modulus of column index: ; Extract all pixels that meet the large deviation threshold from the pixel deviations of the leading and trailing edge regions of the blade. q The deviation values ​​form a matrix S locat_large : ,in S ij_large It is the first in the leaf sensitive area i line, number j The leaf shape deviation value corresponding to the column pixel, and this deviation value is ≥ q. The first one, the second one, the third one... t The leaf shape deviation value corresponding to each pixel in the sensitive area.

[0044] Recognition Matrix S locat_large The element that satisfies the condition for large deviation in leaf shape is denoted as . 1 Elements that do not meet the large deviation condition are marked as 0: ; Statistical matrix S locat_large The number of pixels that satisfy the large deviation condition: ; Calculate the percentage of pixels in the leading and trailing edge regions of the blade that satisfy the large deviation condition relative to the total number of pixels in the corresponding leaf shape deviation image. .

[0045] Step 35, Output Results: Calculate the key parameters characterizing the consistency deviation of blade shape performance: , , .

[0046] Step 4, Efficiency Deviation Prediction: By establishing a quantitative correlation model between leaf shape deviation and efficiency deviation, geometric deviation data is transformed into intuitive performance deviation results, providing accurate and quantifiable basis for subsequent performance consistency judgment. Leaf shape efficiency deviation is the result of the superposition of average leaf shape deviation and large leaf shape deviation; therefore, the formula for calculating deviation efficiency can be written as follows: ; in Represents the average leaf shape deviation Contribution to efficiency impact This represents the maximum leaf shape deviation. and the proportion of large deviations The overall contribution to efficiency.

[0047] Fitting was performed based on a large amount of CFD simulation and experimental data. The impact on efficiency is almost linear, such as Figure 3 As shown, therefore The fitted correlation for the efficiency deviation is: ; Among them, the fitting coefficient A= 0.2537, basic deviation B =-0.0018; After numerous attempts to fit the formula, it was found that... When fitting to the form of independent variables, during adjustment a and b After the gain coefficient, it can be compared with... Maintain a good linear relationship, such as Figure 4 As shown, therefore The fitting correlation is: ; Where the fitting coefficientC =0.0112, additional loss base value D =-0.00002, the gain coefficient of the maximum deviation minus the average deviation. a =0.25, the gain coefficient for the large deviation proportion b =0.15.

[0048] The overall blade efficiency prediction bias is: ; .

[0049] Step 5, performance consistency determination; Based on the efficiency deviation results: If: the impeller performance is determined to be consistent and meets the requirements; If the impeller performance is determined to be inconsistent, it is considered a defective product and must be returned to the production stage for reprocessing or scrapped.

[0050] The description of this invention is given for illustrative and descriptive purposes only and is not intended to be exhaustive or to limit the invention to the forms disclosed. Many modifications and variations will be apparent to those skilled in the art. The embodiments were chosen and described in order to better illustrate the principles and practical application of the invention and to enable those skilled in the art to understand the invention and design various embodiments with various modifications suitable for a particular purpose.

Claims

1. A method for detecting uniformity of efficiency of a compressor impeller, characterized by: The method comprises the following steps: Step 1, sample pretreatment and point cloud scanning; Step 2: color space conversion, converting the RGB color space into the HSV color space and using the characteristics of the hue H dimension to realize accurate quantification of the leaf shape deviation and eliminate background interference; Step 3, deviation parameter calculation, quantifying the overall level of the leaf shape deviation by calculating the average deviation, quantifying the extreme degree of the leaf shape deviation by calculating the maximum deviation, and quantifying the sensitive area distribution density of the leaf shape deviation by calculating the large deviation ratio, so as to construct a comprehensive characterization system of the leaf shape deviation; Step 4, efficiency deviation prediction, converting the geometric deviation data into intuitive performance deviation results by establishing a quantitative correlation model between the leaf shape deviation and the efficiency deviation; Step 5, performance consistency determination, according to the efficiency deviation results.

2. The method of claim 1, wherein: The specific implementation process of step 1 is as follows: The surface of the compressor impeller to be detected is cleaned, and the burrs on the blade root and blade edge are removed by wiping with anhydrous ethanol and sanding; the impeller is fixed on the detection tooling table to ensure that the parallelism of the impeller reference surface and the tooling table plane is ≤0.01 mm, so as to avoid the influence of installation deviation on the scanning accuracy; The pretreated impeller is scanned by a high-precision three-dimensional scanner: the scanning accuracy is set to be higher than 0.005 mm, the scanning point density is not less than 100 points / mm², and the scanning range covers the main blades and splitter blades of the impeller; after scanning, the three-dimensional scanner automatically generates three-dimensional point cloud data of the leaf shape deviation; the color picture representing the leaf shape deviation is adjusted so that the blade leading edge tip and the trailing edge tip are placed at the leftmost and rightmost of the color area respectively, and then the RGB picture is exported.

3. The method for detecting the consistency of compressor impeller efficiency as described in claim 1, characterized in that: The specific implementation process of step 2 comprises the following steps: Step 21, data preprocessing, true color data normalization; ； ； ； Step 22, key parameter calculation; Calculate the maximum value, minimum value and difference value in RGB; Maximum value: ; Minimum value: ; Difference: ; Computing luminance: ; Computing saturation: ; Calculated chrominance: ; Chromaticity correction value: , ; Step 23, background recognition and marking, separate the background and eliminate background interference; The white color of the picture background is marked as H " in the chroma dimension NaN The leaf shape deviation data of the effective area of the leaf blade is retained by excluding the background interference through the parameter definition after the color space conversion. 。 4. The method for detecting the consistency of compressor impeller efficiency as described in claim 3, characterized in that: The specific implementation process of step 2 further comprises the following steps: Step 24, mapping of chroma correction value and leaf shape deviation: ； wherein S d is the maximum lower deviation of the leaf, S m is the maximum upper deviation of the leaf, the if branch in the formula eliminates the background invalid data, and judges whether the chrominance of the current pixel is NaN , if yes NaN , the pixel represents a white background area, and the area does not participate in subsequent calculation because it has no valid leaf deviation data, and the else branch converts the color of the leaf effective area into a specific leaf deviation, and the current pixel is not NaN , the pixel belongs to the leaf effective area; when detecting the main blade, the main blade blade shape deviation wherein S Md is the maximum lower deviation of the main blade, S Mu is the maximum upper deviation of the main blade; When detecting splitter vanes, splitter vane profile deviation wherein S Sd is the maximum lower deviation of the main vane, S Su is the maximum upper deviation of the main vane; Therefore, the matrix of the main blade and the splitter blade leaf shape deviation is: Main blade camber and chroma mapping The chroma values 0°-360° of the main blade region are proportionally mapped to the camber range of the main blade S Md to S Mu The chroma values of the blade region are converted to specific physical camber data, where M is the corrected chroma value of the main blade region; Split vane leaf shape deviation and chroma mapping The chroma values 0°-360° of the split vane region are proportionally mapped to the deviation range of the split vane S Sd to S Su The chroma values of the vane region are converted to specific physical deviation data, wherein S is the corrected chroma value of the split vane region.

5. The method for detecting the consistency of compressor impeller efficiency as described in claim 1, characterized in that: The specific implementation process of step 3 further comprises the following steps: Step 31, after the leaf shape deviation mapping is completed, the following parameters are determined: Impeller wheel diameter: input impeller wheel diameter according to the type of impeller to be detected d ; The blade includes a leading edge region and a trailing edge region, the leading edge region is a region with a chord length of 10% from the blade inlet end, and the trailing edge region is a region with a chord length of 10% from the blade outlet end; Large deviation recognition criterion: ; Step 32, average leaf shape deviation eliminating the influence of impeller size: statistics the deviation values of all effective scanning points of the main blade and the splitter blade, and calculate the arithmetic mean; Step 33, maximum deviation calculation: extract the maximum deviation value in the effective scanning points of the main blade and the splitter blade, and eliminate the influence of the impeller size to quantify the most serious leaf shape deviation in the compressor impeller; ； Step 34, in order to quantify the distribution density of serious leaf shape deviation in the sensitive area, so as to reflect the potential influence degree of leaf shape deviation on aerodynamic performance, the maximum deviation ratio of the blade needs to be calculated, that is, the ratio of the number of large deviation pixels in the area exceeding the set threshold to the total number of effective pixels in the area is calculated; Step 35, output the result, output the key parameters calculated to characterize the consistency deviation of leaf shape performance.

6. The method of claim 5, wherein: The specific process implemented by the step 32 is as follows: Count the number of valid pixels in the corresponding leaf-shaped deviation image: All leaf shape deviations in the image S ij Not for NaN The number of pixels is the total number of effective pixels. Num ,in i , j It is the row and column index of the image pixels. i It's a line number. j These are column numbers, used to locate a specific pixel in the image. m , n These are the total number of rows and columns of the image, representing the pixel scale of the entire scanned image. S ij : No. i line, number j Leaf shape deviation corresponding to column pixels; Sum of leaf deviation in the picture is calculated: , the sum of leaf deviation of all pixels in the effective area of the main blade in the picture is counted; The average leaf shape deviation of the comprehensive main blade and the splitter blade eliminates the influence of the size of the impeller: where Δ S sum_M is the total blade shape deviation of the main blade, Δ S sum_S is the total blade shape deviation of the splitter blade, Num m is the number of effective pixels of the main blade, Num s is the number of effective pixels of the splitter blade.

7. The method for detecting the consistency of compressor impeller efficiency as described in claim 5, characterized in that: The specific process implemented by the step 34 is as follows: Obtain the matrix, the column number of the elements corresponding to the leaf shape deviation, and form the vector matrix of the effective pixels of the blade: ,in J k It is the first in the matrix k One element; Statistics J Number of elements in a vector: ; Identify the row coordinate vector matrix of the pixels corresponding to the leading edge region and the trailing edge region of the blade: The row coordinate vector matrix of the leading edge region is first derived from the blade effective region matrix J In this way, the pixel of the leading edge region with the largest absolute value of deviation is screened out, and the matrix J top10% Mark the positions of these extreme deviations: wherein J r is the matrix J top10% is the element r of the matrix k is the element number of the matrix J is the maximum deviation pixel number of the leading edge region, l is the element number of the matrix r J top10% is the element number of the matrix l ;​ The row coordinate vector matrix of the trailing edge region is first obtained from the blade effective region matrix J The pixel of the trailing edge region with the minimum absolute deviation is screened out, and the matrix J bottom10% The positions of these optimal regions are marked, and the region closest to the standard model in focus is determined. wherein J s is a matrix J bottom10% is the s th element of the matrix s is the element number of the matrix J bottom10% , which takes the value range of N - l +1 to N ; Row indices of the leading edge region and the trailing edge region: ; Modulo of column index: ; Selecting all the deviation values satisfying the large deviation threshold value from the pixel deviation of the leading edge region and the trailing edge region of the blade to constitute a matrix q S locat_large :​ wherein S ij_large is the leaf shape deviation value corresponding to the pixel pair of the i-th row and the j-th column, and the deviation value ≥ i , j the leaf shape deviation value corresponding to the pixel pair of the i-th row and the j-th column; and q , the leaf shape deviation value corresponding to the pixel pair of the i-th row and the j-th column; and t the leaf shape deviation value corresponding to the pixel pair of the i-th row and the j-th column; and Identifying matrices S locat_large The elements satisfying the large deviation condition with the leaf shape are denoted by 1 The elements not satisfying the large deviation condition are denoted by 0. ； Statistical matrix S locat_large Number of pixels satisfying the large deviation condition in the middle ； Calculate the number of pixels in the leading edge region and the trailing edge region of the blade that meet the large deviation condition, that is, calculate the proportion of large deviation in the number of pixels in the picture corresponding to the leaf shape deviation picture: 。 8. The method for detecting the consistency of compressor impeller efficiency as described in claim 7, characterized in that: The specific implementation process of the step 4 is as follows: The leaf shape efficiency deviation is the result of the superposition of the average leaf shape deviation and the large leaf shape deviation, so the deviation efficiency calculation formula can be written as follows: ； wherein represents the average leaf shape deviation the contribution to the efficiency impact, represents the maximum leaf shape deviation and the fraction of large deviations the overall contribution to the efficiency impact.

9. The method for detecting the consistency of compressor impeller efficiency as described in claim 8, characterized in that: The specific implementation process of the step 4 is as follows: The fitted correlation for the efficiency deviation is: ； wherein the fitting coefficients A = 0.2537, the base bias B = -0.0018; The fitting is performed in the form of independent variable Adjusting the gain coefficients of a and b , a good linear relationship with can be maintained, so the fitting correlation of is ； where the fitting coefficients C = 0.0112, the additional loss base value D = -0.00002, the gain coefficient of the maximum deviation minus the average deviation a = 0.25, the gain coefficient of the large deviation proportion b = 0.15; The total blade efficiency prediction deviation is: ； 。 10. The method for detecting the consistency of compressor impeller efficiency as described in claim 9, characterized in that: The specific implementation process of the step 5 is as follows: If: determine that the performance of the impeller is consistent and meets the requirements; If: determine that the performance of the impeller is inconsistent and is unqualified, it needs to be returned to the production link for reprocessing or scrapped.

Images