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 PDFInfo
- 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
Links
- 239000012212 insulator Substances 0.000 title claims abstract description 133
- 238000004364 calculation method Methods 0.000 title claims abstract description 82
- 239000000725 suspension Substances 0.000 title claims abstract description 81
- 230000005540 biological transmission Effects 0.000 claims abstract description 82
- 230000003068 static effect Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 239000004020 conductor Substances 0.000 claims abstract description 18
- 238000010586 diagram Methods 0.000 claims abstract description 15
- 230000005484 gravity Effects 0.000 claims abstract description 12
- 230000009471 action Effects 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000006073 displacement reaction Methods 0.000 claims description 24
- 230000004044 response Effects 0.000 claims description 22
- 230000014509 gene expression Effects 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 4
- 230000033001 locomotion Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 abstract description 14
- 230000000694 effects Effects 0.000 description 14
- 230000007423 decrease Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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
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 setDownwind displacement of point B under action of wire span L and average wind loadA. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angleThe 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:
Wherein (: i) represents the ith column element of the matrix;equivalent background wind pressure;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:
in the formula,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;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:
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:
in the formula,is the downwind peak displacement of the point B under the action of fluctuating wind loadlABA, B is the length of the insulator string between two points;
Gvrespectively taking the average wind load and the vertical gravity load of the suspension insulator string at the target point;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 stringThe calculation formula of (2) is as follows:
in the formula, NcThe number of the split conductors;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 spanPerforming curve integration at a horizontal span; wherein,in the formula,is the load p' and the response yBThe correlation coefficient of (a);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 formulaPerforming 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:Tvl=0;
when the slope of the wire at the geometrical line within the span is not 0:
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw=σo4Ac;
Wherein,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,γwin order to obtain the average wind pressure specific load,
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:
∑Crepresenting summing elements within a computational domain;ccalculating the line length of the wire in the domain;the average wind load is obtained;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 setDownwind displacement of point B under action of wire span L and average wind loadA. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angleHorizontal 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:
wherein (: i) represents the ith column element of the matrix;equivalent background wind pressure;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 angleReaches the maximum valueWith yBFor the target response, it can be calculated according to the above formulaEquivalent static wind load.
The matrix expression of the vibration equation of the lead under the action of wind load is as follows:
in the formula,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;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:
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 ofWhen in useReach extreme value along downwind directionWind deflection angle under the action of incoming wind loadReach extreme value
The calculation formula for calculating the maximum wind drift angle of the suspension insulator string by the equivalent static wind load is as follows:
in the formula,is the downwind peak displacement of the point B under the action of fluctuating wind loadlABA, B is the length of the insulator string between two points;
Gvrespectively taking the average wind load and the vertical gravity load of the suspension insulator string at the target point;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 stringThe calculation formula of (2) is as follows:
in the formula, NcThe number of the split conductors;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 spanPerforming curve integration at a horizontal span;
in the formula,is the load p' and the response yBThe correlation coefficient of (a);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 formulaPerforming 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;
Tvl=0;
when the slope of the wire at the geometrical line within the span is not 0:
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw=σo4Ac;
Wherein,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,ywin order to obtain the average wind pressure specific load,
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:
∑Crepresenting summing elements within a computational domain;ccalculating the line length of the wire in the domain;the average wind load is obtained;equivalent background wind pressure.
The standard value expression of the horizontal wind load of the lead/ground wire of DL/T5154 is as follows:
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,β 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 awayIs 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:
in the formula,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;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:
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:
wherein (: i) represents the ith column element of the matrix;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 angleReaches the maximum valueWith yBFor the target response, it can be calculated according to the above formulaEquivalent 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:
wherein (: i) represents the ith column element of the matrix;equivalent background wind pressure;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:
in the formula,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;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:
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.
Formula Wv=NcPv vCalculating the vertical load W transferred to the suspension insulator string by the wirev=42.082kN;
Formula (II)The order of B points under the action of the calculated average wind loadWind direction displacement is obtained
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 formulaCalculatedFormula (II)Calculated Tvl460.089; formula Wv=Pv l+Tvl+Pv r+TvrCalculated Wv1276.654 kN; formula (II)CalculatedFormula (II)CalculatedFormula (II)CalculatedFormula (II)CalculatedFormula (II)CalculatedFormula (II)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 Increase of WvThe number of the grooves is reduced, and the,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/ |
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/ |
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 setDownwind displacement of point B under action of wire span L and average wind loadA. Length l of insulator chain between two points BABThe hanging point difference h between two ends of the wire and the average wind deflection angleThe 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:
wherein (: i) represents the ith column element of the matrix;equivalent background wind pressure;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:
in the formula,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;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;
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:
in the formula,is the downwind peak displacement of the point B under the action of fluctuating wind loadlABA, B is the length of the insulator string between two points;
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 stringThe calculation formula of (2) is as follows:
in the formula, NcThe number of the split conductors;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 spanIntegration of curves at horizontal span;
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 formulaPerforming 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;
Tvl=0;
when the slope of the wire at the geometrical line within the span is not 0:
in the formula, TwThe calculation formula is the horizontal tension of a single wire in an average wind state: t isw=σo4Ac;
Wherein,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,γwin order to obtain the average wind pressure specific load,
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:
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)
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)
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 |
-
2020
- 2020-03-31 CN CN202010245692.1A patent/CN111651805B/en active Active
Patent Citations (9)
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)
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)
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 |