4.1. Flight Load Calculation of Aircraft Engine
Drawing upon the principles of finite element analysis, this method involves a systematic extraction of incremental trajectory segments within the provided comprehensive flight segment. Subsequently, these diminutive trajectory segments are the focal point for the meticulous computation of aircraft engine flight loads. The procedural framework is delineated as follows:
Step 1 Aircraft Model Specification:
Initial proceedings necessitate the determination of essential design parameters for the aircraft model, encompassing attributes such as the empty weight, the relationship between lift coefficient and angle of attack, zero-lift drag coefficient, lift-induced drag factor, wing area, engine installation angle, and analogous specifications. Employing the least squares method principle, a fitted curve function, predicated on the aircraft’s polar curve and informed by design parameters, is deduced. This mathematical representation facilitates the derivation of lift and drag coefficients across a spectrum of variables, including angle of attack, Mach number, and altitude.
Step 2 Acquisition of Trajectory Data:
Subsequently, requisite data for the computation of aircraft engine flight loads are acquired from measured flight parameters. These encompass fuel weight, airspeed, pitch angle, roll angle, angle of attack, flight path angle, sideslip angle, and other pertinent variables. Retrieval of pertinent flight parameter data is limited to defined temporal intervals.
Step 3 Calculation of Aircraft Drag:
The flight parameter data forms the basis for calculating the aircraft’s lift-to-drag ratio. The computed value is subsequently utilized in Equation (7) to determine the aircraft’s drag within the confines of the current incremental segment.
where
Qi represents the current minute segment of drag,
ρv2/2 is dynamic pressure,
S represents the wing area,
CX is the aircraft’s drag coefficient,
CX0 is the zero-lift drag coefficient,
A is the lift-induced drag factor, Δ
CXh is the drag coefficient’s altitude correction value, and
CYi is the lift coefficient for the current minute segment.
Step 4 Calculate engine thrust
We ascertain the requisite thrust for the aircraft by employing the differential equations governing the aircraft’s center of mass dynamics in the context of the trajectory coordinate system with the contemporaneous flight parameter data. Subsequently, the actual thrust delivered by the engines is determined in accordance with the thrust efficiency conversion ratio.
where
Pi represents the available thrust for the current minor segment of the aircraft,
mi is the mass of the aircraft in the current minor segment;
vi is the flight vector velocity for the current minor segment of the aircraft;
Ri is the engine thrust;
η is the thrust efficiency conversion ratio;
φp is the angle between the engine thrust vector and the aircraft’s axis
Oxb at specific times;
αi is the angle between the projection of airflow axis
Oxa on the aircraft’s symmetry plane and the body axis
Oxb;
βi is the angle between the airflow axis
Oxa and the aircraft’s symmetry plane;
θi is the angle between the trajectory axis
Oxk and the horizontal plane that is positive upwards relative to the terrain.
The ground coordinate system is a fixed system attached to the Earth’s surface. The origin, O, can be set at any location on the ground, typically chosen at the point of aircraft takeoff. The directions of the axes are generally determined using the right-hand rule.
The range for the parameter η typically falls between 90% and 95%. In this study, we opted for a η value of 90%. By selecting a lower efficiency conversion ratio, we intentionally calculate aircraft thrust requirements that are higher than actual values, anticipating a more significant load than reality. This approach ensures that the calculated results lean toward safety.
Step 5 Calculate the lift and side forces
The next step involves the computation of aircraft lift and lateral force by integrating the current flight parameter data with the previously determined aircraft drag values within the current trajectory segment as follows:
where
Yi represents the lift along the
Oya direction;
Zi represents the lateral force along the
Oza direction;
γsi represents the angle between the airflow
Oya-axis and the flight path
Oyk, known as the roll angle; the
Ψai track angle, also referred to as the heading angle, is the angle between the projection of the flight path
Oxk-axis on the horizontal plane and the ground
Oxg-axis, with the convention that rightward is positive.
Step 6 Calculation of load factor
During flight, an aircraft constantly alters its speed, altitude, and heading. The faster these three parameters change, the better the aircraft’s maneuverability. The maneuverability can be characterized by the aircraft’s tangential and normal accelerations, where a larger acceleration corresponds to greater external forces acting on the aircraft. The load factor on the aircraft, defined as the ratio of the sum of all forces on the aircraft, excluding gravity, to the aircraft’s weight, is a critical parameter. The aircraft’s load factor is interconnected with the engine load factor. The load factor is a vector, and its projection on the trajectory coordinate axes is as follows:
here,
nxi,
nyi,
nzi are the engine load factor in the current path coordinate frame, respectively.
Step 7 Check the calculation results
After completing one calculation, check if there are any remaining informational segments. Continue to extract small informational segments and follow steps 2 to 6 to compute engine load information. Ultimately, this process yields engine load information for the entire time interval.
4.2. Engine Flight Load and Thrust Prediction
During the aircraft’s flight load calculations, derivations are conducted using flight parameter information at various instances. When making predictions, only the type of maneuver action is known, and detailed flight parameter data throughout the maneuver action sequence is lacking. Therefore, it is necessary to investigate the parameters governing maneuver actions, establish parameter equations, and provide a parameterized description of the motion.
Step 1 Weight of aircraft
When considering the pattern of aircraft weight variation, the first step involves determining whether the aircraft performs actions such as missile launches or equipment drops during the flight. By excluding such actions, the primary cause of aircraft weight variation becomes fuel consumption. Subsequently, using empirical data, the relationship between the rate of fuel consumption and maneuver actions is calculated. This yields a time-based profile of aircraft weight variation characterized by a nearly steady pattern.
Based on multiple flight profiles, the aircraft’s weight is calculated as 2700 kg, with an initial fuel quantity of 1700 L, and the fuel consumption rates during various time intervals are detailed in
Table 5. Lastly, through the division of missions and task actions, the aircraft’s weight at each moment can be determined.
Step 2 Prediction of pitch angle
When the aircraft performs maneuver actions, there are three types of pitch angle variations: (1) starting from an initial value, maintaining a fixed rate of change until reaching an extreme value, and returning to the initial value at the end of the action; (2) starting from 0°, changing at a fixed rate to reach an extreme value, and maintaining that angle until the end of the action; (3) starting from 0°, changing at a fixed rate to reach an extreme value, and then returning to 0° at a fixed rate (
Figure 17).
Parameter fitting can be achieved by constructing various linear equations. The positional information required during the linear fitting process was acquired during the statistical analysis. This information includes the action’s starting point in time, the time at which the pitch angle reaches its extreme value, the time when the action ends, the initial pitch angle at the start of the action, and the pitch angle at the extreme point of the action. In
Figure 17, the black curve represents the observed pitch angle change of the aircraft during the maneuver, and the red curve (dashed line) represents the linear fitting result.
Step 3 Prediction of roll angle
The variation in the aircraft’s roll angle typically corresponds to turning maneuvers. The aircraft’s roll angle begins at 0° and increases at a specified rate to reach the prescribed roll angle. Depending on the specific maneuver, it maintains the current roll angle for a certain period. It then returns the roll angle to zero at a fixed rate, marking the end of the action, as illustrated in
Figure 18. In most cases, the rate of change in reaching the extreme value of the roll angle and returning from the extreme value to zero is the same. Therefore, similar to the pitch angle, it can be fitted using three linear equations where
x1 corresponds to the start time of the action,
x2 corresponds to the time of the roll angle reaching its extreme value,
x4 corresponds to the end time of the action,
y1 corresponds to the initial roll angle at the start of the action, and
y2 corresponds to the roll angle at the extreme point of the action. In
Figure 18, the black curve represents the observed changes in the aircraft’s roll angle during maneuvers, and the red curve (dashed line) represents the linear fitting result.
Step 4 Altitude and Atmospheric Density
The aircraft’s flight altitude and atmospheric density are crucial parameters in calculating the aircraft’s lift-to-drag ratio. Flight altitude is a function of time based on flight speed and pitch angle. By integrating the previously obtained functions for speed and pitch angle and considering the initial altitude, we can determine the altitude at each moment. The calculation relationship is expressed as Equation (11).
where
represents the flight altitude at time
t1,
H0 is the initial altitude,
v is the flight speed, and
ϑ is the pitch angle.
As shown in
Figure 19, the black curve represents the observed velocity, pitch angle, and altitude profiles during the maneuver action. The red curve (dashed line) represents the results of linear fitting for the velocity and pitch angle profiles and the altitude profile after integration.
An inverse relationship exists between altitude and atmospheric density, as indicated by the relationship in Equation (12). After obtaining information on the altitude variation profile, atmospheric density
ρ at each height
H can be calculated using the following formula:
By solving the parameter histories of each maneuver, the flight trajectory information of the entire predicted profile, and the engine altitude change under the typical profile can be obtained, where the altitude coordinates have been normalized (
Figure 20).
To obtain the engine thrust and load factor of the typical profile, the aircraft performance information and height, pitch angle, roll angle, etc., obtained from steps 1 and 2, are combined with the thrust and load factor calculation method in
Section 4.1. The predicted engine thrust, engine normal load factor, and the tangential load factor of the typical profile can be obtained, as shown in
Figure 21,
Figure 22 and
Figure 23, where each parameter on the vertical axis has been normalized.