A method for calculating average wind load of a super-large aperture radio telescope reflector

Abstract

The application relates to a kind of average wind load calculation methods of super-large caliber radio telescope reflecting surface, for the antenna stiffness and wind speed difference of radio telescope top end and low end due to its large caliber, the radio telescope reflecting surface and wind speed gradient are zoned, the wind force and wind moment size of each reflecting surface partition are obtained, and then the total average wind load of reflecting surface is obtained by summation.This average wind load calculation method of the application, according to the characteristics of radio telescope reflecting surface with high symmetry and the law that wind increases with height, the radio telescope reflecting surface and wind speed gradient are zoned, the wind load of each partition is calculated respectively and then summed, the average wind load of radio telescope reflecting surface is obtained, which can effectively improve the calculation accuracy of average wind load, and provide support for studying the influence of wind load on radio telescope, radio telescope wind-resistant structure design and wind-resistant servo control design.

Description

Technical Field

[0001] This invention belongs to the field of radio telescope technology, specifically relating to a method for calculating the average wind load on the reflector surface of an ultra-large aperture radio telescope. Background Technology

[0002] With the development of radio telescope technology, people's demand for radio telescopes to "see clearly", "see accurately", and "see steadily" has become increasingly strong. Therefore, large aperture, high precision, and high gain have become the inevitable trend in the development of radio telescopes.

[0003] As the aperture of radio telescopes continues to increase, higher demands are placed on their pointing accuracy. Pointing accuracy is easily affected by disturbances such as gusts of wind, gravity, and ice and snow, leading to inaccurate or unstable pointing. With the increase in antenna aperture, its windward area and structural flexibility also increase, making the impact of wind disturbances more significant. Engineers have adopted various measures to overcome the effects of wind loads. Firstly, the Parke radio telescope reduces wind resistance by designing its antenna dish with a mesh structure, but this limits its operating frequency range to a lower range. Secondly, the LMT radio telescope reduces oscillations by increasing the reflective surface of the antenna dish and the weight of the antenna, but this increases its gravitational load. Thirdly, the GBT radio telescope compensates for reflective surface deformation by installing actuators, but these actuators are expensive and easily damaged. Fourthly, the Qinghai Delingha radio telescope isolates itself from wind interference by installing a radome, but this reduces antenna gain, is expensive, and its size cannot meet the requirements for large-aperture antennas. Fifthly, Wodek Gawronski. Modeling and Control of Antennas and Telescopes [M]. Springer, Boston, MA. It is easy to use a servo system for wind interference control by calculating the wind load on the antenna in advance. The key to this is calculating the wind load on the antenna.

[0004] Regarding the calculation of wind load on antennas, in 1964, J.S. Brown and K.M. Keckere conducted wind tunnel experiments on a solid parabolic antenna with a focal diameter ratio of 0.33, obtaining wind load coefficients under different operating conditions, which are applicable to engineering design. my country has also conducted extensive research on the wind load of parabolic antennas, creating scaled-down physical models and simulating actual wind environments and structures in wind tunnels to obtain the antenna's wind load coefficient curve. (References: Gong Zhenbang, Servo Mechanical Transmission Device [M]. National Defense Industry Press, 1980; and Jin Guohua, Antenna Wind Load Analysis and Calculation Software Development [D]. Xi'an University of Electronic Science and Technology, 2013.) Based on the basic characteristics of wind load and dimensional analysis, the formula for calculating the average wind load of a parabolic antenna is derived.

[0005] However, they treat the wind and wind load coefficients experienced by all parts of the antenna reflector as equal, which is within acceptable limits for small-aperture antennas. But for very large-aperture radio telescopes, the wind speed difference between the top and bottom can be as much as two times, and this can be amplified fourfold through dimensional amplification. Furthermore, the significant stiffness difference between the center and edges of a large-aperture antenna leads to a substantial difference in its load coefficient. Therefore, traditional calculation methods for very large-aperture radio telescopes introduce considerable errors, affecting the structural design and wind disturbance control accuracy of the telescope. Summary of the Invention

[0006] To address the aforementioned problems in the existing technology, this invention provides a method for calculating the average wind load on the reflecting surface of an ultra-large aperture radio telescope. The technical problem to be solved by this invention is achieved through the following technical solution:

[0007] This invention provides a method for calculating the average wind load on the reflecting surface of an ultra-large aperture radio telescope, including:

