1. Introduction
As the most basic structural unit of a vehicle, panels have the characteristics of a simple structure, easy loading and unloading, a high load-bearing capacity, etc. They are widely used in the key parts of the wing fuselage skin and air intakes of hypersonic vehicles. However, a dangerous dynamic instability phenomenon, panel flutter [
1,
2], often occurs in actual flight. This is a self-excited oscillation under the coupling of aerodynamic, elastic, and inertial forces, which has strong nonlinear characteristics. Under high dynamic pressure conditions, the structure is prone to Limit Cycle Oscillation (LCO) and dynamic instability, which reduces the fatigue life of the panel and even causes damage. For example, in the late 1950s, NASA’s X-15 test aircraft experienced severe structural vibration in the tail and fairing panels on its first flight, and the conditions of occurrence were much lower than one atmosphere pressure [
3]. Therefore, the study of panel flutter is of great importance for flight safety and aircraft performance assessment.
The systematic study of panel flutter began in the 1950s [
4,
5]. Dowell [
4] analyzed the panel flutter in supersonic airflow using Galerkin’s method, in which the aerodynamic force prediction was based on linear piston theory. The results of the analysis revealed the induced mechanism and evolutionary history of plate aeroelasticity, and gave the critical dynamic pressure of flutter. It was also pointed out that at least sixth-order modes were required to obtain a more accurate solution. These constructive results provided an important reference for the subsequent research on panel flutter problems. Xue et al. [
6] presented a frequency domain method for two-dimensional nonlinear panel flutter with thermal effects obtained from a consistent finite element formulation. The deformation was obtained by von Karman’s nonlinear strain–displacement relation, and the quasi-steady first-order piston theory was employed for aerodynamic loading. The results presented the influence of temperature and dynamic pressure on panel fatigue life. Guo and Mei [
7] investigated the thermo-aeroelastic vibration response of a panel using the structural modal method, which can drastically reduce the number of coupled nonlinear modal equations for the large amplitude nonlinear panel flutter analysis. Using the Ritz method to study the panel aeroelasticity problem, Song and Li [
8] pointed out that the aerodynamic force can significantly change the modal vibration pattern of the structure. Focusing on the flutter problem of composite laminate, Koo et al. [
9] and Singha et al. [
10] established the panel flutter dynamics model based on the finite element method, and discussed the effect of structural and aerodynamic damping. The results indicated that the degree of influence of structural damping on the panel dynamic characteristics is affected by the aerodynamic damping and the fiber layup direction. The effect of structural damping is weak when the aerodynamic damping is large, while the opposite results are obtained when the aerodynamic damping is small. Hamid and Mohammad [
11] presented the nonlinear vibrational frequency analysis of a curved panel under the effects of in-plane compressive and tensile loads, carried out for the first time in the time domain. The first- and third-order piston theories were incorporated for the nonlinear curved panel flutter analysis under the effects of in-plane and thermal loads. Cheng et al. [
12] analyzed the dynamic response of panels in a hypersonic flight environment and pointed out that there are differences in the history of panels to reach chaotic motion at different flight Mach numbers. Xie et al. [
13] considered both aerodynamic nonlinearity and structural nonlinearity in panel flutter analysis, in which third-order piston theory was used for nonlinear aerodynamic loading. The effects of each item of piston theory on panel flutter were studied. The results showed that the geometric term in the quadratic term is the main cause of nonlinearity. Culler et al. [
14] developed a comprehensive aerothermoelastic model for the analysis of panel structures in hypersonic flow, and investigated the impact of fluid–thermal–structural coupling on aerothermoelastic behavior using Galerkin’s method and the third-order piston theory. Chen et al. [
15] analyzed the aeroelastic flutter of a composite panel with functionally graded material. The result illustrated that adding a few amounts of grapheme nanoplatelets can effectually enhance the aeroelastic properties of the plates. The displacement and acceleration feedback control are used to suppress the occurrence of the panel flutter.
During hypersonic flight, the outer surface wall structure of the vehicle is subjected to a rather harsh aerodynamic thermal environment, which seriously affects the performance and flight safety of the vehicle [
16]. The convection cooling system is one of the most commonly used thermal protection strategies [
17,
18,
19,
20]. For instance, Castaldi [
18] presented an effective endothermic fuel platform for regeneratively cooled hypersonic vehicles. The current research on the convection cooling thermal protection system mainly focuses on optimal control [
21], thermal protection efficiency improvement [
16,
22], structure simplification, and mass reduction [
23], but has not yet considered the influence of the coolant force on the stability of the wall plate. Therefore, the dynamics model of panel flutter under the action of coolant is established in this paper based on von Karman’s large deformation theory and the Kirchhoff–Love hypothesis. The kinetic equations are discretized by Galerkin’s method in the space domain and solved by the fourth Runge–Kutta method in the time domain. Based on this, the effects of the head-panel pressure of the coolant, the pressure drop ratio, the coolant injection direction, and the inertial resistance and viscous resistance on panel stability and flight stability were investigated, respectively, and the suppression method of panel flutter was proposed.
This paper is organized as follows.
Section 2 presents the dynamic analytical model of panel flutter under the effect of coolant and verifies the accuracy of the model.
Section 3 analyzes the panel stability under the influence of coolant.
Section 4 analyzes the flight stability of the panel and proposes the suppression method of panel flutter. The conclusions are summarized in
Section 5.
2. Theoretical Analysis
Consider a two-dimensional isotropic simply-supported panel, as shown in
Figure 1, where the panel length is
, the thickness is
, and the transverse vibration displacement is
. The upper surface of the panel is acted on by supersonic airflow, and the Mach number, pressure, and density of the freestream are
, respectively. Below the panel is a convective cooling channel filled with coolant driven by a constant pressure gradient
, ignoring the viscous shear. Based on von Karman’s large deformation thin plate theory and the Kirchhoff–Love hypothesis, the equation of motion for the panel is
The boundary conditions are
where
is transverse vibration displacement,
is material density,
is bending stiffness,
is elastic modulus,
is Poisson ‘s ratio.
is in-plane load caused by assembly or heating, and
is nonlinear membrane load caused by panel deformation.
In order to predict the aerodynamic load under supersonic/hypersonic conditions accurately, the nonlinear modified piston theory [
13] is used.
is the velocity of freestream and is dynamic pressure.
The coolant load on the panel is
The three items on the right side of Equation (5) are coolant pressure, inertial resistance, and viscous resistance. The flow of the coolant in the cooling channel is approximated by the Poiseuille flow. It is called forward flow when the coolant and air flow are in the same direction; otherwise, it is the reverse flow. Take forward flow as an example
is the head-panel pressure. The inertial resistance and viscous resistance are proportional to the second derivative and the first derivative of the deformation to time, respectively.
Through dimensional analysis
where
are the coolant density, the coolant viscosity, the cooling channel thickness, and the thickness of the oscillating boundary layer describing the friction effect, respectively.
The dimensionless variables are defined as
where
denote the in-plane load, air/plate mass ratio, freestream dynamic pressure, head-panel pressure, pressure drop ratio, inertial resistance, and viscous resistance, respectively.
is valued near the atmospheric pressure of the incoming flow, where
is coefficient of head-panel pressure.
cannot be less than 1.0. The non-dimensionalised equation of motion is
Equation (9) is solved using Galerkin’s method. The panel deformation can be approximated by a set of sine functions [
24,
25], which automatically satisfy the simply supported conditions
The dynamic response of the panel is mainly influenced by the low-order mode. Some research work indicated that the first six modes are sufficient to approximate the panel deformation [
4,
14]. Substitute Equation (10) into Equation (9), multiply by
, and then weighted integrate over the panel length. The result is
where
Equations (11) are a set of 2-order nonlinear differential equations in time. Important effecting parameters contain . Setting , the equations become a set of 1-order differential equations, which are solved by the 4 Runge–Kutta method. In this paper, the material properties and geometrical dimensions are . The aircraft flies at an altitude of 30 km. A time step of and initialization value of are adopted. All figures are plots at a typical point .
Figure 2 shows the close agreement between the stability regions given in [
4,
14] and the results of the present model where
. Under the coupling effect of freestream and in-plane load, the dynamic characteristics of the panel are complex, including flat and stable, buckled but dynamically stable, Limit Cycle Oscillation (LCO), anharmonic periodic motion, and chaotic motion. In fact, the regions flat and stable and buckled but dynamically stable belong to panel stabilization, whereas the regions LCO and anharmonic periodic motion and chaotic motion belong to panel flutter. The transition process of the motion state of the panel is called bifurcation. Bifurcation is an important nonlinear phenomenon, which means the global behaviors of nonlinear mathematical systems have changed suddenly.
Figure 3 compares the peak of LCO from the present model to [
4,
6,
14] when in-plane load
, which shows the close agreement with other reference results. More research on panel flutter with coolant effect is carried out in the following sections.
In this section, the dynamic analytical model of panel flutter under the effect of coolant was established. The panel dynamics equations were solved by Galerkin’s method in space domain and the 4th Runge–Kutta method in time domain. The nonlinear modified piston theory was used to predict the unsteady aerodynamic loads. The first six modes were adopted to approximate the panel deformation. Through analysis and derivation, there were eight important effecting parameters: , respectively, which affected the characteristics of the panel flutter. In the end, the accuracy of the solution method was verified.
3. Panel Stability
There are five main factors that contribute to the effect of coolant on panel flutter, including the coefficient of head-panel pressure
,the pressure drop ratio
, the coolant injection direction, inertial resistance
and viscous resistance
. The effect of the pressure drop ratio
on the panel stability, which changes the differential pressure distribution between the upper and lower surfaces of the panel, is discussed first in this section.
Figure 4 shows the time history of the panel response at different pressure drop ratios
when
equals 300.0 and 350.0, respectively. In
Figure 4a, the panel is buckled, but the buckling deformation of the panel decreases as the pressure drop ratio increases. When the pressure drop ratio is infinite, the panel is in a flat and stable state. In
Figure 4b, the motion morph changes with different pressure drop ratios. It is buckled when
equals 1.0 and 2.0 and LCO when
equals 5.0, 10.0, and
. The amplitude of the LCO increases, and the equilibrium position moves towards 0.
Figure 5 shows the curves of the peak of deformation
versus dynamic pressure
under different pressure drop ratios
. It can be seen that the panel bifurcates from the stable state to LCO, and the critical dynamic pressure
at which bifurcation occurs decreases with increasing
. Some cases are shown. The panel bifurcates from the buckling state to LCO at
when
; moreover,
corresponds to
, and
corresponds to
under which the panel bifurcates from flat and stable to LCO.
Figure 6 shows local bifurcation diagrams with several bifurcations of deformation extreme value versus dynamic pressure when
and
. In
Figure 6a, three bifurcations occur: from the buckling state to the anharmonic periodic motion at
; from the anharmonic periodic motion to the buckling state at
; and from the buckling state to the LCO at
. In
Figure 6b, more than three bifurcations occur, and these are not detailed here.
Set
.
Figure 7 shows the time history and phase diagram of the panel deformation when
. The panel is in buckling motion and gathers as a point on the phase plane at
, while the panel performs chaotic motion and the phase diagram is disordered, reflecting strong nonlinearity at
.
Figure 8 shows the results when
, where anharmonic periodic motion occurs. Bifurcation diagrams with several bifurcations of deformation extreme value versus dynamic pressure are shown in
Figure 9. For
, the first bifurcation of the panel occurs at
from the buckling state to the chaotic motion. The second bifurcation occurs at
from chaotic motion to anharmonic periodic motion. The third bifurcation occurs at
from anharmonic periodic motion to LCO. There are other bifurcations, but they are not mentioned here. They all show the complex and strongly nonlinear nature of the panel dynamic response. For
, the bifurcation process is similar but with lower critical dynamic pressure than the result for
. The first three critical dynamic pressures are
, respectively.
Based on the above results, the stability regions of the panel with an in-plane load at different pressure drop ratios are shown in
Figure 10. These curves divide the
plane into two regions: the stable and buckled region and the flutter region. The flutter region has complex motion forms, including LCO, anharmonic periodic motion, chaotic motion, etc. The curves indicate a nonfunctional relationship between
and
. Taking the case with
as an example, three critical dynamic pressures exist on the curve when
, which means three bifurcations and at least four panel response forms, which are buckling motion, flutter (anharmonic periodic motion or chaotic motion), buckling motion, and flutter(LCO) in turn. When
, the panel bifurcates from buckling motion to flutter after which no multiple bifurcations or more complex bifurcation behaviors occur. As the in-plane load increases from 0.0, the critical dynamic pressure of panel flutter decreases as a whole. The stability curves tend to be horizontal with high
and high
. A smaller
results in a larger area of stability and buckling region, which means that the panel is more likely to stabilize.
The coefficient of head-panel pressure
determines the pressure difference on the surface of the panel, which is another important factor affecting the stability of the panel.
Figure 11 plots the peak of deformation versus freestream dynamic pressure
under different coefficients of head-panel pressure
. It can be seen from
Figure 11 that as
increases, the panel bifurcates from buckling motion or stable state to LCO. The buckling deformation and limit cycle amplitude increase gradually and the critical dynamic pressure decreases first and then increases.
Figure 12 shows the stability regions of the panel under an in-plane load for different coefficient of head-panel pressure
. With the increase in
, the area of the stable and buckling region decreases first and then increases. The minimum value occurs at
when
, and at
when
among all cases. When
and
, the upper and lower surfaces of the panel have the same static pressure difference in the opposite direction. However, the area of the stable and buckling region is larger when
than that when
. The same result is concluded when
and
. As the in-plane load increases, the critical dynamic pressure of the panel flutter decreases, and meanwhile, three bifurcations and more complex bifurcation behaviors occur.
Figure 13 shows the relationship between the critical dynamic pressure and the buckling deformation of the panel when the bifurcation occurs. The value
increases from 0.2 to 2.0 along the curves and positive ordinate. The figure mainly reflects two points:
Contradiction between the panel deformation and flutter critical dynamic pressure. It is not difficult to see that the panel deformation is increased with the increase in the critical dynamic pressure. Therefore, not only the flight safety but also the strength safety of the material should be considered in the structural design process.
The safety envelope and flutter envelope of the panel. The
plane is divided into three regions by all curve clusters, labeled I, II, and III in
Figure 13. No flutter occurs regardless of the
and
values in zone I. Flutter always occurs in zone II. Whether flutter occurs or not depends on the values of
. The delimitation curves of zones I and III, and zones I and II form the safety envelope and the flutter envelope of the panel, respectively.
The coolant injection direction is another important factor in the analysis of panel stability.
Figure 14 plots the curves of the peak of deformation versus the freestream dynamic pressure under the condition of the forward injection and reverse injection of coolant. The injection direction mainly affects three aspects of the panel response. First, in the buckling state, the buckling deformation under forward injection is smaller than that under the reverse injection when the dynamic pressure is small, while the magnitude relationship changes when the dynamic pressure is large. There is obviously a transition point. With the increase in
, the values of
corresponding to the transition points are 140, 94, 49.5, and 111.1, respectively. Second, the position at which the bifurcation occurs is also different. The bifurcation occurs first when the coolant is injected reversely when
is equal to 0.6 and 1.0, while opposite results are obtained when
is equal to 1.4 and 1.8. Third, the peak of LCO when the coolant is injected forward is larger than that when the coolant is injected reversely, and the difference between the peaks increases with the increase in
.
Figure 15 shows the stability regions of the panel under the influence of the coolant injection direction. When
, the panel with a small in-plane load is more prone to bifurcation using the reverse injection method. The injection method that leads to earlier bifurcation or multiple bifurcations of the panel with a large in-plane load is correspondingly opposite. The effect of the coolant injection direction on the panel stability when
is opposite to that when
.
Table 1 shows the critical dynamic pressure of the panel without an in-plane load using different injection methods. The results are close when
. The result obtained by the reverse injection is smaller than that of the forward injection when
, and the conclusion is opposite when
.
Considering the influence of inertial and viscous resistances,
Figure 16 shows the convergence process of the buckling response under different inertial and viscous resistances. It indicates that inertial resistance slows down the convergence process, while viscous resistance speeds up the convergence process. However, neither of them change the buckling deformation of the panel.
Figure 17 plots the peak of deformation versus the freestream dynamic pressure under different inertial and viscous resistances, showing that the panel bifurcates from the stable state to LCO. In
Figure 17a, the critical dynamic pressure of the panel decreases, and meanwhile, the peak of LCO gradually increases as inertial resistance increases, but clear changes appear only at small inertial resistances.
Figure 17b shows the influence of viscous resistance, which plays the opposite role to inertial resistance.
Table 2 shows the critical dynamic pressure of the panel under an in-plane load at different inertial resistances. As the in-plane load increases, the critical dynamic pressure first decreases and then tends to be stable. When
, the critical dynamic pressure gradually decreases with increasing inertial resistance. When
, the panel performs multiple bifurcations under the action of large inertial resistances. For instance, the complex multi-bifurcation phenomenon occurs when
under the conditions of
. When
, the critical dynamic pressure is almost constant.
Figure 18 shows the stability regions of the panel under different viscous resistances. As the in-plane load increases, the critical dynamic pressure first decreases and then tends to be stable. As the viscous resistances increase, the critical dynamic pressure increases gradually for small in-plane loads, and remain almost constant for large in-plane loads.
In this section, the stability characteristics of the panel under the action of coolant are presented. The pressure drop ratio with a range of determines the pressure difference distribution on the surfaces of the panel. The critical dynamic pressure under a low pressure drop ratio is larger than that under a high pressure drop ratio. The coefficient of head-panel pressure determines the magnitude of pressure difference. The critical dynamic pressure is large when the coefficient of head-panel pressure is valued far away from 1.0. The effect of the coolant injection direction on panel flutter is affected by the pressure drop ratio and the coefficient of head-panel pressure. When the head-panel pressure is low, the critical dynamic pressure under forward injection is larger than that under reverse injection. Inertial resistance slows down the convergence process of the panel response, while viscous resistance speeds up this process. Inertial resistance reduces the critical dynamic pressure but not obviously, while viscous resistance significantly improves the critical dynamic pressure.
4. Flight Stability and Chatter Suppression
The stability of a panel without an in-plane load at the same height but with different Mach numbers is discussed in this subsection.
Figure 19 shows the variation in the critical dynamic pressure at the panel bifurcation point with a freestream Mach number for different pressure drop ratios. When
is valued less than 20.0, the critical dynamic pressure gradually increases as the Mach number increases. When
, the critical dynamic pressure first decreases and then increases, but the variation is very small. When
, the panel bifurcates several times at a high Mach number. For example, when
, the panel bifurcates from buckling motion to anharmonic periodic motion at
, bifurcates to buckling motion at
, and bifurcates to LCO at
, which greatly reduces the safe dynamic pressure of the flight.
Figure 19 shows that when
, the bifurcation does not occur more than once and the critical dynamic pressure decreases as the pressure drop ratio increases. For
, the largest critical dynamic pressure is 420.1 at
, which is 22.1% higher than the least critical dynamic pressure of 344.0 at
.
Figure 20 shows the variation curves of the critical dynamic pressure of panel flutter with increasing Mach numbers for different coefficients of head-panel pressure. When
, the critical dynamic pressure increases with the increase in the Mach number. When
, the critical dynamic pressure decreases with the increase in the Mach number. When
, the panel bifurcates several times at a high Mach number. When
, the critical dynamic pressure first decreases and then increases as the value of
increases.
Figure 21 and
Figure 22 show the variation curve of the critical dynamic pressure with increasing Mach numbers under inertial resistance and viscous resistance, respectively. It can be seen that all curves in
Figure 21 and
Figure 22 show a gradually decreasing trend as the Mach number increases. As the inertial resistance increases, the critical dynamic pressure decreases gradually, which has an adverse effect on flight safety, but the effect is very small. For example, when
, the critical dynamic pressures are 344.0 at
and 343.06 at
with 0.27% difference. When the inertial resistance takes a value above 200.0, the curve is almost horizontal. As the viscous resistance increases, the critical dynamic pressure when panel bifurcation occurs increases significantly, which leads to an obvious suppression effect on the panel flutter. For
, the critical dynamic pressures are 344.0 at
and 502.8 at
with 50.5% difference.
The panel shows significantly different stability characteristics under the action of coolant. Therefore, this paper proposes the following solutions to suppress fluttering.
(1) Reduce the pressure drop ratio. However, it should be noted that the problem of the material strength limit caused by large flexural deformation of the structure must be considered under a small pressure drop ratio. If the deformation is required to belong and actual flight dynamic pressure is required to reach 380.0, the pressure drop ratio must be less than 2.0.
(2) Increase or decrease the coefficient of head-panel pressure (away from the value of 1.0). The process brings about a large flexural deformation; therefore, the material strength limit problem also needs to be considered.
(3) Adjust the coolant injection direction. For different head-panel pressures, the coolant injection direction has a different effect on the critical dynamic pressure.
(4) Reduce the inertia resistance and improve the viscous resistance. The inertia resistance slows down the convergence process of the panel response, while the viscous resistance speeds up this process. The inertial resistance reduces the critical dynamic pressure, and the viscous resistance increases the critical dynamic pressure. From the theoretical analysis in
Section 2, the inertial resistance is proportional to the coolant concentration, while the viscous resistance is proportional to the coolant viscosity coefficient. Thus, a coolant with a small concentration or a large viscosity coefficient can be selected.
In this section, flight stability of the plate panel is discussed. It can be seen that reducing the pressure drop ratio, and increasing or decreasing the coefficient of head-panel pressure (away from the value of 1.0) can improve the critical dynamic pressure. When the pressure drop ratio is small, the coefficient of head-panel pressure is near 1.0, and the Mach number is high, and the panel is prone to multiple bifurcations occurring. The inertia resistance amplifies the instability of panel flutter, while viscous resistance plays an opposite role. More meaningfully, the flutter suppression strategies are proposed.
5. Discussion and Conclusions
Based on von Karman’s large deformation thin plate theory and the Kirchhoff–Love hypothesis, the dynamics model of panel flutter under the action of coolant is established. The accuracy of the analysis model was verified. On this basis, the effects of the head-panel pressure of coolant, the pressure drop ratio, the coolant injection direction, and the inertial resistance and viscous resistance on panel stability and flight stability were studied, respectively; these are considered to be significant for panel flutter analysis, and the suppression method of panel flutter is proposed. The follow conclusions can be drawn from these results:
(1) The pressure drop ratio determines the pressure difference distribution on the upper and lower surfaces of the panel. A low pressure drop ratio causes multiple bifurcations of the panel with a large in-plane load, and, in general, the critical dynamic pressure when bifurcation occurs under a low pressure drop ratio is larger than that under a high pressure drop ratio. When , the bifurcation occurs no more than one time and the critical dynamic pressure decreases as the pressure drop ratio increases. For , the largest critical dynamic pressure is 420.1 at , which is 22.1% higher than the least critical dynamic pressure valued at 344.0 at . Multiple bifurcations occur at a low Mach number, which greatly reduces the critical dynamic pressure. As the pressure drop ratio decreases, the flexural deformation of the panel when bifurcation occurs increases, so the strength limit of the wall plate material is also considered.
(2) The coefficient of head-panel pressure determines the pressure difference on the upper and lower surfaces of the panel. Values far away from 1.0 can cause a large critical dynamic pressure and avoid multiple bifurcations at high Mach numbers. Moreover, the large panel flexural deformation should be considered.
(3) The coolant injection direction changes the critical dynamic pressure. At low head-panel pressure, the critical dynamic pressure under forward injection is larger than that under reverse injection, while at large head-panel pressure, the result is opposite.
(4) Inertial resistance slows down the convergence process of the panel response, while viscous resistance speeds up this process. Inertial resistance reduces the critical dynamic pressure but not obviously, while viscous resistance significantly improves the critical dynamic pressure. For example, when , the critical dynamic pressures are 344.0 at and 343.06 at with 0.27% difference, while critical dynamic pressures are 344.0 at and 502.8 at with 50.5% difference. The critical dynamic pressure gradually decreases as the Mach number increases.
(5) By analyzing the flutter characteristics of the panel under the action of coolant, the suppression strategy of panel flutter is proposed. Reducing the pressure drop ratio, increasing or decreasing the coefficient of head-panel pressure (away from the value of 1.0), and adjusting the coolant injection direction can play a role. When changing these parameters, the problem of the material strength limit caused by large flexural deformation of the structure must be considered. Reducing inertia resistance and improving viscous resistance can also suppress the flutter. Considering the definitions of inertial resistance and viscous resistance, a coolant with a small concentration or large viscosity coefficient can be selected.