Atmospheric wind field design method for aircraft simulation

Abstract

The invention belongs to the technical field of aircraft simulation, and particularly relates to an atmospheric wind field design method for aircraft simulation. Comprising the following steps: step 1, firstly, determining a latitude and longitude range of a region to be flied by a flight simulator, dividing the region into unit regions according to latitude and longitude coordinates, and then arranging the regions according to latitude and longitude; 2, setting a regional wind speed initial value in each region; 3, calling landform data, and calculating wind power in different height layers in the region; and 4, making the wind field data into a wind field function related to longitude, latitude and height. The method has the beneficial effects that the wind field structure in the aircraft simulator can be more subdivided to fit the real weather condition after the method is adopted.

Description

Technical Field

[0001] The present invention belongs to the technical field of aircraft simulation, and in particular relates to an atmospheric wind field design method for aircraft simulation. Background Art

[0002] In the field of aircraft simulation, aircraft stability is an important research topic. The atmospheric wind field is a key 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 all use a unified single wind field, that is, the wind speed is set on the simulator to create a unified wind field, or several 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 teaching console. This means setting the wind direction on the console, with eight options: east, south, west, north, southeast, northeast, southwest, and northwest. Then, you set the wind force, with options like breeze, strong wind, or gust. Finally, you set the wind speed in meters per second. This setting method can initially meet the wind force requirements of aircraft flight, but it is not precise enough. In real-world scenarios, wind fields are highly variable, so aircraft simulators require more diverse design approaches to ensure that aircraft simulations are more realistic.

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

[0005] Figure 1 Although the single wind direction field in the image can provide wind field data, the wind direction is single and cannot change near the ground. It cannot meet the special wind switching, updraft, downdraft, etc. At the same time, the wind force will not always be a single value for the range of tens to hundreds of kilometers that the aircraft flies.

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

[0007] The purpose of the present invention is to provide an atmospheric wind field design method for aircraft simulation, which can simulate wind field systems of various extreme conditions and can be connected to a real-time meteorological wind field system in real time.

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

[0009] Step 1: First determine the latitude and longitude range of the flight simulator's flight area, divide the area into unit areas according to the latitude and longitude coordinates, and then arrange the areas according to the latitude and longitude;

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

[0011] Step 3: Call the terrain data to calculate the wind speed at different altitudes within the area;

[0012] Step 4: The wind field data is converted into a wind field function related to latitude, longitude and altitude.

[0013] In the step 1, the area is evenly divided into unit areas according to the latitude and longitude coordinates with (0.0001 seconds, 0.0001 seconds) as the unit, and the geographical length of the unit cell is D meters. Then the areas are arranged in ascending order according to the latitude and longitude.

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

[0015] The initial wind speed value in step 2 is preset as follows:

[0016] J 经度下限 is the lower limit of the regional longitude, J is the longitude of the aircraft, and J 经度上限 is the upper limit of regional longitude, W 维度下限 is the lower limit of the regional dimension, W is the dimension of the aircraft, and W 维度上限 is the upper limit of the regional dimension, J 经度下限 <J<J 经度上限 , and W 维度下限 <W<W 维度上限 ,but:

[0017] The regional wind field parameter assignments in the NED coordinate system are:

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

[0019] If the secondary area is a circular area, specify the longitude and latitude of the area center (J 圆形 ,W 圆形 ), where J 圆形 is the longitude of the center of the circle, W 圆形 The radius of the center is L kilometers;

[0020] The radius of the circle is converted into longitude and latitude length as R = 1000L / D, where D is the longitude and latitude interval corresponding to a unit of 1 kilometer;

[0021] When (J 圆形 -(J 区域min+0.00005)) 2 +(W 圆形 -(W 区域min +0.00005)) 2 <R 2 When J 区域min is the lower limit of the regional interval longitude, W 区域min is the lower limit of the regional interval latitude;

[0022] The regional wind field parameter assignments in the NED coordinate system are:

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

[0024] If the secondary area is a square area, specify the minimum longitude of the area J 方形min With the maximum value J 方形max , minimum latitude W 方形min With the maximum value W 方形max ;

[0025] When J 方形min <J 区域min +0.00005 <J 方形max , and W 方形min <W 区域min +0.00005 <W 方形max The regional wind field parameter assignments in the NED coordinate system are:

