A Method and System for Equivalent Loading of Icy Conductor Galloping Load Based on Finite Element Analysis

Abstract

This invention discloses an equivalent loading method and system for galloping loads on icy power transmission lines based on finite element analysis. A conductor model is established based on actual line parameters, and the conductor is modeled and initial tension is applied using engineering simulation software. Three-dimensional aerodynamic parameters are obtained through wind tunnel tests at full wind angle of attack, and an aerodynamic load mathematical model is constructed to calculate the load on the conductor. The calculated load magnitude is then applied as a concentrated load at the nodes of the conductor. This invention utilizes a quasi-static load conversion strategy to equate dynamic aerodynamic forces to concentrated loads at the nodes. Combined with nonlinear conductor form-finding and wind tunnel test data, it achieves efficient and accurate galloping simulation. This solves the problems of low iteration efficiency and high model complexity in dynamic fluid-structure interaction, providing an efficient and reliable technical means for anti-galloping design and disaster prevention assessment of transmission lines.

Description

Technical Field

[0001] This invention belongs to the field of power transmission line technology. It relates to an equivalent loading method for galloping load on icy conductors based on finite element analysis. Background Technology

[0002] Transmission lines are prone to galloping under icing and strong winds, which is a large-amplitude, low-frequency oscillation of conductors caused by the coupling of aerodynamic and structural dynamics. This galloping can cause mechanical damage such as conductor strand breakage and hardware wear, and may also lead to serious accidents such as phase-to-phase flashover and power outages, threatening the safe operation of the power grid. Therefore, accurately simulating the characteristics of conductor galloping and predicting its dynamic response is of great significance for the anti-galloping design and disaster prevention assessment of transmission lines.

[0003] Traditional galloping analysis methods are mostly based on fluid-structure interaction (FSI) dynamic models, using iterative solutions by coupling computational fluid dynamics (CFD) and structural finite element method (FEM). However, such methods have significant limitations: on the one hand, dynamic FSI calculations require repeated iterations of data exchange between the fluid and structural domains, resulting in huge computational resource consumption and excessively long simulation times, making it difficult to meet the timeliness requirements of engineering designs; on the other hand, the complex time-varying nature of aerodynamic forces (such as pulsating wind loads and changes in icing shape) further exacerbates the uncertainty of the model, causing simulation results to often deviate significantly from measured data. In addition, existing equivalent loading methods often use simplified static loads to replace dynamic aerodynamic forces, but lack detailed consideration of the conductor's spatial morphology, initial tension distribution, and aerodynamic nonlinear characteristics, resulting in insufficient simulation accuracy.

[0004] In recent years, scholars have attempted to obtain the three-component force coefficients (drag, lift, and torque) of icing conductors through wind tunnel tests, and based on this, construct a mathematical model of aerodynamic loads, providing a theoretical basis for simplifying simulations. However, how to efficiently and accurately convert dynamic aerodynamic forces into static loads while maintaining the conductor's prestress state and spatial curvature characteristics remains a technical challenge. Furthermore, existing finite element modeling methods often neglect the accurate calculation of conductor form-finding, leading to discrepancies between the initial equilibrium state and actual working conditions, affecting the reliability of subsequent load analyses. Summary of the Invention

[0005] To address the problems of computational complexity, high resource consumption, and insufficient accuracy of existing equivalent loading techniques in traditional fluid-structure interaction methods, the present invention aims to provide an equivalent loading method and system for galloping loads on icy conductors based on finite element analysis. By employing a quasi-static load conversion strategy, dynamic aerodynamic forces are equivalent to nodal concentrated force loading. Combined with nonlinear form-finding of the conductor and wind tunnel test data, efficient and accurate galloping simulation is achieved.

[0006] The technical solution to achieve the purpose of this invention is as follows:

[0007] An equivalent loading method for ice-covered conductor galloping load based on finite element analysis includes the following steps:

[0008] A conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software.

[0009] The aerodynamic parameters of the three forces were obtained by wind tunnel test at full wind angle of attack, and an aerodynamic load mathematical model was constructed to calculate the load on the conductor.

[0010] The calculated load magnitude is then applied to the nodes in the conductor as a concentrated load.

[0011] In the preferred technical solution, establishing a conductor model based on actual line parameters includes:

[0012] Based on the actual transmission line parameter definitions, the physical parameters and cross-sectional characteristics corresponding to the conductor type, the conductor shape-finding method is used to calculate and create the nodes required for the line layout, connect the nodes, and form a conductor model.

[0013] In the preferred technical solution, the nonlinear finite element iterative method is used to find the shape of the conductor, and the equilibrium state of the conductor under its own weight and initial tension is calculated. The required nodes are generated by interpolation algorithm.

[0014] In the preferred technical solution, the equilibrium equation under its own weight is:

[0015]

[0016] Among them, the wire is regarded as having two fixed points ( , )and( , The flexible cable, z is the horizontal coordinate variable. and The horizontal coordinates of the fixed points at both ends of the flexible cable and The vertical coordinates of the fixed points at both ends of the flexible cable are: Let q be the target horizontal tension, q be the mass per unit length, and C be determined by the boundary conditions. Sure.

[0017] In the preferred technical solution, establishing the conductor and applying initial tension includes:

[0018] In a finite element platform, a two-way bar element subjected to only tension / compression is defined, material properties are set, and running tension is applied using the initial strain method to simulate the initial prestress state of the conductor. , Let E be the target horizontal tension, E be the elastic modulus, and A be the conductor area.

[0019] In the preferred technical solution, the construction of the aerodynamic load mathematical model includes:

[0020] Wind tunnel tests were conducted on icing conductors at wind angles of attack from 0° to 360° to obtain the three-component force coefficients and drag coefficients. Lift coefficient Torque coefficient A mathematical model of aerodynamic loads is established based on the curve of angle of attack variation:

[0021] , ，

[0022] in, As resistance, For lift, For torque, ρ air density, V For wind speed, D This is the equivalent diameter.

[0023] In the preferred technical solution, during the loading process of the concentrated load at the nodes in the conductor, quasi-static loading is performed by using the extreme value statistical method to equivalently apply the pulsating wind load. The time-varying fluctuation of the torque is averaged, and its mean value approaches zero. The vertical or lateral displacement of the conductor galloping is then calculated.

[0024] In the preferred technical solution, applying concentrated loads to the nodes in the conductor includes:

[0025] Based on the equivalent loading method of engineering mechanics, the concentrated force load at the nodes at both ends of the conductor is multiplied by... The element length is the equivalent concentrated force load at the node in the middle of the conductor multiplied by one element length;

[0026] The magnitudes of the concentrated lateral and vertical forces at the nodes at both ends of the conductor are:

[0027] ,

[0028] The magnitudes of the concentrated forces in the transverse and vertical directions of the concentrated load at the intermediate node of the conductor are:

[0029] ,

[0030] in, , Where N is the total length of the conductor and N is the number of units.

[0031] This invention also discloses an equivalent loading system for galloping loads of icy conductors based on finite element analysis, comprising:

[0032] The conductor model building module establishes a conductor model based on actual line parameters and uses engineering simulation software to create the conductor and apply initial tension.

[0033] The aerodynamic load model building module obtains the three-part force aerodynamic parameters through a wind tunnel test at full wind angle of attack, builds an aerodynamic load mathematical model, and calculates the load on the conductor.

[0034] The load loading module applies the calculated load to the nodes in the conductor, creating a concentrated load.

[0035] The present invention also discloses a computer storage medium storing a computer program, characterized in that when the computer executes the computer program, it implements the equivalent loading method for ice-covered conductor galloping load based on finite element analysis.

[0036] The present invention also discloses an electronic device, including a memory and a processor, characterized in that the memory stores a computer program, the processor runs the computer program stored in the memory, and when the computer program is executed, it implements the above-mentioned equivalent loading method for ice-covered conductor galloping load based on finite element analysis.

[0037] Compared with the prior art, the significant advantages of this invention are:

[0038] This invention transforms the dynamic fluid-structure interaction model into a quasi-static load model, avoiding the complex iterative calculations of the fluid and structural domains in traditional methods. This significantly reduces computational resource consumption and simulation time, meeting the timeliness requirements of engineering design. An aerodynamic load mathematical model is constructed based on the three-part force coefficients from wind tunnel tests. Combined with conductor nonlinear form-finding technology, the initial tension distribution and spatial curvature characteristics of the conductor are accurately reproduced, ensuring the physical realism of the aerodynamic loading and the accuracy of the dynamic response. An equivalent pulsating wind load is represented using the extreme value statistical method, ignoring secondary factors such as the near-zero mean of torque fluctuations and focusing on the main controlling factors of vertical or lateral displacement. This simplifies the model complexity while ensuring the reliability of the results. Through the equivalent loading strategy of concentrated forces at nodes, combined with time history comparison verification of key displacements, an efficient and reproducible simulation method is provided for the anti-galloping design of transmission lines, assisting in the optimization of disaster prevention measures and operation and maintenance decisions.

[0039] This invention significantly reduces computation time, controls errors to within 10%, balances computational efficiency and accuracy, and also ensures the authenticity of conductor spatial morphology and initial tension distribution. It provides an innovative solution for the analysis of galloping of icy conductors and has significant engineering application value. Attached Figure Description

[0040] Figure 1 The flowchart shows the equivalent loading method for ice-covered conductor galloping load based on finite element analysis.

[0041] Figure 2A graph showing the relationship between the three-component force coefficients and the angle of attack.

[0042] Figure 3 Find the shape diagram for the conductor;

[0043] Figure 4 The initial tension is 22.03 kN, and the lateral displacement at the mid-span node is [value missing].

[0044] Figure 5 The initial tension is 22.03 kN, and the vertical displacement at the mid-span node is [value missing].

[0045] Figure 6 The initial tension is 24 kN, representing the lateral displacement at the mid-span node.

[0046] Figure 7 The initial tension is 24 kN, and the vertical displacement at the mid-span node is given. Detailed Implementation

[0047] The principle of this invention is as follows: In view of the fact that existing simulation methods cannot take fluid-structure interaction into account when load is applied to conductors, this invention integrates nonlinear form finding of conductors, quasi-static equivalent of aerodynamic forces and concentrated force loading at nodes, thereby significantly improving simulation accuracy while ensuring computational efficiency.

[0048] Example:

[0049] like Figure 1 As shown, an equivalent loading method for icing conductor galloping load based on finite element analysis includes the following steps:

[0050] A conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software.

[0051] The aerodynamic parameters of the three forces were obtained by wind tunnel test at full wind angle of attack, and an aerodynamic load mathematical model was constructed to calculate the load on the conductor.

[0052] The calculated load magnitude is then applied to the nodes in the conductor as a concentrated load.

[0053] In a preferred embodiment, the engineering simulation software can be ANSYS or other simulation software, using the ANSYS LINK10 element model to simulate the transmission line conductors and establish conductor elements.

[0054] Establishing a conductor model based on actual line parameters includes:

[0055] Based on the actual transmission line parameter definitions, the physical parameters and cross-sectional characteristics corresponding to the conductor type, the conductor shape-finding method is used to calculate and create the nodes required for the line layout, connect the nodes, and form a conductor model.

[0056] Of course, other methods can also be used to establish the conductor model, which is not limited here.

[0057] In a preferred embodiment, a nonlinear finite element iterative method is used to find the shape of the conductor, calculate the equilibrium state of the conductor under its own weight and initial tension, and generate the required nodes through an interpolation algorithm.

[0058] Of course, other methods can also be used for conductor shape finding, and no specific method is specified here.

[0059] In a preferred embodiment, the equilibrium equation under its own weight is:

[0060]

[0061] Among them, the wire is regarded as having two fixed points ( , )and( , The flexible cable, z is the horizontal coordinate variable. and The horizontal coordinates of the fixed points at both ends of the flexible cable and The vertical coordinates of the fixed points at both ends of the flexible cable are: Let q be the target horizontal tension, q be the mass per unit length, and C be determined by the boundary conditions. Sure.

[0062] In a preferred embodiment, establishing the conductor and applying initial tension includes:

[0063] In a finite element platform, a two-way bar element subjected to only tension / compression is defined, material properties are set, and running tension is applied using the initial strain method to simulate the initial prestress state of the conductor. , Let E be the target horizontal tension, E be the elastic modulus, and A be the conductor area.

[0064] In a preferred embodiment, constructing the aerodynamic load mathematical model includes:

[0065] Wind tunnel tests were conducted on icing conductors at wind angles of attack from 0° to 360° to obtain the three-component force coefficients and drag coefficients. Lift coefficient Torque coefficient A mathematical model of aerodynamic loads is established based on the curve of angle of attack variation:

[0066] , ，

[0067] in, As resistance, For lift, For torque, ρ air density,V For wind speed, D This is the equivalent diameter.

[0068] In a preferred embodiment, during the loading of concentrated loads onto nodes in the conductor, quasi-static loading is performed using an equivalent pulsating wind load through extreme value statistics. The time-varying fluctuations of torque are averaged out, and their mean value approaches zero. The vertical or lateral displacement of the conductor galloping is then calculated.

[0069] In a preferred embodiment, applying concentrated loads to the nodes in the conductor includes:

[0070] Based on the equivalent loading method of engineering mechanics, the concentrated force load at the nodes at both ends of the conductor is multiplied by... The element length is the equivalent concentrated force load at the node in the middle of the conductor multiplied by one element length;

[0071] The magnitudes of the concentrated lateral and vertical forces at the nodes at both ends of the conductor are:

[0072] ,

[0073] The magnitudes of the concentrated forces in the transverse and vertical directions of the concentrated load at the intermediate node of the conductor are:

[0074] ,

[0075] in, , Where N is the total length of the conductor and N is the number of units.

[0076] Another embodiment, an equivalent loading system for icing conductor galloping load based on finite element analysis, includes:

[0077] The conductor model building module establishes a conductor model based on actual line parameters and uses engineering simulation software to create the conductor and apply initial tension.

[0078] The aerodynamic load model building module obtains the three-part force aerodynamic parameters through a wind tunnel test at full wind angle of attack, builds an aerodynamic load mathematical model, and calculates the load on the conductor.

[0079] The load loading module applies the calculated load to the nodes in the conductor, creating a concentrated load.

[0080] The following example illustrates the workflow of the equivalent loading system for icing conductor galloping load based on finite element analysis, including the following steps:

[0081] (1) Based on the actual transmission line parameter definition, the physical parameters and cross-sectional characteristics corresponding to the conductor type, the key nodes required for the line layout are calculated and created using the conductor shape finding method, and then the nodes are connected to form a conductor model;

[0082] (2) Using the LINK10 element model of ANSYS, the conductor of the transmission line is simulated. The conductor element is established by using preset nodes, material properties and cross-sectional characteristics, and the initial running tension is applied to each conductor.

[0083] (3) Obtain the three-part aerodynamic parameters of the conductor under the full wind angle of attack by wind tunnel test of the aerodynamic characteristics of ice splitting, extract the aerodynamic drag coefficient and aerodynamic lift coefficient, and use the aerodynamic load mathematical model to calculate and solve the load applied on the conductor.

[0084] (4) The calculated load magnitude is converted into a static quasi-static load model by mechanical loading method and concentrated load loading method, and applied at each node of the conductor.

[0085] (5) After the specified loading time is completed, the displacement results of the specified nodes are extracted to simulate and reproduce the galloping process of the transmission line. The simulated displacement results are compared with the experimental results obtained by the original fluid-structure interaction iterative method to verify the accuracy and rationality of the equivalent method.

[0086] In step (1), the accurate modeling based on the actual transmission line parameters (conductor type, physical parameters, cross-sectional characteristics) inputs the physical parameters and cross-sectional characteristics (horizontal tension) corresponding to the conductor type. ,density The conductor diameter (D) and conductor area (A) are determined using a nonlinear finite element iterative method for conductor shape finding. The equilibrium state of the conductor under its own weight and initial tension is calculated. Key nodes (suspension points, lowest points) are generated using an interpolation algorithm, and these nodes are connected to form a spatial curve model of the conductor. Fixed points at both ends are considered. , )and( , The equilibrium equation for a flexible cable under its own weight is:

[0087] (1)

[0088] Where z is the horizontal coordinate variable, and The horizontal coordinates of the fixed points at both ends of the flexible cable and The vertical coordinates of the fixed points at both ends of the flexible cable are: Let q be the target horizontal tension, q be the mass per unit length, and C be determined by the boundary conditions. This can be determined by solving the following equations:

[0089] (2)

[0090] in This represents the span length.

[0091] In step (2), define LINK10 elements (tension / compression bidirectional bar elements only) in finite element platforms such as ANSYS and set material properties (elastic modulus E, Poisson's ratio). The initial prestress state of the conductor is simulated by applying operating tension using the initial strain method. The initial strain of the conductor is taken as 0.8 times the stress, i.e.

[0092] (3)

[0093] In step (3), a wind tunnel test with an angle of attack of 0° to 360° is conducted on the icing conductor. The three-component force coefficients are obtained by referring to existing papers on fluid-structure interaction research on aerodynamic coefficients. The drag coefficient, The lift coefficient, Using the torque coefficient as a function of angle of attack, an aerodynamic load formula is established:

[0094] , ， (4)

[0095] in ρ For air density, take 1.25 kg / m³. V represents wind speed, and D represents the equivalent diameter. During the load application process, the quasi-static loading has been equivalent to the pulsating wind load using the extreme value statistical method. The time-varying fluctuations of the torque are averaged out, and its mean approaches zero, resulting in a small overall contribution. Focusing on the vertical or lateral displacement of the conductor galloping, and assuming that torsion does not induce significant coupled vibrations, simplifies the treatment.

[0096] In step (4), the equivalent concentrated force loading of the nodes is considered. Therefore, the equivalent loading method of engineering mechanics is used, multiplying the concentrated force load at the nodes at both ends of the conductor by... The element length is calculated by multiplying the equivalent concentrated force load at the node in the middle of the conductor by one element length, i.e., assuming... Let N be the total length of the conductor and N be the number of units. Then let... Then we have:

[0097] , (5)

[0098] This refers to the magnitude of the concentrated lateral and vertical forces at the two ends of the conductor.

[0099] , (6)

[0100] This refers to the magnitude of the concentrated load in the horizontal and vertical directions at the intermediate node of the conductor.

[0101] In step (5), the displacement results can take into account the displacement time history of key nodes, with the mid-span node as the research focus, and consider vertical or lateral displacement.

[0102] To better illustrate the effectiveness of the method for equivalent loading of galloping loads using ANSYS software provided in this invention, a single-span conductor in a certain project is used as an example for comparison with data obtained through traditional fluid-structure interaction. This method is called the equivalent nodal concentrated load method. The PC configuration for galloping calculations in this invention is as follows: CPU Intel Core i5-1130G7@4.00GHz quad-core, memory 16GB LPDDR4x4266MHz, hard disk 512GB, motherboard Intel Tiger Lake-UP3.

[0103] LGJ300 / 25 steel-cored aluminum stranded wire was selected as the conductor to consider the rationality of mid-span node loading. The variation of the aerodynamic force coefficients of the D-shaped icing-resistant split conductor with wind angle of attack was cited from the relevant paper. Figure 2 As shown, the initial strain of the icing conductor is taken as 0.8 times the stress. The span of a single conductor is 126.77m, the height difference between the two supports is 0m, and the physical parameters of the sub-conductors are shown in Table 1:

[0104] Table 1 Physical Parameters of Icy Sub-conductors

[0105]

[0106] This study selects the maximum aerodynamic coefficient to investigate the equivalent pulsating wind load using the extreme value statistical method, and applies it to the concentrated load. The finite element model diagram and form-finding model for an icing single-span conductor are shown below. Figure 3 As shown in Table 2, the coordinates of the key nodes for shape finding are as follows:

[0107] Table 2 Coordinates of Key Nodes for Shape Finding

[0108]

[0109] Figure 4 and Figure 5This indicates that when considering the equivalent loading of 22.03 kN initial tension in LGJ300 / 25 steel-cored aluminum stranded wire, the lumped mass matrix setting is enabled to activate the prestress effect in the large deformation nonlinear analysis of the conductor. The maximum aerodynamic coefficient is selected for calculation. Half of the concentrated load is applied at both ends of the conductor node, and concentrated loads are applied at the nodes in the conductor. The processing time is 3 seconds, and the lateral and vertical displacements of the mid-span node are obtained. Figure 4 and Figure 5 As can be seen, the lateral amplitude is 0.0231m, slightly larger than the original data of 0.0211m, with an error of 9.47%, which is consistent with the equivalent loading of the lateral nodes; the vertical amplitude is 0.763m, slightly larger than 0.76m, with an error of 0.26%, which is consistent with the equivalent loading of the vertical nodes. The reason for the better equivalent effect of vertical displacement can be attributed to the fact that the conductor usually has higher bending stiffness in the vertical direction due to its own weight; while the lateral stiffness is relatively low, making it more sensitive to lateral displacement under the same loading conditions, and small errors are easily amplified.

[0110] To consider the influence of initial tension on equivalent loading, which affects the sag of the conductor, and to study the impact of conductor sag on galloping, eight initial tension conditions ranging from 10kN to 24kN were selected, with other conditions remaining constant. The galloping amplitude under each condition was calculated. The simulation used 24kN to verify the effectiveness of the equivalent loading method under varying initial tension. The experimental model parameters remained the same: wind speed 4.1m / s, conductor diameter 2.23cm, icing thickness 2.06cm, and span 126.77m. Figure 6 and Figure 7 This indicates that when considering the equivalent load of an initial tension of 24kN for LGJ300 / 25 steel-cored aluminum stranded wire, a concentrated load is applied to the nodes in the conductor, with a processing time of 3 seconds, to obtain the lateral and vertical displacements of the mid-span nodes. Figure 6 and Figure 7 As can be seen, the lateral amplitude is 0.0333m, slightly smaller than the original data of 0.0356m, with an error of 6.46%, which is consistent with the equivalent loading of the lateral node; the vertical amplitude is 0.709m, slightly smaller than 0.72m, with an error of 1.52%, which is consistent with the equivalent loading of the vertical node.

[0111] Simulation analysis shows that the simulation time obtained by fluid-structure interaction iterative conductor galloping is 4000s, but the equivalent method of applying concentrated loads through APDL significantly reduces the time compared to fluid-structure interaction.

[0112] In summary, the equivalent loading method for ice-covered conductor galloping based on finite element analysis proposed in this invention overcomes the incompatibility of traditional fluid-structure interaction load loading methods for studying ice-covered conductor galloping. It seeks a simple and rapid load loading method, achieving the effect of easy implementation, wide applicability, and convenient simulation of equivalent load loading for transmission line galloping based on this simulation through concentrated load loading.

[0113] In another embodiment, a computer storage medium stores a computer program that, when executed, implements the above-described equivalent loading method for icing conductor galloping load based on finite element analysis.

[0114] The specific implementation method is the same as described above, and will not be repeated here.

[0115] In another embodiment, an electronic device includes a memory and a processor. The memory stores the aforementioned computer program, and the processor runs the computer program stored in the memory. When the computer program is executed, it implements the aforementioned equivalent loading method for icing conductor galloping load based on finite element analysis.

[0116] The specific implementation method is the same as described above, and will not be repeated here.

[0117] The above embodiments are preferred embodiments of the present invention, but the embodiments of the present invention are not limited to the above embodiments. Any changes, modifications, substitutions, combinations, or simplifications made without departing from the spirit and principle of the present invention shall be considered equivalent substitutions and shall be included within the protection scope of the present invention.

Claims

1. An ice-coated conductor galloping load equivalent loading method based on finite element analysis, characterized in that, The method comprises the following steps: A conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software; Three-component aerodynamic parameters are obtained through full-attack-angle wind tunnel tests, an aerodynamic load mathematical model is constructed, and loads on the conductor are calculated; The construction of the aerodynamic load mathematical model comprises: The wind tunnel test of 0°-360° wind attack angle is carried out on the iced conductor to obtain the three-component force coefficients, drag coefficient , lift coefficient , and torque coefficient , and the mathematical model of aerodynamic load is established according to the curves of the coefficients changing with attack angle. ， ， wherein, is the drag force, is the lift force, is the torque, ρ is the air density, V is the wind speed, D is the equivalent diameter; The calculated load is applied to the nodes in the conductor as concentrated loads, and the calculated load is applied to the nodes in the conductor as concentrated loads comprise: The method of equivalent loading based on engineering mechanics, the concentrated force load at the nodes at both ends of the wire is multiplied by the unit length, and the equivalent concentrated force load at the nodes in the middle of the wire is multiplied by a unit length.

2. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 1, wherein, The conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software. The conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software.

3. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 2, wherein, The conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software. The conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software.

4. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 3, wherein, The conductor model is established based on actual line parameters, and the conductor is established and initial tension is applied based on engineering simulation software. where the wire is considered as a flexible cable with two fixed points at , ) and ( , ), z is the horizontal coordinate variable, and the horizontal coordinates of the fixed points of the flexible cable, and the vertical coordinates of the fixed points of the flexible cable, is the target horizontal tension, q is the mass per unit length, and C is determined by the boundary conditions .

5. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 1, wherein, The equation of the equilibrium state of the conductor under the action of the self weight is: The two-way bar element subjected to tension / compression only is defined in the finite element platform, material properties are set, initial tension is applied by initial strain method to simulate the initial prestressed state of the conductor, and the initial tension of the conductor , is the target horizontal tension, E is the elastic modulus, and A is the conductor area.

6. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 1, wherein, The conductor is established and initial tension is applied, and the conductor model is established based on actual line parameters.

7. The icing conductor galloping load equivalent loading method based on finite element analysis of claim 1, wherein, In the process of loading the calculated load on the nodes in the conductor as concentrated loads, the equivalent fluctuating wind load is loaded as a quasi-static load through the extreme value statistical method, and the time-varying fluctuation of the torque is homogenized, and the mean value tends to be zero, so that the vertical or lateral displacement of the conductor galloping is calculated. ， The size of the horizontal and vertical concentrated forces of the concentrated loads of the nodes at both ends of the conductor is: ， wherein, , is the total length of the wire, N is the number of units, The size of the horizontal and vertical concentrated forces of the concentrated loads of the nodes at both ends of the conductor is: is the air density, V is the wind speed, D is the equivalent diameter, is the drag coefficient, is the lift coefficient.

8. An icing conductor galloping load equivalent loading system based on finite element analysis, characterized in that, The size of the horizontal and vertical concentrated forces of the concentrated loads of the nodes at both ends of the conductor is: The size of the horizontal and vertical concentrated forces of the concentrated loads of the nodes at both ends of the conductor is: The conductor model construction module, which establishes a conductor model based on actual line parameters and establishes a conductor and applies initial tension based on engineering simulation software; The aerodynamic load model construction module, which obtains three-component aerodynamic parameters through full-attack-angle wind tunnel tests, constructs an aerodynamic load mathematical model, and calculates loads on the conductor; The wind tunnel test of 0°-360° wind attack angle is carried out on the iced conductor to obtain the three-component force coefficients, drag coefficient , lift coefficient , and torque coefficient , and the mathematical model of aerodynamic load is established according to the curves of the coefficients changing with attack angle. ， ， wherein, is the drag force, is the lift force, is the torque, The construction of the aerodynamic load mathematical model comprises: is the air density, V is the wind speed, D is the equivalent diameter; The size of the horizontal and vertical concentrated forces of the concentrated loads of the nodes at both ends of the conductor is: The method of equivalent loading based on engineering mechanics, the concentrated force load at the nodes at both ends of the wire is multiplied by the unit length, and the equivalent concentrated force load at the nodes in the middle of the wire is multiplied by a unit length.

9. A computer storage medium having stored thereon a computer program, characterized in that The computer executes the computer program to realize the equivalent loading method of the icing conductor galloping load based on the finite element analysis in any one of claims 1-7.

10. An electronic device comprising a memory and a processor, characterized in that, The computer executes the computer program to realize the equivalent loading method of the icing conductor galloping load based on the finite element analysis in any one of claims 1-7.

Images