CN113239462A - Simulation method for aircraft turbulent environment simulation - Google Patents

Abstract

The invention discloses a simulation method for simulating a turbulent environment of an airplane, which comprises the following steps: entering a turbulence module, and judging whether the aircraft speed resetting mark is true or not; judging whether the airplane flight freezing mark is true or not; judging whether the vacuum speed VTR of the airplane is greater than 15.0; judging whether the airplane reset directory is 0 or not; judging whether the index of the wind profile of the airplane is not 0; judging whether the airplane micro-storm activation mark is true or not; setting a coarse atmosphere requirement flag of the aircraft to true; judging whether the airplane rough atmosphere mark is true or not; judging whether the type of the aircraft turbulence is 0 or not; judging whether the aircraft turbulence intensity is 0 or not; and judging whether the airplane relocation directory is not 0. Has the advantages that: turbulence models such as rough atmosphere, cobblestone turbulence, discrete vertical gust and the like can be provided, the turbulence intensity is 10-level, the gust intensity can be selected within the range of 0-25 m/s, and therefore the atmosphere disturbance environment is simulated by calculating the turbulence disturbance speed.

Description

Simulation method for aircraft turbulent environment simulation

Technical Field

The invention relates to the technical field of simulation methods, in particular to a simulation method for simulating a turbulent environment of an airplane.

Background

Turbulence, also known as turbulence, is a state of flow of a fluid. When the flow rate is very small, the fluids flow in layers and are not mixed with each other, which is called laminar flow or candy slice; gradually increasing the flow speed, starting to generate wave-shaped oscillation on the streamline of the fluid, wherein the frequency and the amplitude of the oscillation increase along with the increase of the flow speed, and the flow condition is called transition flow; when the flow velocity increases to a great extent, the flow lines are no longer clearly distinguishable and there are many small vortices in the flow field, called turbulence, also called turbulence, flow disturbance or turbulence.

With the continuous improvement of living standard, the airplane riding is more and more popular with people when going out, and the airplane is often easily influenced by atmospheric turbulence when flying in the atmosphere, the interference of the atmospheric turbulence can cause the flying performance of the airplane to be reduced and the airplane to be difficult to operate, and even cause the driver to induce oscillation to endanger the flying safety, therefore, the invention provides a simulation method for simulating the turbulent environment of the airplane.

Disclosure of Invention

Aiming at the problems in the related art, the invention provides a simulation method for simulating a turbulent environment of an airplane, so as to overcome the technical problems in the prior related art.

Therefore, the invention adopts the following specific technical scheme:

a simulation method for aircraft turbulent environment simulation comprises the following steps:

s1, entering a turbulence module, judging whether the speed resetting mark of the airplane is true, if so, executing S19, and if not, executing S2;

s2, judging whether the airplane flight freezing mark is true, if so, executing S19, and if not, executing S3;

s3, judging whether the vacuum speed VTR of the airplane is larger than 15.0, if so, executing S4, and if not, executing S19;

s4, judging whether the airplane reset directory is 0, if so, executing S5, and if not, executing S19;

s5, judging whether the aircraft wind profile index is not 0, if so, executing S7, and if not, executing S6;

s6, judging whether the airplane micro-storm activation mark is true, if so, executing S7, otherwise, setting the airplane rough atmosphere requirement mark as false, and executing S8;

s7, setting the requirement rough atmosphere mark of the airplane to be true, and executing S8;

s8, judging whether the aircraft rough atmosphere mark is true, if so, executing S9, otherwise, executing S10;

s9, judging whether the type of the aircraft turbulence is 0, if so, executing S12, and if not, executing S10;

s10, judging whether the turbulence intensity of the airplane is 0, if so, executing S12, and if not, executing S11;

s11, judging whether the airplane relocation directory is not 0, if so, executing S12, and if not, executing S13;

s12, judging the aircraft turbulence intensity factor eta_TUIf the value is 0, executing S14 if the value is 0, and executing S13 if the value is not 0;

s13, judging whether the type of the turbulent flow of the airplane is not 3, if so, calling a vertical gust removal sub-process and executing S15, and if not, calling a vertical gust calculation sub-process and executing S15;

s14, calling a turbulence eliminating sub-process, then calling a vertical gust eliminating sub-process, and exiting the turbulence module;

s15, judging aircraft gust timing T_JETUPWhether the time is less than 30.0 seconds, if yes, executing S17, and if not, executing S16;

s16, judging whether the type of the aircraft turbulence is 3, if so, executing S17, and if not, executing S18;

s17, calling a turbulence removal sub-process, and exiting from the turbulence module;

s18, sequentially calling and calculating other turbulent flow sub-processes, calling and generating a random number sub-process, calling a filter sub-process, calling and calculating a turbulent flow rate sub-process, and exiting from the turbulent flow module;

s19, judging the time T for repositioning the airplane_REPIf not, executing S21 if yes, and executing S20 if no;

s20, judging whether the vacuum velocity VTR of the airplane is less than 15.0, if so, executing S21, and if not, exiting the turbulence module;

s21, sequentially obtaining the clean turbulence velocity UT = VT = WT =0, UT = VT = WT =0 and eta_TU= η_TUL= η_TUD=0, exit turbulence module, where UT is the draden model output X-direction velocity, VT is the draden model output Y-direction velocity, WT is the draden model output Z-direction velocity, η_TUIs the turbulence intensity factor, η_TULIs the upper-period turbulence intensity factor, η_TUDIs the rough atmospheric turbulence intensity factor.

Further, the process of eliminating the vertical gust in S13 includes the following steps:

s1301, entering a process of removing vertical gusts;

s1302, setting a gust distance factor and a gust timing to be 0;

s1303, setting gust intensity factors and vertical gust intensities to be 0;

s1304, setting gust to require that the cobblestone turbulence intensity and the vertical speed of the gust are both 0;

and S1305, exiting the vertical gust cleaning flow.

Further, the process of calculating the vertical gust in S13 includes the following steps:

s1311, entering a vertical gust calculating sub-process;

s1312, judging whether the turbulence intensity signal is not 0, if so, obtaining a turbulence intensity factor of the vertical gust

Wherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtained_JETLOC=0；

S1313, judging aircraft gust timing T_JETUPIf less than 30.0, if yes, T is obtained_JETUP=T_JETUPΔ t, and execute S1314, wherein t represents the timing interval, if not, the distance factor =0.008Mn + distance factor η_COBLOC=0, and S1315 is performed, where Mn denotes the flight mach number, η_COBLOCRepresenting a cobblestone turbulence intensity factor;

s1314, judging T_JETUPWhether the gust turbulence intensity is larger than 40.0 or not is judged, if yes, the cobble turbulence intensity required by gust turbulence is calculated, S1317 is executed, and if not, S1317 is executed;

s1315, judging whether the distance factor is larger than or equal to 1.0, if so, executing S1316, and if not, executing S1317;

s1316, setting the turbulence type number as 0 and the turbulence intensity number as 0 in sequence, requiring the gust mark as false, the mild turbulence mark as false, the moderate turbulence mark as false and the severe turbulence mark as false, and executing S1317;

s1317, interpolation computation gustHeight factor to vertical gust height factor η_JUM=f(h₁) Wherein h is₁Represents a height;

s1318, obtaining vertical gust speed W_JETUP=f(h₂)∙η_JETLOC∙η_JUMWherein h is₂Represents a distance factor;

and S1319, exiting the vertical gust calculating sub-process.

Further, the turbulence removing sub-process in S14 includes the following steps:

s141, entering a turbulence removal sub-process;

