CN111241731A - Transmission conductor shape finding method based on horizontal tension - Google Patents

Transmission conductor shape finding method based on horizontal tension Download PDF

Info

Publication number
CN111241731A
CN111241731A CN202010014812.7A CN202010014812A CN111241731A CN 111241731 A CN111241731 A CN 111241731A CN 202010014812 A CN202010014812 A CN 202010014812A CN 111241731 A CN111241731 A CN 111241731A
Authority
CN
China
Prior art keywords
transmission conductor
power transmission
tension
conductor
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010014812.7A
Other languages
Chinese (zh)
Other versions
CN111241731B (en
Inventor
吕洪坤
刘孟龙
罗坤
樊建人
汪明军
应明良
池伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN202010014812.7A priority Critical patent/CN111241731B/en
Publication of CN111241731A publication Critical patent/CN111241731A/en
Application granted granted Critical
Publication of CN111241731B publication Critical patent/CN111241731B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a transmission conductor shape-finding method based on horizontal tension, which comprises the following steps: calculating the initial position of a transmission conductor node between two towers by using a catenary formula according to the positions of the hanging line points of the two transmission towers; simulating a power transmission conductor by using a non-linear unit which is only pulled, and setting initial strain; loading a gravity load, and carrying out shape finding analysis on a transmission conductor finite element model under the gravity condition; adjusting the initial strain of the lead according to the difference value of the annual average running horizontal tension and the horizontal tension of the middle node of the lead; and then repeating the shape finding calculation until the difference value between the horizontal tension of the middle node of the conductor and the annual average running horizontal tension is smaller than a threshold value, stopping the calculation, and completing the position shape finding analysis of the power transmission conductor under the gravity condition. According to the invention, the shape finding calculation of the conductor is carried out on the basis of the power transmission conductor model, the iterative calculation difficulty is greatly reduced, the annual average running tension is used as a convergence basis for iterative calculation, the calculation result is more accurate, and the method is more suitable for actual conditions.

Description

Transmission conductor shape finding method based on horizontal tension
Technical Field
The invention relates to a power transmission conductor shape finding method, in particular to a power transmission conductor shape finding method taking horizontal tension in a conductor as an iterative calculation convergence judgment standard.
Background
The transmission line is a lifeline project related to the national civilization, and the transmission conductor is an important component in the transmission line and has direct influence on the overall structural characteristics of the transmission line, so that the study on the structural characteristics of the transmission conductor is an essential part for analyzing the overall transmission line.
The power transmission conductor is a typical flexible cable structure, and has the characteristics of bearing tensile force only and bearing rigidity of external load only when the power transmission conductor is in a stress tension state, so that the premise of analyzing the structural characteristics of the power transmission conductor is to determine the initial position and the initial stress state of the power transmission conductor after being tensioned under the action of gravity.
At present, finite element methods are used for carrying out shape finding analysis on the transmission conductor, and the most common methods are direct iteration and small elastic modulus methods. The direct iteration method is that under the state of very small initial strain, an iteration convergence condition is set, the initial state of the transmission line straight line is directly subjected to iterative calculation, and a finite element model is continuously updated, so that the initial state of the transmission line under the action of gravity is obtained; the small elastic modulus method is characterized in that the elastic modulus of the material of the power transmission conductor is set to be an elastic modulus far smaller than an actual value, the approximate form of the power transmission conductor under the action of gravity is obtained after initial calculation, the actual elastic modulus value is restored under the form to carry out iterative calculation, and finally the final shape finding result of the actual power transmission conductor is obtained through calculation. And the common direct iteration method and the small elastic modulus method fully embody the characteristic of difficult convergence of the iterative computation in the shape finding computation process of complex models such as a power transmission tower line coupling system and the like.
Disclosure of Invention
The invention aims to provide a power transmission conductor shape finding method suitable for a complex coupling model aiming at the defects of the prior art, which is easy to converge and simple to operate, and the solving result can meet the actual precision requirement of engineering.
The purpose of the invention is realized by the following technical scheme: a power transmission conductor form-finding method based on horizontal tension mainly comprises the following steps:
step 1, determining the positions of two power transmission tower line hanging points, equally dividing the horizontal span between the two line hanging points into n parts to obtain n-1 nodes of a power transmission conductor between the two line hanging points, establishing a coordinate system by taking the horizontal projection direction of a connecting line of the two line hanging points as an x axis and the vertical projection direction as a z axis, and calculating the initial coordinate position of each node of the power transmission conductor between the two line hanging points according to line hanging point position data, physical property parameters of the power transmission conductor and a catenary formula;
step 2, constructing a transmission conductor finite element model according to the initial coordinates of each node, and setting the initial strain of the transmission conductor;
step 3, loading gravity on the integral model in the finite element model of the power transmission conductor, setting position constraints at two hanging points of the power transmission conductor, and then performing nonlinear finite element calculation;
step 4, when n is an odd number, extracting the horizontal tension result of the (n + 1)/2-th wire node; when n is even number, the horizontal tension result of the n/2 th wire node is extracted. Performing difference value operation on the horizontal tension result of the node and the annual average running tension of the power transmission conductor, and directly obtaining a power transmission conductor shape finding result under the gravity condition if the difference value is smaller than a threshold value; otherwise, judging the level tension result of the node and the annual average running tension of the power transmission conductor, and reducing the initial strain set by the conductor if the level running tension result of the node is greater than the annual average running tension of the power transmission conductor; if the horizontal tension result of the node is less than the annual average running tension of the transmission conductor, increasing the initial strain set by the conductor;
and 5, repeating the step 3 and the step 4 according to the result of the step 4 until the difference value between the horizontal tension of the node and the annual average running tension of the power transmission conductor in the nonlinear finite element calculation result is smaller than a threshold value, finishing the shape finding calculation of the power transmission conductor, and obtaining the shape finding result of the power transmission conductor under the gravity condition.
Further, the catenary equation for establishing the transmission conductor initial finite element model in step 1 is as follows:
Figure BDA0002358484790000021
Figure BDA0002358484790000022
Figure BDA0002358484790000023
wherein x and z are calculated values of the node coordinates of each wire, and Q is the dead weight of the wire in unit length; h is the horizontal tension of the transmission line, namely the annual average running tension; l is the horizontal span of the wire; and c is the height difference of two ends of the wire.
Further, ANSYS software is adopted for establishing the finite element model in the step 2, only a link10 unit under tension in the finite element software ANSYS is adopted for simulating the power transmission conductor, and the initial strain is set to be epsilon0=1×10-5And sets the elastic modulus, density, poisson's ratio of link10 cell according to actual power conductor parameters.
Further, the initial strain setting is adjusted to an initial strain of one tenth in the increase and decrease amount, i.e., ∈ in step 40/10。
Further, the threshold values in the step 4 and the step 5 are set according to the precision requirement, if the difference value between the node level tension result and the annual average running tension of the power transmission conductor is less than 5% of the annual average running tension, the calculation is finished, and if the difference value is greater than the annual average running tension, the calculation is continued circularly.
The invention has the beneficial effects that: compared with the prior art, the method has the advantages that the initial model of the power transmission conductor is established by the catenary formula, the conductor shape finding calculation is carried out on the basis of the initial model, the iterative calculation difficulty is greatly reduced, the annual average running tension is used as a convergence basis for carrying out iterative calculation, the calculation result is more accurate, and the method is more in line with the actual situation.
Drawings
FIG. 1 is a general flow diagram of the method;
FIG. 2 is a coordinate system used in wire node calculation;
figure 3 is a diagram of the location of the nodes of the transmission conductor;
fig. 4 is a transmission conductor initial finite element model.
Detailed Description
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
As shown in fig. 1, the method for finding the shape of a power transmission conductor according to horizontal tension mainly comprises the following steps:
step 1, determining the positions of two power transmission tower line hanging points, equally dividing the horizontal span between the two line hanging points into n parts to obtain n-1 nodes of a power transmission conductor between the two line hanging points, establishing a coordinate system by taking the horizontal projection direction of a connecting line of the two line hanging points as an x axis and the vertical projection direction as a z axis, and calculating the initial coordinate position of each node of the power transmission conductor between the two line hanging points according to line hanging point position data, physical property parameters of the power transmission conductor and a catenary formula;
the catenary equation is as follows:
Figure BDA0002358484790000031
Figure BDA0002358484790000032
Figure BDA0002358484790000033
wherein x and z are calculated values of the node coordinates of each wire, and Q is the dead weight of the wire in unit length; h is the horizontal tension of the transmission line, namely the annual average running tension; l is the horizontal span of the wire; and c is the height difference of two ends of the wire.
Step 2, constructing a finite element model of the power transmission conductor by adopting ANSYS software according to the initial coordinates of each node, simulating the power transmission conductor by adopting a link10 unit which is only pulled in the finite element software ANSYS, and setting the initial strain epsilon of the power transmission conductor0=1×10-5(ii) a And sets the elastic modulus, density, poisson's ratio of link10 unit according to actual power conductor parameters.
Step 3, loading gravity on the integral model in the finite element model of the power transmission conductor, setting position constraints at two hanging points of the power transmission conductor, and then performing nonlinear finite element calculation;
step 4, when n is an odd number, extracting the horizontal tension result of the (n + 1)/2-th wire node; when n is even number, the horizontal tension result of the n/2 th wire node is extracted. Performing difference value operation on the horizontal tension result of the node and the annual average running tension of the power transmission conductor, and directly obtaining a power transmission conductor shape finding result under the gravity condition if the difference value is smaller than a threshold value; otherwise, judging the level tension result of the node and the annual average running tension of the power transmission conductor, and reducing one tenth of the initial strain of the conductor setting if the level running tension result of the node is greater than the annual average running tension of the power transmission conductor; if the node horizontal tension result is less than the annual average operating tension of the power transmission conductor, increasing one-tenth of the initial strain of the conductor setting; and the threshold is set according to the precision requirement, if the difference value between the node horizontal tension result and the annual average running tension of the power transmission conductor is less than 5% of the annual average running tension, the calculation is finished, and if the difference value is greater than the annual average running tension, the calculation is continued circularly.
And 5, repeating the step 3 and the step 4 according to the result of the step 4 until the difference value between the horizontal tension of the node and the annual average running tension of the power transmission conductor in the nonlinear finite element calculation result is smaller than a threshold value, finishing the shape finding calculation of the power transmission conductor, and obtaining the shape finding result of the power transmission conductor under the gravity condition.
The method for finding the shape of the power transmission conductor is further explained by combining the specific embodiment as follows:
step 1, determining the positions of the hanging wire points of two transmission towers, and calculating the initial position of each node of a transmission conductor between the two towers according to the actual hanging wire point position data and the following catenary formula;
Figure BDA0002358484790000041
Figure BDA0002358484790000042
Figure BDA0002358484790000043
wherein x and z are calculated values of the node coordinates of each wire, and Q is the dead weight of the wire in unit length; h is the horizontal tension of the transmission line, namely the annual average running tension, and the annual average running tension of the lead is obtained by multiplying the breaking force of the ground lead wire by a safety coefficient of 0.95 and an average running tension coefficient of 0.25 in the current engineering; l is the horizontal span of the wire; and c is the height difference of two ends of the wire. The coordinate system uses the horizontal projection direction of the connecting line of the two hanging line points as the x axis, and the vertical projection direction as the z axis, as shown in fig. 2.
And if the span and the height difference between the two wire hanging points are 200m and 20m, L is 200 and c is 20. Other parameters are set according to the parameters of the actual leads JL/G1A-400/35, the self weight Q of the unit lead is 13.2N/m, and the annual average horizontal tension H is 24622N. The initial positions of 199 nodes of the power transmission conductor except for the hanging wire point are calculated by combining the parameters with a catenary formula, and a power transmission conductor finite element model is constructed, as shown in fig. 3 and 4.
Step 2, constructing a finite element model of the transmission conductor by using a nonlinear unit link10 in finite element software ANSYS, setting initial strain in a real constant of a link10 unit, setting the initial strain to be 1E-5, and simultaneously setting other conductor physical parameters, wherein the Poisson ratio is 0.3, and the elastic modulus is 6.5E10 Pa;
step 3, loading gravity on the integral model in the finite element model of the power transmission conductor, setting position constraints at two hanging points of the power transmission conductor, and then performing nonlinear finite element calculation;
step 4, extracting a horizontal tension result of nonlinear finite element analysis of a 100 th node in the middle of the power transmission conductor under the gravity condition, comparing the extracted horizontal tension result of the middle node with a design value of annual average running tension of the power transmission conductor, wherein the calculation result of the horizontal tension of the node is far smaller than the design value of annual average running tension of the power transmission conductor 24622N, and increasing the initial strain 1E-06;
and 5, repeating the step 3 and the step 4 according to the result of the step 4 until the difference value between the horizontal tension of the node and the annual average running tension of the power transmission conductor in the calculation result is smaller than the annual average running tension design value of the power transmission conductor by 5 percent, namely 1231.1N, finishing the conductor form-finding calculation, and obtaining the calculation result of the position of the power transmission conductor under the gravity condition.
The above-described embodiments are intended to illustrate rather than to limit the invention, and any modifications and variations of the present invention are within the spirit of the invention and the scope of the appended claims.

Claims (5)

1. A power transmission conductor form-finding method based on horizontal tension is characterized by mainly comprising the following steps:
step 1, determining the positions of two power transmission tower line hanging points, equally dividing the horizontal span between the two line hanging points into n parts to obtain n-1 nodes of a power transmission conductor between the two line hanging points, establishing a coordinate system by taking the horizontal projection direction of a connecting line of the two line hanging points as an x axis and the vertical projection direction as a z axis, and calculating the initial coordinate position of each node of the power transmission conductor between the two line hanging points according to line hanging point position data, physical property parameters of the power transmission conductor and a catenary formula; .
Step 2, constructing a transmission conductor finite element model according to the initial coordinates of each node, and setting the initial strain of the transmission conductor;
step 3, loading gravity on the integral model in the finite element model of the power transmission conductor, setting position constraints at two hanging points of the power transmission conductor, and then performing nonlinear finite element calculation;
step 4, when n is an odd number, extracting the horizontal tension result of the (n + 1)/2-th wire node; when n is even number, the horizontal tension result of the n/2 th wire node is extracted. Performing difference value operation on the horizontal tension result of the node and the annual average running tension of the power transmission conductor, and directly obtaining a power transmission conductor shape finding result under the gravity condition if the difference value is smaller than a threshold value; otherwise, judging the level tension result of the node and the annual average running tension of the power transmission conductor, and reducing the initial strain set by the conductor if the level running tension result of the node is greater than the annual average running tension of the power transmission conductor; if the horizontal tension result of the node is less than the annual average running tension of the transmission conductor, increasing the initial strain set by the conductor;
and 5, repeating the step 3 and the step 4 according to the result of the step 4 until the difference value between the horizontal tension of the node and the annual average running tension of the power transmission conductor in the nonlinear finite element calculation result is smaller than a threshold value, finishing the shape finding calculation of the power transmission conductor, and obtaining the shape finding result of the power transmission conductor under the gravity condition.
2. A method of horizontal tension transmission conductor form finding as claimed in claim 1, wherein the catenary equation for establishing the initial finite element model of the transmission conductor in step 1 is as follows:
Figure FDA0002358484780000011
Figure FDA0002358484780000012
Figure FDA0002358484780000013
wherein x and z are calculated values of the node coordinates of each wire, and Q is the dead weight of the wire in unit length; h is the horizontal tension of the transmission line, namely the annual average running tension; l is the horizontal span of the wire; and c is the height difference of two ends of the wire.
3. The method of claim 1 wherein the step 2 of creating a finite element model uses ANSYS software and simulating the conductor using the tension-only link10 element of the finite element software ANSYS, the initial strain set to ε0=1×10-5And sets the elastic modulus, density, poisson's ratio of link10 cell according to actual power conductor parameters.
4. The method of claim 1, wherein the initial strain is adjusted in step 4 by a factor of one tenth, i.e., ∈ to0/10。
5. The conductor form-finding method according to the horizontal tension of claim 1, wherein the threshold values in the steps 4 and 5 are set according to the precision requirement, and if the difference value between the node horizontal tension result and the annual average running tension of the power transmission conductor is less than 5% of the annual average running tension, the calculation is finished, and if the difference value is greater than the annual average running tension, the calculation is continued circularly.
CN202010014812.7A 2020-01-07 2020-01-07 Power transmission wire shape finding method based on horizontal tension Active CN111241731B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010014812.7A CN111241731B (en) 2020-01-07 2020-01-07 Power transmission wire shape finding method based on horizontal tension

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010014812.7A CN111241731B (en) 2020-01-07 2020-01-07 Power transmission wire shape finding method based on horizontal tension

Publications (2)

Publication Number Publication Date
CN111241731A true CN111241731A (en) 2020-06-05
CN111241731B CN111241731B (en) 2023-08-18

Family

ID=70872454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010014812.7A Active CN111241731B (en) 2020-01-07 2020-01-07 Power transmission wire shape finding method based on horizontal tension

Country Status (1)

Country Link
CN (1) CN111241731B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112036061A (en) * 2020-08-07 2020-12-04 国网浙江省电力有限公司电力科学研究院 Finite element modeling and dynamic response analysis method for omnibearing long and short leg tower line system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5617326A (en) * 1993-03-15 1997-04-01 Kabushiki Kaisha Toshiba Electronic circuit analyzing method with automatic adjustment of feedback loop effects
CN104166776A (en) * 2014-08-29 2014-11-26 南京理工大学 Transmission line conductor form-finding method based on ANSYS
CN106934147A (en) * 2017-03-09 2017-07-07 北京理工大学 A kind of geometry such as space truss rope net antenna look for shape method
CN109446611A (en) * 2018-10-17 2019-03-08 辽宁工程技术大学 A kind of close coupling tree looks for shape optimum design method
CN110298101A (en) * 2019-06-24 2019-10-01 国网浙江省电力有限公司电力科学研究院 A kind of transmission line of electricity wind-excited responese finite element method coupling wind system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5617326A (en) * 1993-03-15 1997-04-01 Kabushiki Kaisha Toshiba Electronic circuit analyzing method with automatic adjustment of feedback loop effects
CN104166776A (en) * 2014-08-29 2014-11-26 南京理工大学 Transmission line conductor form-finding method based on ANSYS
CN106934147A (en) * 2017-03-09 2017-07-07 北京理工大学 A kind of geometry such as space truss rope net antenna look for shape method
CN109446611A (en) * 2018-10-17 2019-03-08 辽宁工程技术大学 A kind of close coupling tree looks for shape optimum design method
CN110298101A (en) * 2019-06-24 2019-10-01 国网浙江省电力有限公司电力科学研究院 A kind of transmission line of electricity wind-excited responese finite element method coupling wind system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨风利;杨靖波;付东杰;李清华;: "输电线路导线舞动荷载分析" *
谢云云;金颖;黄琳雁;刘旭斐;王昊昊;郁琛;: "适用于有高差线路的覆冰输电导线ANSYS找形方法" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112036061A (en) * 2020-08-07 2020-12-04 国网浙江省电力有限公司电力科学研究院 Finite element modeling and dynamic response analysis method for omnibearing long and short leg tower line system
CN112036061B (en) * 2020-08-07 2024-05-10 国网浙江省电力有限公司电力科学研究院 Finite element modeling and dynamic response analysis method for omnibearing long and short leg tower line system

Also Published As

Publication number Publication date
CN111241731B (en) 2023-08-18

Similar Documents

Publication Publication Date Title
CN108959742B (en) Large-span power transmission tower line system aeroelastic model design method
CN105740548B (en) Method for calculating wind vibration of power transmission line under random wind load
CN105335569B (en) The multi-joint V-type insulator chain mechanical characteristic simulation method of UHV transmission line
CN112784347B (en) Cable-stayed bridge cable force reliability evaluation method based on bridge tower deformation and considering cable breakage
CN111651805B (en) Maximum wind deflection angle and wind vibration coefficient calculation method of suspension insulator string by considering linear shape and linear length influence factors
CN109063393B (en) Method for evaluating stress risk of tangent tower
WO2023168647A1 (en) Pseudo-static force calculation method for power transmission line lateral imbalance force
CN104166776B (en) A kind of transmission line wire based on ANSYS looks for shape method
CN111506946A (en) Calculation method for designing wind load of conventional power transmission tower and line based on effective load method and tower-line separation method and considering tower-line coupling influence
CN109918854A (en) A kind of composite insulation cross arm endurance quality appraisal procedure and device
CN111241731A (en) Transmission conductor shape finding method based on horizontal tension
CN111985018A (en) Calculation method for designing wind load of ultrahigh large-span tower and line based on inertia force method and tower line separation method and considering tower line coupling influence
CN111506944A (en) Tower wind load pulsation reduction coefficient calculation method considering tower line coupling effect
CN110363677A (en) A kind of icing transmission line of electricity looks for shape method
CN111651920A (en) Transmission analysis method for wind-induced mechanical energy of long-span power transmission tower under tower-line system
CN105631169B (en) A kind of cable-rod structure initial strain determines method
CN109784780B (en) Method, device and equipment for evaluating toughness of power system
CN112115616A (en) Dynamic characteristic analysis method and device of power transmission tower
CN111313419A (en) Power distribution network elastic lifting method and system based on dynamic reconfiguration in extreme weather
CN113313402B (en) Line state determining method, device, equipment and medium
CN111507027B (en) Method for judging integral power failure time of steel truss tower structure
CN106126862A (en) Transmission tower finite element modeling method based on GID and ABAQUS
CN109190931B (en) Three-dimensional space coordinate-based hard jumper iterative calculation method
CN113656857A (en) Method for automatically generating stringing construction scheme
CN114172069B (en) Method for determining three-dimensional hard jumper scheme of tension tower

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