CN111651805A - Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors - Google Patents

Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors Download PDF

Info

Publication number
CN111651805A
CN111651805A CN202010245692.1A CN202010245692A CN111651805A CN 111651805 A CN111651805 A CN 111651805A CN 202010245692 A CN202010245692 A CN 202010245692A CN 111651805 A CN111651805 A CN 111651805A
Authority
CN
China
Prior art keywords
wind
insulator string
load
lead
wire
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010245692.1A
Other languages
Chinese (zh)
Other versions
CN111651805B (en
Inventor
赵爽
晏致涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN202010245692.1A priority Critical patent/CN111651805B/en
Publication of CN111651805A publication Critical patent/CN111651805A/en
Application granted granted Critical
Publication of CN111651805B publication Critical patent/CN111651805B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08—Construction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Geometry (AREA)
  • Human Resources & Organizations (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Insulators (AREA)

Abstract

The invention discloses a calculation method for a maximum wind deflection angle and a wind vibration coefficient of a suspended insulator string by considering linear shape and linear length influence factors, which particularly determines the type of a power transmission tower and the arrangement scheme of the power transmission tower, a lead and the insulator string, and obtains physical parameters of the lead and physical parameters of an insulator on the power transmission tower; determining a calculation model of the wind deflection angle of the suspension insulator string by an LRC method by taking the gravity and the average wind load as initial conditions for calculating the lead and the suspension insulator string; calculating to obtain the equivalent static wind load and the component data of the conductor in unit area between the power transmission towers, and correspondingly obtaining the equivalent static wind load and the component distribution diagram of the conductor in unit area between the power transmission towers; calculating the maximum wind deflection angle of the suspension insulator string; and calculating the wind vibration coefficient of the suspension insulator string. Has the advantages that: the method can provide reference for revising the wind load of the wind deflection angle design of the suspension insulator string of the power transmission line, so that the design precision is high.

Description

Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors
Technical Field
The invention relates to the technical field of power transmission tower design, in particular to a calculation method for a maximum wind deflection angle and a wind vibration coefficient of a suspension insulator string by considering linear and line length influence factors.
Background
The windage yaw flashover accident is one of important disasters of a power transmission line, and the windage yaw flashover accident of an insulator string is easily induced by a strong wind load, so that the damage of a power transmission tower component and (or) the occurrence of a large-scale power failure accident can be caused. In engineering design, the wind deflection angle is accurately determined, so that the safe operation of the power transmission line can be ensured, and the construction cost of the power transmission line can be reduced.
In calculating the wind deflection angle of the suspension insulator string, the rigid straight rod method is adopted in the current electric power industry standard in China (for concrete documents, see Zhangsheng, A handbook of design of high-voltage power transmission lines in electric power engineering (2 nd edition) [ M ]. Beijing: Chinese electric power publisher, 2002:66-91 and Shao Tianxiao, and calculation of electric wire mechanics of overhead power transmission lines [ M ]. Beijing: Chinese electric power publisher, 2003: 33-99.). The rigid body straight rod method is used for calculating the wind deflection angle of the suspension insulator string under the action of average wind force through a statics method. The wind deflection angle of the suspension insulator string under the action of fluctuating wind load is calculated on a time domain through a finite element model, and the result shows that the traditional statics method is unsafe (specifically, see the documents: Kondeyi, Li, Longxiahong, and the like. dynamic wind deflection angle finite element analysis of the suspension insulator string [ J ]. electric power construction, 2008,29(9): 5-9.). Because the wind deflection angle of the suspension insulator string is not completely accurately calculated by the existing overhead transmission line design rule, the values of the wind load adjustment coefficients of the transmission line under the conditions of different span, different height difference and different ground-to-ground heights are calculated in the time domain through a finite element model, and the research result provides reference for the revision of the existing design rule. Some scholars also give wind load calculation formulas suitable for design use by root mean square (root mean square) response spectrum expressions of dynamic tension. The traditional rigid body straight rod method does not consider the influences of shielding and pulsating wind between split conductors, analyzes the influences of a hang point height difference, shielding wake flow influence and power amplification action by establishing a finite element model of an insulator string-conductor system, and provides a correction formula of the rigid body straight rod method. The correction formula has reference value and guiding significance for line design and windage yaw prevention, and is particularly described in a reference (Li, Showns, Luo Xian, and the like. windage yaw calculation method for a super high voltage insulator string [ J ] high voltage technology, 2013,39(12): 2924-. Accurate wind load is the premise of effectively calculating the wind deflection angle of the suspension insulator string. Compared with the wind load standard of domestic and foreign power transmission lines, when the wind load of the wire is calculated, the product of the wind pressure non-uniform coefficient and the wind load adjustment coefficient which are standardized in China takes an experience coefficient smaller than 1, and the method is different from the foreign standard and has some defects in consideration of the pulsating wind effect. The dynamic windage yaw time domain calculation is carried out on the continuous multi-span power transmission line, and the result shows that the dynamic amplification effect of the pulsating wind in the actual windage yaw calculation process cannot be ignored. The aerodynamic damping of the wire when wind vibration occurs is increased along with the increase of the average wind speed, the resonance component of the wind vibration response is greatly reduced due to the aerodynamic damping, and the wind vibration response can be ignored in the calculation.
Disclosure of Invention
Aiming at the problems, the invention provides a method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering linear shape and linear length influence factors.
In order to achieve the purpose, the invention adopts the following specific technical scheme:
a method for calculating the maximum wind deflection angle and the wind vibration coefficient of a suspension insulator string by considering linear shape and linear length influence factors has the key technology that the method comprises the following steps:
s1: determining the type of the power transmission tower and the arrangement scheme of the power transmission tower, the lead and the insulator string, and acquiring physical parameters of the lead and physical parameters of the insulator on the power transmission tower;
the physical parameters of the lead at least comprise the linearity and the length of the lead;
s2: determining a calculation model of the wind deflection angle of the suspension insulator string by an LRC method by taking the gravity and the average wind load as initial conditions for calculating the lead and the suspension insulator string;
s3: combining the data of the step S1 and the data of the step S2, calculating to obtain the equivalent static wind load and the component data of the conductor in unit area between the transmission towers, and correspondingly obtaining the equivalent static wind load and the component distribution diagram of the conductor in unit area between the transmission towers;
s4: calculating the maximum wind drift angle of the suspension insulator string according to the equivalent static wind load of the step S3;
s5: and calculating the wind vibration coefficient of the suspension insulator string according to the equivalent static wind load of the step S3.
According to a further technical scheme, the transmission tower types at least comprise a conventional transmission tower and an ultrahigh transmission tower;
the physical parameters of the lead at least comprise the type of the lead, the calculated sectional area of the lead, the elastic modulus of the lead, the linear density and the outer diameter of the lead;
and the physical parameters of the insulator string on the power transmission tower at least comprise the length of the insulator string, the elastic modulus of the insulator string, the quality of the insulator string and the wind shielding area of the insulator string.
In a further technical scheme, in the step S2, when the model of calculating the wind deflection angle of the suspension insulator string is described, a connection point a between the lead and the insulator string, a tail end point B of the suspension insulator string, a motion point B 'of the tail end point of the insulator string in the dynamic state, and a wind deflection angle caused by moving the point B' to the point B "in the dynamic state are set
Figure RE-RE-GDA0002610855590000031
Downwind displacement of point B under action of wire span L and average wind load
Figure RE-RE-GDA0002610855590000032
A. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angle
Figure RE-RE-GDA0002610855590000033
The horizontal distance a 'from the origin of coordinates to the lowest point of the wire, and the horizontal distance b' from the lowest point of the wire to the end of the wire.
The wire is in a catenary configuration under a self-weight state, and shows large geometric deformation under the action of wind load. Previous researches show that the influence of the power transmission tower on the wind vibration response of the wire is small. In order to simplify the calculation, the influence of the tower is ignored, and the hanging point of the insulator on the tower is taken as a fixed hinged support, so that the wind deflection angle research is carried out on the hanging wire suspension insulator string.
In a further technical scheme, in step S3, the equivalent static wind load p per unit area of the conducting wires between the transmission towersESWLIs calculated by the formula：
Figure RE-RE-GDA0002610855590000041
Wherein (: i) represents the ith column element of the matrix;
Figure RE-RE-GDA0002610855590000042
equivalent background wind pressure;
Figure RE-RE-GDA0002610855590000043
the average wind load is obtained;
the incoming wind load can be decomposed into an average wind load and a fluctuating wind load, and the vibration equation matrix expression of the wire under the action of the wind load is as follows:
Figure RE-RE-GDA0002610855590000044
in the formula,
Figure RE-RE-GDA0002610855590000045
y' is the acceleration, the speed and the displacement of the wire node along the wind direction under the action of the pulsating wind load respectively;
Figure RE-RE-GDA0002610855590000046
the displacement of the lead joint along the wind direction under the action of average wind load.
The wire is a light flexible structure, and the performance under the strong wind load is as follows: 1) the structure is greatly deformed, and the geometric nonlinearity is obvious; 2) the structure stress and the displacement do not have a linear relation; 3) under the action of dynamic load, the structure is time-varying rigidity. Therefore, the above equation is a variable coefficient differential equation, and cannot be solved by using the linear superposition principle. The wind vibration response of the wire caused by the incoming wind load can be decomposed into an average response and a pulse response.
M is a quality matrix; c is a damping matrix; a K stiffness matrix; l issIs a node dependent area matrix;
vibration equation matrix expression of lead under action of fluctuating wind loadComprises the following steps:
Figure RE-RE-GDA0002610855590000047
according to the scheme, an LRC method is adopted to achieve equivalent static wind load. The wire is a light flexible structure, and the performance under the strong wind load is as follows: 1) the structure is greatly deformed, and the geometric nonlinearity is obvious; 2) the structure stress and the displacement do not have a linear relation; 3) under the action of dynamic load, the structure is time-varying rigidity.
Therefore, the vibration equation matrix expression of the wire under the action of the pulsating wind load can be obtained by taking the average wind state of the wire as the initial calculation condition.
The vibration equation matrix expression of the lead under the action of wind load cannot be solved by adopting a linear superposition principle. The wind vibration response of the wire caused by the incoming wind load can be decomposed into an average response and a pulse response.
In a further technical scheme, a calculation formula for calculating the maximum wind deflection angle of the suspension insulator string by the equivalent static wind load is as follows:
Figure RE-RE-GDA0002610855590000051
in the formula,
Figure RE-RE-GDA0002610855590000052
is the downwind peak displacement of the point B under the action of fluctuating wind load
Figure RE-RE-GDA0002610855590000053
lABA, B is the length of the insulator string between two points;
Figure RE-RE-GDA0002610855590000054
Figure RE-RE-GDA0002610855590000055
is the downwind displacement of the point B under the action of average wind load,
Figure RE-RE-GDA0002610855590000056
Figure RE-RE-GDA0002610855590000057
is the average wind deflection angle; the specific calculation formula is as follows:
Figure RE-RE-GDA0002610855590000058
Figure RE-RE-GDA0002610855590000059
Gvrespectively taking the average wind load and the vertical gravity load of the suspension insulator string at the target point;
Figure RE-RE-GDA00026108555900000510
Wvrespectively the average wind load and the vertical load transferred to the suspension insulator string by the lead at the target point.
The further technical scheme is as follows: the wire at the target point transmits the average wind load to the suspension insulator string
Figure RE-RE-GDA00026108555900000511
The calculation formula of (2) is as follows:
Figure RE-RE-GDA00026108555900000512
Figure RE-RE-GDA00026108555900000513
in the formula, NcThe number of the split conductors;
Figure RE-RE-GDA00026108555900000514
the uniform average wind load of the unit wire length of a single wire is obtained;hthe calculation mode is a pair formula for the line length of the lead in the horizontal span
Figure RE-RE-GDA00026108555900000515
Performing curve integration at a horizontal span; wherein,
Figure RE-RE-GDA00026108555900000516
in the formula,
Figure RE-RE-GDA00026108555900000517
is the load p' and the response yBThe correlation coefficient of (a);
Figure RE-RE-GDA00026108555900000518
is a response y in the initial conditionBThe influence line of (2).
According to a further technical scheme, when the geometric linear shapes of the left and right cross wires of the target point have points with the slope of 0, the wires transmit vertical loads W to the suspension insulator string according to stress balancevEqual to the total weight of the conductors in the target point vertical span, i.e., the vertical load W transmitted by the conductors at the target point to the suspended insulator string when the transmission tower is a conventional transmission towervThe calculation formula of (2) is as follows: wv＝NcPv v；
Wherein, PvThe gravity of the unit wire length of a single wire;vthe specific calculation formula is that the length of the conducting wire in the vertical span is as follows: within vertical span pair formula
Figure RE-RE-GDA00026108555900000519
Performing curve integration;
in addition, the slope of the geometric linear shape of a certain crossover wire in the left and right spans of the target point is not 0 everywhere, and according to the stress balance, the vertical load transmitted to the suspension insulator string by the crossover wire is equal to the total weight of the crossover wire plus the vertical component of the tension at the lowest point of the wire.
When the power transmission tower is an ultrahigh power transmission tower, the lead at the target point transmits a vertical load W to the suspension insulator stringvThe calculation formula of (2) is as follows: wv＝Pv l+Tvl+Pv r+Tvr；
Wherein,l、rrespectively calculating the lengths of the left span and the right span of the target point; t isvl、TvrAre respectively eyesVertical components of tension at the lowest points of the left and right span wires of the punctuation;
when the slope of the geometric line shape of the wire at a certain point across the wire is 0:
Figure RE-RE-GDA0002610855590000061
Tvl＝0；
when the slope of the wire at the geometrical line within the span is not 0:
Figure RE-RE-GDA0002610855590000062
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw＝σo4Ac；
Wherein,
Figure RE-RE-GDA0002610855590000063
in the formula, subscripts "3" and "4" represent a no-wind state and an average wind state, respectively; a. thecThe stress area of the lead is defined; ecIs the modulus of elasticity of the wire; gamma raycIs the comprehensive specific load of the lead wires,
Figure RE-RE-GDA0002610855590000064
γwin order to obtain the average wind pressure specific load,
Figure RE-RE-GDA0002610855590000065
Figure RE-RE-GDA0002610855590000066
the calculation formula is the average wind load of the unit line length of the lead:
Figure RE-RE-GDA0002610855590000067
lrfor representing span βrTo represent the altitude difference angle.
The further technical scheme is as follows: the equivalenceThe calculation formula for calculating the wind vibration coefficient β of the suspension insulator string by the static wind load is as follows:
Figure RE-RE-GDA0002610855590000071
∑Crepresenting summing elements within a computational domain;ccalculating the line length of the wire in the domain;
Figure RE-RE-GDA0002610855590000072
the average wind load is obtained;
Figure RE-RE-GDA0002610855590000073
equivalent background wind pressure.
And when the height difference between the target point and the adjacent tower wire hanging point is 0, selecting the target point horizontal span as the calculation domain. When the height difference exists, the equivalent static wind load at the position of the target point is more convex, so that the calculation domain is spanned by the left and right sides 1/4 of the selected target point. The calculation of the lines with different engineering design parameters shows that the calculation method of the consistent beta meets the engineering precision requirement. The determined beta is comprehensive in consideration, clear in mechanical significance and more reasonable.
The invention has the beneficial effects that: the invention determines the equivalent static wind load distribution of the wind deflection angle of the suspension insulator string by an LRC method by taking the gravity and the average wind load as initial conditions for the calculation of the lead and the suspension insulator string. And calculating an average wind deflection angle (in a large deformation and nonlinear state) by adopting a rigid body straight rod method model considering linear and linear length influences, and calculating a peak value pulsating wind deflection angle (in a small deformation and approximate linear state) by adopting a linear theory so as to obtain a calculation formula of the maximum wind deflection angle. And deducing a wind vibration coefficient formula of the wind deflection angle of the suspension insulator string according to the physical significance. And providing a calculation model of the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string. And then, determining the maximum wind drift angle and the wind vibration coefficient of the suspension insulator of the conventional power transmission line and the ultra-high large-span power transmission line by adopting a calculation model, comparing the maximum wind drift angle and the wind vibration coefficient with a time domain calculation result and a wind tunnel test measurement result, and verifying the reliability of the model. And finally, researching the influence of the change of engineering parameters (wind speed, altitude difference, breath height, vertical span ratio, span and landform roughness category) on the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by using the calculation model. The research of the invention can provide reference for revising the wind load of the wind deflection angle design of the suspension insulator string of the power transmission line.
Drawings
FIG. 1 is a computational flow diagram of the present invention;
FIG. 2 is a schematic view of a wind deflection angle calculation model of a suspension insulator string according to the present invention;
FIG. 3 is a schematic diagram of a conventional transmission tower calculation circuit of the present invention;
FIG. 4 is a schematic diagram of the distribution of the equivalent static wind load and its components of the conventional transmission line of the present invention;
FIG. 5 is a schematic diagram of the comparison of the numerical simulation wind field characteristics of the present invention with the theoretical wind field characteristics;
FIG. 6 is a schematic diagram of the ultra-high span transmission tower calculation circuit of the present invention;
FIG. 7 is a schematic diagram of the distribution of the equivalent static wind load and its components of the conventional transmission line of the present invention;
FIG. 8 is a schematic diagram comparing the numerical simulation wind field characteristics with the theoretical wind field characteristics of the present invention;
FIG. 9 is a schematic diagram of the distribution of the equivalent static wind load and the component thereof of the ultra-high large-span transmission line of the invention;
FIG. 10 is a schematic illustration of the effect of wind speed on wind yaw angle;
FIG. 11 is a schematic illustration of the effect of head to yaw;
FIG. 12 is a schematic illustration of the effect of pitch on windage yaw;
FIG. 13 is a schematic view of the effect of the droop ratio on the windage yaw;
FIG. 14 is a schematic view of the effect of range on windage yaw;
FIG. 15 is a graphical illustration of the effect of terrain roughness class on wind deflection angle.
Detailed Description
The following provides a more detailed description of the embodiments and the operation of the present invention with reference to the accompanying drawings.
In the embodiment of the invention, a method for calculating the maximum wind deflection angle and the wind vibration coefficient of a suspension insulator string by considering linear and linear length influence factors can be seen by combining figure 1, and the method comprises the following specific steps:
s1: determining the type of the power transmission tower and the arrangement scheme of the power transmission tower, the lead and the insulator string, and acquiring physical parameters of the lead and physical parameters of the insulator on the power transmission tower;
the physical parameters of the lead at least comprise the linearity and the length of the lead;
the transmission tower types include at least conventional transmission towers and ultra-high transmission towers;
the physical parameters of the lead at least comprise the type of the lead, the calculated sectional area of the lead, the elastic modulus of the lead, the linear density and the outer diameter of the lead;
and the physical parameters of the insulator string on the power transmission tower at least comprise the length of the insulator string, the elastic modulus of the insulator string, the quality of the insulator string and the wind shielding area of the insulator string.
S2: determining a calculation model of the wind deflection angle of the suspension insulator string by an LRC method by taking the gravity and the average wind load as initial conditions for calculating the lead and the suspension insulator string, wherein the calculation model is a schematic structural diagram of the calculation model of the wind deflection angle of the suspension insulator string as can be seen by combining FIG. 2;
as can be seen from the combination of FIG. 2, the wind deflection angle caused by the connection point A of the lead and the insulator string, the end point B of the suspended insulator string, the end point movement point B ' of the insulator string in the dynamic state, and the point B ' moving to the point B ' in the dynamic state is set
Figure RE-RE-GDA0002610855590000091
Downwind displacement of point B under action of wire span L and average wind load
Figure RE-RE-GDA0002610855590000092
A. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angle
Figure RE-RE-GDA0002610855590000093
Horizontal distance a' from coordinate origin to lowest point of wire, and lowest point of wireHorizontal distance b' to the end of the wire.
S3: combining the data of the step S1 and the data of the step S2, calculating to obtain the equivalent static wind load and the component data of the conductor in unit area between the transmission towers, and correspondingly obtaining the equivalent static wind load and the component distribution diagram of the conductor in unit area between the transmission towers;
in step S3, the equivalent static wind load p per unit area of the conductor between the transmission towersESWLThe calculation formula of (2) is as follows:
Figure RE-RE-GDA0002610855590000094
wherein (: i) represents the ith column element of the matrix;
Figure RE-RE-GDA0002610855590000095
equivalent background wind pressure;
Figure RE-RE-GDA0002610855590000096
the average wind load is obtained;
in FIG. 2, point B is selected as the target point for calculation, when the downwind displacement y of the target point is reachedBWhen the maximum wind deflection angle is reached, the insulator string wind deflection angle
Figure RE-RE-GDA0002610855590000101
Reaches the maximum value
Figure RE-RE-GDA0002610855590000102
With yBFor the target response, it can be calculated according to the above formula
Figure RE-RE-GDA0002610855590000103
Equivalent static wind load.
The matrix expression of the vibration equation of the lead under the action of wind load is as follows:
Figure RE-RE-GDA0002610855590000104
in the formula,
Figure RE-RE-GDA0002610855590000105
y' is the acceleration, the speed and the displacement of the wire node along the wind direction under the action of the pulsating wind load respectively;
Figure RE-RE-GDA0002610855590000106
is the displacement of the wire node along the wind direction under the action of average wind load
M is a quality matrix; c is a damping matrix; a K stiffness matrix; l issIs a node dependent area matrix;
the matrix expression of the vibration equation of the lead under the action of fluctuating wind load is as follows:
Figure RE-RE-GDA0002610855590000107
s4: calculating the maximum wind drift angle of the suspension insulator string according to the equivalent static wind load of the step S3;
in conjunction with FIG. 2, point B 'is shifted to point B' to cause a wind slip angle of
Figure RE-RE-GDA0002610855590000108
When in use
Figure RE-RE-GDA0002610855590000109
Reach extreme value along downwind direction
Figure RE-RE-GDA00026108555900001010
Wind deflection angle under the action of incoming wind load
Figure RE-RE-GDA00026108555900001011
Reach extreme value
Figure RE-RE-GDA00026108555900001012
The calculation formula for calculating the maximum wind drift angle of the suspension insulator string by the equivalent static wind load is as follows:
Figure RE-RE-GDA00026108555900001013
in the formula,
Figure RE-RE-GDA00026108555900001014
is the downwind peak displacement of the point B under the action of fluctuating wind load
Figure RE-RE-GDA00026108555900001015
lABA, B is the length of the insulator string between two points;
Figure RE-RE-GDA00026108555900001016
Figure RE-RE-GDA00026108555900001017
is the downwind displacement of the point B under the action of average wind load,
Figure RE-RE-GDA00026108555900001018
Figure RE-RE-GDA00026108555900001019
is the average wind deflection angle; the specific calculation formula is as follows:
Figure RE-RE-GDA00026108555900001020
Figure RE-RE-GDA00026108555900001021
Gvrespectively taking the average wind load and the vertical gravity load of the suspension insulator string at the target point;
Figure RE-RE-GDA00026108555900001022
Wvrespectively the average wind load and the vertical load transferred to the suspension insulator string by the lead at the target point.
The wire at the target point transmits the average wind load to the suspension insulator string
Figure RE-RE-GDA00026108555900001023
The calculation formula of (2) is as follows:
Figure RE-RE-GDA0002610855590000111
Figure RE-RE-GDA0002610855590000112
in the formula, NcThe number of the split conductors;
Figure RE-RE-GDA0002610855590000113
the uniform average wind load of the unit wire length of a single wire is obtained;hthe calculation mode is a pair formula for the line length of the lead in the horizontal span
Figure RE-RE-GDA0002610855590000114
Performing curve integration at a horizontal span;
wherein,
Figure RE-RE-GDA0002610855590000115
in the formula,
Figure RE-RE-GDA0002610855590000116
is the load p' and the response yBThe correlation coefficient of (a);
Figure RE-RE-GDA0002610855590000117
is a response y in the initial conditionBThe influence line of (2).
When the power transmission tower is a conventional power transmission tower, the vertical load W transmitted to the suspension insulator string by the lead at the target pointvThe calculation formula of (2) is as follows: wv＝NcPv v；
Wherein, PvThe gravity of the unit wire length of a single wire;vthe specific calculation formula is that the length of the conducting wire in the vertical span is as follows: within vertical span pair formula
Figure RE-RE-GDA0002610855590000118
Performing curve integration;
when the power transmission tower is an ultrahigh power transmission tower, the lead at the target point transmits a vertical load W to the suspension insulator stringvThe calculation formula of (2) is as follows: wv＝Pv l+Tvl+Pv r+Tvr
Wherein,l、rrespectively calculating the lengths of the left span and the right span of the target point; t isvl、TvrThe vertical components of the tension at the lowest points of the left and right two cross-wires of the target point are respectively;
when the slope of the geometric line shape of the wire at a certain point across the wire is 0:
Figure RE-RE-GDA0002610855590000119
Tvl＝0；
when the slope of the wire at the geometrical line within the span is not 0:
Figure RE-RE-GDA00026108555900001110
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw＝σo4Ac；
Wherein,
Figure RE-RE-GDA00026108555900001111
in the formula, subscripts "3" and "4" represent a no-wind state and an average wind state, respectively; a. thecThe stress area of the lead is defined; ecIs the modulus of elasticity of the wire; gamma raycIs the comprehensive specific load of the lead wires,
Figure RE-RE-GDA0002610855590000121
ywin order to obtain the average wind pressure specific load,
Figure RE-RE-GDA0002610855590000122
Figure RE-RE-GDA0002610855590000123
the calculation formula is the average wind load of the unit line length of the lead:
Figure RE-RE-GDA0002610855590000124
lrfor representing span βrTo represent the altitude difference angle.
S5: and calculating the wind vibration coefficient of the suspension insulator string according to the equivalent static wind load of the step S3.
The calculation formula for calculating the wind vibration coefficient beta of the suspension insulator string by the equivalent static wind load is as follows:
Figure RE-RE-GDA0002610855590000125
Figure RE-RE-GDA0002610855590000126
Figure RE-RE-GDA0002610855590000127
∑Crepresenting summing elements within a computational domain;ccalculating the line length of the wire in the domain;
Figure RE-RE-GDA0002610855590000128
the average wind load is obtained;
Figure RE-RE-GDA0002610855590000129
equivalent background wind pressure.
The standard value expression of the horizontal wind load of the lead/ground wire of DL/T5154 is as follows:
Figure RE-RE-GDA00026108555900001210
wherein α' is wind pressure uneven coefficient less than 1scIs coefficient of resistance；βcTaking 1 when calculating the wind deflection angle for adjusting the coefficient of the wind load; dcCalculating the outer diameter of the sub-conductor/ground wire; l ispThe horizontal span of the tower; b islThe coefficient of increase of wind load during ice coating.
α′βcThe average wind load by considering wind pressure non-uniformity is multiplied by βcThus, α' β determines the equivalent static wind load of the lead/ground wirecβ. according to the physical meaning,
Figure RE-RE-GDA00026108555900001211
β calculated using LRC is not constant and is based on p for ease of design and useESWLThe distribution characteristics of (A) are processed in an averaging mode, and the calculation is consistent β. pESWLConvex at the target point position and close to the target point position far away
Figure RE-RE-GDA0002610855590000131
Is non-uniformly distributed. Therefore, a calculation domain is set, and the equivalent static wind load of the target point is averaged in the calculation domain. And when the height difference between the target point and the adjacent tower wire hanging point is 0, selecting the target point horizontal span as the calculation domain. When the height difference exists, the equivalent static wind load at the position of the target point is more convex, so that the calculation domain is spanned by the left and right sides 1/4 of the selected target point.
Specifically, the power transmission lines of the conventional power transmission tower and the ultra-high power transmission tower are calculated and explained in an implementation mode respectively.
The first embodiment: conventional transmission towers and transmission lines.
And calculating the power transmission line of a certain 500kV strain section. The scheme for arranging the power transmission tower, the lead and the insulator string is set to be arranged according to the scheme of a strain tower-a tangent tower-a strain tower. The physical parameters of the insulator string on the strain tower and the linear tower are shown in the table 2. The total length of the strain section is 2200m, no corner is arranged in the section, and a schematic diagram of a calculation circuit is shown in FIG. 3. In the figure, s represents the center vertical arc of the span. And selecting a connection point of the insulator string and the lead at the hanging point 2 as a target point.
The lead type was set to 4 XJLLHA 1/G1A-575/40-45/7, and the physical parameters are shown in Table 1:
TABLE 1 JLHA 1/G1A-575/40-45/steel-cored aluminum strand physical parameters
Calculated cross-sectional area/mm2 Modulus of elasticity/MPa Linear density/(kg km)-1) Outer diameter/mm
621 63 000 1 917 32.4
TABLE 2 physical parameters of insulator chain
Position of Length/m Modulus of elasticity/MPa Mass/kg Area of wind/mm2
Strain tower 8.33 72 000 1 614.6 113 400
Straight line tower 6.832 72 000 1 238.08 101 800
In this example, the wire drag coefficient is 1.1 as specified by the load specification. The wind speed is designed to be 30m/s, and B-type landforms are adopted. The peak factor is valued with reference to the load specification, i.e. gs＝2.5。
The embodiment mainly researches the overall motion of the wires under the action of wind load, the stress of 4 sub-wires is equivalent to a single wire, and the vibration of the secondary span is ignored, so that the purpose of simplifying calculation is achieved. Because the wind speeds of all points of the wire in the whole span are not possible to be the same, in order to consider that the wind load borne by the wire in the whole span is identical with the wind speed of the whole span selected and designed, a wind pressure uneven coefficient is adopted. The wind pressure unevenness coefficient takes into account the unevenness of the distribution of the incoming wind in the span direction, and is an empirical coefficient related to the wind speed and the span. The non-uniformity of wind pressure is considered through the spatial correlation of pulsating wind when simulating a wind field.
1. And equivalent static wind load of the wind deflection angle of the suspension insulator string.
The incoming wind load can be decomposed into an average wind load and a fluctuating wind load, and the vibration equation matrix expression of the wire under the action of the wind load is as follows:
Figure RE-RE-GDA0002610855590000141
in the formula,
Figure RE-RE-GDA0002610855590000142
y' is the acceleration, the speed and the displacement of the wire node along the wind direction under the action of the pulsating wind load respectively;
Figure RE-RE-GDA0002610855590000143
the displacement of the lead joint along the wind direction under the action of average wind load.
The wire is a light flexible structure, and the performance under the strong wind load is as follows: 1) the structure is greatly deformed, and the geometric nonlinearity is obvious; 2) the structure stress and the displacement do not have a linear relation; 3) under the action of dynamic load, the structure is time-varying rigidity. Therefore, the equation (0.1) is a variable coefficient differential equation, and cannot be solved by the linear superposition principle. The wind vibration response of the wire caused by the incoming wind load can be decomposed into an average response and a pulse response. Taking the average wind state of the wire in fig. 2 as the initial condition for calculation, the matrix expression of the vibration equation of the wire under the action of the pulsating wind load is as follows:
Figure RE-RE-GDA0002610855590000144
the stiffness matrix of the formula (0.2) is independent of the function, is a constant coefficient differential equation, and can be solved by adopting a linear superposition principle.
Likewise, the equivalent static wind load of equation (0.2) is calculated using the LRC method. The equivalent static wind load distribution of the unit area of the node is as follows:
Figure RE-RE-GDA0002610855590000151
wherein (: i) represents the ith column element of the matrix;
Figure RE-RE-GDA0002610855590000152
equivalent background wind pressure. The formula (0.3) calculates the equivalent static wind load of the lead under different positions and different responses. In FIG. 2, point B is selected as the calculation target point, which is the target pointDownwind displacement y of punctuationBWhen the maximum wind deflection angle is reached, the insulator string wind deflection angle
Figure RE-RE-GDA0002610855590000153
Reaches the maximum value
Figure RE-RE-GDA0002610855590000154
With yBFor the target response, it can be calculated according to the above formula
Figure RE-RE-GDA0002610855590000155
Equivalent static wind load.
Therefore, in step S3, the equivalent static wind load p per unit area of the conductor between the transmission towersESWLThe calculation formula of (2) is as follows:
Figure RE-RE-GDA0002610855590000156
wherein (: i) represents the ith column element of the matrix;
Figure RE-RE-GDA0002610855590000157
equivalent background wind pressure;
Figure RE-RE-GDA0002610855590000158
the average wind load is obtained;
the matrix expression of the vibration equation of the lead under the action of wind load is as follows:
Figure RE-RE-GDA0002610855590000159
in the formula,
Figure RE-RE-GDA00026108555900001510
y' is the acceleration, the speed and the displacement of the wire node along the wind direction under the action of the pulsating wind load respectively;
Figure RE-RE-GDA00026108555900001511
is the displacement of the wire node along the wind direction under the action of average wind load
M is a quality matrix; c is a damping matrix; a K stiffness matrix; l issIs a node dependent area matrix;
the matrix expression of the vibration equation of the lead under the action of fluctuating wind load is as follows:
Figure RE-RE-GDA00026108555900001512
in the example of the present embodiment, the distribution diagram of the equivalent static wind load per unit area of the conductive line between the transmission towers and the component thereof is shown in fig. 4. In fig. 4, the component of the equivalent static wind load is mainly the average component, and the distribution of the average component is similar to the linear shape in the windless state. The background component peaks at the target point position and approaches 0 at both ends. Due to the symmetry of the structural arrangement and the fact that the target points are located on the symmetry axis, the distribution of the equivalent static wind load and the components thereof has symmetry. And calculating the wind deflection angle and the wind vibration coefficient of the suspension insulator string by adopting the equivalent static wind load determined by the figure 4.
In this embodiment:
Figure RE-RE-GDA00026108555900001513
calculated average wind load transferred to suspension insulator string by wire
Figure RE-RE-GDA00026108555900001514
Formula Wv＝NcPv vCalculating the vertical load W transferred to the suspension insulator string by the wirev＝42.082kN；
Formula (II)
Figure RE-RE-GDA0002610855590000161
Calculated mean declination
Figure RE-RE-GDA0002610855590000162
Formula (II)
Figure RE-RE-GDA0002610855590000163
The order of B points under the action of the calculated average wind loadWind direction displacement is obtained
Figure RE-RE-GDA0002610855590000164
Formula (II)
Figure RE-RE-GDA0002610855590000165
Calculated maximum wind deflection angle
Figure RE-RE-GDA0002610855590000166
Formula (II)
Figure RE-RE-GDA0002610855590000167
Calculated
Figure RE-RE-GDA0002610855590000168
Formula (II)
Figure RE-RE-GDA0002610855590000169
Calculated
Figure RE-RE-GDA00026108555900001610
Formula (II)
Figure RE-RE-GDA00026108555900001611
The calculated wind vibration coefficient β is 1.281.
And performing time domain analysis on the power transmission line, and comparing a time domain result of the wind deflection angle of the target point with a result calculated by adopting the model provided by the text so as to verify the reasonability of the calculation model. A B-class landform wind field is simulated by a harmonic synthesis method, a wind speed spectrum is a Davenport wind speed spectrum used by a load specification, a time interval delta t is 0.125s, and the simulation time is 1024 s. The average wind profile, turbulence profile and wind speed power spectrum at the target point of the simulated wind field are compared to theoretical values as shown in fig. 5. The comparison result shows that the goodness of fit of the analog value and the theoretical value is good.
And establishing an ANSYS finite element model, importing simulation wind field data, and performing time domain analysis on the wind deflection angle of the insulator string. The lead is simulated by using a Link10 unit, and the insulator string is simulated by using a Link8 unit. The wire length is 10m each, and the structural damping ratio is 0.4%, and the aerodynamic damping is considered. The initial operating tension was 55.393 kN. In time domain calculation, the average value of the wind deflection angles of the target points is 58.855 degrees, and the maximum value is 58.375 degrees. Compared with a time domain result, the relative errors of the average wind deflection angle and the maximum wind deflection angle of the calculation model are 0.14% and 0.27% respectively, and the goodness of fit of the calculation model result and the time domain result is good.
As a second embodiment: ultra-high span transmission towers and their lines.
And calculating the super-high large span, wherein the target point corresponds to the measuring point of the test windage yaw angle. The insulator strings are arranged in a duplex manner, the length of each insulator string is 19m, the elastic modulus is 200Gpa, and the wind shielding area is 28m2The mass was 6446.512 kg. The resistance coefficient of the wire is 1.1 according to the specification of the load. Calculation lines as shown in fig. 6, the difference of the wire suspension point height is 328m, and the edge-crossing wire is not 0 at the slope of the geometric line shape in the span. Considering the gradient wind height influence, the distribution of the equivalent static wind load per unit area of the line and the component thereof is determined by the formulas (0.1) to (0.3) as shown in fig. 7. In fig. 7, the component of the equivalent static wind load is dominated by the average component. When the height is lower than the gradient wind height, the average component distribution is similar to the line shape in the windless state. When the height is greater than the gradient wind height, the average component is constant. The background component peaks near the target point and approaches 0 at both ends. The background component is 0 when the height is greater than the gradient wind height. And calculating the wind deflection angle and the wind vibration coefficient of the suspension insulator string by adopting the equivalent static wind load determined by the figure 6. Wherein, the formula
Figure RE-RE-GDA0002610855590000171
Calculated
Figure RE-RE-GDA0002610855590000172
Formula (II)
Figure RE-RE-GDA0002610855590000173
Calculated Tvl460.089; formula Wv＝Pv l+Tvl+Pv r+TvrCalculated Wv1276.654 kN; formula (II)
Figure RE-RE-GDA0002610855590000174
Calculated
Figure RE-RE-GDA0002610855590000175
Formula (II)
Figure RE-RE-GDA0002610855590000176
Calculated
Figure RE-RE-GDA0002610855590000177
Formula (II)
Figure RE-RE-GDA0002610855590000178
Calculated
Figure RE-RE-GDA0002610855590000179
Formula (II)
Figure RE-RE-GDA00026108555900001710
Calculated
Figure RE-RE-GDA00026108555900001711
Formula (II)
Figure RE-RE-GDA00026108555900001712
Calculated
Figure RE-RE-GDA00026108555900001713
Formula (II)
Figure RE-RE-GDA00026108555900001714
Calculated β ═ 1.108.
Considering the influence of the gradient wind height, the maximum wind deflection angle of the target point determined by the wind tunnel test is 49.416 degrees. Compared with the test result considering the gradient wind height correction, the relative error of the maximum wind drift angle of the calculation model is 3.24%, and the goodness of fit between the calculation model result and the test result is good.
In the expression of beta, the main influence parameters are wind speed, wire height, spatial correlation and turbulence, and are divided into wind speed, altitude difference, pitch, sag ratio, span and ground roughness according to the basic parameters of design engineering. And sequentially changing the values of the parameters in the engineering use range, and performing parameter analysis on the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by adopting the calculation model provided by the section.
(1) Wind speed parameter
The design basic wind speed is changed, the wind speed working condition is 10m/s-50m/s, and the interval is 10 m/s. The wind deflection angle distribution corresponding to different wind speeds is shown in fig. 10. The average wind deflection angle and the maximum wind deflection angle are in a nonlinear increasing trend along with the increase of the wind speed, and the increasing amplitude is reduced, which shows that the tangential stiffness of the wind deflection angle of the line is gradually increased.
The influence of the wind speed on the wind vibration coefficient of the wind slip angle is shown in table 3. With the increase of the wind speed, beta shows a nonlinear increasing trend, and the increasing amplitude of beta is gradually increased under the influence of the rigidity change.
TABLE 3 influence of wind speed on the wind vibration coefficient of the wind declination
Basic wind speed/(m/s) 10 20 30 40 50
β 1.246 1.261 1.281 1.344 1.424
(2) Height difference parameter
The effect of the elevation difference was investigated by global translation of the left 1 st transconductor in fig. 3. And the horizontal moving part moves downwards when the height difference is positive, and moves upwards when the height difference is negative. The working condition of the difference between the target point and the hanging point of the wire at the position of the hanging point 1 is minus 60 m-plus 20m, and the interval is 20 m. Along with the reduction of the height difference, the windward load of the conducting wire is increased, and the vertical load transmitted to the target point suspension insulator string by the conducting wire is reduced. According to the formula
Figure RE-RE-GDA0002610855590000181
Figure RE-RE-GDA0002610855590000182
Increase of WvThe number of the grooves is reduced, and the,
Figure RE-RE-GDA0002610855590000183
will increase. The wind deflection angle distributions corresponding to different height differences are shown in fig. 11. The average wind drift angle and the maximum wind drift angle gradually increase along with the reduction of the height difference, and the height difference and the wind drift angle are approximately in a linear relation.
The effect of the height difference on the wind vibration coefficient of the wind deflection angle is shown in table 4. With the increase of the height difference, beta tends to increase first, then decrease and then increase, and beta is larger when the height difference exists than when the height difference does not exist. Within a certain range, the influence of the height difference on the beta is small.
TABLE 4 Effect of height differential on the wind vibration coefficient of the wind deflection angle
Height difference/m -60 -40 -20 0 20
β 1.242 1.258 1.305 1.281 1.286
(3) High call parameter
And changing the size of the designed calling height, and translating the whole line up and down, wherein the working condition of calling height is 25-105 m and the interval is 20 m. It can be seen from fig. 8(a) that as the breath height increases, the mean wind speed increases and the turbulence decreases. The wind deflection angle distributions corresponding to different breathing heights are shown in fig. 12. As the breath height increases, the average wind deflection angle and the maximum wind deflection angle have a nonlinear increasing trend, and the increasing amplitude is reduced.
The effect of the pitch on the wind vibration coefficient of the wind deflection angle is shown in table 5. Under the combined action of average wind and turbulence, the influence of breath height on beta is small, and the influence is not obviously regular.
TABLE 5 influence of breath height on the wind vibration coefficient of the wind declination
Call height/m 25 45 65 85 105
β 1.327 1.281 1.291 1.284 1.305
(4) Sag ratio parameter
The designed vertical span ratio is changed, and the working conditions of the vertical span ratio are 1%, 2.33%, 3%, 4% and 5%. After the vertical span ratio is increased, the wind load of the lead is reduced, the gravity rigidity of the lead is reduced due to the reduction of initial tension, and the vertical span ratio comprehensively influences the windage yaw characteristic of the suspension insulator string from two aspects. The wind deflection angle distributions corresponding to different vertical span ratios are shown in fig. 13. The average wind deflection angle and the maximum wind deflection angle are in a nonlinear decreasing trend along with the increase of the vertical span ratio.
(5) Span length parameter
Keeping the vertical span ratio of the lead unchanged, and changing the span, wherein the working condition of the span is 150-950 m and the interval is 200 m. As the span increases, the spatial dependence of the wire wind load decreases. The wind deflection angle distributions for different spans are shown in fig. 14. With the increase of the span, the average wind deflection angle and the maximum wind deflection angle tend to increase first and then decrease. At a span of 550m, the wind deflection angle reaches a maximum.
The effect of the span on the wind coefficient of windage is shown in table 5. As the range is increased, β decreases, and the magnitude of the decrease gradually decreases.
TABLE 5 influence of span on the wind vibration coefficient of the wind declination
Span/m 150 350 550 750 950
β 1.409 1.313 1.281 1.250 1.246
(6) Roughness parameter of ground
The load-normalized ground roughness categories are classified into 4 categories. Most transmission lines belong to the B-type ground roughness; for a large-span power transmission line, consideration can be given according to A-type ground roughness; a few transmission lines in the city belong to the C-type or D-type ground roughness. The designed ground roughness category is changed, and the working condition is A category-D category. The wind deflection angle distributions for different terrain roughness categories are shown in fig. 15. The ground roughness category is changed from A category to D category, and the average wind deflection angle and the maximum wind deflection angle are gradually reduced.
The effect of the terrain roughness category on the wind coefficient of wind deflection is shown in table 6. The terrain roughness category is changed from A-D, with β gradually increasing and increasing in magnitude.
TABLE 6 influence of the ground roughness class on the wind vibration coefficient of the wind declination
Class of roughness of ground A B C D
β 1.246 1.281 1.414 1.622
It should be noted that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make variations, modifications, additions or substitutions within the spirit and scope of the present invention.

