A method for designing atmospheric wind fields for aircraft simulation

Abstract

The present application belongs to the field of aircraft simulation technology, and particularly relates to a method for designing atmospheric wind field for aircraft simulation. The method comprises the following steps: step 1: firstly, determining the longitude and latitude range of the region to be flown by the flight simulator, dividing the region into unit regions according to the longitude and latitude coordinates, and then arranging the regions according to the longitude and latitude; step 2: setting the initial value of the wind speed in each region; step 3: calling the terrain data to calculate the wind force in different height layers within the region; and step 4: making the wind field data into a wind field function related to the longitude, latitude and height. The present application has the beneficial effect that the wind field structure in the aircraft simulator is more finely subdivided and conforms to the actual weather conditions.

Description

Technical Field

[0001] This invention belongs to the field of aircraft simulation technology, specifically relating to an atmospheric wind field design method for aircraft simulation. Background Technology

[0002] In the field of aircraft simulation, aircraft stability is an important research topic. Atmospheric wind field is an important environmental condition for flight simulation. A changing wind field will affect the simulated control of the aircraft. Currently, the wind field environmental conditions of flight simulators are all uniform single wind field, that is, the wind speed is set on the simulator to create a uniform wind field, or a few random atmospheric turbulences are interspersed in it. These practices cannot well meet the flight experience of flight simulators under various extreme conditions.

[0003] Currently, wind field settings in aircraft simulators use preset values ​​from the training console. This involves setting the wind direction (eight options: East, South, West, North, Southeast, Northeast, Southwest, and Northwest) on the console, followed by setting the wind strength (light breeze, strong wind, gusts, etc.) and finally the wind speed in m / s. While this method can initially meet the wind requirements of aircraft flight, it lacks precision. In real-world scenarios, wind fields are highly variable, so aircraft simulators need to consider more diverse design methods to make the simulation more realistic.

[0004] After using the previous wind field setting method, the wind field in the aircraft simulator will be as shown in the following figure:

[0005] Figure 1 While a single wind direction field can provide wind field data, the wind is monotonous and cannot change near the ground. It cannot adequately meet the needs of special wind transitions, updrafts, downdrafts, etc. Furthermore, the wind force will not always be a single value over the range of tens to hundreds of kilometers in which an aircraft flies.

[0006] Figure 2 It is a large-area meteorological distribution map of wind speed in a certain region. It can be seen that the wind field is distributed very irregularly. At the same time, in mountainous areas, forest areas, and near high mountains, local wind shear will occur. All these situations need to be simulated in the flight simulator. Summary of the Invention

[0007] The purpose of this invention is to provide a method for designing atmospheric wind fields for aircraft simulation, which can simulate wind field systems under various extreme conditions and can be connected to real-time meteorological wind field systems in real-time.

[0008] The technical solution of the present invention is as follows: A method for designing atmospheric wind fields for aircraft simulation, comprising the following steps:

[0009] Step 1: First, determine the latitude and longitude range of the area where the flight simulator will fly, divide the area into unit areas according to latitude and longitude coordinates, and then arrange the units according to latitude and longitude.

[0010] Step 2: Set the initial wind speed value for each area;

[0011] In step 2, the initial wind speed value is preset. The wind speed area with the largest proportion of the same wind speed in the area is set as the global wind speed. Then, values ​​are assigned to the internal small areas by assigning values ​​to elliptical areas, circular areas, and square areas.

[0012] The initial wind speed value for step 2 is preset as follows:

[0013] J 经度下限 J represents the lower limit of the regional longitude, where J is the longitude of the aircraft. 经度上限 W is the upper limit of the region's longitude. 纬度下限 W represents the lower limit of the latitude of the region, where W is the latitude of the aircraft. 纬度上限 J represents the upper limit of the region's latitude. 经度下限 <J<J 经度上限 And W 纬度下限 <W<W 纬度上限 ,but:

[0014] The regional wind field parameters assigned in the NED coordinate system are as follows:

[0015] East axis Xned=x, North axis Yned=y, Earth axis Zned=z

[0016] If the secondary region is a circular region, specify the latitude and longitude of the region center (J). 圆形 W 圆形 ), where J 圆形 W is the longitude of the center of the circle. 圆形 The radius of the circle's center is L kilometers.

[0017] The radius of a circle is converted to latitude and longitude as R = 1000L / D°, where D is the latitude and longitude interval corresponding to 1 kilometer.

[0018] When (J) 圆形 -(J 区域min +0.00005) 2 +(W) 圆形 -(W) 区域min +0.00005) 2 <R 2 At that time, J 区域min W is the lower limit of the longitude of the regional interval. 区域min This represents the lower limit of the latitude of the regional interval;

[0019] The regional wind field parameters assigned in the NED coordinate system are as follows:

[0020] East axis Xned=x, north axis Yned=y, ground axis Zned=z

[0021] If the secondary region is a square region, the minimum value J of the specified region longitude is 方形min The maximum value J 方形max , the minimum value W 方形min of latitude 方形max The maximum value W 方形min ;

[0022] When J 区域min +0.00005<J 方形max , and W 方形min +W 区域min +0.00005<W 方形max

[0023] The NED coordinate system region wind field parameter assignment is respectively:

[0024] East axis Xned=x, north axis Yned=y, ground axis Zned=z

[0025] Step 3: Call topographic data, calculate the wind force in different height layers inside the region;

[0026] Step 4: The wind field data is made into a wind field function related to longitude, latitude and height.

[0027] The initial wind speed value in step 2 is the real-time meteorological data transmission, and the real-time meteorological data is directly assigned to the wind data of the corresponding longitude and latitude region according to the longitude and latitude region.

[0028] The step 3 includes the following:

[0029] The height division is 10 meters apart;

[0030] The geographical terrain flat region:

[0031] The NED coordinate system region wind field parameter assignment is respectively: J is the longitude, w is the latitude, and h is the height above the ground;

[0032] East axis X(J, w, h)=x, x is the east axis wind speed of the flat area;

[0033] North axis Y(J, w, h)=y, y is the north axis wind speed of the flat area;

[0034] Ground axis Z(J, w, h)=z, z is the ground axis wind speed of the flat area;

[0035] There is a windward slope in the region, l is the highest point of the slope, and g is the lowest point of the slope.

[0036] In step 3, h>l+20

[0037] East axis X (J, w, h) = x, x is the east axis wind speed of the windward slope area;

[0038] North axis Y (J, w, h) = y, y is the north axis wind speed of the windward slope area;

[0039] Ground axis Z (J, w, h) = z, z is the ground axis wind speed of the windward slope area.

[0040] In step 3, l < h < l + 20

[0041] East axis X (J, w, h) = x, x is the east axis wind speed of the windward slope area;

[0042] North axis Y (J, w, h) = y, y is the north axis wind speed of the windward slope area;

[0043] Ground axis Z (J, w, h) = z * (h-l) / 20.0, z is the ground axis wind speed of the windward slope area.

[0044] In step 3, g < h < l

[0045] East axis X (J, w, h) = 0, x is the east axis wind speed of the windward slope area;

[0046] North axis Y (J, w, h) = y + x, y is the north axis wind speed of the windward slope area;

[0047] Ground axis Z (J, w, h) = 0.

[0048] In step 3, there is a leeward slope in the region, l is the highest point of the slope, and g is the lowest point of the slope:

[0049] h > l + 20;

[0050] East axis X (J, w, h) = x, x is the east axis wind speed of the leeward slope area;

[0051] North axis Y (J, w, h) = y, y is the north axis wind speed of the leeward slope area;

[0052] Ground axis Z (J, w, h) = z, z is the ground axis wind speed of the leeward slope area.

[0053] In step 3, l < h < l + 20

[0054] X (J, w, h) = x, x is the east axis wind speed of the leeward slope area;

[0055] Y (J, w, h) = y, y is the north axis wind speed of the leeward slope area;

[0056] Z (J, w, h) = z * (h-l) / 20.0, z is the ground axis wind speed of the leeward slope area.

[0057] In g < h < l

[0058] X(J, w, h) = 0;

[0059] Y(J, w, h) = 0;

[0060] Z(J, w, h) = 0.

[0061] The beneficial effect of the present application is that the wind field structure in the aircraft simulator will be more subdivided to fit the real weather conditions after using the method. BRIEF DESCRIPTION OF DRAWINGS

[0062] Figure 1 To set a single wind field situation;

[0063] Figure 2 To a large area of the wind speed of the region weather distribution map;

[0064] Figure 3 To the wind force change graph after the region is divided in the height direction;

[0065] Figure 4 To the effect diagram after the grid wind field;

[0066] Figure 5 To the wind field logic diagram. DETAILED DESCRIPTION

[0067] The present application will be further described in detail below in combination with the drawings and specific embodiments.

[0068] The atmospheric wind field design method for aircraft simulation provided by the present application introduces real-time environmental wind field measurement data and expected setting wind field parameters through setting two input ends, to realize wind field setting input in the aircraft simulator under multiple situations.

[0069] For the existing wind field setting situation, the present application associates the regional wind field division with the geographical latitude and longitude, and divides the wind field into unit blocks in the navigation coordinate system (NED coordinate system). The aircraft only exists wind force change in the height coordinate within the unit block. The unit block can be taken as a region according to the experience of wind field shear, which is 20m to 100m. The initial wind force and wind direction of the wind field can be input according to the detection station data provided by the meteorological department or the data set by the operator according to the needs. When there are mountains and other obstacles in the region, a simple wind force change curve in the height direction is fitted according to the call. Finally, the wind speed in the simulator will be similar to that of Figure 2 The local wind force and wind direction change caused by the change of terrain or weather can be displayed in the local region.

[0070] The wind field structure in the aircraft simulator will be more subdivided to fit the real weather conditions after using the method of the present application, such asFigure 4 The effect diagram after using the method is shown.

[0071] Reference Figure 2 The large-area wind speed weather distribution map in the region is displayed by dividing the region and then displaying the wind speed parameters in the NED coordinate. Then, the wind speed and wind force data in the cell are converted into a function of the processing wind field and latitude and longitude, and the function format is:

[0072] If the cell longitude lower limit < longitude data < cell longitude upper limit, and the cell dimension lower limit < dimension data < cell dimension upper limit, then:

[0073] The NED coordinate system region wind field parameters are respectively:

[0074] The east axis Xned parameter value is:

[0075] The north axis Yned parameter value is:

[0076] The ground axis Zned parameter value is:

[0077] After processing, a wind field function related to latitude, longitude and height is generated.

[0078] The console only needs to directly call the function, and the real-time data of the weather station can be directly imported into the processing program to update the function in real time. The application will be further described in combination with the drawings and embodiments, and the drawings show:

[0079] Figure 5 A logic diagram for generating a new wind field of the technical method is shown, which shows the logic of the processing and calling of the wind field.

[0080] A method for designing an atmospheric wind field for aircraft simulation, comprising the following steps:

[0081] Step 1: First, determine the latitude and longitude range of the region to be flown by the flight simulator, divide the region into unit regions with (0.0001 seconds, 0.0001 seconds) as the unit according to the latitude and longitude coordinates, and then arrange the regions in order from small to large according to the latitude and longitude.

[0082] Step 1: Set the initial wind speed value in each region. The initial wind speed value can be obtained from two sources, one is user pre-setting, such as Figure 2 The user can first set the wind speed region with the largest proportion of the same wind speed in the region as the global wind speed, and then use the assignment method in the elliptical region, circular region and square region to assign values to the small regions inside to cover the wind speed of the small regions with not very large proportion of the region, as follows:

[0083] J 经度下限 J is the longitude of the aircraft, and J 经度上限 is the longitude upper limit of the region, W 维度下限W is the lower bound of the region dimension, and W is the dimension in which the spacecraft is located. 维度上限 J is the upper limit of the region dimension. 经度下限 <J<J 经度上限 And W 维度下限 <W<W 维度上限 ,but:

[0084] The regional wind field parameters assigned in the NED coordinate system are as follows:

[0085] East axis Xned=x, North axis Yned=y, Earth axis Zned=z

[0086] If the secondary region is circular, the user needs to specify the latitude and longitude of the region center (J). 圆形 W 圆形 ), where J 圆形 W is the longitude of the center of the circle. 圆形 The radius of the circle is L kilometers.

[0087] The radius of a circle can be converted to latitude and longitude as R = 1000L / D°, where D is the latitude and longitude interval corresponding to 1 kilometer.

[0088] When (J) 圆形 -(J 区域min +0.00005) 2 +(W) 圆形 -(W) 区域min +0.00005) 2 <R 2 At that time, J 区域min W is the lower limit of the longitude of the regional interval. 区域min This represents the lower limit of the latitude of the regional interval;

[0089] The regional wind field parameters assigned in the NED coordinate system are as follows:

[0090] East axis Xned = x, North axis Yned = y, Earth axis Zned = z

[0091] If the secondary area is a rectangular area, the user needs to specify the minimum longitude J of the area. 方形min With the maximum value J 方形max minimum latitude W 方形min With the maximum value W 方形max .

[0092] When J 方形min <J 区域min +0.00005 <J 方形max And W 方形min <W 区域min +0.00005 <W 方形max

[0093] The NED coordinate system under the regional wind field parameter assignment is respectively:

[0094] The east axis Xned=x, the north axis Yned=y, and the ground axis Zned=z

[0095] The x, y, and z have the same meaning, x represents the wind speed value of the east axis of the region in the NED coordinate system, y represents the wind speed value of the north axis of the region in the NED coordinate system, and z represents the wind speed value of the ground axis of the region in the NED coordinate system. The three values are manually set by the aircraft instructor according to the actual situation or training requirements.

[0096] Second, the real-time meteorological data is sent, and the real-time meteorological data is directly assigned to the wind power data of the corresponding latitude and longitude region according to the latitude and longitude region.

[0097] Step 3: Call the terrain data to calculate the wind power in different height layers in the region.

[0098] The height is divided by 10 meters.

[0099] The geographical terrain is flat:

[0100] The NED coordinate system under the regional wind field parameter assignment is respectively: where J is the longitude, w is the dimension, and h is the space height above the ground.

[0101] East axis X(J, w, h)=x, x is the east axis wind speed of the flat area

[0102] North axis Y(J, w, h)=y, y is the north axis wind speed of the flat area

[0103] Ground axis Z(J, w, h)=z, z is the ground axis wind speed of the flat area

[0104] There is a windward slope in the region, l is the highest point of the slope, and g is the lowest point of the slope:

[0105] h>l+20

[0106] East axis X(J, w, h)=x, x is the east axis wind speed of the windward slope area

[0107] North axis Y(J, w, h)=y, y is the north axis wind speed of the windward slope area

[0108] Ground axis Z(J, w, h)=z, z is the ground axis wind speed of the windward slope area

[0109] l<h<l+20

[0110] East axis X(J, w, h)=x, x is the east axis wind speed of the windward slope area

[0111] North axis Y(J, w, h)=y, y is the north axis wind speed of the windward slope area

[0112] Z(J, w, h) = z * (h - l) / 20.0, z is the ground axis wind speed of the windward slope area

[0113] In g < h < l

[0114] X(J, w, h) = 0

[0115] Y(J, w, h) = y + x, y is the north axis wind speed of the windward slope area

[0116] Z(J, w, h) = 0

[0117] There is a leeward slope in the region, l is the highest point of the slope, g is the lowest point of the slope:

[0118] In h > l + 20

[0119] X(J, w, h) = x, x is the east axis wind speed of the leeward slope area

[0120] Y(J, w, h) = y, y is the north axis wind speed of the leeward slope area

[0121] Z(J, w, h) = z, z is the ground axis wind speed of the leeward slope area

[0122] In l < h < l + 20

[0123] X(J, w, h) = x, x is the east axis wind speed of the leeward slope area

[0124] Y(J, w, h) = y, y is the north axis wind speed of the leeward slope area

[0125] Z(J, w, h) = z * (h - l) / 20.0, z is the ground axis wind speed of the leeward slope area

[0126] In g < h < l

[0127] X(J, w, h) = 0

[0128] Y(J, w, h) = 0

[0129] Z(J, w, h) = 0

[0130] Step 4: Wind field data is made into a wind field function related to latitude, longitude and height, which is called by the teaching console.

[0131] Among them, step 3 is the specific implementation process method of how to use the terrain data (latitude J, longitude w and height h) to make the wind field function (X, Y, Z) based on the wind field data (x, y, z). Finally, the wind field function is made.

Claims

1. A method for atmospheric wind field design in connection with aircraft simulation, characterized in that, It comprises the following steps: Step 1: Firstly, the longitude and latitude range of the flight area to be flown by the flight simulator is determined, the area is divided into unit areas according to the longitude and latitude coordinates, and then the area is arranged according to the longitude and latitude; Step 2: Set the initial wind speed value in each area; The initial wind speed value in step 2 is pre-set, the wind speed region with the largest proportion of the same wind speed in the region is set as the global wind speed, and then the internal small regions are assigned values in the form of elliptical region, circular region and square region; The pre-setting of the initial wind speed value in step 2 is as follows: J 经度下限 is the lower limit of the zone longitude, J is the longitude in which the aircraft is located, J 经度上限 is the upper limit of the zone longitude, W 纬度下限 is the lower limit of the zone latitude, W is the latitude in which the aircraft is located, W 纬度上限 is the upper limit of the zone latitude, J 经度下限 <J<J 经度上限 , and W 纬度下限 <W<W 纬度上限 then: The NED coordinate system under the regional wind field parameter assignment is respectively: East axis Xned=x, north axis Yned=y, and ground axis Zned=z If the secondary area is a circular area, specify the area center longitude and latitude (J 圆形 ,W 圆形 ), where J 圆形 is the circle center longitude, W 圆形 is the circle center latitude, and L is the radius length in kilometers. The radius of the circle is converted into the longitude and latitude length R=1000L / D°, and D is the longitude and latitude interval corresponding to a unit of 1 kilometer; When (J 圆形 - (J 区域min + 0.00005) < 0 2 + (W 圆形 - (W 区域min + 0.00005) < 0 2 < R 2 , J 区域min is the lower limit of the regional interval longitude, and W 区域min is the lower limit of the regional interval latitude; The NED coordinate system under the regional wind field parameter assignment is respectively: East axis Xned=x, north axis Yned=y, and ground axis Zned=z If the secondary area is a square area, the specified area longitude minimum value J 方形min and maximum value J 方形max , the latitude minimum value W 方形min and maximum value W 方形max ; When J 方形min < J 区域min + 0.00005 < J 方形max , and W 方形min < W 区域min + 0.00005 < W 方形max The NED coordinate system under the regional wind field parameter assignment is respectively: East axis Xned=x, north axis Yned=y, and ground axis Zned=z Step 3: Call the terrain data to calculate the wind force in different height layers in the region; Step 4: The wind field data is made into a wind field function related to longitude, latitude and height.

2. A method of atmospheric wind field design for aircraft simulation according to claim 1, wherein: The initial wind speed value in step 2 is real-time meteorological data, and the real-time meteorological data is directly assigned to the wind force data of the corresponding longitude and latitude region according to the longitude and latitude region.

3. A method for designing atmospheric wind fields for aircraft simulation as in claim 1, wherein, The step 3 comprises the following: The height is divided by 10 meters as an interval; The geographical terrain is flat: The NED coordinate system under the regional wind field parameter assignment is respectively: wherein J is the longitude, w is the latitude, and h is the space height above the ground; East axis X(J, w, h)=x, x is the east axis wind speed of the flat area; North axis Y(J, w, h)=y, y is the north axis wind speed of the flat area; Ground axis Z(J, w, h)=z, z is the ground axis wind speed of the flat area; There is a windward slope in the region, l is the highest point of the slope, and g is the lowest point of the slope.

4. A method of atmospheric wind field design for aircraft simulation according to claim 3, wherein: In step 3, h>l+20 East axis X(J, w, h)=x, x is the east axis wind speed of the windward slope area; North axis Y(J, w, h)=y, y is the north axis wind speed of the windward slope area; Ground axis Z(J, w, h)=z, z is the ground axis wind speed of the windward slope area.

5. A method of atmospheric wind field design for aircraft simulation according to claim 3, wherein: In step 3, l<h<l+20 East axis X(J, w, h)=x, x is the east axis wind speed of the windward slope area; North axis Y(J, w, h)=y, y is the north axis wind speed of the windward slope area; Ground axis Z(J, w, h)=z*(h-l) / 20.0, z is the ground axis wind speed of the windward slope area.

6. A method of atmospheric wind field design for aircraft simulation according to claim 3, wherein: In step 3, g<h<l East axis X(J, w, h)=0, x is the east axis wind speed of the windward slope area; North axis Y(J, w, h)=y+x, y is the north axis wind speed of the windward slope area; Ground axis Z(J, w, h)=0.

7. A method of atmospheric wind field design for aircraft simulation according to claim 3, wherein: In step 3, there is a leeward slope in the region, l is the highest point of the slope, and g is the lowest point of the slope: h>l+20; East axis X(J, w, h)=x, x is the east axis wind speed of the leeward slope area; North axis Y(J, w, h) = y, y is the north axis wind speed of the leeward slope area; Ground axis Z(J, w, h) = z, z is the ground axis wind speed of the leeward slope area; In l < h < l + 20 X(J, w, h) = x, x is the east axis wind speed of the leeward slope area; Y(J, w, h) = y, y is the north axis wind speed of the leeward slope area; Z(J, w, h) = z * (h - l) / 20.0, z is the ground axis wind speed of the leeward slope area; In g < h < l X(J, w, h) = 0; Y(J, w, h) = 0; Z(J, w, h) = 0.

Images