Graphical Abstract Figure
Graphical Abstract Figure
Close modal

Abstract

The precision in forming complex double-walled hollow turbine blades significantly influences their cooling efficiency, making the selection of appropriate casting process parameters critical for achieving fine-casting blade formation. However, the high cost associated with real blade casting necessitates strategies to enhance product formation rates and mitigate cost losses stemming from the overshoot phenomenon. We propose a machine learning (ML) data-driven framework leveraging an enhanced whale optimization algorithm (WOA) to estimate product formation under diverse process conditions to address this challenge. Complex double-walled hollow turbine blades serve as a representative case within our proposed framework. We constructed a database using simulation data, employed feature engineering to identify crucial features and streamline inputs, and utilized a whale optimization algorithm-back-propagation neural network (WOA-BP) as the foundational ML model. To enhance WOA-BP’s performance, we introduce an optimization algorithm, the improved chaos whale optimization-back-propagation (ICWOA-BP), incorporating cubic chaotic mapping adaptation. Experimental evaluation of ICWOA-BP demonstrated an average mean absolute error of 0.001995 mm, reflecting a 36.21% reduction in prediction error compared to conventional models, as well as two well-known optimization algorithms (particle swarm optimization (PSO), quantum-based avian navigation optimizer algorithm (QANA)). Consequently, ICWOA-BP emerges as an effective tool for early prediction of dimensional quality in complex double-walled hollow turbine blades.

1 Introduction

The aero-engine is the most critical part of an airplane, and its performance largely depends on the turbine inlet temperature. In the case of limited material temperature resistance, proposing a new blade cooling structure represented by thin-walled, double-walled blades is of great significance for improving engine performance. Currently, research on investment casting of turbine blades has reached a relatively advanced stage, with investment casting being extensively applied in the high-precision molding of intricate structures. Nonetheless, accurately predicting the intricate casting of turbine blades presents a formidable challenge. This challenge arises not only from the process conditions during casting [1], but also from the complex internal structural deformation and external contour. Given the prohibitively high cost of casting double-walled turbine blades, there is a pressing demand for novel methodologies to forecast casting formation. In this work, the modeling of the accuracy prediction of complex structures inside a double-walled turbine blade is used as a case study. Figure 1 shows the complex cooling structure explored.

Fig. 1
Double-walled turbine blade
Fig. 1
Double-walled turbine blade
Close modal

Numerical simulation of the casting process with fluid-thermal-solid coupling is widely used to predict the thermal deformation of various materials [2,3]. Numerical simulation of the casting process is the basis for high-quality, short-cycle, and low-cost manufacturing of new-generation net-forming precision-cast turbine blades. However, due to the long computational cycles required for numerical simulation and the poor access to data sets, it is unsuitable for rapid, real-time evaluation. Therefore, more efficient prediction methods need to be constructed. In today’s era of Industry 4.0 and digitization of manufacturing [4], new communication technologies combined with state-of-the-art algorithms in the field of machine learning (ML) and deep learning (DL) offer a wide range of applications for data-driven smart manufacturing [57]. ML and DL methods are used to automate data-driven quality analytics by enabling them, and their application in this field is known as quality prediction [8,9].

Various types of sensors monitor the casting environment, and ML methods utilize the numerical simulation data of fluid-thermal-solid coupling for casting defect prediction [1012], as well as prediction of product dimensions and quality [1315].

To further reduce the cost of time for a large number of repetitive simulations, small samples of simulation results are used to train neural network models to quickly predict geometric deviations [16] and stress–strain of products under different parametric conditions [17]. Meanwhile, ML can be used in the additive manufacturing stage [18], to predict the quality of porosity [19,20], tensile strength [21], surface roughness [22], microstructural defects [23], and molding quality of the product [2426].

At present, researchers globally try to use intelligent algorithms for regression prediction modeling of the cold-effect design of a turbine blade [27], and regression prediction modeling based on traditional neural networks provides a new solution to the above problems and has achieved specific research results [28,29]. However, researchers have not yet carried out relevant work using intelligent algorithms in the blade casting link, and the simple use of traditional neural networks still has certain deficiencies in the prediction accuracy and endpoint hit rate. To solve the above problems, our research group combined the improved whale swarm optimization algorithm to improve the traditional neural network. The whale optimization algorithm (WOA), as a new swarm intelligent optimization algorithm, still has some deficiencies but has the potential for development. Thus, it is crucial to improve this algorithm and expand its application areas [30,31].

In this paper, an ML-based method for predicting the early dimensional accuracy of the typical structure of a complex double-walled hollow turbine blade is proposed. The framework and key algorithms of the method are investigated, and the effectiveness of the method is demonstrated through a verification using a complex double-walled hollow turbine blade. Figure 2 illustrates the overall technical route of this paper. The main contributions of this paper are as follows:

  • A data-driven framework for predicting the early dimensional accuracy of typical structures of a complex double-walled hollow turbine blade is proposed;

  • The performance of different ML algorithms is investigated, and the results show that the ICWOA algorithm exhibited better prediction accuracy;

  • The potential relationship between the casting conditions with different process parameters and the deformation of the typical structure of a complex double-walled hollow turbine blade is revealed through feature engineering techniques.

Fig. 2
Research for early mass prediction of complex double-walled turbine blades
Fig. 2
Research for early mass prediction of complex double-walled turbine blades
Close modal

2 Whale Optimization Algorithm

Whale optimization Algorithm simulates the humpback whale bubble-net attacking method, which is more straightforward to operate and converges faster than other optimization algorithms. The standard WOA mainly relies on the coefficient vector, A, to select the path to search for the prey and uses the probability, P, to determine the final predation mechanism. The computational flow of the standard WOA is shown in Fig. 3.