s142, calling a first fade-out sub-process (U)_RAL,U_RA) And a second fade-out sub-flow (V)_RAL,V_RA) Wherein, U_RALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine system_RARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine system_RALRepresents the lateral component factor, V, of the turbulent wind speed under the machine system_RARepresenting the lateral component of the turbulent wind speed under the machine system;

s143, calling a third fade-out sub-flow (W)_RAL,W_RA) And a fourth fade-out sub-flow (P)_RAL,P_RA) Wherein W is_RALRepresents the vertical component factor, W, of the turbulent wind speed under the machine system_RARepresents the vertical component, P, of the turbulent wind speed under the machine system_RALPitch angle velocity factor, P, representing turbulent wind velocity generation_RARepresenting the pitch angle velocity resulting from turbulent wind speed;

s144, calling a fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA) Wherein Q is_RALYaw rate factor, Q, representing turbulent wind velocity_RARepresenting yaw rate, R, resulting from turbulent wind speed_RALRoll angular velocity factor, R, representing turbulent wind velocity generation_RARepresenting the roll angular velocity generated by turbulent wind speed;

s145, obtaining eta_TU= η_TUL= η_TUD=0；

S146, clearing all filter outputs F_UI=0、F_VI=0、F_WI=0, I =1,2,3, wherein F_UIRepresenting the longitudinal filter output, F_VIRepresenting the lateral filter output, F_WIRepresenting the vertical filter output;

and S147, exiting the turbulence removing sub-process.

Further, the fade-out sub-process comprises a first fade-out sub-process (U)_RAL,U_RA) And a second decreasing subprocess (V)_RAL,V_RA) And the third fade-out subprocess (W)_RAL,W_RA) And the fourth fade-out sub-process (P)_RAL,P_RA) The fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA）。

Further, the fade-out sub-process comprises the following steps:

entering a fade-out sub-process, and calculating by using a formula RL = RAKFADER RC, wherein RL represents U_RAL、V_RAL、W_RAL、P_RAL、Q_RAL、R_RALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents U_RA、V_RA、W_RA、P_RA、Q_RA、R_RA；

And judging whether the | RL | is smaller than 0.01, if so, obtaining RC =0, and exiting the fade-out sub-process, and if not, obtaining RC = RL, and exiting the fade-out sub-process.

Further, the calculating other turbulent flow sub-flow in S18 includes the following steps:

s1801, entering a sub-process of calculating other turbulence;

s1802, judging whether the rough atmosphere requirement mark is true, if so, calculating a turbulence intensity factor required by wind, and executing S1806, otherwise, executing S1803;

s1803, judging whether the turbulence type is rough atmosphere, if so, obtaining turbulence intensity factor required by the rough atmosphere

And executing S1806, if not, executing S1804;

s1804, judging whether the turbulence type is cobblestone turbulence, if so, obtaining cobblestone atmosphere required turbulence intensity factors

And executing S1805, if not, obtaining eta_TUD=η_COBLOCAnd performs S1806;

s1805, obtaining the intensity of the square wave

And performs S1806;

s1806, judging whether the turbulence intensity signal is 0, if so, executing S1807, and if not, executing 1808;

s1807, setting the rough atmosphere requirement flag to true, and if the rough atmosphere requirement flag is set, obtaining eta_TUD=0, and executing S1808, if not, executing S1808;

s1808 obtaining the increment of the turbulent intensity factor [. eta. ]_TU=η_TUD-η_TULEqual and equal to_TUIf it is greater than 0.01, get [. eta. ]_TUAnd =0.01, and execute S1810, if no, execute S1809;

s1809, Δ η%_TUIf it is less than-0.05, then get Δ η_TU= -0.05, if not, execute S1810;

s1810, obtaining eta in sequence_TU=η_TUL+∆η_TUAnd η_TUL =η_TU；

S1811, calculating a turbulence intensity height factor and a turbulence scale height factor by interpolation, and judging whether the turbulence type is rough atmosphere, if so, executing S1813, and if not, executing S1812;

s1812, judging whether the rough atmosphere requirement mark is true, if so, executing S1813, otherwise, executing S1814;

s1813, sequentially calculating the rough atmospheric turbulence intensity and the turbulence scale to obtain the square wave intensity W of the clean cobblestones_COB=0、A_COB=0, and perform S1821;

s1814, judging whether the turbulence type is cobblestone turbulence, if so, executing S1816, and if not, executing S1815;

s1815, judging whether the turbulence type is vertical gust, if so, executing S1816, and if not, executing S1821;

s1816, obtaining cobblestone timing T_CS=T_CSAnd judging the vacuum speed V of the airplane_TRIf it is greater than 700.0, get η_TU=0.01, and executes S1818, and if not, executes S1817;

s1817, judging vacuum speed V_TRIf it is less than 200.0, if so, the rough atmosphere timing T is obtained_P=0.333333, and execute the next step, if not, get T_P=1/［0.006(V_TR200.0) +3.0 ], and the next step is performed;

s1818, judging T_CSWhether or not greater than T_PIf yes, then T is obtained_CS=0, and executing the next step, if not, executing the next step;

s1819, judging T_CSWhether or not less than T_PIf so, obtaining W_COB=A_COBIf not, obtaining W_COB=-A_COB；

S1820, calculating cobblestone turbulence intensity and turbulence scale;

and S1821, exiting the calculation of other turbulence subprocesses.

Further, the sub-process of generating random numbers in S18 includes the following steps:

s1831, entering into a sub-flow of generating random numbers;

s1832, generating 9 random numbers which are uniformly distributed and are between 0 and 1;

s1833, generating 9 random numbers which are uniformly distributed within-0.5 to + 0.5;

s1834, storing 9 normally distributed random numbers in the last cycle;

s1835, calling external search subroutine and interpolating to solve the normal distribution random number R_U1~R_U3、R_V1~R_V3、R_W1~R_W3Wherein R is_U1~R_U3Denotes a random number, R, normally distributed to U_V1~R_V3Denotes a V-normally distributed random number, R_W1~R_W3Representing that W is normally distributed with random numbers;

s1836, quit the generating random number sub-flow.

Further, the filter subroutine in S18 includes the following steps:

s1841, entering a filter subprocess;

s1842, according to the time constant T₁=2.0L_U/V_TRAnd gain K =2.8264 η_UL_U/V_TRPerforming a calculation, wherein_UIndicates the intensity of the turbulence in the u direction, L_URepresents the scale of the turbulence in the u direction;

s1843, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U1And a random input variable R_f=R_U1Calculating;

s1844, calling EQN200 sub-flow to obtain F_U1L=F_U1、F_U1=N F₀Gain K =1.0, wherein F_UIL(I =1,2, 3) represents the up-cycle vertical filter output, NF₀Represents the filter output;

s1845, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U2And a random input variable R_f=R_U2Calculating;

s1846, calling EQN200 sub-flow to obtain F_U2L=F_U2、F_U2=N F₀Time constant T₁=L_U/V_TRGain, gain

；

S1847, according to EXP = e^-∆t/T1、F₀=F_U3And R_f=R_U3Calculating;

s1848, calling EQN200 sub-flow to obtain F_U3L=F_U3、F_U3=N F₀Time constant T₁=2.0L_V/V_TRGain K = η_U∙2.5133(L_V/V_TR)²Wherein L is_VRepresents the scale of turbulence in the v-direction;

s1849, according to EXP = e^-∆t/T1Filter output value CF₀=F_V1And a random input variable R_f=R_V1Calculating;

s1850, calling EQN200 sub-flow to obtain F_V10L=F_V1、F_V1=N F₀Time constant T₁=L_V/V_TRAnd a gain K =0.15915, wherein F_VIL(I =1,2, 3) represents the upper-period lateral filter output;

