A HIGH-EFFICIENCY SIMULATION METHOD OF 3D WIND FIELD BASED ON DELAY EFFECT

Abstract

Determining coordinates and an initial coordinate system of simulation points according to structural drawings, and transforming the coordinate system so that Y axis is parallel to the wind direction, to obtain a 3D model of N structural simulation points; projecting all the simulation points of the 3D model onto a target 2D plane, and transforming the simulation points into projection points on the 2D plane; calculating the delay time of the wind speed, that is, the time required to move from each point to the projection points on the target plane; using a 2D coherence function to consider the spatial correlation of different simulation points in horizontal and vertical directions; generating fluctuating wind speed by using harmony superposition method; and obtaining the wind speed time history of all the points by using the delay time.

Claims

1. A high-efficiency simulation method of 3D wind field based on a delay effect, comprising the following steps: (1) determining coordinates and an initial coordinate system of simulation points according to structural drawings, and transforming the coordinate system so that Y axis is parallel to a wind direction, to obtain a 3D model of N structural simulation points; selecting a plane perpendicular to y=max(y.sub.1, y.sub.2, . . . , y.sub.N) as a target plane, where y.sub.1, y.sub.2, . . . , y.sub.N are the coordinates of the simulation points along the wind speed direction; projecting all the simulation points of the 3D model onto the target plane; and transforming the simulation points into projection points on a 2D plane; (2) according to the “Taylor's hypothesis”, considering the delay effect of wind speed and the discreteness of the wind speed during the actual simulation process, and calculating the delay time of the simulation point No. i, that is, the time required to move from the simulation points to the projection points at average wind speed: $\begin{array}{cc}{t}_{s1}\left(i\right)=.Math.\frac{\max \left(y_1,y_2,.Math.,y_N\right)-y_i}{V_m\left(z_i\right)*\Delta t}.Math.*\Delta t& \left(2\right)\end{array}$ where ┌*┐ represents rounding up to an integer, Δt is simulated time step, and V.sub.m represents the average wind speed; (3) generating fluctuating wind speed for the projection points in step (1) by using harmony superposition method; using a 2D coherence function to consider the spatial correlation of different simulation points in horizontal and vertical directions; for the longitudinal correlation of simulation points, using temporal correlation of the same simulation point at different times for calculation to obtain the wind speed time history of all the points; and generating the time history by time desynchronization of the fluctuating wind speed, wherein the wind speed of the point at any time t(f) is taken as:
V(x.sub.i,y.sub.i,z.sub.i,t(j))=V.sub.m(z.sub.i)+V.sub.f(x.sub.i,y.sub.i,z.sub.i,t(j)+t.sub.sl(i))  (1) in the formula, x.sub.i, y.sub.i, z.sub.i are an abscissa, an ordinate and a vertical coordinate of the simulation points, and V.sub.f represents the fluctuating wind speed.

Description

DESCRIPTION OF DRAWINGS

[0015] FIG. 1 is a schematic plane view of the structural simulation points projection method of the present invention;

[0016] FIG. 2 is a structural diagram of a transmission tower-line system selected in an implementation case of the present invention, and A, B and C in the figure are target points of axial force extraction after loading;

[0017] FIG. 3 is a comparison diagram of a target simulation point axial force result in an implementation case of the present invention and a axial force time history of a traditional method.

DETAILED DESCRIPTION

[0018] In order to make the invention purposes, features, and advantages of the present invention more obvious and understandable, the technical solutions in the embodiments of the present invention will be described clearly and completely below with reference to the drawings in the embodiments of the present invention. Obviously, the embodiments described below are only some, but not all of embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by those ordinary skilled in the art without contributing creative labor will belong to the protection scope of the present invention.

[0019] Referring to FIGS. 1-3, an embodiment of the present invention proposes a high-efficiency 3D wind field simulation method by taking a transmission tower-line system as an example.

