CN107977492A - Based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain - Google Patents
Based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain Download PDFInfo
- Publication number
- CN107977492A CN107977492A CN201711125066.3A CN201711125066A CN107977492A CN 107977492 A CN107977492 A CN 107977492A CN 201711125066 A CN201711125066 A CN 201711125066A CN 107977492 A CN107977492 A CN 107977492A
- Authority
- CN
- China
- Prior art keywords
- insulator chain
- windage yaw
- stochastic variable
- conducting wire
- wind
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Insulators (AREA)
Abstract
The invention discloses one kind to be based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain, follows the steps below：Transmission line of electricity research object is set, the windage yaw angle value of insulator chain is calculated using the straight rod method of rigidity；According to overhead transmission line specification, shaft tower and insulator chain connection structure are established, and selectes the stochastic variable that the connection structure is related to, distribution function corresponding with stochastic variable is determined according to stochastic variable；Establish geometrical relationship corresponding with bindiny mechanism；Establish the insulator chain windage yaw invalidation functions function with stochastic variable；Insulator chain windage yaw reliability is calculated using Meng Takaluo methods.Beneficial effect：Good reliability, calculating process is simple, can meet the requirement of the prior art.
Description
Technical field
The present invention relates to transmission line of electricity field, specifically a kind of insulator chain based on monte carlo method are non-linear
Windage yaw reliability degree calculation method.
Background technology
Under wind action, the insulator chain of overhead transmission line and its transmission pressure of suspension will produce different cycles
Windage yaw wave.In rocking process, if the distance between powered body portion and shaft tower are less than the electric clearance allowed,
Electric discharge phenomena will occur between power transmission line and shaft tower, that is, windage yaw discharge accident occurs.Insulator chain windage yaw discharge can serious threat
The normal operation of network system, and cause huge economic loss and social influence.
During the angle of wind deflection of domestic aerial high-voltage power transmission line calculates, insulator chain is usually reduced to rigid rod or string is polygon
Shape, the angle of wind deflection of insulator chain is calculated in the case where designing mean wind speed using static method.Engineering also all continue to use in practice this two
Kind method, and the reason for windage yaw accident is analyzed accordingly.In this regard, has the design discipline apoplexy of expert and scholar to China
The calculating of drift angle proposes without considering wind load dynamic effect doubts.Therefore, people begin to use Finite Element Method Simulation have studied
The dynamic windage yaw response of insulator chain under Random Fluctuating Wind Speed effect.Hereafter, many scholars are utilizing finite element soft both at home and abroad
Part establishes the Multi-body model of insulator chain-conducting wire coupling, and is simulated by harmony superposition and consider fluctuating wind space correlation
Property and random pulse wind field of the wind speed along height change, ultimate analysis gone out the change of the parameters such as different spans, height difference, design wind speed
Windage yaw response in the case of change.But finite element model establish it is complicated, for adjacent shaft tower and calculate between shaft tower height difference compared with
Small situation, model is established according to finite element, and process is cumbersome, and for transmission line of electricity, the structure per adjacent shaft tower is not
It is identical, lose time if establishing model using finite element, staff's work efficiency is low, therefore is necessary to propose a kind of simple
Method, to calculate when height difference is smaller between adjacent shaft tower and calculating shaft tower, the non-linear windage yaw reliability of insulator chain.
The content of the invention
In view of the above-mentioned problems, the present invention provides a kind of non-linear windage yaw of the insulator chain based on monte carlo method is reliable
Computational methods are spent, this method considers wind load dynamic effect, and windage yaw angle value is calculated using rigidity straight rod method, and then really
Determine insulator chain windage yaw reliability.
To reach above-mentioned purpose, the concrete technical scheme that the present invention uses is as follows：
One kind is based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain, its key is according to following step
It is rapid to carry out：
S1：Set transmission line of electricity research object, wind load, gravity laod, the gravity of conducting wire according to suffered by insulator chain
Wind load suffered by load, conducting wire, the windage yaw angle value of insulator chain is calculated using the straight rod method of rigidity；
S2：According to overhead transmission line specification, shaft tower and insulator chain connection structure are established, and selectes the connection structure and relates to
And stochastic variable, distribution function corresponding with stochastic variable is determined according to stochastic variable；
S3：Connection structure and stochastic variable are built according to step S2, establish geometrical relationship corresponding with bindiny mechanism；
S4：According to the geometrical relationship of step S3 and the windage yaw angle value of step S1, the insulator chain with stochastic variable is established
Windage yaw invalidation functions function；
S5：According to stochastic variable, the distribution function of stochastic variable and insulator chain windage yaw invalidation functions function, illiteracy is used
Taka Lip river method calculates insulator chain windage yaw reliability.
By above-mentioned design, when height difference is smaller between adjacent shaft tower and calculating shaft tower, using the straight rod method of rigidity to windage yaw
Angle value is calculated, the wind suffered by wind load, gravity laod, the gravity laod of conducting wire, conducting wire suffered by consideration insulator chain
The influence factors such as load, and Meng Takaluo methods are combined, calculate insulator chain windage yaw reliability, good reliability, calculating process letter
It is single, it can meet the requirement of the prior art.
Further, the content of step S1 is：
Wind load suffered by the gravity laod of wind load, gravity laod, conducting wire according to suffered by insulator chain, conducting wire, meter
Calculate insulator chain windage yaw angle value formula be：
In formula (1)：
GhIt is wind load suffered by insulator chain：Gh=W0μzA1；
GvIt is the gravity laod of insulator chain：Gv=p2l；
WvIt is the gravity laod of conducting wire：
WhIt is the wind load suffered by conducting wire：Wh=α W0μzμscβcdLp sin2θ；
Wherein, α is wind evil attacking lung；W0For fundamental wind pressure standard value；μzFor height variation coefficient of wind pressure；μscTo lead
Line Shape Coefficient；βcFor Wind Load Adjustment Coefficients；D is the outside diameter of conducting wire；LpIt is the conducting wire span representated by load node；p1、p2
Gravity laod respectively on conducting wire and insulator chain unit length；l1、 l2The respectively left and right span of shaft tower both sides；h1、h2、
The respectively left and right height difference of the relatively middle hanging point of shaft tower both sides conducting wire hanging point；T is wire tension；A1 bears wind pressure for insulator chain
Area；L is insulator chain length.
Using the above scheme, by using the straight rod method of rigidity, angle of wind deflection is calculated, it is contemplated that wind load dynamic effect.In calculating
Wind evil attacking lung α is taken to be equal to 1.0.The value of other all parameters is fully according to GB50545-2010《110kV-750kV
Overhead transmission line design specification》Regulation value.
Further describe, the stochastic variable variable includes：Wind speed, conducting wire dead weight, insulator chain length, span, absolutely
Height difference, wire diameter, insulator between the dead weight of edge substring, the air gap radius of circle, wire tension, adjacent shaft tower and calculating shaft tower
String outside diameter.
Wherein, the distribution function of stochastic variable and parameter are shown in Table 1.Since the statistical distribution otherness of wind speed is larger, herein
Shared from Chinese meteorological data and the yearly maximum wind speed data of computing electric power line location in the past 50 years has been downloaded in website, and
According to GB50009-2012《Loading code for design of building structures》And GB50545-2010《110kV-750kV overhead transmission lines design
Specification》Regulation calculate the statistical parameter of wind speed.Due to the height difference very little of computing electric power line, wire tension is to angle of wind deflection
Very little is influenced, it is then supposed that tension force and height difference are definite value.
Table 1 designs the statistical parameter of stochastic variable
| Design stochastic variable | Average | The coefficient of variation | Distribution pattern |
| Wind speed v (m/s2) | 23.635 | 0.205 | Extreme I type |
| Conducting wire dead weight p1(N/m) | 10.805 | 0.07 | Normal state |
| Insulator chain length l (m) | 6.832 | 0.05 | Normal state |
| Span LP(m) | 550 | 0.05 | Normal state |
| Insulator chain dead weight p2(N/m) | 1776.2 | 0.07 | Normal state |
| The air gap radius of circle r (m) | 3.7 | - | Definite value |
| Wire tension T (N) | 55393 | - | Definite value |
| Height difference (m) | 3.353 | - | Definite value |
| Wire diameter (mm) | 32.4 | 0.05 | Normal state |
| Insulator chain outside diameter (mm) | 360 | 0.05 | Normal state |
Further describe, step S5 calculates the particular content of insulator chain windage yaw reliability using Meng Takaluo methods
For：
S51：Selected stochastic variable is numbered, 1,2,3 ... I；And determine the distribution function of each stochastic variable
S52：N number of equally distributed random number r in section (0~1) is generated in MATLAB using multiplicative congruential methodi j；
S53：By random number ri jBring into distribution function, obtain random sampling value
S54：The random sampling value that step S53 is obtainedSubstitute into insulator chain windage yaw invalidation functions function Z, and unite
Meter, as Z ＜ 0, the number n of insulator chain windage yaw failure；
S55：According to formulaObtain insulator chain windage yaw failure probability pf；
S56：Pass through formula β=- Φ-1(pf) insulator chain windage yaw RELIABILITY INDEX is calculated, wherein：Φ-1() represents mark
The inverse function of quasi normal distribution, pfFailure probability is represented, β represents RELIABILITY INDEX.
Beneficial effects of the present invention：Suitable for the less situation of height difference between adjacent shaft tower and calculating shaft tower, using rigidity
Straight rod method calculates windage yaw angle value, the gravity laod of wind load, gravity laod, conducting wire suffered by consideration insulator chain,
The influence factors such as the wind load suffered by conducting wire, and Meng Takaluo methods are combined, calculate insulator chain windage yaw reliability, reliability
Good, calculating process is simple, can meet the requirement of the prior art.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention；
Fig. 2 is the shaft tower and suspension insulator schematic diagram of the present invention；
Fig. 3 is research object simulation schematic diagram；
Fig. 4 is the Wind Velocity History schematic diagram at No. 1 hanging point in Fig. 3；
Fig. 5 is simulation point turbulence intensity and B class landforms turbulence intensity contrast schematic diagrams；
Fig. 6 is power spectrum and target power spectrum contrast schematic diagram at characteristic point；
Fig. 7 is the graph of relation of wind speed average and windage yaw RELIABILITY INDEX；
Fig. 8 is graph of relation of the conducting wire from weight-average value and windage yaw RELIABILITY INDEX；
Fig. 9 is the graph of relation of wire diameter average and windage yaw RELIABILITY INDEX；
Embodiment
The embodiment and operation principle of the present invention are described in further detail below in conjunction with the accompanying drawings.
As seen in Figure 1, a kind of to be based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain, it is special
Sign is to follow the steps below：
S1：Set transmission line of electricity research object, wind load, gravity laod, the gravity of conducting wire according to suffered by insulator chain
Wind load suffered by load, conducting wire, the windage yaw angle value of insulator chain is calculated using the straight rod method of rigidity；
Wherein, the windage yaw angle value formula for calculating insulator chain is：
In above-mentioned formula：
GhIt is wind load suffered by insulator chain：Gh=W0μzA1；
GvIt is the gravity laod of insulator chain：Gv=p2l；
WvIt is the gravity laod of conducting wire：
WhIt is the wind load suffered by conducting wire：Wh=α W0μzμscβcdLp sin2θ；
Wherein, α is wind evil attacking lung；W0For fundamental wind pressure standard value；μzFor height variation coefficient of wind pressure；μscTo lead
Line Shape Coefficient；βcFor Wind Load Adjustment Coefficients；D is the outside diameter of conducting wire；LpIt is the conducting wire span representated by load node；p1、p2
Gravity laod respectively on conducting wire and insulator chain unit length；l1、 l2The respectively left and right span of shaft tower both sides；h1、h2、
The respectively left and right height difference of the relatively middle hanging point of shaft tower both sides conducting wire hanging point；T is wire tension；A1Wind pressure is born for insulator chain
Area；L is insulator chain length.
S2：According to overhead transmission line specification, shaft tower and insulator chain connection structure are established, and selectes the connection structure and relates to
And stochastic variable, distribution function corresponding with stochastic variable is determined according to stochastic variable；
The stochastic variable variable includes：Wind speed, conducting wire dead weight, the dead weight of insulator chain length, span, insulator chain, sky
Height difference, wire diameter, insulator chain outside diameter between gas gap radius of circle, wire tension, adjacent shaft tower and calculating shaft tower.
S3：Connection structure and stochastic variable are built according to step S2, establish geometrical relationship corresponding with bindiny mechanism；At this
In embodiment, illustrated by taking pole tower forms as an example, be specifically shown in Fig. 2.
S4：According to the geometrical relationship of step S3 and the windage yaw angle value of step S1, the insulator chain with stochastic variable is established
Windage yaw invalidation functions function；
B represents cross-arm length；L represents insulator chain length；R represents the air gap radius of circle；α is represented between cross-arm and tower body
Inclined angle, shaft tower α herein are 115 °；Represent windage yaw angle.
Wherein, the statistical parameter value of the stochastic variable in insulator chain windage yaw invalidation functions function is shown in Table 1.Due to we
The height difference very little of computing electric power line in case, influence very little of the wire tension to angle of wind deflection, it is then supposed that tension force and height difference are
Definite value.
S5：According to stochastic variable, the distribution function of stochastic variable and insulator chain windage yaw invalidation functions function, illiteracy is used
Taka Lip river method calculates insulator chain windage yaw reliability.
S51：Selected stochastic variable is numbered, 1,2,3 ... I；And determine the distribution function of each stochastic variable
S52：N number of equally distributed random number r in section (0~1) is generated in MATLAB using multiplicative congruential methodi j；
S53：By random number ri jBring into distribution function, obtain random sampling value
S54：The random sampling value that step S53 is obtainedSubstitute into insulator chain windage yaw invalidation functions function Z, and unite
Meter, as Z ＜ 0, the number n of insulator chain windage yaw failure；
S55：According to formulaObtain insulator chain windage yaw failure probability pf；
S56：Pass through formula β=- Φ-1(pf) insulator chain windage yaw RELIABILITY INDEX is calculated, wherein：Φ-1() represents mark
The inverse function of quasi normal distribution, pfFailure probability is represented, β represents RELIABILITY INDEX.
Through the above scheme, the reasonability of angle of wind deflection is calculated using the straight rod method of rigidity in order to explain.Have below by introducing
Meta-model is limited, is compared and analyzed：
In the present embodiment, research object is：Using the 3 cross-line roads of certain 500kV as research object, share 2 anchor supports and
2 tangent towers, every section of span length 550m, total length 1650m are whole across interior no corner.It is specifically shown in Fig. 3.Research object steel-cored aluminium strand thing
Reason parameter is shown in Table 2.The physical parameter of anchor support and the insulator chain on tangent tower is shown in Table 3.
2 steel-cored aluminium strand physical parameter of table
| Calculate sectional area/mm2 | Coefficient of elasticity/MPa | Line density/(kgkm-1) | Initial tension/kN |
| 621 | 630 00 | 191 7 | 55.4 |
3 insulator chain physical parameter of table
| Position | Length/m | Coefficient of elasticity/MPa | Quality/kg |
| Anchor support | 8.330 | 72000 | 1614 |
| Tangent tower | 6.832 | 72000 | 1238 |
The simulation of fluctuating wind：
S is composed using Davenportu(f)：
In formula：X=1200 ω/U10, ω is circular frequency, and n is frequency, and k is and the relevant constant of landforms, it is assumed that calculating it is defeated
Electric line is located in B class landforms, takes k=0.005.U10 is the wind speed that 10m highly locates.
Since transmission line of electricity is larger along conducting wire direction scale, fluctuating wind between conducting wire diverse location is considered herein
The spatial coherence of speed, by the cross-spectral density function of fluctuating wind speed come i in representation space, j point-to-point transmission fluctuating wind speeds
Degree of correlation, concrete form are as follows：
In formulaFor phase angle, Coh (ω) is the spatial correlation function that Davenport recommends.
Fluctuating wind field is simulated using harmony superposition：It is 1 024s to simulate time-histories total length, and time step takes 0.25s, frequency
The interception scope of rate takes 0.00rad/s-6.28rad/s, and sample frequency points are 4 096.Conducting wire upper edge line direction every
10m selects a wind speed simulation point, and symbiosis is into 166 random wind speed points.When mean wind speed at 10m height is 30m/s,
In the research object for simulating generation, wind speed time history curve at No. 1 hanging point, specific wind speed time history curve figure is shown in Fig. 4.
To examine the accuracy of wind speed simulation, by the turbulence intensity of modeling wind field characteristic point and power spectral density and target wind
Speed spectrum is compared as shown in Figure 5, Figure 6.It can be seen that from Fig. 5, Fig. 6：Wind field turbulence intensity and the analogue value of wind speed and theory
It is higher to be worth the goodness of fit.
After wind speed is simulated, wind load time-history is calculated using following formula, is applied on numerical model.
In formula：ρ is atmospheric density, and standard value is 1.225 5kg/m3；D is the equivalent diameter windward of conducting wire；L grows for node
Degree；CDFor resistance coefficient；UzAnd VzMean wind speed and fluctuating wind speed respectively at height z.The damping ratio of conducting wire takes 0.4%, and
Influence of the air damping to conducting wire.
Based on the studies above object, the straight rod method result of rigidity is contrasted with finite element result.
In table 4：Angle of wind deflection time-histories statistical value of No. 1 hanging point under various criterion wind speed and calculate according to the straight rod method of rigidity
Angle of wind deflection numerical value.Wherein：Represent the angle of wind deflection numerical value that the straight rod method of rigidity calculates；Represent the angle of wind deflection that ANSYS is calculated
The average of time-history curves；Represent the maximum for the angle of wind deflection time-history curves that ANSYS is calculated；δ represents the straight rod method meter of rigidity
Relative error between the angle of wind deflection time-history curves maximum that the angle of wind deflection numerical value and ANSYS of calculation are calculated.
In table 4, under different standard wind speeds, air damping effect, multi gear are considered in finite element model calculating
The anchor support at circuit model both ends is the work that intercouples " pullling " for producing effect, each shelves circuit are swung outside power transmission line road surface
Under precondition, the angle of wind deflection result that the straight rod method of rigidity calculates is also in engineering acceptable error range.Then
Angle of wind deflection that the straight rod method of rigidity calculates the result is that relatively safety, it is rational.
Insulator chain windage yaw reliability influence factors are analyzed：
According to GB 50068-2001《Building structure RELIABILITY DESIGN is sought unity of standard》To serviceability limit state reliability
Index is taken 0~1.5 regulation by degree of reversibility, is chosen insulator chain windage yaw target reliability index herein and is checked for 1.5.
In order to check in the case of different design stochastic variable statistical parameters, according to Current specifications design shaft tower, its
Whether windage yaw RELIABILITY INDEX reaches requirement.The cross-arm length B for assuming first that shaft tower is the unknown parameter of structure, in known table 4
In other specification in the case of.Follow the steps below calculating：
1) the design stochastic variable in table 5 determines the standard value of each parameter.Wherein the design load of wind speed according to
GB50545-2010《110kV-750kV overhead transmission line design specifications》Regulation：The transmission line of electricity circuit of 500kV takes the return period
For the design wind speed value of 50 years.Other design stochastic variable standard values are according to GB50009-2012《Loading code for design of building structures》
Regulation carry out value.
2) combined using serviceability limit state, the standard value of each design stochastic variable is substituted into formulaIn, to calculate cross-arm length B.
3) after cross-arm length B is determined, its corresponding RELIABILITY INDEX is calculated with MC methods.
4) change the statistical parameter value of some design stochastic variable, repeat 1)~3) step, wind speed average can be obtained
With the relation curve of windage yaw RELIABILITY INDEX, Fig. 7 is referred to；Conducting wire is from the relation curve of weight-average value and windage yaw RELIABILITY INDEX, in detail
See Fig. 8；The relation curve of wire diameter average and windage yaw RELIABILITY INDEX, refers to Fig. 9.
As can see from Figure 7 with the increase of wind speed average, according to the windage yaw reliability of the shaft tower of current specifications design
Index can be reduced constantly.This means in the higher area of wind speed average, shaft tower is designed according to current method for normalizing
When, it is possible that its insulator chain windage yaw RELIABILITY INDEX is less than the situation of target reliability index.And in wind speed average
Less area, the insulator chain windage yaw RELIABILITY INDEX being designed according to current specifications method to shaft tower is completely full
The requirement of foot-eye RELIABILITY INDEX.
, can according to the windage yaw of the shaft tower of current specifications design as can see from Figure 7 as conducting wire is from the reduction of weight-average value
Can constantly it be reduced by degree index.
As can see from Figure 8 with the increase of wire diameter average, the shaft tower designed according to current specifications, its insulator
String windage yaw RELIABILITY INDEX can be reduced constantly.This means for outside diameter is larger or the less conducting wire of dead weight, according to existing rule
When model is designed shaft tower, it is possible that its insulator chain windage yaw RELIABILITY INDEX is less than target reliability index
Situation.
(1) smaller in height difference, in the case that wind evil attacking lung takes 1.0, the straight rod method of rigidity calculates angle of wind deflection in specification
Resultant error be within the acceptable range, it is believed that the straight rod method of rigidity in specification is rational.
(2) it is non-linear due to insulator chain windage yaw invalidation functions function in the larger area of wind speed average, according to existing
The shaft tower of method for normalizing design, its insulator chain windage yaw RELIABILITY INDEX may be less than target reliability index.It is therefore proposed that
In the area that wind speed average is larger, more than 50 year return period, higher design wind speed value is taken to refer to improve its windage yaw reliability
Mark.
(3) it is reliable in the less area of wind speed average, the shaft tower designed according to current specifications method, its insulator chain windage yaw
Degree index fully achieves the requirement of target reliability index, it is proposed that windage yaw does not have to choose a minimum design wind in calculating
Speed value, but directly choose local air speed value.
(4) it can improve its insulator chain windage yaw reliability using larger, the less conducting wire of outside diameter of conducting oneself with dignity in the design phase and refer to
Mark.
It should be pointed out that it is limitation of the present invention that described above, which is not, the present invention is also not limited to the example above,
What those skilled in the art were made in the essential scope of the present invention changes, is modified, adds or replaces, and also should
Belong to protection scope of the present invention.
Claims (4)
1. one kind is based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain, it is characterised in that according to following steps
Carry out：
S1：Setting transmission line of electricity research object, the gravity laod of wind load, gravity laod, conducting wire according to suffered by insulator chain,
Wind load suffered by conducting wire, the windage yaw angle value of insulator chain is calculated using the straight rod method of rigidity；
S2：According to overhead transmission line specification, shaft tower and insulator chain connection structure are established, and selectes what the connection structure was related to
Stochastic variable, distribution function corresponding with stochastic variable is determined according to stochastic variable；
S3：Connection structure and stochastic variable are built according to step S2, establish geometrical relationship corresponding with bindiny mechanism；
S4：According to the geometrical relationship of step S3 and the windage yaw angle value of step S1, the insulator chain windage yaw with stochastic variable is established
Invalidation functions function；
S5：According to stochastic variable, the distribution function of stochastic variable and insulator chain windage yaw invalidation functions function, illiteracy Taka is used
Lip river method calculates insulator chain windage yaw reliability.
2. the non-linear windage yaw reliability degree calculation method of the insulator chain according to claim 1 based on Monte Carlo, it is special
Sign is that the content of step S1 is：
Wind load suffered by the gravity laod of wind load, gravity laod, conducting wire according to suffered by insulator chain, conducting wire, calculates exhausted
The windage yaw angle value formula of edge substring is：
In formula (1)：
GhIt is wind load suffered by insulator chain：Gh=W0μzA1；
GvIt is the gravity laod of insulator chain：Gv=p2l；
WvIt is the gravity laod of conducting wire：
WhIt is the wind load suffered by conducting wire：Wh=α W0μzμscβcdLpsin2θ；
Wherein, α is wind evil attacking lung；W0For fundamental wind pressure standard value；μzFor height variation coefficient of wind pressure；μscFor wire body
Type coefficient；βcFor Wind Load Adjustment Coefficients；D is the outside diameter of conducting wire；LpIt is the conducting wire span representated by load node；p1、p2Respectively
For the gravity laod on conducting wire and insulator chain unit length；l1、l2The respectively left and right span of shaft tower both sides；h1、h2, respectively
For the left and right height difference of the relatively middle hanging point of shaft tower both sides conducting wire hanging point；T is wire tension；A1The face of wind pressure is born for insulator chain
Product；L is insulator chain length.
3. the non-linear windage yaw reliability degree calculation method of the insulator chain according to claim 1 based on Monte Carlo, it is special
Sign is that the stochastic variable variable includes：Wind speed, conducting wire dead weight, the dead weight of insulator chain length, span, insulator chain, air
Height difference, wire diameter, insulator chain outside diameter between gap radius of circle, wire tension, adjacent shaft tower and calculating shaft tower.
4. the non-linear windage yaw reliability degree calculation method of the insulator chain according to claim 1 based on Monte Carlo, it is special
Sign is that the particular content of step S5 is：
S51：Selected stochastic variable is numbered, 1,2,3 ... I；And determine the distribution function of each stochastic variable
S52：N number of equally distributed random number in section (0~1) is generated in MATLAB using multiplicative congruential method
S53：By random numberBring into distribution function, obtain random sampling value
S54：The random sampling value that step S53 is obtainedSubstitute into insulator chain windage yaw invalidation functions function Z, and count, work as Z
During ＜ 0, the number n of insulator chain windage yaw failure；
S55：According to formulaObtain insulator chain windage yaw failure probability pf；
S56：Pass through formula β=- Φ-1(pf) insulator chain windage yaw RELIABILITY INDEX is calculated, wherein：Φ-1() is representing standard just
The inverse function of state distribution, pfFailure probability is represented, β represents RELIABILITY INDEX.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711125066.3A CN107977492B (en) | 2017-11-14 | 2017-11-14 | Monte Carlo insulator string based nonlinear windage yaw reliability calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711125066.3A CN107977492B (en) | 2017-11-14 | 2017-11-14 | Monte Carlo insulator string based nonlinear windage yaw reliability calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107977492A true CN107977492A (en) | 2018-05-01 |
| CN107977492B CN107977492B (en) | 2020-12-04 |
Family
ID=62013549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711125066.3A Active CN107977492B (en) | 2017-11-14 | 2017-11-14 | Monte Carlo insulator string based nonlinear windage yaw reliability calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107977492B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108876194A (en) * | 2018-07-16 | 2018-11-23 | 国网福建省电力有限公司 | Power distribution network methods of risk assessment under a kind of typhoon disaster scene |
| CN109271751A (en) * | 2018-11-16 | 2019-01-25 | 重庆科技学院 | A kind of maximum dynamic angle of wind deflection of suspension insulator determines method |
| CN109471997A (en) * | 2018-10-19 | 2019-03-15 | 海南电网有限责任公司电力科学研究院 | A kind of insulator chain windage yaw distance calculating method |
| CN109902351A (en) * | 2019-01-28 | 2019-06-18 | 浙江大学 | A kind of simplified calculation method of ice coating wire dynamic windage yaw |
| CN110378050A (en) * | 2019-07-25 | 2019-10-25 | 重庆科技学院 | A kind of power transmission tower down wind air damping compares calculation method |
| CN111157017A (en) * | 2019-11-29 | 2020-05-15 | 湖南大学 | Efficient reliability analysis method for bias value of accelerometer |
| CN114757077A (en) * | 2022-04-21 | 2022-07-15 | 内蒙古电力（集团）有限责任公司内蒙古电力经济技术研究院分公司 | Construction method of wind deflection angle prediction model of double-split line suspension insulator string |
| WO2024198505A1 (en) * | 2023-03-30 | 2024-10-03 | 贵州电网有限责任公司 | Wire wind-induced swaying load monitoring method, apparatus and device, and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105046581A (en) * | 2015-07-09 | 2015-11-11 | 广西电网有限责任公司电力科学研究院 | Multi-parameter wind load-based power transmission line strong wind tripping risk assessment method |
| CN105335569A (en) * | 2015-11-16 | 2016-02-17 | 国家电网公司 | Ultrahigh voltage transmission line multi-connected V-shaped insulator string mechanical property analogue simulation method |
| CN106096110A (en) * | 2016-06-03 | 2016-11-09 | 国网电力科学研究院武汉南瑞有限责任公司 | Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast |
| CN106096105A (en) * | 2016-06-02 | 2016-11-09 | 浙江大学 | Power transmission circuit caused by windage transient response computational methods |
| CN107092766A (en) * | 2017-06-27 | 2017-08-25 | 国网江苏省电力公司电力科学研究院 | A kind of wire windage yaw insulation check method based on finite element |
-
2017
- 2017-11-14 CN CN201711125066.3A patent/CN107977492B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105046581A (en) * | 2015-07-09 | 2015-11-11 | 广西电网有限责任公司电力科学研究院 | Multi-parameter wind load-based power transmission line strong wind tripping risk assessment method |
| CN105335569A (en) * | 2015-11-16 | 2016-02-17 | 国家电网公司 | Ultrahigh voltage transmission line multi-connected V-shaped insulator string mechanical property analogue simulation method |
| CN106096105A (en) * | 2016-06-02 | 2016-11-09 | 浙江大学 | Power transmission circuit caused by windage transient response computational methods |
| CN106096110A (en) * | 2016-06-03 | 2016-11-09 | 国网电力科学研究院武汉南瑞有限责任公司 | Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast |
| CN107092766A (en) * | 2017-06-27 | 2017-08-25 | 国网江苏省电力公司电力科学研究院 | A kind of wire windage yaw insulation check method based on finite element |
Non-Patent Citations (4)
| Title |
|---|
| BO YAN 等: "numerical study on dynamics swing of suspension insulator string in overhead transmission line under wind load", 《IEEE TRANSACTIONS ON POWERED DELIVERY》 * |
| MARCELO M.ROCHA等: "A comparison of proper orthogonal decomposition and Monte Carlo simulation of wind pressure data", 《JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS》 * |
| 严波 等: "绝缘子串风偏角风荷载调整系数的研究", 《工程力学》 * |
| 肖刚 等: "计算高可靠性系统失效概率的统计估计蒙特卡洛方法", 《核科学与工程》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108876194A (en) * | 2018-07-16 | 2018-11-23 | 国网福建省电力有限公司 | Power distribution network methods of risk assessment under a kind of typhoon disaster scene |
| CN109471997A (en) * | 2018-10-19 | 2019-03-15 | 海南电网有限责任公司电力科学研究院 | A kind of insulator chain windage yaw distance calculating method |
| CN109471997B (en) * | 2018-10-19 | 2023-06-13 | 海南电网有限责任公司电力科学研究院 | Wind deflection distance calculation method for insulator string |
| CN109271751A (en) * | 2018-11-16 | 2019-01-25 | 重庆科技学院 | A kind of maximum dynamic angle of wind deflection of suspension insulator determines method |
| CN109271751B (en) * | 2018-11-16 | 2022-08-26 | 重庆科技学院 | Method for determining maximum dynamic wind deflection angle of suspension insulator string |
| CN109902351A (en) * | 2019-01-28 | 2019-06-18 | 浙江大学 | A kind of simplified calculation method of ice coating wire dynamic windage yaw |
| CN110378050A (en) * | 2019-07-25 | 2019-10-25 | 重庆科技学院 | A kind of power transmission tower down wind air damping compares calculation method |
| CN110378050B (en) * | 2019-07-25 | 2022-03-22 | 重庆科技学院 | Calculation method for downwind pneumatic damping ratio of power transmission tower |
| CN111157017A (en) * | 2019-11-29 | 2020-05-15 | 湖南大学 | Efficient reliability analysis method for bias value of accelerometer |
| CN114757077A (en) * | 2022-04-21 | 2022-07-15 | 内蒙古电力（集团）有限责任公司内蒙古电力经济技术研究院分公司 | Construction method of wind deflection angle prediction model of double-split line suspension insulator string |
| CN114757077B (en) * | 2022-04-21 | 2024-04-16 | 内蒙古电力（集团）有限责任公司内蒙古电力经济技术研究院分公司 | Construction method of wind deflection angle prediction model of double-split-line suspension insulator string |
| WO2024198505A1 (en) * | 2023-03-30 | 2024-10-03 | 贵州电网有限责任公司 | Wire wind-induced swaying load monitoring method, apparatus and device, and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107977492B (en) | 2020-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107977492A (en) | Based on the non-linear windage yaw reliability degree calculation method of Monte Carlo insulator chain | |
| CN106096105B (en) | Power transmission circuit caused by windage transient response calculation method | |
| Zhou et al. | Methodology that improves water utilization and hydropower generation without increasing flood risk in mega cascade reservoirs | |
| Elshaer et al. | LES evaluation of wind-induced responses for an isolated and a surrounded tall building | |
| CN102663251B (en) | Physical prediction method for wind power station power based on computational fluid mechanics model | |
| CN101661530B (en) | Method for acquiring steady-state equivalent wind speed and generated power in wind power station based on correlation analysis | |
| CN107705032A (en) | A kind of line ice-melting method based on the offline risk evaluation model of power network icing disaster | |
| Tao et al. | Spectral modelling of typhoon winds considering nexus between longitudinal and lateral components | |
| CN106709641A (en) | Monte-Carlo simulation based small interference probability risk analysis and simulation method | |
| Vigueras-Rodríguez et al. | Spectral coherence model for power fluctuations in a wind farm | |
| CN106121913B (en) | A kind of wind turbine arrangement method of marine wind electric field | |
| Zhou et al. | Influence of inflow conditions on turbine loading and wake structures predicted by large eddy simulations using exact geometry | |
| Jeong et al. | Wake impacts on aerodynamic and aeroelastic behaviors of a horizontal axis wind turbine blade for sheared and turbulent flow conditions | |
| CN105069236A (en) | Generalized load joint probability modeling method considering node spatial correlation of wind power plant | |
| Zhang et al. | Failure analysis of large‐scale wind power structure under simulated typhoon | |
| Ke et al. | Non‐Gaussian characteristics and extreme distribution of fluctuating wind pressures on large cylindrical–conical steel cooling towers | |
| CN103248049A (en) | Method of obtaining electrical power system voltage stability domain tangent plane containing DFIG (Doubly Fed Induction Generator) wind power plant | |
| Uchida et al. | New Evaluation Technique for WTG Design Wind Speed Using a CFD‐Model‐Based Unsteady Flow Simulation with Wind Direction Changes | |
| CN105281371A (en) | Telescopic active static safety domain taking wind power generation into account | |
| Olauson et al. | Simulating intra‐hourly wind power fluctuations on a power system level | |
| Fu et al. | LES simulation of hemispherical dome's aerodynamic characteristics in smooth and turbulence boundary layer flows | |
| CN105740516A (en) | Method for determining electric insulation reliability of power transmission line | |
| Rivera‐Arreba et al. | Effect of atmospheric stability on the dynamic wake meandering model applied to two 12 MW floating wind turbines | |
| Liu et al. | Improved local time-stepping schemes for storm surge modeling on unstructured grids | |
| Haouas et al. | Wind farm power forecasting |
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 |