s1851, according to EXP = e^-∆t/T1Filter output value CF₀=F_V2Upper period filter output value DF₀=F_V2LRandom input R_f=R_V1And an upper period random input R_fL=R_V2LCalculating;

s1852, calling EQN300 subflow to obtain F_V2L=F_V2、F_V2=N F₀Time constant 1T₁=L_V/V_TRGain K = η_VL_V/2V_TRAnd time constant 2T₂= 3 ∙ T1, wherein η_VRepresenting the turbulence intensity of the turbulence in the v direction;

s1853, according to EXP = e^-∆t/T1Filter output value CF₀=F_V3Upper period filter output value DF₀=F_V3LRandom input R_f=R_V3And an upper period random input R_fL=R_V3LCalculating;

s1854, calling EQN400 sub-flow to obtain F_V3L=F_V3、F_V3=N F₀Time constant T₁=2.0L_W/V_TRAnd gain K = η_W∙2.5133(L_W/V_TR)²Wherein η_WIndicating the intensity of the turbulence in the w direction, L_WRepresents the scale of the turbulence in the w direction;

s1855, according to EXP = e^-∆t/T1Filter output CF₀=F_W1And a random number R_f=R_W1Calculating;

s1856, calling EQN200 sub-flow to obtain F_W1L=F_W1、F_W1=N F₀Time and time constantNumber T₁=2.0L_W/V_TRAnd a gain K =0.15915, wherein F_WIL(I =1,2, 3) represents the upper-cycle vertical filter output;

s1857, according to EXP = e^-∆t/T1Filter output value CF₀=F_W2Upper period filter output value DF₀=F_W2LRandom input R_f=R_W2And an upper period random input R_fL=R_W2LCalculating;

s1858, calling EQN300 subflow to obtain F_W2L=F_W2、F_W2=N F₀Time constant 1T₁=L_W/V_TRGain K = η_WL_W/2V_TRAnd time constant 2T₂=√3∙T1；

S1859, according to EXP = e^-∆t/T1Filter output value CF₀=F_W3Upper period filter output DF₀=F_W3LA random number R_f=R_W3And an upper periodic random number R_fL=R_W3LCalculating;

s1860, call EQN400 sub-flow to get F_W3L=F_W3And F_W3=N F₀；

S1861, quitting the filter subprocess.

Further, the turbulence rate calculating sub-process in S18 includes the following steps:

s1871, entering a sub-process for calculating the turbulent flow rate;

s1872, judging whether the turbulence type is vertical gust, and if so, obtaining a turbulence factor eta_THL=1.0, and execute S1878, if not, execute S1873;

s1873 and judging the flying height h_RWhether or not less than the turbulent low-height limit h_L-0.5, if yes, perform S1875, if no, perform S1874;

s1874 and judging h_RWhether or not it is greater than h_U+0.5, if yes, execute S1875, if no, execute S1876;

s1875 obtaining a turbulence factor eta_THL=0, and perform S1878;

s1876 and judging h_RWhether or not h is less than or equal to_LIf yes, get η_THL=1.0-(h_L-h_R) ∙ 2.0.0, and executing S1878, if not, executing S1877;

s1877 and judging h_RWhether or not the height of the turbulent flow is greater than or equal to the limit h_UIf yes, get η_THL=1.0+(h_U-h_R) ∙ 2.0.0, and executing S1878, if not, obtaining eta_THL=1.0, and perform S1878;

s1878, calculating U in sequence_RAL=F_U1∙F_U2+F_U3、U_RAL=F_V1∙F_V2+F_V3、w_RAL= F_W1∙F_W2+F_W3∙W_JETUP+W_COB、P_RAL=0.025l(W_RAL-W_JETUP)/L_W、Q_RAL=0.025(W_RAL-W_JETUP)81.5/L_W、R_RAL=0.025∙81.5V_RAL/L_W、U_T=U_RAL∙η_THL、V_T=V_RAL∙η_THL、W_T=W_RAL∙η_THL、P_T=P_RAL∙η_THL、Q_T=Q_RAL∙η_THL、R_T=R_RAL∙η_THL、

、

、

、U_TOL=U_T、V_TOL=V_TAnd W_TOL=W_TWherein P is_TRepresenting the DRYDEN model output pitch angular velocity, Q_TRepresenting DRYDEN model output yaw rate, R_TRepresenting the DRYDEN model output roll angular velocity, U_TLRepresents the X-direction speed, V, of the DRYDEN model output in the upper period_TLRepresents the Y-direction speed, W of the DRYDEN model output in the upper period_TLTo representThe upper period DRYDEN model outputs Z-direction speed,

represents;

represents the acceleration of the output Y direction of the DRYDEN model,

represents the Z-direction acceleration of the DRYDEN model output; (ii) a

S1879 and exiting the turbulence rate calculation sub-process.

The invention has the beneficial effects that: by using the invention, turbulence models such as rough atmosphere, cobblestone turbulence, discrete vertical gust and the like can be provided, the turbulence intensity is 10-level, and the gust intensity can be selected within the range of 0-25 m/s, so that the atmosphere disturbance environment can be simulated by calculating the turbulence disturbance speed. In addition, the invention can be divided into vertical gust turbulence, rough atmosphere and cobblestone turbulence according to the type of the turbulence, and can be divided into nine types of simulated turbulence of light turbulence, moderate turbulence and severe turbulence according to the intensity of the turbulence.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.

FIG. 1 is a schematic flow chart of a simulation method for aircraft turbulent environment simulation to eliminate vertical gusts according to an embodiment of the present invention;

FIG. 2 is a schematic view of a sub-flow of calculating a vertical gust in a simulation method for simulating an aircraft turbulent environment according to an embodiment of the present invention;

FIG. 3 is a schematic diagram of a turbulence removal sub-flow of a simulation method for aircraft turbulent environment simulation according to an embodiment of the present invention;

FIG. 4 is a schematic diagram of a fade-out sub-flow in a simulation method for aircraft turbulent environment simulation according to an embodiment of the present invention;

FIG. 5 is a schematic diagram of a sub-flow of generating random numbers in a simulation method for simulating an aircraft turbulent environment according to an embodiment of the invention;

FIG. 6 is a schematic view of a sub-flow EQN200 of a simulation method for simulating a turbulent environment of an aircraft according to an embodiment of the present invention;

FIG. 7 is a schematic view illustrating a sub-flow EQN300 of a simulation method for simulating a turbulent environment of an aircraft according to an embodiment of the present invention;

FIG. 8 is a schematic view of a sub-flow EQN400 in a simulation method for simulating an aircraft turbulent environment according to an embodiment of the present invention.

Detailed Description

For further explanation of the various embodiments, the drawings which form a part of the disclosure and which are incorporated in and constitute a part of this specification, illustrate embodiments and, together with the description, serve to explain the principles of operation of the embodiments, and to enable others of ordinary skill in the art to understand the various embodiments and advantages of the invention, and, by reference to these figures, reference is made to the accompanying drawings, which are not to scale and wherein like reference numerals generally refer to like elements.

According to an embodiment of the invention, a simulation method for aircraft turbulent environment simulation is provided.

The present invention will be further described with reference to the accompanying drawings and the detailed description, and as shown in fig. 1 to 8, a simulation method for aircraft turbulent environment simulation according to an embodiment of the present invention includes the following steps:

s1, entering a turbulence module, judging whether the speed resetting mark of the airplane is true, if so, executing S19, and if not, executing S2;

s2, judging whether the airplane flight freezing mark is true, if so, executing S19, and if not, executing S3;

