High-voltage line wind vibration simulation analysis method, system and equipment based on Modelica language and medium

Abstract

The invention relates to high-voltage line wind vibration analysis, in particular to a high-voltage line wind vibration simulation analysis method, system and device based on the Modelica language and a medium. Constructing a parameterized sub-model of the high-voltage line system based on the physical parameters; packaging each parameterized sub-model into a model assembly by using a Modelica language, integrating the corresponding model assembly based on an actual topological structure of the high-voltage line system, defining an initial state of the system, and establishing a kinetic equation for each cable section to form a high-voltage line system model; acquiring wind vector data of the current environment in real time; inputting the wind vector data of the current environment into the high-voltage line system model for simulation analysis, solving the kinetic equation of each cable section through numerical integration, and calculating the movement speed of each cable section at a specified moment in the future and the position of the midpoint of the cable section. According to the invention, high-precision and high-efficiency multi-physical-field wind vibration simulation can be realized.

Description

TECHNICAL FIELD

[0001] The application relates to high-voltage line wind vibration analysis, in particular to a high-voltage line wind vibration simulation analysis method, system, device and medium based on a Modelica language. BACKGROUND

[0002] As a key component of the backbone network of the power system, the high-voltage line undertakes the important mission of large-capacity and long-distance power transmission, and the safe and stable operation of the high-voltage line has an increasingly prominent impact on the reliability of the entire power grid. The high-voltage line is exposed to the natural environment for a long time and is easily affected by the continuous action of wind load, which may cause wind vibration phenomena such as wind deviation and dancing. These vibrations not only cause conductor fatigue damage and fitting wear, but also may cause serious accidents such as inter-phase flashover, wire breakage and tower collapse, which directly threaten the safe operation of the power grid. Therefore, it is of great significance to carry out high-voltage line wind vibration analysis to ensure the safety of the power grid and prevent accidents.

[0003] At present, for the analysis of high-voltage line wind vibration, the engineering field mainly uses finite element analysis software to establish a refined model for simulation calculation, and simulates the dynamic response of the line under wind load through meshing and numerical iteration. In addition, some studies also use simplified methods such as rigid straight rod method or statics single pendulum model to achieve rapid calculation and evaluation by appropriately simplifying the line system.

[0004] However, the existing simulation analysis methods for high-voltage lines have obvious deficiencies: the traditional simulation method often simplifies the line system too much, ignoring key dynamic characteristics such as the elastic deformation of the conductor and the time-varying effect of aerodynamic force, resulting in limited simulation accuracy; the finite element-based analysis method has high accuracy, but requires a large amount of meshing and iterative calculation, and has low calculation efficiency, making it difficult to meet the rapid evaluation needs of multiple working conditions and multiple scenarios; in addition, the existing simulation tools have limited capabilities in multi-physical field coupling, making it difficult to efficiently realize model integration and collaborative simulation in different disciplines such as aerodynamics, multi-body dynamics and structural mechanics. SUMMARY

[0005] In view of the technical problems that the existing simulation analysis methods for high-voltage lines have the problem of excessive simplification of traditional methods, resulting in the lack of dynamic characteristics, the low calculation efficiency of finite element analysis methods, making it difficult to meet the rapid evaluation of multiple working conditions, and the difficulty of multi-physical field coupling, the application provides a high-voltage line wind vibration simulation analysis method, system, device and medium based on a Modelica language, which simulates the dynamic characteristics of the cable by using a parameterized sub-model after discretizing the cable, encapsulates components based on Modelica and constructs a high-voltage line system model, and performs simulation analysis by using a numerical integration method, realizing high-precision and high-efficiency multi-physical field wind vibration simulation.

[0006] In a first aspect, the application provides a high-voltage line wind vibration simulation analysis method based on Modelica language, comprising the following steps: S1. Obtain the physical parameters of the high-voltage line system, including cable parameters and tower parameters; S2. Construct a parameterized submodel of the high-voltage line system based on the physical parameters, comprising: Discretize the continuous cable into multiple cable segments; Construct a rigid body model, a gravity model, an elastic property model, and an aerodynamic property model of a single cable segment, and construct an equivalent tower model; The gravity model calculates the gravity acting on the cable segment through the mass of the single cable segment; The elastic property model calculates the elastic deformation variable based on the positions of the two ends of the cable segment and the preset length of the cable segment, and calculates the tension of the adjacent cable segment based on the elastic deformation variable; The aerodynamic property model calculates the aerodynamic force acting on the cable segment based on the current wind vector and the motion speed of the cable segment; The equivalent tower model calculates the tension of the adjacent cable segment on the tower based on the displacement of the tower node through the equivalent stiffness matrix; S3. Encapsulate each parameterized submodel into a model component using Modelica language, integrate the corresponding model components based on the actual topological structure of the high-voltage line system, define the initial state of the system, and establish the dynamic equation for each cable segment to form a high-voltage line system model; S4. Real-time acquisition of the current environmental wind vector data; S5. Input the current environmental wind vector data into the high-voltage line system model for simulation analysis, solve the dynamic equation of each cable segment through numerical integration, and calculate the motion speed of each cable segment and the position of the midpoint of the cable segment at a specified future time, wherein: The motion speed of each cable segment at the initial time is a preset value, and in each subsequent integration time step, the aerodynamic force acting on the cable segment is calculated based on the motion speed of the cable segment calculated in the previous time step, and then the dynamic equation of each cable segment is solved.

[0007] Further, in step S1, the cable parameters include the line span, span distribution, conductor cross-sectional area, elastic modulus, cable expansion coefficient, and unit length mass of the cable, and the tower parameters include the structural size, tower height, cross arm position, tower body segmented length, and material density of the tower.

[0008] Further, in step S2, the length of each cable segment is ≤50 times the diameter of the cable.

[0009] Further, in step S2, the gravity model of the first The model expression of the gravity model of the first cable segment is:

[0010] in Representing the The weight of each cable segment; Representing the The quality of each cable segment ; Mass per unit length of high-voltage line cable; Representing the The length of each cable segment; It represents gravitational acceleration.

[0011] It should be further explained that in step S2, the first The elastic characteristic model expression for each cable segment is:

