CN113177343A - Power transmission line sag calculation method, equipment and computer medium - Google Patents
Power transmission line sag calculation method, equipment and computer medium Download PDFInfo
- Publication number
- CN113177343A CN113177343A CN202110581297.5A CN202110581297A CN113177343A CN 113177343 A CN113177343 A CN 113177343A CN 202110581297 A CN202110581297 A CN 202110581297A CN 113177343 A CN113177343 A CN 113177343A
- Authority
- CN
- China
- Prior art keywords
- transmission line
- wire
- power transmission
- sag
- model
- 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.)
- Pending
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 65
- 238000004364 calculation method Methods 0.000 title claims description 21
- 239000004020 conductor Substances 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000012805 post-processing Methods 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 5
- 238000002940 Newton-Raphson method Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 9
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000004590 computer program Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a method, equipment and computer medium for calculating sag of a power transmission line, wherein the method comprises the following steps: establishing a transmission line conductor initial position model according to a catenary equation; carrying out shape finding on the initial position model of the transmission line conductor, and determining a conductor finite element model; applying a temperature condition to the wire finite element model to obtain a first sag value of the wire; and applying different temperature boundary conditions to different section units in the axial direction of the wire, and extracting a second sag value of the wire through post-processing. The method is based on a finite element simulation model of a catenary formula method, calculates the line sag change condition under the factors of temperature and stress matching and the like, applies an uneven temperature distribution field near the power transmission line under the condition of the forest fire to the finite element model of the wire, calculates the sag change when the temperature unevenness occurs on the large-span power transmission line wire along the length direction of the large-span power transmission line wire, and can more accurately predict the forest fire tripping risk according to the solved sag change information.
Description
Technical Field
The invention relates to the technical field of data simulation, in particular to a method and equipment for calculating sag of a power transmission line and a computer medium.
Background
Due to the influence of natural and human factors such as extreme weather, burning wasteland, sacrifice and the like, large-scale mountain fire disasters nearby the corridor of the overhead transmission line sometimes occur. When a mountain fire occurs below the overhead transmission line, the gap insulation strength of the transmission line can be obviously reduced, the transmission line is easy to trip and stop, and the safety operation of a power grid is seriously threatened. However, the existing mountain fire trip risk prediction method does not consider the effect of high-temperature factors on the sag of the power transmission line, and the accuracy of mountain fire trip risk assessment is influenced.
At present, the traditional catenary method is mostly selected to calculate the sag length of the power transmission line at home and abroad, but the method does not consider the influence of the temperature of a wire transition point on the sag and the influence of the forest fire on the temperature of the wire under the condition of the forest fire, the influence of the temperature change on the sag of the wire cannot be accurately calculated, and the accurate forest fire tripping risk prediction result is difficult to obtain.
Disclosure of Invention
The invention provides a method, equipment and computer medium for calculating sag of a power transmission line.
In a first aspect, the present invention provides a method for calculating sag of a power transmission line, including:
establishing a transmission line conductor initial position model according to a catenary equation;
carrying out shape finding on the initial position model of the transmission line conductor, and determining a conductor finite element model;
applying a temperature condition to the wire finite element model to obtain a first sag value of the wire;
and applying different temperature boundary conditions to different section units in the axial direction of the wire, and extracting a second sag value of the wire through post-processing.
The power transmission line sag calculation method provided by the invention is based on a finite element simulation model of a catenary formula method, calculates the line sag change condition under the factors of temperature and stress matching and the like, applies an uneven temperature distribution field near the power transmission line under the condition of the forest fire to the finite element model of the wire, calculates the sag change when the temperature unevenness occurs on the wire of the large-span power transmission line along the length direction of the wire, and can more accurately predict the forest fire tripping risk according to the solved sag change information.
Optionally, before the establishing the initial position model of the transmission line conductor according to the catenary equation, the method further includes: selecting a unit according to the type of the lead; and defining material parameters according to the model of the lead.
Optionally, the unit is provided as a LINK10 unit; the material parameters comprise the diameter of the wire, the sectional area of the wire, the elastic coefficient, the linear expansion coefficient, the mass per unit length, the calculated tension and the span of the power transmission line.
Optionally, the catenary equation is:
wherein σ0The axial stress of the lowest point of the sag is expressed in MPa, and the gamma represents the dead weight specific load of the power transmission line and is expressed in MPa/m.
Optionally, the shape finding is performed on the initial position model of the transmission line conductor, specifically: applying boundary conditions and initial strain to the electric transmission line conductor initial position model; and setting large deformation and stress rigidization and carrying out nonlinear solution to obtain a finite element model of the wire, wherein the finite element model of the wire meets the geometric shape and theoretical error conditions of the wire.
Optionally, the nonlinear solving includes: and (4) performing iterative convergence calculation by adopting a Newton-Raphson method.
Optionally, the shape finding of the model of the initial position of the conductor of the power transmission line further includes: applying initial strain and a gravity field to the lead, wherein the calculation equation of the initial strain is as follows:
in the formula,
indicating initial response of probeChange, sigma0The axial stress of the lowest point of the sag is expressed in Mpa, gamma represents the dead weight specific load of the power transmission line and is expressed in Mpa/m,
representing the initial strain of the wire in the horizontal direction after stringing is completed.
Optionally, the
The method specifically comprises the following steps:
wherein P is the tension of the wire, A is the sectional area of the wire, and E is the Young's modulus of the wire.
The process of shape finding of the lead by applying initial strain and gravity can ensure the accuracy of shape finding calculation, and compared with a theoretical value, the simulation value considers the nonlinear condition of the structure and can more accurately simulate the actual condition of the engineering.
Meanwhile, the shape finding is carried out on the lead by directly establishing a catenary equation of the lead, and the elastic modulus of the lead does not need to be reset after the shape finding, so that the shape finding speed is higher.
In a second aspect, the present invention further provides a data processing apparatus, including a processor, coupled with a memory, where the memory stores a program, and the program is executed by the processor, so that the data processing apparatus executes the power transmission line sag calculation method according to the first aspect.
In a third aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the method for calculating the sag of the power transmission line according to the first aspect.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for calculating sag of a power transmission line according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be understood that the step numbers used herein are for convenience of description only and are not intended as limitations on the order in which the steps are performed.
It is to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the specification of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of the described features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The term "and/or" refers to and includes any and all possible combinations of one or more of the associated listed items.
As shown in fig. 1, in a first aspect, an embodiment of the present invention provides a method for calculating a sag of a power transmission line, including the following steps.
S11: and establishing a transmission line conductor initial position model according to a catenary equation.
In the embodiment, before the initial position model of the transmission line conductor is established according to the catenary equation, unit selection is carried out according to the type of the conductor, and then material parameters are defined according to the type of the conductor.
Specifically, a proper unit can be selected from finite element software ANSYS, and a power transmission wire model is established according to the wire model definition material parameters.
Wherein, the wire model definition material parameters mainly include: the diameter of the wire, the sectional area of the wire, the elastic coefficient, the linear expansion coefficient, the mass per unit length, the calculated tension, the span of the power transmission line and other parameters; and the wire of the power transmission line is selected to be subjected to simulation by a LINK10 unit.
In this embodiment, the geometric model of the initial position of the wire is established according to the catenary equation of the wire, which is:
wherein σ0The axial stress of the lowest point of the sag is expressed in MPa, and the gamma represents the dead weight specific load of the power transmission line and is expressed in MPa/m.
S12: and carrying out shape finding on the initial position model of the transmission line conductor, and determining a conductor finite element model.
In this embodiment, boundary conditions and initial strain are applied to the initial position model of the transmission line conductor, large deformation and stress rigidization are set, and nonlinear solution is performed to obtain a conductor finite element model satisfying the geometric shape and theoretical error conditions of the conductor.
Specifically, firstly, a boundary condition and an initial strain are applied to an initial position model of a transmission line lead, and an initial strain and a gravity field are applied to the lead, wherein a calculation equation of the initial strain is as follows:
in the formula,
representing the tentative initial strain, σ0The axial stress of the lowest point of the sag is expressed in Mpa, gamma represents the dead weight specific load of the power transmission line and is expressed in Mpa/m,
representing the initial strain of the wire in the horizontal direction after stringing is completed.
Wherein,
in particular to
P is the tension of the wire, A is the cross-sectional area of the wire, and E is the Young's modulus of the wire.
And then setting large deformation and stress rigidization effects, performing iterative convergence calculation by adopting a Newton-Raphson method, outputting the displacement and the internal force of the wire, adding the obtained wire displacement value into the initial position model of the wire of the power transmission line, updating the node coordinates of the model until the geometric shape and the theoretical error value of the wire in the model meet the preset requirements, and finishing shape finding.
The invention carries out shape finding on the lead by applying initial strain and gravity, compares the arc sag analog value after shape finding with a theoretical value, and sets the comparison error within an allowable range so as to ensure the accuracy of shape finding calculation. Meanwhile, the shape finding is carried out on the lead by directly establishing a catenary equation of the lead, and the elastic modulus of the lead does not need to be reset after the shape finding, so that the shape finding speed is higher.
S13: and applying a temperature condition to the finite element model of the wire to obtain a first sag value of the wire.
Specifically, in the finite element model of the conductor, a temperature condition is applied to the whole conductor to simulate the sag change of the transmission conductor when the conductor is uniformly heated, and a first sag value of the conductor is obtained by post-processing extraction, namely: sag values when the wire temperature is uniformly distributed.
S14: and applying different temperature boundary conditions to different section units in the axial direction of the wire, and extracting a second sag value of the wire through post-processing.
Specifically, an uneven temperature distribution field is applied to the initial model of the wire, and different temperature boundary conditions are applied to units on different sections of the wire to obtain a second sag value of the wire, so that the effect of uneven temperature distribution of the wire under the high-temperature condition of the mountain fire is better simulated.
The method provided by the invention considers the influence caused by the rise of the operating temperature of the wire under the condition of the forest fire, calculates the sag value in a sectional mode when calculating the temperature field, considers the temperature of the wire transfer point of the power transmission wire, and takes the condition that the sag of the power transmission line suddenly rises due to the sudden stress change into account.
In a second aspect, an embodiment of the present invention further provides a data processing device, which includes a processor, where the processor is coupled to a memory, and the memory stores a program, where the program is executed by the processor, so that the data processing device executes the power transmission line sag calculation method according to any one of the above embodiments.
In a third aspect, an embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the method for calculating the sag of the power transmission line according to any of the above embodiments is implemented.
It will be understood by those skilled in the art that all or part of the processes of the methods of the above embodiments may be implemented by a computer program, which may be stored in a computer-readable storage medium, and may include the processes of the embodiments of the methods when executed. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (10)
1. A sag calculation method for a power transmission line is characterized by comprising the following steps:
establishing a transmission line conductor initial position model according to a catenary equation;
carrying out shape finding on the initial position model of the transmission line conductor, and determining a conductor finite element model;
applying a temperature condition to the wire finite element model to obtain a first sag value of the wire;
and applying different temperature boundary conditions to different section units in the axial direction of the wire, and extracting a second sag value of the wire through post-processing.
2. The sag calculation method for the power transmission line according to claim 1, wherein before the establishing of the initial position model of the power transmission line conductor according to the catenary equation, the sag calculation method further comprises:
selecting a unit according to the type of the lead;
and defining material parameters according to the model of the lead.
3. The power transmission line sag calculation method according to claim 2, further comprising:
the cell is set as a LINK10 cell;
the material parameters comprise the diameter of the wire, the sectional area of the wire, the elastic coefficient, the linear expansion coefficient, the mass per unit length, the calculated tension and the span of the power transmission line.
4. The method for calculating the sag of the power transmission line according to claim 1, further comprising:
the catenary equation is:
wherein σ0The axial stress of the lowest point of the sag is expressed in MPa, and the gamma represents the dead weight specific load of the power transmission line and is expressed in MPa/m.
5. The method for calculating the sag of the power transmission line according to claim 1, wherein the shaping of the initial position model of the wire of the power transmission line comprises:
applying boundary conditions and initial strain to the electric transmission line conductor initial position model;
and setting large deformation and stress rigidization and carrying out nonlinear solution to obtain a finite element model of the wire, wherein the finite element model of the wire meets the geometric shape and theoretical error conditions of the wire.
6. The power transmission line sag calculation method according to claim 5, wherein the nonlinear solving comprises:
and (4) performing iterative convergence calculation by adopting a Newton-Raphson method.
7. The power transmission line sag calculation method according to claim 5, further comprising:
applying initial strain and a gravity field to the lead, wherein the calculation equation of the initial strain is as follows:
in the formula,
representing the tentative initial strain, σ0The axial stress of the lowest point of the sag is expressed in Mpa, gamma represents the dead weight specific load of the power transmission line and is expressed in Mpa/m,
representing the initial strain of the wire in the horizontal direction after stringing is completed.
8. The method according to claim 7, wherein the sag calculation is performed by using a linear interpolation method
The method specifically comprises the following steps:
wherein P is the tension of the wire, A is the sectional area of the wire, and E is the Young's modulus of the wire.
9. A data processing apparatus, characterized by comprising:
a processor coupled to a memory, the memory storing a program that is executed by the processor to cause the data processing apparatus to perform the power transmission line sag calculation method of any one of claims 1 to 8.
10. A computer storage medium storing computer instructions for performing the method of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110581297.5A CN113177343A (en) | 2021-05-26 | 2021-05-26 | Power transmission line sag calculation method, equipment and computer medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110581297.5A CN113177343A (en) | 2021-05-26 | 2021-05-26 | Power transmission line sag calculation method, equipment and computer medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113177343A true CN113177343A (en) | 2021-07-27 |
Family
ID=76927571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110581297.5A Pending CN113177343A (en) | 2021-05-26 | 2021-05-26 | Power transmission line sag calculation method, equipment and computer medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113177343A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113983983A (en) * | 2021-10-27 | 2022-01-28 | 深圳飞赛精密钣金技术有限公司 | Wire temperature and sag actual measurement method and system |
| CN114840958A (en) * | 2022-06-30 | 2022-08-02 | 国网湖北省电力有限公司经济技术研究院 | Calculation method for three-dimensional space attitude of construction block |
| CN115357965A (en) * | 2022-06-27 | 2022-11-18 | 中建三局第一建设工程有限责任公司 | Self-anchored suspension bridge and method for determining bridge forming line shape thereof |
| CN117932898A (en) * | 2023-12-28 | 2024-04-26 | 国家电网有限公司华东分部 | Overhead line induced electricity determining method and device, storage medium and computer equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108021785A (en) * | 2017-12-15 | 2018-05-11 | 中国电建集团河南省电力勘测设计院有限公司 | One kind string weight Tensile Sag algorithm |
| CN108955609A (en) * | 2018-05-28 | 2018-12-07 | 广东电网有限责任公司 | Aerial cable sag computing method and device |
-
2021
- 2021-05-26 CN CN202110581297.5A patent/CN113177343A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108021785A (en) * | 2017-12-15 | 2018-05-11 | 中国电建集团河南省电力勘测设计院有限公司 | One kind string weight Tensile Sag algorithm |
| CN108955609A (en) * | 2018-05-28 | 2018-12-07 | 广东电网有限责任公司 | Aerial cable sag computing method and device |
Non-Patent Citations (1)
| Title |
|---|
| 张世龙: "输电导线增容过程中温升和力学响应的有限元分析", 万方数据库, 30 July 2015 (2015-07-30), pages 37 - 48 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113983983A (en) * | 2021-10-27 | 2022-01-28 | 深圳飞赛精密钣金技术有限公司 | Wire temperature and sag actual measurement method and system |
| CN115357965A (en) * | 2022-06-27 | 2022-11-18 | 中建三局第一建设工程有限责任公司 | Self-anchored suspension bridge and method for determining bridge forming line shape thereof |
| CN115357965B (en) * | 2022-06-27 | 2023-10-27 | 中建三局第一建设工程有限责任公司 | Self-anchored suspension bridge and bridge forming line shape determining method thereof |
| CN114840958A (en) * | 2022-06-30 | 2022-08-02 | 国网湖北省电力有限公司经济技术研究院 | Calculation method for three-dimensional space attitude of construction block |
| CN117932898A (en) * | 2023-12-28 | 2024-04-26 | 国家电网有限公司华东分部 | Overhead line induced electricity determining method and device, storage medium and computer equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113177343A (en) | Power transmission line sag calculation method, equipment and computer medium | |
| CN106096105B (en) | Power transmission circuit caused by windage transient response calculation method | |
| CN105004504B (en) | The starting of oscillation wind speed appraisal procedure and selection method of composite insulator | |
| CN110160492B (en) | Method and device for monitoring inclination of power transmission tower | |
| CN105740549A (en) | Wind vibration simulation method for 220kV power transmission line tower wire coupling system | |
| He et al. | A method for analyzing stability of tower-line system under strong winds | |
| CN106651140B (en) | Module difference evaluation method and device for power transmission line risk in typhoon area | |
| CN104535233B (en) | A kind of icing transmission line of electricity stress monitoring system | |
| CN106934096B (en) | Method for solving steel core temperature based on surface temperature of overhead conductor | |
| CN107358020B (en) | Method for calculating line tightening length of distribution network wire with overlarge sag | |
| CN106482848A (en) | One kind is based on the generalized inverse three-core cable conductor temperature dynamic acquisition method of M P | |
| CN110781589B (en) | Method for detecting overheating fault of lap joint of gas insulated metal closed power transmission line | |
| CN105205202A (en) | Current carrying capacity calculation method | |
| CN113128080A (en) | Method and apparatus for predicting temperature of wire harness | |
| CN116222466A (en) | Method, device, equipment and medium for monitoring icing thickness of tension tower wire | |
| CN110224335A (en) | A kind of aerial condutor erection method and system considering creep compliance | |
| CN115224692A (en) | Probabilistic load flow calculation method and system considering electro-thermal mechanics dynamics of overhead conductor | |
| KR101658748B1 (en) | Reliability design assistance device and reliability design assistance method | |
| CN114880828A (en) | Cable quality scoring method and system | |
| CN114172069B (en) | Method for determining three-dimensional hard jumper scheme of tension tower | |
| CN117932898B (en) | Overhead line induced electricity determining method and device, storage medium and computer equipment | |
| CN117313409B (en) | Arc mathematical model parameter optimization method, system, terminal and medium | |
| CN118378461B (en) | Support structure checking simulation method and equipment based on mechanical numerical simulation | |
| CN112380669A (en) | Method for quickly acquiring maximum temperature of steady-state core of four-loop groove cable | |
| CN118114409A (en) | Static equilibrium state analysis method and system for continuous multi-span transmission line |
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 |