Fig. 3
WOA algorithm flow
Fig. 3
WOA algorithm flow
Close modal

WOA includes two mechanisms, shrinking encircling and spiraling updating, to model the unique bubble-net hunting strategy of humpback whales; in addition, humpback whales also randomly search for prey. The following section briefly introduces the mathematical model of WOA.

2.1 Encircling Prey.

Humpback whales can identify the location of their prey and surround them. The location of the true optimal solution in the target space is not a priori. Moreover, WOA assumes that the optimal solution for the current population is the location of the prey, denoted by X*, and that after X* has been determined, the other individuals move toward X* in the following manner:
(1)
(2)
where t is the number of generations of the current iteration, A and C are the coefficients, X*(t) is the position vector of the current optimal solution (i.e., the location of the prey), and X(t) is the position vector of the current solution. X* needs to be updated each time a better solution appears during each iteration, and A and C are calculated as follows:
(3)
(4)
where α is a coefficient that decreases linearly with iterative generations from 2 to 0, α=22tTmax, t is the current number of iterations, Tmax is the maximum number of iterations, and r1 and r2 are random number between [0,1].

2.2 Bubble-Net Attacking Method.

WOA includes two movement mechanisms to simulate the bubble-net hunting strategy of the humpback whale population. The first is a shrinking circle encirclement, which linearly decreases the value of α in Eq. (3); the another is a spiral update, which is mathematically modeled as follows:
(5)
(6)
where D is the distance between the current search individual and the current optimal solution, b is the spiral shape parameter, and k is a [1,1] uniformly distributed random number.
Since two types of predation behaviors exist in the process of approaching prey, WOA chooses either bubble-net predation or contraction encirclement, according to the probability p. The position update formula is as follows:
(7)
where p is the probability of the predation mechanism and the random number with the value range [0,1]. As the number of iterations t increases, the parameter A and the convergence factor α gradually decrease, and if |A|<1, each whale gradually surrounds the current optimal solution, which belongs to the local optimal-seeking stage in WOA.

2.3 Searching for Prey.

WOA simulates the behavior of humpback whale pods searching for prey locations in addition to their approach to prey, as expressed by the following equation:
(8)
(9)
where Xrand is a random selection of individuals in the current population, and in this movement mode, humpback whales search randomly according to each other’s position. WOA uses the numerical magnitude of |A| to determine whether the population of whales is conducting a shrinking circle encirclement or searching for prey; when |A|<1, the population conducts a shrinking circle encirclement, and when |A|1, the population conducts a search for prey.

3 Improvements to Whale Optimization Algorithm

Although WOA has good optimization-seeking capability, the basic WOA still has the disadvantages of low solution accuracy, slow convergence, and a tendency to fall into the local optimum. In the population intelligence optimization algorithm, randomly generated arrays are required; however, the results are often poor. Since chaotic sequences are nonlinear, ergodic, and unpredictable, random arrays can be replaced by chaotic mappings. To overcome these shortcomings, we improve the WOA in three aspects: population initialization, position update strategy, and prevention of falling into the local optimum. The improved algorithm is ICWOA, and the overall flowchart is shown in Fig. 4.

Fig. 4
ICWOA algorithm flow
Fig. 4
ICWOA algorithm flow
Close modal

3.1 Cubic Map.

The standard Cubic mapping [32] function can be expressed as shown in Eq. (10)
(10)
Here b and c are chaotic influence factors, and the range of cubic mapping is different for different values of b and c. Generally, when c(2.3,3), the sequence generated by cubic mapping is chaotic. In addition, xn(2,2) when b=1, and xn(1,1) when b=4. The maximum Lyapunov exponent of 16 commonly used 1D chaotic mappings was calculated and analyzed in the literature [33], and the cubic mapping expression is shown in Eq. (11)
(11)
where xn(0,1); ρ is the control parameter, and the chaotic property of the cubic mapping has a strong relationship with the value of parameter ρ. From Ref. [33], the chaotic property of cubic mapping is known to be similar to the maximum Lyapunov exponent of the logistic mapping and tent mapping, and it is better than 1D mappings such as sine, circle, Singer, and Kent mapping. An improved cubic mapping with a maximum Lyapunov exponent of 1.0984 was proposed in the literature [34], which is better than the Lyapunov exponent of 0.6930 for the logistic mapping.

In this paper, the distribution curves of the values taken after 103 iterations of the cubic chaos mapping are given. The initial value of xn(0,1), and the cubic mapping has good chaos traversal when x0=0.3 and ρ=2.59. The simulation results are shown in Fig. 5.

3.2 Improved Nonlinear Convergence Factor.

Since the convergence factor for linear variation does not well regulate the global search ability and local exploitation ability, a nonlinear convergence factor [3537] is proposed in this paper as follows:
(12)
where afirst and afinal are the initial and final values of the control parameters, respectively; max_iter is the maximum number of iterations; t is the current number of iterations; μ and φ are the regulation parameters.

3.3 Adaptive Weighting Strategy and Stochastic Difference Method Variation Strategy.

WOA is prone to falling into the local optimum and early convergence in late local development. To enable the algorithm to maintain the population diversity and to jump out of local optimum in time, we propose an adaptive weight strategy [38,39] and random difference variation strategy [40,41].

The mathematical expressions of the adaptive weighting strategy are as follows:

(13)
(14)
The stochastic difference variance strategy is as follows:
(15)
where r1 and r2 are both random numbers of [0,1]; X(t) is a randomly selected individual in the population.

Each individual must go through the encircling predation, spiral update, and prey search stages, and the adaptive weighting strategy is used to update the position when the individual performs encircling predation or spiral update. Afterward, the individual needs to update it again by the random difference variation strategy to take the optimal position before and after its change, which speeds up the convergence of the population and effectively prevents the population from falling into the local optimum. The population works in concert by these two strategies, giving the algorithm a better optimality-finding effect.

4 Data Processing of Critical Parts of a Double-Walled Blade

For data-driven prediction algorithms, the quality of the dataset is critical. In this study, the deformation of structures such as impact holes with different aperture sizes, rotary section, trailing-edge cutback, and pin fin were the key structures affecting the cooling efficiency of the blade [4244]. From the perspective of double-wall blade casting, among the effects of different process parameters on the deformation of single-crystal nickel-based high-temperature alloys during the solidification stage, the pouring temperature [45], preheating temperature, and drawing speed had a great influence on the deformation of casting complex structures of the blade. Therefore, we processed the data for the numerical simulation casting deformation of these structures under different process parameters and used them as the initial input features for the prediction model.

4.1 Impact-Hole Error Statistics.

Since the procast software exports the model as an stereolithography (STL) point cloud file, the point cloud was extracted locally for the impact-hole area, a circle was fitted to the 3D discrete points, with known coordinates of the three-dimensional spatial discrete points (xi,yi,zi), and a spatial circle was constructed so that the spatial points were as close as possible to the fitted spatial circle.

First, all the discrete points were on a plane as much as possible, and the plane equation can be expressed as
(16)
The above can be written in matrix form as
(17)
where

Calculation of the error can be done using the Euclidean paradigm: MAL1=i=1n(MAL1)2.

This is a hyper-deterministic equation to solve for A=(MTM)1MTL1, which is the normal vector of the plane, according to the method of least squares.

Assuming that all discrete points are on the circle, the middle of the line connecting any two points had to pass through the center of the circle. Let the center of the circle be C(x0,y0,z0). We took two points, P1(x1,y1,z1) and P2(x2,y2,z2). Then, the vector1 of the line connecting P1 and P2 is expressed as (x2x1,y2y1,z2z1), and the coordinates P12 of the midpoint of the line connecting P1 and P2 are (x1+x22,y1+y22,z1+z22).

The center of the circle C was connected to P12 with vector2 as (x1+x22x0,y1+y22y0,z1+z22z0). For P1 and P2 to be on the circle, vector1*vector2=0 had to be satisfied, i.e.,
(18)
Thus, it was possible to obtain
(19)
where Δx12=x2x1, Δy12=y2y1, Δz12=z2z1,l1=x22+y22+z22x12y12z122.
All points were on the circle; then, we had
(20)
where the above equation is also an overdetermined equation. Calculation of the error can be done using the Euclidean paradigm: BCL2=i=1n(BCL2)2. Since the center of the circle C had to be in the plane of the aforementioned control, it satisfied ax0+by0+cz01=0, i.e.,
(21)
where A is the normal vector of the plane, which was found by A=(MTM)1MTL1. Therefore, it could be solved by constructing an optimization problem under the constraints of Eq. (20)
(22)
where λ is the Lagrangian multiplier. When we took the derivative concerning C with respect to λ and made the derivative value 0, the transformation yielded
(23)
The radius of a circle could be determined from the average of the distances from all points to the center of the circle
(24)

Therefore, using the STL point cloud data, we constructed a spatial circle so that the spatial points were as close as possible to the fitted spatial circle, and according to the least squares method, we could derive the normal vector of the plane. According to the point coordinates, to find the center of the spatial circle and, at the same time, get the fitted radius, we coded and solved the derivation process of the equation of the fitted circle curve using matlab, as shown in Fig. 6, to obtain the simulation mode with l 1.0-mm, 1.2-mm, 1.3-mm, and 1.5-mm impact-hole circle fitting.

Fig. 6
Impact-hole circle fitting: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), and (d) impact hole (1.5 mm)
Fig. 6
Impact-hole circle fitting: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), and (d) impact hole (1.5 mm)
Close modal

Figure 7 shows the error probability density function (PDF) of the variation amount of four different orifice impact holes. The dynamic leaf impact-hole apertures were 1.0 mm, 1.2 mm, 1.3 mm, and 1.5 mm, and the variation of 1.0-mm impact-hole aperture was 0.10–0.15 mm; the 1.2-mm impact-hole aperture variation was 0.06–0.20 mm; the 1.3-mm impact-hole aperture variation was 0.11–0.21 mm, and the 1.5-mm impact-hole aperture variation was 0.16–0.30 mm.

Fig. 7
Four types of impact-hole aperture variation
Fig. 7
Four types of impact-hole aperture variation
Close modal

Although the statistical mean and standard deviation of the overall geometric deviation parameters on the impact holes of different bore diameters differed, their probability density functions were close to the Gaussian distribution, and the statistical results are shown in Fig. 8. The mean value of the 1.0-mm impact holes was 0.12142, with a standard deviation of 0.00919. The mean value of the 1.2-mm impact holes was 0.13033, with a standard deviation of 0.02939. The mean value of the 1.3-mm impact holes was 0.16109, with a standard deviation of 0.02097. Further, the mean value of 1.5-mm impact holes was 0.21790, with a standard deviation of 0.02985. The statistical mean value of the variation of different impact-hole diameters gradually increased, and the standard deviation exhibited a decreasing trend when the hole diameter was 1.3 mm, and it started to increase again when the hole diameter was 1.5 mm. The variation of the impact-hole diameter directly affected the cooling efficiency in its working condition.

Fig. 8
Statistical results for different impact holes
Fig. 8
Statistical results for different impact holes
Close modal

4.2 Rotary Section Error Statistics.

Figure 9 shows the error PDF of the radius of curvature of the rotary section. The original design model, chamfer, of the front chamber rotary section of the moving leaf was 0.4–0.8 mm. Because the model shrinks irregularly after the simulation casting, it does not form a complete chamfer angle, so the radius of curvature was used here instead of the chamfer to indicate the change law. Moreover, the radius of curvature of the chamfer of the front chamber rotary section increased after the shrinkage deformation, the deformation of the front chamber rotary section area was larger, and the radius of curvature was mainly concentrated in 8–9.5 mm. The original design model chamfer of the dynamic lobe middle chamber rotary section was 0.4 mm, the radius of curvature of the chamfer of the front chamber rotary section increased after shrinkage deformation, and the radius of curvature was concentrated in 8.5–10 mm.

Fig. 9
Curvature of rotary section
Fig. 9
Curvature of rotary section
Close modal

The statistical mean and standard deviation of the geometric deviation parameters of the front chamber rotary section and the middle chamber rotary section were different, and the statistical results are shown in Fig. 10. The mean value of the anterior chamber was 8.99706, with a standard deviation of 0.69427; the mean value of the middle chamber was 9.29326, with a standard deviation of 0.62595. The mean value of the middle chamber rotary section was higher than that of the front chamber rotary section. Still, the standard deviation of the middle chamber rotary section was lower than that of the front chamber rotary section, revealing that the curvature variation of the middle chamber rotary section was more concentrated.

Fig. 10
Statistical results of the front and middle chamber rotary section
Fig. 10
Statistical results of the front and middle chamber rotary section
Close modal

4.3 Trailing-Edge Cutback and Pin Fin Error Statistics.

Figure 11 shows the error PDF of the trailing-edge cutback and pin fin. Here, the maximum spacing of the trailing-edge cutback position was taken to indicate its variation law. The spacing size of the trailing-edge cutback in the original design model was 0.4 mm, and the spacing increased at the trailing-edge cutback after shrinkage deformation, mainly concentrated in 0.1–0.18 mm. A Gaussian distribution for the trailing-edge cutback was fitted to the distribution, where the mean was 0.14839, and the standard deviation was 0.0302. The diameter of the original design model of the pin fin was 1.2 mm, and here, the same circle fitting as the impact hole was used to fit the cylindrical cross section of the pin fin to obtain the diameter of the pin fin after casting deformation. The specific method is introduced in Sec. 4.1, so we will not repeat it here. The radius of the pin fin after shrinkage deformation changed in the range of 0.02–0.20 mm, which was mainly concentrated in 0.04–0.16 mm. A Gaussian distribution for pin fin was fitted to the distribution, where the mean was 0.09754, and the standard deviation was 0.03934.

Fig. 11
Curvature of trailing-edge cutback and pin fin
Fig. 11
Curvature of trailing-edge cutback and pin fin
Close modal
After processing the initial feature input data, and considering the volatility of the investment casting process, as well as the fact that the deformation deviation follows the characteristics of a normal distribution, the initial data are normalized to the database using the Zscore method, as shown below:
(25)
where x is the original data, μ is the mean of x, and σ is the standard deviation. In this way, the bias index was scaled to zero mean and unit variance.

5 Comparison of Prediction Model Results

The dataset with normalization done in Sec. 4 is used as input to train the early neural network model. To illustrate the performance of the ICWOA algorithm, in addition to comparing the ICWOA algorithm with the WOA algorithm, the ICWOA algorithm is also compared with the PSO algorithm. To better compare the differences between the different algorithms, each algorithm has the same setup conditions, where the population size and the maximum number of iterations are set to 30 and 1000, respectively. Figure 12 depicts the critical structure prediction modeling framework.

Fig. 12
Modeling process of critical structure prediction modeling with ICWOA
Fig. 12
Modeling process of critical structure prediction modeling with ICWOA
Close modal

Convergence performance is one of the important indicators of optimization algorithms. To evaluate the ICWOA algorithm, we show the convergence curves of different thickness prediction models in Fig. 13. It can be seen that the ICWOA-BP model converged faster than the WOA-BP model, and at the same time, the prediction error was better. This indicates that the WOA algorithm may fall into local optimal solutions. In addition, the fast iteration speed was also proved to be another advantage of ICWOA. From these results, it can be concluded that ICWOA outperforms WOA in global optimization for the same amount of data.

Fig. 13
Convergence curves of the ICWOA and WOA in different test sets
Fig. 13
Convergence curves of the ICWOA and WOA in different test sets
Close modal

A detailed comparison between ICWOA-BP, WOA-BP, and PSO-BP predicted deformations and simulation results is shown in Fig. 14, and it is obvious that the WOA-BP prediction accuracy was significantly improved after the improvement. From Figs. 14(b)14(d), WOA-BP as well as PSO-BP were caught in the phenomenon of the local optimal solution, and large deviations appeared locally. In contrast, the optimized ICWOA-BP did not have a large deviation point, and the maximal difference between the actual value and the predicted value was stabilized at ±0.005mm under the data of the deformation of four different impact holes, except that the improved ICWOA-BP performed slightly better than PSO-BP. In addition, the fast iteration speed was another advantage of ICWOA. In the validation set, the average MAE of the optimized ICWOA-BP model was 1.955×103, a 36.21% reduction in the error compared with the original model WOA-BP. The result indicates that the prediction accuracy of the ICWOA model is satisfactory.

Fig. 14
Predicted deformation bias of ICWOA, WOA, and PSO in different ML models: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), (d) impact hole (1.5 mm), (e) front chamber rotary section, (f) middle chamber rotary section, (g) trailing edge cutback, and (h) pin fin
Fig. 14
Predicted deformation bias of ICWOA, WOA, and PSO in different ML models: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), (d) impact hole (1.5 mm), (e) front chamber rotary section, (f) middle chamber rotary section, (g) trailing edge cutback, and (h) pin fin
Close modal

To further demonstrate the advantages and usefulness of the algorithms, the ICWOA is compared with the quantum-based avian navigation optimizer algorithm (QANA) [46]. Inspired by the navigation mechanism of migratory birds, Zamani et al. proposed a robust and scalable differential evolutionary algorithm. QANA is very competitive state-of-the-art the swarm intelligence (SI) and the differential evolution (DE) algorithms and was awarded the best theory paper in Engineering Applications of Artificial Intelligence from 2020 to 2022. It is also proven to be accurately applicable to solving real-world engineering optimization problems. The QANA algorithm is due to the WOA algorithm in terms of accuracy and convergence speed [47], so it is compared with the ICWOA algorithm to verify the optimization effect. A detailed comparison between ICWOA-BP and QANA-BP predicted deformations and simulation results is shown in Fig. 15. In some feature dimensions, such as Figs. 15(a), 15(c), and 15(d), the QANA algorithm and the ICWOA algorithm do not differ much in their predictions. In other feature dimensions, the ICWOA algorithm is more stable and has a higher prediction accuracy, and the QANA algorithm appears to have a large deviation.

Fig. 15
ICWOA and QANA comparison: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), (d) impact hole (1.5 mm), (e) front chamber rotary section, (f) middle chamber rotary section, (g) trailing edge cutback, and (h) pin fin
Fig. 15
ICWOA and QANA comparison: (a) impact hole (1.0 mm), (b) impact hole (1.2 mm), (c) impact hole (1.3 mm), (d) impact hole (1.5 mm), (e) front chamber rotary section, (f) middle chamber rotary section, (g) trailing edge cutback, and (h) pin fin
Close modal

In addition, the improved ICWOA model and the comparison of MAE results of different models are shown in Table 1. It can be seen that there are some effects of prediction bias in the data with different features, but the prediction accuracy of the ICWOA-BP model is the highest among all the models, which is improved by about 22.05% on average compared with PSO-BP, and the performance of ICWOA-BP is better than PSO-BP. Meanwhile, the average MAE improved by 14.29% compared to QANA-BP, and ICWOA-BP was superior to QANA-BP.

Table 1

Performance evaluation of the ML models

MAE in testing phase
Key parts of the bladeBPPSO-BPWOA-BPICWOA-BPQANA-BP
Impact hole (1.0 mm)5.696 × 10−31.217 × 10−31.976 × 10−31.124 × 10−31.249 × 10−3
Impact hole (1.2 mm)6.162 × 10−32.954 × 10−33.078 × 10−31.791 × 10−32.083 × 10−3
Impact hole (1.3 mm)6.420 × 10−32.479 × 10−32.726 × 10−31.406 × 10−31.616 × 10−3
Impact hole (1.5 mm)6.507 × 10−32.511 × 10−33.229 × 10−32.353 × 10−32.801 × 10−3
Front chamber rotary section6.485 × 10−32.177 × 10−33.856 × 10−31.940 × 10−32.180 × 10−3
Middle chamber rotary section7.462 × 10−33.658 × 10−33.214 × 10−32.567 × 10−33.209 × 10−3
Trailing-edge cutback6.629 × 10−32.211 × 10−33.056 × 10−32.131 × 10−32.449 × 10−3
Pin fin9.789 × 10−32.870 × 10−33.387 × 10−32.331 × 10−32.664 × 10−3
Average6.894 × 10−32.510 × 10−33.065 × 10−31.955 × 10−32.281 × 10−3
MAE in testing phase
Key parts of the bladeBPPSO-BPWOA-BPICWOA-BPQANA-BP
Impact hole (1.0 mm)5.696 × 10−31.217 × 10−31.976 × 10−31.124 × 10−31.249 × 10−3
Impact hole (1.2 mm)6.162 × 10−32.954 × 10−33.078 × 10−31.791 × 10−32.083 × 10−3
Impact hole (1.3 mm)6.420 × 10−32.479 × 10−32.726 × 10−31.406 × 10−31.616 × 10−3
Impact hole (1.5 mm)6.507 × 10−32.511 × 10−33.229 × 10−32.353 × 10−32.801 × 10−3
Front chamber rotary section6.485 × 10−32.177 × 10−33.856 × 10−31.940 × 10−32.180 × 10−3
Middle chamber rotary section7.462 × 10−33.658 × 10−33.214 × 10−32.567 × 10−33.209 × 10−3
Trailing-edge cutback6.629 × 10−32.211 × 10−33.056 × 10−32.131 × 10−32.449 × 10−3
Pin fin9.789 × 10−32.870 × 10−33.387 × 10−32.331 × 10−32.664 × 10−3
Average6.894 × 10−32.510 × 10−33.065 × 10−31.955 × 10−32.281 × 10−3

6 Conclusion

In this paper, a data-driven framework based on ML techniques was developed for estimating and screening early products. A complex double-walled hollow turbine blade was analyzed as a typical case of the proposed framework. Owing to the extremely high cost of real blade casting, this paper adopts WOA-BP neural network as the basic model of the machine learning algorithm, and proposes a cubic chaos mapping adaptive optimization algorithm ICWOA-BP based on the compilation of the database using the simulation data, and adopts the characteristic structure of the complex double-wall hollow turbine blade as the initial input, and simplifies the data processing. The main findings of this paper are as follows:

  1. A data-driven framework for early dimensional accuracy prediction of complex double-walled hollow turbine blades is proposed, and the effectiveness of the method is illustrated with an example of a critical part of a complex double-walled hollow turbine blade;

  2. We found that six feature dimensions provided the best performance. The selection of features in different critical areas had a significant impact on the performance of the predictive model;

  3. The optimized WOA method was used for the hyperparameter optimization of BP, and the prediction results in the small dataset could still maintain good prediction accuracy, showing good improvement. In the validation set, the average MAE of the optimized ICWOA-BP model was 1.955×103, reducing the error by 36.21% compared with the original model. The prediction errors of the improved model on the test set were affected by different feature data, the average improvement was about 22.09% compared with PSO-BP, and the performance of ICWOA-BP was better than that of PSO-BP. Meanwhile, the average MAE improved by 14.29% compared to QANA-BP, and ICWOA-BP was superior to QANA-BP.

Acknowledgment

This work was supported in part by the National Natural Science Foundation of China (Grant Nos. 52475491 and 51705440), in part by the Natural Science Foundation of Fujian Province, China (Grant No. 2019J01044), in part by the Aeronautical Science Foundation of China (Grant Nos. 20170368001 and 20230003068002), in part by the Foundation of State Key Laboratory of Intelligent Manufacturing Equipment and Technology (Grant No. IMETKF2024013), in part by the CAS-Fujian STS Project (Grant No. 2022T3071), and in part by the National Science and Technology Major Project (Grant Nos. J2019-II-0022-0043 and J2019-VII-0013-0153).

Declaration of Competing Interest

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

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Ren
,
S.
,
Bu
,
K.
,
Mou
,
S.
,
Zhang
,
R.
, and
Bai
,
B.
,
2023
, “
Control of Dimensional Accuracy of Hollow Turbine Blades During Investment Casting
,”
J. Manuf. Process.
,
99
(
8
), pp.
548
562
.
2.
Wang
,
D.
,
He
,
B.
,
Liu
,
S.
,
Liu
,
C.
, and
Fei
,
L.
,
2016
, “
Dimensional Shrinkage Prediction Based on Displacement Field in Investment Casting
,”
Int. J. Adv. Manuf. Technol.
,
85
(
1–4
), pp.
201
208
.
3.
Dong
,
Y.
,
Li
,
X.
,
Zhao
,
Q.
,
Yang
,
J.
, and
Dao
,
M.
,
2017
, “
Modeling of Shrinkage During Investment Casting of Thin-Walled Hollow Turbine Blades
,”
J. Mater. Process. Technol.
,
244
(
6
), pp.
190
203
.
4.
Fordal
,
J. M.
,
Schjølberg
,
P.
,
Helgetun
,
H.
,
Skjermo
,
T. Ø.
,
Wang
,
Y.
, and
Wang
,
C.
,
2023
, “
Application of Sensor Data Based Predictive Maintenance and Artificial Neural Networks to Enable Industry 4.0
,”
Adv. Manuf.
,
11
(
2
), pp.
248
263
.
5.
Mohamed
,
O. A.
,
Masood
,
S. H.
, and
Bhowmik
,
J. L.
,
2021
, “
Modeling, Analysis, and Optimization of Dimensional Accuracy of FDM-Fabricated Parts Using Definitive Screening Design and Deep Learning Feedforward Artificial Neural Network
,”
Adv. Manuf.
,
9
(
1
), pp.
115
129
.
6.
Ferreira
,
C.
, and
Gonçalves
,
G.
,
2022
, “
Remaining Useful Life Prediction and Challenges: A Literature Review on the Use of Machine Learning Methods
,”
J. Manuf. Syst.
,
63
(
4
), pp.
550
562
.
7.
Wang
,
J.
,
Ma
,
Y.
,
Zhang
,
L.
,
Gao
,
R. X.
, and
Wu
,
D.
,
2018
, “
Deep Learning for Smart Manufacturing: Methods and Applications
,”
J. Manuf. Syst.
,
48
(
7
), pp.
144
156
.
8.
Deng
,
F.
,
Rui
,
X.
,
Lu
,
S.
,
Liu
,
Z.
,
Sun
,
H.
, and
Volk
,
W.
,
2023
, “
Deep Learning-Based Inline Monitoring Approach of Mold Coating Thickness for Al-Si Alloy Permanent Mold Casting
,”
Int. J. Adv. Manuf. Technol.
,
130
(
11
), pp.
1
9
.
9.
Yang
,
F.
,
Zhao
,
W.
,
Ru
,
Y.
,
Pei
,
Y.
,
Li
,
S.
,
Gong
,
S.
, and
Xu
,
H.
,
2023
, “
Deep Learning Accelerates the Development of Ni-Based Single Crystal Superalloys: A Physical-Constrained Neural Network for Creep Rupture Life Prediction
,”
Mater. Des.
,
232
(
8
), p.
112174
.
10.
Lee
,
J.
,
Lee
,
Y. C.
, and
Kim
,
J. T.
,
2020
, “
Fault Detection Based on One-Class Deep Learning for Manufacturing Applications Limited to an Imbalanced Database
,”
J. Manuf. Syst.
,
57
(
10
), pp.
357
366
.
11.
Ke
,
K.-C.
, and
Huang
,
M.-S.
,
2020
, “
Quality Prediction for Injection Molding by Using a Multilayer Perceptron Neural Network
,”
Polymers
,
12
(
8
), p.
1812
.
12.
Liu
,
Z.
,
Song
,
Y.
,
Tang
,
R.
,
Duan
,
G.
, and
Tan
,
J.
,
2023
, “
Few-Shot Defect Recognition of Metal Surfaces Via Attention-Embedding and Self-Supervised Learning
,”
J. Intell. Manuf.
,
34
(
8
), pp.
3507
3521
.
13.
Obregon
,
J.
,
Hong
,
J.
, and
Jung
,
J.-Y.
,
2021
, “
Rule-Based Explanations Based on Ensemble Machine Learning for Detecting Sink Mark Defects in the Injection Moulding Process
,”
J. Manuf. Syst.
,
60
(
7
), pp.
392
405
.
14.
Ge
,
Z.
,
Song
,
Z.
, and
Gao
,
F.
,
2013
, “
Self-Training Statistical Quality Prediction of Batch Processes With Limited Quality Data
,”
Ind. Eng. Chem. Res.
,
52
(
2
), pp.
979
984
.
15.
Lyu
,
Y.
,
Wen
,
Z.
, and
Chen
,
A.
,
2023
, “
A Novel Transfer Learning Approach Based on Deep Degradation Feature Adaptive Alignment for Remaining Useful Life Prediction With Multi-condition Data
,”
J. Intell. Manuf.
12
), pp.
1
19
.
16.
Zhu
,
Z.
,
Ferreira
,
K.
,
Anwer
,
N.
,
Mathieu
,
L.
,
Guo
,
K.
, and
Qiao
,
L.
,
2020
, “
Convolutional Neural Network for Geometric Deviation Prediction in Additive Manufacturing
,”
Proc. CIRP
,
91
(
8
), pp.
534
539
.
17.
Li
,
L.
, and
Anand
,
S.
,
2020
, “
Hatch Pattern Based Inherent Strain Prediction Using Neural Networks for Powder Bed Fusion Additive Manufacturing
,”
J. Manuf. Process.
,
56
(
8
), pp.
1344
1352
.
18.
Ciccarelli
,
M.
,
Corradini
,
F.
,
Germani
,
M.
,
Menchi
,
G.
,
Mostarda
,
L.
,
Papetti
,
A.
, and
Piangerelli
,
M.
,
2023
, “
SPECTRE: A Deep Learning Network for Posture Recognition in Manufacturing
,”
J. Intell. Manuf.
,
34
(
8
), pp.
3469
3481
.
19.
O’Donnell
,
J.
, and
Yoon
,
H.-S.
,
2024
, “
Determination of Multi-component Failure in Automotive System Using Deep Learning
,”
ASME J. Comput. Inf. Sci. Eng.
,
24
(
2
), p.
021005
.
20.
Zhang
,
B.
,
Liu
,
S.
, and
Shin
,
Y. C.
,
2019
, “
In-Process Monitoring of Porosity During Laser Additive Manufacturing Process
,”
Addit. Manuf.
,
28
, pp.
497
505
.
21.
Zhang
,
J.
,
Wang
,
P.
, and
Gao
,
R. X.
,
2019
, “
Deep Learning-Based Tensile Strength Prediction in Fused Deposition Modeling
,”
Comput. Ind.
,
107
(
5
), pp.
11
21
.
22.
Li
,
Z.
,
Zhang
,
Z.
,
Shi
,
J.
, and
Wu
,
D.
,
2019
, “
Prediction of Surface Roughness in Extrusion-Based Additive Manufacturing With Machine Learning
,”
Rob. Comput.-Integr. Manuf.
,
57
(
6
), pp.
488
495
.
23.
Bartlett
,
J. L.
,
Jarama
,
A.
,
Jones
,
J.
, and
Li
,
X.
,
2020
, “
Prediction of Microstructural Defects in Additive Manufacturing From Powder Bed Quality Using Digital Image Correlation
,”
Mater. Sci. Eng. A.
,
794
(
9
), p.
140002
.
24.
Liu
,
Y.
,
Wu
,
Q.
,
Wang
,
P.
, and
Zhuang
,
W.
,
2023
, “
Forming Process Prediction of a Self-Piercing Riveted Joint in Carbon Fibre Reinforced Composites and Aluminium Alloy Based on Deep Learning
,”
J. Manuf. Process.
,
106
(
11
), pp.
453
464
.
25.
Wu
,
Y.
,
Meng
,
Y.
, and
Shao
,
C.
,
2022
, “
End-to-End Online Quality Prediction for Ultrasonic Metal Welding Using Sensor Fusion and Deep Learning
,”
J. Manuf. Process.
,
83
(
11
), pp.
685
694
.
26.
Mongan
,
P. G.
,
Hinchy
,
E. P.
,
O’Dowd
,
N. P.
, and
McCarthy
,
C. T.
,
2021
, “
Quality Prediction of Ultrasonically Welded Joints Using a Hybrid Machine Learning Model
,”
J. Manuf. Process.
,
71
(
11
), pp.
571
579
.
27.
Zhou
,
H.
,
Yu
,
K.
,
Luo
,
Q.
,
Luo
,
L.
,
Du
,
W.
, and
Wang
,
S.
,
2022
, “
Design Methods and Strategies for Forward and Inverse Problems of Turbine Blades Based on Machine Learning
,”
J. Thermal Sci.
,
31
(
1
), pp.
82
95
.
28.
Wang
,
Z.
,
Yin
,
Y.
,
Wang
,
Y.
,
Sun
,
T.
, and
Luan
,
Y.
,
2022
, “
Similarity Characteristics of Geometric Scaling Matrix Cooling Channels in Turbine Blade
,”
Appl. Therm. Eng.
,
212
(
7
), p.
118601
.
29.
Kamath
,
R.
,
Venkobarao
,
V.
,
Kopold
,
R.
, and
Subramaniam
,
C.
,
2021
, “
Machine Learning-Based Turbine Vane Position Estimation for Advanced Engine Airpath Control
,”
SAE Int. J. Eng.
,
14
(
6
), pp.
833
852
.
30.
Guiying
,
N.
,
Dunqian
,
C.
, and
Yongquan
,
Z.
,
2019
, “
Improved Differential Evolution Algorithm for Solving 0-1 Programming Problems
,”
J. Syst. Sci. Math. Sci.
,
39
(
1
), pp.
120
132
. .
in Chinese
.
31.
Zhong
,
M.
, and
Long
,
W.
,
2017
, “
Whale Optimization Algorithm Based on Stochastic Adjustment Control Parameter
,”
Sci. Technol. Eng.
,
17
(
12
), pp.
68
73
. .
in Chinese
.
32.
Palacios
,
A.
,
2002
, “
Cycling Chaos in One-Dimensional Coupled Iterated Maps
,”
Int. J. Bifurcat. Chaos
,
12
(
08
), pp.
1859
1868
.
33.
Feng
,
J.
,
Zhang
,
J.
,
Zhu
,
X.
, and
Lian
,
W.
,
2017
, “
A Novel Chaos Optimization Algorithm
,”
Multim. Tools Appl.
,
76
(
16
), pp.
17405
17436
.
34.
Wang
,
G.-Y.
, and
Yuan
,
F.
,
2013
, “
Cascade Chaos and Its Dynamic Characteristics
,”
Acta Phys. Sinica
,
62
(
2
), p.
020506
. .
in Chinese
.
35.
Zhang
,
C.-X.
,
Zhou
,
K.-Q.
,
Ye
,
S.-Q.
, and
Zain
,
A. M.
,
2021
, “
An Improved Cuckoo Search Algorithm Utilizing Nonlinear Inertia Weight and Differential Evolution for Function Optimization Problem
,”
IEEE Access
,
9
(
11
), pp.
161352
161373
.
36.
Ding
,
H.
,
Wu
,
Z.
, and
Zhao
,
L.
,
2020
, “
Whale Optimization Algorithm Based on Nonlinear Convergence Factor and Chaotic Inertial Weight
,”
Concurr. Comput.: Practice Exp.
,
32
(
24
), p.
e5949
.
37.
Long
,
W.
,
Cai
,
S.
,
Jiao
,
J.
,
Tang
,
M.
, and
Wu
,
T.
,
2017
, “
Improved Whale Optimization Algorithm for Large Scale Optimization Problems
,”
Syst. Eng.-Theory Practice
,
37
(
11
), pp.
2983
2994
. .
in Chinese
.
38.
Yue
,
Y.
,
You
,
H.
,
Wang
,
S.
, and
Cao
,
L.
,
2021
, “
Improved Whale Optimization Algorithm and Its Application in Heterogeneous Wireless Sensor Networks
,”
Int. J. Distrib. Sensor Netw.
,
17
(
5
), p.
155014772110181
.
39.
Cao
,
D.
,
Xu
,
Y.
,
Yang
,
Z.
,
Dong
,
H.
, and
Li
,
X.
,
2023
, “
An Enhanced Whale Optimization Algorithm With Improved Dynamic Opposite Learning and Adaptive Inertia Weight Strategy
,”
Complex Intell. Syst.
,
9
(
1
), pp.
767
795
.
40.
Mostafa Bozorgi
,
S.
, and
Yazdani
,
S.
,
2019
, “
IWOA: An Improved Whale Optimization Algorithm for Optimization Problems
,”
J. Comput. Des. Eng.
,
6
(
3
), pp.
243
259
.
41.
Lin
,
M.
,
Wang
,
Z.
,
Chen
,
D.
, and
Zheng
,
W.
,
2022
, “
Particle Swarm-Differential Evolution Algorithm With Multiple Random Mutation
,”
Appl. Soft. Comput.
,
120
(
5
), p.
108640
.
42.
Bai
,
B.
,
Li
,
Z.
,
Hao
,
M.
,
Li
,
Y.
,
Zhang
,
K.
,
Li
,
Z.
, and
Li
,
J.
,
2023
, “
Turbine Blade Endwall Heat Transfer and Film Cooling Performance With Multigap Jets and Film Holes
,”
ASME J. Eng. Gas Turbines Power.
,
145
(
11
), p.
111020
.
43.
Xu
,
Q.
,
Wang
,
P.
, and
Du
,
Q.
,
2022
, “
DDES Investigation of the Cooling Performance and Turbulence Characteristics of the Turbine Vane Trailing Edge Cutback Slot With Fan-Shaped Holes
,”
Int. Commun. Heat Mass Transf.
,
138
(
11
), p.
106405
.
44.
Yeranee
,
K.
,
Rao
,
Y.
,
Yang
,
L.
, and
Li
,
H.
,
2022
, “
Enhanced Thermal Performance of a Pin-Fin Cooling Channel for Gas Turbine Blade by Density-Based Topology Optimization
,”
Int. J. Therm. Sci.
,
181
(
11
), p.
107783
.
45.
Dong
,
Y.-W.
,
Shao
,
P.-F.
,
Guo
,
X.
,
Xu
,
B.
,
Yin
,
C.-P.
, and
Tan
,
Z.-y.
,
2023
, “
Deformation Characterization Method of Typical Double-Walled Turbine Blade Structure During Casting Process
,”
J. Iron Steel Res. Int.
,
30
(
2
), pp.
1
20
.
46.
Zamani
,
H.
,
Nadimi-Shahraki
,
M. H.
, and
Gandomi
,
A. H.
,
2021
, “
QANA: Quantum-Based Avian Navigation Optimizer Algorithm
,”
Eng. Appl. Artif. Intell.
,
104
(
9
), p.
104314
.
47.
Nadimi-Shahraki
,
M. H.
,
2023
, “
An Effective Hybridization of Quantum-Based Avian Navigation and Bonobo Optimizers to Solve Numerical and Mechanical Engineering Problems
,”
J. Bionic Eng.
,
20
(
3
), pp.
1361
1385
.