[0012] in, Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; k To increase the tensile stiffness of the cable; Representing the The elastic deformation of a cable segment; b The cable damping coefficient; represent The rate of change; represent The unit direction; Each cable segment's elasticity model is equipped with a left-side interface and a right-side interface.

[0013] It should be further explained that, The calculation method is as follows: Calculate the current cable segment Left adjacent cable segment to current cable segment When applying tension, the left interface is used to input the real-time position of the adjacent cable segment to the left of the current cable segment. The interface on the right is used to input the real-time position of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the left of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to its left. ; Calculate the current cable segment Right adjacent cable segment to current cable segment When tension is applied, the left interface is used to input the real-time position of the current cable segment. The right-hand interface is used to input the real-time position of the cable segment adjacent to the right of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the right of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to the right of the current cable segment. .

[0014] It should be further explained that in step S2, the first The aerodynamic characteristic model expression for each cable segment is as follows:

[0015] in, For the first The aerodynamic forces acting on each cable segment; ρ Represents air density; Represents the average wind vector in the ground coordinate system; Represents the normal force coefficient; To be Transformation matrix to the cable coordinate system; S This refers to the cross-sectional area of ​​the cable. For the first The speed of movement of each cable segment; For the first The relative speed of motion between the cable segment and the air .

[0016] It should be further explained that the method for constructing the equivalent model of the iron tower in step S2 is as follows: For each iron tower, the tower body is divided into several rectangular cross-section beams according to the mechanical force transmission path of the iron tower; A parameter calculation module is set up for each I-beam cross-section to calculate the moment of inertia of each I-beam cross-section. Cross-section beam quality and cross-section beam stiffness The expressions include:

[0017] Where A and B are the outer and inner side lengths of the rectangular section of the beam; The length of the cross-section beam; Density of the steel in the cross-section beam; The elastic modulus of the beam material; Based on the parameter calculation module for each cross-section beam, a preliminary equivalent model of each iron tower is built using the beam component of the Modelica model library; Modal analysis was performed on the preliminary model of each tower to calculate its first-order bending frequency. Then compare The first-order bending frequency of an actual iron tower f 1: like This reduces the stiffness of the beam section. and moment of inertia of cross section ; like This increases the stiffness of the beam section. and moment of inertia of cross section ; Then repeat the above calculations until the first-order bending frequency of the equivalent model of the tower is reached. Compared to the actual iron tower f1 If the error is ≤5%, the equivalent model of the corresponding iron tower is obtained.

[0018] It should be further noted that in step S3, the first The dynamic equations for each cable segment are expressed as follows:

[0019] in, Representing the The quality of each cable segment; Representing the The acceleration of each cable segment; Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; Representative and the All adjacent cable segments of the first cable segment are related to the second cable segment. The sum of the tensile forces on each cable segment; Representative and the The tower connected to the cable segment is the first The tensile force on each cable segment; Representing the The weight of each cable segment; Representing the The aerodynamic forces acting on each cable segment.

[0020] It should be further noted that step S3 also includes: encapsulating the high-voltage line system model into a dynamic link library based on the FMI standard to support co-simulation with third-party simulation tools, specifically including: The high-voltage line system model is compiled into C code, and dynamic link libraries, output state variables, input interfaces and exported functions are generated according to the FMI2.0 standard. It enables co-simulation with third-party simulation tools through the FMI interface, supporting collaborative analysis of multiphysics fields. Third-party simulation tools include Simulink and Dymola.

[0021] It should be further noted that in step S5, a CVODE variable-order, variable-step integrator is used to solve the dynamic equations for each cable segment. Combined with high-stiffness damping and zero-crossing event shielding technology, computational efficiency is improved. Specifically, this includes: CVODE variable order and step size integration: The CVODE solver is used to dynamically adjust the integration order and step size according to the local truncation error, balancing computational accuracy and efficiency. High stiffness damping treatment: Introducing a large damping coefficient into the elastic property model ( b =104~105Ns / m), suppressing the influence of high-frequency longitudinal vibration on the low-frequency response of wind vibration, and avoiding the solver having too small a step size due to "illness" problem; Zero-crossing event shielding: To address the discontinuity issue when the cable stretch crosses zero (ΔL=0), Modelica's noEvent mechanism is used to skip the precise moment search, retaining only the critical time step of low-frequency wind vibration, thereby improving calculation speed.

[0022] Secondly, this application provides a high-voltage line wind vibration simulation analysis system based on the Modelica language, used to implement the above-mentioned high-voltage line wind vibration simulation analysis method, including: The physical parameter acquisition module is used to acquire the physical parameters of the high-voltage line system, including cable parameters and tower parameters; The parameterized sub-model construction module is used to construct parameterized sub-models of high-voltage line systems based on physical parameters. This includes discretizing continuous cables into multiple cable segments, constructing rigid body models, gravity models, elastic property models, and aerodynamic property models for individual cable segments, and constructing equivalent models of towers. The high-voltage line system model building module is used to encapsulate various parameterized sub-models into model components using the Modelica language. Based on the actual topology of the high-voltage line system, it integrates the corresponding model components, defines the initial state of the system, and establishes dynamic equations for each cable segment to form a high-voltage line system model. The wind vector data acquisition module is used to acquire wind vector data of the current environment in real time. The simulation analysis module is used to input the wind vector data of the current environment into the high-voltage line system model for simulation analysis. It solves the dynamic equation of each cable segment through numerical integration and calculates the velocity of each cable segment and the position of the midpoint of the cable segment at a specified time in the future.

[0023] Thirdly, this application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps of the above-mentioned high-voltage line wind vibration simulation analysis method.

[0024] Fourthly, this application provides a storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described high-voltage line wind vibration simulation analysis method.

[0025] As can be seen from the above technical solutions, this application has the following advantages: 1. This application discretizes a continuous cable into multiple cable segments and constructs a parametric sub-model that includes a rigid body model, a gravity model, an elastic characteristic model, an aerodynamic characteristic model, and an equivalent tower model. This enables a more refined simulation of the dynamic characteristics and aerodynamic response of the cable, effectively overcoming the problem of oversimplification in existing models and significantly improving the accuracy and reliability of wind vibration simulation.

[0026] 2. This application encapsulates each parametric sub-model into a standardized model component using the Modelica language, integrates them based on the actual topology of the high-voltage line system, and establishes dynamic equations for each cable segment. This achieves unified description and efficient integration of multi-disciplinary models, solves the problem of multi-physics coupling in existing technologies, and enhances the flexibility and applicability of the simulation system.

[0027] 3. This application acquires wind field data in real time and inputs it into the system model for simulation analysis. It uses the numerical integration method to solve the dynamic equation of the cable segment and combines it with an efficient solution strategy, which greatly improves the calculation efficiency and overcomes the limitation of the traditional finite element method in terms of calculation time. It can meet the rapid simulation needs of multiple scenarios and multiple working conditions in engineering practice. Attached Figure Description

[0028] To more clearly illustrate the technical solution of this application, the accompanying drawings used in the description will be briefly introduced below. Obviously, the accompanying 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.

[0029] Figure 1 This is a flowchart of a high-voltage line wind vibration simulation analysis method based on the Modelica language in one embodiment of this application.

[0030] Figure 2 This is the first embodiment of the present invention. A schematic diagram of the stress analysis of each cable segment.

[0031] Figure 3 This is a schematic diagram of the X-direction displacement simulation curve in the model verification of this invention.

[0032] Figure 4 This is a schematic diagram of the Y-direction displacement simulation curve during the model verification of this invention.

[0033] Figure 5 This is a schematic diagram of the Z-direction displacement simulation curve in the model verification of this invention.

[0034] Figure 6 This is a schematic flowchart illustrating the simulation curve of the synchronous change of displacement in the xz direction of the high-voltage line after gravity balance in the model verification of this invention.

[0035] Figure 7 This is a schematic block diagram of a high-voltage line wind vibration simulation analysis system based on the Modelica language in one embodiment of this application.

[0036] Figure 8 This is a schematic diagram of the hardware structure of an electronic device in one embodiment of this application. Detailed Implementation

[0037] To make the purpose, features, and advantages of this application more apparent and understandable, specific embodiments and accompanying drawings will be used to clearly and completely describe the technical solution protected by this application. Obviously, the embodiments described below are only some embodiments of this application, and not all embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0038] The following will describe in detail the high-voltage line wind vibration simulation analysis method involved in this application. Specific details such as particular system structures and technologies are presented for illustrative purposes rather than limiting, in order to provide a thorough understanding of the embodiments of this application. However, those skilled in the art will understand that this application can also be implemented in other embodiments without these specific details.

[0039] In the high-voltage line wind vibration simulation analysis method involved in this application, the term "comprising" indicates the presence of the described feature, whole, step, operation, element and / or component, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or collections thereof. The terms "comprising," "including," "having," and variations thereof all mean "including but not limited to," unless otherwise specifically emphasized.

[0040] To facilitate a clear description of the technical solutions of this application, the terms "first" and "second" are used to distinguish identical or similar items with essentially the same function and effect. Those skilled in the art will understand that the terms "first" and "second" do not limit the quantity or execution order, and that the terms "first" and "second" do not necessarily imply that they are different.

[0041] The terms "one embodiment" or "some embodiments" used in this application mean that one or more embodiments of this application include the specific features, structures, or characteristics described in that embodiment. Therefore, the phrases "in one embodiment," "in some embodiments," "in other embodiments," "in still other embodiments," etc., appearing in different parts of this application do not necessarily refer to the same embodiment, but rather mean "one or more, but not all, embodiments," unless otherwise specifically emphasized.

[0042] The technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings.

[0043] The high-voltage line wind vibration simulation analysis method provided in this application embodiment is executed by computer equipment. Correspondingly, the high-voltage line wind vibration simulation analysis system based on the Modelica language runs on the computer equipment.

[0044] Figure 1 This is a flowchart of a high-voltage line wind vibration simulation analysis method based on the Modelica language, according to one embodiment of this application. Figure 1 The implementing entity can be a high-voltage line wind vibration simulation analysis system. Depending on different requirements, the order of steps in this flowchart can be changed, and some steps can be omitted.

[0045] like Figure 1As shown, the high-voltage line wind vibration simulation analysis method based on Modelica language includes: Step S1: Obtain the physical parameters of the high-voltage line system, including cable parameters and tower parameters.

[0046] By obtaining detailed physical parameters of the high-voltage line system, an accurate and comprehensive data foundation is provided for the subsequent construction of a high-fidelity parametric simulation model, ensuring the authenticity and reliability of the model input conditions. In some specific embodiments, cable parameters include cable span, span distribution, conductor cross-sectional area, elastic modulus, coefficient of linear expansion, and mass per unit length; tower parameters include tower structural dimensions, tower height, crossarm position, tower segment length, and material density.

[0047] By clearly defining cable parameters, including line span, span distribution, conductor cross-sectional area, elastic modulus, coefficient of linear expansion, and mass per unit length, as well as tower parameters, including structural dimensions, tower height, crossarm position, tower segment length, and material density, the physical parameters obtained in step S1 can comprehensively cover the key factors affecting the mechanical and wind-induced vibration behavior of the line, providing complete parameter definitions for building a high-precision model.

[0048] Step S2, constructing a parameterized sub-model of the high-voltage line system based on physical parameters, including: Discretize the continuous cable into multiple cable segments; Construct rigid body models, gravity models, elastic property models, and aerodynamic property models for a single cable segment, as well as construct an equivalent model of the tower; The gravity model calculates the gravity acting on a cable segment by using the mass of a single cable segment. The elastic characteristic model calculates the elastic deformation by using the positions of both ends of the cable segment and the preset length of the cable segment, and calculates the tension of adjacent cable segments based on the elastic deformation. The aerodynamic characteristic model calculates the aerodynamic forces acting on the cable segment based on the current wind vector and the cable segment's velocity. The equivalent model of the tower is based on the displacement of the tower nodes, and the tension of the tower on the adjacent cable segments is calculated by the equivalent stiffness matrix.

[0049] By discretizing the continuous cable into multiple cable segments and constructing rigid body models, gravity models, elastic characteristic models, and aerodynamic characteristic models for each individual cable segment, as well as constructing an equivalent model of the tower, a refined mathematical description of the line mass distribution, elastic deformation, aerodynamic load, and tower boundary constraints is achieved, laying a model foundation for accurately simulating the complex dynamic behavior of high-voltage lines under wind load.

[0050] In some specific embodiments, the length of each cable segment is ≤50 times the cable diameter.

[0051] By limiting the length of each cable segment in step S2 to no more than 50 times the cable diameter, the length of the discretized cable segments is ensured to be small enough. This ensures computational efficiency while accurately capturing the mechanical properties and aerodynamic loads distributed along the line, effectively balancing model accuracy and computational cost.

[0052] In some specific embodiments, the first The gravity model expression for each cable segment is:

[0053] in Representing the The weight of each cable segment; Representing the The quality of each cable segment ; Mass per unit length of high-voltage line cable; Representing the The length of each cable segment; It represents gravitational acceleration.

[0054] By providing the first The explicit mathematical expression of the gravity model of a cable segment, that is, the gravity acting on it is directly calculated based on the mass and gravitational acceleration of the cable segment, realizes the accurate quantification of the basic and constant load of the line self-weight, and provides a reliable static load input for solving the dynamic equation.

[0055] In some specific embodiments, the first The elastic characteristic model expression for each cable segment is:

[0056] in, Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; k To increase the tensile stiffness of the cable; Representing the The elastic deformation of a cable segment; b The cable damping coefficient; represent The rate of change; represent The unit direction; Each cable segment's elasticity model is equipped with a left-side interface and a right-side interface.

[0057] By providing the first The explicit mathematical expression of the elastic characteristic model of each cable segment is obtained. Based on the elastic deformation and its rate of change of the cable segment, combined with the cable tensile stiffness and damping coefficient, the tensile force between adjacent cable segments is calculated. With the configuration of left and right side interfaces, the model achieves an accurate description and interconnection of the elastic interaction and energy dissipation mechanism inside the cable.

[0058] In some specific embodiments, The calculation method is as follows: Calculate the current cable segment Left adjacent cable segment to current cable segment When applying tension, the left interface is used to input the real-time position of the adjacent cable segment to the left of the current cable segment. The interface on the right is used to input the real-time position of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the left of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to its left. ; Calculate the current cable segment Right adjacent cable segment to current cable segment When tension is applied, the left interface is used to input the real-time position of the current cable segment. The right-hand interface is used to input the real-time position of the cable segment adjacent to the right of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the right of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to the right of the current cable segment. .

[0059] By defining specific methods for calculating the unit direction vector and elastic deformation in the elastic characteristic model, and using real-time location information of adjacent cable segments to dynamically calculate the direction and magnitude of the tension, the accurate calculation of the tension vector as the line shape changes in real time is achieved, ensuring the vector accuracy of the elastic force calculation and the spatial consistency of the model.

[0060] In some specific embodiments, the first The aerodynamic characteristic model expression for each cable segment is as follows:

[0061] in, For the first The aerodynamic forces acting on each cable segment; ρ Represents air density; Represents the average wind vector in the ground coordinate system; Represents the normal force coefficient; To be Transformation matrix to the cable coordinate system; S This refers to the cross-sectional area of ​​the cable. For the first The speed of movement of each cable segment; For the first The relative speed of motion between the cable segment and the air .

[0062] By providing the first The model of the aerodynamic characteristics of a cable segment is given an explicit mathematical expression. The relative velocity is calculated based on the current wind vector and the moving speed of the cable segment. The aerodynamic force is solved by considering air density, normal force coefficient, coordinate transformation and cross-sectional area. This achieves a high-fidelity simulation of the complex aerodynamic interaction between wind and moving cable.

[0063] In some specific embodiments, the method for constructing the equivalent model of the iron tower is as follows: For each iron tower, the tower body is divided into several rectangular cross-section beams according to the mechanical force transmission path of the iron tower; A parameter calculation module is set up for each I-beam cross-section to calculate the moment of inertia of each I-beam cross-section. Cross-section beam quality and cross-section beam stiffness The expressions include:

[0064] Where A and B are the outer and inner side lengths of the rectangular section of the beam; The length of the cross-section beam; Density of the steel in the cross-section beam; The elastic modulus of the beam material; Based on the parameter calculation module of each cross-section beam, a preliminary equivalent model of each iron tower is built using Modelica's beam component; Modal analysis was performed on the preliminary model of each tower to calculate its first-order bending frequency. Then compare The first-order bending frequency of an actual iron towerf 1: like This reduces the stiffness of the beam section. and moment of inertia of cross section ; like This increases the stiffness of the beam section. and moment of inertia of cross section ; Then repeat the above calculations until the first-order bending frequency of the equivalent model of the tower is reached. Compared to the actual iron tower f1 If the error is ≤5%, the equivalent model of the corresponding iron tower is obtained.

[0065] By disassembling the tower into a U-shaped cross-section beam according to the mechanical force transmission path, calculating its cross-sectional properties, and building a preliminary equivalent model based on the Modelica library, and then iteratively adjusting the model parameters to make the error between its first-order bending frequency and the actual tower less than or equal to 5%, high-precision equivalence of the tower's dynamic characteristics is achieved, which significantly simplifies the tower model and ensures the accuracy of the simulation coupled with the line.

[0066] Step S3: Use Modelica language to encapsulate each parameterized sub-model into a model component, integrate the corresponding model component based on the actual topology of the high-voltage line system, define the initial state of the system, and establish dynamic equations for each cable segment to form a high-voltage line system model.

[0067] By using the Modelica language to encapsulate each parametric sub-model into corresponding reusable model components, and integrating these components based on the actual topology of the high-voltage line system, defining the initial state of the system, and establishing dynamic equations for each cable segment, a complete and flexible high-voltage line system model is formed, supporting unified modeling and simulation of multibody dynamics and multiphysics.

[0068] In some embodiments, the rigid body model components from the Modelica model library are used to encapsulate the rigid body models of each cable segment into rigid body model components. These rigid body model components (which have both graphical and 3D interfaces) can simulate cable models with a circular cross-section and include a dynamics interface and a kinematics interface. The dynamics interface is represented by a white connector, and kinematic variables are obtained by setting dynamic variables. The kinematics interface is represented by a yellow connector, and dynamic variables are obtained by setting motion variables. Dynamic variables include rigid body position, velocity, and acceleration, while kinematic variables include force and torque. Furthermore, by configuring the model parameters of the rigid body model in the Modelica model library, coordinate systems, mass properties, and cross-sectional dimension properties can be set.