s3, judging whether the vacuum speed VTR of the airplane is larger than 15.0, if so, executing S4, and if not, executing S19;

s4, judging whether the airplane reset directory is 0, if so, executing S5, and if not, executing S19;

s5, judging whether the aircraft wind profile index is not 0, if so, executing S7, and if not, executing S6;

s6, judging whether the airplane micro-storm activation mark is true, if so, executing S7, otherwise, setting the airplane rough atmosphere requirement mark as false, and executing S8;

s7, setting the requirement rough atmosphere mark of the airplane to be true, and executing S8;

s8, judging whether the aircraft rough atmosphere mark is true, if so, executing S9, otherwise, executing S10;

s9, judging whether the type of the aircraft turbulence is 0, if so, executing S12, and if not, executing S10;

s10, judging whether the turbulence intensity of the airplane is 0, if so, executing S12, and if not, executing S11;

s11, judging whether the airplane relocation directory is not 0, if so, executing S12, and if not, executing S13;

s12, judging the aircraft turbulence intensity factor eta_TUIf the value is 0, executing S14 if the value is 0, and executing S13 if the value is not 0;

s13, judging whether the type of the turbulent flow of the airplane is not 3, if so, calling a vertical gust removal sub-process and executing S15, and if not, calling a vertical gust calculation sub-process and executing S15;

the process for eliminating the vertical gust comprises the following steps:

s1301, entering a process of removing vertical gusts;

s1302, setting a gust distance factor and a gust timing to be 0;

s1303, setting gust intensity factors and vertical gust intensities to be 0;

s1304, setting gust to require that the cobblestone turbulence intensity and the vertical speed of the gust are both 0;

and S1305, exiting the vertical gust cleaning flow.

The process for calculating the vertical gust comprises the following steps:

s1311, entering a vertical gust calculating sub-process;

s1312, judging whether the turbulence intensity signal is not 0, if so, obtaining a turbulence intensity factor of the vertical gust

Wherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtained_JETLOC=0；

S1313, judging aircraft gust timing T_JETUPIf less than 30.0, if yes, T is obtained_JETUP=T_JETUPΔ t, and execute S1314, wherein t represents the timing interval, if not, the distance factor =0.008Mn + distance factor η_COBLOC=0, and S1315 is performed, where Mn denotes the flight mach number, η_COBLOCRepresenting a cobblestone turbulence intensity factor;

s1314, judging T_JETUPWhether the gust turbulence intensity is larger than 40.0 or not is judged, if yes, the cobble turbulence intensity required by gust turbulence is calculated, S1317 is executed, and if not, S1317 is executed;

s1315, judging whether the distance factor is larger than or equal to 1.0, if so, executing S1316, and if not, executing S1317;

s1316, setting the turbulence type number as 0 and the turbulence intensity number as 0 in sequence, requiring the gust mark as false, the mild turbulence mark as false, the moderate turbulence mark as false and the severe turbulence mark as false, and executing S1317;

s1317, calculating gust height factor to vertical gust height factor eta by interpolation_JUM=f(h₁) Wherein h is₁Represents a height;

s1318, obtaining vertical gust speed W_JETUP=f(h₂)∙η_JETLOC∙η_JUMWherein h is₂Represents a distance factor;

and S1319, exiting the vertical gust calculating sub-process.

S14, calling a turbulence eliminating sub-process, then calling a vertical gust eliminating sub-process, and exiting the turbulence module;

wherein, the turbulence eliminating sub-process comprises the following steps:

s141, entering a turbulence removal sub-process;

s142, calling a first fade-out sub-process (U)_RAL,U_RA) And a second fade-out sub-flow (V)_RAL,V_RA) Wherein, U_RALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine system_RARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine system_RALRepresents the lateral component factor, V, of the turbulent wind speed under the machine system_RARepresenting the lateral component of the turbulent wind speed under the machine system;

s143, calling a third fade-out sub-flow (W)_RAL,W_RA) And a fourth fade-out sub-flow (P)_RAL,P_RA) Wherein W is_RALRepresents the vertical component factor, W, of the turbulent wind speed under the machine system_RARepresents the vertical component, P, of the turbulent wind speed under the machine system_RALPitch angle velocity factor, P, representing turbulent wind velocity generation_RARepresenting the pitch angle velocity resulting from turbulent wind speed;

s144, calling a fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA) Wherein Q is_RALYaw rate factor, Q, representing turbulent wind velocity_RARepresenting yaw rate, R, resulting from turbulent wind speed_RALRoll angular velocity factor, R, representing turbulent wind velocity generation_RARepresenting the roll angular velocity generated by turbulent wind speed;

s145, obtaining eta_TU= η_TUL= η_TUD=0；

S146, clearing all filter outputs F_UI=0、F_VI=0、F_WI=0, I =1,2,3, wherein F_UIRepresenting the longitudinal filter output, F_VIRepresenting the lateral filter output, F_WIRepresenting the vertical filter output;

and S147, exiting the turbulence removing sub-process.

Specifically, the fade-out sub-process includes a first fade-out sub-process (U)_RAL,U_RA) And a second decreasing subprocess (V)_RAL,V_RA) And the third fade-out subprocess (W)_RAL,W_RA) And the fourth fade-out sub-process (P)_RAL,P_RA) The fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA）。

The fade-out sub-process comprises the following steps:

entering a fade-out sub-process, and calculating by using a formula RL = RAKFADER RC, wherein RL represents U_RAL、V_RAL、W_RAL、P_RAL、Q_RAL、R_RALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents U_RA、V_RA、W_RA、P_RA、Q_RA、R_RA；

And judging whether the | RL | is smaller than 0.01, if so, obtaining RC =0, and exiting the fade-out sub-process, and if not, obtaining RC = RL, and exiting the fade-out sub-process.

S15, judging aircraft gust timing T_JETUPWhether the time is less than 30.0 seconds, if yes, executing S17, and if not, executing S16;

s16, judging whether the type of the aircraft turbulence is 3, if so, executing S17, and if not, executing S18;

s17, calling a turbulence removal sub-process, and exiting from the turbulence module;

s18, sequentially calling and calculating other turbulent flow sub-processes, calling and generating a random number sub-process, calling a filter sub-process, calling and calculating a turbulent flow rate sub-process, and exiting from the turbulent flow module;

wherein, the calculation of other turbulent flow sub-processes comprises the following steps:

s1801, entering a sub-process of calculating other turbulence;

s1802, judging whether the rough atmosphere requirement mark is true, if so, calculating a turbulence intensity factor required by wind, and executing S1806, otherwise, executing S1803;

s1803, judging whether the turbulence type is rough atmosphere, if so, obtaining turbulence intensity factor required by the rough atmosphere

And executing S1806, if not, executing S1804;

s1804, judging whether the turbulence type is presentObtaining the cobble atmospheric required turbulence intensity factor if the cobble atmospheric required turbulence intensity factor is cobble turbulence

And executing S1805, if not, obtaining eta_TUD=η_COBLOCAnd performs S1806;

s1805, obtaining the intensity of the square wave

And performs S1806;

s1806, judging whether the turbulence intensity signal is 0, if so, executing S1807, and if not, executing 1808;

s1807, setting the rough atmosphere requirement flag to true, and if the rough atmosphere requirement flag is set, obtaining eta_TUD=0, and executing S1808, if not, executing S1808;

s1808 obtaining the increment of the turbulent intensity factor [. eta. ]_TU=η_TUD-η_TULEqual and equal to_TUIf it is greater than 0.01, get [. eta. ]_TUAnd =0.01, and execute S1810, if no, execute S1809;

s1809, Δ η%_TUIf it is less than-0.05, then get Δ η_TU= -0.05, if not, execute S1810;