[0026] The east axis Xned=x, the north axis Yned=y, and the earth axis Zned=z.

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

[0028] The step 3 includes the following:

[0029] The height is divided into intervals of 10 meters;

[0030] Geographically flat areas:

[0031] The regional wind field parameter assignments in the NED coordinate system are: where J is the longitude, w is the latitude, and h is the spatial height above the ground;

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

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

[0034] Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in 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, when h>1+20

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

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

[0039] Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in the windward slope area.

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

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

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

[0043] Earth axis Z(J, w, h) = z*(hl) / 20.0, where z is the earth axis wind speed in the windward slope area.

[0044] In step 3, <h<l

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

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

[0047] The earth's 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: when h>l+20;

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

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

[0051] Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in the leeward slope area;

[0052] In l <h<l+20

[0053] X(J, w, h) = x, where x is the east-axis wind speed in the leeward slope area;

[0054] Y(J, w, h) = y, where y is the north-axis wind speed in the leeward slope area;

[0055] Z(J, w, h) = z*(hl) / 20.0, where z is the axial wind speed in the leeward slope area;

[0056] In g <h<l

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

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

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

[0060] The beneficial effect of the present invention is that after adopting the method, the wind field structure in the aircraft simulator will be more detailed and fit the real meteorological conditions. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0062] Figure 2 A large-area meteorological distribution map of wind speed in a certain area;

[0063] Figure 3 This is the wind force variation diagram in the height direction after regional division;

[0064] Figure 4 This is the effect diagram after gridding the wind field;

[0065] Figure 5 This is the wind farm logic diagram. DETAILED DESCRIPTION

[0066] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0067] The present invention provides an atmospheric wind field design method for aircraft simulation, which respectively introduces real-time environmental wind field measurement data and expected wind field parameters by setting two input terminals, thereby realizing wind field setting input in various situations in the aircraft simulator.

[0068] In view of the existing wind field settings, the present invention associates the regional wind field division with the geographical longitude and latitude, and uses the navigation coordinate system (NED coordinate system) to divide the wind field into unit blocks. The aircraft only experiences wind force changes under the altitude coordinates within the unit blocks. The unit blocks can be taken as an area of ​​20m to 100m according to the experience of wind field shear. The initial wind force and 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 encountering obstacles such as mountains in the area, a simple wind force change curve in the altitude direction is fitted according to the call. Finally, the wind speed in the simulator will be the same as Figure 2 In local areas, it can show changes in local wind force and direction caused by terrain or meteorological changes.

[0069] After adopting the method of the present invention, the wind field structure in the aircraft simulator will be more detailed and fit the real meteorological conditions, such as Figure 4 Shown is the effect diagram after using this method.

[0070] refer to Figure 2 The large-area wind speed meteorological distribution map in the image is divided into regions and then the wind speed parameters are displayed in NED coordinates. The wind speed and wind force data in the cell are then converted into a function that processes the wind field and longitude and latitude. The function format is:

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

[0072] The regional wind field parameters in the NED coordinate system are:

[0073] East axis Xned parameter value =

[0074] North axis Yned parameter value =

[0075] Earth axis Zned parameter value =

[0076] After processing, a wind field function related to latitude, longitude and altitude will be generated.

[0077] The teaching control console only needs to call the function directly, 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 present invention will be further described below with reference to the accompanying drawings and implementations. In the accompanying drawings:

[0078] Figure 5 A new wind farm logic diagram is generated for this technical method, which shows the logic of the wind farm's processing call.

[0079] A method for designing an atmospheric wind field for aircraft simulation includes the following steps:

[0080] Step 1: First determine the latitude and longitude range of the area where the flight simulator will fly, and divide the area into unit areas evenly according to the latitude and longitude coordinates with (0.0001 seconds, 0.0001 seconds) as the unit. The geographical length of the unit cell is D meters. Then arrange the areas in ascending order of latitude and longitude.

[0081] Step 1: Set the initial wind speed value in each area. The initial wind speed value can come from two sources. One is pre-set by the user, such as Figure 2 , the user can first set the wind speed area with the largest proportion of the same wind speed in the area as the global wind speed, and then use the elliptical area, circular area, and square area assignment method to assign values ​​to the internal small areas to cover the wind speed of small areas that are not very large. The details are as follows:

[0082] J 经度下限is the lower limit of the regional longitude, J is the longitude of the aircraft, and J 经度上限 is the upper limit of regional longitude, W 维度下限 is the lower limit of the regional dimension, W is the dimension of the aircraft, and W 维度上限 is the upper limit of the regional dimension, J 经度下限 <J<J 经度上限 , and W 维度下限 <W<W 维度上限 ,but:

[0083] The regional wind field parameter assignments in the NED coordinate system are:

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

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

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

[0087] When (J 圆形 -(J 区域min +0.00005)) 2 +(W 圆形 -(W 区域min +0.00005)) 2 <R 2 When J 区域min is the lower limit of the regional interval longitude, W 区域min is the lower limit of the regional interval latitude;

[0088] The regional wind field parameter assignments in the NED coordinate system are:

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

[0090] If the secondary area is a square 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 .

[0091] When J 方形min <J 区域min +0.00005 <J 方形max , and W 方形min <W 区域min +0.00005 <W方形max The regional wind field parameter assignments in the NED coordinate system are:

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

[0093] Here, x, y, and z have the same meaning: x represents the wind speed about the east axis of the region in the NED coordinate system, y represents the wind speed about the north axis of the region in the NED coordinate system, and z represents the wind speed about the earth axis of the region in the NED coordinate system. These three values ​​are manually set by the aircraft instructor based on actual conditions or training requirements.

[0094] The second is the sending of real-time meteorological data. The real-time meteorological data is directly assigned to the wind data of the corresponding longitude and latitude areas according to the longitude and latitude areas.

[0095] Step 3: Call the terrain data to calculate the wind force at different height layers within the area.

[0096] The height is divided into intervals of 10 meters.

[0097] Geographically flat areas:

[0098] The regional wind field parameter assignments in the NED coordinate system are: where J is the longitude, w is the latitude, and h is the spatial height above the ground.

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

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

[0101] Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in the flat area

[0102] 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:

[0103] When h>l+20

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

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

[0106] Earth axis Z(J, w, h) = z, z is the earth axis wind speed in the windward slope area

[0107] In l <h<l+20

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

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

[0110] Earth axis Z(J, w, h) = z*(hl) / 20.0, where z is the earth axis wind speed in the windward slope area

[0111] In g <h<l

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

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

[0114] Earth axis Z(J, w, h) = 0

[0115] 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:

[0116] When h>l+20

[0117] East axis X(J, w, h) = x, where x is the east axis wind speed in the leeward slope area

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

[0119] Earth axis Z(J, w, h) = z, z is the earth axis wind speed in the leeward slope area

[0120] In l <h<l+20

[0121] X(J, w, h) = x, where x is the east-axis wind speed in the leeward slope area

[0122] Y(J, w, h) = y, where y is the north-axis wind speed in the leeward slope area

[0123] Z(J, w, h) = z*(hl) / 20.0, where z is the axial wind speed in the leeward slope area

[0124] In g <h<l

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

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

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

[0128] Step 4: The wind field data is converted into a wind field function related to latitude, longitude and altitude for use by the teaching control console.

[0129] Step 3 is the specific execution process of how to use topographic data (latitude and longitude J, w and height h) to create a wind field function (X, Y, Z) based on the wind field data (x, y, z). The final result is the wind field function.

Claims

1. A method for designing an atmospheric wind field for aircraft simulation, characterized in that: The steps include: 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 the latitude and longitude coordinates, and then arrange the areas according to the latitude and longitude; Step 2: Set the initial value of regional wind speed in each area; Step 3: Call the terrain data to calculate the wind force at different height layers within the area; Step 4: The wind field data is made into a wind field function related to latitude, longitude and altitude.

2. The atmospheric wind field design method for aircraft simulation according to claim 1, characterized in that: In the step 1, the area is evenly divided into unit areas according to the longitude and latitude coordinates with (0.0001 seconds, 0.0001 seconds) as the unit, and the geographical length of the unit cell is D meters. Then the areas are arranged in the order of longitude and latitude from small to large.

3. The atmospheric wind field design method for aircraft simulation according to claim 1, characterized in that: In the step 2, the initial wind speed value is preset, and the wind speed area with the largest proportion of the same wind speed in the area is set as the global wind speed, and then the internal small area is assigned a value by assigning values ​​in the elliptical area, circular area, and square area.

4. The atmospheric wind field design method for aircraft simulation according to claim 3, characterized in that: The initial wind speed value of step 2 is preset as follows: J 经度下限 is the lower limit of the regional longitude, J is the longitude of the aircraft, and J 经度上限 is the upper limit of the regional longitude, W 维度下限 is the lower limit of the regional dimension, W is the dimension where the aircraft is located, and W 维度上限 is the upper limit of the regional dimension, J 经度下限 <J<J 经度上限 , and W 维度下限 <W<W 维度上限 ,but: The regional wind field parameter assignments in the NED coordinate system are: East axis Xned = x, North axis Yned = y, Earth axis Zned = z If the secondary area is a circular area, specify the longitude and latitude of the area center (J 圆形 ,W 圆形 ), where J 圆形 is the longitude of the center of the circle, W 圆形 The radius of the center is L kilometers; The radius of a circle is converted into longitude and latitude length as R = 1000L / D, where D is the longitude and latitude interval corresponding to 1 kilometer; When (J 圆形 -(J 区域min +0.00005)) 2 +(W 圆形 -(W 区域min +0.00005)) 2 <R 2 When J 区域min is the lower limit of the regional interval longitude, W 区域min is the lower limit of the regional interval latitude; The regional wind field parameter assignments in the NED coordinate system are: East axis Xned = x, North axis Yned = y, Earth axis Zned = z If the secondary region is a square region, specify the minimum longitude of the region J 方形min With the maximum value J 方形max , minimum latitude W 方形min With the maximum value W 方形max ; When J 方形min <J 区域min +0.00005 <J 方形max , and W 方形min <W 区域min +0.00005 <W 方形max The regional wind field parameter assignments in the NED coordinate system are: The east axis Xned=x, the north axis Yned=y, and the earth axis Zned=z.

5. The atmospheric wind field design method for aircraft simulation according to claim 1, characterized in that: The initial wind speed value in step 2 is sent by real-time meteorological data, and the real-time meteorological data is directly assigned to the wind data of the corresponding longitude and latitude area according to the longitude and latitude area.

6. The atmospheric wind field design method for aircraft simulation according to claim 1, characterized in that: The step 3 includes the following: The height is divided into intervals of 10 meters; Geographically flat areas: The regional wind field parameter assignments in the NED coordinate system are: where J is longitude, w is latitude, and h is the spatial height above the ground; East axis X(J, w, h) = x, where x is the east axis wind speed in the flat area; North axis Y(J, w, h) = y, y is the north axis wind speed in the flat area; Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in 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.

7. The atmospheric wind field design method for aircraft simulation according to claim 6, characterized in that: In step 3, when h>1+20 East axis X(J, w, h) = x, where x is the east axis wind speed in the windward slope area; North axis Y(J, w, h) = y, y is the north axis wind speed in the windward slope area; Earth axis Z(J, w, h) = z, where z is the earth axis wind speed in the windward slope area.

8. The atmospheric wind field design method for aircraft simulation according to claim 6, characterized in that: In step 3 described above, <h<l+20 East axis X(J, w, h) = x, where x is the east axis wind speed in the windward slope area; North axis Y(J, w, h) = y, y is the north axis wind speed in the windward slope area; Earth axis Z(J, w, h) = z*(hl) / 20.0, where z is the earth axis wind speed in the windward slope area.

9. The atmospheric wind field design method for aircraft simulation according to claim 6, characterized in that: In step 3 described above, <h<l East axis X(J, w, h) = 0, x is the east axis wind speed in the windward slope area; North axis Y(J, w, h) = y + x, y is the north axis wind speed in the windward slope area; The earth's axis Z(J, w, h)=0.

10. The atmospheric wind field design method for aircraft simulation according to claim 6, characterized in that: 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: At h>l+20; East axis X(J, w, h) = x, where x is the east axis wind speed in the leeward slope area; North axis Y(J, w, h) = y, y is the north axis wind speed in the leeward slope area; Earth axis Z(J, w, h) = z, z is the earth axis wind speed in the leeward slope area; In l <h<l+20 X(J, w, h) = x, where x is the east-axis wind speed in the leeward slope area; Y(J, w, h) = y, y is the north axis wind speed in the leeward slope area; Z(J, w, h) = z*(hl) / 20.0, z is the axial wind speed in the leeward slope area; In g <h<l X(J,w,h)=0; Y(J,w,h)=0; Z(J,w,h)=0.

Images