Claims (8)

1. A method for calculating the maximum wind deflection angle and the wind vibration coefficient of a suspension insulator string by considering linear shape and linear length influence factors is characterized by comprising the following specific steps of:
s1: determining the type of the power transmission tower and the arrangement scheme of the power transmission tower, the lead and the insulator string, and acquiring physical parameters of the lead and physical parameters of the insulator on the power transmission tower;
the physical parameters of the lead at least comprise the linearity and the length of the lead;
s2: determining a calculation model of the wind deflection angle of the suspension insulator string by an LRC method by taking the gravity and the average wind load as initial conditions for calculating the lead and the suspension insulator string;
s3: combining the data of the step S1 and the data of the step S2, calculating to obtain the equivalent static wind load and the component data of the conductor in unit area between the transmission towers, and correspondingly obtaining the equivalent static wind load and the component distribution diagram of the conductor in unit area between the transmission towers;
s4: calculating the maximum wind drift angle of the suspension insulator string according to the equivalent static wind load of the step S3;
s5: and calculating the wind vibration coefficient of the suspension insulator string according to the equivalent static wind load of the step S3.
2. The method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 1, wherein the method comprises the following steps:
the transmission tower types include at least conventional transmission towers and ultra-high transmission towers;
the physical parameters of the lead at least comprise the type of the lead, the calculated sectional area of the lead, the elastic modulus of the lead, the linear density and the outer diameter of the lead;
and the physical parameters of the insulator string on the power transmission tower at least comprise the length of the insulator string, the elastic modulus of the insulator string, the quality of the insulator string and the wind shielding area of the insulator string.
3. The method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 1, wherein the method comprises the following steps: in the step S2, when the model of calculating the wind deflection angle of the suspension insulator string is described, a connection point a between the lead and the insulator string, a tail end point B of the suspension insulator string, a motion point B 'at the tail end point of the insulator string in a dynamic state, and a wind deflection angle caused by moving the point B' to B ″ in the dynamic state are set
Figure FDA0002433917420000021
Downwind displacement of point B under action of wire span L and average wind load
Figure FDA0002433917420000022
A. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angle
Figure FDA0002433917420000023
The horizontal distance a 'from the origin of coordinates to the lowest point of the wire, and the horizontal distance b' from the lowest point of the wire to the end of the wire.
4. The method for calculating the maximum wind drift angle and the wind vibration coefficient of the suspended insulator string by considering the linear shape and the linear length influence factors according to claim 3, wherein in the step S3, the equivalent static wind load p per unit area of the conducting wire between the power transmission towersESWLThe calculation formula of (2) is as follows:
Figure FDA0002433917420000024
wherein (: i) represents the ith column element of the matrix;
Figure FDA0002433917420000025
equivalent background wind pressure;
Figure FDA0002433917420000026
the average wind load is obtained;
the matrix expression of the vibration equation of the lead under the action of wind load is as follows:
Figure FDA0002433917420000027
in the formula,
Figure FDA0002433917420000028
y' is the acceleration, the speed and the displacement of the wire node along the wind direction under the action of the pulsating wind load respectively;
Figure FDA0002433917420000029
for the following of the wire node under the action of average wind loadDisplacement of wind direction
M is a quality matrix; c is a damping matrix; a K stiffness matrix; l issIs a node dependent area matrix;
the matrix expression of the vibration equation of the lead under the action of fluctuating wind load is as follows:
Figure FDA00024339174200000210
5. the method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 3, wherein the method comprises the following steps: the calculation formula for calculating the maximum wind drift angle of the suspension insulator string by the equivalent static wind load is as follows:
Figure FDA00024339174200000211
in the formula,
Figure FDA00024339174200000212
is the downwind peak displacement of the point B under the action of fluctuating wind load
Figure FDA00024339174200000213
lABA, B is the length of the insulator string between two points;
Figure FDA0002433917420000031
Figure FDA0002433917420000032
is the downwind displacement of the point B under the action of average wind load,
Figure FDA0002433917420000033
Figure FDA0002433917420000034
is the average wind deflection angle; the specific calculation formula is as follows:
Figure FDA0002433917420000035
Figure FDA0002433917420000036
Gvrespectively taking the average wind load and the vertical gravity load of the suspension insulator string at the target point;
Figure FDA0002433917420000037
Wvrespectively the average wind load and the vertical load transferred to the suspension insulator string by the lead at the target point.
6. The method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 5, wherein the method comprises the following steps: the wire at the target point transmits the average wind load to the suspension insulator string
Figure FDA0002433917420000038
The calculation formula of (2) is as follows:
Figure FDA0002433917420000039
Figure FDA00024339174200000310
in the formula, NcThe number of the split conductors;
Figure FDA00024339174200000311
the uniform average wind load of the unit wire length of a single wire is obtained;hthe calculation mode is a pair formula for the line length of the lead in the horizontal span
Figure FDA00024339174200000312
Integration of curves at horizontal span；
Wherein,
Figure FDA00024339174200000313
in the formula,
Figure FDA00024339174200000314
is the load p' and the response yBThe correlation coefficient of (a);
Figure FDA00024339174200000315
is a response y in the initial conditionBThe influence line of (2).
7. The method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 5, wherein the method comprises the following steps: when the power transmission tower is a conventional power transmission tower, the vertical load W transmitted to the suspension insulator string by the lead at the target pointvThe calculation formula of (2) is as follows: wv＝NcPv v；
Wherein, PvThe gravity of the unit wire length of a single wire;vthe specific calculation formula is that the length of the conducting wire in the vertical span is as follows: within vertical span pair formula
Figure FDA00024339174200000316
Performing curve integration;
when the power transmission tower is an ultrahigh power transmission tower, the lead at the target point transmits a vertical load W to the suspension insulator stringvThe calculation formula of (2) is as follows: wv＝Pv l+Tvl+Pv r+Tvr
Wherein,l、vrespectively calculating the lengths of the left span and the right span of the target point; t isvl、TvrThe vertical components of the tension at the lowest points of the left and right two cross-wires of the target point are respectively;
when the slope of the geometric line shape of the wire at a certain point across the wire is 0:
Figure FDA0002433917420000041
Tvl＝0；
when the slope of the wire at the geometrical line within the span is not 0:
Figure FDA0002433917420000042
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw＝σo4Ac；
Wherein,
Figure FDA0002433917420000043
in the formula, subscripts "3" and "4" represent a no-wind state and an average wind state, respectively; a. thecThe stress area of the lead is defined; ecIs the modulus of elasticity of the wire; gamma raycIs the comprehensive specific load of the lead wires,
Figure FDA0002433917420000044
γwin order to obtain the average wind pressure specific load,
Figure FDA0002433917420000045
Figure FDA0002433917420000046
the calculation formula is the average wind load of the unit line length of the lead:
Figure FDA0002433917420000047
lrfor representing span βrTo represent the altitude difference angle.
8. The method for calculating the maximum wind deflection angle and the wind vibration coefficient of the suspension insulator string by considering the linear shape and the linear length influence factors according to claim 4, wherein the method comprises the following steps: the calculation formula for calculating the wind vibration coefficient beta of the suspension insulator string by the equivalent static wind load is as follows:
Figure FDA0002433917420000051
Figure FDA0002433917420000052
Figure FDA0002433917420000053
∑Crepresenting summing elements within a computational domain;ccalculating the line length of the wire in the domain;
Figure FDA0002433917420000054
the average wind load is obtained;
Figure FDA0002433917420000055
equivalent background wind pressure.
CN202010245692.1A 2020-03-31 2020-03-31 Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors Active CN111651805B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010245692.1A CN111651805B (en) 2020-03-31 2020-03-31 Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010245692.1A CN111651805B (en) 2020-03-31 2020-03-31 Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors

Publications (2)

Publication Number Publication Date
CN111651805A true CN111651805A (en) 2020-09-11
CN111651805B CN111651805B (en) 2022-12-27

Family

ID=72346457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010245692.1A Active CN111651805B (en) 2020-03-31 2020-03-31 Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors

Country Status (1)

Country Link
CN (1) CN111651805B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112100724A (en) * 2020-09-16 2020-12-18 广东电网有限责任公司电力科学研究院 Wind resistance evaluation method and related device for power transmission tower
CN112484639A (en) * 2020-11-23 2021-03-12 国网北京市电力公司 Method and device for determining windage yaw position of wire, storage medium and processor
CN113109587A (en) * 2021-03-15 2021-07-13 南方电网科学研究院有限责任公司 Method, system and device for measuring wind speed of power transmission line and storage medium
CN113704937A (en) * 2021-08-06 2021-11-26 国网江苏省电力有限公司 Method for displaying postures of wires under different working conditions of distribution network
CN114595609A (en) * 2022-03-10 2022-06-07 大连理工大学 Pseudo-static force calculation method for lateral unbalanced force of power transmission line

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676667A (en) * 2013-12-04 2014-03-26 国家电网公司 Three-dimensional simulation, measurement and control system for transformer substation project flexible conductor assembling
CN104463353A (en) * 2014-11-20 2015-03-25 广东电网有限责任公司河源供电局 Intelligent prediction method for ultra-high plant safe distance of overhead transmission line
CN104655171A (en) * 2015-01-14 2015-05-27 西安工程大学 Electric transmission line simulation de-icing method based on concentrated load
CN104715105A (en) * 2015-02-11 2015-06-17 福州大学 Maximum sag modeling method of substation flexible conductor on basis of catenary
CN106597233A (en) * 2016-12-06 2017-04-26 武汉大学 Power transmission line windage yaw discharge early warning method considering dynamic air density
DE102015221553A1 (en) * 2015-11-03 2017-05-04 Innogy Se Overhead transmission system and method for extending the electrical transmission capacity of an overhead transmission system
CN109033666A (en) * 2018-08-07 2018-12-18 国网河南省电力公司电力科学研究院 Ice coating wire based on beam element waves ANSYS method for numerical simulation
CN109271751A (en) * 2018-11-16 2019-01-25 重庆科技学院 A kind of maximum dynamic angle of wind deflection of suspension insulator determines method
CN109614659A (en) * 2018-11-16 2019-04-12 重庆科技学院 A kind of suspension insulator windage yaw load inversion coefficient determines method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676667A (en) * 2013-12-04 2014-03-26 国家电网公司 Three-dimensional simulation, measurement and control system for transformer substation project flexible conductor assembling
CN104463353A (en) * 2014-11-20 2015-03-25 广东电网有限责任公司河源供电局 Intelligent prediction method for ultra-high plant safe distance of overhead transmission line
CN104655171A (en) * 2015-01-14 2015-05-27 西安工程大学 Electric transmission line simulation de-icing method based on concentrated load
CN104715105A (en) * 2015-02-11 2015-06-17 福州大学 Maximum sag modeling method of substation flexible conductor on basis of catenary
DE102015221553A1 (en) * 2015-11-03 2017-05-04 Innogy Se Overhead transmission system and method for extending the electrical transmission capacity of an overhead transmission system
CN106597233A (en) * 2016-12-06 2017-04-26 武汉大学 Power transmission line windage yaw discharge early warning method considering dynamic air density
CN109033666A (en) * 2018-08-07 2018-12-18 国网河南省电力公司电力科学研究院 Ice coating wire based on beam element waves ANSYS method for numerical simulation
CN109271751A (en) * 2018-11-16 2019-01-25 重庆科技学院 A kind of maximum dynamic angle of wind deflection of suspension insulator determines method
CN109614659A (en) * 2018-11-16 2019-04-12 重庆科技学院 A kind of suspension insulator windage yaw load inversion coefficient determines method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HOOMAN KEYHAN，等: "On the influence of wind-conductor interactions in stress analysis of overhead transmission line towers", 《ADVANCES IN CIVIL, ENVIRONMENTAL, AND MATERIALS RESEARCH (ACEM’ 12)》 *
曹树森，等: "接触线风振疲劳可靠性分析", 《机械科学与技术》 *
程唯，等: "绝缘子串非线性风偏可靠度研究*", 《建筑结构》 *
赵爽，等: "基于风洞试验的苏通大跨越输电塔风振系数研究", 《建筑结构学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112100724A (en) * 2020-09-16 2020-12-18 广东电网有限责任公司电力科学研究院 Wind resistance evaluation method and related device for power transmission tower
CN112100724B (en) * 2020-09-16 2021-10-19 广东电网有限责任公司电力科学研究院 Wind resistance evaluation method and related device for power transmission tower
CN112484639A (en) * 2020-11-23 2021-03-12 国网北京市电力公司 Method and device for determining windage yaw position of wire, storage medium and processor
CN113109587A (en) * 2021-03-15 2021-07-13 南方电网科学研究院有限责任公司 Method, system and device for measuring wind speed of power transmission line and storage medium
CN113109587B (en) * 2021-03-15 2023-02-28 南方电网科学研究院有限责任公司 Method, system and device for measuring wind speed of power transmission line and storage medium
CN113704937A (en) * 2021-08-06 2021-11-26 国网江苏省电力有限公司 Method for displaying postures of wires under different working conditions of distribution network
CN114595609A (en) * 2022-03-10 2022-06-07 大连理工大学 Pseudo-static force calculation method for lateral unbalanced force of power transmission line
CN114595609B (en) * 2022-03-10 2024-02-09 大连理工大学 Pseudo static force calculation method for lateral unbalance force of power transmission line

Also Published As

Publication number Publication date
CN111651805B (en) 2022-12-27

Similar Documents

Publication Publication Date Title
CN111651805B (en) Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors
CN109614659B (en) Method for determining wind deflection load adjustment coefficient of suspension insulator string
CN111506946B (en) Calculation method for designing wind load of conventional power transmission tower and line based on effective load method and tower-line separation method and considering tower-line coupling influence
CN107977492B (en) Monte Carlo insulator string based nonlinear windage yaw reliability calculation method
CN106096105B (en) Power transmission circuit caused by windage transient response calculation method
CN101866395B (en) Method for calculating load of hanging point of power transmission line tower during conductor galloping
CN111985018A (en) Calculation method for designing wind load of ultrahigh large-span tower and line based on inertia force method and tower line separation method and considering tower line coupling influence
CN109271751B (en) Method for determining maximum dynamic wind deflection angle of suspension insulator string
CN110287562B (en) High-speed pantograph structure parameter optimization method
CN111523172A (en) Bridge forming linear analysis method for main cable of spatial special-shaped cable surface suspension bridge
CN111651804B (en) Conventional power transmission line design wind load calculation method considering tower line coupling influence
CN112287424B (en) Calculation method for designing wind load of ultrahigh large-span tower and line based on effective load method and tower line separation method and considering tower line coupling influence
CN105719071B (en) Method for evaluating safety performance of power transmission line under action of strong wind
CN111651806B (en) Conventional power transmission tower design wind load calculation method based on inertial force method
CN111506944B (en) Tower wind load pulsation reduction coefficient calculation method considering tower line coupling effect
CN114218699A (en) Rigid contact net modeling method based on absolute node coordinates and additional constraint method
CN111985019B (en) Wind vibration coefficient simplified calculation method based on effective load method for conventional single tower with cantilever cross arm
CN112580243A (en) Power transmission line deicing jump dynamic response simulation analysis method
CN110968935A (en) Algorithm for maximum jumping height of deicing of transmission line with random height difference
CN111651808B (en) Simplified calculation method for wind vibration coefficient of ultrahigh power transmission tower considering gradient wind and concrete quality influence
CN111506945B (en) Equivalent damping coefficient calculation method for power transmission tower resonance response based on tower line coupling influence factor
CN109190931B (en) Three-dimensional space coordinate-based hard jumper iterative calculation method
CN111651807B (en) Simplified calculation method of ultrahigh single-tower wind vibration coefficient based on effective load method
CN103050909A (en) Method for installing interphase spacers and system thereof
CN117406022B (en) Travelling wave speed dynamic correction method suitable for high-voltage line length change

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant