CN113239462A - Simulation method for aircraft turbulent environment simulation - Google Patents
Simulation method for aircraft turbulent environment simulation Download PDFInfo
- Publication number
- CN113239462A CN113239462A CN202110568308.6A CN202110568308A CN113239462A CN 113239462 A CN113239462 A CN 113239462A CN 202110568308 A CN202110568308 A CN 202110568308A CN 113239462 A CN113239462 A CN 113239462A
- Authority
- CN
- China
- Prior art keywords
- turbulence
- executing
- ral
- sub
- judging whether
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 123
- 238000004088 simulation Methods 0.000 title claims abstract description 41
- 230000004913 activation Effects 0.000 claims abstract description 4
- 230000008014 freezing Effects 0.000 claims abstract description 4
- 238000007710 freezing Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000004140 cleaning Methods 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
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 etaTUIf 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 TJETUPWhether 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 airplaneREPIf 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 etaTU= η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 gustWherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtainedJETLOC=0;
S1313, judging aircraft gust timing TJETUPIf less than 30.0, if yes, T is obtainedJETUP=TJETUPΔ 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 TJETUPWhether 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(h1) Wherein h is1Represents a height;
s1318, obtaining vertical gust speed WJETUP=f(h2)∙ηJETLOC∙ηJUMWherein h is2Represents 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,URA) And a second fade-out sub-flow (V)RAL,VRA) Wherein, URALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine systemRARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine systemRALRepresents the lateral component factor, V, of the turbulent wind speed under the machine systemRARepresenting the lateral component of the turbulent wind speed under the machine system;
s143, calling a third fade-out sub-flow (W)RAL,WRA) And a fourth fade-out sub-flow (P)RAL,PRA) Wherein W isRALRepresents the vertical component factor, W, of the turbulent wind speed under the machine systemRARepresents the vertical component, P, of the turbulent wind speed under the machine systemRALPitch angle velocity factor, P, representing turbulent wind velocity generationRARepresenting the pitch angle velocity resulting from turbulent wind speed;
s144, calling a fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA) Wherein Q isRALYaw rate factor, Q, representing turbulent wind velocityRARepresenting yaw rate, R, resulting from turbulent wind speedRALRoll angular velocity factor, R, representing turbulent wind velocity generationRARepresenting the roll angular velocity generated by turbulent wind speed;
s145, obtaining etaTU= ηTUL= ηTUD=0;
S146, clearing all filter outputs FUI=0、FVI=0、FWI=0, I =1,2,3, wherein FUIRepresenting the longitudinal filter output, FVIRepresenting the lateral filter output, FWIRepresenting 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,URA) And a second decreasing subprocess (V)RAL,VRA) And the third fade-out subprocess (W)RAL,WRA) And the fourth fade-out sub-process (P)RAL,PRA) The fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA)。
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 URAL、VRAL、WRAL、PRAL、QRAL、RRALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents URA、VRA、WRA、PRA、QRA、RRA;
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 atmosphereAnd executing S1806, if not, executing S1804;
s1804, judging whether the turbulence type is cobblestone turbulence, if so, obtaining cobblestone atmosphere required turbulence intensity factorsAnd executing S1805, if not, obtaining etaTUD=ηCOBLOCAnd 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 etaTUD=0, and executing S1808, if not, executing S1808;
s1808 obtaining the increment of the turbulent intensity factor [. eta. ]TU=ηTUD-ηTULEqual and equal toTUIf 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 sequenceTU=η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 cobblestonesCOB=0、ACOB=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 TCS=TCSAnd judging the vacuum speed V of the airplaneTRIf it is greater than 700.0, get ηTU=0.01, and executes S1818, and if not, executes S1817;
s1817, judging vacuum speed VTRIf it is less than 200.0, if so, the rough atmosphere timing T is obtainedP=0.333333, and execute the next step, if not, get TP=1/[0.006(VTR200.0) +3.0 ], and the next step is performed;
s1818, judging TCSWhether or not greater than TPIf yes, then T is obtainedCS=0, and executing the next step, if not, executing the next step;
s1819, judging TCSWhether or not less than TPIf so, obtaining WCOB=ACOBIf not, obtaining WCOB=-ACOB;
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 RU1~RU3、RV1~RV3、RW1~RW3Wherein R isU1~RU3Denotes a random number, R, normally distributed to UV1~RV3Denotes a V-normally distributed random number, RW1~RW3Representing 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 T1=2.0LU/VTRAnd gain K =2.8264 ηULU/VTRPerforming a calculation, whereinUIndicates the intensity of the turbulence in the u direction, LURepresents the scale of the turbulence in the u direction;
s1843, according to EXP = e-∆t/T1Initial value CF of filter0=FU1And a random input variable Rf=RU1Calculating;
s1844, calling EQN200 sub-flow to obtain FU1L=FU1、FU1=N F0Gain K =1.0, wherein FUIL(I =1,2, 3) represents the up-cycle vertical filter output, NF0Represents the filter output;
s1845, according to EXP = e-∆t/T1Initial value CF of filter0=FU2And a random input variable Rf=RU2Calculating;
S1847, according to EXP = e-∆t/T1、F0=FU3And Rf=RU3Calculating;
s1848, calling EQN200 sub-flow to obtain FU3L=FU3、FU3=N F0Time constant T1=2.0LV/VTRGain K = ηU∙2.5133(LV/VTR)2Wherein L isVRepresents the scale of turbulence in the v-direction;
s1849, according to EXP = e-∆t/T1Filter output value CF0=FV1And a random input variable Rf=RV1Calculating;
s1850, calling EQN200 sub-flow to obtain FV10L=FV1、FV1=N F0Time constant T1=LV/VTRAnd a gain K =0.15915, wherein FVIL(I =1,2, 3) represents the upper-period lateral filter output;
s1851, according to EXP = e-∆t/T1Filter output value CF0=FV2Upper period filter output value DF0=FV2LRandom input Rf=RV1And an upper period random input RfL=RV2LCalculating;
s1852, calling EQN300 subflow to obtain FV2L=FV2、FV2=N F0Time constant 1T1=LV/VTRGain K = ηVLV/2VTRAnd time constant 2T2= 3 ∙ T1, wherein ηVRepresenting the turbulence intensity of the turbulence in the v direction;
s1853, according to EXP = e-∆t/T1Filter output value CF0=FV3Upper period filter output value DF0=FV3LRandom input Rf=RV3And an upper period random input RfL=RV3LCalculating;
s1854, calling EQN400 sub-flow to obtain FV3L=FV3、FV3=N F0Time constant T1=2.0LW/VTRAnd gain K = ηW∙2.5133(LW/VTR)2Wherein ηWIndicating the intensity of the turbulence in the w direction, LWRepresents the scale of the turbulence in the w direction;
s1855, according to EXP = e-∆t/T1Filter output CF0=FW1And a random number Rf=RW1Calculating;
s1856, calling EQN200 sub-flow to obtain FW1L=FW1、FW1=N F0Time and time constantNumber T1=2.0LW/VTRAnd a gain K =0.15915, wherein FWIL(I =1,2, 3) represents the upper-cycle vertical filter output;
s1857, according to EXP = e-∆t/T1Filter output value CF0=FW2Upper period filter output value DF0=FW2LRandom input Rf=RW2And an upper period random input RfL=RW2LCalculating;
s1858, calling EQN300 subflow to obtain FW2L=FW2、FW2=N F0Time constant 1T1=LW/VTRGain K = ηWLW/2VTRAnd time constant 2T2=√3∙T1;
S1859, according to EXP = e-∆t/T1Filter output value CF0=FW3Upper period filter output DF0=FW3LA random number Rf=RW3And an upper periodic random number RfL=RW3LCalculating;
s1860, call EQN400 sub-flow to get FW3L=FW3And FW3=N F0;
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 etaTHL=1.0, and execute S1878, if not, execute S1873;
s1873 and judging the flying height hRWhether or not less than the turbulent low-height limit hL-0.5, if yes, perform S1875, if no, perform S1874;
s1874 and judging hRWhether or not it is greater than hU+0.5, if yes, execute S1875, if no, execute S1876;
s1875 obtaining a turbulence factor etaTHL=0, and perform S1878;
s1876 and judging hRWhether or not h is less than or equal toLIf yes, get ηTHL=1.0-(hL-hR) ∙ 2.0.0, and executing S1878, if not, executing S1877;
s1877 and judging hRWhether or not the height of the turbulent flow is greater than or equal to the limit hUIf yes, get ηTHL=1.0+(hU-hR) ∙ 2.0.0, and executing S1878, if not, obtaining etaTHL=1.0, and perform S1878;
s1878, calculating U in sequenceRAL=FU1∙FU2+FU3、URAL=FV1∙FV2+FV3、wRAL= FW1∙FW2+FW3∙WJETUP+WCOB、PRAL=0.025l(WRAL-WJETUP)/LW、QRAL=0.025(WRAL-WJETUP)81.5/LW、RRAL=0.025∙81.5VRAL/LW、UT=URAL∙ηTHL、VT=VRAL∙ηTHL、WT=WRAL∙ηTHL、PT=PRAL∙ηTHL、QT=QRAL∙ηTHL、RT=RRAL∙ηTHL、、、、UTOL=UT、VTOL=VTAnd WTOL=WTWherein P isTRepresenting the DRYDEN model output pitch angular velocity, QTRepresenting DRYDEN model output yaw rate, RTRepresenting the DRYDEN model output roll angular velocity, UTLRepresents the X-direction speed, V, of the DRYDEN model output in the upper periodTLRepresents the Y-direction speed, W of the DRYDEN model output in the upper periodTLTo 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 etaTUIf 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 gustWherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtainedJETLOC=0;
S1313, judging aircraft gust timing TJETUPIf less than 30.0, if yes, T is obtainedJETUP=TJETUPΔ 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 TJETUPWhether 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 interpolationJUM=f(h1) Wherein h is1Represents a height;
s1318, obtaining vertical gust speed WJETUP=f(h2)∙ηJETLOC∙ηJUMWherein h is2Represents 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,URA) And a second fade-out sub-flow (V)RAL,VRA) Wherein, URALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine systemRARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine systemRALRepresents the lateral component factor, V, of the turbulent wind speed under the machine systemRARepresenting the lateral component of the turbulent wind speed under the machine system;
s143, calling a third fade-out sub-flow (W)RAL,WRA) And a fourth fade-out sub-flow (P)RAL,PRA) Wherein W isRALRepresents the vertical component factor, W, of the turbulent wind speed under the machine systemRARepresents the vertical component, P, of the turbulent wind speed under the machine systemRALPitch angle velocity factor, P, representing turbulent wind velocity generationRARepresenting the pitch angle velocity resulting from turbulent wind speed;
s144, calling a fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA) Wherein Q isRALYaw rate factor, Q, representing turbulent wind velocityRARepresenting yaw rate, R, resulting from turbulent wind speedRALRoll angular velocity factor, R, representing turbulent wind velocity generationRARepresenting the roll angular velocity generated by turbulent wind speed;
s145, obtaining etaTU= ηTUL= ηTUD=0;
S146, clearing all filter outputs FUI=0、FVI=0、FWI=0, I =1,2,3, wherein FUIRepresenting the longitudinal filter output, FVIRepresenting the lateral filter output, FWIRepresenting 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,URA) And a second decreasing subprocess (V)RAL,VRA) And the third fade-out subprocess (W)RAL,WRA) And the fourth fade-out sub-process (P)RAL,PRA) The fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA)。
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 URAL、VRAL、WRAL、PRAL、QRAL、RRALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents URA、VRA、WRA、PRA、QRA、RRA;
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 TJETUPWhether 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 atmosphereAnd 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 turbulenceAnd executing S1805, if not, obtaining etaTUD=ηCOBLOCAnd 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 etaTUD=0, and executing S1808, if not, executing S1808;
s1808 obtaining the increment of the turbulent intensity factor [. eta. ]TU=ηTUD-ηTULEqual and equal toTUIf 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 sequenceTU=η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 cobblestonesCOB=0、ACOB=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 TCS=TCSAnd judging the vacuum speed V of the airplaneTRIf it is greater than 700.0, get ηTU=0.01, and executes S1818, and if not, executes S1817;
s1817, judging vacuum speed VTRIf it is less than 200.0, if so, the rough atmosphere timing T is obtainedP=0.333333, and execute the next step, if not, get TP=1/[0.006(VTR200.0) +3.0 ], and the next step is performed;
s1818, judging TCSWhether or not greater than TPIf yes, then T is obtainedCS=0, and executing the next step, if not, executing the next step;
s1819, judging TCSWhether or not less than TPIf so, obtaining WCOB=ACOBIf not, obtaining WCOB=-ACOB;
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 RU1~RU3、RV1~RV3、RW1~RW3Wherein R isU1~RU3Denotes a random number, R, normally distributed to UV1~RV3Denotes a V-normally distributed random number, RW1~RW3Representing 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 T1=2.0LU/VTRAnd gain K =2.8264 ηULU/VTRPerforming a calculation, whereinUIndicates the intensity of the turbulence in the u direction, LURepresents the scale of the turbulence in the u direction;
s1843, according to EXP = e-∆t/T1Initial value CF of filter0=FU1And a random input variable Rf=RU1Calculating;
s1844, calling EQN200 sub-flow to obtain FU1L=FU1、FU1=N F0Gain K =1.0, wherein FUIL(I =1,2, 3) represents the up-cycle vertical filter output, NF0Represents the filter output;
s1845, according to EXP = e-∆t/T1Initial value CF of filter0=FU2And a random input variable Rf=RU2Calculating;
S1847, according to EXP = e-∆t/T1、F0=FU3And Rf=RU3Calculating;
s1848, calling EQN200 sub-flow to obtain FU3L=FU3、FU3=N F0Time constant T1=2.0LV/VTRGain K = ηU∙2.5133(LV/VTR)2Wherein L isVRepresents the scale of turbulence in the v-direction;
s1849, according to EXP = e-∆t/T1Filter output value CF0=FV1And a random input variable Rf=RV1Calculating;
s1850, calling EQN200 sub-flow to obtain FV10L=FV1、FV1=N F0Time constant T1=LV/VTRAnd a gain K =0.15915, wherein FVIL(I =1,2, 3) represents the upper-period lateral filter output;
s1851, according to EXP = e-∆t/T1Filter output value CF0=FV2Upper period filter output value DF0=FV2LRandom input Rf=RV1And an upper period random input RfL=RV2LCalculating;
s1852, calling EQN300 subflow to obtain FV2L=FV2、FV2=N F0Time constant 1T1=LV/VTRGain K = ηVLV/2VTRAnd time constant 2T2= 3 ∙ T1, wherein ηVRepresenting the turbulence intensity of the turbulence in the v direction;
s1853, according to EXP = e-∆t/T1Filter output value CF0=FV3Upper period filter output value DF0=FV3LRandom input Rf=RV3And an upper period random input RfL=RV3LCalculating;
s1854, calling EQN400 sub-flow to obtain FV3L=FV3、FV3=N F0Time constant T1=2.0LW/VTRAnd gain K = ηW∙2.5133(LW/VTR)2Wherein ηWIndicating the intensity of the turbulence in the w direction, LWRepresents the scale of the turbulence in the w direction;
s1855, according to EXP = e-∆t/T1Filter output CF0=FW1And a random number Rf=RW1Calculating;
s1856, calling EQN200 sub-flow to obtain FW1L=FW1、FW1=N F0Time constant T1=2.0LW/VTRAnd a gain K =0.15915, wherein FWIL(I =1,2, 3) represents the upper-cycle vertical filter output;
s1857, according to EXP = e-∆t/T1Filter output value CF0=FW2Upper period filter output value DF0=FW2LRandom input Rf=RW2And an upper period random input RfL=RW2LCalculating;
s1858, calling EQN300 subflow to obtain FW2L=FW2、FW2=N F0Time constant 1T1=LW/VTRGain K = ηWLW/2VTRAnd time constant 2T2=√3∙T1;
S1859, according to EXP = e-∆t/T1Filter output value CF0=FW3Upper period filter output DF0=FW3LA random number Rf=RW3And an upper periodic random number RfL=RW3LCalculating;
s1860, call EQN400 sub-flow to get FW3L=FW3And FW3=N F0;
S1861, quitting the filter subprocess.
Specifically, the EQN200 sub-process includes the following steps:
proceed to EQN200 subroutine;
calculate NF0=CF0 e-∆t/T1+RfK(1- e-∆t/T1);
Judgment | NF0If | is less than 0.00001, if yes, NF is obtained0=0, if not, executing the next step;
the subroutine EQN200 is exited.
EQN300 the subflow includes the following steps:
proceed to EQN300 subroutine;
calculate NF0=(2CF0-DF0∙EXP)∙EXP+∆t∙[K∙EXP(Rf-RfL)/T1 2];
Judgment | NF0If | is less than 0.00001, if yes, NF is obtained0=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 F1=(2.0∙CF0-DF0∙EXP)∙EXP,F21=KT2/T1 2,F22= T1 2/T2,F23=(∆tT1/T2) -∆t,F31=Rf∙F21∙[F22+(-F23-F22)∙EXP],F32=RfL∙F21∙[F22(EXP-1.0)+F23]∙EXP,NF0=F1+F31+F32;
Judgment | NF0If | is less than 0.00001, if yes, NF is obtained0=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 etaTHL=1.0, and execute S1878, if not, execute S1873;
s1873 and judging the flying height hRWhether or not less than the turbulent low-height limit hL-0.5, if yes, perform S1875, if no, perform S1874;
s1874 and judging hRWhether or not greater than the turbulent high height limit hU+0.5, if yes, execute S1875, if no, execute S1876;
s1875 obtaining a turbulence factor etaTHL=0, and perform S1878;
s1876 and judging hRWhether or not h is less than or equal toLIf yes, get ηTHL=1.0-(hL-hR) ∙ 2.0.0, and executing S1878, if not, executing S1877;
s1877 and judging hRWhether or not h is greater than or equal toUIf yes, get ηTHL=1.0+(hU-hR) ∙ 2.0.0, and executing S1878, if not, obtaining etaTHL=1.0, and perform S1878;
s1878, calculating U in sequenceRAL=FU1∙FU2+FU3、URAL=FV1∙FV2+FV3、wRAL= FW1∙FW2+FW3∙WJETUP+WCOB、PRAL=0.025l(WRAL-WJETUP)/LW、QRAL=0.025(WRAL-WJETUP)81.5/LW、RRAL=0.025∙81.5VRAL/LW、UT=URAL∙ηTHL、VT=VRAL∙ηTHL、WT=WRAL∙ηTHL、PT=PRAL∙ηTHL、QT=QRAL∙ηTHL、RT=RRAL∙ηTHL、、、、UTOL=UT、VTOL=VTAnd WTOL=WTWherein P isTRepresenting the DRYDEN model output pitch angular velocity, QTRepresenting DRYDEN model output yaw rate, RTRepresenting the DRYDEN model output roll angular velocity, UTLRepresents the X-direction speed, V, of the DRYDEN model output in the upper periodTLRepresents the Y-direction speed, W of the DRYDEN model output in the upper periodTLRepresents 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 airplaneREPIf 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 etaTU= η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 isTUIs 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 etaTUIf 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 TJETUPWhether 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 airplaneREPIf 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 etaTU= η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 isTUIs 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 gustWherein, IAXTURIN represents the turbulence intensity, if not, the turbulence intensity factor eta of the vertical gust is obtainedJETLOC=0;
S1313, judging aircraft gust timing TJETUPIf less than 30.0, if yes, T is obtainedJETUP=TJETUPΔ 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 TJETUPWhether 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 interpolationJUM=f(h1) Wherein h is1Represents a height;
s1318, obtaining vertical gust speed WJETUP=f(h2)∙ηJETLOC∙ηJUMWherein h is2Represents 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,URA) And a second fade-out sub-flow (V)RAL,VRA) Wherein, URALRepresents the longitudinal component factor, U, of the turbulent wind speed under the machine systemRARepresenting the longitudinal component, V, of the turbulent wind velocity under the machine systemRALRepresents the lateral component factor, V, of the turbulent wind speed under the machine systemRARepresenting the lateral component of the turbulent wind speed under the machine system;
s143, calling a third fade-out sub-flow (W)RAL,WRA) And a fourth fade-out sub-flow (P)RAL,PRA) Wherein W isRALRepresents the vertical component factor, W, of the turbulent wind speed under the machine systemRARepresents the vertical component, P, of the turbulent wind speed under the machine systemRALPitch angle velocity factor, P, representing turbulent wind velocity generationRARepresenting the pitch angle velocity resulting from turbulent wind speed;
s144, calling a fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA) Wherein Q isRALIndicating deviation of turbulent wind velocity generationAngular velocity factor, QRARepresenting yaw rate, R, resulting from turbulent wind speedRALRoll angular velocity factor, R, representing turbulent wind velocity generationRARepresenting the roll angular velocity generated by turbulent wind speed;
s145, obtaining etaTU= ηTUL= ηTUD=0;
S146, clearing all filter outputs FUI=0、FVI=0、FWI=0, I =1,2,3, wherein FUIRepresenting the longitudinal filter output, FVIRepresenting the lateral filter output, FWIRepresenting 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,URA) And a second decreasing subprocess (V)RAL,VRA) And the third fade-out subprocess (W)RAL,WRA) And the fourth fade-out sub-process (P)RAL,PRA) The fifth fade-out sub-process (Q)RAL,QRA) And a sixth fade-out sub-flow (R)RAL,RRA)。
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 URAL、VRAL、WRAL、PRAL、QRAL、RRALRaffader =0.95, a factor fixed value representing the calculation speed, RC represents URA、VRA、WRA、PRA、QRA、RRA;
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 atmosphereAnd executing S1806, if not, executing S1804;
s1804, judging whether the turbulence type is cobblestone turbulence, if so, obtaining cobblestone atmosphere required turbulence intensity factorsAnd executing S1805, if not, obtaining etaTUD=ηCOBLOCAnd 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 etaTUD=0, and executing S1808, if not, executing S1808;
s1808 obtaining the increment of the turbulent intensity factor [. eta. ]TU=ηTUD-ηTULEqual and equal toTUIf 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 sequenceTU=η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 cobblestonesCOB=0、ACOB=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 TCS=TCSAnd judging the vacuum speed V of the airplaneTRIf it is greater than 700.0, get ηTU=0.01, and executes S1818, and if not, executes S1817;
s1817, judging vacuum speed VTRIf it is less than 200.0, if so, the rough atmosphere timing T is obtainedP=0.333333, and execute the next step, if not, get TP=1/[0.006(VTR200.0) +3.0 ], and the next step is performed;
s1818, judging TCSWhether or not greater than TPIf yes, then T is obtainedCS=0, and executing the next step, if not, executing the next step;
s1819, judging TCSWhether or not less than TPIf so, obtaining WCOB=ACOBIf not, obtaining WCOB=-ACOB;
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 RU1~RU3、RV1~RV3、RW1~RW3Wherein R isU1~RU3Denotes a random number, R, normally distributed to UV1~RV3Denotes a V-normally distributed random number, RW1~RW3Representing 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 T1=2.0LU/VTRAnd gain K =2.8264 ηULU/VTRPerforming a calculation, whereinUIndicates the intensity of the turbulence in the u direction, LURepresents the scale of the turbulence in the u direction;
s1843, according to EXP = e-∆t/T1Initial value CF of filter0=FU1And a random input variable Rf=RU1Calculating;
s1844, calling EQN200 sub-flow to obtain FU1L=FU1、FU1=N F0Gain K =1.0, wherein,FUIL(I =1,2, 3) represents the up-cycle vertical filter output, NF0Represents the filter output;
s1845, according to EXP = e-∆t/T1Initial value CF of filter0=FU2And a random input variable Rf=RU2Calculating;
S1847, according to EXP = e-∆t/T1、F0=FU3And Rf=RU3Calculating;
s1848, calling EQN200 sub-flow to obtain FU3L=FU3、FU3=N F0Time constant T1=2.0LV/VTRGain K = ηU∙2.5133(LV/VTR)2Wherein L isVRepresents the scale of turbulence in the v-direction;
s1849, according to EXP = e-∆t/T1Filter output value CF0=FV1And a random input variable Rf=RV1Calculating;
s1850, calling EQN200 sub-flow to obtain FV10L=FV1、FV1=N F0Time constant T1=LV/VTRAnd a gain K =0.15915, wherein FVIL(I =1,2, 3) represents the upper-period lateral filter output;
s1851, according to EXP = e-∆t/T1Filter output value CF0=FV2Upper period filter output value DF0=FV2LRandom input Rf=RV1And an upper period random input RfL=RV2LCalculating;
s1852, calling EQN300 subflow to obtain FV2L=FV2、FV2=N F0Time constant 1T1=LV/VTRGain K = ηVLV/2VTRAnd time constant 2T2= 3 ∙ T1, wherein ηVRepresenting the turbulence intensity of the turbulence in the v direction;
s1853, according to EXP = e-∆t/T1Filter output value CF0=FV3Upper period filter output value DF0=FV3LRandom input Rf=RV3And an upper period random input RfL=RV3LCalculating;
s1854, calling EQN400 sub-flow to obtain FV3L=FV3、FV3=N F0Time constant T1=2.0LW/VTRAnd gain K = ηW∙2.5133(LW/VTR)2Wherein ηWIndicating the intensity of the turbulence in the w direction, LWRepresents the scale of the turbulence in the w direction;
s1855, according to EXP = e-∆t/T1Filter output CF0=FW1And a random number Rf=RW1Calculating;
s1856, calling EQN200 sub-flow to obtain FW1L=FW1、FW1=N F0Time constant T1=2.0LW/VTRAnd a gain K =0.15915, wherein FWIL(I =1,2, 3) represents the upper-cycle vertical filter output;
s1857, according to EXP = e-∆t/T1Filter output value CF0=FW2Upper period filter output value DF0=FW2LRandom input Rf=RW2And an upper period random input RfL=RW2LCalculating;
s1858, calling EQN300 subflow to obtain FW2L=FW2、FW2=N F0Time constant 1T1=LW/VTRGain K = ηWLW/2VTRAnd time constant 2T2=√3∙T1;
S1859, according to EXP = e-∆t/T1Filter output value CF0=FW3Upper period filter output DF0=FW3LA random number Rf=RW3And an upper periodic random number RfL=RW3LCalculating;
s1860, call EQN400 sub-flow to get FW3L=FW3And FW3=N F0;
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 etaTHL=1.0, and execute S1878, if not, execute S1873;
s1873 and judging the flying height hRWhether or not less than the turbulent low-height limit hL-0.5, if yes, perform S1875, if no, perform S1874;
s1874 and judging hRWhether or not it is greater than hU+0.5, if yes, execute S1875, if no, execute S1876;
s1875 obtaining a turbulence factor etaTHL=0, and perform S1878;
s1876 and judging hRWhether or not h is less than or equal toLIf yes, get ηTHL=1.0-(hL-hR) ∙ 2.0.0, and executing S1878, if not, executing S1877;
s1877 and judging hRWhether or not the height of the turbulent flow is greater than or equal to the limit hUIf yes, get ηTHL=1.0+(hU-hR) ∙ 2.0.0, and executing S1878, if not, obtaining etaTHL=1.0, and perform S1878;
s1878, calculating U in sequenceRAL=FU1∙FU2+FU3、URAL=FV1∙FV2+FV3、wRAL= FW1∙FW2+FW3∙WJETUP+WCOB、PRAL=0.025l(WRAL-WJETUP)/LW、QRAL=0.025(WRAL-WJETUP)81.5/LW、RRAL=0.025∙81.5VRAL/LW、UT=URAL∙ηTHL、VT=VRAL∙ηTHL、WT=WRAL∙ηTHL、PT=PRAL∙ηTHL、QT=QRAL∙ηTHL、RT=RRAL∙ηTHL、、、、UTOL=UT、VTOL=VTAnd WTOL=WTWherein P isTRepresenting the DRYDEN model output pitch angular velocity, QTRepresenting DRYDEN model output yaw rate, RTRepresenting the DRYDEN model output roll angular velocity, UTLRepresents the X-direction speed, V, of the DRYDEN model output in the upper periodTLRepresents the Y-direction speed, W of the DRYDEN model output in the upper periodTLRepresents 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.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110568308.6A CN113239462A (en) | 2021-05-25 | 2021-05-25 | Simulation method for aircraft turbulent environment simulation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110568308.6A CN113239462A (en) | 2021-05-25 | 2021-05-25 | Simulation method for aircraft turbulent environment simulation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113239462A true CN113239462A (en) | 2021-08-10 |
Family
ID=77138649
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110568308.6A Pending CN113239462A (en) | 2021-05-25 | 2021-05-25 | Simulation method for aircraft turbulent environment simulation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113239462A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116337396A (en) * | 2023-05-30 | 2023-06-27 | 中国航空工业集团公司哈尔滨空气动力研究所 | Method for actively simulating wind tunnel test by using high-altitude atmospheric turbulence |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5860807A (en) * | 1996-04-22 | 1999-01-19 | The United States Of America As Represented By The National Aeronautics And Space Administration | System and method for finite element simulation of helicopter turbulence |
CN101650883A (en) * | 2009-02-13 | 2010-02-17 | 中国人民解放军空军航空大学 | Simulation method of atmospheric turbulence on flight simulator |
CN105354387A (en) * | 2015-11-26 | 2016-02-24 | 哈尔滨工业大学 | Improved gaussian white noise based three-dimensional atmospheric turbulence generating method |
CN205175659U (en) * | 2015-12-02 | 2016-04-20 | 同济大学 | A wind -tunnel turbulent flow analogue means |
CN106874535A (en) * | 2016-12-28 | 2017-06-20 | 南京航空航天大学 | A kind of method for designing and emulation mode of high-order air turbulence formed filter |
CN107590353A (en) * | 2017-10-19 | 2018-01-16 | 郑州云海信息技术有限公司 | The cluster of the server of air turbulence field stimulation method and use KNL processors |
CN107729638A (en) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | Anisotropy In The Atmospheric Turbulent Field method for numerical simulation |
US20180232473A1 (en) * | 2017-02-14 | 2018-08-16 | Corvid Technologies LLC | Methods, systems, and computer program products for determining the effect of turbulence on an entity |
CN111859529A (en) * | 2020-06-11 | 2020-10-30 | 北京航空航天大学 | Multi-grid disturbance domain updating acceleration method for aircraft streaming numerical simulation |
CN112800639A (en) * | 2021-04-13 | 2021-05-14 | 江苏普旭科技股份有限公司 | Simulation method for aircraft wind environment simulation |
-
2021
- 2021-05-25 CN CN202110568308.6A patent/CN113239462A/en active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5860807A (en) * | 1996-04-22 | 1999-01-19 | The United States Of America As Represented By The National Aeronautics And Space Administration | System and method for finite element simulation of helicopter turbulence |
CN101650883A (en) * | 2009-02-13 | 2010-02-17 | 中国人民解放军空军航空大学 | Simulation method of atmospheric turbulence on flight simulator |
CN105354387A (en) * | 2015-11-26 | 2016-02-24 | 哈尔滨工业大学 | Improved gaussian white noise based three-dimensional atmospheric turbulence generating method |
CN205175659U (en) * | 2015-12-02 | 2016-04-20 | 同济大学 | A wind -tunnel turbulent flow analogue means |
CN106874535A (en) * | 2016-12-28 | 2017-06-20 | 南京航空航天大学 | A kind of method for designing and emulation mode of high-order air turbulence formed filter |
US20180232473A1 (en) * | 2017-02-14 | 2018-08-16 | Corvid Technologies LLC | Methods, systems, and computer program products for determining the effect of turbulence on an entity |
CN107729638A (en) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | Anisotropy In The Atmospheric Turbulent Field method for numerical simulation |
CN107590353A (en) * | 2017-10-19 | 2018-01-16 | 郑州云海信息技术有限公司 | The cluster of the server of air turbulence field stimulation method and use KNL processors |
CN111859529A (en) * | 2020-06-11 | 2020-10-30 | 北京航空航天大学 | Multi-grid disturbance domain updating acceleration method for aircraft streaming numerical simulation |
CN112800639A (en) * | 2021-04-13 | 2021-05-14 | 江苏普旭科技股份有限公司 | Simulation method for aircraft wind environment simulation |
Non-Patent Citations (1)
Title |
---|
王化东: "大气紊流建模与仿真", 飞行力学, vol. 32, no. 01, 28 February 2014 (2014-02-28) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116337396A (en) * | 2023-05-30 | 2023-06-27 | 中国航空工业集团公司哈尔滨空气动力研究所 | Method for actively simulating wind tunnel test by using high-altitude atmospheric turbulence |
CN116337396B (en) * | 2023-05-30 | 2023-07-21 | 中国航空工业集团公司哈尔滨空气动力研究所 | Method for actively simulating wind tunnel test by using high-altitude atmospheric turbulence |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Robert | A semi-Lagrangian and semi-implicit numerical integration scheme for the primitive meteorological equations | |
Menut et al. | CHIMERE 2013: a model for regional atmospheric composition modelling | |
Vogel et al. | Influence of topography and biogenic volatile organic compounds emission in the state of Baden‐Württemberg on ozone concentrations during episodes of high air temperatures | |
Lorenz | Maximum simplification of the dynamic equations | |
Schoeberl et al. | The zonally averaged circulation of the middle atmosphere | |
Bao et al. | Large-eddy simulation over complex terrain using an improved immersed boundary method in the Weather Research and Forecasting Model | |
Krishnamurti et al. | Downstream amplification and formation of monsoon disturbances | |
CN113239462A (en) | Simulation method for aircraft turbulent environment simulation | |
Berry et al. | An analysis of cloud drop growth by collection: Part III. Accretion and self-collection | |
Wang et al. | Top-hat representation of turbulence statistics in cloud-topped boundary layers: A large eddy simulation study | |
Somerville | The forgiving air: understanding environmental change | |
CN104794311A (en) | Global optimal double-pulse orbit transfer method under time free conditions | |
Petry et al. | Chemical conversion of subsonic aircraft emissions in the dispersing plume: Calculation of effective emission indices | |
Byrd | A Computer Model of the Origin of H i Radial Velocity Asymmetries in M31 | |
Granier et al. | Ozone and other trace gases in the Arctic and Antarctic regions: Three‐dimensional model simulations | |
CN110807270B (en) | Method for inverting engine parameters and estimating trajectory based on tail flame radiation line type | |
Lamraoui et al. | WRF hindcasts of cold front passages over the ARM eastern North Atlantic site: A sensitivity study | |
Kärcher | Transport of exhaust products in the near trail of a jet engine under atmospheric conditions | |
Ruzmaikin | Climate as a game of chance | |
Danilin et al. | Evolution of the concentrations of trace species in an aircraft plume: Trajectory study | |
CN108647419B (en) | Low-altitude flame infrared radiation characteristic estimation method and device along with height change | |
Nuissier et al. | A numerical simulation of Hurricane Bret on 22–23 August 1999 initialized with airborne Doppler radar and dropsonde data | |
CN102915391A (en) | Method for correcting numerical simulation of atmospheric turbulences | |
Pickersgill et al. | An examination of the formation of linear lee waves generated by giant Martian volcanoes | |
Kanari | On the study of numerical experiments of two layer Lake Biwa |
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 |