OOA-Optimized Kriging-RBF Method for Radial Deformation Reliability Analysis of Compressor Blade Angle Crack

1. Introduction

Compressor blade radial deformation poses a critical design challenge in modern aeroengines, where precise tip clearance control directly influences aerodynamic efficiency, operational safety, and engine performance [1,2,3]. Under operational conditions, blades experience complex multi-physics loading, including centrifugal forces, aerodynamic pressures, thermal gradients, and vibratory excitations, resulting in intricate radial displacement patterns that must be carefully managed to prevent blade-casing contact while minimizing tip leakage losses [4,5]. Angle crack defects, which typically occur at blade tip corners due to stress concentration and high-cycle fatigue, significantly complicate this deformation behavior by introducing local compliance reductions, altered load paths, and time-dependent structural property changes [6,7]. These crack-induced effects create highly nonlinear relationships between operational parameters and radial clearance, necessitating sophisticated reliability assessment methodologies to ensure safe and efficient engine operation under uncertainty.

Traditional structural reliability assessment approaches for compressor blade systems predominantly rely on first-order reliability methods (FORM), second-order reliability methods (SORM), and Monte Carlo simulation techniques [8,9,10]. While FORM and SORM offer computational efficiency through linearization or quadratic approximation of limit state functions, they often fail to capture the complex nonlinear coupling effects inherent in crack-damaged blade structures, particularly when multiple uncertainty sources interact simultaneously [11,12]. Monte Carlo simulation provides accurate reliability estimates but requires extensive computational resources due to the need for repeated high-fidelity finite element analyses, making it impractical for design optimization and real-time monitoring applications. Furthermore, conventional reliability methods struggle with the multi-modal response surfaces and discontinuous behavior introduced by crack propagation, leading to inaccurate failure probability estimates and potentially unsafe design margins [13]. These limitations highlight the critical need for advanced computational approaches that can balance accuracy and efficiency in modeling crack-structure interactions, motivating the development of surrogate modeling techniques for reliability assessment of damaged blade systems.

Surrogate modeling techniques have emerged as promising alternatives to address computational challenges in structural reliability assessment, with Kriging, Gaussian Process (GP) models and Radial Basis Function (RBF) networks gaining particular attention for their interpolation capabilities and uncertainty quantification properties [14,15,16]. Su et al. proposed a GP-based dynamic surrogate model for complex engineering structural reliability analysis [17]. However, individual surrogate models exhibit inherent limitations when applied to crack-influenced blade deformation problems. GP models excel in uncertainty quantification and global trend capture but may struggle with sharp response transitions and discontinuities introduced by crack boundaries [18]. Conversely, RBF networks demonstrate superior nonlinear mapping capabilities and local approximation performance but provide limited uncertainty information essential for reliability analysis [19,20]. Chen et al. developed a multi-fidelity data aggregation using convolutional neural networks for multi-fidelity modeling [21]. Recent research efforts have explored hybrid surrogate modeling approaches to leverage complementary strengths of different metamodeling techniques. Wang et al. introduced a new adaptive extreme response surface approach for time-variant reliability problems [22]. Xia et al. developed a hybrid approach to seismic reliability assessment [23]. Additionally, Rabhi et al. proposed a combined a dynamic reliability method and a me-ta-model (reduced model) to obtain good results in the reliability and optimization of physical systems [24]. However, most existing hybrid schemes focus on simple ensemble averaging or sequential correction methods without addressing the specific challenges posed by crack-induced response characteristics.

