Dynamic evaluation method and system for aeolian vibration of power transmission line

Abstract

The invention relates to the technical field of power transmission line evaluation, and discloses a power transmission line aeolian vibration dynamic evaluation method and system.The method comprises the steps that when environmental meteorological data meet meteorological sudden change conditions, relative errors between simulation data and actually-measured vibration data and between simulation data and actually-measured tension data are calculated respectively, and if the errors exceed a preset threshold value, the power transmission line is evaluated; if yes, optimizing the initial coupling simulation model, then calculating fatigue damage distribution and residual life of the wire along the span, and generating a fatigue life evaluation result; and finally, under a plurality of constraint conditions, carrying out iterative solution on a preset damper installation optimization model, and determining a damper installation optimization scheme. By constructing a complete technical link of'multi-source data input-coupling simulation calculation-dynamic calibration optimization-fatigue evaluation-anti-vibration optimization ', the damper installation optimization method has the advantages that the anti-vibration performance of the damper is improved, and the anti-vibration performance of the damper is improved. The problems that the dynamic coupling effect of a wind field-wire-damper system is ignored, and the dynamic influence of wind speed and wind direction changes of the wind field on the wire vibration response cannot be captured in real time in the prior art are effectively solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of power transmission line evaluation, in particular to a power transmission line aeolian vibration dynamic evaluation method and system. BACKGROUND

[0002] In the natural environment, the power transmission line is prone to periodic aeolian vibration caused by the micro wind. Although the vibration amplitude is small, the vibration frequency is high and the duration is long, which will make the conductor repeatedly bear alternating stress. The connection parts of the conductor, the damper and the insulator, and the aluminum strand layer of the conductor itself will be damaged due to long-term alternating stress accumulation. With the passage of time, the damage will gradually intensify, and eventually may lead to conductor breakage, wire breakage accident, and even line outage, causing huge economic losses.

[0003] However, the traditional method is a point monitoring method using a single acceleration sensor. The vibration data is collected by installing an acceleration sensor at a local point of the conductor, and then the fatigue damage is calculated. However, the above method is limited by single-point measurement characteristics, cannot cover the entire conductor, is difficult to reflect the spatial distribution law of the conductor vibration mode, and is prone to blind areas, leading to misjudgment of the fatigue danger section.

[0004] Therefore, the existing technology usually uses an evaluation method of a static empirical model to estimate the conductor stress and fatigue life by presetting the empirical formula of wind conditions and vibration response. However, the above method uses a static empirical formula, ignores the dynamic coupling effect of the wind field-conductor-damper system, cannot capture the dynamic influence of the change of wind speed and direction of the wind field on the vibration response of the conductor in real time, and leads to large deviation of the evaluation result from the actual working condition. SUMMARY

[0005] The present application provides a power transmission line aeolian vibration dynamic evaluation method and system, which solves the technical problem that the existing technology ignores the dynamic coupling effect of the wind field-conductor-damper system, cannot capture the dynamic influence of the change of wind speed and direction of the wind field on the vibration response of the conductor in real time, and leads to large deviation of the evaluation result from the actual working condition.

[0006] The first aspect of the present application provides a power transmission line aeolian vibration dynamic evaluation method, comprising:

[0007] inputting the vibration data, tension data and environmental meteorological data of the conductor of the power transmission line into a preset initial coupling simulation model for simulation, and outputting simulation data;

[0008] under the meteorological mutation condition of the environmental meteorological data, calibrating the initial coupling simulation model according to the error value between the simulation data and the vibration data and the tension data, respectively;

[0009] Based on the target simulation data corresponding to the calibrated coupling simulation model, fatigue life of the conductor of the power transmission line is evaluated, and a fatigue life evaluation result is generated.

[0010] Under multiple constraint conditions, an iteration solution is performed on the preset damper installation optimization model, a damper installation optimization scheme is determined according to a solution result, and the damper installation optimization scheme and the fatigue life evaluation result are output to an evaluation platform.

[0011] Optionally, the inputting of the vibration data, the tension data and the environmental meteorological data of the conductor of the power transmission line into the preset initial coupling simulation model and the outputting of simulation data include:

[0012] The vibration data, the tension data and the environmental meteorological data of the conductor of the power transmission line are acquired.

[0013] The vibration data, the tension data and the environmental meteorological data are preprocessed.

[0014] Based on the preprocessed vibration data, tension data and environmental meteorological data, an initial coupling simulation model is constructed; wherein the initial coupling simulation model includes a fluid dynamics model and a finite element model.

[0015] The preprocessed vibration data, tension data and environmental meteorological data are input into the initial coupling simulation model for simulation, and point wind load, displacement data, stress data, simulation amplitude and simulation frequency of the conductor are output.

[0016] The simulation data are generated in combination with the point wind load, the displacement data, the stress data, the simulation amplitude and the simulation frequency of the conductor.

[0017] Optionally, the inputting of the preprocessed vibration data, tension data and environmental meteorological data into the initial coupling simulation model for simulation and the outputting of the point wind load, the displacement data, the stress data, the simulation amplitude and the simulation frequency of the conductor include:

[0018] The preprocessed environmental meteorological data are input into the fluid dynamics model of the initial coupling simulation model for simulation, and the point wind load of the conductor is output.

[0019] The point wind load, the preprocessed vibration data and tension data of the conductor are input into the finite element coupling model of the initial coupling simulation model for simulation, and the displacement data, the stress data, the simulation amplitude and the simulation frequency are output.

[0020] Optionally, the initial coupled simulation model is calibrated according to error values between the simulation data and the vibration data and the tension data respectively under the weather mutation condition of the environmental meteorological data, including:

[0021] calculating a first error value between a simulation amplitude of the simulation data and a measured amplitude of the vibration data;

[0022] determining a dynamic component of a support reaction force of a finite element model of the initial coupled simulation model based on displacement data and mechanical properties of the simulation data;

[0023] calculating a second error value between the dynamic component of the support reaction force and a dynamic tension increment of the tension data;

[0024] when the first error value or the second error value is greater than a preset error threshold and the weather of the environmental meteorological data is in a weather mutation condition, calibrating the initial coupled simulation model.

[0025] Optionally, a fatigue life of the conductor of the power transmission line is evaluated based on target simulation data corresponding to the calibrated coupled simulation model, and a fatigue life evaluation result is generated, including:

[0026] vibration data, tension data and environmental meteorological data of the conductor of the power transmission line are input into the calibrated coupled simulation model for simulation to generate new simulation data;

[0027] jumping to execute the step of calculating the first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data until the first error value and the second error value are both less than or equal to the preset error threshold, and then determining the current simulation data as the target simulation data;

[0028] extracting a stress amplitude of stress data of the target simulation data and a cycle number corresponding to the stress amplitude;

[0029] based on a conductor stress-life curve, the stress amplitude and the cycle number corresponding to the stress amplitude are used to calculate a damage distribution of the conductor;

[0030] based on the damage distribution of the conductor, a maximum damage degree of the conductor is determined;

[0031] the fatigue life of the conductor is evaluated using the maximum damage degree of the conductor and a preset damage rate, and a fatigue life evaluation result is generated.

[0032] Optionally, the method comprises: under multiple constraint conditions, iteratively solving a preset anti-vibration hammer installation optimization model, determining an anti-vibration hammer installation optimization scheme according to a solving result, and outputting the anti-vibration hammer installation optimization scheme and the fatigue life evaluation result to an evaluation platform.

[0033] A first constraint condition is constructed based on a distance between a clamp of the conductor and an anti-vibration hammer installation position.

[0034] A second constraint condition is constructed according to a distance between two adjacent anti-vibration hammers.

[0035] A third constraint condition is constructed according to the number of the anti-vibration hammers.

[0036] An anti-vibration hammer installation optimization model is constructed under the first constraint condition, the second constraint condition and the third constraint condition, with minimization of unevenness of damage distribution of the conductor as a target.

