Three-degree-of-freedom fluid-structure interaction simulation method for icing conductor galloping

Abstract

The application discloses a three-degree-of-freedom fluid-structure coupling simulation method for icing conductor galloping, comprising: generating fluctuating wind speed time history, simulating the flow field around the conductor and calculating the supercooled water droplet collision distribution and icing growth rate, dynamically updating the circumferential and spanwise non-uniform distribution of the icing shape, and then establishing vertical, horizontal and torsional three-degree-of-freedom dynamic equations containing mass eccentric inertia coupling effect, adopting a partitioned two-way fluid-structure coupling strategy, alternately executing flow field solving and time domain integration in each time step, updating the conductor motion boundary and iterating until the simulation is terminated, and finally determining the galloping critical condition and dynamic characteristics based on the time history response. The application overcomes the precision defects of the traditional uniform icing and steady wind field assumptions, and can provide a high-precision numerical simulation tool for transmission line galloping risk assessment and anti-galloping design.

Description

Technical Field

[0001] This application relates to the field of fluid dynamics coupling simulation technology, and more specifically, to a three-degree-of-freedom fluid-structure interaction simulation method for the dancing of an icy conductor. Background Technology

[0002] As a critical component of the power system, the operating status of power conductors directly affects the reliability and economy of the power grid. Icing galloping, a serious natural disaster, is one of the main factors threatening the safety of transmission lines, easily causing accidents such as flashover, conductor strand breakage, and even tower collapse, posing a severe challenge to the stable operation of the power grid.

[0003] However, existing methods for predicting galloping still have significant shortcomings: Regarding icing morphology, current studies often assume uniform or regular geometric shapes, but actual icing is controlled by atmospheric turbulence, and the collision efficiency of supercooled water droplets exhibits a significant non-uniform distribution along the circumference and axial direction of the conductor. The assumption of uniform icing cannot reflect its true impact on conductor mass eccentricity and aerodynamic characteristics. Regarding wind input, most models treat the wind field as a steady-state uniform flow, failing to reflect the continuous modulation of the effective angle of attack by pulsating wind and its influence on the galloping initiation threshold. Regarding multi-degree-of-freedom coupling description, existing models often decouple the vertical, horizontal, and torsional degrees of freedom. However, iced conductor galloping is essentially a self-excited vibration with strong three-degree-of-freedom coupling. Torsional motion changes the effective angle of attack, thus modulating aerodynamic forces, while vertical and horizontal vibrations, in turn, affect torsional excitation. The lack of a two-way feedback mechanism severely restricts the model's prediction accuracy. Regarding fluid-structure interaction, traditional methods simplify aerodynamic forces to steady functions only related to the incoming flow, ignoring the reaction of the conductor's motion to the surrounding flow field, and failing to capture the phase relationship between aerodynamic forces and motion.

[0004] The aforementioned shortcomings lead to significant deviations in the prediction of galloping response by existing models in real complex environments. There is an urgent need for a high-precision galloping simulation method that can comprehensively consider the dynamic evolution of icing morphology, the unsteady characteristics of pulsating wind fields, and three-degree-of-freedom two-way fluid-structure interaction. Summary of the Invention

[0005] The purpose of this application is to provide a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of icy conductors, which aims to solve the problem that the existing technology uses the assumptions of uniform icing and steady-state wind field, resulting in low accuracy of galloping simulation.

[0006] The first aspect of this invention provides a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor, comprising: Generating fluctuating wind speed time histories based on wind speed power spectrum; The flow field around the conductor is simulated based on the time history of the pulsating wind speed. The collision distribution and icing growth rate of supercooled water droplets on the surface of the conductor are calculated. The icing geometry of the conductor surface is updated based on the collision distribution and the icing growth rate to obtain a non-uniform icing morphology. Based on the non-uniformly distributed icing morphology, a dynamic equation containing inertial coupling and three degrees of freedom is constructed; The three-degree-of-freedom dynamic equations are coupled with the flow field solution. In each time step, the flow field is first solved to obtain transient aerodynamic forces and torques. The transient aerodynamic forces and torques are then substituted into the three-degree-of-freedom dynamic equations to obtain the displacement response. The conductor boundary is updated based on the displacement response, and the time history response of the icy conductor is output after the time step iteration terminates. The critical conditions and dynamic characteristics of icing conductor galloping are determined based on the time history response.