Optimization algorithms play a crucial role in surrogate model parameter tuning and hybrid weight determination, with evolutionary algorithms, swarm intelligence methods, and gradient-based techniques being widely employed [25,26]. Traditional genetic algorithms and particle swarm optimization approaches often suffer from premature convergence and local optima entrapment when dealing with high-dimensional parameter spaces typical of hybrid surrogate models. Liu et al. developed an efficient reliability analysis framework for OST using the kriging model of optimal linear unbiased estimation [27]. More recent nature-inspired algorithms such as whale optimization, grey wolf optimization, and artificial bee colony algorithms have shown improved global search capabilities but may lack sufficient exploitation ability for fine-tuning complex hyperparameter configurations [28,29,30]. Liu et al. developed a method combining the shear strength reduction (SSR) technique, the surrogate model and the adaptive pool-based sampling strategy for efficient slope reliability analysis [31]. Kumar et al introduced a data-driven model approximating the relationship between the inputs and outputs by using an adaptive sparse polynomial chaos expansion approach [32]. Wong et al. developed a novel method considering the equality constraints, inequality and physical constraints of the hydro-thermal systems, which combined with simulated annealing, genetic algorithms, and two hybrid optimization techniques [33]. However, these conventional optimization algorithms still face significant limitations in handling the complex multi-modal landscapes characteristic of crack-damaged systems. The Osprey Optimization Algorithm (OOA) offers distinct advantages over traditional methods through its unique dual-phase search mechanism that combines aggressive exploration via high-altitude soaring behavior with precise exploitation through targeted diving strategies, enabling superior performance in navigating complex, multi-modal parameter spaces while maintaining computational efficiency [34]. Current optimization-based surrogate modeling approaches face two critical limitations: (1) poor exploration-exploitation balance in high-dimensional spaces, resulting in suboptimal configurations; and (2) insufficient efficiency and accuracy in optimization processes.

To address these challenges, this study proposes a novel OOA-optimized Kriging-RBF (OOA-KR) method for efficient reliability assessment of compressor blade radial clearance considering angle crack defects. The proposed approach integrates the uncertainty quantification capabilities of Kriging with the nonlinear mapping strengths of RBF neural networks through an adaptive weighting scheme optimized by the Osprey Optimization Algorithm (OOA). The OOA algorithm combines global exploration abilities inspired by osprey hunting behavior with local exploitation mechanisms derived from fish positioning dynamics, providing superior parameter optimization performance for hybrid surrogate model construction. Multiple uncertainty sources including crack geometry variations, operational temperature fluctuations, loading condition changes, and material property scatter are systematically incorporated to evaluate their effects on blade radial deformation distribution and reliability assessment accuracy.

The paper is organized as follows: Section 2 introduces the OOA-optimized Kriging-RBF methodology. Section 3 presents the structural reliability assessment framework. Section 4 describes the experimental setup and results. Section 5 concludes the study.

2. OOA-Optimized Kriging-RBF Method (OOA-KR)

The proposed OOA-KR method represents a novel hybrid metamodeling approach that synergistically integrates the global modeling capability of Kriging with the local approximation advantages of RBF neural networks through an adaptive weighted fusion strategy. This hybrid framework leverages Kriging’s proficiency in uncertainty quantification and global trend modeling, while capitalizing on RBF’s superior performance in capturing nonlinear local features and handling high-dimensional scattered data. The key innovation lies in the simultaneous optimization of fusion weights and hyperparameters of both base models using the Osprey Optimization Algorithm (OOA), which combines the global exploration ability of eagle search behavior with the local exploitation capability of fish swarm dynamics. The overall framework of the proposed method is illustrated in Figure 1, demonstrating the complete workflow from data preprocessing to performance function modeling.

The mathematical formulation of the OOA-KR method can be expressed as:

$$Y(x)={\lambda }_1{g}_1(x)+{\lambda }_2{g}_2(x)$$

(1)

where g₁ (x) represents the Kriging surrogate model, g₂(x) denotes the RBF neural network model, and λ₁, λ₂ are the adaptive fusion weights subject to the constraint λ₁+λ₂=1.

The OOA-KR method follows a systematic approach consisting of the following steps:

Step 1: Collect and preprocess the training dataset D={(xᵢ, yᵢ), i=1, 2, ..., n}, including data normalization, outlier detection, and quality assessment to ensure data integrity.

Step 2: Construct the Kriging surrogate model g₁(x) based on the training data, establishing optimal spatial correlation functions and estimating model parameters through maximum likelihood estimation.

Step 3: Build the RBF neural network surrogate model g₂(x) using the same training dataset, determining the optimal network architecture, basis function types, and center locations.

Step 4: Design the weighted fusion strategy Y(x) = λ₁g₁(x) + λ₂g₂(x) with the constraint λ₁+λ₂ = 1 to combine the predictions from both base models while maintaining prediction consistency.

Step 5: Define the multi-objective optimization function incorporating cross-validation error metrics (RMSE, R²) and model complexity penalties as the fitness function for OOA.

Step 6: Initialize OOA parameters including population size, maximum iterations, search space bounds, and algorithm control parameters for both osprey hunting and fish positioning phases.

Step 7: Execute OOA iterative optimization to simultaneously determine optimal fusion weights {λ₁, λ₂} and hyperparameters of both base models through bi-phase population update mechanisms.

Step 8: Construct the final optimized fusion surrogate model using the obtained optimal parameters and evaluate its performance on independent test datasets through comprehensive validation metrics.

2.1. Kriging

Kriging, originally developed by Krige for geostatistical applications, has emerged as a powerful metamodeling technique for approximating expensive computational models in engineering design and reliability analysis. Unlike polynomial-based response surface methods that rely on global fitting, kriging provides an interpolative approach that exactly reproduces the observed data while offering optimal predictions at unsampled locations with quantified uncertainty bounds.

The kriging metamodel treats the deterministic computer response as a realization of a stochastic process, expressed mathematically as [14]:

$$g_1(x)=f{(x)}^T\beta +Z(x)$$

(2)

where g₁(x) represents the unknown function of interest, f(x)=[f₁(x), f₂(x), …, f_p(x)]^T denotes a vector of known regression functions, β=[β₁, β₂, …, β_p] is the vector of regression coefficients, and Z(x) represents a zero-mean Gaussian random process with covariance structure:

$$\mathrm{Cov}\text{}\left[Z(x_i),Z(x_j)\right]={\alpha }^{2}R(x_i,x_j)$$

(3)

where term f^T(x)β captures the global trend of the response surface, while Z(x) models the local deviations to ensure interpolation through the training data points.

Given a set of n training points {x⁽¹⁾, x⁽²⁾, …, x⁽ⁿ⁾} and their corresponding responses {y⁽¹⁾, y⁽²⁾, …, y⁽ⁿ⁾}, the unknown parameters are estimated using maximum likelihood estimation:

$$\widehat{\beta }={(F^T{R}^{-1}F)}^{-1}{F}^T{R}^{-1}y$$

(4)

$${\widehat{\sigma }}^2={\displaystyle \frac{1}{n}}{(y-F\widehat{\beta })}^T{R}^{-1}(y-F\widehat{\beta })$$

(5)

where F is n×p regression matrix with elements F_ij=f(x_i), R is the n×n correlation matrix with elements R=R(x_i, x_j), and y=[ y⁽¹⁾, y⁽²⁾, …, y⁽ⁿ⁾].

The correlation function R(x_i, x_j) governs the smoothness and local behavior of the metamodel. Commonly employed correlation functions include the exponential, Gaussian, and Matérn families. For multivariate inputs, a separable correlation structure is typically assumed:

$$R(x_i,x_j)={\displaystyle \prod _{k=1}^{d}R}\text{}\left(x_{ik},x_{jk}\right)$$

(6)

where d is the input dimension and R represents the correlation function.

The kriging predictor at an arbitrary point x is given by:

$$\begin{array}{l}g(x)={\mathit{f}}^T(x)\beta +{\mathit{r}}^T(x)\alpha \\ s.t.\left\{\begin{array}{l}r(x)={[R(x,x^{(1)}),R(x,x^{(2)}),\dots ,R(x,x^{(n)})]}^T\\ \alpha =R^{-1}(y-F\widehat{\beta })\end{array}\right.\end{array}$$

(7)

2.2. RPF

Radial basis function (RBF) networks represent a class of neural network architectures that have gained considerable attention in engineering metamodeling due to their universal approximation capabilities and relatively simple mathematical formulation. Unlike traditional feedforward neural networks, RBF networks employ radially symmetric basis functions that provide localized responses, making them particularly suitable for interpolating scattered data in high-dimensional spaces.

The RBF metamodel approximates an unknown function y(x) as a linear combination of radial basis functions [16]:

$$g_2(x)={\displaystyle \sum _{i=1}^n{w}_i}\varphi ‖(x_j-c_{ij})‖+p(x)$$

(8)

where w_i represents the weight associated with the i-th basis function, ϕ(⋅) denotes the radial basis function, c_i are the center locations(typically coinciding with training points), ‖⋅‖ represents the Euclidean norm, and p(x) is an optional polynomial term to ensure well-posedness.

Common types of RBF include the Gaussian function:

$${\varphi }_r=\left\{\begin{array}{l}\exp (-(\epsilon r{)}^2)\\ \sqrt{1+{(\epsilon r)}^2}\\ 1/\sqrt{1+{(\epsilon r)}^2}\\ r^2\ln (r),wherer=‖x-c_i‖\end{array}\right.$$

(9)

where ε is a shape parameter controlling the width of the basic functions.

2.3. OOA -Optimized Weighted Fusion and Hyper-Parameter Solution

The OOA-KR method formulates the hyperparameter optimization as a multi-dimensional problem where the decision variables Θ include fusion weights {λ₁, λ₂}, Kriging parameters {θ_k}, and RBF parameters {ε, c_i, w_i}. The objective function balances prediction accuracy and model generalization:

$$\min \Theta$$

(10)

Similar to other optimization algorithms, the population is randomly initialized in the optimization space:

$$x_{i,j}=l{b}_j+r\cdot (u{b}_j-l{b}_j)$$

(11)

where x_i,j is index individuals; lb_j is the lower bound of the search space for variable j; ub_j is the upper bound; r is a random number uniformly distributed in the interval [0,1].

The osprey is a powerful hunter. Due to its strong eyesight, it can detect the position of underwater fish. After determining the fish’s location, they attacked it and hunted for it underwater. The first stage of population update in OOA is modeled based on the simulation of the natural behavior of ospreys. Modeling the fish attacked by ospreys can lead to significant changes in the ospreys’ positions in the search space, which enhances the exploration ability of OOA in identifying the best areas and escaping local optima.

In OOA design, for each osprey, the positions of other ospreys with better objective function values in the search space are regarded as underwater fish. The fish group for each osprey is specified by

$$F{P}_i=\left\{X_k\right|k\in \{1,2,\dots ,N\}\wedge F_k<F_i\}\cup \{X_{best}\}$$

(12)

where FP_i is the fish set of the i-th osprey, and X_best is the best position for the osprey.

The osprey randomly detected the position of one of the fish and attacked it. Based on the simulation of the osprey’s movement towards the fish, the new position of the corresponding osprey is calculated using Eq. (13). If the new position is better, replace the osprey’s previous position according to Eq. (14)

$$x_{i,j}^{P1}=x_{i,j}+r_{i,j}\cdot (S{F}_{i,j}-I_{i,j}\cdot x_{i,j})s.t.x_{i,j}^{P1}=\left\{\begin{array}{l}{x}_{i,j}^{P1},l{b}_j\le x_{i,j}^{P1}\le u{b}_j\\ l{b}_j,x_{i,j}^{P1}<l{b}_j\\ u{b}_j,x_{i,j}^{P1}>u{b}_j\end{array}\right.$$

(13)

$$X_i=\left\{\begin{array}{l}{X}_i^{P1},F_i^{P1}<F_i\\ X_i,else\end{array}\right.$$

(14)

wherein SF is the fish selected by the osprey, r is a random number between 0 and 1, and the value of I is one of {1,2}.

After hunting a fish, an osprey will take it to a suitable (safe) location where it eats. The second stage of population update in OO 8 is based on modeling the simulation of the natural behavior of ospreys. Modeling that brings the fish to the right position leads to minor changes in the osprey’s position in the search space, which increases the utilization ability of OOA in local search and conversions to better solutions near the discovered ones.

In the design of OOA, to simulate this natural behavior of ospreys, first of all, for each member of the population, a new random position is calculated using Eq. (15) as the “position suitable for eating fish”. Then, if the value of the objective function is improved at this new position, it replaces the previous position of the corresponding osprey according to Eq. (16).

$$\begin{array}{l}{x}_{i,j}^{P2}=x_{i,j}+{\displaystyle \frac{l{b}_j+r\cdot (u{b}_j-l{b}_j)}{t}},i=1,2,\dots ,N;j=1,2,\dots ,m;t=1,2,\dots ,T\\ s.t.x_{i,j}^{P2}=\left\{\begin{array}{l}{x}_{i,j}^{P2}.l{b}_j\le x_{i,j}^{P2}\le u{b}_j\\ l{b}_j,x_{i,j}^{P2}<l{b}_j\\ u{b}_j.x_{i,j}^{P2}>u{b}_j\end{array}\right.\end{array}$$

(15)

$$X_i=\left\{\begin{array}{l}{X}_i^{P2},F_i^{P2}<F_i\\ X_i,else\end{array}\right.$$

(16)

Here, t is the number of iterations, and T is the maximum number of iterations.

The OOA algorithm optimizes the fusion weights and hyperparameters through the following systematic steps:

Step 1: Initialize osprey population with random parameter vectors within predefined bounds, ensuring each individual represents a complete set of fusion weights and model hyperparameters.

Step 2: Evaluate fitness of each individual by constructing Kriging and RBF models with corresponding hyperparameters, computing fusion predictions, and calculating cross-validation error metrics.

Step 3: Identify the fish population for each osprey based on fitness ranking, where better-performing individuals serve as potential targets for position updates.

Step 4: Hunting phase by updating osprey positions through simulated fish attack behavior, promoting global exploration of the parameter space.

Step 5: Feeding phase by fine-tuning osprey positions through local search mechanisms, enhancing exploitation of promising regions.

Step 6: Update population by replacing inferior solutions with improved ones, maintaining diversity while progressing toward optimal parameter configurations.

Step 7: Check convergence criteria (maximum iterations, fitness threshold, or population stagnation) and terminate if satisfied, otherwise return to Step 2.

Step 8: Extract optimal fusion weights and hyperparameters from the best individual, construct the final OOA-KR model, and validate performance on test data.

The algorithm of OOA-KR is shown as follows:

3. Structural Reliability Assessment of Compressor Blade Angle Crack Based on OOA-KR Method

Reliability theory is an important tool to deal with uncertainty in the design and evaluation of engineering structures. The core of structural reliability analysis is to quantify the failure probability, that is, the probability that the structure fails to meet the performance requirements under the given random variable conditions, as shown in Figure 2.

For compressor blade radial clearance deformation, the failure probability P_f is expressed by

$$P_f={\displaystyle {\int }_{g\left(X\right)\le 0}f\left(x\right)}dx$$

(17)

where g(X) is the limit state function; g(X)≤0 stands for structural security boundary; f(x) is the joint probability density function of random variables.

The limit state function for compressor blade radial clearance deformation is defined as [1,3]:

$$g(x)={\delta }_{allow}-{\delta }_{radial}(x)$$

(18)

where δ_allow represents the maximum allowable radial clearance deformation and δ_radial(x) is the actual radial clearance deformation, where g(x)≥0 indicates that the blade radial clearance is within the safe range and g(x)<0 indicates structural failure.

Since the limit state function (LSF) is usually highly nonlinear, the direct calculation of failure probability will consume a lot of computing resources. Monte Carlo (MC) simulation estimates the probability of failure by random sampling, and its basic idea is to calculate failure events through a large number of samples. Specifically, the failure probability can be expressed as:

$${\widehat{P}}_f\approx {\displaystyle \frac{1}{N}}{\displaystyle {\sum }_{i=1}^{N}I\left[g\left(X_i\right)\le 0\right]}$$

(19)

In this paper, MC simulation is used to determine the radial clearance deformation distribution of compressor blades, and hence, the reliability degree is expressed as

$$\left\{\begin{array}{l}R=P(X\ge 0)={\displaystyle \int I}(g(x))g(x)dx={\displaystyle \frac{1}{N}}{\displaystyle \sum I}(g(x))={\displaystyle \frac{N_r}{N}}\hfill \\ \mathrm{s}.\mathrm{t}.I(g(x))=\left\{\begin{array}{l}0,g(x)\le 0\hfill \\ 1,g(x)\ge 0\hfill \end{array}\right.\hfill \end{array}\right.$$

(20)

in which R is the degree of reliability; I(·) is the indicator function of the margin of safety; N_r is the number of samples located in the margin of safety; and N indicates the number of total samples.

The sensitivity index S_i and importance index I_i are defined as

$$\left\{\begin{array}{l}{S}_i={\displaystyle \raisebox{1ex}{${\left(\frac{x_i-{\mu }_x}{{\sigma }_x}\right)}^2$}\!\left/ \!\raisebox{-1ex}{${\beta }^2$}\right.}\hfill \\ I_i={\displaystyle \raisebox{1ex}{$\left|S_{d,i}\right|$}\!\left/ \!\raisebox{-1ex}{${\displaystyle \sum _{i=1}^m\left|S_{d,i}\right|}$}\right.}\hfill \end{array}\right.$$

(21)

where x_i is the value of the i-th input parameter; µ_i is the mean value of the i-th samples; β is a normalizing factor for the sensitivities; S_d,i is the sensitivity index of the i-th samples; and m is the total number of input samples.

4. Experiments and Results

4.1. Experiment Environment

The experimental setup utilized a workstation featuring an Intel 13000KF processor, NVIDIA GeForce RTX 4090 graphics card, 64 GB memory, and 2 TB NVMe solid-state drive for rapid data retrieval. The computing environment operated on Windows 11 with Python 3.9 deployed through Anaconda distribution.

4.2. Deterministic Simulation of Blade with Angle Crack

A comprehensive three-dimensional finite element model of the compressor blade with angle crack defects was established. The angle crack is located at the intersection between the blade tip edge and the suction, where complex multiaxial stress states occur during operation, including centrifugal tension, aerodynamic bending loads, and thermal gradients. Cracks of varying lengths are simulated by adjusting the contact mode between the diminutive triangular blade block and the remaining blade structure. The spatial relationship between these small triangular blade blocks and the intact blade segments is illustrated in Figure 3(a), while the mesh configuration employed in the simulation is depicted in Figure 3(b). For blade model with angle crack, the number of nodes and cells for remaining blade and small triangular blades are, 12,315 and 8,321, and 1,273 and 963, respectively.

The simulation of the angular crack is achieved by manipulating the contact approach. A binding contact is established on a no-crack surface, while a frictionless contact is implemented directly at the crack surface. The length and position of the crack are controlled by altering the relative placement of binding and frictionless. To enhance the simulation efficiency, the crack surface is subdivided into some cells, as shown in Figure 4. Frictionless contact is applied between some cells surface on the triangular blade and the upper contact surface of remaining blade, while binding contact is employed for the other cell surfaces of the triangular blade.

This study examines a turbine compressor blade with constrained root conditions, accounting for centrifugal forces generated by high-speed rotation. Figure 5 shows the geometric structure and boundary conditions of the turbine blade. Convective heat transfer boundaries are applied across four distinct regions: root, lower, middle, upper sections, and tip area. Each region is characterized by specific temperature values (T₁, T₂, T₃, T₄) and corresponding heat transfer coefficients (a₁, a₂, a₃, a₄). Nine stochastic parameters follow normal distribution patterns, detailed in Table 1, with numerical designation 1-4 representing exterior-to-interior positioning. Besides, the distribution of crack lengths is included in Table 1. Finite element modeling incorporates these boundary specifications, yielding results presented in Figure 6 that serve as baseline data for subsequent reliability assessment.

4.3. Structural Reliability Assessment of Compressor Blade Angle Crack

As shown in Table 2, it is the relationship between Monte Carlo simulation sample sizes and reliability degrees. As the MCS increases from 10² to 10⁵, the reliability improves from 0.9512 to 0.9713, indicating convergence around 0.9713 when sample sizes reach 10⁴ and above. This demonstrates that larger sample sizes yield more stable reliability estimates.

Based on the probability distribution of input parameters, MC simulation is used to input 10⁴ sets of randomized parameters into the established radial deformation analysis model to obtain clearance distribution. The obtained radial deformation is displayed in Figure 7, which follows a normal distribution with a mean value of 0.735 mm and a standard deviation of 0.037 mm. Figure 8 shows the distribution characteristics of the samples.

Table 3 quantifies the sensitivity indices (S) and importance indices (I) for ten critical parameters. Parameter 1 (ω) exhibits the highest sensitivity (S=0.42) and importance (I = 0.37), confirming rotational torque as the most influential factor. Parameters 2, 6, and 10 demonstrate moderate-to-high sensitivity values (S=0.36-0.39), while parameters 5, 8, and 9 show relatively low sensitivity (S=0.09-0.17). Figure 9 visualizes the sensitivity index distribution across all parameters, with color-coded bars representing different parameter categories. The results establish a clear hierarchy: high-sensitivity parameters (ω, T₁), moderate-sensitivity parameters (T₂, T₃, σ₁, σ₂), and low-sensitivity parameters (ε₃, ε₄). This quantitative ranking provides essential guidance for design optimization and manufacturing tolerance allocation in fatigue life prediction applications.

4.4. Modeling Efficiency and Precision Analysis

Training data were generated using Latin Hypercube Sampling (LHS) to ensure uniform coverage across the multi-dimensional uncertainty parameter space. A total of 300 training samples were generated, with each sample corresponding to a complete multi-physics finite element analysis calculating the blade’s radial deformation response under specific parameter combinations. Data preprocessing included normalization of all input variables to [0,1] interval, outlier detection using 3σ criterion, and dataset partitioning with 7:2:1 ratio for training, validation, and testing sets.

The error metrics are used in this paper, including Root Mean Square Error (RMSE), Mean Absolute Error (MAE) and Coefficient of Determination (R²), which are defined as follows [3]:

$$\mathbf{RMSE}=\sqrt{{\displaystyle \frac{1}{n}}{\displaystyle \sum _{i=1}^n{\left({\widehat{y}}_i-y_i\right)}^2}}$$

(22)

$$\mathbf{MAE}={\displaystyle \frac{1}{n}}{\displaystyle \sum _{i=1}^n\left|y_i-{\widehat{y}}_i\right|}$$

(23)

$${\mathbf{R}}^2=1-{\displaystyle \frac{{\displaystyle \sum _{i=1}^n{(y_i-{\widehat{y}}_i)}^2}}{{\displaystyle \sum _{i=1}^n{(y_i-\overline{y})}^2}}}$$

(24)

where y_i is the i-th radial deformation of true samples; n is total number of samples; ŷ_i is predicted value of the response variable for the i-th observation; ȳ is arithmetic mean of all measured values.

Table 4 presents a comprehensive performance comparison of different surrogate models across multiple evaluation metrics. The proposed OOA-KR method demonstrates superior performance with the highest R² value of 0.8842 and the lowest prediction errors (RMSE: 0.568, MAE: 0.412). Among the baseline methods, PSO-Kriging achieves the best performance (R²: 0.8616), followed by GA-Kriging (R²: 0.8445) and Kriging (R²: 0.8313), while ANN shows the poorest accuracy (R²: 0.6723). Regarding computational efficiency, basic methods (RBF, Kriging) require minimal training time (6.28-8.94s), whereas optimization-based approaches demand significantly longer training periods (52.82-88.31s) due to their iterative search processes. All models exhibit fast prediction times (0.013-0.043s), making them suitable for real-time applications. The results indicate that the proposed OOA-KR method achieves an optimal balance between prediction accuracy and computational efficiency.

Table 5 demonstrates the simulation precision of reliability estimation across different sample sizes. The proposed OOA-KR method achieves the highest simulation precisions of 91.9%, 96.3%, and 97.6% for sample sizes of 10², 10³, and 10⁴, respectively, significantly outperforming baseline methods. Among conventional approaches, PSO-Kriging exhibits the best performance (87.9%-90.7%), while ANN shows the lowest accuracy (69.7%-76.4%). All methods demonstrate improved precision with increasing sample sizes, with the proposed OOA-KR method maintaining superior performance across all sample configurations, indicating its effectiveness for reliability assessment applications.

5. Conclusions

This study proposes a novel OOA-optimized Kriging-RBF (OOA-KR) method for efficient reliability assessment of compressor blade radial clearance considering angle crack defects. The comprehensive investigation leads to the following key conclusions:

(1) Angle cracks significantly alter blade structural behavior, with crack simulation achieved through contact manipulation between triangular blade blocks and remaining structures, effectively capturing crack-induced compliance changes.

(2) The proposed OOA-KR hybrid model achieves superior prediction accuracy (RMSE=0.568, R²=0.8842), outperforming individual Kriging, RBF, and other optimization-based surrogate models by significant margins.

(3) Under normal operational conditions, the blade achieves a reliability of 0.9713 with mean radial deformation of 0.735 mm and standard deviation of 0.037 mm based on Monte Carlo simulation with 10⁴ samples.

(4) Global sensitivity analysis identifies rotational speed as the most critical factor (S=0.42, I=0.37), followed by temperature parameters and heat transfer coefficients, providing quantitative guidance for design priorities.

(5) The OOA-KR method demonstrates excellent computational efficiency with training time of 52.82s and prediction time of 0.025s, achieving optimal balance between accuracy and computational cost.

(6) The methodology achieves reliability assessment precision up to 97.6% for large sample sizes, significantly superior to conventional approaches and suitable for practical engineering applications.

The proposed approach demonstrates superior performance in capturing nonlinear multi-physics effects while maintaining computational efficiency, making it valuable for aerospace reliability applications. Future work will focus on extending the methodology to dynamic loading conditions and multiple crack interaction scenarios.