s1810, obtaining eta in sequence_TU=η_TUL+∆η_TUAnd η_TUL =η_TU；

S1811, calculating a turbulence intensity height factor and a turbulence scale height factor by interpolation, and judging whether the turbulence type is rough atmosphere, if so, executing S1813, and if not, executing S1812;

s1812, judging whether the rough atmosphere requirement mark is true, if so, executing S1813, otherwise, executing S1814;

s1813, sequentially calculating the rough atmospheric turbulence intensity and the turbulence scale to obtain the square wave intensity W of the clean cobblestones_COB=0、A_COB=0, and perform S1821;

s1814, judging whether the turbulence type is cobblestone turbulence, if so, executing S1816, and if not, executing S1815;

s1815, judging whether the turbulence type is vertical gust, if so, executing S1816, and if not, executing S1821;

s1816, obtaining cobblestone timing T_CS=T_CSAnd judging the vacuum speed V of the airplane_TRIf it is greater than 700.0, get η_TU=0.01, and executes S1818, and if not, executes S1817;

s1817, judging vacuum speed V_TRIf it is less than 200.0, if so, the rough atmosphere timing T is obtained_P=0.333333, and execute the next step, if not, get T_P=1/［0.006(V_TR200.0) +3.0 ], and the next step is performed;

s1818, judging T_CSWhether or not greater than T_PIf yes, then T is obtained_CS=0, and executing the next step, if not, executing the next step;

s1819, judging T_CSWhether or not less than T_PIf so, obtaining W_COB=A_COBIf not, obtaining W_COB=-A_COB；

S1820, calculating cobblestone turbulence intensity and turbulence scale;

and S1821, exiting the calculation of other turbulence subprocesses.

The sub-process for generating random numbers comprises the following steps:

s1831, entering into a sub-flow of generating random numbers;

s1832, generating 9 random numbers which are uniformly distributed and are between 0 and 1;

s1833, generating 9 random numbers which are uniformly distributed within-0.5 to + 0.5;

s1834, storing 9 normally distributed random numbers in the last cycle;

s1835, calling external search subroutine and interpolating to solve the normal distribution random number R_U1~R_U3、R_V1~R_V3、R_W1~R_W3Wherein R is_U1~R_U3Denotes a random number, R, normally distributed to U_V1~R_V3Denotes a V-normally distributed random number, R_W1~R_W3Representing that W is normally distributed with random numbers;

s1836, quit the generating random number sub-flow.

The filter sub-process includes the following steps:

s1841, entering a filter subprocess;

s1842, according to the time constant T₁=2.0L_U/V_TRAnd gain K =2.8264 η_UL_U/V_TRPerforming a calculation, wherein_UIndicates the intensity of the turbulence in the u direction, L_URepresents the scale of the turbulence in the u direction;

s1843, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U1And a random input variable R_f=R_U1Calculating;

s1844, calling EQN200 sub-flow to obtain F_U1L=F_U1、F_U1=N F₀Gain K =1.0, wherein F_UIL(I =1,2, 3) represents the up-cycle vertical filter output, NF₀Represents the filter output;

s1845, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U2And a random input variable R_f=R_U2Calculating;

s1846, calling EQN200 sub-flow to obtain F_U2L=F_U2、F_U2=N F₀Time constant T₁=L_U/V_TRGain, gain

；

S1847, according to EXP = e^-∆t/T1、F₀=F_U3And R_f=R_U3Calculating;

s1848, calling EQN200 sub-flow to obtain F_U3L=F_U3、F_U3=N F₀Time constant T₁=2.0L_V/V_TRGain K = η_U∙2.5133(L_V/V_TR)²Wherein L is_VRepresents the scale of turbulence in the v-direction;

s1849, according to EXP = e^-∆t/T1Filter output value CF₀=F_V1And a random input variable R_f=R_V1Calculating;

s1850, calling EQN200 sub-flow to obtain F_V10L=F_V1、F_V1=N F₀Time constant T₁=L_V/V_TRAnd a gain K =0.15915, wherein F_VIL(I =1,2, 3) represents the upper-period lateral filter output;

s1851, according to EXP = e^-∆t/T1Filter output value CF₀=F_V2Upper period filter output value DF₀=F_V2LRandom input R_f=R_V1And an upper period random input R_fL=R_V2LCalculating;

s1852, calling EQN300 subflow to obtain F_V2L=F_V2、F_V2=N F₀Time constant 1T₁=L_V/V_TRGain K = η_VL_V/2V_TRAnd time constant 2T₂= 3 ∙ T1, wherein η_VRepresenting the turbulence intensity of the turbulence in the v direction;

s1853, according to EXP = e^-∆t/T1Filter output value CF₀=F_V3Upper period filter output value DF₀=F_V3LRandom input R_f=R_V3And an upper period random input R_fL=R_V3LCalculating;

s1854, calling EQN400 sub-flow to obtain F_V3L=F_V3、F_V3=N F₀Time constant T₁=2.0L_W/V_TRAnd gain K = η_W∙2.5133(L_W/V_TR)²Wherein η_WIndicating the intensity of the turbulence in the w direction, L_WRepresents the scale of the turbulence in the w direction;

s1855, according to EXP = e^-∆t/T1Filter output CF₀=F_W1And a random number R_f=R_W1Calculating;

s1856, calling EQN200 sub-flow to obtain F_W1L=F_W1、F_W1=N F₀Time constant T₁=2.0L_W/V_TRAnd a gain K =0.15915, wherein F_WIL(I =1,2, 3) represents the upper-cycle vertical filter output;

s1857, according to EXP = e^-∆t/T1Filter output value CF₀=F_W2Upper period filter output value DF₀=F_W2LRandom input R_f=R_W2And an upper period random input R_fL=R_W2LCalculating;

s1858, calling EQN300 subflow to obtain F_W2L=F_W2、F_W2=N F₀Time constant 1T₁=L_W/V_TRGain K = η_WL_W/2V_TRAnd time constant 2T₂=√3∙T1；

S1859, according to EXP = e^-∆t/T1Filter output value CF₀=F_W3Upper period filter output DF₀=F_W3LA random number R_f=R_W3And an upper periodic random number R_fL=R_W3LCalculating;

s1860, call EQN400 sub-flow to get F_W3L=F_W3And F_W3=N F₀；

S1861, quitting the filter subprocess.

Specifically, the EQN200 sub-process includes the following steps:

proceed to EQN200 subroutine;

calculate NF₀=CF₀ e^-∆t/T1+R_fK(1- e^-∆t/T1)；

Judgment | NF₀If | is less than 0.00001, if yes, NF is obtained₀=0, if not, executing the next step;

the subroutine EQN200 is exited.

EQN300 the subflow includes the following steps:

proceed to EQN300 subroutine;

calculate NF₀=(2CF₀-DF₀∙EXP)∙EXP+∆t∙［K∙EXP(R_f-R_fL)/T₁ ²］;

Judgment | NF₀If | is less than 0.00001, if yes, NF is obtained₀=0, if not, executing the next step；

The subroutine EQN300 is exited.

EQN400 the subroutine includes the following steps:

proceed to EQN400 subroutine;