[0007] In one optional implementation, the fluctuating wind speed time history is generated based on the wind speed power spectrum, specifically as follows: The wind speed power spectrum is discretized by frequency, and a random phase is assigned to each frequency. The signal is then converted into a time-domain signal by inverse Fourier transform. After zero-mean and standardization processing of the time-domain signal, the time history of the pulsating wind speed is obtained.

[0008] In one optional implementation, the flow field around the conductor is simulated based on the pulsating wind speed time history to calculate the collision distribution of supercooled water droplets on the conductor surface and the icing growth rate, specifically: Using the fluctuating wind speed time history as the inlet boundary condition of the flow field, the continuous phase flow field around the conductor is solved; Introducing a discrete phase of supercooled water droplets into a continuous flow field, tracking the motion trajectory of the water droplets and statistically analyzing the positions and fluxes of the impact on the conductor surface, the collision distribution on the conductor surface is obtained. Based on the collision distribution and the mass conservation relationship of the water droplet freezing process, the ice growth rate at each location on the conductor surface is determined.

[0009] In one optional implementation, the icing geometry of the conductor surface is updated based on the collision distribution and the icing growth rate, specifically as follows: Divide the surface into several micro-elements along the circumference of the conductor. Based on the collision distribution and ice growth rate of each micro-element, accumulate the ice layer thickness and update the outer contour coordinates of the ice layer by micro-element. The geometric boundary of the conductor section is reconstructed based on the updated outer contour coordinates, forming a non-uniformly distributed icing pattern.

[0010] In one optional implementation, a dynamic equation comprising inertial coupling three degrees of freedom is constructed based on the non-uniformly distributed icing morphology, specifically: The mass eccentricity parameter and rotational inertia parameter of the conductor are determined based on the non-uniform icing morphology. The unit-length iced conductor is transformed into a mass-spring-damped system containing vertical, horizontal, and torsional degrees of freedom. Introducing mass eccentricity and moment of inertia parameters, the three-degree-of-freedom dynamic equations are obtained, expressed as follows:

[0011] In the formula, The mass per unit length of the conductor including ice. , and These are vertical acceleration, horizontal acceleration, and torsional angular acceleration, respectively. This represents the offset between the ice-covered centroid and the conductor centroid. , and These are vertical velocity, horizontal velocity, and torsional angular velocity, respectively. , and These are vertical displacement, horizontal displacement, and torsional angular displacement, respectively. , and These are the structural damping coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are the stiffness coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are instantaneous lift, drag, and aerodynamic torque, respectively. This is the moment of inertia of a unit length of conductor about its centroid.

[0012] In one optional implementation, the three-degree-of-freedom dynamic equations are coupled with the flow field solution, and the transient aerodynamic forces and torques are obtained by solving the flow field at each time step, specifically as follows: The flow field is solved first in each time step. The fluid control equations are solved by computational fluid dynamics and a turbulence model is introduced to obtain the transient lift, drag and aerodynamic torque acting on the surface of the conductor. By substituting the transient lift, drag, and aerodynamic torque into the three-degree-of-freedom dynamic equations and using the time-domain integration method for discrete solution, the vertical displacement, horizontal displacement, and torsional angular displacement at the current time step are obtained.

[0013] In one optional implementation, the conductor boundary is updated based on the displacement response, and the time history response of the icy conductor is output after the time step iteration terminates, specifically as follows: Based on the displacement response obtained at the current time step, the conductor boundary is driven to synchronously perform vertical translation, horizontal translation, and torsional rotation, and the positions of the mesh nodes in the computational domain are updated synchronously. Determine whether the current time step has reached the preset total simulation duration. If not, proceed to the next time step to continue iterating. If it has reached the preset total simulation duration, terminate the iteration calculation and output the time history data of displacement, velocity, acceleration and aerodynamic force corresponding to each time step of the icy conductor, which constitute the time history response of the icy conductor.

[0014] A second aspect of the present invention provides an electronic device, characterized in that it includes: a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor.

[0015] A third aspect of the present invention provides a computer-readable storage medium, characterized in that a computer program is stored on the computer-readable storage medium, and the computer program, when executed by a processor, is a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor.

[0016] This application has at least the following advantages or beneficial effects: (1) This invention obtains the real ice shape with non-uniform distribution in the circumferential and spanwise directions by simulating the dynamic evolution of ice morphology, which overcomes the problem that the regular cross section assumption deviates greatly from the actual ice morphology, and provides more reliable geometric input for glaciation analysis.