[0020] Data source of implementation case: see “Fu X and L H N, Dynamic analysis of transmission tower-line system subjected to wind and rain loads, Journal of Wind Engineering and Industrial Aerodynamics, 2016, 157, 95-103”.

[0021] In the embodiment of the present invention, the establishment of a numerical model of a transmission tower can use self-programming or related commercial software. In the present embodiment, the widely used finite element analysis software ANSYS is used as an example to implement application of a time delay method of wind speed in applying wind loads to a transmission tower structure. In conjunction with the process shown in FIG. 1 and the technical solution of the present invention, specific illustration is as follows:

[0022] (1) The embodiment has three transmission towers with a height of 99.9 m, which are made of Q235 and Q345 equilateral angle steel, and the distance between every two towers is 500 m. The transmission towers are connected by wires. For information on tower structure and conductors, see related introduction in “Section 5” and “Section 6” of “Fu X and L H N, Dynamic analysis of transmission tower-line system subjected to wind and rain loads, Journal of Wind Engineering and Industrial Aerodynamics, 2016, 157, 95-103”. A transmission tower-line system which contains three towers and four lines finite element model is established by using ANSYS software. The ends of the conductors on both sides are in rigid connection. The BEAM188 element is selected to simulate the rods of the transmission towers. Rigid joints are used to simplify the connection between the members. An elasto-perfectly plastic model is adopted for a steel constitutive model.

[0023] (2) In the present embodiment, Davenport wind spectrum and Davenport coherence function models are used to simulate the spatial correlation of the fluctuating wind. The attenuation coefficients of the coherence function in x, y, and z directions are 16, 8 and 10 respectively, and an exponential law is used to simulate an average wind profile. The ground roughness is taken from Class B (i.e., α=0.15) of the Chinese standard “DL/T5551-2018. Load code for the design of overhead Transmission Line. National Energy Administration. 2018”, and the basic wind speed at 10 m height is 16 m/s.

[0024] (3) In the present embodiment, the traditional harmony superposition method considering 3D correlation and the time delay method of wind speed are used to generate wind loads. The time delay method of wind speed uses a 2D coherence function to consider the horizontal and longitudinal spatial correlation of different simulation points. For the spatial correlation of different longitudinal simulation points, according to the “Taylor's hypothesis”, temporal correlation of the same simulation point at different times is used for calculation to obtain the wind speed time history of all the points, and the simulation time is realized by time desynchronization of the fluctuating wind speed. The wind speed of the point at any time t can be calculated according to formula (4):

V(x.sub.i,y.sub.i,z.sub.i,t)=V.sub.m(z.sub.i)+V.sub.f(x.sub.i,y.sub.i,z.sub.i,t+t.sub.sl)  (4)

[0025] (4) The statistical characteristics of the force at target points A, B and C of the transmission tower-line system in the two cases are extracted and compared to illustrate the rationality of the time delay method of wind speed.

[0026] When the present invention issued, it should be noted that in actual simulation, the wind speed is discretized, and the formula (4) should be converted into the following formula (3) for calculation:

[00003]$\begin{array}{cc}V\left(x_i,y_i,z_i,t\left(j\right)\right)=V_m\left(z_i\right)+V_f\left(x_i,y_i,z_i,t\left(j\right)+t_{s1}\left(i\right)\right)\mathrm{where}{t}_{s1}\left(i\right)=.Math.\frac{\max \left(y_1,y_2,.Math.,y_N\right)-y_i}{V_m\left(z_i\right)*\Delta t}.Math.*\Delta t;.Math.*.Math.& \left(3\right)\end{array}$

represents rounding up to an integer; and Δt is simulated time step.

[0027] The above embodiments are only used to illustrate, but not to limit, the technical solutions of the present invention. Although the present invention has been described in detail with reference to the above embodiments, it should be understood by those ordinary skilled in the art that: the technical solutions recorded in the above embodiments can also be amended, or part of the technical features can be equivalently replaced. These amendments or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.