sequentially calculate F₁=(2.0∙CF₀-DF₀∙EXP)∙EXP，F₂₁=KT₂/T₁ ²，F₂₂= T₁ ²/T₂，F₂₃=(∆tT₁/T₂₎ -∆t，F₃₁=R_f∙F₂₁∙［F₂₂+(-F₂₃-F₂₂)∙EXP］，F₃₂=R_fL∙F₂₁∙［F₂₂(EXP-1.0)+F₂₃］∙EXP，NF₀=F₁+F₃₁+F₃₂；

Judgment | NF₀If | is less than 0.00001, if yes, NF is obtained₀=0, if not, executing the next step;

the subroutine EQN400 is exited.

The turbulence rate calculating sub-process comprises the following steps:

s1871, entering a sub-process for calculating the turbulent flow rate;

s1872, judging whether the turbulence type is vertical gust, and if so, obtaining a turbulence factor eta_THL=1.0, and execute S1878, if not, execute S1873;

s1873 and judging the flying height h_RWhether or not less than the turbulent low-height limit h_L-0.5, if yes, perform S1875, if no, perform S1874;

s1874 and judging h_RWhether or not greater than the turbulent high height limit h_U+0.5, if yes, execute S1875, if no, execute S1876;

s1875 obtaining a turbulence factor eta_THL=0, and perform S1878;

s1876 and judging h_RWhether or not h is less than or equal to_LIf yes, get η_THL=1.0-(h_L-h_R) ∙ 2.0.0, and executing S1878, if not, executing S1877;

s1877 and judging h_RWhether or not h is greater than or equal to_UIf yes, get η_THL=1.0+(h_U-h_R) ∙ 2.0.0, and executing S1878, if not, obtaining eta_THL=1.0, and perform S1878;

s1878, calculating U in sequence_RAL=F_U1∙F_U2+F_U3、U_RAL=F_V1∙F_V2+F_V3、w_RAL= F_W1∙F_W2+F_W3∙W_JETUP+W_COB、P_RAL=0.025l(W_RAL-W_JETUP)/L_W、Q_RAL=0.025(W_RAL-W_JETUP)81.5/L_W、R_RAL=0.025∙81.5V_RAL/L_W、U_T=U_RAL∙η_THL、V_T=V_RAL∙η_THL、W_T=W_RAL∙η_THL、P_T=P_RAL∙η_THL、Q_T=Q_RAL∙η_THL、R_T=R_RAL∙η_THL、

、

、

、U_TOL=U_T、V_TOL=V_TAnd W_TOL=W_TWherein P is_TRepresenting the DRYDEN model output pitch angular velocity, Q_TRepresenting DRYDEN model output yaw rate, R_TRepresenting the DRYDEN model output roll angular velocity, U_TLRepresents the X-direction speed, V, of the DRYDEN model output in the upper period_TLRepresents the Y-direction speed, W of the DRYDEN model output in the upper period_TLRepresents the Z-direction speed of the DRYDEN model output in the upper period,

represents;

represents the acceleration of the output Y direction of the DRYDEN model,

represents the Z-direction acceleration of the DRYDEN model output;

s1879 and exiting the turbulence rate calculation sub-process.

S19, judging the time T for repositioning the airplane_REPIf not, executing S21 if yes, and executing S20 if no;

s20, judging whether the vacuum velocity VTR of the airplane is less than 15.0, if so, executing S21, and if not, exiting the turbulence module;

s21, sequentially obtaining the clean turbulence velocity UT = VT = WT =0, UT = VT = WT =0 and eta_TU= η_TUL= η_TUD=0, exiting the turbulence module;

wherein UT is the X-direction speed of the DRYDEN model output, VT is the Y-direction speed of the DRYDEN model output, WT is the Z-direction speed of the DRYDEN model output, and eta is_TUIs the turbulence intensity factor, η_TULIs the upper-period turbulence intensity factor, η_TUDIs the rough atmospheric turbulence intensity factor.

For the convenience of understanding the above technical solution of the present invention, the following detailed description will be made about the relevant calculation formula in the turbulence module of the present invention.

The turbulence module has the following related calculation formula:

wherein, the turbulent module formula variable corresponding table is as follows:

in summary, by means of the technical scheme of the invention, turbulence models such as rough atmosphere, cobblestone turbulence, discrete vertical gust and the like can be provided through the use of the invention, the turbulence intensity is 10-level, the gust intensity can be selected within the range of 0-25 m/s, and the atmospheric disturbance environment is simulated by calculating the turbulence disturbance speed. In addition, the invention can be divided into vertical gust turbulence, rough atmosphere and cobblestone turbulence according to the type of the turbulence, and can be divided into nine types of simulated turbulence of light turbulence, moderate turbulence and severe turbulence according to the intensity of the turbulence.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (10)

1. A simulation method for simulating a turbulent environment of an aircraft is characterized by comprising the following steps:

s1, entering a turbulence module, judging whether the speed resetting mark of the airplane is true, if so, executing S19, and if not, executing S2;

s2, judging whether the airplane flight freezing mark is true, if so, executing S19, and if not, executing S3;

s3, judging whether the vacuum speed VTR of the airplane is larger than 15.0, if so, executing S4, and if not, executing S19;

s4, judging whether the airplane reset directory is 0, if so, executing S5, and if not, executing S19;

s5, judging whether the aircraft wind profile index is not 0, if so, executing S7, and if not, executing S6;

s6, judging whether the airplane micro-storm activation mark is true, if so, executing S7, otherwise, setting the airplane rough atmosphere requirement mark as false, and executing S8;

s7, setting the requirement rough atmosphere mark of the airplane to be true, and executing S8;

s8, judging whether the aircraft rough atmosphere mark is true, if so, executing S9, otherwise, executing S10;

s9, judging whether the type of the aircraft turbulence is 0, if so, executing S12, and if not, executing S10;

s10, judging whether the turbulence intensity of the airplane is 0, if so, executing S12, and if not, executing S11;

s11, judging whether the airplane relocation directory is not 0, if so, executing S12, and if not, executing S13;

s12, judging the aircraft turbulence intensity factor eta_TUIf the value is 0, executing S14 if the value is 0, and executing S13 if the value is not 0;

s13, judging whether the type of the turbulent flow of the airplane is not 3, if so, calling a vertical gust removal sub-process and executing S15, and if not, calling a vertical gust calculation sub-process and executing S15;

s14, calling a turbulence eliminating sub-process, then calling a vertical gust eliminating sub-process, and exiting the turbulence module;

s15, judging aircraft gust timing T_JETUPWhether the time is less than 30.0 seconds, if yes, executing S17, and if not, executing S16;

s16, judging whether the type of the aircraft turbulence is 3, if so, executing S17, and if not, executing S18;

s17, calling a turbulence removal sub-process, and exiting from the turbulence module;

s18, sequentially calling and calculating other turbulent flow sub-processes, calling and generating a random number sub-process, calling a filter sub-process, calling and calculating a turbulent flow rate sub-process, and exiting from the turbulent flow module;

s19, judging the time T for repositioning the airplane_REPIf not, executing S21 if yes, and executing S20 if no;

s20, judging whether the vacuum velocity VTR of the airplane is less than 15.0, if so, executing S21, and if not, exiting the turbulence module;

s21, sequentially obtaining the clean turbulence velocity UT = VT = WT =0, UT = VT = WT =0 and eta_TU= η_TUL= η_TUD=0, exiting the turbulence module;

wherein UT is the X-direction speed of the DRYDEN model output, VT is the Y-direction speed of the DRYDEN model output, WT is the Z-direction speed of the DRYDEN model output, and eta is_TUIs the turbulence intensity factor, η_TULIs the upper-period turbulence intensity factor, η_TUDIs the rough atmospheric turbulence intensity factor.

2. The simulation method for simulating the turbulent environment of the airplane as claimed in claim 1, wherein the process of eliminating the vertical gust in S13 includes the following steps:

s1301, entering a process of removing vertical gusts;

s1302, setting a gust distance factor and a gust timing to be 0;

s1303, setting gust intensity factors and vertical gust intensities to be 0;

s1304, setting gust to require that the cobblestone turbulence intensity and the vertical speed of the gust are both 0;

and S1305, exiting the vertical gust cleaning flow.

3. The simulation method for simulating the turbulent environment of the airplane as claimed in claim 1, wherein the step of calculating the vertical gust sub-flow in S13 includes the following steps:

s1311, entering a vertical gust calculating sub-process;

s1312, judging whether the turbulence intensity signal is not 0, if so, obtaining a turbulence intensity factor of the vertical gust

Wherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtained_JETLOC=0；

S1313, judging aircraft gust timing T_JETUPIf less than 30.0, if yes, T is obtained_JETUP=T_JETUPΔ t, and execute S1314, wherein t represents the timing interval, if not, the distance factor =0.008Mn + distance factor η_COBLOC=0, and S1315 is performed, where Mn denotes the flight mach number, η_COBLOCRepresenting a cobblestone turbulence intensity factor;

s1314, judging T_JETUPWhether the gust turbulence intensity is larger than 40.0 or not is judged, if yes, the cobble turbulence intensity required by gust turbulence is calculated, S1317 is executed, and if not, S1317 is executed;

s1315, judging whether the distance factor is larger than or equal to 1.0, if so, executing S1316, and if not, executing S1317;

s1316, setting the turbulence type number as 0 and the turbulence intensity number as 0 in sequence, requiring the gust mark as false, the mild turbulence mark as false, the moderate turbulence mark as false and the severe turbulence mark as false, and executing S1317;

s1317, calculating gust height factor to vertical gust height factor eta by interpolation_JUM=f(h₁) Wherein h is₁Represents a height;

s1318, obtaining vertical gust speed W_JETUP=f(h₂)∙η_JETLOC∙η_JUMWherein h is₂Represents a distance factor;

and S1319, exiting the vertical gust calculating sub-process.

4. The simulation method for aircraft turbulent environment simulation according to claim 1, wherein the turbulence removal sub-process in S14 includes the following steps:

s141, entering a turbulence removal sub-process;

s142, calling a first fade-out sub-process (U)_RAL,U_RA) And a second fade-out sub-flow (V)_RAL,V_RA) Wherein, U_RALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine system_RARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine system_RALRepresents the lateral component factor, V, of the turbulent wind speed under the machine system_RARepresenting the lateral component of the turbulent wind speed under the machine system;

s143, calling a third fade-out sub-flow (W)_RAL,W_RA) And a fourth fade-out sub-flow (P)_RAL,P_RA) Wherein W is_RALRepresents the vertical component factor, W, of the turbulent wind speed under the machine system_RARepresents the vertical component, P, of the turbulent wind speed under the machine system_RALPitch angle velocity factor, P, representing turbulent wind velocity generation_RARepresenting the pitch angle velocity resulting from turbulent wind speed;

s144, calling a fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA) Wherein Q is_RALIndicating deviation of turbulent wind velocity generationAngular velocity factor, Q_RARepresenting yaw rate, R, resulting from turbulent wind speed_RALRoll angular velocity factor, R, representing turbulent wind velocity generation_RARepresenting the roll angular velocity generated by turbulent wind speed;

s145, obtaining eta_TU= η_TUL= η_TUD=0；

S146, clearing all filter outputs F_UI=0、F_VI=0、F_WI=0, I =1,2,3, wherein F_UIRepresenting the longitudinal filter output, F_VIRepresenting the lateral filter output, F_WIRepresenting the vertical filter output;

and S147, exiting the turbulence removing sub-process.

5. A simulation method for aircraft turbulence environment simulation according to claim 4, characterized in that the fade-out sub-routine comprises a first fade-out sub-routine (U)_RAL,U_RA) And a second decreasing subprocess (V)_RAL,V_RA) And the third fade-out subprocess (W)_RAL,W_RA) And the fourth fade-out sub-process (P)_RAL,P_RA) The fifth fade-out sub-process (Q)_RAL,Q_RA) And a sixth fade-out sub-flow (R)_RAL,R_RA）。

6. The simulation method for aircraft turbulent environment simulation according to claim 5, wherein the fade-out sub-process comprises the following steps:

entering a fade-out sub-process, and calculating by using a formula RL = RAKFADER RC, wherein RL represents U_RAL、V_RAL、W_RAL、P_RAL、Q_RAL、R_RALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents U_RA、V_RA、W_RA、P_RA、Q_RA、R_RA；

And judging whether the | RL | is smaller than 0.01, if so, obtaining RC =0, and exiting the fade-out sub-process, and if not, obtaining RC = RL, and exiting the fade-out sub-process.

7. The simulation method for aircraft turbulent environment simulation according to claim 1, wherein the step of calculating other turbulent sub-processes in S18 comprises the following steps:

s1801, entering a sub-process of calculating other turbulence;

s1802, judging whether the rough atmosphere requirement mark is true, if so, calculating a turbulence intensity factor required by wind, and executing S1806, otherwise, executing S1803;

s1803, judging whether the turbulence type is rough atmosphere, if so, obtaining turbulence intensity factor required by the rough atmosphere

And executing S1806, if not, executing S1804;

s1804, judging whether the turbulence type is cobblestone turbulence, if so, obtaining cobblestone atmosphere required turbulence intensity factors

And executing S1805, if not, obtaining eta_TUD=η_COBLOCAnd performs S1806;

s1805, obtaining the intensity of the square wave

And performs S1806;

s1806, judging whether the turbulence intensity signal is 0, if so, executing S1807, and if not, executing 1808;

s1807, setting the rough atmosphere requirement flag to true, and if the rough atmosphere requirement flag is set, obtaining eta_TUD=0, and executing S1808, if not, executing S1808;

s1808 obtaining the increment of the turbulent intensity factor [. eta. ]_TU=η_TUD-η_TULEqual and equal to_TUIf it is greater than 0.01, get [. eta. ]_TUAnd =0.01, and execute S1810, if no, execute S1809;

s1809, Δ η%_TUWhether or not less than-0.05, ifIf it is, then Δ η_TU= -0.05, if not, execute S1810;

s1810, obtaining eta in sequence_TU=η_TUL+∆η_TUAnd η_TUL =η_TU；

S1811, calculating a turbulence intensity height factor and a turbulence scale height factor by interpolation, and judging whether the turbulence type is rough atmosphere, if so, executing S1813, and if not, executing S1812;

s1812, judging whether the rough atmosphere requirement mark is true, if so, executing S1813, otherwise, executing S1814;

s1813, sequentially calculating the rough atmospheric turbulence intensity and the turbulence scale to obtain the square wave intensity W of the clean cobblestones_COB=0、A_COB=0, and perform S1821;

s1814, judging whether the turbulence type is cobblestone turbulence, if so, executing S1816, and if not, executing S1815;

s1815, judging whether the turbulence type is vertical gust, if so, executing S1816, and if not, executing S1821;

s1816, obtaining cobblestone timing T_CS=T_CSAnd judging the vacuum speed V of the airplane_TRIf it is greater than 700.0, get η_TU=0.01, and executes S1818, and if not, executes S1817;

s1817, judging vacuum speed V_TRIf it is less than 200.0, if so, the rough atmosphere timing T is obtained_P=0.333333, and execute the next step, if not, get T_P=1/［0.006(V_TR200.0) +3.0 ], and the next step is performed;

s1818, judging T_CSWhether or not greater than T_PIf yes, then T is obtained_CS=0, and executing the next step, if not, executing the next step;

s1819, judging T_CSWhether or not less than T_PIf so, obtaining W_COB=A_COBIf not, obtaining W_COB=-A_COB；