[0017] (2) This invention introduces the time history of fluctuating wind speed as the boundary condition of the flow field, and combines it with the partitioned two-way fluid-structure coupling solution to make up for the defect that the quasi-static aerodynamic model cannot capture unsteady aerodynamic effects, and is more in line with the actual atmospheric boundary layer conditions.

[0018] (3) The present invention establishes a three-degree-of-freedom dynamic equation with mass eccentric inertial coupling, which fully describes the strong nonlinear coupling in the vertical, horizontal and torsional directions, and can more accurately reflect the self-excited vibration mechanism of the galloping. Attached Figure Description

[0019] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] Figure 1 This is a flowchart of a three-degree-of-freedom fluid-structure interaction simulation method for ice-covered conductor galloping provided by the present invention; Figure 2This is a schematic diagram of an electronic device according to this application. Detailed Implementation

[0021] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0022] Please refer to Figure 1 , Figure 1 This is a flowchart of a three-degree-of-freedom fluid-structure interaction simulation method for icy conductor galloping, proposed in an embodiment of this application, including the following steps: Step S1: Generate the time history of fluctuating wind speed based on the wind speed power spectrum; Specifically, the time history of fluctuating wind speed is generated based on the wind speed power spectrum, as follows: The wind speed power spectrum is discretized by frequency, and a random phase is assigned to each frequency. The signal is then converted into a time-domain signal by inverse Fourier transform. After zero-mean and standardization processing of the time-domain signal, the time history of the pulsating wind speed is obtained.

[0023] In this embodiment, the time history of fluctuating wind speed is generated based on the von Kármán wind speed spectrum, wherein the expression for the power spectral density of wind speed is:

[0024] In the formula, The power spectral density is the wind speed. For turbulence intensity, The average wind speed, For turbulent integral scale, For frequency; Furthermore, frequency discretization of the power spectral density is performed to set the frequency resolution. The amplitude of the frequency component is calculated using the following expression:

[0025] In the formula, The amplitude of the frequency component; Each frequency component is assigned a random phase, which follows an interval... The uniform distribution is expressed as:

[0026] In the formula, For the complex wind speed spectrum, The imaginary unit, It is a random phase; For real-valued signals, the spectrum must satisfy conjugate symmetry, meaning the negative frequency components are the complex conjugates of the positive frequency components. After generating the complex spectrum of all frequency components in the frequency domain, it is converted into a time-domain signal using an inverse Fourier transform. To ensure the signal's standard deviation meets a specified value, the signal is zero-meaned and its standard deviation is calibrated to obtain the pulsating wind speed time history.

[0027] Step S2: Simulate the flow field around the conductor based on the pulsed wind speed time history, calculate the collision distribution and icing growth rate of supercooled water droplets on the conductor surface, update the icing geometry on the conductor surface based on the collision distribution and the icing growth rate, and obtain a non-uniformly distributed icing morphology. Existing studies typically assume that icing is uniformly distributed, but under the influence of pulsating wind fields, conductor icing exhibits a significant non-uniform distribution.

[0028] In this embodiment, the generated fluctuating wind speed time history is used as the inlet boundary condition for the flow field to numerically simulate the gas-liquid two-phase flow field around the conductor. In the continuous phase (gas phase) solution, the Reynolds-averaged method (RANS) is used to numerically simulate the flow field around the conductor.

[0029] Among them, the Reynolds-averaged method decomposes transient flow into average and fluctuating quantities by performing time averaging on the Navier-Stokes equations.

[0030] In this embodiment, the shear stress transport (SST) k-ω turbulence model is preferably used to solve the governing equations in a closed loop.

[0031] The shear stress transport (SST) k-ω turbulence model includes the turbulent kinetic energy k equation and the specific dissipation rate ω equation. For example, the expression for the turbulent kinetic energy k equation is:

[0032] The expression for the specific dissipation rate ω transport equation is:

[0033] In the formula, In this embodiment, the trajectory of supercooled water droplets in the flow field is tracked using a discrete phase model (DPM). As discrete phase particles, the motion equations of supercooled water droplets take into account inertial forces, drag forces, gravity, and turbulent diffusion effects.

[0034] Furthermore, based on the intersection of the particle trajectory and the conductor wall, the collision positions and collision flux of each particle with the conductor surface are statistically analyzed, and the local collection coefficient and collision efficiency on each surface micro-element are calculated, thereby obtaining the collision distribution of supercooled water droplets on the conductor surface.

[0035] Based on the obtained collision distribution, the mass flow rate of water droplets captured at each location on the conductor surface is determined. Then, the mass conservation equations for the icing surface are solved simultaneously to determine the icing growth rate at each location.

[0036] The mass conservation equation is used to determine the ratio of the frozen portion to the overflow portion of the mass of water droplets captured at various locations on the conductor surface.

[0037] After determining the icing growth rate at each location, several surface micro-elements are divided along the circumference of the conductor. Based on the collision distribution and icing growth rate corresponding to each micro-element, the ice layer thickness is accumulated element by element and the outer contour coordinates of the ice layer are updated. Based on the updated outer contour coordinates, the geometric boundary of the conductor section is reconstructed to form an icing morphology that is non-uniformly distributed along the circumference and axis of the conductor.

[0038] Step S3: Construct a dynamic equation with inertial coupling of three degrees of freedom based on the non-uniformly distributed icing morphology; Based on the Euler-Bernoulli beam theory, the icy conductor per unit length is equivalent to a three-degree-of-freedom mass-spring-damped system containing vertical, horizontal, and torsional degrees of freedom, and vertical displacement, horizontal displacement, and torsional angular displacement are selected as generalized coordinates.

[0039] Determine the mass eccentricity parameter of the conductor cross-section based on the non-uniform icing pattern. and rotational inertia parameters .

[0040] Based on Newton's second law and the law of rigid body rotation, the three-degree-of-freedom dynamic equations of the icing conductor under the action of the flow field are established, and the expression is:

[0041] In the formula, The mass per unit length of the conductor including ice. , and These are vertical acceleration, horizontal acceleration, and torsional angular acceleration, respectively. This represents the offset between the ice-covered centroid and the conductor centroid. , and These are vertical velocity, horizontal velocity, and torsional angular velocity, respectively. , and These are vertical displacement, horizontal displacement, and torsional angular displacement, respectively. , and These are the structural damping coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are the stiffness coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are instantaneous lift, drag, and aerodynamic torque, respectively. This is the moment of inertia of a unit length of conductor about its centroid.

[0042] When the ice is evenly distributed The vertical, horizontal, and torsional degrees of freedom are independent of each other, and the vertical, horizontal, and torsional vibrations of the system do not affect each other. At this point, even if there is external aerodynamic excitation, there is no inertial coupling between the three degrees of freedom, and their respective responses are determined only by their corresponding stiffness and damping. When the ice cover is not uniformly distributed The three degrees of freedom are bidirectionally correlated through an inertial coupling term. Specifically, the torsional angular acceleration... pass and The term generates additional inertial forces in the vertical and horizontal equations, and the components of this force in the vertical and horizontal directions excite vertical and horizontal vibrations, respectively. Simultaneously, vertical acceleration... and horizontal acceleration pass The term generates an additional inertial torque in the torsion equation, which excites torsional vibration. A two-way, real-time inertial coupling closed loop is formed between the three degrees of freedom.

[0043] Step S4: Couple the three-degree-of-freedom dynamic equations with the flow field solution, and first solve the flow field in each time step to obtain transient aerodynamic forces and torques. Substitute the transient aerodynamic forces and torques into the three-degree-of-freedom dynamic equations to obtain the displacement response. In this embodiment, the geometric parameters, material parameters and initial conditions of the transmission line are input, an initial calculation grid is established, and the calculation time step and total simulation duration are set.

[0044] Within the current time step, the current boundary position of the conductor is used as the boundary of the computational domain. The Reynolds average method combined with the turbulence model is used to solve the flow field around the conductor, and the aerodynamic force distribution acting on the surface of the conductor is obtained, including instantaneous lift, drag and aerodynamic torque.

[0045] After the flow field calculation converges at the current time step via the first interface in the user-defined function, the aerodynamic distribution is read and substituted as the right-hand excitation term into the three-degree-of-freedom dynamic equation. The motion equation is then discretized using the time-domain integration method to obtain the displacement, velocity, and acceleration responses of the conductor in the vertical, horizontal, and torsional directions within the current time step.

[0046] In this embodiment, the fourth-order Runge-Kutta method is preferred as a time-domain integration method. Taking vertical motion as an example, its iterative format is as follows: Let the state variable (speed), (Acceleration), then the velocity and displacement are finally updated:

[0047] In the formula, and The first Time step and the Vertical velocity of the time step, and The first Time step and the Vertical displacement at the time step, , , and These represent the first, second, third, and fourth order slopes of the velocity increment, respectively. , , and These represent the first, second, third, and fourth order slopes of the displacement increment, respectively.

[0048] Step S5: Update the conductor boundary according to the displacement response, and output the time history response of the icy conductor after the time step iteration terminates; Specifically, the conductor boundary is updated based on the displacement response, and the time history response of the icy conductor is output after the time step iteration terminates, as follows: Based on the displacement response obtained at the current time step, the conductor boundary is driven to synchronously perform vertical translation, horizontal translation, and torsional rotation, and the positions of the mesh nodes in the computational domain are updated synchronously. Determine whether the current time step has reached the preset total simulation duration. If not, proceed to the next time step to continue iterating. If it has reached the preset total simulation duration, terminate the iteration calculation and output the time history data of displacement, velocity, acceleration and aerodynamic force corresponding to each time step of the icy conductor, which constitute the time history response of the icy conductor.

[0049] Step S6: Determine the critical conditions and dynamic characteristics of the ice-covered conductor galloping based on the time history response.

[0050] Specifically, the critical conditions include the critical wind speed for the galloping and the initial angle of attack range, and the dynamic characteristics include the amplitude and frequency of the galloping.

[0051] The method involves progressively increasing the average wind speed within a set range, and repeatedly applying the aforementioned fluid-structure interaction simulation for each wind speed condition to obtain the steady-state vibration amplitude of the conductor at each wind speed. A wind speed-amplitude curve is then plotted, and the wind speed at which the amplitude abruptly increases is identified as the critical wind speed for galloping.

[0052] In this process, simulations were performed for each angle of attack by changing the incoming wind direction angle, and the occurrence of significant conductor galloping was recorded at each angle of attack. The range of angles of attack that could induce galloping was defined as the initial angle of attack interval.

[0053] Among them, the steady-state vibration amplitude obtained under various wind speed conditions is plotted as a curve of wind speed variation to obtain the amplitude-wind speed relationship and obtain the galloping amplitude frequency.

[0054] This disclosure also provides an electronic device, please refer to... Figure 2 , Figure 2 This is a schematic diagram of an electronic device illustrated in an embodiment of this disclosure. For example... Figure 2 As shown, the electronic device 100 includes a memory 110 and a processor 120. The memory 110 and the processor 120 are connected via a bus for communication. The memory 110 stores a computer program that can run on the processor 120 to implement the steps in the three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor disclosed in this embodiment.

[0055] The disclosed embodiments also provide a computer-readable storage medium that, when executed by a processor of a computer device, enables the computer device to perform steps in a three-degree-of-freedom fluid-structure interaction simulation method for icing conductor galloping, as described in the present disclosure embodiments.

[0056] This application describes embodiments with reference to flowchart illustrations and / or block diagrams of methods, apparatuses, electronic devices, and computer program products according to embodiments of this application. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0057] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing terminal device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0058] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal equipment, causing a series of operational steps to be performed on the computer or other programmable terminal equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable terminal equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0059] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.

[0060] Although preferred embodiments of the present application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present application.

[0061] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only 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 thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.

[0062] The above provides a detailed description of a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the method and its core ideas. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor, characterized in that, include: Generating fluctuating wind speed time histories based on wind speed power spectrum; The flow field around the conductor is simulated based on the time history of the pulsating wind speed. The collision distribution and icing growth rate of supercooled water droplets on the surface of the conductor are calculated. The icing geometry of the conductor surface is updated based on the collision distribution and the icing growth rate to obtain a non-uniform icing morphology. Based on the non-uniformly distributed icing morphology, a dynamic equation containing inertial coupling and three degrees of freedom is constructed; The three-degree-of-freedom dynamic equations are coupled with the flow field solution. In each time step, the flow field is first solved to obtain transient aerodynamic forces and torques. The transient aerodynamic forces and torques are then substituted into the three-degree-of-freedom dynamic equations to obtain the displacement response. The conductor boundary is updated based on the displacement response, and the time history response of the icy conductor is output after the time step iteration terminates. The critical conditions and dynamic characteristics of icing conductor galloping are determined based on the time history response.

2. The three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The time history of fluctuating wind speed is generated based on the wind speed power spectrum, specifically as follows: The wind speed power spectrum is discretized by frequency, and a random phase is assigned to each frequency. The signal is then converted into a time-domain signal by inverse Fourier transform. After zero-mean and standardization processing of the time-domain signal, the time history of the pulsating wind speed is obtained.

3. The three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The flow field around the conductor is simulated based on the time history of the fluctuating wind speed, and the collision distribution of supercooled water droplets on the conductor surface and the icing growth rate are calculated, specifically: Using the fluctuating wind speed time history as the inlet boundary condition of the flow field, the continuous phase flow field around the conductor is solved; Introducing a discrete phase of supercooled water droplets into a continuous flow field, tracking the motion trajectory of the water droplets and statistically analyzing the positions and fluxes of the impact on the conductor surface, the collision distribution on the conductor surface is obtained. Based on the collision distribution and the mass conservation relationship of the water droplet freezing process, the ice growth rate at each location on the conductor surface is determined.

4. The three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The icing geometry of the conductor surface is updated based on the collision distribution and the icing growth rate, specifically as follows: Divide the surface into several micro-elements along the circumference of the conductor. Based on the collision distribution and ice growth rate of each micro-element, accumulate the ice layer thickness and update the outer contour coordinates of the ice layer by micro-element. The geometric boundary of the conductor section is reconstructed based on the updated outer contour coordinates, forming a non-uniformly distributed icing pattern.

5. A three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, Based on the non-uniformly distributed icing morphology, a dynamic equation involving inertial coupling and three degrees of freedom is constructed, specifically as follows: The mass eccentricity parameter and rotational inertia parameter of the conductor are determined based on the non-uniform icing morphology. The unit-length iced conductor is transformed into a mass-spring-damped system containing vertical, horizontal, and torsional degrees of freedom. Introducing mass eccentricity and moment of inertia parameters, the three-degree-of-freedom dynamic equations are obtained, expressed as follows: In the formula, The mass per unit length of the conductor including ice. , and These are vertical acceleration, horizontal acceleration, and torsional angular acceleration, respectively. This represents the offset between the ice-covered centroid and the conductor centroid. , and These are vertical velocity, horizontal velocity, and torsional angular velocity, respectively. , and These are vertical displacement, horizontal displacement, and torsional angular displacement, respectively. , and These are the structural damping coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are the stiffness coefficients in the vertical, horizontal, and torsional directions, respectively. , and These are instantaneous lift, drag, and aerodynamic torque, respectively. This is the moment of inertia of a unit length of conductor about its centroid.

6. The three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The three-degree-of-freedom dynamic equations are coupled with the flow field solution, and the transient aerodynamic forces and moments are obtained by solving the flow field in each time step. Specifically: The flow field is solved first in each time step. The fluid control equations are solved by computational fluid dynamics and a turbulence model is introduced to obtain the transient lift, drag and aerodynamic torque acting on the surface of the conductor. By substituting the transient lift, drag, and aerodynamic torque into the three-degree-of-freedom dynamic equations and using the time-domain integration method for discrete solution, the vertical displacement, horizontal displacement, and torsional angular displacement at the current time step are obtained.

7. The three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The conductor boundary is updated based on the displacement response. After the time step iteration terminates, the time history response of the icy conductor is output, specifically: Based on the displacement response obtained at the current time step, the conductor boundary is driven to synchronously perform vertical translation, horizontal translation, and torsional rotation, and the positions of the mesh nodes in the computational domain are updated synchronously. Determine whether the current time step has reached the preset total simulation duration. If not, proceed to the next time step to continue iterating. If it has reached the preset total simulation duration, terminate the iteration calculation and output the time history data of displacement, velocity, acceleration and aerodynamic force corresponding to each time step of the icy conductor, which constitute the time history response of the icy conductor.

8. A three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor according to claim 1, characterized in that, The critical conditions include the critical wind speed for the galloping and the initial angle of attack range, and the dynamic characteristics include the amplitude and frequency of the galloping.

9. An electronic device, characterized in that, include: A processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor as described in any one of claims 1 to 8.

10. A computer-readable storage medium, characterized in that, A computer program is stored on the computer-readable storage medium, which, when executed by a processor, implements a three-degree-of-freedom fluid-structure interaction simulation method for the galloping of an icy conductor as described in any one of claims 1 to 8.

Images