[0008] Step 1: Determine the structural parameters and operating conditions of the radio telescope;

[0009] Step 2: Based on the structural parameters and operating conditions, obtain the height from the center of the equivalent aperture of the radio telescope's reflecting surface to the ground.

[0010] Step 3: Determine the wind speed gradient and divide the wind speed gradient according to the wind field distribution at the site;

[0011] Step 4: Divide the vertical projection circle of the radio telescope's reflecting surface into uniform areas according to the wind speed gradient;

[0012] Step 5: Divide the reflective surface of the radio telescope into a ring shape according to the wind speed gradient;

[0013] Step 6: Based on the division results of the projected circular surface and the annular division results of the radio telescope's reflecting surface, the reflecting surface of the radio telescope is divided into sections, and each sub-region is numbered.

[0014] Step 7: Classify the sub-regions according to their shape features and calculate the proportion coefficient of each sub-region;

[0015] Step 8: Calculate the wind force and wind moment corresponding to each sub-region based on the scaling factor of the sub-region and the center distance between each sub-region and the reflecting surface of the radio telescope;

[0016] Step 9: Calculate the average wind load on the radio telescope reflector based on the wind force and wind moment corresponding to each sub-region.

[0017] In one embodiment of the present invention, the structural parameters of the radio telescope include: the focal-to-diameter ratio of the reflecting surface of the radio telescope, the aperture area and aperture diameter of the reflecting surface, and the distance from the center of the reflecting surface to the ground; the operating conditions of the radio telescope include the antenna elevation angle.

[0018] In one embodiment of the present invention, the height from the center of the equivalent aperture of the radio telescope's reflecting surface to the ground is calculated according to the following formula:

[0019] ;

[0020] In the formula, This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground. This indicates the distance from the center of the reflecting surface to the ground. Indicates the focal length of the reflecting surface. This indicates the antenna elevation angle.

[0021] In one embodiment of the present invention, step 3 includes:

[0022] Based on the topography and climate of the site, determine the topography type of the site;

[0023] Based on the terrain type of the site and the measured wind speed data, the wind speed gradient was determined and divided.

[0024] In one embodiment of the present invention, step 4 includes:

[0025] Using the diameter of the projected circular surface, the projected circular surface is divided into two symmetrical semicircles.

[0026] According to unit length The two semicircular surfaces are divided evenly, wherein, , Indicates the wind speed gradient. This indicates the antenna elevation angle.

[0027] In one embodiment of the present invention, step 5 includes:

[0028] The radio telescope's reflective surface is divided into rings from the center outwards. The radius of the innermost circle is the wind speed gradient, and the width of the ring formed by two adjacent circles is the wind speed gradient.

[0029] In one embodiment of the present invention, step 6 includes:

[0030] Based on the division result of the upper semicircle of the projected circular surface, the divided areas are divided from bottom to top according to... Number them sequentially;

[0031] Based on the division result of the lower semicircle of the projected circular surface, the divided areas are divided from top to bottom according to... Number them sequentially;

[0032] Based on the annular division of the radio telescope's reflecting surface, the divided areas are arranged from the inside out according to... Number them sequentially;

[0033] The division result of the projection circle is projected onto the reflecting surface of the radio telescope and combined with the annular division result of the reflecting surface of the radio telescope to achieve the partitioning of the reflecting surface of the radio telescope.

[0034] Each sub-region is numbered according to the region division number, and denoted as . ,in, This indicates the partition region number in the annular partitioning result of the reflecting surface, where the sub-region is located. m This indicates the number of the partitioned region in the result of the projection circle.

[0035] In one embodiment of the present invention, in step 7, the proportional coefficient of each sub-region is calculated according to the following formula, wherein,

[0036] when hour,

[0037] ;

[0038] In the formula, Subregion The corresponding proportionality coefficient;

[0039] when hour,

[0040] ;

[0041] when hour,

[0042] ;

[0043] when hour,

[0044] .

[0045] In one embodiment of the present invention, step 8 includes:

[0046] The height of each sub-region from the ground is calculated using the following formula:

[0047] ;

[0048] In the formula, Subregion Height to the ground This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground.

[0049] Based on the scaling factor of each sub-region and the height of each sub-region from the ground, the wind force corresponding to each sub-region is calculated according to the following formula:

