CN110991017A - Flight/propulsion system/jet noise comprehensive real-time model modeling method - Google Patents

Flight/propulsion system/jet noise comprehensive real-time model modeling method Download PDF

Info

Publication number
CN110991017A
CN110991017A CN201911137583.1A CN201911137583A CN110991017A CN 110991017 A CN110991017 A CN 110991017A CN 201911137583 A CN201911137583 A CN 201911137583A CN 110991017 A CN110991017 A CN 110991017A
Authority
CN
China
Prior art keywords
noise
model
engine
jet
flight
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911137583.1A
Other languages
Chinese (zh)
Other versions
CN110991017B (en
Inventor
孙庆彪
黄向华
赵晓春
冉鹏宇
王佳雯
张子豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201911137583.1A priority Critical patent/CN110991017B/en
Publication of CN110991017A publication Critical patent/CN110991017A/en
Application granted granted Critical
Publication of CN110991017B publication Critical patent/CN110991017B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a flight/propulsion system/jet noise comprehensive real-time model modeling method, which comprises the following steps: establishing a jet flow noise real-time prediction model; establishing a turbofan engine component level model with a large bypass ratio; establishing a dynamic model and a kinematic model of the double-engine conveyor; flight/engine/jet noise model synthesis and modification. Aiming at the problems that in the process of researching and designing the aircraft noise, a flight system, a propulsion system and noise calculation are usually carried out separately by a traditional method, the coupling relation among the systems is difficult to be fully considered, the noise condition in the whole flight process cannot be covered, the calculation amount is huge, and the size of jet noise cannot be calculated in real time according to the flight state and the engine state, a jet noise real-time prediction model is established and is combined with a flight model and an engine nonlinear model, so that the real-time simulation of the flight state, the engine performance and the jet noise is realized.

Description

Flight/propulsion system/jet noise comprehensive real-time model modeling method
Technical Field
The invention discloses a flight/propulsion system/jet noise comprehensive real-time model modeling method, and belongs to the technical field of modeling and simulation of aero-engines.
Background
Turbofan engine noise can be classified into jet noise, fan noise, combustion chamber noise, and turbine noise, where jet noise can be viewed approximately as being proportional to the higher order power of the jet velocity, and jet noise increases sharply with increasing jet velocity, which is the dominant noise source for the engine. The existing passive jet noise control technology has limited potential and often has negative influence on the performance of an engine, a real-time simulation model of jet noise needs to be established for actively controlling far-field noise of jet, and the establishment of a flight/propulsion system/jet noise comprehensive model has very important significance in consideration of the influence of the jet noise on the working state and the flight state of the engine.
The simulation technology is an important means for supporting the autonomous research and development of the aero-engine, the research and development efficiency and quality of the aero-engine can be greatly improved, the research and development period is shortened, and the research and development cost is reduced. At present, a plurality of engine complete machine models are known at home and abroad, and most of the models utilize the component characteristics of an engine to establish a nonlinear component-level model of the engine, and Newton-Raphson and other methods are used for solving a nonlinear equation set, or a volume dynamics method is used for replacing a classical iterative algorithm, so that the real-time model simulation of the engine is realized. However, when the noise of the engine is designed, the calculation is usually directly carried out according to an unsteady process changing along with time, for example, various CFD software is used, the calculation amount of the method is huge, only a plurality of sample points are usually calculated, the noise condition in the whole flight process cannot be covered, and the requirement of real-time simulation cannot be met far, and the efficiency problem of the noise simulation of the engine is the main bottleneck of the technologies of improving the noise design efficiency of the engine, developing a noise airworthiness evaluation system, actively controlling noise, optimizing multidisciplinary design and the like at present. Based on the problems, no turbofan engine real-time model with jet flow noise prediction function is published in China.
Disclosure of Invention
The invention provides a flight/propulsion system/jet noise comprehensive real-time model modeling method, aiming at overcoming the defects in the prior art, the invention provides a flight/propulsion system/jet noise comprehensive real-time model modeling method, aiming at the problems that the traditional method has huge calculation amount, cannot cover the noise condition in the whole flight process, cannot meet the real-time simulation requirement and the like in the engine noise design process, a jet noise real-time prediction model is established, and is combined with an engine nonlinear model and a flight model, so that the real-time simulation of the engine jet noise by using engine performance characteristic parameters and flight state parameters is realized.
The technical solution of the invention is as follows: a flight/propulsion system/jet noise comprehensive real-time model, the comprehensive real-time model comprises a jet noise real-time prediction model, a high bypass ratio turbofan engine component level model, and a dynamics and kinematics model of a double-engine transporter; the dynamics and kinematics model of the double-engine transporter calculates to obtain flight height, Mach number and attack angle, and the parameters are used as the input of a component-level model of the turbofan engine with large bypass ratio and a jet noise real-time prediction model; the high bypass ratio turbofan engine component level model calculates to obtain jet flow speed, temperature and air flow of an inner bypass and an outer bypass, the parameters are used as input of a jet flow noise real-time prediction model, and engine thrust is calculated to be used as input of a dynamics and kinematics model of the double-engine transporter; the jet flow noise real-time prediction model ignores the influence of viscous force and heat conduction on the basis of a sound wave equation of a quadrupole sound source sound field caused by turbulence stress tensor, takes the longitudinal length of a jet flow mixing area as a linear function of the diameter of a nozzle, obtains a semi-empirical model of far-field noise, and improves the operation speed. The method comprises the following steps:
1) establishing a jet flow noise real-time prediction model;
2) establishing a turbofan engine component level model with a large bypass ratio;
3) establishing a dynamics and kinematics model of the double-engine conveyor;
4) flight/engine/jet noise model synthesis and modification.
The step 1) comprises the following specific steps:
1.1) partition calculation noise: by analyzing the characteristics of the jet noise of the engine and according to the difference of the types of noise sources, the jet flow field of the turbofan engine which exhausts separately is divided into four areas for analysis and calculation, wherein the four areas are respectively a jet flow core area, a sufficient mixing area, an outer ring shear layer and an inner ring shear layer, and a noise solution equation is simplified aiming at different areas; respectively calculating outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing area noise and caudal vertebra separation noise through a noise solving equation, and obtaining the total sound pressure level of the noise source through the noise superposition of each part;
1.2) correcting the noise model: and correcting the noise model by using the flight speed, the attack angle and the geometric dimension of the engine spray pipe, and calculating the corrected jet flow Mach number according to the following formula:
Figure RE-GDA0002371915560000031
in the formula voAs absolute velocity of jet, MfMach number of flight, αjFor angle of flight, nc,oTaking 0.62 as an adjustment coefficient;
for large bypass ratio turbofan engines, nc,oIn relation to the internal and external culvert flow rate ratio, the following formula is used for correction:
Figure RE-GDA0002371915560000041
and correcting the inner ring shearing and mixing noise, the fully mixed region noise and the caudal vertebra separation noise by adopting the same method.
The step 1.1) of calculating the noise in a partition mode comprises the following specific steps:
1.1.1) calculate the outer ring shear blending noise:
the basic calculation formula of the sound power level of the outer ring shear mixing noise is as follows:
Figure RE-GDA0002371915560000042
Figure RE-GDA0002371915560000043
in the formula, NoIs a speed index, Ve,oThe absolute velocity of the jet flow of the outer duct, CambAt ambient speed of sound, poFor bypass jet density, ρambIs atmospheric density, WoIs a density index, Mc,oMach number is calibrated for jet velocity, as a correction to jet velocity in flight, [ theta ]oIs the pointing angle between the jet noise source and the observation point; when the correction is not considered, directly using the distance and the angle between the nozzle core position and the observation point as the pointing angle at the moment;
at the moment, the basic calculation formula of the total sound pressure level of the shearing and mixing noise of the outer ring of the engine is as follows:
Figure RE-GDA0002371915560000044
in the formula Aj,outFor an ideal jet area, i.e., the jet area when the jet expands fully isentropically to ambient air pressure, CFM56-5B is a convergent nozzle, so the ideal jet area can be found from its nozzle exit area and flow velocity, and the formula is:
Figure RE-GDA0002371915560000045
in the formula, Ath,oFor outer culvert spray pipe throat area, because be the shrink spray pipe, so its area is the nozzle area:
Figure RE-GDA0002371915560000051
1.1.2) calculating the inner ring shear blending noise:
the basic calculation formula of the sound pressure level of the inner ring shearing mixing noise is as follows:
Figure RE-GDA0002371915560000052
the calculation of the ideal jet flow area is the same as step 1.1, wherein:
Figure RE-GDA0002371915560000053
1.1.3) fully-mixed region noise:
the basic calculation formula of the total sound pressure level of the noise in the fully mixed region is as follows:
Figure RE-GDA0002371915560000054
the ideal jet flow area of the sufficient mixing area is calculated according to the mixed airflow of the inner and outer culvert jet flows, and the formula is as follows:
Aj,mix=Aj,in+Aj,out
1.1.4) caudal vertebra separation noise:
the basic calculation formula of the total sound pressure level of the caudal vertebra separation noise is as follows:
Figure RE-GDA0002371915560000055
Figure RE-GDA0002371915560000061
1.1.5) total output noise:
the operation of the sound pressure level is carried out according to a logarithmic rule, the superposition result of n different sound pressure level noise sources is the logarithm of the energy sum, the total jet noise sound pressure level of the engine is the superposition of outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing noise and caudal vertebra separation noise, and the formula is as follows:
Figure RE-GDA0002371915560000063
step 1.2) geometric dimension correction of engine nozzle in noise model correctionNow, the constant term is modified, taking the fully-doped region noise as an example, according to the Lighthill acoustical formula, the intensity of the noise is proportional to the eighth power of the velocity, i.e. N at this timemIs 8, but in practice, as the flow rate increases, N ismThe total noise sound pressure level is mainly determined by the maximum component, in the noise of the sufficient mixing area, the external culvert flow field is the maximum influence factor, the higher the internal culvert flow velocity-velocity ratio is, the lower the total sound pressure level is, the smaller the same area of the culvert jet flow field is, the higher the total sound pressure level is, and the noise of the sufficient mixing area is corrected as follows by combining test data:
Figure RE-GDA0002371915560000062
in the above formula, ve,m=(vomo+vimi)/(mo+mi) Wherein m iso,miRespectively, the culvert flow and the culvert flow.
Step 2) establishing a component-level model of the turbofan engine with a large bypass ratio, treating gas in the engine according to one-dimensional flow, neglecting combustion delay, and establishing the component-level model of the turbofan engine by utilizing the characteristics of the components of the turbofan engine; selecting high-low pressure shaft rotating speed, fan internal pressure ratio, fan external pressure ratio, booster stage pressure ratio and high-low pressure turbine pressure drop ratio, selecting a balance equation of flow of each section and a high-low pressure shaft power matching equation, and performing model solution by adopting a Newton-Raphson method, wherein the concrete steps are as follows:
and (3) modeling the fan by adopting the blade root and the blade tip separately, and selecting 8 initial guess values in the steady-state common working category: low spool speed NlHigh pressure shaft rotational speed NhTip pressure ratio coefficient Z of fanCl,tipRoot pressure ratio coefficient Z of fanCl,corePressure ratio coefficient Z of booster stageCl,ipPressure ratio coefficient Z of high-pressure compressorChHigh pressure turbine pressure ratio coefficient ZThLow pressure turbine pressure ratio coefficient ZTlStructure x ═ x1x2x3x4x5x6x7x8]T= [NlNhZcl,tipZcl,coreZcl,ipZchZThZTl]T
According to the conditions of flow continuity, pressure balance and power balance met in the working process of the engine, the following 8 same working equations are selected:
(1) the flow of the fan outer duct outlet and the inlet of the outer duct tail spray pipe is balanced:
Figure RE-GDA0002371915560000071
(2) the flow of the fan inner duct outlet and the booster stage inlet is balanced:
Figure RE-GDA0002371915560000072
(3) flow balance between the outlet of the booster stage and the inlet of the high-pressure compressor:
Figure RE-GDA0002371915560000073
(4) the flow of the outlet of the combustion chamber is balanced with the inlet of the high-pressure turbine:
Figure RE-GDA0002371915560000074
(5) the flow of the high-pressure turbine outlet and the low-pressure turbine inlet are balanced:
Figure RE-GDA0002371915560000075
(6) the flow of the low-pressure turbine outlet and the tail nozzle of the outer duct is balanced:
Figure RE-GDA0002371915560000076
(7) and (3) power balance of the high-pressure rotor:
Figure RE-GDA0002371915560000081
(8) and (3) low-voltage rotor power balance:
Figure RE-GDA0002371915560000082
performing iterative operation by using a Newton-Raphson method; in the dynamic process, the power difference between the turbine of the engine and the fan and the compressor generates the rotating acceleration, at the moment, the power balance equations (7) and (8) are not established any more, the rotating acceleration is calculated by the following rotor dynamics equation, in the formula, J1,J2The moment of inertia of the low-pressure shaft and the high-pressure shaft of the engine.
Figure RE-GDA0002371915560000083
Figure RE-GDA0002371915560000084
Dynamic iteration is carried out by adopting a one-pass method, so that both real-time performance and model convergence are considered.
And 3) establishing a dynamics and kinematics model of the double-engine conveyor:
the established double-engine transport plane model meets the simulation functions of airplane motion simulation and ground running state during air flight, the airplane during air flight is regarded as a rigid body with six degrees of freedom, and a six-degree-of-freedom motion model is established; the ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; (2) the nose landing gear lifts off the ground until the lift force is larger than the neutral position of the airplane; when the magnitude of the bearing reaction force of the ground to the nose landing gear becomes 0, the stage (1) can be judged to enter the stage (2), and the motion model is calculated by utilizing the fourth-order Runge Kutta algorithm.
And 3) establishing a dynamics and kinematics model of the double-engine conveyer, establishing a six-degree-of-freedom motion model and a ground running model of the double-engine conveyer, and assuming as follows: (1) the airplane is a rigid body, and the mass change in the flying process is ignored; (2) the surface of the earth is approximately a plane, and the influence of the rotation and revolution of the earth is neglected; (3) the gravitational acceleration is constant; (4) the plane of the airplane is symmetrical in left-right quality and geometry according to a plane body coordinate system xoz; (5) neglecting the effect of the spoiler, comprising the following specific steps:
3.1) six-freedom-degree motion model:
the airplane is regarded as a rigid body with six degrees of freedom when flying in the air, and the speed and the angular speed change in all directions of the airplane depend on aerodynamic force and aerodynamic moment acting on the airplane, and thrust and a thrust distance of an engine; the main pneumatic control surfaces are set as follows: the system comprises inner and outer sections of ailerons, flaps, upper and lower sections of rudders, a horizontal stabilizer, inner and outer sections of elevators and spoilers, wherein the rolling moment is obtained through the differential motion of the ailerons, the yawing moment is provided through the rudders, the pitching moment is obtained through the elevators, and each aerodynamic surface has a deflection limit angle and a limit deflection rate; establishing a pneumatic force and pneumatic moment model, and obtaining the pneumatic force and the pneumatic moment according to the deflection angle of the pneumatic control surface and the dimensionless coefficient; the engine thrust is obtained by an engine model; six dimensionless coefficients such as lift coefficient are utilized to calculate the force and moment of the plane along the coordinate system of the plane in all directions in the air, and the six-freedom-degree motion model is as follows:
Figure RE-GDA0002371915560000091
3.2) ground running motion model: the ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; (2) the nose landing gear lifts off to lift force which is larger than the neutral position of the airplane. When the magnitude of the bearing reaction force of the ground to the nose landing gear becomes 0, the stage (1) can be judged to enter the stage (2);
the equation of motion for the first stage of jogging is:
Figure RE-GDA0002371915560000101
after entering the second stage, the two wheels run, the reaction force of the front wheel support is 0, and the motion equation is as follows:
Figure RE-GDA0002371915560000102
and solving the model by using a Longge Kutta method.
And 4) integrating and correcting the flight/engine/jet noise model:
determining a parameter transfer relation among an aircraft, an engine and a jet flow noise model, wherein a dynamics and kinematics model of a double-engine transporter and a jet flow noise real-time prediction model both need a turbofan engine component level model with a large bypass ratio to provide flight altitude and flight speed, and the input of the jet flow noise real-time prediction model also needs a flight attack angle and the distance and angle between a measuring point and an engine spray pipe in consideration of the directivity problem of noise radiation; the jet flow speed, area, temperature and flow of the inner and outer ducts are calculated through a turbofan engine component level model with a large duct ratio and are used as the input of a jet flow noise real-time prediction model. The method comprises the following specific steps:
4.1) calculating according to the flying height to obtain the gas property of the atmospheric environment where the airplane is located, then obtaining the flight state parameters of the airplane, including the flying speed and the attack angle, and simultaneously inputting the geometric parameters of the engine;
4.2) calculating the engine performances such as the thrust oil consumption rate of the engine and the altitude speed of the airplane through a nonlinear model of the turbofan engine, the airplane and a jet flow noise prediction model to preliminarily obtain the noise radiation size of a noise source;
4.3) correcting the noise radiation by utilizing the flying speed, the flying height, the flying attack angle and the like, and correcting the sound source size by utilizing the geometric parameters of the engine; and then, calculating attenuation in the noise propagation process according to the distance angle of the measuring point, and finally obtaining various performances and a noise radiation field of the engine.
The invention has the beneficial effects that:
the model is corrected by utilizing the flight speed, the height, the flight attack angle and the geometrical parameters of the jet pipe, so that the model is more in line with the practical application requirements, the three models are integrated to obtain a flight/propulsion/jet flow noise model meeting the real-time simulation requirements, the performance characteristic parameters, the flight state parameters and the jet flow noise radiation of the engine can be simulated in real time, the model fully considers the mutual relations among the aircraft, the engine and the jet flow noise, can be used for the early work of active jet flow noise control and noise airworthiness certification, and has a great application prospect.
Drawings
FIG. 1 is a sectional view of a jet flow field.
FIG. 2 is a flow of jet noise model calculations.
Fig. 3 is a schematic view of a position parameter.
FIG. 4 is a schematic view of the structural parameters of the nozzle.
FIG. 5 is a turbofan engine/jet noise ensemble model parameter transfer relationship.
FIG. 6 is a view of a turbofan engine/jet noise integration model architecture.
Detailed Description
A flight/propulsion system/jet noise comprehensive real-time model modeling method comprises the following steps:
step 1) establishing a jet flow noise real-time prediction model:
step 1.1) analyzing and calculating the jet flow field of the separately exhausted turbofan engine into a plurality of areas according to different types of noise sources by analyzing the characteristics of the jet flow noise of the engine, wherein the jet flow field of the separately exhausted turbofan engine is mainly divided into four parts, namely a jet flow core area, a sufficient mixing area, an outer ring shear layer and an inner ring shear layer, and the noise solution equation is simplified aiming at different areas. And respectively calculating outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing region noise and tail cone separation noise through a noise solving equation, and then obtaining the total sound pressure level of the noise source through the noise superposition of all parts.
Step 1.2) the noise model is corrected. Since the correction of the jet flow velocity requires the influence of the flight velocity and the flight attack angle to be added to the absolute velocity of the jet flow velocity, the corrected jet flow mach number can be calculated according to the following formula:
Figure RE-GDA0002371915560000121
in the formula voAs absolute velocity of jet, MfMach number of flight, αjFor angle of flight, nc,oFor adjusting the coefficient, refer to the foreign test data, n under the condition of subsonic velocityc,oIt is preferable to take 0.62. Based on analysis of the test data, for a turbofan engine with a large bypass ratio, nc,oMainly related to the ratio of the inner and outer culvert flow rates, the following formula can be approximately used for correction:
Figure RE-GDA0002371915560000122
for the inner ring shearing and mixing noise, the fully mixed region noise and the caudal vertebra separation noise, the correction ideas are similar and are not repeated here.
Step 2) establishing a component level model of the turbofan engine with a large bypass ratio:
the gas is processed in the engine according to one-dimensional flow, combustion delay is ignored, and a turbofan engine component level model is built by utilizing characteristics of the turbofan engine components. The method comprises the steps of selecting high-low pressure shaft rotating speed, fan internal pressure ratio, fan external pressure ratio, booster stage pressure ratio and high-low pressure turbine pressure drop ratio, selecting a balance equation of flow of each section and a high-low pressure shaft power matching equation, and performing model solution by adopting a Newton-Raphson method.
Step 3) establishing a dynamic model and a kinematic model of the double-engine conveyer:
the established model meets the simulation functions of airplane motion simulation and ground running state of air flight. The airplane flying in the air is regarded as a rigid body with six degrees of freedom, and a six-degree-of-freedom motion model is established. The ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; (2) the nose landing gear lifts off to lift force which is larger than the neutral position of the airplane. When the magnitude of the reaction force of the ground to the nose gear becomes 0, it is determined that the stage 1 is shifted to the stage 2. And calculating the motion model by utilizing a fourth-order Runge Kutta algorithm.
Step 4), flight/engine/jet noise model synthesis and correction:
and defining parameter transfer relations among the aircraft, the engine and a jet flow noise model, wherein the engine model and the noise model both need the flight model to provide flight altitude and flight speed, and the input of the noise model also needs a flight attack angle and the distance and the angle of a measuring point from an engine nozzle in consideration of the directivity problem of noise radiation. And calculating the jet flow speed, area, temperature and flow of an internal and external culvert through a turbofan engine model, and taking the jet flow speed, area, temperature and flow as the input of a noise model.
Example 1
In order to facilitate understanding of those skilled in the art, the present invention will be further described with reference to the following examples and drawings, which are not intended to limit the present invention.
In the embodiment, by taking the establishment of a CFM56-5B turbofan engine/jet noise comprehensive real-time model as an example, on the basis of an existing CFM56-5B turbofan engine nonlinear model established through public data, a jet noise model is established and is combined with the engine nonlinear model to obtain a CFM56-5B turbofan engine real-time model with a jet noise prediction function.
Step 1) establishing a jet flow noise prediction model:
through characteristic analysis of engine jet flow noise, a turbofan engine jet flow field for separately exhausting is divided into a plurality of areas for analysis and calculation according to different types of noise sources, the areas are mainly divided into four parts, as shown in fig. 1, a jet flow core area, a sufficient mixing area, an outer ring shear layer and an inner ring shear layer are respectively arranged, and noise solution equations are simplified aiming at different areas.
The jet noise model calculation process is as shown in fig. 2, and the outer ring shear blending noise, the inner ring shear blending noise, the fully blended region noise and the caudal vertebra separation noise are calculated respectively, and then the total sound pressure level of the noise source is obtained through the noise superposition of each part.
Step 1.1) calculating the outer ring shear blending noise:
the basic calculation formula of the sound power level of the outer ring shear mixing noise is as follows:
Figure RE-GDA0002371915560000141
in the formula, NoIs a speed index, Ve,oThe absolute velocity of the jet flow of the outer duct, CambAt ambient speed of sound, poFor bypass jet density, ρambIs atmospheric density, WoIs a density index, Mc,oThe mach number is calibrated for jet velocity as a correction to jet velocity in flight. ThetaoThe pointing angle between the source of the jet noise and the observation point, see fig. 3 for details. When the correction is not considered, the distance and the angle of the nozzle core position from the observation point are directly used as the pointing angle at the moment.
At the moment, the basic calculation formula of the total sound pressure level of the shearing and mixing noise of the outer ring of the engine is as follows:
Figure RE-GDA0002371915560000142
in the formula Aj,outFor an ideal jet area, i.e., the jet area when the jet expands fully isentropically to ambient air pressure, CFM56-5B is a convergent nozzle, so the ideal jet area can be found from its nozzle exit area and flow velocity, and the formula is:
Figure RE-GDA0002371915560000151
in the formula, Ath,oFor outer culvert spray pipe throat area, because be the shrink spray pipe, so its area is the nozzle area:
Figure RE-GDA0002371915560000152
the geometrical parameters used in the formula are shown in fig. 4.
Step 1.2) calculating the inner ring shearing and mixing noise:
similar to the outer ring shear blending noise, the basic calculation formula of the sound pressure level of the inner ring shear blending noise is as follows:
Figure RE-GDA0002371915560000153
the ideal jet area calculation is similar to step 1.1, where:
Figure RE-GDA0002371915560000154
step 1.3) fully mixing the noise in the region:
the basic calculation formula of the total sound pressure level of the noise in the fully mixed region is as follows:
Figure RE-GDA0002371915560000155
Figure RE-GDA0002371915560000161
the ideal jet flow area of the sufficient mixing area is calculated according to the mixed airflow of the inner and outer culvert jet flows, and the formula is as follows:
Aj,mix=Aj,in+Aj,out
step 1.4) separating noise of the caudal vertebra:
the basic calculation formula of the total sound pressure level of the caudal vertebra separation noise is as follows:
Figure RE-GDA0002371915560000162
step 1.5) total output noise:
the operation of sound pressure level is carried out according to logarithmic law (energy law), the superposition result of n different sound pressure level noise sources is the logarithm of the energy sum, the total jet noise sound pressure level of the engine is the superposition of outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing noise and caudal vertebra separation noise, and the formula is as follows:
Figure RE-GDA0002371915560000164
step 1.6) correcting the noise model by using the flight speed, the attack angle and the geometrical size of the engine nozzle:
since the correction of the jet flow velocity requires the influence of the flight velocity and the flight attack angle to be added to the absolute velocity of the jet flow velocity, the corrected jet flow mach number can be calculated according to the following formula:
Figure RE-GDA0002371915560000163
in the formula voAs absolute velocity of jet, MfMach number of flight, αjFor angle of flight, nc,oFor adjusting the coefficient, refer to the foreign test data, n under the condition of subsonic velocityc,oIt is preferable to take 0.62. Based on analysis of the test data, for a turbofan engine with a large bypass ratio, nc,oMainly related to the ratio of the inner and outer culvert flow rates, the following formula can be approximately used for correction:
Figure RE-GDA0002371915560000171
the corrected jet velocity at this time is:
Figure RE-GDA0002371915560000172
for the inner ring shearing and mixing noise, the fully mixed region noise and the caudal vertebra separation noise, the correction ideas are similar and are not repeated here.
For the correction of the engine size, mainly the correction of a constant term is reflected in the correction of a constant term, taking the fully-mixed zone noise as an example, according to the Lighthill acoustic formula, the intensity of the noise is proportional to the eighth power of the speed, namely N at the momentmIs 8, but in practice, as the flow rate increases, N ismWill reduce gradually, the total noise sound pressure level depends on the biggest component mainly, in fully mixing the noise of district, the external culvert flow field is the biggest influence factor, the higher the velocity ratio of internal culvert flow velocity and external culvert flow velocity, the lower the total sound pressure level, the same internal culvert jet flow fieldThe smaller the area, the higher the total sound pressure level. In combination with the experimental data, the noise in the fully-doped region was corrected as follows:
Figure RE-GDA0002371915560000173
in the above formula, ve,m=(vomo+vimi)/(mo+mi) Wherein m iso,miRespectively, the culvert flow and the culvert flow.
Step 2) establishing a component level model of the turbofan engine with a large bypass ratio:
the engine modeling object is a CFM56-5B civil turbofan engine. Usually, the airflow at the outlet of the fan is divided into two paths, and the two paths enter an inner bypass compressor and an outer bypass according to an inner bypass, an outer bypass and a component-level working equation, but for a CFM56-5B turbofan engine with a large bypass ratio, the diameter of a fan blade is large, the characteristic difference between the blade root and the blade tip of the fan is obvious, and the blade root and the blade tip of the fan need to be separately modeled. Selecting 8 initial guess values in a steady-state common working category: low spool speed NlHigh pressure shaft rotational speed NhTip pressure ratio coefficient Z of fanCl,tipRoot pressure ratio coefficient Z of fanCl,corePressure ratio coefficient Z of booster stageCl,ipPressure ratio coefficient Z of high-pressure compressorChHigh pressure turbine pressure ratio coefficient ZThLow pressure turbine pressure ratio coefficient ZTlStructure x ═ x1x2x3x4x5x6x7x8]T= [NlNhZcl,tipZcl,coreZcl,ipZchZThZTl]T. According to the conditions of flow continuity, pressure balance and power balance met in the working process of the engine, the following 8 same working equations are selected:
(9) the flow of the fan outer duct outlet and the inlet of the outer duct tail spray pipe is balanced:
Figure RE-GDA0002371915560000181
(10) the flow of the fan inner duct outlet and the booster stage inlet is balanced:
Figure RE-GDA0002371915560000182
(11) flow balance between the outlet of the booster stage and the inlet of the high-pressure compressor:
Figure RE-GDA0002371915560000183
(12) the flow of the outlet of the combustion chamber is balanced with the inlet of the high-pressure turbine:
Figure RE-GDA0002371915560000184
(13) the flow of the high-pressure turbine outlet and the low-pressure turbine inlet are balanced:
Figure RE-GDA0002371915560000185
(14) the flow of the low-pressure turbine outlet and the tail nozzle of the outer duct is balanced:
Figure RE-GDA0002371915560000186
(15) and (3) power balance of the high-pressure rotor:
Figure RE-GDA0002371915560000191
(16) and (3) low-voltage rotor power balance:
Figure RE-GDA0002371915560000192
iterative operations were performed using the Newton-Raphson method. In the dynamic process, the power difference between the turbine of the engine and the fan and the compressor generates the rotating acceleration, at the moment, the power balance equations (7) and (8) are not established, and the rotating acceleration is moved by the following rotorCalculation by the mechanical equation, wherein J1,J2The moment of inertia of the low-pressure shaft and the high-pressure shaft of the engine.
Figure RE-GDA0002371915560000193
Figure RE-GDA0002371915560000194
Dynamic iteration is carried out by adopting a one-pass method, so that both real-time performance and model convergence are considered.
Step 3) establishing a dynamic model and a kinematic model of the double-engine conveyer:
an airbus A320 airplane is taken as a research object, and a six-degree-of-freedom motion model and a ground running model are established. The A320 airplane adopts two CFM56-5B engines, is a 150-seat-class medium-short range airplane, and has the empty weight of 41t, the maximum takeoff weight of 73.5t, the maximum cruising speed of 0.82Ma and the cruising power of 5000km when being fully loaded. The following assumptions were made: (1) the airplane is a rigid body, and the mass change in the flying process is ignored; (2) the surface of the earth is approximately a plane, and the influence of the rotation and revolution of the earth is neglected; (3) the gravitational acceleration is constant; (4) the plane of the airplane is symmetrical in left-right quality and geometry according to a plane body coordinate system xoz; (5) the effect of the spoiler is ignored.
Step 3.1) six-degree-of-freedom motion model
When flying in the air, an airplane can be regarded as a rigid body with six degrees of freedom, and the speed and the angular speed of the airplane in all directions change depending on aerodynamic force and aerodynamic moment acting on the airplane, and thrust and a thrust moment of an engine. The main pneumatic control surfaces of air passenger a320 are: the wing comprises inner and outer sections of ailerons, a flap, upper and lower sections of rudders, a horizontal stabilizer, inner and outer sections of elevators and a spoiler. Wherein, the rolling moment is obtained through the differential motion of the ailerons, the yawing moment is provided through the rudder, the pitching moment is obtained through the elevator, and each aerodynamic surface has a deflection limit angle and a limit deflection rate. Aerodynamic force and aerodynamic moment models can be established, and the aerodynamic force and the aerodynamic moment can be obtained according to the deflection angle of the aerodynamic control surface and the dimensionless coefficient. The engine thrust is obtained from an engine model. By utilizing six dimensionless coefficients such as the disclosed lift coefficient, the force and moment of the plane along the coordinate system of the plane in all directions borne by the plane in the air can be calculated, and the six-freedom-degree motion model is as follows:
Figure RE-GDA0002371915560000201
step 3.2) ground running motion model
The ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; and (2) lifting the nose landing gear off the ground until the lift force is larger than the neutral position of the airplane, and lifting off the ground. When the magnitude of the reaction force of the ground to the nose gear becomes 0, it is determined that the stage 1 is shifted to the stage 2.
The equation of motion for the first stage of jogging is:
Figure RE-GDA0002371915560000202
after entering the second stage, the two wheels run, the reaction force of the front wheel support is 0, and the motion equation is as follows:
Figure RE-GDA0002371915560000211
and solving the model by using a Longge Kutta method.
Step 4), flight/engine/jet noise model synthesis and correction:
the object to be modeled is a comprehensive model of engine performance and jet noise radiation, and the parameter transfer relation among the flight environment, the engine model and the jet noise model needs to be fully considered. The engine is a strong nonlinear model, and the performance of the engine is related to the altitude and Mach number of the flight besides the characteristics of the engine. The problem of noise radiation in the flight process of an airplane is a complex unsteady process, and three main factors influencing the noise radiation of the airplane change along with the working states of an engine and the airplane, which are mainly represented by the following three aspects: (1) the jet noise source strength depends on the engine power, flight speed and flight attitude; (2) the measurement noise depends on the distance from the airplane to the measurement point, the polar direction angle and the azimuth direction angle; (3) jet velocity and flight velocity produce convection amplification and doppler shift phenomena. In view of the above relationships, the parameter transfer relationships of the integrated model are shown in FIG. 5.
The engine model and the noise model both need to input the flying height and the flying speed, the directivity problem of noise radiation is considered, and the input of the noise model also needs the flying attack angle and the distance and the angle between the measuring point and the engine nozzle. And calculating the jet flow speed, area, temperature and flow of an internal and external culvert through a turbofan engine model, and taking the jet flow speed, area, temperature and flow as the input of a noise model.
The overall structure of the turbofan engine/jet noise model is shown in FIG. 6.
And 4.1) calculating according to the flying height to obtain the gas property of the atmospheric environment where the airplane is located, then obtaining the flight state parameters of the airplane, including the flying speed, the attack angle and the like, and simultaneously inputting the geometric parameters of the engine.
And 4.2) calculating the engine performances such as the thrust oil consumption rate of the engine and the altitude speed of the airplane through a nonlinear model of the turbofan engine, the airplane and a jet flow noise prediction model, and preliminarily obtaining the noise radiation size of the noise source.
And 4.3) correcting the noise radiation by using the flight speed, the flight height, the flight attack angle and the like, and correcting the sound source size by using the geometric parameters of the engine. And then, calculating attenuation in the noise propagation process according to the distance angle of the measuring point, and finally obtaining various performances and a noise radiation field of the engine.
The invention aims at the problems that the traditional method has huge calculation amount, cannot cover the noise condition in the whole flight process, cannot meet the real-time simulation requirement and the like in the engine noise design process, establishes a jet flow noise real-time prediction model, a flight model and a turbofan engine component level model, takes the influence of the flight state on the noise radiation into consideration, corrects the model by utilizing the flight speed, the height, the flight attack angle and the nozzle geometrical parameters to make the model more accord with the practical application requirement, synthesizes the three models to obtain the flight/propulsion/jet flow noise model meeting the real-time simulation requirement, can carry out real-time simulation on the performance characteristic parameters, the flight state parameters and the jet flow noise radiation of the engine, fully considers the interrelation among an aircraft, the engine and the jet flow noise, and can be used for the early work of active jet flow noise control and noise airworthiness certification, has wide application prospect.