S1820, calculating cobblestone turbulence intensity and turbulence scale;

and S1821, exiting the calculation of other turbulence subprocesses.

8. The simulation method for simulating the turbulent environment of the airplane as claimed in claim 1, wherein the step of generating the random number subprocess in S18 includes the following steps:

s1831, entering into a sub-flow of generating random numbers;

s1832, generating 9 random numbers which are uniformly distributed and are between 0 and 1;

s1833, generating 9 random numbers which are uniformly distributed within-0.5 to + 0.5;

s1834, storing 9 normally distributed random numbers in the last cycle;

s1835, calling external search subroutine and interpolating to solve the normal distribution random number R_U1~R_U3、R_V1~R_V3、R_W1~R_W3Wherein R is_U1~R_U3Denotes a random number, R, normally distributed to U_V1~R_V3Denotes a V-normally distributed random number, R_W1~R_W3Representing that W is normally distributed with random numbers;

s1836, quit the generating random number sub-flow.

9. The simulation method for simulating the turbulent environment of an aircraft according to claim 1, wherein the filter subprocess in S18 includes the following steps:

s1841, entering a filter subprocess;

s1842, according to the time constant T₁=2.0L_U/V_TRAnd gain K =2.8264 η_UL_U/V_TRPerforming a calculation, wherein_UIndicates the intensity of the turbulence in the u direction, L_URepresents the scale of the turbulence in the u direction;

s1843, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U1And a random input variable R_f=R_U1Calculating;

s1844, calling EQN200 sub-flow to obtain F_U1L=F_U1、F_U1=N F₀Gain K =1.0, wherein,F_UIL(I =1,2, 3) represents the up-cycle vertical filter output, NF₀Represents the filter output;

s1845, according to EXP = e^-∆t/T1Initial value CF of filter₀=F_U2And a random input variable R_f=R_U2Calculating;

s1846, calling EQN200 sub-flow to obtain F_U2L=F_U2、F_U2=N F₀Time constant T₁=L_U/V_TRGain, gain

；

S1847, according to EXP = e^-∆t/T1、F₀=F_U3And R_f=R_U3Calculating;

s1848, calling EQN200 sub-flow to obtain F_U3L=F_U3、F_U3=N F₀Time constant T₁=2.0L_V/V_TRGain K = η_U∙2.5133(L_V/V_TR)²Wherein L is_VRepresents the scale of turbulence in the v-direction;

s1849, according to EXP = e^-∆t/T1Filter output value CF₀=F_V1And a random input variable R_f=R_V1Calculating;

s1850, calling EQN200 sub-flow to obtain F_V10L=F_V1、F_V1=N F₀Time constant T₁=L_V/V_TRAnd a gain K =0.15915, wherein F_VIL(I =1,2, 3) represents the upper-period lateral filter output;

s1851, according to EXP = e^-∆t/T1Filter output value CF₀=F_V2Upper period filter output value DF₀=F_V2LRandom input R_f=R_V1And an upper period random input R_fL=R_V2LCalculating;

s1852, calling EQN300 subflow to obtain F_V2L=F_V2、F_V2=N F₀Time constant 1T₁=L_V/V_TRGain K = η_VL_V/2V_TRAnd time constant 2T₂= 3 ∙ T1, wherein η_VRepresenting the turbulence intensity of the turbulence in the v direction;

s1853, according to EXP = e^-∆t/T1Filter output value CF₀=F_V3Upper period filter output value DF₀=F_V3LRandom input R_f=R_V3And an upper period random input R_fL=R_V3LCalculating;

s1854, calling EQN400 sub-flow to obtain F_V3L=F_V3、F_V3=N F₀Time constant T₁=2.0L_W/V_TRAnd gain K = η_W∙2.5133(L_W/V_TR)²Wherein η_WIndicating the intensity of the turbulence in the w direction, L_WRepresents the scale of the turbulence in the w direction;

s1855, according to EXP = e^-∆t/T1Filter output CF₀=F_W1And a random number R_f=R_W1Calculating;

s1856, calling EQN200 sub-flow to obtain F_W1L=F_W1、F_W1=N F₀Time constant T₁=2.0L_W/V_TRAnd a gain K =0.15915, wherein F_WIL(I =1,2, 3) represents the upper-cycle vertical filter output;

s1857, according to EXP = e^-∆t/T1Filter output value CF₀=F_W2Upper period filter output value DF₀=F_W2LRandom input R_f=R_W2And an upper period random input R_fL=R_W2LCalculating;

s1858, calling EQN300 subflow to obtain F_W2L=F_W2、F_W2=N F₀Time constant 1T₁=L_W/V_TRGain K = η_WL_W/2V_TRAnd time constant 2T₂=√3∙T1；

S1859, according to EXP = e^-∆t/T1Filter output value CF₀=F_W3Upper period filter output DF₀=F_W3LA random number R_f=R_W3And an upper periodic random number R_fL=R_W3LCalculating;

s1860, call EQN400 sub-flow to get F_W3L=F_W3And F_W3=N F₀；

S1861, quitting the filter subprocess.

10. The simulation method for aircraft turbulent environment simulation according to claim 1, wherein the step of calculating the turbulent velocity subprocess in S18 comprises the following steps:

s1871, entering a sub-process for calculating the turbulent flow rate;

s1872, judging whether the turbulence type is vertical gust, and if so, obtaining a turbulence factor eta_THL=1.0, and execute S1878, if not, execute S1873;

s1873 and judging the flying height h_RWhether or not less than the turbulent low-height limit h_L-0.5, if yes, perform S1875, if no, perform S1874;

s1874 and judging h_RWhether or not it is greater than h_U+0.5, if yes, execute S1875, if no, execute S1876;

s1875 obtaining a turbulence factor eta_THL=0, and perform S1878;

s1876 and judging h_RWhether or not h is less than or equal to_LIf yes, get η_THL=1.0-(h_L-h_R) ∙ 2.0.0, and executing S1878, if not, executing S1877;

s1877 and judging h_RWhether or not the height of the turbulent flow is greater than or equal to the limit h_UIf yes, get η_THL=1.0+(h_U-h_R) ∙ 2.0.0, and executing S1878, if not, obtaining eta_THL=1.0, and perform S1878;

s1878, calculating U in sequence_RAL=F_U1∙F_U2+F_U3、U_RAL=F_V1∙F_V2+F_V3、w_RAL= F_W1∙F_W2+F_W3∙W_JETUP+W_COB、P_RAL=0.025l(W_RAL-W_JETUP)/L_W、Q_RAL=0.025(W_RAL-W_JETUP)81.5/L_W、R_RAL=0.025∙81.5V_RAL/L_W、U_T=U_RAL∙η_THL、V_T=V_RAL∙η_THL、W_T=W_RAL∙η_THL、P_T=P_RAL∙η_THL、Q_T=Q_RAL∙η_THL、R_T=R_RAL∙η_THL、

、

、

、U_TOL=U_T、V_TOL=V_TAnd W_TOL=W_TWherein P is_TRepresenting the DRYDEN model output pitch angular velocity, Q_TRepresenting DRYDEN model output yaw rate, R_TRepresenting the DRYDEN model output roll angular velocity, U_TLRepresents the X-direction speed, V, of the DRYDEN model output in the upper period_TLRepresents the Y-direction speed, W of the DRYDEN model output in the upper period_TLRepresents the Z-direction speed of the DRYDEN model output in the upper period,

represents;

represents the acceleration of the output Y direction of the DRYDEN model,

represents the Z-direction acceleration of the DRYDEN model output;

s1879 and exiting the turbulence rate calculation sub-process.

Images