[0050] ;

[0051] In the formula, Subregion Corresponding wind force Subregion The corresponding wind load factor, Indicates air density, This represents the average wind speed at the reference height. Indicates reference altitude. This represents the wind speed profile index. Indicates the antenna elevation angle;

[0052] Based on the projection relationship of the radio telescope's reflecting surface in the vertical direction, the distance from each sub-region to the midpoint of the radio telescope's reflecting surface is calculated as follows:

[0053] ;

[0054] In the formula, Subregion The distance to the midpoint of the radio telescope's reflecting surface. Indicates the focal length of the reflecting surface. This indicates the diameter of the reflecting surface. Represents a sub-region in the projected circular surface Distance to the center point of the projected circular surface;

[0055] Based on the wind force corresponding to each sub-region and the distance from each sub-region to the midpoint of the radio telescope's reflecting surface, the wind moment corresponding to each sub-region is calculated using the following formula:

[0056] ;

[0057] In the formula, Subregion The corresponding wind torque.

[0058] In one embodiment of the present invention, step 9 includes:

[0059] The wind force and wind moment corresponding to each sub-region are summed according to the following formula to obtain the average wind load on the radio telescope reflector.

[0060] ;

[0061] .

[0062] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0063] 1. The method for calculating the average wind load on the reflector surface of a large-aperture radio telescope of the present invention divides the reflector surface and wind speed gradient into regions based on the highly symmetrical characteristics of the radio telescope reflector surface and the law that wind increases with height. The wind load of each region is calculated separately and then summed to obtain the average wind load on the reflector surface of the radio telescope, which can effectively improve the calculation accuracy of the average wind load.

[0064] 2. The method for calculating the average wind load on the reflector surface of a large-aperture radio telescope of the present invention has high accuracy in calculating the average wind load, providing support for studying the influence of wind load on radio telescopes, the design of wind-resistant structures for radio telescopes, and the design of wind-resistant servo control.

[0065] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of the present invention more apparent and understandable, preferred embodiments are described in detail below with reference to the accompanying drawings. Attached Figure Description

[0066] Figure 1 This is a flowchart of a method for calculating the average wind load on the reflector surface of an ultra-large aperture radio telescope, provided in an embodiment of the present invention.

[0067] Figure 2 This is a schematic diagram of the structure of an ultra-large aperture radio telescope provided in an embodiment of the present invention;

[0068] Figure 3 This is a schematic diagram of the region division of the projection circle of the reflecting surface of a radio telescope in the vertical direction, according to an embodiment of the present invention.

[0069] Figure 4 This is a schematic diagram of the wind load calculation results of the reflecting surface of a radio telescope provided in an embodiment of the present invention. Detailed Implementation

[0070] To further illustrate the technical means and effects adopted by the present invention to achieve the intended purpose, the following describes in detail, with reference to the accompanying drawings and specific embodiments, a method for calculating the average wind load on the reflector surface of an ultra-large aperture radio telescope proposed according to the present invention.

[0071] The foregoing and other technical contents, features, and effects of the present invention will be clearly presented in the following detailed description of specific embodiments in conjunction with the accompanying drawings. Through the description of the specific embodiments, a more in-depth and concrete understanding can be gained of the technical means and effects adopted by the present invention to achieve its intended purpose. However, the accompanying drawings are for reference and illustration only and are not intended to limit the technical solutions of the present invention.

[0072] Example 1

[0073] Please see Figure 1 , Figure 1 This is a flowchart illustrating a method for calculating the average wind load on the reflector surface of a very large aperture radio telescope, as provided in this embodiment of the invention. The method includes:

[0074] Step 1: Determine the structural parameters and operating conditions of the radio telescope;

[0075] The structural parameters of a radio telescope include: the focal-to-diameter ratio of the reflector, the aperture area and diameter of the reflector, and the distance from the center of the reflector to the ground; the operating conditions of the radio telescope include the antenna elevation angle.

[0076] It should be noted that the aperture and focal diameter ratio of the radio telescope's reflecting surface are determined during the initial design phase and cannot be changed after construction. The focal diameter ratio affects the observation performance of the radio telescope. , Indicates the focal length of the reflecting surface. This indicates the diameter of the reflecting surface.

[0077] Step 2: Based on the structural parameters and operating conditions, obtain the height from the center of the equivalent aperture of the radio telescope's reflecting surface to the ground.

[0078] In one alternative implementation, such as Figure 2 The schematic diagram of the ultra-large aperture radio telescope shown is given. The height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground is calculated according to equation (1):

[0079] (1);

[0080] In the formula, This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground. This indicates the distance from the center of the reflecting surface to the ground. Indicates the focal length of the reflecting surface. This indicates the antenna elevation angle.

[0081] Step 3: Determine the wind speed gradient and divide the wind speed gradient according to the wind field distribution at the site;

[0082] In an optional implementation, step 3 includes:

[0083] Step 3.1: Determine the terrain type of the site based on its topography and climate;

[0084] The site, or installation location of the Very Large Aperture Radio Telescope, is determined according to national standards based on the site's topography and climate. Different topographic types correspond to different wind speed profile indices.

[0085] Step 3.2: Determine the wind speed gradient and divide the wind speed gradient according to the terrain type of the site and the measured wind speed data.

[0086] Optionally, a wind speed gradient of 10m can be selected according to national standards.

[0087] Step 4: Divide the vertical projection circle of the radio telescope's reflecting surface into uniform sections based on the wind speed gradient.

[0088] In an optional implementation, step 4 includes:

[0089] Step 4.1: Using the diameter of the projected circular surface, divide the projected circular surface into two symmetrical semicircles;

[0090] Step 4.2: According to unit length Divide the two semicircles evenly, where, , Indicates the wind speed gradient. This indicates the antenna elevation angle.

[0091] Step 5: Divide the reflector surface of the radio telescope into rings according to the wind speed gradient;

[0092] Because the reflector of a radio telescope has a highly symmetrical structure, the reflector of the radio telescope is divided into ring-shaped sections.

[0093] In an optional implementation, step 5 includes:

[0094] Divide the radio telescope's reflecting surface into rings from the center outwards. The radius of the innermost circle is the wind speed gradient, and the width of the ring formed by two adjacent circles is the wind speed gradient. Therefore, for any ring, its radius is a positive integer multiple of the radius of the innermost circle.

[0095] Step 6: Based on the division results of the projected circular surface and the annular division results of the radio telescope's reflecting surface, the reflecting surface of the radio telescope is divided into regions, and each sub-region is numbered.

[0096] See also Figure 3The illustrated embodiment of the present invention provides a schematic diagram of the region division of the projection circle of the reflecting surface of a radio telescope in the vertical direction. In an optional embodiment, step 6 includes:

[0097] Step 6.1: Based on the division result of the upper semicircle of the projected circular surface, divide the divided areas from bottom to top according to... Number them sequentially;

[0098] Step 6.2: Based on the division result of the lower semicircle of the projected circular surface, divide the divided area from top to bottom according to... Number them sequentially;

[0099] Step 6.3: Based on the annular division of the radio telescope's reflecting surface, divide the area from the inside out according to... Number them sequentially;

[0100] Step 6.4: Project the division result of the projection circle onto the reflecting surface of the radio telescope, and combine it with the annular division result of the reflecting surface of the radio telescope to achieve the partitioning of the reflecting surface of the radio telescope.

[0101] Step 6.5: Number each sub-region according to the region division number, denoted as ,in, This indicates the partition region number in the annular partitioning result of the reflecting surface, where the sub-region is located. m This indicates the number of the partitioned region in the result of the projection circle.

[0102] In this embodiment, , , , All are integers. For any subregion The distance from the center point of the projection circle to the region within the projection circle, which is also the radius of the region within the projection circle, is:

[0103] (2);

[0104] (3);

[0105] in, Represents a sub-region in the projected circular surface The distance to the center point of the projected circle, which is the sub-region The radius in the projected circular plane.

[0106] Step 7: Classify the sub-regions according to their shape characteristics and calculate the proportion coefficient of each sub-region;

[0107] In one optional implementation, all sub-regions are divided into four categories based on shape characteristics, including: , , and .

[0108] In this embodiment, the scaling factor for each sub-region is calculated according to equations (4)-(7), where,

[0109] when hour,

[0110] (4);

[0111] In the formula, Subregion The corresponding proportionality coefficient;

[0112] when hour,

[0113] (5);

[0114] when hour,

[0115] (6);

[0116] when hour,

[0117] (7).

[0118] Step 8: Calculate the wind force and wind moment corresponding to each sub-region based on the scale factor of the sub-region and the center distance between each sub-region and the reflector surface of the radio telescope;

[0119] Specifically, the formula for calculating wind load is as follows:

[0120] (8);

[0121] (9);

[0122] in, F Indicates wind force. M Indicates wind torque, Indicates the wind load factor. Indicates dynamic pressure head. It is the aperture area. It is the aperture diameter. , Indicates air density, Indicates distance from the ground The average wind speed at that location.

[0123] In this embodiment, a wind load calculation formula for each sub-region is constructed based on the wind load calculation formula, and the wind force and wind moment corresponding to each sub-region are calculated.

[0124] In an optional implementation, step 8 includes:

[0125] Step 8.1: Calculate the height of each sub-region from the ground according to equation (10):

[0126] (10);

[0127] In the formula, Subregion Height to the ground This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground.

[0128] Step 8.2: Based on the scaling factor of each sub-region and the height of each sub-region from the ground, calculate the wind force corresponding to each sub-region according to equation (11):

[0129] (11);

[0130] In the formula, Subregion Corresponding wind force Subregion The corresponding wind load factor, Indicates air density, This represents the average wind speed at the reference height. Indicates reference altitude. This represents the wind speed profile index. Indicates the antenna elevation angle;

[0131] Step 8.3: Based on the projection relationship of the radio telescope's reflecting surface in the vertical direction, the distance from each sub-region to the midpoint of the radio telescope's reflecting surface is calculated as follows:

[0132] (12);

[0133] In the formula, Subregion The distance to the midpoint of the radio telescope's reflecting surface. Indicates the focal length of the reflecting surface. This indicates the diameter of the reflecting surface. Represents a sub-region in the projected circular surface Distance to the center point of the projected circular surface;

[0134] Step 8.4: Based on the wind force corresponding to each sub-region and the distance from each sub-region to the midpoint of the radio telescope's reflecting surface, calculate the wind moment corresponding to each sub-region according to equation (13):

[0135] (13);

[0136] In the formula, Subregion The corresponding wind torque.

[0137] It should be noted that, in this embodiment, for ease of simulation, wind speed is simulated using an exponential law, and the simulation formula is as follows:

[0138] (14);

[0139] in, Indicates distance from the ground The average wind speed at that location This represents the average wind speed at the reference height. Indicates reference altitude. This indicates the wind speed profile index.

[0140] Generally, a reference height of 10m is taken, and the average wind speed at the reference height is the average wind speed at 10m above the ground.

[0141] Step 9: Calculate the average wind load on the radio telescope reflector based on the wind force and wind moment corresponding to each sub-region.

[0142] In an optional implementation, step 9 includes:

[0143] The wind force and wind moment corresponding to each sub-region are summed according to equations (15)-(16) to obtain the average wind load on the radio telescope reflector.

[0144] (15);

[0145] (16).

[0146] It should be noted that, for the reflector of the radio telescope, since the wind torque on the upper half of the reflector is opposite in direction to the wind torque on the lower half of the reflector, the total wind torque on the radio telescope is obtained by subtracting the wind torque on the upper half of the reflector from the wind torque on the lower half of the reflector in equation (16).

[0147] The present invention provides a method for calculating the average wind load on the reflector surface of a large-aperture radio telescope. Based on the highly symmetrical characteristics of the radio telescope reflector surface and the law that wind speed increases with altitude, the reflector surface and wind speed gradient of the radio telescope are divided into regions. The wind load of each region is calculated separately and then summed to obtain the average wind load on the reflector surface of the radio telescope. This method can effectively improve the accuracy of the average wind load calculation and provide support for studying the influence of wind load on radio telescopes, the design of wind-resistant structures for radio telescopes, and the design of wind-resistant servo control.

[0148] Example 2

[0149] This embodiment illustrates the specific process and effect of the method for calculating the average wind load on the reflector surface of the ultra-large aperture radio telescope in Embodiment 1 through simulation experiments.

[0150] The first step is to determine the structural parameters and operating conditions of the radio telescope, and to calculate the height of the equivalent aperture center from the ground. This example uses the QTT radio telescope, currently under construction, with an aperture of 110m, a focal diameter ratio (P / D) of 0.3, an antenna elevation angle of 0°, a focal length of 33m, and a height of 60m from the midpoint of the radio telescope's reflecting surface to the ground.

[0151] The second step involves determining the wind speed gradients to be divided based on the wind field distribution at the site, and then partitioning the vertical projection of the radio telescope's reflecting surface. The site's topography is analyzed to be of type B, with a wind speed profile index... Taking a value of 0.16, the local air density in Qitai is 1.0197 kg / m³, and the wind speed gradient is 5 m. Therefore, the reflector surface of the radio telescope is divided into 22 sections, each with a height of 5 m.

[0152] The third step involves dividing the radio telescope's reflector surface into annular zones with a minimum radius of 5m. Combined with wind speed divisions, this results in 220 sub-regions. The scaling factor for each sub-region is shown in Table 1. The distances from each sub-region to the ground are 5m, 10m, 15m, 20m, 25m, 30m, 35m, 40m, 45m, 50m, 55m, 60m, 65m, 70m, 75m, 80m, 85m, 90m, 95m, 100m, 105m, and 110m. For example, based on measured data, the average annual wind speed at a height of 10m at Qitai is 4m / s. The calculated wind speed at 110m is 8m / s, obtained from wind tunnel data. .

[0153] Table 1. Proportional coefficients for each sub-region

[0154]

[0155] Fourth step, substituting each parameter into equation (11) yields the wind force experienced by any sub-region as follows: The distance from any sub-region to the midpoint of the antenna reflector is The forces acting on the reflector sections at each wind speed gradient are shown in Table 2. The total wind torque acting on the radio telescope is: .

[0156] Table 2. Forces on antenna reflector sections at different wind speed gradients

[0157]

[0158] As can be seen from the above experiments, the method for calculating the average wind load on the reflector surface of the ultra-large aperture radio telescope of this invention can significantly improve the accuracy of the calculation by dividing the radio telescope reflector surface into sections and considering various factors such as wind field characteristics, antenna structure, antenna operating conditions, and site wind field characteristics. The calculation process is relatively simple and the physical meaning is relatively clear, providing support for further research on the impact of wind load on antennas, antenna wind-disruption-resistant structural design, and antenna wind-disruption-resistant control.

[0159] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations are intended to cover non-exclusive inclusion, such that an article or device comprising a list of elements includes not only those elements but also other elements not expressly listed. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the article or device comprising said element. Terms such as "connected" or "linked" are not limited to physical or mechanical connections but can include electrical connections, whether direct or indirect. The orientations or positional relationships indicated by terms such as "upper," "lower," "left," and "right" are based on the orientations or positional relationships shown in the accompanying drawings and are used only for the convenience of describing the invention and for simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as limiting the invention.

[0160] The above description, in conjunction with specific preferred embodiments, provides a further detailed explanation of the present invention. It should not be construed that the specific implementation of the present invention is limited to these descriptions. For those skilled in the art, various simple deductions or substitutions can be made without departing from the concept of the present invention, and all such modifications and substitutions should be considered within the scope of protection of the present invention.

Claims

1. A method for calculating the average wind load on the reflecting surface of an ultra-large aperture radio telescope, characterized in that, include: Step 1: Determine the structural parameters and operating conditions of the radio telescope; Step 2: Based on the structural parameters and operating conditions, obtain the height from the center of the equivalent aperture of the radio telescope's reflecting surface to the ground. Step 3: Determine the wind speed gradient and divide the wind speed gradient according to the wind field distribution at the site; Step 4: Divide the vertical projection circle of the radio telescope's reflecting surface into uniform areas according to the wind speed gradient; Step 5: Divide the reflective surface of the radio telescope into a ring shape according to the wind speed gradient; Step 6: Based on the division results of the projected circular surface and the annular division results of the radio telescope's reflecting surface, the reflecting surface of the radio telescope is divided into zones, and each sub-region is numbered; Step 6 includes: Based on the division result of the upper semicircle of the projected circular surface, the divided areas are divided from bottom to top according to... Number them sequentially; Based on the division result of the lower semicircle of the projected circular surface, the divided areas are divided from top to bottom according to... Number them sequentially; Based on the annular division of the radio telescope's reflecting surface, the divided areas are arranged from the inside out according to... Number them sequentially; The division result of the projection circle is projected onto the reflecting surface of the radio telescope and combined with the annular division result of the reflecting surface of the radio telescope to achieve the partitioning of the reflecting surface of the radio telescope. Each sub-region is numbered according to the region division number, and denoted as . ,in, This indicates the partition region number in the annular partitioning result of the reflecting surface, where the sub-region is located. m This indicates the number of the partitioned region in the result of the projection circle; Step 7: Classify the sub-regions according to their shape features and calculate the proportion coefficient of each sub-region; Step 8: Calculate the wind force and wind moment corresponding to each sub-region based on the scaling factor of the sub-region and the center distance between each sub-region and the reflecting surface of the radio telescope; Step 9: Calculate the average wind load on the radio telescope reflector based on the wind force and wind moment corresponding to each sub-region.

2. The method for calculating the average wind load on the reflector surface of a large-aperture radio telescope according to claim 1, characterized in that, The structural parameters of the radio telescope include: the focal-to-diameter ratio of the reflector, the aperture area and diameter of the reflector, and the distance from the center of the reflector to the ground; the operating conditions of the radio telescope include the antenna elevation angle.

3. The method for calculating the average wind load on the reflecting surface of a large-aperture radio telescope according to claim 1, characterized in that, The height from the center of the equivalent aperture of the radio telescope's reflecting surface to the ground is calculated using the following formula: ； In the formula, This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground. This indicates the distance from the center of the reflecting surface to the ground. Indicates the focal length of the reflecting surface. This indicates the antenna elevation angle.

4. The method for calculating the average wind load on the reflector surface of a large-aperture radio telescope according to claim 1, characterized in that, Step 3 includes: Based on the topography and climate of the site, determine the topography type of the site; Based on the terrain type of the site and the measured wind speed data, the wind speed gradient was determined and divided.

5. The method for calculating the average wind load on the reflecting surface of a large-aperture radio telescope according to claim 1, characterized in that, Step 4 includes: Using the diameter of the projected circular surface, the projected circular surface is divided into two symmetrical semicircles. According to unit length The two semicircular surfaces are divided evenly, wherein, , Indicates the wind speed gradient. This indicates the antenna elevation angle.

6. The method for calculating the average wind load on the reflecting surface of a large-aperture radio telescope according to claim 5, characterized in that, Step 5 includes: The radio telescope's reflective surface is divided into rings from the center outwards. The radius of the innermost circle is the wind speed gradient, and the width of the ring formed by two adjacent circles is the wind speed gradient.

7. The method for calculating the average wind load on the reflecting surface of a super-large aperture radio telescope according to claim 1, characterized in that, In step 7, the scaling factor for each sub-region is calculated according to the following formula, where, when hour, ； In the formula, Subregion The corresponding proportionality coefficient; when hour, ； when hour, ； when hour, 。 8. The method for calculating the average wind load on the reflecting surface of a large-aperture radio telescope according to claim 7, characterized in that, Step 8 includes: The height of each sub-region from the ground is calculated using the following formula: ； In the formula, Subregion Height to the ground This indicates the height of the center of the equivalent aperture of the radio telescope's reflecting surface from the ground. Indicates the wind speed gradient; Based on the scaling factor of each sub-region and the height of each sub-region from the ground, the wind force corresponding to each sub-region is calculated according to the following formula: ； In the formula, Subregion Corresponding wind force Subregion The corresponding wind load factor, Indicates air density, This represents the average wind speed at the reference height. Indicates reference altitude. This represents the wind speed profile index. Indicates the antenna elevation angle; Based on the projection relationship of the radio telescope's reflecting surface in the vertical direction, the distance from each sub-region to the midpoint of the radio telescope's reflecting surface is calculated as follows: ； In the formula, Subregion The distance to the midpoint of the radio telescope's reflecting surface. Indicates the focal length of the reflecting surface. This indicates the diameter of the reflecting surface. Represents a sub-region in the projected circular surface Distance to the center point of the projected circular surface; Based on the wind force corresponding to each sub-region and the distance from each sub-region to the midpoint of the radio telescope's reflecting surface, the wind moment corresponding to each sub-region is calculated using the following formula: ； In the formula, Subregion The corresponding wind torque.

9. The method for calculating the average wind load on the reflecting surface of a super-large aperture radio telescope according to claim 8, characterized in that, Step 9 includes: The wind force and wind moment corresponding to each sub-region are summed according to the following formula to obtain the average wind load on the radio telescope reflector. ； 。

Images