Claims (10)

1. A flight/propulsion system/jet noise comprehensive real-time model is characterized in that the comprehensive real-time model comprises a jet noise real-time prediction model, a turbofan engine component level model with a large bypass ratio and a dynamics and kinematics model of a double-engine transporter; the dynamics and kinematics model of the double-engine transporter calculates to obtain flight height, Mach number and attack angle, and the parameters are used as the input of a component-level model of the turbofan engine with large bypass ratio and a jet noise real-time prediction model; the high bypass ratio turbofan engine component level model calculates to obtain jet flow speed, temperature and air flow of an inner bypass and an outer bypass, the parameters are used as input of a jet flow noise real-time prediction model, and engine thrust is calculated to be used as input of a dynamics and kinematics model of the double-engine transporter; the jet flow noise real-time prediction model ignores the influence of viscous force and heat conduction on the basis of a sound wave equation of a quadrupole sound source sound field caused by turbulence stress tensor, takes the longitudinal length of a jet flow mixing area as a linear function of the diameter of a nozzle, obtains a semi-empirical model of far-field noise, and improves the operation speed.
2. A method of modeling a flight/propulsion system/jet noise integrated real time model as claimed in claim 1, comprising the steps of:
1) establishing a jet flow noise real-time prediction model;
2) establishing a turbofan engine component level model with a large bypass ratio;
3) establishing a dynamics and kinematics model of the double-engine conveyor;
4) flight/engine/jet noise model synthesis and modification.
3. A method for modeling a flight/propulsion system/jet noise integrated real-time model according to claim 2, wherein said step 1) comprises the following steps:
1.1) partition calculation noise: by analyzing the characteristics of the jet noise of the engine and according to the difference of the types of noise sources, the jet flow field of the turbofan engine which exhausts separately is divided into four areas for analysis and calculation, wherein the four areas are respectively a jet flow core area, a sufficient mixing area, an outer ring shear layer and an inner ring shear layer, and a noise solution equation is simplified aiming at different areas; respectively calculating outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing area noise and caudal vertebra separation noise through a noise solving equation, and obtaining the total sound pressure level of the noise source through the noise superposition of each part;
1.2) correcting the noise model: and correcting the noise model by using the flight speed, the attack angle and the geometric dimension of the engine spray pipe, and calculating the corrected jet flow Mach number according to the following formula:
Figure RE-FDA0002371915550000021
in the formula voAs absolute velocity of jet, MfMach number of flight, αjFor angle of flight, nc,oTaking 0.62 as an adjustment coefficient;
for large bypass ratio turbofan engines, nc,oIn relation to the internal and external culvert flow rate ratio, the following formula is used for correction:
Figure RE-FDA0002371915550000022
and correcting the inner ring shearing and mixing noise, the fully mixed region noise and the caudal vertebra separation noise by adopting the same method.
4. A modeling method of flight/propulsion system/jet noise integrated real-time model according to claim 3, characterized in that said step 1.1) of calculating noise by partition comprises the following concrete steps:
1.1.1) calculate the outer ring shear blending noise:
the basic calculation formula of the sound power level of the outer ring shear mixing noise is as follows:
Figure RE-FDA0002371915550000023
Figure RE-FDA0002371915550000031
in the formula, NoIs a speed index, Ve,oThe absolute velocity of the jet flow of the outer duct, CambAt ambient speed of sound, poFor bypass jet density, ρambIs atmospheric density, WoIs a density index, Mc.oMach number is calibrated for jet velocity, as a correction to jet velocity in flight, [ theta ]oIs the pointing angle between the jet noise source and the observation point; when the correction is not considered, directly using the distance and the angle between the nozzle core position and the observation point as the pointing angle at the moment;
at the moment, the basic calculation formula of the total sound pressure level of the shearing and mixing noise of the outer ring of the engine is as follows:
Figure RE-FDA0002371915550000032
in the formula Aj,outFor an ideal jet area, i.e., the jet area when the jet expands fully isentropically to ambient air pressure, CFM56-5B is a convergent nozzle, so the ideal jet area can be found from its nozzle exit area and flow velocity, and the formula is:
Figure RE-FDA0002371915550000033
in the formula, Ath,oFor outer culvert spray pipe throat area, because be the shrink spray pipe, so its area is the nozzle area:
Figure RE-FDA0002371915550000034
1.1.2) calculating the inner ring shear blending noise:
the basic calculation formula of the sound pressure level of the inner ring shearing mixing noise is as follows:
Figure RE-FDA0002371915550000035
the calculation of the ideal jet flow area is the same as step 1.1, wherein:
Figure RE-FDA0002371915550000041
1.1.3) fully-mixed region noise:
the basic calculation formula of the total sound pressure level of the noise in the fully mixed region is as follows:
Figure RE-FDA0002371915550000042
the ideal jet flow area of the sufficient mixing area is calculated according to the mixed airflow of the inner and outer culvert jet flows, and the formula is as follows:
Aj,mix=Aj,in+Aj,out
1.1.4) caudal vertebra separation noise:
the basic calculation formula of the total sound pressure level of the caudal vertebra separation noise is as follows:
Figure RE-FDA0002371915550000043
1.1.5) total output noise:
the operation of the sound pressure level is carried out according to a logarithmic rule, the superposition result of n different sound pressure level noise sources is the logarithm of the energy sum, the total jet noise sound pressure level of the engine is the superposition of outer ring shearing mixing noise, inner ring shearing mixing noise, full mixing noise and caudal vertebra separation noise, and the formula is as follows:
Figure RE-FDA0002371915550000044
5. the modeling method of flight/propulsion system/jet noise integrated real-time model as claimed in claim 4, wherein said step 1.2) of engine nozzle geometry modification in noise model modification is embodied in the modification of constant term, such as fully blended region noise, according to Lighthill Acoustic formula, the intensity of noise is proportional to the eighth power of velocity, that is, N at this timemIs 8, but in practice, as the flow rate increases, N ismThe total noise sound pressure level is mainly determined by the maximum component, in the noise of the sufficient mixing area, the external culvert flow field is the maximum influence factor, the higher the internal culvert flow velocity-velocity ratio is, the lower the total sound pressure level is, the smaller the same area of the culvert jet flow field is, the higher the total sound pressure level is, and the noise of the sufficient mixing area is corrected as follows by combining test data:
Figure RE-FDA0002371915550000051
in the above formula, ve,m=(vomo+vimi)/(mo+mi) Wherein m iso,miRespectively, the culvert flow and the culvert flow.
6. The modeling method of a flight/propulsion system/jet noise integrated real-time model according to claim 2, wherein said step 2) builds a turbofan engine component level model with large bypass ratio, gas is processed in the engine according to one-dimensional flow, combustion delay is ignored, and turbofan engine component level model is built by utilizing characteristics of turbofan engine components; selecting high-low pressure shaft rotating speed, fan internal pressure ratio, fan external pressure ratio, booster stage pressure ratio and high-low pressure turbine pressure drop ratio, selecting a balance equation of flow of each section and a high-low pressure shaft power matching equation, and performing model solution by adopting a Newton-Raphson method, wherein the concrete steps are as follows:
to the fan adoptAnd modeling the blade root and the blade tip separately, and selecting 8 initial guess values in a steady-state common working category: low spool speed NlHigh pressure shaft rotational speed NhTip pressure ratio coefficient Z of fanCl,tipRoot pressure ratio coefficient Z of fanCl,corePressure ratio coefficient Z of booster stageCl,ipPressure ratio coefficient Z of high-pressure compressorChHigh pressure turbine pressure ratio coefficient ZThLow pressure turbine pressure ratio coefficient ZTlStructure x ═ x1x2x3x4x5x6x7x8]T=[NlNhZcl,tipZcl,coreZcl,ipZchZThZTl]T
According to the conditions of flow continuity, pressure balance and power balance met in the working process of the engine, the following 8 same working equations are selected:
(1) the flow of the fan outer duct outlet and the inlet of the outer duct tail spray pipe is balanced:
Figure RE-FDA0002371915550000061
(2) the flow of the fan inner duct outlet and the booster stage inlet is balanced:
Figure RE-FDA0002371915550000062
(3) flow balance between the outlet of the booster stage and the inlet of the high-pressure compressor:
Figure RE-FDA0002371915550000063
(4) the flow of the outlet of the combustion chamber is balanced with the inlet of the high-pressure turbine:
Figure RE-FDA0002371915550000064
(5) the flow of the high-pressure turbine outlet and the low-pressure turbine inlet are balanced:
Figure RE-FDA0002371915550000065
(6) the flow of the low-pressure turbine outlet and the tail nozzle of the outer duct is balanced:
Figure RE-FDA0002371915550000066
(7) and (3) power balance of the high-pressure rotor:
Figure RE-FDA0002371915550000067
(8) and (3) low-voltage rotor power balance:
Figure RE-FDA0002371915550000068
performing iterative operation by using a Newton-Raphson method; in the dynamic process, the power difference between the turbine of the engine and the fan and the compressor generates the rotating acceleration, at the moment, the power balance equations (7) and (8) are not established any more, the rotating acceleration is calculated by the following rotor dynamics equation, in the formula, J1,J2The moment of inertia of the low-pressure shaft and the high-pressure shaft of the engine.
Figure RE-FDA0002371915550000071
Figure RE-FDA0002371915550000072
Dynamic iteration is carried out by adopting a one-pass method, so that both real-time performance and model convergence are considered.
7. The modeling method of flight/propulsion system/jet noise integrated real-time model of claim 2, characterized in that said step 3) establishes a dynamics and kinematics model of the dual-engine transport:
the established double-engine transport plane model meets the simulation functions of airplane motion simulation and ground running state during air flight, the airplane during air flight is regarded as a rigid body with six degrees of freedom, and a six-degree-of-freedom motion model is established; the ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; (2) the nose landing gear lifts off the ground until the lift force is larger than the neutral position of the airplane;
when the magnitude of the bearing reaction force of the ground to the nose landing gear becomes 0, the stage (1) can be judged to enter the stage (2), and the motion model is calculated by utilizing the fourth-order Runge Kutta algorithm.
8. The modeling method of flight/propulsion system/jet noise integrated real-time model as claimed in claim 7, wherein said step 3) is to build dynamics and kinematics model of the dual-engine transport plane, build six-degree-of-freedom motion model and ground run model thereof, and make the following assumptions: (1) the airplane is a rigid body, and the mass change in the flying process is ignored; (2) the surface of the earth is approximately a plane, and the influence of the rotation and revolution of the earth is neglected; (3) the gravitational acceleration is constant; (4) the plane of the airplane is symmetrical in left-right quality and geometry according to a plane body coordinate system xoz; (5) neglecting the effect of the spoiler, comprising the following specific steps:
3.1) six-freedom-degree motion model:
the airplane is regarded as a rigid body with six degrees of freedom when flying in the air, and the speed and the angular speed change in all directions of the airplane depend on aerodynamic force and aerodynamic moment acting on the airplane, and thrust and a thrust distance of an engine; the main pneumatic control surfaces are set as follows: the system comprises inner and outer sections of ailerons, flaps, upper and lower sections of rudders, a horizontal stabilizer, inner and outer sections of elevators and spoilers, wherein the rolling moment is obtained through the differential motion of the ailerons, the yawing moment is provided through the rudders, the pitching moment is obtained through the elevators, and each aerodynamic surface has a deflection limit angle and a limit deflection rate; establishing a pneumatic force and pneumatic moment model, and obtaining the pneumatic force and the pneumatic moment according to the deflection angle of the pneumatic control surface and the dimensionless coefficient; the engine thrust is obtained by an engine model; six dimensionless coefficients such as lift coefficient are utilized to calculate the force and moment of the plane along the coordinate system of the plane in all directions in the air, and the six-freedom-degree motion model is as follows:
Figure RE-FDA0002371915550000081
3.2) ground running motion model: the ground running is divided into two stages, namely (1) from a starting position to the liftoff of a nose landing gear; (2) the nose landing gear lifts off to lift force which is larger than the neutral position of the airplane. When the magnitude of the bearing reaction force of the ground to the nose landing gear becomes 0, the stage (1) can be judged to enter the stage (2);
the equation of motion for the first stage of jogging is:
Figure RE-FDA0002371915550000091
after entering the second stage, the two wheels run, the reaction force of the front wheel support is 0, and the motion equation is as follows:
Figure RE-FDA0002371915550000092
and solving the model by using a Longge Kutta method.
9. The modeling method of flight/propulsion system/jet noise integrated real-time model according to claim 2, characterized in that said step 4) flight/engine/jet noise model integration and modification:
determining a parameter transfer relation among an aircraft, an engine and a jet flow noise model, wherein a dynamics and kinematics model of a double-engine transporter and a jet flow noise real-time prediction model both need a turbofan engine component level model with a large bypass ratio to provide flight altitude and flight speed, and the input of the jet flow noise real-time prediction model also needs a flight attack angle and the distance and angle between a measuring point and an engine spray pipe in consideration of the directivity problem of noise radiation; the jet flow speed, area, temperature and flow of the inner and outer ducts are calculated through a turbofan engine component level model with a large duct ratio and are used as the input of a jet flow noise real-time prediction model.
10. The modeling method of flight/propulsion system/jet noise integrated real-time model according to claim 9, wherein said step 4) of flight/engine/jet noise model integration and modification comprises the following steps:
4.1) calculating according to the flying height to obtain the gas property of the atmospheric environment where the airplane is located, then obtaining the flight state parameters of the airplane, including the flying speed and the attack angle, and simultaneously inputting the geometric parameters of the engine;
4.2) calculating the engine performances such as the thrust oil consumption rate of the engine and the altitude speed of the airplane through a nonlinear model of the turbofan engine, the airplane and a jet flow noise prediction model to preliminarily obtain the noise radiation size of a noise source;
4.3) correcting the noise radiation by utilizing the flying speed, the flying height, the flying attack angle and the like, and correcting the sound source size by utilizing the geometric parameters of the engine; and then, calculating attenuation in the noise propagation process according to the distance angle of the measuring point, and finally obtaining various performances and a noise radiation field of the engine.
CN201911137583.1A 2019-11-19 2019-11-19 Modeling method for flight and propulsion system and jet flow noise comprehensive real-time model Active CN110991017B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911137583.1A CN110991017B (en) 2019-11-19 2019-11-19 Modeling method for flight and propulsion system and jet flow noise comprehensive real-time model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911137583.1A CN110991017B (en) 2019-11-19 2019-11-19 Modeling method for flight and propulsion system and jet flow noise comprehensive real-time model

