Abstract

With the development of the high performance and high reliability of aeroengine, the blade-tip radial running clearance (BTRRC) of high pressure turbine seriously influences the reliability and performance of aeroengine, wherein the radial deformation control of turbine casing has to be concerned in BTRRC design. To improve BTRRC design, the optimum control-based probabilistic optimization of turbine casing radial deformation was implemented using time-varying least square support vector machine (T-LSSVM) by considering nonlinear material properties and dynamic thermal load. First the T-LSSVM method was proposed and its mathematical model was established. And then the nonlinear dynamic optimal control model of casing radial deformation was constructed with T-LSSVM. Thirdly, through the numerical experiments, the T-LSSVM method is demonstrated to be a promising approach in reducing additional design samples and improving computational efficiency with acceptable computational precision. Through the optimum control-based probabilistic optimization for nonlinear dynamic radial turbine casing deformation, the optimum radial deformation is 7.865 × 10−4 m with acceptable reliability degree 0.995 6, which is reduced by 7.86 × 10−5 m relative to that before optimization. These results validate the effectiveness and feasibility of the proposed T-LSSVM method, which provides a useful insight into casing radial deformation, BTRRC control, and the development of gas turbine with high performance and high reliability.

1. Introduction

Blade-tip radial running clearance (BTRRC) of high pressure turbine (HPT) seriously influences the performance and reliability of gas turbine [1]. Large BTRRC causes the decline of aeroengine performance just like low working efficiency and thrust, high specific fuel consumption, and so forth, while small BTRRC results in the friction fault between blade-tip and casing [2, 3]. Hence, BTRRC has to hold a reasonable value under operating conditions. Hereinto, casing radial deformation is one important factor of directly effecting on BTRRC [4]. The optimum control of casing radial deformation should be done to make BTRRC more reasonable. In advanced propulsive systems, turbine casing bears vast thermal load and pressure difference under working conditions due to high temperature gradient and high rotational velocity [5, 6]. A significant objective of the probabilistic optimal control of casing radial deformation is to seek an optimum radial deformation to offset the radial deformations of turbine disk and blade under satisfying an acceptable reliability degree. Therefore, an effective optimization method-based probabilistic analysis is required for the optimum control of turbine casing radial deformation.

With the growth in the computing power of current computers, finite element (FE) method has become a common and important technique in product development process, which has been applied to the design analysis of aeroengine components like stress analysis, thermal analysis, vibration analysis, and fatigue life estimate [710]. Excessive computational cost and impractical runtime are confronted when FE method is adopted in the probabilistic optimal design of complex structure like turbine casing because of too large computational loads for nonlinear analysis, dynamic analysis, and multiple loop computations. To reduce the runtime, one alternative is to construct a simple surrogate model (also called as response surface method, RSM) to approximate the response of FE solver [5]. The surrogate model expresses the relationship between the objective or constraint functions and design variables by simple equations, which only cost a small number of FE analyses in order to significantly reduce computing time. The typical surrogate models include polynomial response surface method [1116] and support vector machine (SVM) [1721], and so forth. Polynomial response surface model uses least squares regression analysis to fit low-order polynomials to a set of experimental data [15, 16]. Currently, SVM as an intelligent statistical learning method and implicit performance function has been employed to the reliability analysis and optimal design of complex structure instead of the polynomial response surface method, because SVM holds the advantages of small training sample, high computational accuracy, and efficiency. To further improve the computing speed for building SVM model, the least square SVM (LSSVM) was developed and proved to be an effective method for structure reliability design-based steady analysis [18]. However, it is difficult that the LSSVM is directly applied to the nonlinear dynamic optimal analysis of complex structure with time-varying characteristics [22, 23].

To accomplish the optimum control for nonlinear dynamic radial deformation of turbine casing with LSSVM, the time-varying LSSVM (T-LSSVM) method with high efficiency and high precision is presented. And the mathematical model of T-LSSVM is established and verified to be effective and feasible by numerical experimentation. Finally, the optimum control-based probabilistic optimization for nonlinear dynamic radial deformation of turbine casing was performed with the proposed T-LSSVM, subject to constraints on reliability degree and other practical conditions with nonlinear material property and dynamic thermal load.

2. Time-Varying LSSVM (T-LSSVM)

During aeroengine working, thermal stress and pressure difference are of variation. Hence, the probabilistic optimization of casing radial deformation needs a time-varying (or dynamic) optimization method. In this section, the T-LSSVM method is proposed based on LSSVM for dynamic probabilistic optimization.

2.1. Least Square Support Vector Machine

SVM is an important surrogate model based on intelligent statistical learning theory [24] and holds high computational precision and efficiency. LSSVM adopts linear equation in high-dimensional space replacing quadratic programming problem in low-dimensional space to improve computational efficiency by taking the term of squared error as the objective function and regarding an equation as the constraints [17, 18]. The mathematical model of LSSVM is summarized as follows.

From a certain kind of assumed distribution ( and ), the sampling points are generated, and a set of functions which maps a point in the space onto the space is denoted by in which is a set of parameters and is an undetermined parameter vector.

The regression object is to find a function which makes (1) have the lowest expected risk: in which is an error function [21] . Function can be determined by the following method.

For the nonlinear regression, each sampling point is mapped by a nonlinear function onto the high-dimensional space to conduct the linear regression, and then the original space nonlinear regression effect is attained. Thus, the function is rewritten as

Obviously, the problem of solving the regression function can be transformed into obtaining the following optimal solution: Considering the permissible errors, two slack variables and are introduced, and the optimization function is where is a penalty coefficient.

For obtaining the solution of this quadratic program, the Lagrange function is introduced: where , .

In the optimization process, a kernel function is used to replace the inner product in higher dimensional space, where the Lagrange duality problem is expressed by

After getting the optimized solution , , and , the regression estimating function (LSSVM model) is where SV is a set of support vectors for a given sample set; , is the number of sample values; is the number of random variables.

2.2. Time-Varying LSSVM

For time-varying probabilistic optimization problem, the response of each calculation is a stochastic process so that it is difficult to finish the probabilistic analysis of complex structure with response surface method. In the face of this situation, conventional approaches fit a lot of response surface models in this time-varying process to select one reasonable response at one time point as the computational point of probabilistic optimization. However, in all of the loop computations of probabilistic optimization, it is difficult to make this computational point reasonable and feasible, so that the computational accuracy of probabilistic design is low. To resolve this issue using LSSVM, this paper advances the T-LSSVM method for the dynamic probabilistic optimization of turbine casing radial deformation. T-LSSVM is used to calculate a single extreme value rather than a series of dynamic output responses under different input samples within a time domain , which is equivalent to transform a stochastic process of output response into a random variable. In each stochastic analysis, the random variable as the response can be ensured to be the most reasonable and effective, which can pledge the analytical precision. And then the T-LSSVM model is established based on these data for probabilistic optimization. The procedure of probabilistic optimization with T-LSSVM is drawn as follows.(1)Establish the FE model of structure and select reasonable parameters (dynamic load, constraint conditions, time domain, etc.).(2)Extract the samples of input and output variables and define the maximum output responses for each stochastic analysis within time domain.(3)Fit the response surface function of T-LSSVM to these samples.(4)Apply the T-LSSVM function of to accomplish structural dynamic probabilistic optimization.

Obviously, T-LSSVM is promising to reduce computing cost and enhance calculation efficiency. For probabilistic analysis of turbine casing radial deformation, the maximum deformation may be considered as the output response of T-LSSVM for optimum control during the time domain , because the time-varying deformation of turbine casing is secure in as long as the maximum response is safe at the time point.

According to the basic thought of T-LSSVM, the data set with the maximum output response and input random variables is obtained through a number of stochastic analyses in time domain . With the th input samples , the extremum of output response is within time domain . The data set consisting of the maximum output responses is used to fit the extremum response curve : where is a response surface function, called T-LSSVM function (or model) here.

In fact, for dynamic probabilistic design the T-LSSVM function is a key factor because a valid T-LSSVM function is conductive to enhance the efficiency and precision of probabilistic design. If the T-LSSVM function replacing FE model is applied to structural dynamic probabilistic design, this method is called T-LSSVM method, which belongs to the global response surface method. Based on (8), the T-LSSVM function is built as follows: where , is the number of fitted LSSVM functions of in one stochastic analysis process within .

3. Problem Formulation

Casing radial deformation under gas turbine operation is an important component influencing BTRRC variation [13, 25, 26]. In BTRRC control, small casing radial deformation produces a small BTRRC, which degrades the reliability of gas turbine, while large casing radial deformation increases BTRRC and lessens the performance of gas turbine. Clearly, casing radial deformation should be potentially controlled-based probabilistic optimization subject to a reasonable reliability degree.

3.1. FE Model of Turbine Casing

An aeroengine high pressure turbine is studied to build its simplified model in Figure 1(a). Detailed studies of turbine stator system have been conducted [1, 3, 2527]. The casing structure is modeled as a ring-like structure made up of a series of arcs joined together [27]. The lining ring of casing is a sensitive component so that its expansion and contraction lead to the variation of casing radial deformation and blade-tip clearance. Hence, the inner lining ring is only studied and simplified as a cylindrical structure. For the convenient calculation, an axial cross-section of the axisymmetric FE model is extracted as an interesting submodel as shown in Figure 1(b). The left axial degree of freedom of the submodel is restrained to prevent the axial motion. The inner surface of FE model is exposed to turbine inlet hot gas, and the outside surface is exposed to compressor discharge air.

3.2. Parameters Selection

Under operating condition, rotor speed and gas temperature change with time (time-varying), and the flight profile parameter of which is selected for aeroengine load spectrums [13, 27] as shown in Figure 2. Thermal expansion coefficients and thermal conductivities possess obvious nonlinearity as listed in Table 1. Hence, the inner and outer surface temperatures of casing are time-varying and are predetermined by gas temperature and corresponding convection heat transfer coefficients [2830]. The geometric parameters of casing are predetermined in manufacturing processes, which possesses obvious randomness. Because of only analyzing the radial deformation of casing, the thickness is regarded in casing optimization design. For analysis convenience, the maximum values of aforementioned nonlinear and dynamic parameters were selected as random variables by the extremum selection method [31]. Therefore, the distribution characteristics of these parameters were listed in Table 2.

3.3. Probabilistic Analysis Approach

For turbine casing radial deformation analysis, output response and input random variables are denoted by and , respectively. If the predetermined minimum deformation is , the limit state function with (10) is From (11), indicates the system failure while denotes that the system is secure. According to the limit state function, the fact that the maximum casing radial deformation within time domain is greater than the permissible maximum deformation is able to ensure the safety of aeroengine. This is also the reason why the presented T-LSSVM method is applied in the dynamic probabilistic optimization of turbine casing radial deformation.

When , and the random variables , are mutually independent, , , , and meet where and are the functions of mean value and variance function, respectively; and are the mean value and variance function functions of the limit state function ; ; is the number of random variables.

When obeys a normal distribution, its reliability index and reliability degree , respectively, are

3.4. Optimization Model of Turbine Casing Radial Deformation

Constraints for turbine casing optimization consist of design variable constraints and performance constraints. The design variables constraint change ranges. We consider the reliability degree of casing radial deformation as the performance constraints to minimize the maximum value of casing radial deformation. The reliability degree is a function of random parameters. Hence, under the permissible reliability degree and the limited variable ranges, the reliability optimization model is formulated as in which is the th circular computation; is the design variable vector of the th loop; is the time; and are the allowable mean and standard deviation of after the th loop, respectively; and are the mean and standard deviation of output response after the th loop; is the inverse normal distribution function; is the equality constraint of the th loop that the design variables should meet; and are the lower bound and upper bound of random variables; is the predetermined reliability degree, respectively.

In (14), the correlative parameters () are denoted by random design variables and the expression of output response is gained based on these design variables. In line with Monte Carlo (MC) method, the mean value and standard deviation (, is the variation coefficient of design variables) of design variable vector rather than the mean value and standard deviation are inputted into (14) to obtain the mean value model of structural reliability optimal design.

3.5. Basic Idea of Optimization Design for Casing Radial Deformation

The FE method and T-LSSVM are to be adopted in the probabilistic optimal design of casing radial deformation. In this analysis process, the optimization design needs many times loop calculations. The flowchart of structural probabilistic optimization using T-LSSVM is shown in Figure 3.

4. Experimental Study

The experimental study for assessing the proposed method consists of two stages. During the first stage, two numerical test functions with nonlinearity or dynamic are employed in assessing the accuracy and efficiency of T-LSSVM through the comparison of the LSSVM [18]. The second stage is the reliability optimization and verification of casing radial deformation with the T-LSSVM. All calculations are performed on the same computer.

4.1. Numerical Experiment

The first test problem, denoted by P1, is the Haupt function with nonlinearity [5], defined as where , . It has three local optimal solutions. One of them is the global solution with a value of −1.4565 as .

The second test problem, denoted by P2, is a nonlinear dynamic function with semi-infinite and multivariable constraints, defined as where the semi-infinite time-varying constraints () are

It has an optimal solution with a function value of 157.7923.

LSSVM and T-LSSVM are used to conduct the two problems. The comparisons of results are listed in Table 3.

As shown in Table 3, the objective values of LSSVM are −1.319 72 and 149.687 81 for two problems and their corresponding relative errors are 9.391% and 5.136%, respectively. However, the objective values of T-LSSVM are −1.455 91 and 158.183 75, respectively, which are quite close to the theoretical optimal values (the relative errors are only 0.041% and 0.567%). The number of needed samples and elapsed time for two test problems with T-LSSVM are consistent with LSSVM. Obviously, T-LSSVM has higher computational accuracy with keeping computational efficiency in nonlinear dynamic optimization. Hence, T-LSSVM provides a useful insight into the optimum control-based nonlinear dynamic probabilistic optimization of turbine casing radial deformation.

4.2. Optimization Design of Turbine Casing Radial Deformation
4.2.1. Dynamic Deterministic Analysis

Base on the FE model of casing and MC method, the dynamic deterministic analysis of casing radial deformation was fulfilled by the thermosetting coupling analysis. Casing time-varying radial deformation is shown in Figure 4. Obviously, turbine casing deformation reaches the maximum at  s approximately, and the maximum radial deformation is 8.651 × 10−4 m.

In this dynamic analysis, the maximum value of casing radial deformation determines whether the system is safe during the whole time domain [ s], because the limit state function is secure in the whole analysis process if the extremum of casing radial deformation in [ s] is smaller than the predetermined deformation . Therefore, the ERSM is suitable for the probabilistic optimal design of casing radial deformation.

4.2.2. Fit and Validate ERSFs

A small quantity of samples of each random variable was extracted and afterwards inputted into the FE model of turbine casing to gain the corresponding extremum output responses. The obtained many groups of samples are called as fitting samples with input samples and output samples and are used to construct LSSVM and T-LSSVM models. The required sample number and fitting time of LSSVM and T-LSSVM models are drawn in Table 4. Afterwards 200 groups of testing samples were extracted and were used to test LSSVM and T-LSSVM models. The analysis results are shown in Table 4.

As shown in Table 4, in the stage of fitting response surface, T-LSSVM costs 638.2 s, which consumes less runtime than LSSVM (815.3 s) although the same number of samples is needed to fit response surface models. From the test results of two response surface models by 200 times simulations, it is revealed that T-LSSVM saves computing time 1.411 s and improves computational precision 3.611% relative to LSSVM, which demonstrates the high computational efficiency and accuracy of T-LSSVM.

4.2.3. Probabilistic Optimal Design for Turbine Casing Radial Deformation

In the optimization model (14) of casing radial deformation, it is presumed that the reliability degree is 0.99 and the minimum allowable value is 8 × 10−4 m. The reliability optimal design of casing radial deformation is enforced with the proposed T-LSSVM through 10 000 times simulations in each iteration. The probabilistic distributions of maximum output response before and after optimization were shown in Figure 5. From the optimization, the optimal solution of T-LSSVM is = (3 × 10−3, 1 034.3, 5 910, 1.851, and 27.73). The objective values, computational accuracy, optimal time, iterations, and function counter are drawn in Table 5. The distribution of casing radial deformation before and after optimization is shown in Figure 6.

As shown in the Figure 5, the variation of casing radial deformation (or extremum output response) obeys a normal distribution. The mean and standard deviation obtained by T-LSSVM are 7.865 × 10−4 m and 1.834 ×10−5 m.

As shown in Table 5 and Figures 5 and 6, the casing radial deformation after optimization with T-LSSVM reduces by 7.86 × 10−5 m with the reliability degree = 0.995 6 to that before optimization, when the maximum permissible radial deformation = 8 × 10−4 m, which pledges the reliability and work efficiency of blade-tip clearance. The solutions are promising to preferably control blade-tip clearance and develop the high performance and high reliability of gas turbine. The computational efficiency and precision of T-LSSVM are also acceptable for structural nonlinear dynamic probabilistic optimization.

5. Conclusions

To seek the optimum control of casing radial deformation with thermal stress, pressure differential, and nonlinear material properties, a probabilistic optimization method (T-LSSVM) is proposed. Some conclusions are summarized as follows.(1)The mathematical model of T-LSSVM method is established for structural dynamic probabilistic optimization.(2)The optimization model of turbine casing radial deformation is constructed with T-LSSVM.(3)The proposed T-LSSVM is promising to improve computational precision with the acceptable calculation efficiency.(4)The establishment of T-LSSVM is promising to significantly reduce the total computational cost due to a small number of expensive FE analyses.(5)Desirable radial deformation and reasonable reliability are obtained for turbine casing design; meanwhile, T-LSSVM is proved to be feasible and effective in nonlinear dynamic optimization.

Despite current work demonstrates the potential of T-LSSVM, more studies should be performed in our future work to validate T-LSSVM in other complex design optimization problems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to thank the National Natural Science Foundation of China which cosupported the paper (Grant nos. 51175017 and 51275024).