[0069] In this application, the rigid body model of a single cable segment is causal, and the connection direction is unidirectional. To achieve identical model structures in both directions, FreeMotion is used. Adjacent ports are used for force calculation, and the dynamic ports connected to FreeMotion are used to calculate the six-degree-of-freedom results. The FreeMotion module simulates free motion with six spatial degrees of freedom, meaning the motion between the two interfaces of the MBS connector is completely unconstrained. Depending on the user's state selection, it generates either four state variables for the relative Euler symmetry parameter (quaternion) qRel or three state variables for the corresponding rotation angle phiRel. By default, the relative Euler symmetry parameter is used. Additionally, three state variables are generated for the relative angular velocity omRel, relative displacement xRel, and relative velocity vRel. A total of up to 13 potential state variables are generated.

[0070] In one embodiment of this application, the first cable segment is connected to both the left and right ends. The stress analysis diagram of each cable segment is shown below. Figure 2 As shown, where, Indicates the relationship with the first The leftmost cable segment adjacent to the first cable segment is paired with the second cable segment. The tensile force on each cable segment Indicates the relationship with the first The rightmost cable segment adjacent to the first cable segment is paired with the second cable segment. The tensile force of each cable segment.

[0071] In some embodiments of this application, the globalForce component from the Modelica model library is used to encapsulate the gravity model of each cable segment into a gravity model component.

[0072] In some embodiments of this application, the elasticity model of each cable segment is encapsulated into an elasticity model component using the elasticModel component in the Modelica model library. The model can be easily constructed based on non-causal modeling, and engineers only need to focus on the modeling characteristics theory related to elasticity.

[0073] In some embodiments of this application, the FreeMotion component in the Modelica model library is used to encapsulate the aerodynamic characteristic models of each cable segment into an aerodynamic characteristic model component. Based on non-causal modeling, the model can be easily constructed, and engineers only need to focus on the aerodynamic-related modeling characteristic theory.

[0074] In some specific embodiments, in step S3, the first The dynamic equations for each cable segment are expressed as follows:

[0075] in, Representing the The quality of each cable segment; Representing the The acceleration of each cable segment; Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; Representative and the All adjacent cable segments of the first cable segment are related to the second cable segment. The sum of the tensile forces on each cable segment; Representative and the The tower connected to the cable segment is the first The tensile force on each cable segment; Representing the The weight of each cable segment; Representing the The aerodynamic forces acting on each cable segment.

[0076] By explicitly giving the first The dynamic equations for each cable segment are expressed such that the product of its mass and acceleration equals the sum of the tension of the adjacent cable segment, the tension of the tower, gravity, and aerodynamic force. This establishes a complete mathematical description of Newton's second law for each discrete cable segment, forming the core mathematical foundation for the simulation and solution of the entire system.

[0077] In some specific embodiments, step S3 further includes: encapsulating the high-voltage line system model into a dynamic link library based on the FMI standard to support co-simulation with third-party simulation tools, specifically including: The high-voltage line system model is compiled into C code, and dynamic link libraries, output state variables, input interfaces and exported functions are generated according to the FMI2.0 standard. It enables co-simulation with third-party simulation tools through the FMI interface, supporting collaborative analysis of multiphysics fields. Third-party simulation tools include Simulink and Dymola.

[0078] By encapsulating the high-voltage line system model into a dynamic link library based on the FMI standard, providing C code, state variables, input interfaces, and exported functions, the model achieves standard interfaces and co-simulation capabilities with third-party simulation tools (such as Simulink and Dymola), thereby improving the model's interoperability and the application scope of multiphysics collaborative analysis.

[0079] Step S4: Obtain wind vector data of the current environment in real time.

[0080] By acquiring real-time wind vector data of the current environment, accurate and time-varying external excitation inputs are provided for the high-voltage line system model, enabling the simulation analysis to closely match the actual wind field conditions, thereby effectively evaluating the dynamic response and safety status of the line under specific wind conditions.

[0081] Step S5: Input the current environmental wind vector data into the high-voltage line system model for simulation analysis. Solve the dynamic equations of each cable segment through numerical integration, and calculate the velocity of each cable segment and the position of the midpoint of the cable segment at a specified future time. At the initial moment, the velocity of each cable segment is a preset value. In each subsequent integral time step, the aerodynamic force on the cable segment is calculated based on the velocity of the cable segment calculated in the previous time step, and then the dynamic equation of each cable segment is solved.

[0082] By inputting real-time wind vector data into the high-voltage line system model for simulation analysis, and using the numerical integration method to solve the dynamic equations of each cable segment to calculate its future velocity and position, efficient numerical simulation and prediction of the dynamic evolution process of the line system from the initial state under the coupling of aerodynamic loads and internal forces are achieved.

[0083] In some specific embodiments, a CVODE variable-order, variable-step integrator is used to solve the dynamic equations for each cable segment. Combined with high-stiffness damping and zero-crossing event shielding technology, computational efficiency is improved. Specifically, this includes: CVODE variable order and variable step size integration: The CVODE solver is used to dynamically adjust the integration order (1-5) and step size (1ms-1s) according to the local truncation error, balancing computational accuracy and efficiency. High stiffness damping treatment: Introducing a large damping coefficient into the elastic property model ( b =104~105Ns / m), suppressing the influence of high-frequency longitudinal vibration on the low-frequency response of wind vibration, and avoiding the solver having too small a step size due to "illness" problem; Zero-crossing event shielding: To address the discontinuity issue when the cable stretch crosses zero (ΔL=0), Modelica's noEvent mechanism is used to skip the precise moment search, retaining only the critical time step of low-frequency wind vibration (such as 1ms) to improve calculation speed.

[0084] By employing the CVODE variable-order, variable-step integrator to solve the dynamic equations, and combining high-stiffness damping treatment with Zero-crossing event shielding technology, the integration order and step size are dynamically adjusted to suppress high-frequency vibration and discontinuity issues. This significantly improves the computational efficiency and stability of numerical solutions while ensuring the accuracy of low-frequency response calculations for wind-induced vibrations.

[0085] In one specific embodiment, the steps of the high-voltage line wind vibration simulation analysis method based on the Modelica language include: Step S1: Obtain the physical parameters of the high-voltage line system, including cable parameters and tower parameters; The cable parameters include the cable's line span, span distribution, conductor cross-sectional area, elastic modulus, coefficient of linear expansion, and mass per unit length. The tower parameters include the tower's structural dimensions, tower height, crossarm position, tower body segment length, and material density.

[0086] Step S2, constructing a parameterized sub-model of the high-voltage line system based on physical parameters, including: Discretize the continuous cable into multiple cable segments, with each cable segment having a length ≤ 50 times the cable diameter; Construct rigid body models, gravity models, elastic property models, and aerodynamic property models for a single cable segment, as well as construct an equivalent model of the tower; The gravity model calculates the gravity acting on a cable segment based on the mass of a single cable segment. The gravity model expression for each cable segment is:

[0087] in Representing the The weight of each cable segment; Representing the The quality of each cable segment ; Mass per unit length of high-voltage line cable; Representing the The length of each cable segment; Represents gravitational acceleration; The elastic characteristic model calculates the elastic deformation by considering the positions of both ends of the cable segment and the preset length of the cable segment, and then calculates the tensile force of adjacent cable segments based on the elastic deformation. The elastic characteristic model expression for each cable segment is:

[0088] in, Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; k To increase the tensile stiffness of the cable; Representing the The elastic deformation of a cable segment; b The cable damping coefficient; represent The rate of change; represent The unit direction; Each cable segment's elasticity model is equipped with a left-side interface and a right-side interface; The calculation method is as follows: Calculate the current cable segment Left adjacent cable segment to current cable segment When applying tension, the left interface is used to input the real-time position of the adjacent cable segment to the left of the current cable segment. The interface on the right is used to input the real-time position of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the left of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to its left. ; Calculate the current cable segment Right adjacent cable segment to current cable segment When tension is applied, the left interface is used to input the real-time position of the current cable segment. The right-hand interface is used to input the real-time position of the cable segment adjacent to the right of the current cable segment. , The value is the Euclidean distance between the real-time position of the adjacent cable segment to the right of the current cable segment and the real-time position of the current cable segment. Subtract the initial distance between the current cable segment and the adjacent cable segment to the right of the current cable segment. ; The aerodynamic characteristic model calculates the aerodynamic forces acting on the cable segment based on the current wind vector and the cable segment's velocity. The aerodynamic characteristic model expression for each cable segment is as follows:

[0089] in, For the first The aerodynamic forces acting on each cable segment; ρ Represents air density; Represents the average wind vector in the ground coordinate system; Represents the normal force coefficient; To be Transformation matrix to the cable coordinate system; S This refers to the cross-sectional area of ​​the cable. For the first The speed of movement of each cable segment; For the first The relative speed of motion between the cable segment and the air ; The equivalent model of the tower is based on the displacement of the tower nodes. The tensile force exerted by the tower on adjacent cable segments is calculated through the equivalent stiffness matrix. The method for constructing the equivalent model of the tower is as follows: For each iron tower, the tower body is divided into several rectangular cross-section beams according to the mechanical force transmission path of the iron tower; A parameter calculation module is set up for each I-beam cross-section to calculate the moment of inertia of each I-beam cross-section. Cross-section beam quality and cross-section beam stiffness The expressions include:

[0090] Where A and B are the outer and inner side lengths of the rectangular section of the beam; The length of the cross-section beam; Density of the steel in the cross-section beam; The elastic modulus of the beam material; Based on the parameter calculation module for each cross-section beam, a preliminary equivalent model of each iron tower is built using the beam component of the Modelica model library; Modal analysis was performed on the preliminary model of each tower to calculate its first-order bending frequency. Then compare The first-order bending frequency of an actual iron tower f 1: like This reduces the stiffness of the beam section. and moment of inertia of cross section ; like This increases the stiffness of the beam section. and moment of inertia of cross section ; Then repeat the above calculations until the first-order bending frequency of the equivalent model of the tower is reached. Compared to the actual iron tower f1 If the error is ≤5%, the equivalent model of the corresponding iron tower is obtained.

[0091] Step S3: Using Modelica, each parametric sub-model is encapsulated into a model component. Based on the actual topology of the high-voltage line system, the corresponding model components are integrated, the initial state of the system is defined, and dynamic equations are established for each cable segment to form a high-voltage line system model. The dynamic equations for each cable segment are expressed as follows:

[0092] in, Representing the The quality of each cable segment; Representing the The acceleration of each cable segment; Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; Representative and the All adjacent cable segments of the first cable segment are related to the second cable segment. The sum of the tensile forces on each cable segment; Representative and the The tower connected to the cable segment is the first The tensile force on each cable segment; Representing the The weight of each cable segment; Representing the The aerodynamic forces acting on each cable segment; Based on the FMI standard, the high-voltage line system model is encapsulated into a dynamic link library, supporting co-simulation with third-party simulation tools, specifically including: The high-voltage line system model is compiled into C code, and dynamic link libraries, output state variables, input interfaces and exported functions are generated according to the FMI2.0 standard. It enables co-simulation with third-party simulation tools through the FMI interface, supporting collaborative analysis of multiphysics fields. Third-party simulation tools include Simulink and Dymola.

[0093] Step S4: Obtain wind vector data of the current environment in real time.

[0094] Step S5: Input the current wind vector data into the high-voltage line system model for simulation analysis. Use the CVODE variable-order variable-step integrator to solve the dynamic equations of each cable segment, calculate the velocity of each cable segment at a specified future time, and the position of the midpoint of the cable segment. Combine high-stiffness damping and zero-crossing event shielding technology to improve computational efficiency. At the initial moment, the velocity of each cable segment is a preset value. In each subsequent integral time step, the aerodynamic force on the cable segment is calculated based on the velocity of the cable segment calculated in the previous time step, and then the dynamic equation of each cable segment is solved. The CVODE solver is used to dynamically adjust the integration order (1-5th order) and step size (1ms-1s) based on the local truncation error, balancing computational accuracy and efficiency. A large damping coefficient is introduced into the elastic property model, and a setting is made. b =105Ns / m, suppressing the influence of high-frequency longitudinal vibration on the low-frequency response of wind vibration, and avoiding the solver having too small a step size due to "illness" problem; To address the discontinuity issue when the cable stretch crosses zero (ΔL=0), Modelica's noEvent mechanism is used to skip the precise moment search, retaining only the critical time step of low-frequency wind vibration, set to 1ms, to improve calculation speed.

[0095] The accuracy of the high-voltage line wind vibration simulation analysis method in this embodiment is verified by measuring wind vibration data. This method is then applied to the line response analysis under different wind field conditions and vibration suppression measures, specifically including: Model Validation: Measured wind vibration data under a periodic galloping wind condition of 10 m / s (frequency 0.04 Hz) were input into the model of this embodiment. The simulation outputs the three-directional displacement curves of the line (X-direction tower displacement 5 mm, Y-direction tilt due to self-weight, Z-direction synchronous fluctuation). The schematic diagram of the X-direction displacement simulation curve is shown below. Figure 3 The simulation curve of displacement in the Y direction is shown in the figure. Figure 4 The simulated displacement curve in the Z direction is shown in the figure. Figure 5 As shown.

[0096] from Figures 3 to 5 The simulation results show that the X direction is the wind direction, and the tower displacement is 5mm under the action of the cable swinging and pulling the tower. The Y direction is the cable pulling direction. The tower in the model only has cables on one side, so the tower is tilted in this direction due to the gravity of the cables. The tower displacement in the Z direction is particularly small, so the difference between the initial Z-direction and X-direction displacement changes is quite obvious.

[0097] A schematic flowchart of the simulation curve of synchronous displacement change in the XZ direction of the high-voltage line after gravity equilibrium is shown below. Figure 6 As shown, it can be seen that after gravity reaches equilibrium, the Z-axis displacement changes are basically synchronized with the X-axis displacement changes.

[0098] The simulation data obtained after model verification showed an error of less than 5% compared with the measured wind vibration data, indicating that the high-voltage line wind vibration simulation analysis method of this application supports rapid simulation of different wind fields and line parameters.

[0099] The following are embodiments of the high-voltage line wind vibration simulation analysis system based on the Modelica language provided in this application. This high-voltage line wind vibration simulation analysis system based on the Modelica language belongs to the same inventive concept as the high-voltage line wind vibration simulation analysis method in the above embodiments. For details not described in detail in the embodiments of the high-voltage line wind vibration simulation analysis system, please refer to the embodiments of the high-voltage line wind vibration simulation analysis method based on the Modelica language described above.

[0100] like Figure 7 As shown, the high-voltage line wind vibration simulation analysis system based on the Modelica language includes: The physical parameter acquisition module is used to acquire the physical parameters of the high-voltage line system, including cable parameters and tower parameters; The parameterized sub-model construction module is used to construct parameterized sub-models of high-voltage line systems based on physical parameters. This includes discretizing continuous cables into multiple cable segments, constructing rigid body models, gravity models, elastic property models, and aerodynamic property models for individual cable segments, and constructing equivalent models of towers. The high-voltage line system model building module is used to encapsulate various parameterized sub-models into model components using the Modelica language. Based on the actual topology of the high-voltage line system, it integrates the corresponding model components, defines the initial state of the system, and establishes dynamic equations for each cable segment to form a high-voltage line system model. The wind vector data acquisition module is used to acquire wind vector data of the current environment in real time. The simulation analysis module is used to input the wind vector data of the current environment into the high-voltage line system model for simulation analysis. It solves the dynamic equation of each cable segment through numerical integration and calculates the velocity of each cable segment and the position of the midpoint of the cable segment at a specified time in the future.

[0101] The high-voltage line wind vibration simulation analysis system in this embodiment is used to implement a high-voltage line wind vibration simulation analysis method based on the Modelica language.

[0102] This application also provides an electronic device for implementing the various embodiments of this application. Figure 8To illustrate the hardware structure of an electronic device according to various embodiments of this application, as shown in the following diagram... Figure 8 As shown, the electronic device includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor.

[0103] Those skilled in the art will understand that the electronic device structure involved in the embodiments of this application does not constitute a limitation on the electronic device. The electronic device may include more or fewer components than shown in the figure, or combine certain components, or have different component arrangements.

[0104] In embodiments of this application, electronic devices include, but are not limited to, laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. Electronic devices may also represent various forms of mobile devices and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the embodiments of this application described and / or claimed herein.

[0105] In this application embodiment, the processor can be implemented using at least one of an Application-Specific Integrated Circuit (ASIC), a Digital Signal Processor (DSP), a Digital Signal Processing Device (DSPD), a processor, a controller, a microcontroller, a microprocessor, or an electronic unit designed to perform the functions described herein. In some cases, such implementations can be implemented within a controller. For software implementations, implementations such as processes or functions can be implemented with separate software modules that allow the performance of at least one function or operation. The software code can be implemented by a software application (or program) written in any suitable programming language, and the software code can be stored in memory and executed by the controller.

[0106] In addition, the electronic device includes some functional modules not shown, which will not be described in detail here.

[0107] Those skilled in the art will understand that the various aspects of the electronic device provided in this application can be implemented as a system, method, or program product. Therefore, the various aspects of this application can be specifically implemented in the following forms: a completely hardware implementation, a completely software implementation (including firmware, microcode, etc.), or a combination of hardware and software aspects, collectively referred to herein as a "circuit," "module," or "system."

[0108] This application also provides a storage medium storing a program product capable of implementing a high-voltage line wind vibration simulation analysis method based on the Modelica language. In some possible implementations, various aspects of this application can also be implemented as a program product comprising program code that, when run on a terminal device, causes the terminal device to perform the steps described in the "Exemplary Methods" section of this specification according to various exemplary embodiments of this application.

[0109] The storage medium may be any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example,, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples (a non-exhaustive list) of readable storage media include: electrical connections having one or more wires, portable disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0110] The above description of the disclosed embodiments enables those skilled in the art to make or use this application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of this application. Therefore, this application is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A simulation analysis method for wind-induced vibration of high-voltage lines based on Modelica language, characterized in that, include: S1. Obtain the physical parameters of the high-voltage line system, including cable parameters and tower parameters; S2. Constructing a parameterized sub-model of the high-voltage line system based on physical parameters, including: Discretize the continuous cable into multiple cable segments; Construct rigid body models, gravity models, elastic property models, and aerodynamic property models for a single cable segment, as well as construct an equivalent model of the tower; The gravity model calculates the gravity acting on a cable segment by using the mass of a single cable segment. The elastic characteristic model calculates the elastic deformation by using the positions of both ends of the cable segment and the preset length of the cable segment, and calculates the tension of adjacent cable segments based on the elastic deformation. The aerodynamic characteristic model calculates the aerodynamic forces acting on the cable segment based on the current wind vector and the cable segment's velocity. The equivalent model of the tower is based on the displacement of the tower nodes, and the tension of the tower on adjacent cable segments is calculated by the equivalent stiffness matrix. S3. Using the Modelica language, each parameterized sub-model is encapsulated into a model component. Based on the actual topology of the high-voltage line system, the corresponding model components are integrated, the initial state of the system is defined, and dynamic equations are established for each cable segment to form a high-voltage line system model. S4. Acquire real-time wind vector data of the current environment; S5. Input the current wind vector data into the high-voltage line system model for simulation analysis. Solve the dynamic equations of each cable segment through numerical integration, and calculate the velocity of each cable segment and the position of the midpoint of the cable segment at a specified future time. At the initial moment, the velocity of each cable segment is a preset value. In each subsequent integral time step, the aerodynamic force on the cable segment is calculated based on the velocity of the cable segment calculated in the previous time step, and then the dynamic equation of each cable segment is solved.

2. The high-voltage line wind vibration simulation analysis method as described in claim 1, characterized in that, In step S2, the first The gravity model expression for each cable segment is: in Representing the The weight of each cable segment; Representing the The quality of each cable segment ; Mass per unit length of high-voltage line cable; Representing the The length of each cable segment; It represents gravitational acceleration.

3. The high-voltage line wind vibration simulation analysis method as described in claim 1, characterized in that, In step S2, the first The elastic characteristic model expression for each cable segment is: in, Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; k To increase the tensile stiffness of the cable; Representing the The elastic deformation of a cable segment; b The cable damping coefficient; represent The rate of change; represent The unit direction; Each cable segment's elasticity model is equipped with a left-side interface and a right-side interface.

4. The high-voltage line wind vibration simulation analysis method as described in claim 3, characterized in that, In step S2, the first The aerodynamic characteristic model expression for each cable segment is: in, For the first The aerodynamic forces acting on each cable segment; ρ Represents air density; Represents the average wind vector in the ground coordinate system; Represents the normal force coefficient; To be Transformation matrix to the cable coordinate system; S This refers to the cross-sectional area of ​​the cable. For the first The speed of movement of each cable segment; For the first The relative speed of motion between the cable segment and the air .

5. The high-voltage line wind vibration simulation analysis method as described in claim 1, characterized in that, In step S2, the method for constructing the equivalent model of the iron tower is as follows: For each iron tower, the tower body is divided into several rectangular cross-section beams according to the mechanical force transmission path of the iron tower; A parameter calculation module is set up for each I-beam cross-section to calculate the moment of inertia of each I-beam cross-section. Cross-section beam quality and cross-section beam stiffness The expressions include: Where A and B are the outer and inner side lengths of the rectangular section of the beam; The length of the cross-section beam; Density of the steel in the cross-section beam; The elastic modulus of the beam material; Based on the parameter calculation module for each cross-section beam, a preliminary equivalent model of each iron tower is built using the beam component of the Modelica model library; Modal analysis was performed on the preliminary model of each tower to calculate its first-order bending frequency. Then compare The first-order bending frequency of an actual iron tower f 1: like This reduces the stiffness of the beam section. and moment of inertia of cross section ; like This increases the stiffness of the beam section. and moment of inertia of cross section ; Then repeat the above calculations until the first-order bending frequency of the equivalent model of the tower is reached. Compared to the actual iron tower f1 If the error is ≤5%, the equivalent model of the corresponding iron tower is obtained.

6. The high-voltage line wind vibration simulation analysis method as described in claim 1, characterized in that, In step S3, the first The dynamic equations for each cable segment are expressed as follows: in, Representing the The quality of each cable segment; Representing the The acceleration of each cable segment; Representing the The adjacent cable segments of the first cable segment are paired with the second cable segment. The tensile force on each cable segment; Representative and the All adjacent cable segments of the first cable segment are related to the second cable segment. The sum of the tensile forces on each cable segment; Representative and the The tower connected to the cable segment is the first The tensile force on each cable segment; Representing the The weight of each cable segment; Representing the The aerodynamic forces acting on each cable segment.

7. The high-voltage line wind vibration simulation analysis method as described in claim 1, characterized in that, Step S3 also includes: encapsulating the high-voltage line system model into a dynamic link library based on the FMI standard to support co-simulation with third-party simulation tools.

8. A high-voltage line wind vibration simulation analysis system based on Modelica language, characterized in that, The method for implementing the high-voltage line wind vibration simulation analysis method as described in any one of claims 1-7 includes: The physical parameter acquisition module is used to acquire the physical parameters of the high-voltage line system, including cable parameters and tower parameters; The parameterized sub-model construction module is used to construct parameterized sub-models of high-voltage line systems based on physical parameters. This includes discretizing continuous cables into multiple cable segments, constructing rigid body models, gravity models, elastic property models, and aerodynamic property models for individual cable segments, and constructing equivalent models of towers. The high-voltage line system model building module is used to encapsulate various parameterized sub-models into model components using the Modelica language. Based on the actual topology of the high-voltage line system, it integrates the corresponding model components, defines the initial state of the system, and establishes dynamic equations for each cable segment to form a high-voltage line system model. The wind vector data acquisition module is used to acquire wind vector data of the current environment in real time. The simulation analysis module is used to input the wind vector data of the current environment into the high-voltage line system model for simulation analysis. It solves the dynamic equation of each cable segment through numerical integration and calculates the velocity of each cable segment and the position of the midpoint of the cable segment at a specified time in the future.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes a computer program, it implements the steps of the high-voltage line wind vibration simulation analysis method as described in any one of claims 1-7.

10. A storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the high-voltage line wind vibration simulation analysis method as described in any one of claims 1-7.

Images