[0037] A genetic algorithm is used to iteratively solve the anti-vibration hammer installation optimization model, to obtain an optimal number of the anti-vibration hammers and respective anti-vibration hammer installation positions.

[0038] An anti-vibration hammer installation optimization scheme is generated in combination with the optimal number of the anti-vibration hammers and the respective anti-vibration hammer installation positions.

[0039] The anti-vibration hammer installation optimization scheme and the fatigue life evaluation are output to the evaluation platform.

[0040] The second aspect of the application provides a power transmission line aeolian vibration dynamic evaluation system, comprising:

[0041] A simulation module is configured to input vibration data, tension data and environmental meteorological data of a conductor of a power transmission line into a preset initial coupling simulation model for simulation, and output simulation data.

[0042] A calibration module is configured to calibrate the initial coupling simulation model according to error values between the simulation data and the vibration data and the tension data under a meteorological mutation condition of the environmental meteorological data.

[0043] An evaluation module is configured to perform fatigue life evaluation on the conductor of the power transmission line based on target simulation data corresponding to the calibrated coupling simulation model, and generate a fatigue life evaluation result.

[0044] An output module is configured to, under multiple constraint conditions, iteratively solve a preset anti-vibration hammer installation optimization model, determine an anti-vibration hammer installation optimization scheme according to a solving result, and output the anti-vibration hammer installation optimization scheme and the fatigue life evaluation result to an evaluation platform.

[0045] The third aspect of the present application provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer program is executed by the processor to make the processor execute the steps of the power transmission line aeolian vibration dynamic evaluation method according to any one of the above aspects.

[0046] The fourth aspect of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed to implement the virtual impedance control parameter optimization method according to any one of the above aspects.

[0047] The fifth aspect of the present application provides a computer program product, characterized in that the computer program product comprises a computer program stored on a non-transitory computer readable storage medium, and the computer program comprises program instructions, wherein when the program instructions are executed by a computer, the computer executes the power transmission line aeolian vibration dynamic evaluation method according to any one of the above aspects.

[0048] From the above technical solutions, the present application has the following advantages:

[0049] The present application obtains conductor vibration data, tension data and environmental meteorological data, inputs the three types of data into a preset initial coupling simulation model for simulation calculation, and outputs simulation data. Then, when the environmental meteorological data meets the meteorological mutation condition, the relative errors of the simulation data and the measured vibration data and the measured tension data are calculated respectively. If the error exceeds the preset threshold, the key parameters of the initial coupling simulation model are optimized, the model is calibrated, and the target simulation data is output. Then, the fatigue damage distribution and the remaining life of the conductor along the span are calculated, and the fatigue life evaluation result is generated. Finally, under multiple constraint conditions, the preset damper installation optimization model is iteratively solved to determine the damper installation optimization scheme, and the damper installation optimization scheme and the fatigue life evaluation result are synchronously output to the evaluation platform. The present application constructs a complete technical link of "multi-source data input-coupling simulation calculation-dynamic calibration optimization-fatigue evaluation-damper optimization", effectively solves the problem that the prior art ignores the dynamic coupling effect of the wind field-conductor-damper system, and cannot capture the dynamic influence of the change of wind speed and direction of the wind field on the conductor vibration response in real time. BRIEF DESCRIPTION OF DRAWINGS

[0050] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments or the prior art description will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present application, and those skilled in the art can also obtain other drawings according to these drawings without creative labor.

[0051] Figure 1A step flow chart of a dynamic evaluation method of a power transmission line aeolian vibration is provided for the embodiment one of the present application.

[0052] Figure 2 A schematic diagram of a distributed fiber layout cross section is provided for the embodiment one of the present application.

[0053] Figure 3 A schematic diagram of a distributed fiber layout span view is provided for the embodiment one of the present application.

[0054] Figure 4 A structural block diagram of a dynamic evaluation system of a power transmission line aeolian vibration is provided for the embodiment two of the present application.

[0055] Figure 5 A structural block diagram of a computer device is provided for the embodiment three of the present application.

[0056] Wherein, the meaning of the reference signs is as follows:

[0057] 1, conductor; 2, optical fiber; 3, damper installation position; 11, steel core; 12, aluminum strand; 21, optical fiber sensing node. DETAILED DESCRIPTION

[0058] The embodiment of the present application provides a dynamic evaluation method and system of a power transmission line aeolian vibration, which is used for solving the technical problem that the dynamic coupling effect of a wind field-conductor-damper system is ignored in the prior art, the dynamic influence of the change of wind speed and direction of a wind field on the vibration response of a conductor cannot be captured in real time, and the evaluation result is greatly deviated from the actual working condition.

[0059] In order to make the invention purpose, features and advantages of the present application more obvious and easy to understand, the technical solutions in the embodiments of the present application will be described clearly and completely in combination with the drawings in the embodiments of the present application. Obviously, the following described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the protection scope of the present application.

[0060] Please refer to Figures 1 to 3 , Figure 1 A step flow chart of a dynamic evaluation method of a power transmission line aeolian vibration is provided for the embodiment one of the present application.

[0061] The dynamic evaluation method of a power transmission line aeolian vibration provided by the present application comprises the following steps.

[0062] In step 101, the vibration data, tension data and environmental meteorological data of the conductor 1 of the power transmission line are input into a preset initial coupling simulation model for simulation, and simulation data is output.

[0063] In the embodiments of the present application, the vibration data refers to the dynamic vibration response data generated by the conductor 1 of the distributed transmission line under the excitation of wind and the like, which is collected by a distributed optical fiber sensor and the like, and contains key parameters such as vibration amplitude, vibration frequency and dynamic strain time history of the conductor 1 at each point along the span.

[0064] The tension data refers to the tension data acting on the conductor 1 of the transmission line, which is collected by a tension sensor installed at the hanging point of the insulator string, and contains static tension and dynamic tension increment.

[0065] The environmental meteorological data refers to the meteorological parameters of the environment where the transmission line is located, which is collected by a multi-point weather station in the span, and mainly includes wind speed, wind direction, air density, turbulence intensity and the like.

[0066] The initial coupling simulation model refers to a coupling simulation system that is pre-constructed and combines a computational fluid dynamics (CFD, Computational Fluid Dynamics) wind field model and a finite element (FEM, Finite Element Method) structure vibration model, and is used to simulate the dynamic interaction of “wind field-conductor 1-damper”.

[0067] The simulation data refers to the result data output by the initial coupling simulation model, and contains vibration amplitude, frequency, stress time history and dynamic component of the support reaction force of the conductor 1.

[0068] The vibration data, tension data and environmental meteorological data of the conductor 1 of the transmission line are collected by the multi-source monitoring module, and the three types of data are input into the preset initial coupling simulation model; the initial coupling simulation model is composed of a CFD wind field simulation unit and a conductor 1-damper-insulator FEM coupled vibration model, first takes the wind speed distribution in the environmental meteorological data as the inlet boundary condition, solves the Navier-Stokes equation through the CFD wind field simulation unit, outputs the wind load data on the surface of the conductor 1 and transmits it to the FEM coupled vibration model, and then combines the static tension in the tension data as the initial boundary constraint of the FEM model, solves the system motion differential equation to obtain the simulation data such as displacement, stress and vibration response of the conductor 1.

[0069] Further, the step 101 includes the following sub-steps:

[0070] S11, obtaining the vibration data, tension data and environmental meteorological data of the conductor 1 of the transmission line.

[0071] In the embodiment of the present application, the vibration data, tension data and environmental meteorological data of the conductor 1 of the transmission line are collected by the multi-source monitoring module, wherein the distributed optical fiber sensor is arranged with optical fiber sensing nodes 21 at an interval of 5m along the span of the conductor 1, the dynamic strain at each point of the conductor 1 is collected in real time and converted into vibration data, the spoke type tension sensor is installed at the hanging point of the insulator string to collect the total tension of the conductor 1 and separate the tension data, and the three-point weather station in the span collects the environmental meteorological data such as wind speed, wind direction, temperature and humidity.

[0072] Specifically, referring to Figure 2 and Figure 3 , the conductor 1 is provided with a steel core 11 and an aluminum strand 12; the optical fiber 2 is embedded in the second layer of the aluminum strand 12; and the optical fiber sensing nodes 21 are arranged along the direction of the conductor 1 at an interval of 5m.

[0073] S12, the vibration data, tension data and environmental meteorological data are preprocessed.

[0074] In the embodiment of the present application, preprocessing refers to wavelet filtering, Fourier transform, low-pass filtering and other data processing methods for vibration data, tension data and environmental meteorological data.

[0075] The vibration data is filtered by wavelet to remove environmental electromagnetic interference and sensor noise, then converted into acceleration signal by twice integration, and finally the standardized vibration data such as vibration amplitude and dominant frequency is extracted by Fourier transform.

[0076] The tension data is filtered by low-pass filter to separate static tension and dynamic tension increment, and the jump value caused by instantaneous disturbance of the sensor is removed to form standardized tension data.

[0077] For the collected raw meteorological signals such as wind speed, wind direction, temperature and humidity, the three-point wind speed data is first interpolated and mapped to each position of the conductor 1 according to the exponential wind profile formula, the non-uniformity of the wind field caused by the difference in terrain is corrected, and the abnormal data with instantaneous fluctuation of wind speed exceeding 3m / s and sudden change of wind direction exceeding 60° is removed to obtain standardized environmental meteorological data.

[0078] S13, based on the preprocessed vibration data, tension data and environmental meteorological data, an initial coupled simulation model is constructed; wherein the initial coupled simulation model includes a fluid dynamics model and a finite element model.

[0079] In the embodiment of the present application, the fluid dynamics model refers to a wind field simulation model constructed based on the theory of computational fluid dynamics, which simulates the real wind field distribution and the wind load on the surface of the conductor 1 by solving the fluid motion equation (Navier-Stokes equation).

[0080] Finite element model refers to a structural vibration model constructed based on finite element method, which discretizes components such as conductor 1, damping hammer, insulator and the like into units, and calculates vibration response by solving system motion differential equation.

[0081] The normalized wind speed distribution in the pretreated meteorological data is taken as the inlet boundary condition of the fluid dynamics model, the fluid physical parameters are determined in combination with the measured air density, temperature and humidity, the k-omega SST turbulence model is used to define the turbulence characteristics of the wind field, the upstream 15D x downstream 20D x lateral 10D calculation domain is constructed with conductor 1 as the center, and the line terrain data is imported to correct the wind field flow boundary, and the fluid dynamics model is built; then, the finite element model is built relying on the tension data and the vibration data, the static tension in the pretreated tension data is taken as the initial tensioning condition of the conductor 1, the material parameters of the conductor 1 are set, the LINK180 unit is used to discretize the conductor 1, the COMBIN14 spring-damping unit is selected for the damping hammer, and the BEAM188 beam unit is used to simulate the bending characteristics of the insulator, and the vibration amplitude and frequency in the pretreated vibration data are taken as the initial verification benchmark of the finite element model, so that the model grid division and boundary constraint are ensured to meet the actual vibration scene; finally, the coupling correlation of the two models is realized, that is, the fluid dynamics model solves the Navier-Stokes equation to obtain the wind load data of the conductor 1 surface, which is mapped to the node force load of the finite element model through radial basis function interpolation, the finite element model solves the motion differential equation to output the displacement response of the conductor 1, and the displacement data is fed back to the fluid dynamics model to correct the wind field flow boundary to simulate the dynamic interaction of the wind field and the conductor 1, and finally an initial coupled simulation model including the fluid dynamics model and the finite element model is formed.

[0082] S14, input the pretreated vibration data, tension data and environmental meteorological data into the initial coupled simulation model for simulation, and output the wind load, displacement data, stress data, simulation amplitude and simulation frequency of each point of the conductor 1.

[0083] In the embodiment of the application, the wind load at each point refers to the wind load acting on each discrete point on the surface of the conductor 1 calculated by the fluid dynamics model, which includes the load size and the action direction.

[0084] The displacement data refers to the spatial position change data of the conductor 1 at each point along the span, which includes the time history of the transverse and longitudinal displacement.

[0085] The stress data refers to the internal force state data of the conductor 1 derived based on the displacement data, which includes the dynamic stress time history and the maximum stress value, which is obtained through geometric equation and physical equation.

[0086] The simulation amplitude refers to the displacement peak value extracted after the displacement data output by the finite element model is subjected to Fourier transform.

[0087] Simulation frequency refers to the dominant frequency of vibration identified after frequency domain analysis of displacement data.

[0088] The pre-processed environmental meteorological data is taken as the driving input of the fluid dynamics model, parameters such as the normalized wind speed distribution, air density, and turbulence intensity are substituted into the model, the inlet is set as a velocity boundary, the outlet is set as a pressure boundary, and the surface of the conductor 1 is set as a no-slip wall surface, the Navier-Stokes equation is solved to obtain the global wind pressure distribution on the surface of the conductor 1, and then the wind load of each point of the conductor 1 is calculated by integration; then the wind load is mapped to the node force load of the finite element model through radial basis function interpolation, and the pre-processed tension data is imported, the static tension is taken as the initial boundary constraint of the conductor 1, the dynamic tension increment is taken as the model verification benchmark, the simulation time step is set according to the frequency range in the pre-processed vibration data, the motion differential equation of the finite element model is solved, and the displacement data of the conductor 1 at each point along the span is obtained; based on the displacement data, the unit strain is calculated through the geometric equation, for example, strain = displacement gradient, and then the stress data of each point is output according to the physical equation, for example, stress = elastic matrix x strain, according to the material elastic matrix of the conductor 1; finally, the displacement data is subjected to Fourier transform, the peak displacement of vibration of each point is extracted as the simulation amplitude, and the dominant peak frequency in the frequency spectrum is identified as the simulation frequency.

[0089] Further, the step S14 includes the following sub-steps:

[0090] S141, input the pre-processed environmental meteorological data into the fluid dynamics model of the initial coupled simulation model for simulation, and output the wind load of each point of the conductor 1.

[0091] In the embodiment of the application, the model boundary and physical property configuration are completed according to the key parameters in the pre-processed meteorological data, the normalized wind speed distribution is taken as the inlet velocity boundary of the fluid dynamics model, the fluid physical parameters are set in combination with the pre-processed air density p and molecular dynamic viscosity μ, the line terrain data is imported to define the terrain boundary of the calculation domain, the k-ω SST turbulence model is used to simulate the turbulence characteristics of the wind field, and the model is ensured to be consistent with the actual atmospheric environment; then the fluid dynamics model solving process is started, the Navier-Stokes equation is numerically dispersed based on the finite volume method, the wind field flow field distribution is iteratively solved with a time step of 0.02 s, the pressure gradient change on the surface of the conductor 1 is calculated, the force size and direction of each element on the surface of the conductor 1 are obtained by pressure integration of the discrete elements, and finally the wind load of each point of the conductor 1 along the span is integrated and output.

[0092] S142, input the wind load of each point of the conductor 1, the pre-processed vibration data and the tension data into the finite element coupling model of the initial coupled simulation model for simulation, and output the displacement data, the stress data, the simulation amplitude and the simulation frequency.

[0093] In the embodiment of the present application, the initial state configuration of the model is completed based on the pretreated tension data, the static tension of the tension data is taken as the initial tensioning condition of the conductor 1, the initial strain of the conductor unit is iteratively adjusted through the statics analysis module of the finite element software, so that the deviation between the static tension calculated by the model and the measured static tension is ≤5%, and the dynamic tension increment of the tension data is taken as the verification benchmark of the dynamic response of the model; then the wind load of each point of the conductor 1 is imported, the discrete wind load data is mapped to the node force load of the conductor unit in the finite element model through radial basis function interpolation, and it is ensured that the spatial distribution of the wind load is consistent with the output of the fluid dynamics model; then the frequency range in the pretreated vibration data is combined, the transient analysis time step of the finite element model is set to 0.05s, the Newmark-β integration method is called to solve the system motion differential equation, and the displacement data of the conductor 1 at each point along the span is obtained; based on the displacement data, the unit strain is calculated through the geometric equation (strain = displacement gradient), and then the stress data of each point is output according to the physical equation (stress = elastic matrix * strain) combined with the elastic matrix of the conductor material; finally, the Fourier transform is performed on the displacement data, the peak value of the displacement time history is extracted as the simulation amplitude, and the frequency component with the maximum energy in the frequency domain spectrum is identified as the simulation frequency.

[0094] S15, combining the wind load, displacement data, stress data, simulation amplitude and simulation frequency of each point of the conductor 1, simulation data is generated.

[0095] In the embodiment of the present application, the wind load at each point is indexed according to the span position z and the time step t of the conductor 1, and the displacement data and stress data at the corresponding time step are synchronously associated to ensure the one-to-one correspondence between the wind load and the structure response data in space and time; then the simulation amplitude and the simulation frequency are supplemented with annotations, the simulation amplitude and the simulation frequency of each spatial position z are associated with the maximum stress value of the position to form a mapping relationship of “position-amplitude-frequency-stress”; at the same time, a data validity checking mechanism is introduced to eliminate abnormal data caused by unstable numerical iteration, and finally the integrated wind load at each point, displacement time history, stress time history, localized simulation amplitude and simulation frequency are packaged into a simulation data set according to a unified data format.

[0096] Step 102, under the meteorological mutation condition of the environmental meteorological data, the initial coupled simulation model is calibrated according to the error values between the simulation data and the vibration data and the tension data respectively.

[0097] In the embodiment of the present application, the meteorological mutation condition refers to the determination standard for representing the drastic change of the wind field in the environmental meteorological data, specifically the wind speed change rate ≥1.5m / s², the wind direction change rate ≥30° / s and the duration ≥10s.

[0098] The error value refers to a quantitative index for measuring the deviation of the simulation data from the measured data, including the relative error of the simulation amplitude and the measured amplitude in the vibration data, and the deviation of the simulation frequency and the measured frequency.

[0099] Under the condition that the environmental meteorological data satisfies the meteorological mutation condition of a wind speed change rate ≥ 1.5 m / s², a wind direction change rate ≥ 30° / s, and a duration ≥ 10s, the simulation data and the measured vibration data and tension data in the corresponding period of the meteorological mutation condition are obtained; then the error value is calculated, for the vibration data, the relative error of the simulation amplitude and the measured amplitude is taken as the core index, and the deviation of the simulation frequency and the measured frequency is compared; for the tension data, the absolute error of the simulation dynamic tension increment and the measured dynamic tension increment is mainly calculated; if any error value exceeds a preset threshold value, the calibration process of the initial coupled simulation model is started, that is, for the fluid dynamics model, the turbulence intensity correction coefficient of the k-ω SST turbulence model is adjusted based on the error feedback to optimize the calculation accuracy of the wind load; for the finite element coupled model, the damper coefficient of the anti-vibration hammer is corrected based on the tension error, and the conductor 1 unit stiffness matrix is fine-tuned combined with the vibration frequency deviation; after calibration, the coupled simulation is re-run, the error calculation and parameter adjustment are repeated, and until all error values fall within the threshold range, the model calibration is completed and the calibrated coupled simulation model is output.

[0100] Further, step 102 comprises the following sub-steps:

[0101] S21, calculating a first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data.

[0102] In the embodiment of the application, the first error value refers to a deviation quantitative index between the simulation amplitude in the simulation data and the measured amplitude in the vibration data.

[0103] The amplitude data of the corresponding space-time dimension is extracted from the simulation data and the vibration data, and the simulation amplitude and the measured amplitude at the same time period and the same conductor 1 position z are screened out to ensure that the two groups of data are completely aligned in space-time; then the relative error formula is used to calculate the first error value, and the specific formula is:

[0104] The first error value = |simulation amplitude-measured amplitude| / measured amplitude*100%.

[0105] S22, determining the dynamic component of the support reaction force of the finite element model of the initial coupled simulation model based on the displacement data and the mechanical properties of the simulation data.

[0106] In the embodiment of the application, the mechanical properties refer to the core parameters in the finite element model that affect the structural dynamic response, mainly including the mass matrix, the damping matrix and the stiffness matrix.

[0107] The dynamic component of the support reaction force refers to a reaction force component generated at a boundary constraint position of the finite element model and used for balancing the dynamic load of the conductor 1 vibration, and is a high-frequency fluctuation part of the support reaction force over time.

[0108] According to displacement data in simulation data, a structure dynamics basic equation is called through a finite element post-processing module, second-order derivation is performed on the displacement data to obtain an acceleration vector, first-order derivation is performed to obtain a velocity vector, and then the system motion differential equation is substituted into the mechanical characteristic parameters of the finite element model, and the support reaction force vector is derived through equation transformation; then, the support reaction force vector is separated into a static component and a dynamic component, a low-pass filter is used to remove the slowly changing static component, and the remaining high-frequency fluctuation part is the dynamic component of the support reaction force; finally, effectiveness verification is performed to ensure that the peak value is synchronized with the peak time of the displacement data, and the numerical range conforms to the common sense of mechanics, and finally the dynamic component of the support reaction force along the time sequence is output.

[0109] S23, a second error value between the dynamic component of the support reaction force and the dynamic tension increment of the tension data is calculated.

[0110] In the embodiment of the application, the dynamic tension increment refers to a periodic fluctuation component generated by conductor 1 vibration in the tension data, which is obtained after the static component is separated from the total tension.

[0111] The second error value refers to a deviation quantization index between the dynamic component of the support reaction force output by the finite element model and the dynamic tension increment actually measured by the tension sensor.

[0112] The time sequence of the dynamic component of the support reaction force is extracted from the finite element model output, and the dynamic tension increment of the corresponding period is called from the preprocessed tension data, so that the two groups of data are compared under the same time dimension; then, the second error value is calculated by using the absolute error formula, and the formula is:

[0113] Second error value = | ΔR(t) - ΔF(t) | / ΔF(t) * 100%

[0114] In the formula, ΔR(t) is the dynamic component of the support reaction force, and ΔF(t) is the dynamic tension increment.

[0115] Taking the absolute value can eliminate the interference of positive and negative deviations on error judgment, and directly reflect the difference between the two support values at a certain time; considering that instantaneous fluctuations may cause individual abnormal errors, sliding average processing is required for the calculation results to smooth the influence of short-term disturbances, and the average second error value along the time sequence is obtained.

[0116] S24, when the first error value or the second error value is greater than a preset error threshold value, and the weather of the environmental weather data is in a weather mutation condition, the initial coupled simulation model is calibrated.

[0117] In the embodiment of the present application, the preset error threshold refers to an error critical value set based on engineering precision requirements and reliability of measured data, including a first error threshold and a second error threshold.

[0118] If the first error value exceeds the first error threshold, it is determined that the error is mainly caused by insufficient calculation accuracy of wind load or deviation of vibration response simulation of the conductor 1; if the second error value exceeds the second error threshold, it is focused on deviation of boundary constraint or damping parameter of the finite element model; if both exceed the threshold, the coupling parameter of the fluid dynamics model and the finite element model is simultaneously associated. Then, a hierarchical calibration process is started based on the error source: for the case that the first error value exceeds the threshold, the turbulence intensity correction coefficient in the k-ω SST turbulence model of the fluid dynamics model is adjusted, and the wind load at each point of the conductor 1 is recalculated, while the damping matrix of the conductor unit in the finite element model is fine-tuned to suppress excessive vibration response; for the case that the second error value exceeds the threshold, the damping coefficient of the anti-hammer is corrected based on the dynamic tension increment error, and the stiffness parameter of the insulator unit is adjusted to optimize the boundary support force transmission characteristics. During the calibration process, the input of the meteorological mutation condition is kept stable, the coupled simulation is re-run after each parameter adjustment, and the new first error value and the second error value are calculated until both sets of errors fall within the preset threshold, and finally the calibrated coupled simulation model is output.

[0119] In step 103, the fatigue life of the conductor 1 of the power transmission line is evaluated based on the target simulation data corresponding to the calibrated coupled simulation model, and a fatigue life evaluation result is generated.

[0120] In the embodiment of the present application, the target simulation data refers to the simulation data output by the calibrated coupled simulation model and verified by the measured data.

[0121] The fatigue life evaluation result refers to a structured report integrating the damage distribution along the span of the conductor 1, the residual life at each point, the fatigue danger section and the evaluation accuracy.

[0122] The dynamic stress time history data of the conductor 1 at each point along the span is extracted from the target simulation data, which is derived from the displacement data calculated by the calibrated finite element model and has been verified by the measured vibration and tension data with an error meeting the requirements; at the same time, the S-N fatigue curve (conductor stress-life curve) of the conductor material is called to serve as the basis for fatigue life calculation. Then, the evaluation is carried out according to the Miner linear cumulative damage theory: first, the rainflow counting method is used for the dynamic stress time history to extract different stress amplitudes Δσ i and corresponding cycle numbers n i , which converts the complex stress sequence into discrete "stress amplitude-cycle number" data pairs; second, the fatigue life N i at each stress amplitude is calculated based on the S-N curve, that is, N i =C / Δσ m, and then the cumulative damage degree D of each point of the conductor 1 is calculated by formula D = ∑(n i / N i ), wherein D = 1 indicates fatigue failure; in the third step, the residual life is calculated based on the cumulative damage degree distribution, the damage rate D rate =D total / T total in the monitoring period is determined, wherein D total is the total damage in the period, T total is the monitoring period, and then formula T remain =(1-D max ) / D rate is used, wherein D max is the maximum damage degree along the span, so that the residual life of each point is obtained, and the life prediction result is corrected in combination with the 90% confidence interval, so that the influence of extreme working conditions is avoided; finally, the fatigue life evaluation result is generated by integrating the damage degree distribution along the span, the residual life and the dangerous section information.

[0123] Further, the step 103 comprises the following sub-steps:

[0124] S31, input the vibration data, tension data and environmental meteorological data of the conductor of the power transmission line into the calibrated coupled simulation model for simulation to generate new simulation data.

[0125] In the embodiment of the present application, the calibrated coupled simulation model refers to the coupled model corrected by the first error value and the second error value, which comprises the optimized fluid dynamics sub-model and the finite element coupled sub-model.

[0126] The new simulation data refers to the structured data set output by the calibrated model, which comprises multidimensional parameters such as wind load, displacement, stress, amplitude and frequency.

[0127] The preprocessed environmental meteorological data is converted into parameters such as wind speed, wind direction and air density according to the format required by the fluid dynamics model, and the wind field boundary condition is reconstructed in combination with the optimized turbulence intensity correction coefficient in the calibrated model; at the same time, the static tension in the preprocessed tension data is taken as the initial tension constraint of the model, and the dynamic tension increment is taken as the real-time verification benchmark in the simulation process, and the measured amplitude and frequency in the vibration data are taken together to form the verification dimension of the model input. After starting the simulation, the calibrated fluid dynamics model calculates the wind load of each point of the conductor 1 based on the optimized parameters, and the load mapping technology is used to transfer the wind load to the finite element coupled model, the model calls the calibrated damping matrix and stiffness matrix, and solves the displacement and stress response of the conductor vibration, and the deviation between the simulation amplitude and the measured amplitude and the deviation between the simulation dynamic tension and the measured dynamic tension are compared in real time during the process to ensure the stability of the simulation process; so that the new simulation data is obtained, which comprises the wind load, displacement time history, stress time history, amplitude and frequency aligned in time and space, and is accompanied by a data confidence label.

[0128] S32, jump to execute the step of calculating the first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data, until the first error value and the second error value are both less than or equal to the preset error threshold, and then the current simulation data is determined as the target simulation data.

[0129] In the embodiment of the application, jump to execute refers to returning to the previous error calculation and model calibration link when the error value does not meet the threshold requirement, and adjusting the parameters through the iterative adjustment and re-simulation cycle process.

[0130] The simulation amplitudes of each time step and each conductor position in the latest simulation data are matched with the measured amplitudes of the corresponding space-time points in the preprocessed vibration data, and the first error value is calculated according to the first error formula; at the same time, the dynamic component of the support reaction force and the dynamic tension increment data are calculated according to the second error formula to form two error value sequences. Then, the two error values are compared with the preset error threshold point by point, and if all the error values meet the threshold, the current simulation data is determined as the target simulation data; if any error value exceeds the threshold, return to step S21, adjust the parameters according to the calibration strategy for iterative adjustment, run the coupled simulation again to generate new simulation data, and repeat the error calculation and threshold comparison process. This cycle continues until both error values are less than or equal to the preset threshold, ensuring that the model outputs the target simulation data that meets the accuracy requirement in the stable state.

[0131] S33, extract the stress amplitude of the stress data of the target simulation data and the cycle number corresponding to the stress amplitude.

[0132] In the embodiment of the application, the stress amplitude refers to half of the difference between the maximum stress and the minimum stress in a complete cycle in the dynamic stress time history.

[0133] The cycle number refers to the total number of times a certain stress amplitude appears in the dynamic stress time history.

[0134] The dynamic stress time history curve of conductor 1 along the span is extracted from the target simulation data, the cycle of the time history curve is identified by the rain flow counting method, the stress time history is regarded as a fluctuating "waveform", the rainwater is simulated to flow along the waveform trajectory, the complete cycle of the stress from the peak value to the valley value is tracked, and the maximum stress value and the minimum stress value of each closed cycle are distinguished, and each stress amplitude is calculated by the maximum stress value and the minimum stress value, and the calculation formula is:

[0135] Stress amplitude = (maximum stress value - minimum stress value) / 2

[0136] Simultaneously, the number of times each stress amplitude occurs throughout the entire time history is counted to obtain the correspondence between "stress amplitude - number of cycles". During this period, small cycles with stress amplitudes less than the material fatigue limit are removed, and consecutively occurring identical stress amplitudes are merged to obtain the stress amplitude and corresponding number of cycles at each point along conductor 1.

[0137] S34. Based on the conductor stress-life curve, calculate the damage distribution of conductor 1 using the stress amplitude and the number of cycles corresponding to the stress amplitude.

[0138] In this embodiment of the invention, the conductor stress-life curve refers to a curve that describes the fatigue life relationship of the conductor material under different stress amplitudes.

[0139] Damage distribution refers to the spatial distribution of the total cumulative damage at each point along a span of one conductor.

[0140] The stress-life curve corresponding to the conductor model is determined by the conductor model. Then, each stress amplitude Δσ is retrieved one by one from the stress amplitude-cycle number data pairs. i and its corresponding number of iterations n i The fatigue life under this stress amplitude is calculated by combining the conductor stress-life curve. Subsequently, according to Miner's linear cumulative damage theory, through D... i =n i / N i Calculate the damage contribution of each stress amplitude to conductor 1, and sum the damage contributions of all stress amplitudes at the same location to obtain the total cumulative damage D at that location. total The total cumulative damage is calculated sequentially at each point along the first span of the conductor to form a continuous damage numerical sequence from the tower connection point to the midpoint of the span. Then, a damage distribution cloud map is generated by interpolation.

[0141] S35. Based on the damage distribution of conductor 1, determine the maximum damage degree of conductor 1.

[0142] In this embodiment of the invention, the maximum damage degree refers to the total cumulative damage degree with the largest value in the damage distribution.

[0143] The total cumulative damage D at each point along the first span of the conductor is obtained based on the damage distribution data. total By eliminating outliers caused by calculation errors, the remaining effective damage data are arranged in order of the spatial location of the conductor. The damage values ​​at each location are compared point by point, and the damage data with the largest value is identified as the maximum damage value of conductor 1.

[0144] S36. Using the maximum damage degree and preset damage rate of conductor 1, evaluate the fatigue life of conductor 1 and generate fatigue life evaluation results.

[0145] In the embodiment of the present application, the preset damage rate refers to the preset cumulative speed of fatigue damage of the conductor 1 per unit time.

[0146] Residual fatigue life = (1-maximum damage degree) / preset damage rate;

[0147] The conductor 1 can be calculated by the above formula to determine the safe operation time. If the maximum damage degree is 0.6 and the preset damage rate is 0.1 / year, the residual life is 4 years. In combination with the current operation period of the conductor 1, the residual proportion of the total design fatigue life can be further calculated, and finally the residual life value, the maximum damage degree position and the risk level (such as the residual life <5 years for high risk) are integrated to generate the fatigue life evaluation result including the evaluation basis, the calculation process and the recommended measures.

[0148] Step 104, under multiple constraint conditions, the preset vibration damper installation optimization model is iteratively solved, the vibration damper installation optimization scheme is determined according to the solving result, and the vibration damper installation optimization scheme and the fatigue life evaluation result are output to the evaluation platform.

[0149] In the embodiment of the present application, the constraint condition refers to the limiting condition that needs to be met in the vibration damper installation optimization process, including the installation quantity, the position range, the dynamic tension control and the cost, etc.

[0150] Iterative solution refers to the process of generating and verifying the installation scheme multiple times by an optimization algorithm.

[0151] The vibration damper installation optimization scheme refers to the optimal scheme determined after iterative solution, including the installation position, quantity, model and optimization effect data of the vibration damper.

[0152] The evaluation platform refers to a visualization platform for displaying the vibration damper optimization scheme and the fatigue life evaluation result.

[0153] The upper limit of the damper installation quantity, the installation position range, the dynamic tension control threshold and the cost constraint are constructed into multiple constraint conditions, and the above conditions are set as the boundary limit parameters of the model; the conductor damage distribution, the fatigue life evaluation result and the dynamic parameters of the damper are input into the preset damper installation optimization model, and the iterative solution is started, the initial installation scheme is generated by using the genetic algorithm optimization, the optimized damage distribution and tension response are calculated by substituting into the calibrated coupled simulation model, whether all constraint conditions are satisfied is checked, if not, the scheme is adjusted and regenerated, if satisfied, the objective function value is calculated, and the scheme with the highest objective function value is screened out through multiple rounds of iteration; the damper installation optimization scheme is determined according to the solving result, wherein the damper installation optimization scheme includes the specific installation position, model and quantity of each group of dampers, and the damage degree comparison data before and after optimization is attached; finally, the damper installation optimization scheme and the fatigue life evaluation result are output to the evaluation platform. The damper installation position 3 is provided with a damper.

[0154] Further, step 104 includes the following sub-steps:

[0155] S41, based on the distance between the clamp of the conductor 1 and the damper installation position 3, a first constraint condition is constructed.

[0156] In the embodiment of the application, the first constraint condition refers to the constraint set based on the distance between the clamp and the damper installation position 3.

[0157] The minimum distance between the clamp and the damper installation position 3 is not less than 1.5m, and the maximum distance is not more than 5m, and the first constraint condition is constructed.

[0158] S42, according to the distance between each two adjacent dampers, a second constraint condition is constructed.

[0159] In the embodiment of the application, the second constraint condition refers to the constraint set based on the distance between each two adjacent dampers.

[0160] According to the wavelength characteristics of the conductor 1 vibration, the reasonable distance range of adjacent dampers is determined, the minimum distance is set to be not less than 5m, the maximum distance is set to be not more than 15m, if the number of dampers in the span exceeds 2 groups, the difference between all adjacent distances is additionally constrained to be not more than 3m, that is, 5m≤L≤15m, and |L i -L j |≤3m, the second constraint condition is constructed, wherein L is the distance between each two adjacent dampers, L i , and L j are different adjacent distances.

[0161] S43, according to the number of dampers, a third constraint condition is constructed.

[0162] In the embodiment of the present application, the third constraint condition refers to a constraint set based on the number of damping hammers.

[0163] The number of damping hammers is determined according to the span length of the conductor 1, historical damping effect data and cost budget. If the span length is less than or equal to 200 m, the lower limit of the number of damping hammers is set to 1 group and the upper limit is set to 3 groups. If the span length is greater than 200 m, the lower limit of the number of damping hammers is adjusted to 2 groups and the upper limit is adjusted to 5 groups. At the same time, the matching relationship between the number of damping hammers and the damage distribution of the conductor 1 needs to be additionally constrained, that is, at least one group of damping hammers is configured in each 30 m range of the high damage section, and the third constraint condition is constructed according to the above number requirement.

[0164] S44, under the first constraint condition, the second constraint condition and the third constraint condition, an optimization model for installing damping hammers is constructed with the minimization of the unevenness of the damage distribution of the conductor 1 as the target.

[0165] In the embodiment of the present application, the minimization of the unevenness of the damage distribution refers to an index quantifying the damage difference of each section of the conductor 1, which is calculated by the difference between the maximum damage degree and the average damage degree.

[0166] The unevenness of the damage distribution is taken as the calculation index, and the calculation formula of the unevenness of the damage distribution is:

[0167] The unevenness of the damage distribution = (the maximum damage degree - the average damage degree) / the average damage degree * 100%.

[0168] The smaller the value of the unevenness of the damage distribution obtained in the above formula represents the smaller the damage difference of each section of the conductor 1, which avoids local excessive damage. Then the decision variable of the model is determined as the installation position coordinates and the number of damping hammers, and the first constraint condition, the second constraint condition and the third constraint condition are converted into boundary restrictions of the decision variable. Then the objective function is constructed with the minimization of the unevenness of the damage distribution, and the calibrated coupled simulation model module is embedded to obtain the optimization model for installing damping hammers.

[0169] S45, the optimization model for installing damping hammers is iteratively solved by using a genetic algorithm to obtain the optimal number of damping hammers and the installation positions of the damping hammers.

[0170] In the embodiment of the present application, the genetic algorithm refers to an optimization algorithm simulating the biological evolution process, which generates an initial population by coding, iteratively optimizes through selection, crossover and mutation operations, and finally finds the optimal solution of the objective function.

[0171] The optimal number and the installation positions of the damping hammers refer to the configuration parameters of the damping hammers obtained after iteration convergence.

[0172] The number and installation position of the damper are coded as a chromosome, and an initial population is generated based on the third constraint condition; then, the fitness value of each individual is calculated, that is, the unevenness of damage distribution is output by the damper installation optimization model, and "1 / unevenness" is taken as the fitness value, while the individual is checked whether it meets the first constraint condition and the second constraint condition, and the individual that does not meet the constraint is given a very low fitness value; the iteration link is entered, the individual with high fitness value is selected by roulette selection, the offspring population is generated by using single-point crossover and random mutation, and the optimal individual with a fitness value in the top 10% of the parent population is reserved to maintain excellent characteristics; the selection, crossover and mutation process is repeated until the population fitness tends to be stable, at this time, the individual solution with the highest fitness value is decoded to obtain the optimal number of dampers and the specific installation position of each damper.

[0173] S46, in combination with the optimal number of dampers and the installation position of each damper 3, an optimal damper installation scheme is generated.

[0174] In the embodiment of the present application, the optimal number is associated with the corresponding installation position, and the specific coordinates and damper model of each position are labeled; then, it is verified whether the combination meets the first, second and third constraint conditions, and the comparison data of damage distribution before and after optimization is attached; then, the installation construction points are supplemented, including the clamp fastening torque, the direction of the damper and the spacing requirement from other fittings, and finally the optimal damper installation scheme is obtained, which includes the configuration parameters, the constraint verification result, the damping effect and the construction specification.

[0175] S47, the optimal damper installation scheme and the fatigue life evaluation are output to the evaluation platform.

[0176] In the embodiment of the present application, the optimal damper installation scheme and the fatigue life evaluation are transmitted to the evaluation platform, and the platform automatically analyzes and classifies the display after receiving, wherein the optimal damper installation scheme is visualized in a three-dimensional span model, and the fatigue life data is displayed in the form of a trend chart, and at the same time, the platform generates a comprehensive evaluation report, including the scheme compliance judgment, the life improvement percentage and the operation and maintenance priority suggestion.

[0177] Please refer to Figure 4 , Figure 4 A structural block diagram of a power transmission line aeolian vibration dynamic evaluation system provided in Embodiment Two of the present application.

[0178] The power transmission line aeolian vibration dynamic evaluation system provided by the present application comprises:

[0179] The simulation module 201 is configured to input the vibration data, tension data and environmental meteorological data of the conductor 1 of the power transmission line into a preset initial coupling simulation model for simulation, and output simulation data.

[0180] The calibration module 202 is configured to calibrate the initial coupled simulation model according to error values between the simulation data and the vibration data and the tension data respectively under the meteorological mutation condition of the environmental meteorological data;

[0181] The evaluation module 203 is configured to perform fatigue life evaluation on the conductor 1 of the power transmission line based on the target simulation data corresponding to the calibrated coupled simulation model, and generate a fatigue life evaluation result;

[0182] The output module 204 is configured to iteratively solve the preset damper installation optimization model under a plurality of constraint conditions, determine a damper installation optimization scheme according to a solution result, and output the damper installation optimization scheme and the fatigue life evaluation result to an evaluation platform.

[0183] Further, the simulation module 201 comprises:

[0184] The acquisition sub-module is configured to acquire vibration data, tension data and environmental meteorological data of the conductor 1 of the power transmission line;

[0185] The preprocessing sub-module is configured to preprocess the vibration data, the tension data and the environmental meteorological data;

[0186] The construction sub-module is configured to construct an initial coupled simulation model based on the preprocessed vibration data, tension data and environmental meteorological data; wherein the initial coupled simulation model comprises a fluid dynamics model and a finite element model;

[0187] The first simulation sub-module is configured to input the preprocessed vibration data, tension data and environmental meteorological data into the initial coupled simulation model for simulation, and output point wind load, displacement data, stress data, simulation amplitude and simulation frequency of the conductor 1;

[0188] The combination sub-module is configured to combine the point wind load, displacement data, stress data, simulation amplitude and simulation frequency of the conductor 1 to generate simulation data.

[0189] Further, the combination sub-module comprises:

[0190] The input sub-module is configured to input the preprocessed environmental meteorological data into the fluid dynamics model of the initial coupled simulation model for simulation, and output point wind load of the conductor 1;

[0191] The output sub-module is configured to input the point wind load of the conductor 1, the preprocessed vibration data and the tension data into the finite element coupled model of the initial coupled simulation model for simulation, and output displacement data, stress data, simulation amplitude and simulation frequency.

[0192] Further, the calibration module 202 comprises:

[0193] a first calculation submodule, configured to calculate a first error value between a simulation amplitude of the simulation data and a measured amplitude of the vibration data;

[0194] a dynamic component submodule, configured to determine a dynamic component of the reaction force of the finite element model of the initial coupled simulation model based on displacement data of the simulation data and the mechanical properties;

[0195] a second calculation submodule, configured to calculate a second error value between the dynamic component of the reaction force and a dynamic tension increment of the tension data;

[0196] a calibration submodule, configured to calibrate the initial coupled simulation model when the first error value or the second error value is greater than a preset error threshold and a weather of the environmental weather data is in a weather mutation condition.

[0197] Further, the evaluation module 203 comprises:

[0198] a second simulation submodule, configured to input the vibration data, the tension data and the environmental weather data of the conductor of the power transmission line into the calibrated coupled simulation model to perform simulation, and generate new simulation data;

[0199] a jump submodule, configured to jump to execute the step of calculating the first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data until the first error value and the second error value are both less than or equal to the preset error threshold, and then determine the current simulation data as target simulation data;

[0200] an extraction submodule, configured to extract a stress amplitude of the stress data of the target simulation data and a cycle number corresponding to the stress amplitude;

[0201] a third calculation submodule, configured to calculate a damage distribution of the conductor 1 based on a conductor stress-life curve and the stress amplitude and the cycle number corresponding to the stress amplitude;

[0202] a damage submodule, configured to determine a maximum damage degree of the conductor 1 based on the damage distribution of the conductor 1;

[0203] an evaluation submodule, configured to evaluate a fatigue life of the conductor 1 based on the maximum damage degree of the conductor 1 and a preset damage rate, and generate a fatigue life evaluation result.

[0204] Further, the output module 204 comprises:

[0205] a first constraint submodule, configured to construct a first constraint condition based on a distance between the clamp of the conductor 1 and the installation position 3 of the damper;

[0206] a second constraint submodule, configured to construct a second constraint condition according to distances between two adjacent dampers;

[0207] The third constraint sub-module is configured to construct a third constraint condition according to the number of the anti-vibration hammers;

[0208] The model construction sub-module is configured to construct an anti-vibration hammer installation optimization model under the first constraint condition, the second constraint condition and the third constraint condition, with the minimum damage distribution unevenness of the conductor 1 as the target.

[0209] The iteration sub-module is configured to solve the anti-vibration hammer installation optimization model by using a genetic algorithm to obtain the optimal number of the anti-vibration hammers and the respective anti-vibration hammer installation positions 3.

[0210] The optimization sub-module is configured to generate an anti-vibration hammer installation optimization scheme in combination with the optimal number of the anti-vibration hammers and the respective anti-vibration hammer installation positions 3.

[0211] The output platform sub-module is configured to output the anti-vibration hammer installation optimization scheme and the fatigue life assessment to an evaluation platform.

[0212] Please refer to Figure 5 , Figure 5 A structural block diagram of a computer device provided for the third embodiment of the present application.

[0213] The electronic device of the embodiment of the present application comprises a memory 301 and a processor 302, the memory 301 stores a computer program, and the computer program is executed by the processor 302 to make the processor 302 execute the power transmission line breeze vibration dynamic evaluation method of any one of the above embodiments.

[0214] The memory 301 can be an electronic memory such as a flash memory, an EEPROM (Electrically Erasable Programmable Read-Only Memory), an EPROM, a hard disk or a ROM. The memory 301 has a storage space 303 for program codes 313 for performing any of the method steps in the above described methods. For example, the storage space 303 for program codes can comprise individual program codes 313 for implementing the various steps in the above described methods, respectively. These program codes can be read from or written to one or more computer program products. These computer program products include program code carriers such as a hard disk, a compact disk (CD), a memory card or a floppy disk. The program codes can be compressed, for example, in a suitable form. These codes, when run by a computing processing device, cause the computing processing device to perform the individual steps in the above described methods. These program codes can be read from or written to one or more computer program products. These computer program products include program code carriers such as a hard disk, a compact disk (CD), a memory card or a floppy disk. The program codes can be compressed, for example, in a suitable form. These codes, when run by a computing processing device, cause the computing processing device to perform the individual steps in the above described power line aeolian vibration dynamic assessment method.

[0215] Embodiment four of the present application provides a computer readable storage medium, having stored thereon a computer program, which when executed by a computer implements the power line aeolian vibration dynamic assessment method of any embodiment of the present application.

[0216] Embodiment five of the present application provides a computer program product, comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions, wherein the program instructions, when executed by a computer, cause the computer to perform the power line aeolian vibration dynamic assessment method of any embodiment described above.

[0217] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above described system, device and unit can refer to the corresponding processes in the foregoing method embodiments, which will not be described herein.

[0218] In several embodiments provided in the present application, it should be understood that the disclosed system, device and method can be implemented in other manners. For example, the described device embodiments are merely schematic. For example, the division of the units is only a logical function division. There can be another division manner for the actual implementation, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between different units, can be indirect couplings or communication connections through some interfaces, devices or units, and can be electrical, mechanical or in other forms.

[0219] The units described as separate components can or can not be physically separate, and the components shown as units can or can not be physical units, i.e., can be located in one place, or can be distributed on multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the embodiment.

[0220] In addition, each functional unit in each embodiment of the present application can be integrated into a processing unit, or each unit can exist physically, or two or more units can be integrated into one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0221] The integrated unit, if realized in the form of a software functional unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application, essentially or the part that contributes to the prior art, or all or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a number of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0222] The above embodiments are merely used to describe the technical solutions of the present application, rather than limit them. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that the technical solutions recorded in the foregoing embodiments can be modified, or some technical features can be replaced by equivalent ones; and these modifications or replacements do not cause the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for dynamic assessment of aeolian vibrations of a power transmission line, characterized in that, The method comprises the following steps: inputting vibration data, tension data and environmental meteorological data of a conductor of a power transmission line into a preset initial coupling simulation model for simulation, and outputting simulation data; under the condition of a meteorological mutation of the environmental meteorological data, calibrating the initial coupling simulation model according to error values between the simulation data and the vibration data and the tension data respectively; based on target simulation data corresponding to the calibrated coupling simulation model, performing fatigue life assessment on the conductor of the power transmission line, and generating fatigue life assessment results; under multiple constraint conditions, iteratively solving a preset anti-vibration hammer installation optimization model, determining an anti-vibration hammer installation optimization scheme according to a solution result, and outputting the anti-vibration hammer installation optimization scheme and the fatigue life assessment results to an evaluation platform.

2. The method for galloping dynamic assessment of a power transmission line according to claim 1, characterized in that, The method comprises the following steps: obtaining vibration data, tension data and environmental meteorological data of a conductor of a power transmission line; preprocessing the vibration data, the tension data and the environmental meteorological data; based on the preprocessed vibration data, tension data and environmental meteorological data, constructing an initial coupling simulation model; wherein the initial coupling simulation model comprises a fluid dynamics model and a finite element model; inputting the preprocessed vibration data, tension data and environmental meteorological data into the initial coupling simulation model for simulation, and outputting point wind load, displacement data, stress data, simulation amplitude and simulation frequency of the conductor; combining the point wind load of the conductor, the displacement data, the stress data, the simulation amplitude and the simulation frequency, and generating simulation data.

3. The galloping dynamic assessment method of a power transmission line according to claim 2, characterized in that, The method comprises the following steps: inputting the preprocessed environmental meteorological data into the fluid dynamics model of the initial coupling simulation model for simulation, and outputting point wind load of the conductor; inputting the point wind load of the conductor, the preprocessed vibration data and tension data into the finite element coupling model of the initial coupling simulation model for simulation, and outputting the displacement data, the stress data, the simulation amplitude and the simulation frequency.

4. The method for galloping dynamic assessment of a power transmission line according to claim 1, characterized in that, The method comprises the following steps: calculating a first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data; based on the displacement data and the mechanical properties of the simulation data, determining a dynamic component of a support reaction force of the finite element model of the initial coupling simulation model; calculating a second error value between the dynamic component of the support reaction force and a dynamic tension increment of the tension data; When the first error value or the second error value is greater than a preset error threshold, and a weather of the environmental weather data is in a weather mutation condition, the initial coupled simulation model is calibrated.

5. The method for galloping dynamic assessment of a power transmission line according to claim 4, characterised in that, The fatigue life of the conductor of the power transmission line is evaluated based on the target simulation data corresponding to the calibrated coupled simulation model, and a fatigue life evaluation result is generated, including: The vibration data, tension data and environmental weather data of the conductor of the power transmission line are input into the calibrated coupled simulation model for simulation, and new simulation data are generated; The step of calculating the first error value between the simulation amplitude of the simulation data and the measured amplitude of the vibration data is jumped to be executed until the first error value and the second error value are both less than or equal to the preset error threshold, and the current simulation data is determined as the target simulation data; The stress amplitude of the stress data of the target simulation data and the cycle number corresponding to the stress amplitude are extracted; Based on the conductor stress-life curve, the stress amplitude and the cycle number corresponding to the stress amplitude are used to calculate the damage distribution of the conductor; Based on the damage distribution of the conductor, the maximum damage degree of the conductor is determined; The fatigue life of the conductor is evaluated based on the maximum damage degree of the conductor and a preset damage rate, and a fatigue life evaluation result is generated.

6. The method for galloping dynamic assessment of a power transmission line according to claim 1, characterised in that, The preset damper installation optimization model is iteratively solved under multiple constraint conditions, a damper installation optimization scheme is determined according to a solving result, and the damper installation optimization scheme and the fatigue life evaluation result are output to an evaluation platform, including: A first constraint condition is constructed based on the distance between the clamp of the conductor and the damper installation position; A second constraint condition is constructed according to the distance between two adjacent dampers; A third constraint condition is constructed according to the number of dampers; A damper installation optimization model is constructed with the minimum damage distribution unevenness of the conductor as the target under the first constraint condition, the second constraint condition and the third constraint condition; A genetic algorithm is used to iteratively solve the damper installation optimization model to obtain the optimal number of dampers and the installation position of each damper; A damper installation optimization scheme is generated in combination with the optimal number of dampers and the installation position of each damper; The damper installation optimization scheme and the fatigue life evaluation are output to the evaluation platform.

7. A dynamic evaluation system for micro-wind vibration of transmission lines, characterized in that, including: A simulation module is configured to input the vibration data, tension data and environmental weather data of the conductor of the power transmission line into a preset initial coupled simulation model for simulation, and output simulation data; A calibration module is configured to calibrate the initial coupled simulation model according to the error values between the simulation data and the vibration data and the tension data under the weather mutation condition of the environmental weather data; An evaluation module is configured to evaluate the fatigue life of the conductor of the power transmission line based on the target simulation data corresponding to the calibrated coupled simulation model, and generate a fatigue life evaluation result. An output module is configured to iteratively solve the preset installation optimization model of the anti-vibration hammer under multiple constraint conditions, determine an installation optimization scheme of the anti-vibration hammer according to a solving result, and output the installation optimization scheme of the anti-vibration hammer and the fatigue life evaluation result to an evaluation platform.

8. An electronic device, comprising: The computer program is stored in the memory and executed by the processor, and the processor executes the method for dynamic evaluation of aeolian vibration of a power transmission line according to any one of claims 1-6.

9. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed to implement the method for dynamic evaluation of aeolian vibration of a power transmission line according to any one of claims 1-6.

10. A computer program product, characterised in that, The computer program product comprises a computer program stored on a non-transitory computer-readable storage medium, and the computer program comprises program instructions, wherein when the program instructions are executed by a computer, the computer executes the method for dynamic evaluation of aeolian vibration of a power transmission line according to any one of claims 1-6.

Images