Publications (2)

Publication Number Publication Date
CN110991017A true CN110991017A (en) 2020-04-10
CN110991017B CN110991017B (en) 2022-05-20

Family

ID=70085092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911137583.1A Active CN110991017B (en) 2019-11-19 2019-11-19 Modeling method for flight and propulsion system and jet flow noise comprehensive real-time model

Country Status (1)

Country Link
CN (1) CN110991017B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111833836A (en) * 2020-07-02 2020-10-27 中国航空工业集团公司上海航空测控技术研究所 Civil aircraft cabin or cockpit environmental noise simulator
CN112182765A (en) * 2020-09-30 2021-01-05 中车长春轨道客车股份有限公司 Noise data acquisition method and device
CN112507476A (en) * 2020-11-03 2021-03-16 南京航空航天大学 Integrated modeling method for variable-geometry air inlet and engine
CN112580222A (en) * 2020-12-28 2021-03-30 中国航空工业集团公司沈阳飞机设计研究所 Spray pipe improvement method and system for solving problem of caudal vertebra ablation
CN112924180A (en) * 2021-01-15 2021-06-08 南京航空航天大学 Jet flow noise active control test platform based on engine spray pipe scaling model
CN113378487A (en) * 2021-01-27 2021-09-10 中国民用航空飞行学院 CFD simulation-based aircraft engine jet flow prediction method and system
CN113550840A (en) * 2020-04-26 2021-10-26 中国航发商用航空发动机有限责任公司 Fan shock wave noise prediction and suppression method, device and storage medium
US20220121787A1 (en) * 2020-05-07 2022-04-21 Nanjing University Of Aeronautics And Astronautics Method for component-level non-iterative construction of airborne real-time model of variable-cycle engine
CN114608789A (en) * 2022-04-07 2022-06-10 中国空气动力研究与发展中心低速空气动力研究所 Test method for studying jet flow noise and sound transmission
WO2022126472A1 (en) * 2020-12-17 2022-06-23 大连理工大学 Multiple geometric parameters-adjustable intake/exhaust/engine integrated aviation propulsion system modeling method
CN117382889A (en) * 2023-12-12 2024-01-12 中国航空工业集团公司金城南京机电液压工程研究中心 Power system of high-speed aircraft and transition stage mode switching method thereof

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090265049A1 (en) * 2008-04-22 2009-10-22 Honeywell International, Inc. Aircraft system emissions and noise estimation mechanism
CN106875756A (en) * 2017-04-10 2017-06-20 南京航空航天大学 A kind of continuously declining near aircraft noise optimization method based on ant group algorithm
CN106934074A (en) * 2015-12-29 2017-07-07 中航商用航空发动机有限责任公司 Global optimum's fanjet air intake duct noise abatement method for designing
CN108520134A (en) * 2018-03-30 2018-09-11 天津大学 A kind of Engine Noise weight analysis method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090265049A1 (en) * 2008-04-22 2009-10-22 Honeywell International, Inc. Aircraft system emissions and noise estimation mechanism
CN106934074A (en) * 2015-12-29 2017-07-07 中航商用航空发动机有限责任公司 Global optimum's fanjet air intake duct noise abatement method for designing
CN106875756A (en) * 2017-04-10 2017-06-20 南京航空航天大学 A kind of continuously declining near aircraft noise optimization method based on ant group algorithm
CN108520134A (en) * 2018-03-30 2018-09-11 天津大学 A kind of Engine Noise weight analysis method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HABIBI K 等: "A robust numerical approach for the prediction of turbofan engine noise", 《MEETINGS ON ACOUSTICS. ACOUSTICAL SOCIETY OF AMERICA》 *
刘楠 等: "应用改进粒子群算法的涡轴发动机性能寻优", 《南京航空航天大学学报》 *
刘沛清等: "大型飞机增升装置气动噪声研究进展", 《空气动力学学报》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113550840A (en) * 2020-04-26 2021-10-26 中国航发商用航空发动机有限责任公司 Fan shock wave noise prediction and suppression method, device and storage medium
CN113550840B (en) * 2020-04-26 2023-02-03 中国航发商用航空发动机有限责任公司 Fan shock wave noise prediction and suppression method, device and storage medium
US20220121787A1 (en) * 2020-05-07 2022-04-21 Nanjing University Of Aeronautics And Astronautics Method for component-level non-iterative construction of airborne real-time model of variable-cycle engine
CN111833836A (en) * 2020-07-02 2020-10-27 中国航空工业集团公司上海航空测控技术研究所 Civil aircraft cabin or cockpit environmental noise simulator
CN112182765A (en) * 2020-09-30 2021-01-05 中车长春轨道客车股份有限公司 Noise data acquisition method and device
CN112507476A (en) * 2020-11-03 2021-03-16 南京航空航天大学 Integrated modeling method for variable-geometry air inlet and engine
CN112507476B (en) * 2020-11-03 2024-05-24 南京航空航天大学 Integrated modeling method for variable geometry air inlet channel and engine
WO2022126472A1 (en) * 2020-12-17 2022-06-23 大连理工大学 Multiple geometric parameters-adjustable intake/exhaust/engine integrated aviation propulsion system modeling method
CN112580222A (en) * 2020-12-28 2021-03-30 中国航空工业集团公司沈阳飞机设计研究所 Spray pipe improvement method and system for solving problem of caudal vertebra ablation
CN112580222B (en) * 2020-12-28 2024-02-23 中国航空工业集团公司沈阳飞机设计研究所 Spray pipe improvement method and system for solving tail cone ablation problem
CN112924180A (en) * 2021-01-15 2021-06-08 南京航空航天大学 Jet flow noise active control test platform based on engine spray pipe scaling model
CN112924180B (en) * 2021-01-15 2022-04-05 南京航空航天大学 Jet flow noise active control test platform based on engine spray pipe scaling model
CN113378487B (en) * 2021-01-27 2022-04-12 中国民用航空飞行学院 CFD simulation-based aircraft engine jet flow prediction method and system
CN113378487A (en) * 2021-01-27 2021-09-10 中国民用航空飞行学院 CFD simulation-based aircraft engine jet flow prediction method and system
CN114608789A (en) * 2022-04-07 2022-06-10 中国空气动力研究与发展中心低速空气动力研究所 Test method for studying jet flow noise and sound transmission
CN114608789B (en) * 2022-04-07 2023-03-21 中国空气动力研究与发展中心低速空气动力研究所 Test method for studying jet flow noise and sound transmission
CN117382889A (en) * 2023-12-12 2024-01-12 中国航空工业集团公司金城南京机电液压工程研究中心 Power system of high-speed aircraft and transition stage mode switching method thereof
CN117382889B (en) * 2023-12-12 2024-02-09 中国航空工业集团公司金城南京机电液压工程研究中心 Power system of high-speed aircraft and transition stage mode switching method thereof

Also Published As

Publication number Publication date
CN110991017B (en) 2022-05-20

Similar Documents

Publication Publication Date Title
CN110991017B (en) Modeling method for flight and propulsion system and jet flow noise comprehensive real-time model
Zhang et al. Review on ducted fans for compound rotorcraft
Hooker et al. Over wing nacelle installations for improved energy efficiency
Collier et al. Environmentally responsible aviation-real solutions for environmental challenges facing aviation
CN112613119B (en) Modeling method for multi-geometric parameter adjustable in-line/out integrated aviation propulsion system
Houacine et al. Gauss pseudospectral method for less noise and fuel consumption of aircraft operations
CN111079325B (en) Turbofan engine jet noise real-time calculation and prediction method based on proxy model
Arend et al. Generation after next propulsor research: Robust design for embedded engine systems
Qing et al. Kriging assisted integrated rotor-duct optimization for ducted fan in hover
Mi et al. Numerical simulation on rigid foreign object exclusion in the turboprop engine intake system with a bypass duct
CN113051662B (en) Pneumatic modeling and performance evaluation method for folding wingtip variant aircraft based on CFD and DATCOM
Takahashi et al. Inlet Diffusor Buoyancy-An Overlooked Term in the Thrust Equation
Shea et al. Turbine powered simulator calibration and testing for hybrid wing body powered airframe integration
Kritikos et al. Prediction of certification noise levels generated by contra-rotating open rotor engines
Rostami et al. Multidisciplinary Analysis Program for Light Aircraft (MAPLA)
Baltman et al. A methodology for determining the interdependence of fuel burn and LTO noise of a commercial supersonic transport
Prudden Rotor aerodynamic interaction effects for multirotor unmanned aircraft systems in forward flight
Liang et al. Propeller Installation Effect and Surrogate Based Optimization of a Tail-sitter Drone
Tonti Development of a Flight Dynamics Modelof a Flying Wing Configuration
Misté et al. Assessment of the Engine Installation Performance of a Redesigned Tilt-rotor Intake System
Gibertini et al. Wind-tunnel tests of the ERICA tiltrotor optimised air-intake
de Bruin et al. Aircraft noise simulation environment, suitable for development of noise abatement procedures
Fisher et al. Aeroacoustics research in Europe: The ceas-asc report on 2001 highlights
Jirasek et al. Computational design of S-duct intakes for the NATO AVT-251 multi-disciplinary configuration
Nikolaeva et al. Construction of a digital model and analysis of the characteristics of the air intake of the advanced supersonic passenger aircraft

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant