CN101859996A - Method for iron tower anti-lead galloping of coupled system of iron tower, lead and anti-galloping device - Google Patents
- ️Wed Oct 13 2010
Info
-
Publication number
- CN101859996A CN101859996A CN201010187045A CN201010187045A CN101859996A CN 101859996 A CN101859996 A CN 101859996A CN 201010187045 A CN201010187045 A CN 201010187045A CN 201010187045 A CN201010187045 A CN 201010187045A CN 101859996 A CN101859996 A CN 101859996A Authority
- CN
- China Prior art keywords
- lead
- iron tower
- galloping
- unit
- conductor Prior art date
- 2010-05-25 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for iron tower anti-lead galloping of a coupled system of an iron tower, leads and an anti-galloping device. The method comprises the following steps: establishing an iron tower model; setting a lead galloping parameter; establishing a coupled system module of the iron tower, the leads and the anti-galloping device; adding an iron tower node load; setting tower foot node restraining; calculating and solving by utilizing a core calculation module; utilizing a post processing module to directly calculate to output an internal force-time course curve of an iron tower unit; checking the bearing capacity of the iron tower unit according to the calculated output internal force-time course curve of the iron tower unit; and judging whether the iron tower unit needs to be reinforced. The invention considers the coupling effect of the iron tower, the leads and anti-galloping device, accurately studies the lead galloping regulation after the weapon is additionally arranged, can calculate and analyze the influence of lead galloping on the iron tower, studies the bearing capacity of iron tower galloping resistance, provides a basis for reinforcing the iron tower unit, and can satisfy the requirements of scientific research and engineering design.
Description
Technical field
The present invention relates to the method for a kind of iron tower, lead, the anti-conductor galloping of anti-dance device coupled system iron tower.
Background technology
Conductor galloping is the solid coupling nonlinear vibration of stream that ice coating wire produces under the excitation of relatively stable wind.
It is big that ice coating wire is waved amplitude, longer duration, and harm is big.Gently then alternate flashover, damage lead, gold utensil, heavy then cause accidents such as line tripping power failure, wire breaking and tower falling, cause heavy economic losses.Southern ice damage in 2008 just the accident that many iron tower of power transmission line destroy because of conductor galloping occurs.Conductor galloping has contingency on time and space, be difficult to observation.The country that takes place frequently is waved by China, and conductor galloping becomes one of key factor that threatens the overhead transmission line safe and stable operation.
Since the thirties in 20th century, some scholars to conductor galloping mechanism, wave rule and done a large amount of research, proposed the multiple mechanism theory of waving, comprised that eccentric inertia coupling instability Mechanism, the fitful wind that reversing of mechanism, O.Nigol wave mechanism, P.Yu of vertically waving of Den Hartog brings out mechanism etc.Chinese scholars has also been done a large amount of experimental studies to waving of lead, has proposed some anti-dance measures that effectively also found application in engineering practice.
The research of waving in the past is conceived to conductor galloping itself more, and iron tower is approximately certain spring support border, and this research method has following shortcoming: (1), do not consider the coupling effect between iron tower, lead, the anti-dance device; (2), do not consider adjacent shelves lead intercoupling, influencing each other in the process of waving; (3), ice coating wire waves to the influence of iron tower of power transmission line research seldom, the safety research under the condition of waving of iron tower is seldom mentioned.
Therefore, in the process of research conductor galloping, consider the coupling effect between iron tower, lead, the anti-dance device, promote the design theory and the level of the anti-conductor galloping of electric power pylon, improve the anti-dance bearing capacity of iron tower, for conveniently calculating and engineering design service become a very important job.
Summary of the invention
The research that the objective of the invention is to overcome above-mentioned existing conductor galloping is not considered the coupling effect of lead and iron tower and anti-dance device, only iron tower is approximately the deficiency that certain spring supports the border; And provide the method for a kind of iron tower, lead, the anti-conductor galloping of anti-dance device coupled system iron tower.
The present invention utilizes three-dimensional big displacement, big rotation, Nonlinear dynamic explicit Finite Element Method to calculate the process of waving of ice coating wire, set up iron tower, lead, anti-dance device coupling model, consider the coupling effect of iron tower and lead and anti-dance device, what lead behind the anti-dance device was installed in computational analysis additional waves rule and the stress of iron tower under the conductor galloping effect, so that the anti-dance of iron tower bearing capacity is carried out safety evaluation.Can realize iron tower, lead, anti-dance device Fourier Series expansion technique rapid modeling, wave the output of computational analysis and result of calculation, but animation shows the overall process of conductor galloping, with figure, the time history curve of textual form output lead end tension force and the time history curve of iron tower component internal force, as the foundation of the anti-dance design of iron tower.
The objective of the invention is to reach by following measure: the method for iron tower, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is characterized in that it comprises the steps:
(1), set up iron tower model: import the three-dimensional pylon line model of .dxf that vertical type iron tower internal force analysis software TTA generates, with among the .dxf of pylon wired information be converted into the unit information of FEM (finite element) model, with the position of determining unit in the space; The size in each cross section, unit in the FEM (finite element) model that has generated is carried out assignment or directly imported the .out destination file that TTA generates, and the definition unit type is beam element type or bar unit type; In the unit of pylon FEM (finite element) model sectional dimension and after the position in space is determined, again to the cross section of beam element towards adjusting, make that beam element cross-wise direction and actual conditions are identical, promptly obtain an iron tower model; Set up a plurality of iron tower model according to actual span, the corner number of degrees, tower position elevation according to the method described above again;
(2), the conductor galloping parameter is set: the conductor galloping parameter is set, comprises lead division number, division circular section radius, strain insulator string length, lead parameter, icing parameter, wind speed, initial angle of attack θ and aerodynamic coefficient on iron tower model; The cell type and the material properties of definition wire insulation substring, bundle conductor conductor spacer, wherein insulator string is defined as bar unit, and conductor spacer is defined as beam element;
Described lead parameter comprises section radius, steel core section area, aluminium stranded conductor area of section, Unit Weight, modulus of elasticity, the icing parameter is the ice covering thickness of semiellipse type long axis direction, and described aerodynamic coefficient is determined according to icing parameter, conductor cross-section, wind speed and initial angle of attack θ by wind tunnel experiment;
(3), set up iron tower, lead, anti-dance device coupling system model: on the iron tower hanging point, add lead, and the stringing tension force of input lead and dividing elements quantity; Adopt two fens iterative methods that the original shape of lead is calculated, make lead behind weight application and initial tension, arrive the stable equilibrium; Described initial tension T=0.25Tp, Tp are the Fracture Force of lead; Conductor spacer in bundle conductor, adding mutually according to actual pitch; On bundle conductor, add anti-dance device according to engineering reality, wave by GTA and calculate emulation platform input anti-dance device size, weight and installation site;
(4), after setting up iron tower, lead, anti-dance device coupling system model, add the iron tower joint load, the column foot joint constraint is set, translational degree of freedom and the rotational freedom of column foot in X-axis, Y-axis, three directions of Z axle retrained the definition load curve; Time parameter is set: calculate concluding time and result of calculation output time at interval, generate calculation document, utilize the core calculations module to calculate and find the solution;
(5), after calculating finishes, utilize post-processing module, the animation of waving that directly shows lead and pylon by the track of waving in the arbitrary cross section of figure output lead, utilizes the lead modal displacement-time history opisometer of output to calculate waving frequency, waving amplitude of lead; Directly calculate the internal force-time history curve of output iron tower unit;
(6), check the bearing capacity of iron tower unit, judge in view of the above whether the iron tower unit needs to strengthen according to the iron tower element force that calculates output-time history curve.
In technique scheme, the method that the calculating of core calculations module is found the solution in described (4) step is:
The first step: read iron tower, lead, anti-dance device coupling system model;
Second step: in 0＜t＜=200 second, on lead, slowly add carrying force and quiet wind load;
In 200＜t＜=300 second, keep gravity and quiet wind load;
In 300＜t＜T second, add dynamic load, wherein T is the concluding time;
The 3rd step: obtain lead node coordinate, rotational angle, speed, acceleration, angular speed, angular acceleration by the dynamic calculation module;
The 4th step: according to data computation lead unit speed, acceleration, angular speed, angular acceleration and the angle of attack variation in the 3rd step;
The 5th step: according to the data in the 4th step try to achieve lead unit suffered wave load;
The 6th step: the data according to the 5th step are divided load equally the lead unit node;
The 7th step: according to the data in the 6th step will divide equally load on the lead unit node add in second step at 300＜t＜T dynamic load of second;
The 8th step: the dynamic load that obtains according to the 7th step constantly carries out dynamic calculation again to next, enters for the 3rd step; Entered for the 4th step again, the 5th step, the 6th step, the 7th step, the so circulation of the 8th step;
The 9th step: with data record and the output that conforms to the output time of setting in the 3rd step.
In technique scheme, two fens iterative methods of described definite lead original shape are:
Step 1: the shape of differentiate line under the balance of gravity and stringing tension force;
Step 2: the shape of being tried to achieve is carried out conversion, the coordinate difference of all nodes of lead and lead hitch point is multiplied by one less than 1.0 coefficient η, form new wire shape;
Step 3: new wire shape is carried out gravity load calculating, and obtain the lead Horizontal Tension under the gravity;
Step 4: the contrast Horizontal Tension and the stringing tension force that calculate, if Horizontal Tension less than stringing tension force, reduction ratio η, on the contrary increase, this coefficient is determined according to dichotomy;
Step 5: repeating
step2 equates with stringing tension force up to Horizontal Tension to
step4, when difference absolute value during less than 50N, thinks equal, and the shape that obtain this moment is the lead original shape.
The present invention is by setting up iron tower, lead, anti-dance device coupling system model; Lead applied wave load and carry out finite element method (fem) analysis; And utilize post-processing module, with animation, figure and text formatting output result of calculation.
The method of iron tower of the present invention, lead, the anti-conductor galloping of anti-dance device coupled system iron tower, considered the coupling effect of iron tower, lead, anti-dance device, the research of conductor galloping rule has more accuracy behind the anti-dance device to installing additional, also the computable analysis conductor galloping is to the influence of iron tower, the anti-bearing capacity of waving of research iron tower, the foundation of iron tower unit reinforcement is provided, can satisfies the needs of scientific research, engineering design.
Description of drawings
Fig. 1 calculates emulation platform master interface for GTA waves;
Fig. 2 is iron tower, lead, anti-dance device coupling system model figure;
Fig. 3 is the schematic cross-section of ice coating wire;
Fig. 4 is an ice coating wire force diagram schematic diagram;
Fig. 5 is the flow chart of core calculations module of the present invention;
Fig. 6 is the curve chart of lead stretched wire among the present invention;
Fig. 7 is six-multiple conductor middle part node movement locus figure in iron tower, lead, the anti-dance device coupling system model;
Fig. 8 is six-multiple conductor middle part modal displacement-time history curve in iron tower, lead, the anti-dance device coupling system model;
Fig. 9 arranges schematic diagram for the six-multiple conductor anti-dance device.
Embodiment
Describe performance of the present invention in detail below in conjunction with accompanying drawing, but they do not constitute limitation of the invention, only do for example.By explanation, advantage of the present invention will become clear more and understand easily simultaneously.
The method of iron tower of the present invention, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is to set up iron tower, lead, anti-dance device coupling system model earlier, the load that system is subjected to is as excitation then, response by the central difference method solving system in the Nonlinear Dynamic Finite Element Method, wave the post-processing module output result of calculation of calculating emulation platform by GTA at last, whole process is all waved to calculate in the emulation platform (applying for the computer software registration in addition) at GTA and carried out, and is convenient, fast.
Fig. 1 calculates emulation platform master interface for GTA waves, the major function that GTA waves the calculating emulation platform is as follows: import the pylon model from the .dxf file; From the .out file, import iron tower unit sectional dimension; The cross section that the iron tower unit is set towards; Set up single conductor, bundle conductor model; Add panel load and joint constraint, add the curve of load and time control parameter is set; Calculate the original shape of single conductor, bundle conductor; The icing of single conductor, bundle conductor is waved calculating; The time history curve output of iron tower element force and text output; The time history curve output of modal displacement and text output; The animation of whole pylon displacement shows; The animation that track is waved in the arbitrary cross section of lead shows.
Fig. 2 is iron tower, lead, anti-dance device coupling system model figure, and two ends are strain tower, hangs a phase six-multiple conductor, span 400m.
Fig. 3 is the schematic cross-section of ice coating wire, has shown the ice covering thickness 4.0mm of lead steel core diameter 7.04mm, wire diameter 23.94mm, semiellipse type long axis direction among the figure.
Fig. 4 is an ice coating wire force diagram schematic diagram, and among the figure: e is an eccentric throw, G 1Be icing gravity, G 2Be lead gravity, the θ angle is the initial angle of attack,
Be lead y direction speed,
Be lead z direction speed, U is a wind speed, and β is relative wind velocity and y angular separation, U rBe the relative wind velocity that causes because of the lead translation, α is the angle of attack of relative wind velocity, F L(x) be aerodynamic lift, F D(x) be aerodynamic drag, M A(x) be aerodynamic moment.
Fig. 6 is the curve chart of lead stretched wire among the present invention, and A, B are contour hitch point among the figure, and L is a span, and the O point is a lead catenary curve minimum point, and x is the horizontal range that the curve any point is ordered from O, h 0Distance for O point and horizontal reference line.
Fig. 7 is six-multiple conductor middle part node movement locus figure in iron tower, lead, the anti-dance device coupling system model, curve shown in the figure is the movement locus of six-multiple conductor middle part node, comprise since 0 constantly until the y, the z coordinate information that calculate the interior nodes finish time, as can be seen, conductor galloping is based on vertical direction after adding anti-dance device, and horizontal direction is auxilliary.
Fig. 8 is six-multiple conductor middle part modal displacement-time history curve in iron tower, lead, the anti-dance device coupling system model.Among the figure, U x, U y, U zBe respectively node along x, y, z direction of principal axis with close the displacement of direction with U, as can be seen during conductor galloping modal displacement based on the vertical direction displacement.
Fig. 9 arranges schematic diagram for the six-multiple conductor anti-dance device.
The method of iron tower of the present invention, lead, the anti-conductor galloping of anti-dance device coupled system iron tower, it comprises the steps:
(1), set up iron tower model: import the three-dimensional pylon line model of .dxf that vertical type iron tower internal force analysis software TTA generates, with among the .dxf of pylon wired information be converted into the unit information of FEM (finite element) model, with the position of determining unit in the space; The size in each cross section, unit in the FEM (finite element) model that has generated is carried out assignment or directly imported the .out destination file that TTA generates, and the definition unit type is beam element type or bar unit type; In the unit of pylon FEM (finite element) model sectional dimension and after the position in space is determined, again to the cross section of beam element towards adjusting, make that beam element cross-wise direction and actual conditions are identical, promptly obtain an iron tower model; Set up a plurality of iron tower model according to actual span, the corner number of degrees, tower position elevation according to the method described above again;
(2), the conductor galloping parameter is set: the conductor galloping parameter is set, comprises lead division number, division circular section radius, strain insulator string length, lead parameter, icing parameter, wind speed, initial angle of attack θ and aerodynamic coefficient on iron tower model; The cell type and the material properties of definition wire insulation substring, bundle conductor conductor spacer, wherein insulator string is defined as bar unit, and conductor spacer is defined as beam element;
Described lead parameter comprises section radius, steel core section area, aluminium stranded conductor area of section, Unit Weight, modulus of elasticity, the icing parameter is the ice covering thickness of semiellipse type long axis direction, and described aerodynamic coefficient is determined according to icing parameter, conductor cross-section, wind speed and initial angle of attack θ by wind tunnel experiment;
(3), set up iron tower, lead, anti-dance device coupling system model: on the iron tower hanging point, add lead, and the stringing tension force of input lead and dividing elements quantity; Adopt two fens iterative methods that the original shape of lead is calculated, make lead behind weight application and initial tension, arrive the stable equilibrium; Described initial tension T=0.25Tp, Tp are the Fracture Force of lead; Conductor spacer in bundle conductor, adding mutually according to actual pitch; On bundle conductor, add anti-dance device according to engineering reality, wave by GTA and calculate emulation platform input anti-dance device size, weight and installation site;
(4), after setting up iron tower, lead, anti-dance device coupling system model, add the iron tower joint load, the column foot joint constraint is set, translational degree of freedom and the rotational freedom of column foot in X-axis, Y-axis, three directions of Z axle retrained the definition load curve; Time parameter is set: calculate concluding time and result of calculation output time at interval, generate calculation document, utilize the core calculations module to calculate and find the solution;
(5), after calculating finishes, utilize post-processing module, the animation of waving that directly shows lead and pylon by the track of waving in the arbitrary cross section of figure output lead, utilizes the lead modal displacement-time history opisometer of output to calculate waving frequency, waving amplitude of lead; Directly calculate the internal force-time history curve of output iron tower unit;
(6), check the bearing capacity of iron tower unit, judge in view of the above whether the iron tower unit needs to strengthen according to the iron tower element force that calculates output-time history curve.
The method that the calculating of core calculations module is found the solution in described (4) step is:
The first step: read iron tower, lead, anti-dance device coupling system model;
Second step: in 0＜t＜=200 second, on lead, slowly add carrying force and quiet wind load;
In 200＜t＜=300 second, keep gravity and quiet wind load;
In 300＜t＜T second, add dynamic load, wherein T is the concluding time;
The 3rd step: obtain lead node coordinate, rotational angle, speed, acceleration, angular speed, angular acceleration by the dynamic calculation module;
The 4th step: according to data computation lead unit speed, acceleration, angular speed, angular acceleration and the angle of attack variation in the 3rd step;
The 5th step: according to the data in the 4th step try to achieve lead unit suffered wave load;
The 6th step: the data according to the 5th step are divided load equally the lead unit node;
The 7th step: according to the data in the 6th step will divide equally load on the lead unit node add in second step at 300＜t＜T dynamic load of second;
The 8th step: the dynamic load that obtains according to the 7th step constantly carries out dynamic calculation again to next, enters for the 3rd step; Entered for the 4th step again, the 5th step, the 6th step, the 7th step, the so circulation of the 8th step;
The 9th step: with data record and the output that conforms to the output time of setting in the 3rd step.
In technique scheme, two fens iterative methods of described definite lead original shape are:
Step 1: the shape of differentiate line under the balance of gravity and stringing tension force;
Step 2: the shape of being tried to achieve is carried out conversion, the coordinate difference of all nodes of lead and lead hitch point is multiplied by one less than 1.0 coefficient η, form new wire shape;
Step 3: new wire shape is carried out gravity load calculating, and obtain the lead Horizontal Tension under the gravity;
Step 4: the contrast Horizontal Tension and the stringing tension force that calculate, if Horizontal Tension less than stringing tension force, reduction ratio η, on the contrary increase, this coefficient is determined according to dichotomy;
Step 5: repeating
step2 equates with stringing tension force up to Horizontal Tension to step 4, when difference absolute value during less than 50N, thinks equal, and the shape that obtain this moment is the lead original shape.
Accuracy and necessity below by the calculating of computational analysis comparative illustration original shape:
According to the method for the present invention's proposition original shape of determining and the original shape that is assumed to straight line, degree of speeding after adding carrying force and parabola theories are separated and are compared, the stringing tension force of lead Horizontal Tension and setting compares, and the accuracy and the necessity of the inventive method is described.
Stretched wire (lead) accounting equation based on parabola theories is:
y = h 0 + x 2 2 ! h 0
S = x + x 3 3 ! h 0 2
In the formula, y is the ordinate of arbitrfary point, and S is that span mid point O is the lead physical length (as shown in Figure 5) between the point of x to abscissa,
H is the lead Horizontal Tension, and q is uniform load value (being the unit length wire weight when only being subjected to gravity).
Get single LGJ-300/40 lead (national standard), 400 meters spans, the two strain towers in two ends, stringing tension force and be 9.80665 meters of 24000 Ns, acceleration of gravity 2/ second (acceleration of gravity in the dimension 20-40 ° zone, the corresponding Northern Hemisphere).
1, parabola theories
Separate according to parabola theories, the degree of speeding of above-mentioned lead is: 9.258 meters;
2, straight line original shape
Suppose that original shape is a straight line, add the degree of speeding that obtains lead behind the carrying force and be: be 7.72 meters, littler by 16.57% than theoretical solution; The Horizontal Tension that obtains is: 28774.1 Ns, bigger by 19.89% than the stringing tension force of setting;
3, the original shape that calculates according to the inventive method
According to the inventive method, the original shape that calculates is for straight line, adds the degree of speeding that obtains lead behind the carrying force to be: 9.278 meters, bigger by 0.216% than theoretical solution; The Horizontal Tension that obtains is: 23947.3 Ns, littler by 0.22% than the stringing tension force of setting;
From being not difficult to find out, if the supposition original shape is a straight line, not only degree of speeding is inaccurate after adding carrying force, and the Horizontal Tension of lead does not meet with stringing tension force yet; If carry out determining of original shape according to the inventive method, after adding carrying force, lead degree of speeding and theoretical solution error are only 0.216%, and Horizontal Tension that obtains and stringing tension error only are 0.22%.Contrast as can be seen, the original shape computational methods that the present invention proposes not only are necessary, and its calculating accuracy rate is also high.
Below the method for ice coating wire equivalence for single a kind of material and sectional dimension described:
The equivalent area of ice coating wire and equivalent density calculation: actual ice coating wire is formed (lead also can adopt other metal to make) by aluminium, steel, three kinds of materials of icing, therefore needing the ice coating wire equivalence is single a kind of material and sectional dimension, calculates its equivalent area and equivalent density.
The lead equivalent area is:
The equivalent density of ice coating wire is:
ρ＝(M calc+A iceρ ice)/A
In the formula, A IceBe the area of icing on the cross section, ρ IceDensity for ice is taken as 897.6kg/m 3, M CalcCalculated weight for the unit length lead.
Below the beam element that adopts in the core calculations module is described:
Use the discrete ice coating wire of Belytschko beam element in the processing, this beam element uses the corotation technology to handle big rotation, the distortion of beam is divided into rigid body translation, rotation and true strain, rigid body displacement does not produce strain, have only true strain just to produce strain, the big rotation of big displacement that is suitable for the process of waving of lead is handled.
Below the static load in the core calculations module is described:
The static load of lead includes gravity and wind carries the two large divisions.The deadweight of lead can be given by the whole acceleration along gravity direction, and the quiet wind load of lead can obtain by following formula.
Quiet wind load expression formula:
F L ( x ) = 1 2 ρ air U 2 d C L ( α ) F D ( x ) = 1 2 ρ air U 2 d C D ( α ) M A ( x ) = 1 2 ρ air U 2 d 2 C M ( α )
In the formula, U is wind speed, and α is the angle of attack, and aerodynamic coefficient is got the value of wind tunnel experiment, ρ in the formula AirBe atmospheric density, d is a diameter of wire, F L(x), F D(x) and M A(x) be respectively aerodynamic lift, aerodynamic drag and the aerodynamic moment of lead at the x place.
Because it is dynamic explicit analytical method that GTA calculates emulation platform, for as far as possible accurately analog conducting wire carry at quiet wind, poised state under the action of gravity, by slow loaded load and the final method of load a period of time that keeps the static balance state is calculated.
After calculating quiet wind under the certain wind speed and the initial angle of attack according to following formula and carrying, quiet wind is carried on the node that loads on lead by the method for slow loading.This loading procedure is calculated 200s (s represents second) altogether, and the every 1.0s of load increases by 0.5%, becomes 100% in 200s moment load; Kept for 100% load 100s time afterwards, make the scope of lead fluctuation within the acceptable scope; At 300s constantly, cancel Jing Fengzai, and replace, begin dynamically to wave calculating afterwards with dynamically waving load.
Below the load of dynamically waving in the core calculations module is described:
Dynamically wave load is meant that ice coating wire bears when waving load, comprise lead gravity, icing gravity and inertia force, the inertia torque that produces when eccentric torque that eccentric icing causes and icing motion, and dynamic wind load (comprising the aerodynamic lift, aerodynamic drag and the aerodynamic moment that produce by wind).
Under wind action, lead can be subjected to vertically, level and twisting action power (showing as Fig. 4), and the size of these active forces and conductor length are linear relationship, and the aerodynamic lift of supposing unit length is F L(x), aerodynamic drag is F D(x) and aerodynamic moment be M A(x), they are approximate along the distribution of x axle, and supposition is certainly normal usually, is expressed as:
F L ( x ) = 1 2 ρ air ( U - w · ) 2 d [ C L ( α ) cos β - C D ( α ) tan β cos β ] F D ( x ) = 1 2 ρ air ( U - w · ) 2 d [ C L ( α ) tan β cos β + C D ( α ) cos β ] M A ( x ) = 1 2 ρ air ( U - w · ) 2 d 2 C M ( α ) 2 cos 2 β
ρ in the formula AirBe atmospheric density, d is a diameter of wire,
Be lead z direction speed, U is a wind speed, and β is relative wind velocity and y angular separation, and α is the angle of attack of relative wind velocity, aerodynamic coefficient C L(x), C D(x) and C M(x) can determine by the wind tunnel experiment of corresponding ice type.
Comparison example: have, the comparison of no anti-dance device
In order to check the anti-dance effect of anti-dance device, establish, do not have the two strain towers-single-phase six-multiple conductor model (6 division LGJ-300/40 leads, 400 meters spans) of anti-dance device, calculated the wave situation of two kinds of models under 8m/s wind speed, 30 ° of initial angle of attack conditions.
Result of calculation:
A. time of experiencing in the lead stage of dancing is increased to 5000s from 800s, illustrates the formation that two parallel wielded-proof device can effectively suppress conductor galloping is installed;
B. the maximum amplitude of conductor galloping is reduced to 1.63m from 3.84m, reduces degree and reaches 60%, illustrates that two parallel wielded-proof device is installed has effectively reduced the amplitude of waving;
C. the frequency of conductor galloping became 20 times/minute from 21-22 time/minute, illustrated the frequency of waving that two parallel wielded-proof device can change lead is installed.
By have in this example, the contrast of the conductor galloping rule of no anti-dance device, the method for iron tower of the present invention, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is safe and reliable, convenient and practical.
Claims (3)
1. the method for iron tower, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is characterized in that it comprises the steps:
(1), set up iron tower model: import the three-dimensional pylon line model of .dxf that vertical type iron tower internal force analysis software TTA generates, with among the .dxf of pylon wired information be converted into the unit information of FEM (finite element) model, with the position of determining unit in the space; The size in each cross section, unit in the FEM (finite element) model that has generated is carried out assignment or directly imported the .out destination file that TTA generates, and the definition unit type is beam element type or bar unit type; In the unit of pylon FEM (finite element) model sectional dimension and after the position in space is determined, again to the cross section of beam element towards adjusting, make that beam element cross-wise direction and actual conditions are identical, promptly obtain an iron tower model; Set up a plurality of iron tower model according to actual span, the corner number of degrees, tower position elevation according to the method described above again;
(2), the conductor galloping parameter is set: the conductor galloping parameter is set, comprises lead division number, division circular section radius, strain insulator string length, lead parameter, icing parameter, wind speed, initial angle of attack θ and aerodynamic coefficient on iron tower model; The cell type and the material properties of definition wire insulation substring, bundle conductor conductor spacer, wherein insulator string is defined as bar unit, and conductor spacer is defined as beam element;
Described lead parameter comprises section radius, steel core section area, aluminium stranded conductor area of section, Unit Weight, modulus of elasticity, the icing parameter is the ice covering thickness of semiellipse type long axis direction, and described aerodynamic coefficient is determined according to icing parameter, conductor cross-section, wind speed and initial angle of attack θ by wind tunnel experiment;
(3), set up iron tower, lead, anti-dance device coupling system model: on the iron tower hanging point, add lead, and the stringing tension force of input lead and dividing elements quantity; Adopt two fens iterative methods that the original shape of lead is calculated, make lead behind weight application and initial tension, arrive the stable equilibrium; Described initial tension T=0.25Tp, Tp are the Fracture Force of lead; Conductor spacer in bundle conductor, adding mutually according to actual pitch; On bundle conductor, add anti-dance device according to engineering reality, wave by GTA and calculate emulation platform input anti-dance device size, weight and installation site;
(4), after setting up iron tower, lead, anti-dance device coupling system model, add the iron tower joint load, the column foot joint constraint is set, translational degree of freedom and the rotational freedom of column foot in X-axis, Y-axis, three directions of Z axle retrained the definition load curve; Time parameter is set: calculate concluding time and result of calculation output time at interval, generate calculation document, utilize the core calculations module to calculate and find the solution;
(5), after calculating finishes, utilize post-processing module, the animation of waving that directly shows lead and pylon by the track of waving in the arbitrary cross section of figure output lead, utilizes the lead modal displacement-time history opisometer of output to calculate waving frequency, waving amplitude of lead; Directly calculate the internal force-time history curve of output iron tower unit;
(6), check the bearing capacity of iron tower unit, judge in view of the above whether the iron tower unit needs to strengthen according to the iron tower element force that calculates output-time history curve.
2. the method for iron tower according to claim 1, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is characterized in that the method that the calculating of core calculations module is found the solution in described (4) step is:
The first step: read iron tower, lead, anti-dance device coupling system model;
Second step: in 0＜t＜=200 second, on lead, slowly add carrying force and quiet wind load;
In 200＜t＜=300 second, keep gravity and quiet wind load;
In 300＜t＜T second, add dynamic load, wherein T is the concluding time;
The 3rd step: obtain lead node coordinate, rotational angle, speed, acceleration, angular speed, angular acceleration by the dynamic calculation module;
The 4th step: according to data computation lead unit speed, acceleration, angular speed, angular acceleration and the angle of attack variation in the 3rd step;
The 5th step: according to the data in the 4th step try to achieve lead unit suffered wave load;
The 6th step: the data according to the 5th step are divided load equally the lead unit node;
The 7th step: according to the data in the 6th step will divide equally load on the lead unit node add in second step at 300＜t＜T dynamic load of second;
The 8th step: the dynamic load that obtains according to the 7th step constantly carries out dynamic calculation again to next, enters for the 3rd step; Entered for the 4th step again, the 5th step, the 6th step, the 7th step, the so circulation of the 8th step;
The 9th step: with data record and the output that conforms to the output time of setting in the 3rd step.
3. the method for iron tower according to claim 1, lead, the anti-conductor galloping of anti-dance device coupled system iron tower is characterized in that two fens iterative methods of described definite lead original shape are:
Step 1: the shape of differentiate line under the balance of gravity and stringing tension force;
Step 2: the shape of being tried to achieve is carried out conversion, the coordinate difference of all nodes of lead and lead hitch point is multiplied by one less than 1.0 coefficient η, form new wire shape;
Step 3: new wire shape is carried out gravity load calculating, and obtain the lead Horizontal Tension under the gravity;
Step 4: the contrast Horizontal Tension and the stringing tension force that calculate, if Horizontal Tension less than stringing tension force, reduction ratio η, on the contrary increase, this coefficient is determined according to dichotomy;
Step 5: repeating step 2 equates with stringing tension force up to Horizontal Tension to step 4, when difference absolute value during less than 50N, thinks equal, and the shape that obtain this moment is the lead original shape.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010187045 CN101859996B (en) | 2010-05-25 | 2010-05-25 | Method for iron tower anti-lead galloping of coupled system of iron tower, lead and anti-galloping device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010187045 CN101859996B (en) | 2010-05-25 | 2010-05-25 | Method for iron tower anti-lead galloping of coupled system of iron tower, lead and anti-galloping device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101859996A true CN101859996A (en) | 2010-10-13 |
| CN101859996B CN101859996B (en) | 2013-06-19 |
Family
ID=42945724
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010187045 Active CN101859996B (en) | 2010-05-25 | 2010-05-25 | Method for iron tower anti-lead galloping of coupled system of iron tower, lead and anti-galloping device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101859996B (en) |
Cited By (12)
* Cited by examiner, † Cited by third party| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101975565A (en) * | 2010-10-15 | 2011-02-16 | 合肥工业大学 | GPS-based transmission conductor waving monitoring system and monitoring method |
| CN102623939A (en) * | 2012-02-06 | 2012-08-01 | 广东电网公司电力科学研究院 | A wind-induced vibration control method for transmission lines based on simulated random gust field calculation |
| CN104504225A (en) * | 2015-01-20 | 2015-04-08 | 华北电力大学（保定） | Single-column stay wire tower torsion frequency estimation method based on single degree of freedom model |
| CN104638559A (en) * | 2013-11-07 | 2015-05-20 | 国家电网公司 | Method for installing 500kV compact transmission line phase-to-phase spacers |
| CN104833393A (en) * | 2015-05-12 | 2015-08-12 | 西安工程大学 | Transmission line iced conductor galloping aerodynamic parameter monitoring device and monitoring method |
| CN104913900A (en) * | 2015-06-03 | 2015-09-16 | 西安工程大学 | Measurement method of power transmission line iced conductor galloping aerodynamic parameters |
| CN106199232A (en) * | 2015-11-09 | 2016-12-07 | 全玉生 | The positive sequence impedance method of high pressure overhead power line dynamic compatibilization |
| CN106370373A (en) * | 2016-11-03 | 2017-02-01 | 中国电力科学研究院 | Compact transmission line wire anti-dance test device |
| CN109033666A (en) * | 2018-08-07 | 2018-12-18 | 国网河南省电力公司电力科学研究院 | Ice coating wire based on beam element waves ANSYS method for numerical simulation |
| CN110375949A (en) * | 2019-08-08 | 2019-10-25 | 南方电网科学研究院有限责任公司 | Wire wind tunnel test model and manufacturing method thereof |
| CN111753456A (en) * | 2020-06-30 | 2020-10-09 | 国网河南省电力公司电力科学研究院 | A method and system for using a spacer bar suitable for preventing wire overturning |
| CN114279495A (en) * | 2021-12-16 | 2022-04-05 | 国网江苏省电力有限公司盐城供电分公司 | A monitoring system for galloping overhead lines |
Families Citing this family (1)
* Cited by examiner, † Cited by third party| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107526867A (en) * | 2017-07-17 | 2017-12-29 | 北京道亨时代科技有限公司 | A kind of processing method of steel tower auxiliary material |
Citations (2)
* Cited by examiner, † Cited by third party| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005213822A (en) * | 2004-01-28 | 2005-08-11 | Chugoku Electric Power Co Inc:The | Design-assist method for uneven leg in transmission line steel tower, design-assist apparatus for uneven leg, computer program and program storage medium |
| CN101577408A (en) * | 2009-03-25 | 2009-11-11 | 中国电力科学研究院 | Comprehensive analysis method of reliability of transmission line tower structure |
-
2010
- 2010-05-25 CN CN 201010187045 patent/CN101859996B/en active Active
Patent Citations (2)
* Cited by examiner, † Cited by third party| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005213822A (en) * | 2004-01-28 | 2005-08-11 | Chugoku Electric Power Co Inc:The | Design-assist method for uneven leg in transmission line steel tower, design-assist apparatus for uneven leg, computer program and program storage medium |
| CN101577408A (en) * | 2009-03-25 | 2009-11-11 | 中国电力科学研究院 | Comprehensive analysis method of reliability of transmission line tower structure |
Non-Patent Citations (1)
* Cited by examiner, † Cited by third party| Title |
|---|
| 杨靖波 等: "输电线路铁塔导线耦合振动动力特性研究", 《振动 测试与诊断》, vol. 28, no. 2, 30 June 2008 (2008-06-30), pages 147 - 150 * |
Cited By (19)
* Cited by examiner, † Cited by third party| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101975565B (en) * | 2010-10-15 | 2013-04-03 | 合肥工业大学 | GPS-based transmission conductor waving monitoring system and monitoring method |
| CN101975565A (en) * | 2010-10-15 | 2011-02-16 | 合肥工业大学 | GPS-based transmission conductor waving monitoring system and monitoring method |
| CN102623939A (en) * | 2012-02-06 | 2012-08-01 | 广东电网公司电力科学研究院 | A wind-induced vibration control method for transmission lines based on simulated random gust field calculation |
| CN102623939B (en) * | 2012-02-06 | 2015-09-09 | 广东电网公司电力科学研究院 | A kind of power transmission line transportation work style calculated based on simulation stochastic matrix wind field is shaken administering method |
| CN104638559B (en) * | 2013-11-07 | 2017-02-22 | 国家电网公司 | Method for installing 500kV compact transmission line phase-to-phase spacers |
| CN104638559A (en) * | 2013-11-07 | 2015-05-20 | 国家电网公司 | Method for installing 500kV compact transmission line phase-to-phase spacers |
| CN104504225A (en) * | 2015-01-20 | 2015-04-08 | 华北电力大学（保定） | Single-column stay wire tower torsion frequency estimation method based on single degree of freedom model |
| CN104504225B (en) * | 2015-01-20 | 2017-09-22 | 华北电力大学（保定） | Single-column Guywire tower torsion frequency evaluation method based on one degree of freedom modeling |
| CN104833393A (en) * | 2015-05-12 | 2015-08-12 | 西安工程大学 | Transmission line iced conductor galloping aerodynamic parameter monitoring device and monitoring method |
| CN104913900B (en) * | 2015-06-03 | 2017-10-31 | 西安工程大学 | The assay method of powerline ice-covering conductor galloping aerodynamic parameter |
| CN104913900A (en) * | 2015-06-03 | 2015-09-16 | 西安工程大学 | Measurement method of power transmission line iced conductor galloping aerodynamic parameters |
| CN106199232A (en) * | 2015-11-09 | 2016-12-07 | 全玉生 | The positive sequence impedance method of high pressure overhead power line dynamic compatibilization |
| CN106370373A (en) * | 2016-11-03 | 2017-02-01 | 中国电力科学研究院 | Compact transmission line wire anti-dance test device |
| CN109033666A (en) * | 2018-08-07 | 2018-12-18 | 国网河南省电力公司电力科学研究院 | Ice coating wire based on beam element waves ANSYS method for numerical simulation |
| CN109033666B (en) * | 2018-08-07 | 2022-09-09 | 国网河南省电力公司电力科学研究院 | Beam-unit-based ANSYS numerical simulation method for galloping of ice-coated wire |
| CN110375949A (en) * | 2019-08-08 | 2019-10-25 | 南方电网科学研究院有限责任公司 | Wire wind tunnel test model and manufacturing method thereof |
| CN111753456A (en) * | 2020-06-30 | 2020-10-09 | 国网河南省电力公司电力科学研究院 | A method and system for using a spacer bar suitable for preventing wire overturning |
| CN111753456B (en) * | 2020-06-30 | 2023-02-03 | 国网河南省电力公司电力科学研究院 | A method and system for using a spacer bar suitable for preventing wire overturning |
| CN114279495A (en) * | 2021-12-16 | 2022-04-05 | 国网江苏省电力有限公司盐城供电分公司 | A monitoring system for galloping overhead lines |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101859996B (en) | 2013-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101859996B (en) | 2013-06-19 | Method for iron tower anti-lead galloping of coupled system of iron tower, lead and anti-galloping device |
| CN101847855B (en) | 2013-08-21 | Conductor galloping resistant method of iron towers of overhead power transmission lines |
| Ma et al. | 2021 | Component-based fragility analysis of transmission towers subjected to hurricane wind load |
| Zhou et al. | 2016 | Study on galloping behavior of iced eight bundle conductor transmission lines |
| CN103942417B (en) | 2017-01-04 | Ice coating wire ice-shedding emulation test method |
| Wardlaw et al. | 1975 | Wind tunnel and analytical investigations into the aeroelastic behaviour of bundled conductors |
| Hamada et al. | 2017 | The response of a guyed transmission line system to boundary layer wind |
| Wei et al. | 2012 | Modal analysis of a cable-stayed bridge |
| CN115391881B (en) | 2023-04-18 | Bridge tower wake region sling wind-induced vibration numerical value prediction method |
| CN104504225B (en) | 2017-09-22 | Single-column Guywire tower torsion frequency evaluation method based on one degree of freedom modeling |
| Chen et al. | 2018 | Experimental and theoretical investigation of galloping of transversely inclined slender prisms |
| Xie et al. | 2019 | Wind tunnel test of an aeroelastic model of a catenary system for a high-speed railway in China |
| Cai et al. | 2020 | Investigation on galloping of D-shape iced 6-bundle conductors in transmission tower line |
| Asgari et al. | 2015 | Optimization of pre-tensioning cable forces in highly redundant cable-stayed bridges |
| Sun et al. | 2024 | Fatigue of suspender anchorages under axial and bending loads of suspension bridges |
| Konstantakopoulos et al. | 2010 | A mathematical model for a combined cable system of bridges |
| CN104636603B (en) | 2017-08-29 | A kind of two-freedom computational methods of extra-high voltage single-column Guywire tower torsion frequency |
| Shevchenko et al. | 2015 | Refinement of wind loads on lattice support structures of the intersystem overhead power transmission lines 750 kV |
| CN103050909B (en) | 2016-02-10 | A kind of inter-phase spacer installation method and system thereof |
| Templin et al. | 1981 | Design and performance of a multi-degree-of-freedom aeroelastic building model |
| Dua et al. | 2015 | Dynamic analysis of overhead transmission line under turbulent wind loading |
| Shklyarchuk et al. | 2007 | Nonlinear Aeroelastic vibrations and galloping of iced conductor lines under wind. |
| Li et al. | 2024 | Theoretical algorithm for maximum jump height of the conductor after ice-shedding considering elevation angle and nonlinearity |
| CN112036115A (en) | 2020-12-04 | Icing conductor galloping simulation test method |
| Wang et al. | 2017 | Numerical computation for vibration characteristics of long-span bridges with considering vehicle-wind coupling excitations based on finite element and neural network models |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 2010-10-13 | C06 | Publication | |
| 2010-10-13 | PB01 | Publication | |
| 2010-11-24 | C10 | Entry into substantive examination | |
| 2010-11-24 | SE01 | Entry into force of request for substantive examination | |
| 2013-06-19 | C14 | Grant of patent or utility model | |
| 2013-06-19 | GR01 | Patent grant | |
| 2015-02-11 | C56 | Change in the name or address of the patentee |
Owner name: CENTRAL SOUTHERN CHINA ELECTRIC POWER DESIGN INSTI Free format text: FORMER NAME: CENTRAL SOUTHERN CHINA ELECTRIC POWER DESIGN INSTITUTE, CHINA POWER ENGINEERING CONSULTING GROUP COR |
| 2015-02-11 | CP03 | Change of name, title or address |
Address after: 430071 Wuchang District, Hubei, South Central Road No. two, No. 12, No. Patentee after: Co., Ltd of Central Southern China Electric Power Design Institute, China Power Engineering Consulting Group Corporation Address before: 430071 No. 668 Democratic Road, Wuchang District, Hubei, Wuhan (Central South Electric Power Design Institute) Patentee before: Central Southern China Electric Power Design Institute, China Power Engineering Consulting Group Cor |