# Fatigue Strength Evaluation for Remanufacturing Impeller of Centrifugal Compressor Based on Modified Grey Relational Model

Authors: Qingchao Sun; Bowen Shi; Xiaokai Mu; Kepeng Sun
Journal: Advances in Materials Science and Engineering
 (2020)

---

## Abstract

The fatigue strength, as the essential basis of residual life evaluation, is required to be obtained timely for remanufacturing. Since impeller damage is characterized with very-high-cycle fatigue (VHCF), it is difficult to directly test the strength data. The transformation method of multisource strength data is proposed to predict fatigue strength for impeller based on grey relational theory. The multisource strength data, as factor space, primarily include available existing experimental data and operating data, while the strength data of the remanufacturing impeller are taken as target data. The fatigue strength model of material and component are presented to analyze the influence factors of remanufacturing target strength. And similar material provides a theoretical basis for selecting reference data reasonably. Considering the correlation and difference between available data and target data, the grey relational function is established, and the correction function of the target residual is brought forward to reduce the transformation deviation. The entropy-weight theory is implemented to determine the different impacts of multisource data on target strength. A test case, predicting the unknown impeller fatigue strength with various impellers, is applied to validate the proposed transformation method, and the results show that the predicted strength data are consistent with the experimental data well.


---

## Body

## 1. Introduction

Remanufacturing is one of the effective ways to solve resource shortage and environmental protection problem restricting to manufacturing development. It can save 50% in cost, 60% in energy, and 70% in material compared to new products [1]. Remanufacturability assessment, aiming at estimating whether products or parts are worth to be remanufactured, is the first step towards remanufacturing. The key issue of remanufacturability assessment is to evaluate residual life of products for ensuring safety in service in the next life cycle, where the fatigue strength is the essential basis for residual life evaluation [2, 3].Centrifugal compressor, as a kind of high value-added product, is widely used in petrochemical industry, natural gas, and so on. Since the centrifugal compressor is generally produced with one piece, it has diversity in material, structure, and working condition. Moreover, the service time of a centrifugal compressor is long, up to 20–30 years, the fatigue damage typically belongs to high-cycle fatigue (HCF) or very high-cycle fatigue (VHCF) [4, 5]. The fatigue tests of HCF or VHCF for the centrifugal compressor impeller is very complicated and time-consuming. However, in order to reduce losses caused by compressor breakdown, the fatigue strength is required to be obtained timely for remanufacturability assessment.Due to considerable divergence of the fatigue strength data for the centrifugal compressor impeller, it is difficult to directly obtain precise fatigue strength data based on existing theoretical models. Garwood and Nishijima et al. [6–8] found that, for low- or medium-strength steel with Vickers hardness (HV) less than 400, there is a linear relationship between rotating-bending or tension-compression fatigue strength σw,MPa and tensile strength Rm or Vickers hardness, as the fatigue fracture is mainly caused by surface damage. However, the linear relationship is not suitable to the material with higher hardness or tensile strength, which shows obvious divergence [5, 9]. The reason for this phenomenon is that the damage mechanism of the centrifugal compressor impeller includes both surface damage and internal damage. Surface damage, the primary form of which is fatigue crack, mainly starts from the discontinuous part of geometry and material [6]. Internal damage under low stress is mainly caused by internal defects, such as nonmetallic inclusion and the GBF (granular bright facet) area [3]. Therefore, the data transformation, based on the available existing strength data, is a significant way to obtain target strength data for centrifugal compressor impeller remanufacturing.The prerequisite of strength data transformation is to determine the impacts of material, structure, surface condition, and working condition on fatigue strength. In recent years, some researchers studied the influence of internal defect on fatigue strength and put forward several models based on micromechanics, such as the parallel layer model by Tanaka and Mura [10, 11] and the dislocation pileup model by Chang et al. [12, 13]. Among all models, the inclusion equivalent projected area model by Murakami et al. [14–16] is the most representative model. All these models provide theoretical basis to predict the fatigue strength of the centrifugal compressor, but it is quite difficult to transform the existing data to target data directly because high-strength steel is sensitive to small defects and nonmetallic inclusions [9] and the size and distribution of inclusion are random obviously [17–19].For fatigue strength of centrifugal compressor impeller, the available existing experimental data and actual operating data are typically small-sample discrete data. The statistical deduction of small-sample discrete data is used to obtain the transformed data by multisource data acquirement and fusion. The common data fusion algorithm includes the grey relational method, Kalman filter method, adaptive weighted method, D-S evidence theory, and neural network theory [20–24]. These algorithms were mainly applied to signal analysis and fault diagnosis at first and recently used in the assessment of fatigue strength and residual life [25–27]. The fatigue strength of the centrifugal compressor impeller is influenced by multiple impact factors, such as design parameters, material, and structure. Since it is difficult to obtain all of the data on the impact factors, the existing available data of high-strength steel are often incomplete and scattered. Grey relational analysis [28, 29], assessing related data with the grey relational grade of sample data, has great advantage in mapping relationship analysis of incomplete and scattered data, but the problems of qualitative analysis and lack of precision still exist at the same time [30]. So, the grey relational theory must be modified in data transformation.This paper proposes a data transformation method, according to available existing experimental data and actual operating data, to timely obtain fatigue strength data of the centrifugal compressor impeller. Firstly, impact factors are determined by analyzing the material and component fatigue strength characteristics of high-strength steel. Then, the transformation model of multisource strength data is presented based on the modified grey relational theory. The grey analysis function of factor space and target data is established; meanwhile, the residual modification on target data and the parameter correction and optimization ensure the rationality of transformation data.


## 2. Fatigue Strength Model

The prerequisite of establishing impact factor space and transforming fatigue strength data is to determine the impact factors of fatigue strength and the influence of each factor on fatigue strength. The fatigue strength models for material FV520B-I for the centrifugal compressor impeller are presented to analyze the impact factors of fatigue strength in two aspects: (1) the material fatigue strength model, which focuses on the influence of material hardness and yield strength; (2) the component fatigue strength model, which focuses on the influence of stress gradient, structure size, surface condition, working condition, etc.


### 2.1. Fatigue Strength Model of Material

Material fatigue strength is generally evaluated with material static performance. There are many impact factors on material fatigue strength; however, there is no an authoritative theory, with a quantitative and comprehensive analysis, containing all kinds of impact factors. Therefore, instead of all kinds of factors, the primary factors that have significant impacts on material fatigue strength include tensile strength, material hardness, and inclusion size [8, 14, 31].For low-strength steel with Vickers hardness (HV) less than 400, there is the linear relationship between rotating-bending or tension-compression fatigue strengthσw,MPa and material Vickers hardness [6–8]:(1)σw≈1.6HV±0.1HVHV≤400.There is also a linear relationship between fatigue strength and tensile strengthσb,MPa for structural steel with tensile strength less than 1200 MPa:(2)σw≈0.5σbσb≤1200MPa.The essence of fatigue strength is the microplastic deformation resistance of metal material, which is different from tensile strength in physical characteristics. Thus fatigue strength is only approximately estimated with tensile strength.The mechanical properties of material FV520B used in the centrifugal compressor impeller are similar to high-strength steel (tensile strengthσb>1200MPa). The above linear relationship between fatigue strength and tensile strength does not exist when tensile strength reaches up to 1250 MPa. Inclusion has been proved to be the key factor affecting fatigue strength of high-strength steel, and Murakami and Endo [32] proposed the expression of the impacts of steel matrix hardness HV and inclusion size Ain on material fatigue strength:(3)σw=1.56HV+120Ain1/6.The fatigue strength can also be determined by substituting inclusion sizeAin with the size of GBF area AGBF [33], and the size of GBF is formulated as follows:(4)AGBF=2HV+120σw6.The influence of stress ratior and hydrogen on fatigue strength of high-strength steel is also significant. Considering the influence of stress ratio and hydrogen, the relationship between material hardness, inclusion size, and fatigue strength can be expressed as [34](5)σw=2.7HV+12015/16Ain3/161−r2α,where α is the material parameter related to the material Vickers hardness, α=0.226+HV×10−4.


### 2.2. Material Similarity

Due to the great difference of fatigue property between different materials, the similar materials with similar fatigue strength are selected in data transformation, which ensures the effectiveness and accuracy of the transformed data. Generalized tangent modulus theory, estimating generality and similarity of different materials by the tangent modulus factorΦt, plays a great role in fatigue life prediction [35–37]. It can be used to determine fatigue property of unknown material according to the known material. Φt, a product factor of tangent modulus Et, is used to express the equilibrium state of elastomer. The tangent modulus factor curve σ¯−Φt can be expressed as(6)Φt=1σ¯dσ¯dε¯=E¯tσ¯=1σ¯EtE.Combined withS-N curve, the n¯−Φt curve can be obtained by the relationship of S=σ¯ and mean fatigue life n¯=N/106 as well as the relationship of S=σ¯ and Φt:(7)n¯=NN0=1−Φ1−Φ21−Φt−Φ21/s−11/t,where σ¯andn¯ are, respectively, the dimensionless values of S and N obtained from the S-N curve and σ¯−Φtandn¯−Φt curves are the dimensionless curve taking tangent modulus factor Φt as the parametric variable. Furthermore, the σ¯−n¯−Φt curve can be obtained by putting the two curves into a diagram. Figure 1 shows the σ¯−n¯−Φt curve for four kinds of material FV520B-I, X12Cr13-I, KMN-I, and FV520B-S. The σ¯Φt curve is derived from the stress-strain curve of material, while the n¯Φt curve is derived from the S-N curve of material. Thus, the σ¯−n¯−Φt curve reflects the working performance of material in the condition of static load and symmetrical cyclic load, respectively.Figure 1
σ
¯
−
n
¯
−
Φ
t curve of these four materials.As shown in Figure1, the above four kinds of materials have a high similarity for both σ¯−Φt curve and n¯−Φt curve, and especially when Φt≥1.4, the curves for these material gradually tend to be unanimous. So the four kinds of material can be regarded as similar material, and their fatigue properties are similar. Therefore, the obtained fatigue strength data are more scientific and reliable using these materials.


### 2.3. Fatigue Strength of Components

Material fatigue strength is often obtained by experiments, and the component fatigue strength is obtained by transforming material fatigue strength, taking into account the impacts of structure shape, stress gradient, surface condition, and temperature The relationship between material fatigue strength and component fatigue strength can be expressed as [38](8)σaf,C=σA,tsc⋅ftot,af,(9)ftot,af=1fST,af⋅fGR,af2−1.0+f1,af2fm,af⋅fTI,af⋅fTP,af⋅fGS,af,where σaf,C is the endurance stress limit of the local component, σA,tsc is the material alternating stress limit, ftot,af is the overall impact factor, fST,af presents the statistical impact factor, fGR,af is the stress gradient impact factor, f1,af is the combined surface roughness and degree of forging impact factor, fm,af is the mean stress impact factor, fTI,af is the thermal impact factor, fTP,af is the technological parameter impact factor, and fGS,af is the general surface factor.The statistical impact ist assumed as a log-normal distribution of the strength parameters:(10)fST,af=σaf,90σaf,i,where σaf,90 is the material alternating stress limit with a survival probability of 90% and σaf,i is the endurance stress limit value with the required survival probability i.The stress gradient impact factor on endurance stress limit of the local component can be expressed as below:(11)fGR,af=1+σA,b/σA,tsc−12/bv⋅χ′v,where χ′ is the relative stress gradient, σA,b is the material alternating stress limit for bending, σA,tsc is the material alternating stress limit for tension or compression, ν is the material parameter, and b is the diameter of the smooth sample under flexural load.If the material sample has a mean roughness depthRZ,s=1, the respective surface condition of a component acts fully on the endurance stress limit:(12)fSR,af=1−0.22lgRt,C0.64⋅lgσUTS+0.45⋅lgRt,C0.531−0.22⋅lgRt,M0.64⋅lgσUTS+0.45⋅lgRt,M0.53,where Rt,M is the maximum roughness depth of the material sample and Rt,C is the maximum roughness depth of the component.The mean stress impact on the endurance stress limit can be expressed as(13)fm,af=σA,tscσAσm,where σA,tsc is the material alternating stress limit and σAσm is the material alternating stress limit at the mean stress σm.If the temperature is greater than100°C, the appropriate impact factor is determined as(14)fTE,af=1−aT,af⋅10−3⋅T°C−100,where aT,af is the impact coefficient of temperature and T is the temperature in centigrade.Taking into account the differing strengths of materials as a function of the effective diameter of the semi or the unfinished casting, the type of material and the technological treatment, e.g., tempering, the technological parameter impact factor can be calculated as(15)fTP,af=1Kddeff,S⋅minKddeff,Kddeff,P,where deff is the effective diameter of the semi or the unfinished casting, deff,P represents deff up to which no technological parameter influence is considered, and deff,S represents deff according to the respective material standard.General surface factorfGS,af can be defined according to special surface properties, such as shot peening, rolling, carburizing, nitriding, induction hardening, and flame hardening, and these impact factors are primarily based on professional experience.


### 2.4. Grey Relational Model of Multisource Strength Data considering Incomplete Factor

The influence of design parameters, working conditions, and transmission modes on fatigue strength is smaller than material parameters and stress ratio, especially in the normal working mode. In this paper, design parameters, working conditions, and transmission modes are in the range of normal working parameters. Combined with the fatigue strength models above, the influence of these parameters on fatigue strength follows a roughly linear law.Grey relational analysis for data with linear laws can be used to obtain a relatively precise result by considering the correlation of source data in factor space and establishing grey relational functions of target data and source data [28, 39]. However, it is difficult to obtain the strength data of the centrifugal compressor impellor, thus the source data are often incomplete. In addition, source data, formed and obtained in different periods, have obvious multisource characteristic. Thus, it is necessary to correct the transformation result acquired from limited data.


### 2.5. Grey Relational Grade for Fatigue Strength Data

Based on the point relational coefficient model, the relational grade of fatigue strength can be determined, with taking the influence of HV, tensile strength, inclusion size, etc., on material fatigue strength, as well as the influence of structure shape, stress gradient, surface condition, mean stress, temperature, etc., on component fatigue strength into account.Supposing the measured arrayX0=x01,x02,…,x0n is the characteristic array, and the reference array is denoted as(16)Xi=xi1,xi2,...,xin.The relational coefficientαik between the characteristic array and the reference array is defined as [39](17)αik=mini=1mmaxk=1nx0k−xik+ξmaxi=1mmaxk=1nx0k−xikx0k−xik+ξmaxi=1mmaxk=1nx0k−xik,where ξ is the resolution coefficient and the value is in the range of 0～1.DenoteγX0,Xi with the following formula:(18)γX0,Xi=1n∑k=1nαik.IfγX0,Xi meets the axiom of normalization, integrality, duality, and proximity [40], γX0,Xi is defined as the grey relational grade between X0andXi.The grey relational analysis is an approximate analysis method. Since the influence of each factor on fatigue strength follows an approximately linear law and source data are always incomplete, the calculation result has some deviation compared with the real data. It is necessary to reduce the deviation by residual error modification [29].For any reference arrayXi=xi1,xi2,…,xin, if X0andXi have an increment Δxi, the variation of target value fi is denoted with Δfk and the variation of first data is Δf0. The modified residual error of target value is(19)βX0,Xi=εiΔf0,where(20)εi=1n∑k=1nΔfkΔf0.Furthermore, the modified transformation model based on grey relational grade can be established as follows:(21)fi′=γX0,Xifi+βX0,Xi.In equation (21), since the original target value is modified with residual error considering individual difference of data in factor space, data deviation caused by data transformation can be decreased compared with the original target value.


### 2.6. Weight Calculation of Multisource Strength Data

Generally speaking, the factors in factor space and target data have the same parameter name but have the difference in material, construction, or working condition data. Due to the multisource data in factor space, it is necessary to weigh the influence of different types of data on transformation result, to distinguish the influence grade of each factor. Entropy-weight theory has an advantage on simple calculation and little subjective interference to calculate weight coefficient [41, 42].Based on the comprehensive evaluation method [43, 44], the entropy-weight coefficient decision matrix of strength data for centrifugal compressor is established to express the correlation among three kinds of factors that include working condition, design parameter, and construction parameter. Thus, the decision matrix is(22)Qm×n=q11⋯q1n⋮⋮⋮qm1⋯qmn.Inm×n matrix, the number of influence factors are n and each factor is assessed m times independently. The treated decision matrix Qm×n reduces subjective interference to some extent.DefineXij as the membership degree of the jth factor, and QM=maxqij,Qm=minqij,i=1,2,…,m. Then the membership degree of each factor in Qm×n is(23)Xij=Qij−QmQM−Qm,positive factors,QM−QijQM−Qm,negative factors.Normalizing each factor of the membership degree:(24)xij=Xij∑i=1m∑j=1nXij.Thus, the decision matrix of membership degreexm×n will be established. The entropy value of the jth factor is(25)ej=−1lnm∑i=1myijlnyij.To avoid the failure of logarithm in equation (25) due to possible zero element in decision matrix xm×n, yij is expressed as(26)yij=1+xij∑i=1m1+xij.Finally, the weight valueωj, which reflects the influence of jth factor on strength data, can be calculated with(27)ωj=1−ejn−∑j=1nej.Data transformation is determined not only by the geometric correlation of factor space but also by the logic correlation between element of factor space and target value, which is consistent with perceptual knowledge. It provides the mathematical basis for the following multidata transformation. In this circumstance, the strength data could be transformed equivalently by the modified grey relational model, as long as the space factor is established and the element of space is quantified.


## 2.1. Fatigue Strength Model of Material

Material fatigue strength is generally evaluated with material static performance. There are many impact factors on material fatigue strength; however, there is no an authoritative theory, with a quantitative and comprehensive analysis, containing all kinds of impact factors. Therefore, instead of all kinds of factors, the primary factors that have significant impacts on material fatigue strength include tensile strength, material hardness, and inclusion size [8, 14, 31].For low-strength steel with Vickers hardness (HV) less than 400, there is the linear relationship between rotating-bending or tension-compression fatigue strengthσw,MPa and material Vickers hardness [6–8]:(1)σw≈1.6HV±0.1HVHV≤400.There is also a linear relationship between fatigue strength and tensile strengthσb,MPa for structural steel with tensile strength less than 1200 MPa:(2)σw≈0.5σbσb≤1200MPa.The essence of fatigue strength is the microplastic deformation resistance of metal material, which is different from tensile strength in physical characteristics. Thus fatigue strength is only approximately estimated with tensile strength.The mechanical properties of material FV520B used in the centrifugal compressor impeller are similar to high-strength steel (tensile strengthσb>1200MPa). The above linear relationship between fatigue strength and tensile strength does not exist when tensile strength reaches up to 1250 MPa. Inclusion has been proved to be the key factor affecting fatigue strength of high-strength steel, and Murakami and Endo [32] proposed the expression of the impacts of steel matrix hardness HV and inclusion size Ain on material fatigue strength:(3)σw=1.56HV+120Ain1/6.The fatigue strength can also be determined by substituting inclusion sizeAin with the size of GBF area AGBF [33], and the size of GBF is formulated as follows:(4)AGBF=2HV+120σw6.The influence of stress ratior and hydrogen on fatigue strength of high-strength steel is also significant. Considering the influence of stress ratio and hydrogen, the relationship between material hardness, inclusion size, and fatigue strength can be expressed as [34](5)σw=2.7HV+12015/16Ain3/161−r2α,where α is the material parameter related to the material Vickers hardness, α=0.226+HV×10−4.


## 2.2. Material Similarity

Due to the great difference of fatigue property between different materials, the similar materials with similar fatigue strength are selected in data transformation, which ensures the effectiveness and accuracy of the transformed data. Generalized tangent modulus theory, estimating generality and similarity of different materials by the tangent modulus factorΦt, plays a great role in fatigue life prediction [35–37]. It can be used to determine fatigue property of unknown material according to the known material. Φt, a product factor of tangent modulus Et, is used to express the equilibrium state of elastomer. The tangent modulus factor curve σ¯−Φt can be expressed as(6)Φt=1σ¯dσ¯dε¯=E¯tσ¯=1σ¯EtE.Combined withS-N curve, the n¯−Φt curve can be obtained by the relationship of S=σ¯ and mean fatigue life n¯=N/106 as well as the relationship of S=σ¯ and Φt:(7)n¯=NN0=1−Φ1−Φ21−Φt−Φ21/s−11/t,where σ¯andn¯ are, respectively, the dimensionless values of S and N obtained from the S-N curve and σ¯−Φtandn¯−Φt curves are the dimensionless curve taking tangent modulus factor Φt as the parametric variable. Furthermore, the σ¯−n¯−Φt curve can be obtained by putting the two curves into a diagram. Figure 1 shows the σ¯−n¯−Φt curve for four kinds of material FV520B-I, X12Cr13-I, KMN-I, and FV520B-S. The σ¯Φt curve is derived from the stress-strain curve of material, while the n¯Φt curve is derived from the S-N curve of material. Thus, the σ¯−n¯−Φt curve reflects the working performance of material in the condition of static load and symmetrical cyclic load, respectively.Figure 1
σ
¯
−
n
¯
−
Φ
t curve of these four materials.As shown in Figure1, the above four kinds of materials have a high similarity for both σ¯−Φt curve and n¯−Φt curve, and especially when Φt≥1.4, the curves for these material gradually tend to be unanimous. So the four kinds of material can be regarded as similar material, and their fatigue properties are similar. Therefore, the obtained fatigue strength data are more scientific and reliable using these materials.


## 2.3. Fatigue Strength of Components

Material fatigue strength is often obtained by experiments, and the component fatigue strength is obtained by transforming material fatigue strength, taking into account the impacts of structure shape, stress gradient, surface condition, and temperature The relationship between material fatigue strength and component fatigue strength can be expressed as [38](8)σaf,C=σA,tsc⋅ftot,af,(9)ftot,af=1fST,af⋅fGR,af2−1.0+f1,af2fm,af⋅fTI,af⋅fTP,af⋅fGS,af,where σaf,C is the endurance stress limit of the local component, σA,tsc is the material alternating stress limit, ftot,af is the overall impact factor, fST,af presents the statistical impact factor, fGR,af is the stress gradient impact factor, f1,af is the combined surface roughness and degree of forging impact factor, fm,af is the mean stress impact factor, fTI,af is the thermal impact factor, fTP,af is the technological parameter impact factor, and fGS,af is the general surface factor.The statistical impact ist assumed as a log-normal distribution of the strength parameters:(10)fST,af=σaf,90σaf,i,where σaf,90 is the material alternating stress limit with a survival probability of 90% and σaf,i is the endurance stress limit value with the required survival probability i.The stress gradient impact factor on endurance stress limit of the local component can be expressed as below:(11)fGR,af=1+σA,b/σA,tsc−12/bv⋅χ′v,where χ′ is the relative stress gradient, σA,b is the material alternating stress limit for bending, σA,tsc is the material alternating stress limit for tension or compression, ν is the material parameter, and b is the diameter of the smooth sample under flexural load.If the material sample has a mean roughness depthRZ,s=1, the respective surface condition of a component acts fully on the endurance stress limit:(12)fSR,af=1−0.22lgRt,C0.64⋅lgσUTS+0.45⋅lgRt,C0.531−0.22⋅lgRt,M0.64⋅lgσUTS+0.45⋅lgRt,M0.53,where Rt,M is the maximum roughness depth of the material sample and Rt,C is the maximum roughness depth of the component.The mean stress impact on the endurance stress limit can be expressed as(13)fm,af=σA,tscσAσm,where σA,tsc is the material alternating stress limit and σAσm is the material alternating stress limit at the mean stress σm.If the temperature is greater than100°C, the appropriate impact factor is determined as(14)fTE,af=1−aT,af⋅10−3⋅T°C−100,where aT,af is the impact coefficient of temperature and T is the temperature in centigrade.Taking into account the differing strengths of materials as a function of the effective diameter of the semi or the unfinished casting, the type of material and the technological treatment, e.g., tempering, the technological parameter impact factor can be calculated as(15)fTP,af=1Kddeff,S⋅minKddeff,Kddeff,P,where deff is the effective diameter of the semi or the unfinished casting, deff,P represents deff up to which no technological parameter influence is considered, and deff,S represents deff according to the respective material standard.General surface factorfGS,af can be defined according to special surface properties, such as shot peening, rolling, carburizing, nitriding, induction hardening, and flame hardening, and these impact factors are primarily based on professional experience.


## 2.4. Grey Relational Model of Multisource Strength Data considering Incomplete Factor

The influence of design parameters, working conditions, and transmission modes on fatigue strength is smaller than material parameters and stress ratio, especially in the normal working mode. In this paper, design parameters, working conditions, and transmission modes are in the range of normal working parameters. Combined with the fatigue strength models above, the influence of these parameters on fatigue strength follows a roughly linear law.Grey relational analysis for data with linear laws can be used to obtain a relatively precise result by considering the correlation of source data in factor space and establishing grey relational functions of target data and source data [28, 39]. However, it is difficult to obtain the strength data of the centrifugal compressor impellor, thus the source data are often incomplete. In addition, source data, formed and obtained in different periods, have obvious multisource characteristic. Thus, it is necessary to correct the transformation result acquired from limited data.


## 2.5. Grey Relational Grade for Fatigue Strength Data

Based on the point relational coefficient model, the relational grade of fatigue strength can be determined, with taking the influence of HV, tensile strength, inclusion size, etc., on material fatigue strength, as well as the influence of structure shape, stress gradient, surface condition, mean stress, temperature, etc., on component fatigue strength into account.Supposing the measured arrayX0=x01,x02,…,x0n is the characteristic array, and the reference array is denoted as(16)Xi=xi1,xi2,...,xin.The relational coefficientαik between the characteristic array and the reference array is defined as [39](17)αik=mini=1mmaxk=1nx0k−xik+ξmaxi=1mmaxk=1nx0k−xikx0k−xik+ξmaxi=1mmaxk=1nx0k−xik,where ξ is the resolution coefficient and the value is in the range of 0～1.DenoteγX0,Xi with the following formula:(18)γX0,Xi=1n∑k=1nαik.IfγX0,Xi meets the axiom of normalization, integrality, duality, and proximity [40], γX0,Xi is defined as the grey relational grade between X0andXi.The grey relational analysis is an approximate analysis method. Since the influence of each factor on fatigue strength follows an approximately linear law and source data are always incomplete, the calculation result has some deviation compared with the real data. It is necessary to reduce the deviation by residual error modification [29].For any reference arrayXi=xi1,xi2,…,xin, if X0andXi have an increment Δxi, the variation of target value fi is denoted with Δfk and the variation of first data is Δf0. The modified residual error of target value is(19)βX0,Xi=εiΔf0,where(20)εi=1n∑k=1nΔfkΔf0.Furthermore, the modified transformation model based on grey relational grade can be established as follows:(21)fi′=γX0,Xifi+βX0,Xi.In equation (21), since the original target value is modified with residual error considering individual difference of data in factor space, data deviation caused by data transformation can be decreased compared with the original target value.


## 2.6. Weight Calculation of Multisource Strength Data

Generally speaking, the factors in factor space and target data have the same parameter name but have the difference in material, construction, or working condition data. Due to the multisource data in factor space, it is necessary to weigh the influence of different types of data on transformation result, to distinguish the influence grade of each factor. Entropy-weight theory has an advantage on simple calculation and little subjective interference to calculate weight coefficient [41, 42].Based on the comprehensive evaluation method [43, 44], the entropy-weight coefficient decision matrix of strength data for centrifugal compressor is established to express the correlation among three kinds of factors that include working condition, design parameter, and construction parameter. Thus, the decision matrix is(22)Qm×n=q11⋯q1n⋮⋮⋮qm1⋯qmn.Inm×n matrix, the number of influence factors are n and each factor is assessed m times independently. The treated decision matrix Qm×n reduces subjective interference to some extent.DefineXij as the membership degree of the jth factor, and QM=maxqij,Qm=minqij,i=1,2,…,m. Then the membership degree of each factor in Qm×n is(23)Xij=Qij−QmQM−Qm,positive factors,QM−QijQM−Qm,negative factors.Normalizing each factor of the membership degree:(24)xij=Xij∑i=1m∑j=1nXij.Thus, the decision matrix of membership degreexm×n will be established. The entropy value of the jth factor is(25)ej=−1lnm∑i=1myijlnyij.To avoid the failure of logarithm in equation (25) due to possible zero element in decision matrix xm×n, yij is expressed as(26)yij=1+xij∑i=1m1+xij.Finally, the weight valueωj, which reflects the influence of jth factor on strength data, can be calculated with(27)ωj=1−ejn−∑j=1nej.Data transformation is determined not only by the geometric correlation of factor space but also by the logic correlation between element of factor space and target value, which is consistent with perceptual knowledge. It provides the mathematical basis for the following multidata transformation. In this circumstance, the strength data could be transformed equivalently by the modified grey relational model, as long as the space factor is established and the element of space is quantified.


## 3. Transformation of Multisource Strength Data of Centrifugal Compressor Impeller

The strength of the centrifugal compressor impeller can be determined by transforming the multisource data, including design parameter, working condition, material, and structure. The difference of the known strength data and unknown data of the centrifugal compressor impeller can be divided into the following types: (1) the same material and structure, but different working conditions, (2) the same structure and working condition, but different materials, (3) the similar material, structure, and working condition, but different design parameters, (4) the similar material and blade structure, but different structure of the compressor impeller, and (5) the comprehensive differences involving many impact factors.


### 3.1. Parametric Treatment of Each Impact Factor

The influence of the stress ratio on fatigue strength is discussed with material parameters, and these parameters have the largest influence on fatigue strength of the compressor impeller. The difference between these parameters for different materials is significant, which cannot be transformed according to approximate grey relational grade. Therefore, the fatigue strength of target material is obtained by analogically scaling tensile strength and Vickers hardness of material in factor space based on equation (5). The scaling factor κMi, used as the whole scaling factor of multisource data transformation, is formulated as follows:(28)κMi=σbiσb0HV0+120HVi+12016/15⋅1−r0α01−riαi⋅2αi−α0,where α=0.226+HV×10−4.The design parameters and quantifiable working condition parameters related to fatigue strength of the compressor impeller, such as blade inlet width, blade inlet diameter, blade thickness, and other typical impeller blade structure parameters, can be determined according to equations (8)–(15). Assuming that other parameters are constant, the influence of these parameters on the fatigue strength conforms to a linear law approximately. Bringing source data into factor space, the corresponding grey relational grade can be obtained by equations (17)–(21).There are some unquantifiable working condition factors, including lubrication condition, wear degree, and cleanness of fluid, also affecting fatigue strength of the compressor impeller. And it is difficult to directly qualify these factors. Thus, the evaluation mechanism of “good” or “bad” is used to quantify these factors. The rating scales are shown in Table1.Table 1
Rating scales for unquantifiable working condition factors.
Parameter name
Symbol
Good
Bad
Lubrication condition
L
1
5
Wear degree
W
1
5
Cleanness of fluid
C
1
5The impeller transmission mode affects the fatigue strength, and the actual carrying capacity is not the same even if the same impeller is in different levels. Therefore, the impeller transmisson mode, including the impeller position leveln, blade form N, impeller form M, rotor form P, and other qualitative parameters, should be quantified by the evaluation mechanism as shown in Table 1 and bring them into factor space in data transformation.


### 3.2. Comprehensive Transformation Model of Multisource Strength Data

In the above factors, the material factor acted as the independent correction factorκMi, while the other three kinds of factors have different influence on data translation in different conditions. Therefore, the weights of the three kinds of impacts factors on strength data translation should be taken into account for coordinating these factors, and the general translation model is(29)Si′=κMiωWiSWi+ωPiSPi+ωSiSSi,where Si′ is the predicted fatigue strength data of the unknown centrifugal compressor impeller, κMi is the material correction factor, SWi is the transformation data of working condition, SPi is the transformation data of the design parameter, SSi is the transformation data of the transmission mode, ωWi is the influence weight of the working condition, ωPi is the influence weight of design parameters, and ωSi is the influence weight of the transmission mode.


## 3.1. Parametric Treatment of Each Impact Factor

The influence of the stress ratio on fatigue strength is discussed with material parameters, and these parameters have the largest influence on fatigue strength of the compressor impeller. The difference between these parameters for different materials is significant, which cannot be transformed according to approximate grey relational grade. Therefore, the fatigue strength of target material is obtained by analogically scaling tensile strength and Vickers hardness of material in factor space based on equation (5). The scaling factor κMi, used as the whole scaling factor of multisource data transformation, is formulated as follows:(28)κMi=σbiσb0HV0+120HVi+12016/15⋅1−r0α01−riαi⋅2αi−α0,where α=0.226+HV×10−4.The design parameters and quantifiable working condition parameters related to fatigue strength of the compressor impeller, such as blade inlet width, blade inlet diameter, blade thickness, and other typical impeller blade structure parameters, can be determined according to equations (8)–(15). Assuming that other parameters are constant, the influence of these parameters on the fatigue strength conforms to a linear law approximately. Bringing source data into factor space, the corresponding grey relational grade can be obtained by equations (17)–(21).There are some unquantifiable working condition factors, including lubrication condition, wear degree, and cleanness of fluid, also affecting fatigue strength of the compressor impeller. And it is difficult to directly qualify these factors. Thus, the evaluation mechanism of “good” or “bad” is used to quantify these factors. The rating scales are shown in Table1.Table 1
Rating scales for unquantifiable working condition factors.
Parameter name
Symbol
Good
Bad
Lubrication condition
L
1
5
Wear degree
W
1
5
Cleanness of fluid
C
1
5The impeller transmission mode affects the fatigue strength, and the actual carrying capacity is not the same even if the same impeller is in different levels. Therefore, the impeller transmisson mode, including the impeller position leveln, blade form N, impeller form M, rotor form P, and other qualitative parameters, should be quantified by the evaluation mechanism as shown in Table 1 and bring them into factor space in data transformation.


## 3.2. Comprehensive Transformation Model of Multisource Strength Data

In the above factors, the material factor acted as the independent correction factorκMi, while the other three kinds of factors have different influence on data translation in different conditions. Therefore, the weights of the three kinds of impacts factors on strength data translation should be taken into account for coordinating these factors, and the general translation model is(29)Si′=κMiωWiSWi+ωPiSPi+ωSiSSi,where Si′ is the predicted fatigue strength data of the unknown centrifugal compressor impeller, κMi is the material correction factor, SWi is the transformation data of working condition, SPi is the transformation data of the design parameter, SSi is the transformation data of the transmission mode, ωWi is the influence weight of the working condition, ωPi is the influence weight of design parameters, and ωSi is the influence weight of the transmission mode.


## 4. Case Study

### 4.1. Fatigue Strength Prediction for Impeller

It is unlike to obtain failure strength data of impeller directly due to the shortage of failure data. Therefore, it should select the exciting source data with similar centrifugal compressor impeller to transform equivalently. The No. 5 and the predicted No. 0 impellers are used in the centrifugal coal gas compressor. The No. 1 and No. 4 impellers are used in the centrifugal ethylene compressor, and the No. 2 and No. 3 impellers are used in the centrifugal air compressor. The appearance of No. 3 and No. 5 impellers is shown in Figure2. The five reference impellers in this paper have similar material of FV520B-I, X12Cr13-I, KMN-I, and FV520B-S, respectively, as well similar structure and working condition, with predicted impeller. Meanwhile, these five reference impellers have similar stress distribution with the predicted No. 0 impeller, just like Figure 3. The reference impellers data are listed in Tables 2 and 3.Figure 2
No. 3 (a) and No. 5 (b) reference impellers.
(a)
(b)Figure 3
(a) The stress distribution of the No. 0 impeller and (b) its blade.
(a)
(b)Table 2
Design parameters of the reference impeller.
Parameter name
No. 0
No. 1
No. 2
No. 3
No. 4
No. 5
Width of blade inlet
b
1 (mm)
38.9
29
72
57
45
59
Diameter of blade inlet
D
1 (mm)
357
225
480
340
245
324
Blade thickness
t (mm)
2.5
2.8
3.8
6
5.6
5
Velocity of blade inlet
c
1 (m/s)
73
134.5
155.5
89.5
81.2
120
Installation angle of blade inlet
β
1A (°)
33
30
25
25
20
28
Slope of wheel cover
θ (°)
9.12
9
7.5
0
6
9.1
Curvature radius of blade
R (mm)
434.5
325
576
413
397
512Table 3
Parameter of working condition, transmission mode, and material property of the reference impeller.
Parameter name
No. 0
No. 1
No. 2
No. 3
No. 4
No. 5
Input speed
n
0 (r/min)
8877
22500
6667
5035
11200
8200
Inlet pressure
P
in (kg/cm2)
1.02
1.025
1.1005
0.955
1.008
1.048
Inlet temperature
t
in (°C)
27
21
30
37
29
32
Outlet pressure
P
out (kg/cm2)
1.78
2.954
3.533
3.793
2.223
1.898
Lubrication condition
L
1
1
1
1
1
1
Wear degree
W
0
1
0
0
0
1
Cleanness of fluid
C
1
1
2
3
3
2
Position of impeller
N
1
1
1
1
1
1
Impeller form
N
1
1
2
1
1
1
Blade form
F
2
2
2
3
2
3
Rotor form
W
R
1
1
1
3
1
2
Stress ratio
R
0.82
0.72
0.85
0.77
0.81
0.79
Material
M
FV520B-I
FV520B-I
X12Cr13-I
KMN-I
FV520B-S
FV520B-S
Tensile strength
σ
b
1078
1078
735
931
931
887
Vickers hardness
HV
370
370
280
320
320
300
Fatigue strength value
S (MPa)
—
662
445
529
476
579The four kinds of material can be regarded as similar material, and their fatigue properties are similar. Meanwhile, these six impellers have roughly similar structure, shape, and operating mode, so they have similar stress distribution and fatigue damage positions. Therefore, it is a reasonable assumption that all the impellers have a similar fatigue damage mechanism, and similar nondimensional fatigue damage versus the nondimensional number of cycles.In order to show multisource data transformation with the method presented above, the No. 1 reference impeller is exampled as follows.


### 4.2. Transformation of Multisource Design Parameter Data

The factor space of design parameters for the No. 0 compressor impeller isX0=38.93572.573339.12434.5, and the factor space for the No. 1 impeller is X1=292252.8134.5309325. If ξ=0.5, according to equations (17) and (18), the relational grade of two factor spaces is(30)γPX0,X1=1n∑k=1nαik=0.7155.If the same proportion of incrementΔxi=0.1xi is added to each element of X0andX1, the residual error for fatigue strength data affected by design parameters is βPX0,X1=97.3MPa according to equations (19)–(21).Thus, the transformed fatigue strength for the No. 1 impeller in the condition of multisource design parameters is(31)SP1′=γPX0,X1S1+βPX0,X1=570.9MPa.


### 4.3. Transformation of Multisource Structure Parameter Data

According to Table3, the factor space of the structure for the No. 0 compressor impeller and No. 1 compressor impeller is(32)XS0=ni,Ni,Fi,WRi=1121,XS1=n1,N1,F1,WR1=1121.The relational grade of two factor spaces isγSX0,X1=1, and the transformed fatigue strength is 662 MPa.


### 4.4. Transformation of Multisource Working Condition Data

The working conditions, such as lubrication, wear, and cleanliness, are classified and are evaluated by the professional engineers in factory according to Table1. The factor space is established as follows.(33)
X
W
0
=
P
in
0
,
P
out
0
,
n
00
,
t
in
0
,
L
0
,
W
0
,
C
0
=
1.02
,
1.78
,
8877
,
27
,
1
,
0
,
1
,
X
W
1
=
P
in
1
,
P
out
1
,
n
01
,
t
in
1
,
L
1
,
W
1
,
C
1
=
1.025
,
2.954
,
22500
,
21
,
1
,
1
,
1
.The relational grade of two factor spaces isγWX0,X1=0.7325, and the residual error of fatigue strength data affected by working conditions is βWX0,X1=75MPa. So the transformed fatigue strength value is(34)SW1′=γWX0,X1S1+βWX0,X1=559.1MPa.


### 4.5. Transformation of Multisource Material Parameter Data

According to equation (28), the scaling coefficient of data transformation κMi=0.89.


### 4.6. Weight Calculation of Multisource Data

Table4 confirms the weight of each element in factor space. It is available to confirm the influence of each subfactor on impeller fatigue strength by grading.Table 4
Factor set, subfactor set, and the rating scale.
Factor sets
Subfactor sets
Score
Parameter (U1)
u
11
Width of blade inlet
4
u
12
Diameter of blade inlet
3
u
13
Installation angle of blade inlet
1
u
14
Curvature radius of blade
1
u
15
Blade thickness
5
u
16
Velocity of blade inlet
3
u
17
Slope of wheel cover
1
Working condition (U2)
U
21
Inlet pressure
2
U
22
Outlet pressure
3
U
23
Input speed
4
U
24
Inlet temperature
2
u
31
Lubrication condition
2
Working condition (U3)
u
32
Wear degree
2
u
33
Cleanness of fluid
1According to Table4 and judging one by one, based on the comprehensive evaluation method, the decision matrix of entropy weight is(35)Qm×n=3.32.64.53.22.44.14.22.54.83.62.14.3.According to the above decision matrix, the membership degree matrixxm×n can be obtained by equations (23) and (27) and the weight value of the corresponding three impact factors is ωPi=0.345, ωSi=0.267, and ωWi=0.363 respectively.Taking the weight value into equation (29), the fatigue strength value for the No. 1 impeller after transformation is(36)S1′=0.89×0.345×570.9+0.267×662+0.363×559.1=513MPa.Similarly, the equivalent data of the fatigue strength for the else four impellers can also be obtained as shown in Table5.Table 5
Equivalent data for fatigue strength after transformation.
Impeller number
No. 1
No. 2
No. 3
No. 4
No. 5
Equivalent data (MPa)
513
472
622
507
598


### 4.7. Effectiveness Discussion for the Transformed Fatigue Strength Data

The transformed fatigue strength data with the proposed method in this paper and the theoretical strength data of centrifugal compressor impeller are presented as below:(37)S1′=513,S2′=472,S3′=622,S4′=507,S5′=598,S1=624,S2=420,S3=592,S4=467,S5=505.Taking the distribution of theoretical fatigue strength data as reference, the effectiveness of the transformation data can be discussed by judging whether the transformation data and theoretical data obey the approximate population distribution. According to hypothesis testing rules, compatibility and consistency for the transformation data and theoretical data are solved, that the mean value of strength and the variance of strength are tested byt-test and F-test, respectively. The critical value is t0.058 in the condition of significant level α=0.1. The test statistics is calculated as follows:(38)t=542.4−521.675.43×1/5+1/5=0.4360≤t0.058=1.8595.Similarly,(39)F0.954,4=0.1565≤F=4125.37255.29=0.5686≤F0.054,4=6.39.There is no significant difference between the transformation data and theoretical data in the mean value and the variance. Therefore, these two sets of data are consistent and approximately obey the same population distribution. It shows that the transformation data have certain credibility and the transformed strength can be regarded as the fatigue strength of the No. 0 impeller.In Figure3, the maximum stress of the impeller, as well as the most likely to fatigue damage position, is located near the blade root connecting with wheel cover. Since the tensile strength of material FV520B-I for the No. 0 centrifugal compressor impeller is close to the high-strength steel up to 1100 MPa, ultrasonic fatigue testing is applied to verify the failure stress state of this position and test the rationality of transformation strength data. The S-N curve of FV520B-I is shown in Figure 4.Figure 4
Experimental result of fatigue property for material FV520B-I.In Figure4, there is no fatigue limit at 107 cycles and the fatigue strength is about 550 MPa corresponding to the stress amplitude at 4×108∼109 cycles. By comparing the transformed fatigue strength data obtained by the proposed method in this paper and the material strength data of the centrifugal compressor impeller, we see that the maximum difference in the 5 studied cases is 15%, which meets the requirements of safety margin (100% or 150%) according to codes. Therefore, it shows the transformed strength data are effective.


## 4.1. Fatigue Strength Prediction for Impeller

It is unlike to obtain failure strength data of impeller directly due to the shortage of failure data. Therefore, it should select the exciting source data with similar centrifugal compressor impeller to transform equivalently. The No. 5 and the predicted No. 0 impellers are used in the centrifugal coal gas compressor. The No. 1 and No. 4 impellers are used in the centrifugal ethylene compressor, and the No. 2 and No. 3 impellers are used in the centrifugal air compressor. The appearance of No. 3 and No. 5 impellers is shown in Figure2. The five reference impellers in this paper have similar material of FV520B-I, X12Cr13-I, KMN-I, and FV520B-S, respectively, as well similar structure and working condition, with predicted impeller. Meanwhile, these five reference impellers have similar stress distribution with the predicted No. 0 impeller, just like Figure 3. The reference impellers data are listed in Tables 2 and 3.Figure 2
No. 3 (a) and No. 5 (b) reference impellers.
(a)
(b)Figure 3
(a) The stress distribution of the No. 0 impeller and (b) its blade.
(a)
(b)Table 2
Design parameters of the reference impeller.
Parameter name
No. 0
No. 1
No. 2
No. 3
No. 4
No. 5
Width of blade inlet
b
1 (mm)
38.9
29
72
57
45
59
Diameter of blade inlet
D
1 (mm)
357
225
480
340
245
324
Blade thickness
t (mm)
2.5
2.8
3.8
6
5.6
5
Velocity of blade inlet
c
1 (m/s)
73
134.5
155.5
89.5
81.2
120
Installation angle of blade inlet
β
1A (°)
33
30
25
25
20
28
Slope of wheel cover
θ (°)
9.12
9
7.5
0
6
9.1
Curvature radius of blade
R (mm)
434.5
325
576
413
397
512Table 3
Parameter of working condition, transmission mode, and material property of the reference impeller.
Parameter name
No. 0
No. 1
No. 2
No. 3
No. 4
No. 5
Input speed
n
0 (r/min)
8877
22500
6667
5035
11200
8200
Inlet pressure
P
in (kg/cm2)
1.02
1.025
1.1005
0.955
1.008
1.048
Inlet temperature
t
in (°C)
27
21
30
37
29
32
Outlet pressure
P
out (kg/cm2)
1.78
2.954
3.533
3.793
2.223
1.898
Lubrication condition
L
1
1
1
1
1
1
Wear degree
W
0
1
0
0
0
1
Cleanness of fluid
C
1
1
2
3
3
2
Position of impeller
N
1
1
1
1
1
1
Impeller form
N
1
1
2
1
1
1
Blade form
F
2
2
2
3
2
3
Rotor form
W
R
1
1
1
3
1
2
Stress ratio
R
0.82
0.72
0.85
0.77
0.81
0.79
Material
M
FV520B-I
FV520B-I
X12Cr13-I
KMN-I
FV520B-S
FV520B-S
Tensile strength
σ
b
1078
1078
735
931
931
887
Vickers hardness
HV
370
370
280
320
320
300
Fatigue strength value
S (MPa)
—
662
445
529
476
579The four kinds of material can be regarded as similar material, and their fatigue properties are similar. Meanwhile, these six impellers have roughly similar structure, shape, and operating mode, so they have similar stress distribution and fatigue damage positions. Therefore, it is a reasonable assumption that all the impellers have a similar fatigue damage mechanism, and similar nondimensional fatigue damage versus the nondimensional number of cycles.In order to show multisource data transformation with the method presented above, the No. 1 reference impeller is exampled as follows.


## 4.2. Transformation of Multisource Design Parameter Data

The factor space of design parameters for the No. 0 compressor impeller isX0=38.93572.573339.12434.5, and the factor space for the No. 1 impeller is X1=292252.8134.5309325. If ξ=0.5, according to equations (17) and (18), the relational grade of two factor spaces is(30)γPX0,X1=1n∑k=1nαik=0.7155.If the same proportion of incrementΔxi=0.1xi is added to each element of X0andX1, the residual error for fatigue strength data affected by design parameters is βPX0,X1=97.3MPa according to equations (19)–(21).Thus, the transformed fatigue strength for the No. 1 impeller in the condition of multisource design parameters is(31)SP1′=γPX0,X1S1+βPX0,X1=570.9MPa.


## 4.3. Transformation of Multisource Structure Parameter Data

According to Table3, the factor space of the structure for the No. 0 compressor impeller and No. 1 compressor impeller is(32)XS0=ni,Ni,Fi,WRi=1121,XS1=n1,N1,F1,WR1=1121.The relational grade of two factor spaces isγSX0,X1=1, and the transformed fatigue strength is 662 MPa.


## 4.4. Transformation of Multisource Working Condition Data

The working conditions, such as lubrication, wear, and cleanliness, are classified and are evaluated by the professional engineers in factory according to Table1. The factor space is established as follows.(33)
X
W
0
=
P
in
0
,
P
out
0
,
n
00
,
t
in
0
,
L
0
,
W
0
,
C
0
=
1.02
,
1.78
,
8877
,
27
,
1
,
0
,
1
,
X
W
1
=
P
in
1
,
P
out
1
,
n
01
,
t
in
1
,
L
1
,
W
1
,
C
1
=
1.025
,
2.954
,
22500
,
21
,
1
,
1
,
1
.The relational grade of two factor spaces isγWX0,X1=0.7325, and the residual error of fatigue strength data affected by working conditions is βWX0,X1=75MPa. So the transformed fatigue strength value is(34)SW1′=γWX0,X1S1+βWX0,X1=559.1MPa.


## 4.5. Transformation of Multisource Material Parameter Data

According to equation (28), the scaling coefficient of data transformation κMi=0.89.


## 4.6. Weight Calculation of Multisource Data

Table4 confirms the weight of each element in factor space. It is available to confirm the influence of each subfactor on impeller fatigue strength by grading.Table 4
Factor set, subfactor set, and the rating scale.
Factor sets
Subfactor sets
Score
Parameter (U1)
u
11
Width of blade inlet
4
u
12
Diameter of blade inlet
3
u
13
Installation angle of blade inlet
1
u
14
Curvature radius of blade
1
u
15
Blade thickness
5
u
16
Velocity of blade inlet
3
u
17
Slope of wheel cover
1
Working condition (U2)
U
21
Inlet pressure
2
U
22
Outlet pressure
3
U
23
Input speed
4
U
24
Inlet temperature
2
u
31
Lubrication condition
2
Working condition (U3)
u
32
Wear degree
2
u
33
Cleanness of fluid
1According to Table4 and judging one by one, based on the comprehensive evaluation method, the decision matrix of entropy weight is(35)Qm×n=3.32.64.53.22.44.14.22.54.83.62.14.3.According to the above decision matrix, the membership degree matrixxm×n can be obtained by equations (23) and (27) and the weight value of the corresponding three impact factors is ωPi=0.345, ωSi=0.267, and ωWi=0.363 respectively.Taking the weight value into equation (29), the fatigue strength value for the No. 1 impeller after transformation is(36)S1′=0.89×0.345×570.9+0.267×662+0.363×559.1=513MPa.Similarly, the equivalent data of the fatigue strength for the else four impellers can also be obtained as shown in Table5.Table 5
Equivalent data for fatigue strength after transformation.
Impeller number
No. 1
No. 2
No. 3
No. 4
No. 5
Equivalent data (MPa)
513
472
622
507
598


## 4.7. Effectiveness Discussion for the Transformed Fatigue Strength Data

The transformed fatigue strength data with the proposed method in this paper and the theoretical strength data of centrifugal compressor impeller are presented as below:(37)S1′=513,S2′=472,S3′=622,S4′=507,S5′=598,S1=624,S2=420,S3=592,S4=467,S5=505.Taking the distribution of theoretical fatigue strength data as reference, the effectiveness of the transformation data can be discussed by judging whether the transformation data and theoretical data obey the approximate population distribution. According to hypothesis testing rules, compatibility and consistency for the transformation data and theoretical data are solved, that the mean value of strength and the variance of strength are tested byt-test and F-test, respectively. The critical value is t0.058 in the condition of significant level α=0.1. The test statistics is calculated as follows:(38)t=542.4−521.675.43×1/5+1/5=0.4360≤t0.058=1.8595.Similarly,(39)F0.954,4=0.1565≤F=4125.37255.29=0.5686≤F0.054,4=6.39.There is no significant difference between the transformation data and theoretical data in the mean value and the variance. Therefore, these two sets of data are consistent and approximately obey the same population distribution. It shows that the transformation data have certain credibility and the transformed strength can be regarded as the fatigue strength of the No. 0 impeller.In Figure3, the maximum stress of the impeller, as well as the most likely to fatigue damage position, is located near the blade root connecting with wheel cover. Since the tensile strength of material FV520B-I for the No. 0 centrifugal compressor impeller is close to the high-strength steel up to 1100 MPa, ultrasonic fatigue testing is applied to verify the failure stress state of this position and test the rationality of transformation strength data. The S-N curve of FV520B-I is shown in Figure 4.Figure 4
Experimental result of fatigue property for material FV520B-I.In Figure4, there is no fatigue limit at 107 cycles and the fatigue strength is about 550 MPa corresponding to the stress amplitude at 4×108∼109 cycles. By comparing the transformed fatigue strength data obtained by the proposed method in this paper and the material strength data of the centrifugal compressor impeller, we see that the maximum difference in the 5 studied cases is 15%, which meets the requirements of safety margin (100% or 150%) according to codes. Therefore, it shows the transformed strength data are effective.


## 5. Conclusions

To obtain the fatigue strength of the centrifugal compressor impeller quickly, a fatigue strength data transformation model based on modified grey relational theory is established, combining with the impact of tensile strength, structure size, mean stress, temperature, etc., on fatigue strength. This model has the following characteristics:(1)
The grey relational model is established to analyze the impact of each factor on fatigue strength, and it is modified with residual error to reduce the deviation by source strength data incompletion.(2)
The multisource strength data transformation model and the decision matrix of factor space based on entropy-weight theory are presented, to evaluate the impact of multisource data, such as structure parameter, working condition, and transmission mode, on fatigue strength comprehensively.(3)
Taking the influence of nonquantitative parameters into account, the working condition and transmission mode are quantified with the related standard. By analyzing the influences of quantitative and nonquantitative parameters comprehensively, the transformed strength data are more scientific.The very-high-cycle fatigue test of FV520B-I is applied to verify the rationality of the fatigue strength transformation model, which agrees with multisource data transformation result well; therefore, the above model and method are effective.



---

*Source: 1236130-2020-01-30.xml*