CN108133089B - Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line - Google Patents
Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line Download PDFInfo
- Publication number
- CN108133089B CN108133089B CN201711310584.2A CN201711310584A CN108133089B CN 108133089 B CN108133089 B CN 108133089B CN 201711310584 A CN201711310584 A CN 201711310584A CN 108133089 B CN108133089 B CN 108133089B
- Authority
- CN
- China
- Prior art keywords
- shielding ring
- support rod
- connecting plate
- lead
- damper
- 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.)
- Active
Links
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]
Abstract
A load application method for finite element numerical calculation of an electrostatic field of a shielding ring of a power transmission line is characterized in that a geometric model is established according to the structural dimensions of the shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, a tower cross arm and a damper; establishing two air domain models, namely an air domain I model and an air domain II model; defining dielectric constants of a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, a tower cross arm and a damper material, and defining the dielectric constant of an air domain; dividing a calculation domain into finite element meshes with tetrahedral structures; carrying out effective grid division; and (4) applying a load. The invention does not divide the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring and the damper, only divides the air around the shielding ring into grids, applies electric potential on the surfaces of the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring and the damper, and then carries out solution calculation, thereby conveniently obtaining the electric field strength value of the surface of the shielding ring, and the maximum value appears on the surface of the shielding ring.
Description
Technical Field
The invention relates to a finite element numerical calculation method of an electrostatic field, in particular to a load application method for the finite element numerical calculation of the electrostatic field of a shielding ring of a power transmission line.
Background
Finite element numerical calculation is to disperse a continuous body structure into finite unit bodies which are hinged at nodes, load loads on the nodes and calculate the electric field of each node under the action of the loads. The original continuum solution is replaced by an approximation of the solution for the discrete body, which is close to the true value when the cells are divided sufficiently densely.
At present, the finite element numerical calculation method of the electrostatic field is mature, and the load application method is to divide the whole shielding ring into grids, apply electric potential to each node, and then perform solution calculation. However, since the shielding ring made of the metal material is an equipotential body, the internal electric field intensity is 0, after the shielding ring is wholly divided, the surface grid nodes belong to a part of the calculation domain made of the metal material, and the electric field intensity of the surface nodes is also 0, so that the actual electric field intensity value cannot be observed.
Disclosure of Invention
The invention provides a load applying method for electric transmission line shielding ring electrostatic field finite element numerical calculation, which does not need to carry out mesh subdivision on a shielding ring, only carries out mesh division on air around the shielding ring, thus, nodes on the surface of the shielding ring belong to an air calculation domain, the calculated node electric field strength is the shielding ring surface field strength, and the electric field strength value on the surface of the shielding ring can be conveniently obtained.
The technical scheme adopted by the invention is as follows:
a load applying method for finite element numerical calculation of an electrostatic field of a shielding ring of a power transmission line comprises the following steps:
step 1: establishing a geometric model for the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper by utilizing finite element software according to the structural dimensions of the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper;
step 2: establishing two air domain models, wherein the air domain model I is 5% larger than the shielding ring model, and the size of the air domain model II is three times of that of the model of the overall calculation region;
and step 3: defining dielectric constants of a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, a tower cross arm and a damper material, and defining the dielectric constant of an air domain;
and 4, step 4: dividing a calculation domain into finite element meshes with tetrahedral structures;
and 5: effective grid division is carried out, a reasonable unit shape is ensured, and the node shape is a three-dimensional node electrostatic entity unit; establishing the minimum value of the unit size during the whole grid division, carrying out the unit grid division on the established cross arm and the two air domain models, and enabling all nodes on the divided unit grid bodies to be equal and to have the same value; and the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring and the damper are not subjected to grid division, so that the metal surfaces of the shielding ring, the shielding ring support rod, the lead and the shielding ring are at the same potential.
Step 6: applying load, selecting a cross arm entity, applying zero potential on the cross arm, and loading the zero potential on an air boundary; selecting a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring and the outermost layer of the metal surface of the damper, wherein the potential of the outermost layer is high potential, and loading U/H
And kV voltage.
The invention relates to a load applying method for electrostatic field finite element numerical calculation of a transmission line shielding ring, which is compared with the applying method for electrostatic field finite element numerical calculation load in the prior art:
the method of the invention is provided with: the shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper are not divided, only the air around the shielding ring is subjected to grid division, electric potentials are applied to the surfaces of the shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper, then solution calculation is carried out, the electric field strength value of the surface of the shielding ring can be obtained conveniently, and the maximum value appears on the surface of the shielding ring.
Drawings
Fig. 1 is a diagram showing an electric field intensity distribution of a surface of a shield ring in the prior art.
FIG. 2 is a cross arm geometric model diagram.
FIG. 3 is a geometric model diagram of a shield ring body, a shield ring support rod, a connecting plate of a lead and a shield ring, and an anti-vibration hammer;
the vibration damper comprises a cross arm 1, a shielding ring 2, a shielding ring support rod 3, a connecting plate of a wire 4 and a shielding ring and a vibration damper 5.
FIG. 4 is a diagram of an air space model;
wherein, 6-air domain one face.
FIG. 5 is a diagram of an air domain two model;
among them, 7-the second plane of the air domain.
Fig. 6 is a diagram of a split air domain one and air domain two model.
Wherein, 8-the subdivision grid of air area one, 9-the subdivision grid of air area two.
Fig. 7 is a model diagram of a shielding ring body, a shielding ring support rod, a connecting plate of a lead and a shielding ring, and a damper.
FIG. 8 is a distribution diagram of electric field intensity on the surface of the shielding ring according to the present invention.
Detailed Description
A load applying method for finite element numerical calculation of an electrostatic field of a shielding ring of a power transmission line comprises the following steps:
step 1: according to the structural dimensions of the shielding ring, the shielding ring support rod, the connecting plate of the conducting wire and the shielding ring, the tower cross arm and the damper, a geometric model is established for the shielding ring, the shielding ring support rod, the connecting plate of the conducting wire and the shielding ring, the tower cross arm and the damper by using finite element software ANSYS 12.0.
Step 2: two air domain models are established, the air domain one model is 5% larger than the shielding ring model, according to the calculation experience, the splitting and the calculation are convenient, for example, the size of the shielding ring model is 1 meter in length, 1 meter in width and 1 meter in height, and the size of the air domain one model can be set to be 1.05 meter in length, 1.05 meter in width and 1.05 meter in height. The air domain two model size is three times the model of the overall calculation region. Strictly speaking theoretically, the electric field intensity at the infinite boundary is zero, but if a very large air domain model is established, a great calculation amount is caused, and according to the calculation experience, when the air domain model is 3 times of the whole calculation region model, the boundary electric field intensity is attenuated to zero, which can be approximately considered as zero.
And step 3: defining the relative dielectric constants of a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, a tower cross arm and a damper material, wherein the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper material are made of metal, and the dielectric constant of the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper material relative to air is set to be 1; the relative dielectric constant defining the air domain is 1.
And 4, step 4: the computational domain is divided into finite element meshes of tetrahedral structure.
And 5: in order to ensure the calculation precision and improve the calculation efficiency, effective grid division is carried out, and a reasonable unit shape is ensured, wherein the node shape is a three-dimensional node electrostatic entity unit; establishing the minimum value of the unit size during the whole grid division, carrying out the unit grid division on the established cross arm and the two air domain models, and enabling all nodes on the divided unit grid bodies to be equal and to have the same value; and the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring and the damper are not subjected to grid division, so that the metal surfaces of the shielding ring, the shielding ring support rod, the lead and the shielding ring are at the same potential.
Step 6: applying load, selecting a cross arm entity, applying zero potential on the cross arm, and loading the zero potential on an air boundary; selecting a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, and the outermost layer of the metal surface of the damper, wherein the potential of the outermost layer is high potential, loading load on the high potential for calculating the electric field intensity, and the numerical value of the outermost layer is equal to the line voltage grade divided by the line voltage grade
For example: for a 1000kV line, a high potential is applied 1000-
I.e. a potential of 577.37 kV.
The shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper are subjected to grid division, potentials are applied to all nodes, and then the electric field intensity on the surface of the shielding ring obtained after solution calculation is shown as 1.
In the method for applying the finite element numerical calculation load of the electrostatic field in the prior art, a shielding ring support rod, a connecting plate of a lead and the shielding ring and a damper are subjected to grid division, a potential is applied to each node, then solution calculation is carried out, and the electric field intensity on the surface of the shielding ring is shown in figure 1. From FIG. 1, the values of the electric field strength of the surface of the shielding ring were observed, and the length in the finite element software was m, and the voltage in the finite element software was kV, so that the maximum value of the surface field strength was 0.106 × 10-14kV/m, almost close to zero, and if the electric field intensity value of the surface of the shielding ring is desired, it is inconvenient to take the value from the air adjacent to the shielding ring.
The method of the invention is provided with: the shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper are not subdivided, only the air around the shielding ring is subjected to grid division, electric potentials are applied to the surfaces of the shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper, then solution calculation is carried out, and the strength of the surface of the shielding ring is shown in a graph 8. From the observation of the electric field intensity on the surface of the shield ring in fig. 8, it is convenient to obtain the electric field intensity on the surface of the shield ring, the maximum value of which appears on the surface of the shield ring, and the maximum electric field intensity on the surface of the shield ring is 2816.9 kV/m.
Claims (1)
1. A load application method for finite element numerical calculation of an electrostatic field of a shielding ring of a power transmission line is characterized by comprising the following steps:
step 1: establishing a geometric model for the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper by utilizing finite element software according to the structural dimensions of the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper;
step 2: establishing two air domain models, wherein the air domain model I is 5% larger than the shielding ring model, and the size of the air domain model II is three times of that of the model of the overall calculation region;
and step 3: defining dielectric constants of a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, a tower cross arm and a damper material, and defining the dielectric constant of an air domain, wherein the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper material are made of metal, and the dielectric constant of the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring, the tower cross arm and the damper material relative to air is set to be 1; defining a relative dielectric constant of 1 for the air domain;
and 4, step 4: dividing a calculation domain into finite element meshes with tetrahedral structures;
and 5: effective grid division is carried out, a reasonable unit shape is ensured, and the node shape is a three-dimensional node electrostatic entity unit; establishing the minimum value of the unit size during the whole grid division, carrying out the unit grid division on the established cross arm and the two air domain models, and enabling all nodes on the divided unit grid bodies to be equal and to have the same value; the shielding ring, the shielding ring support rod, the connecting plate of the lead and the shielding ring and the damper are not subjected to grid division, so that the metal surfaces of the shielding ring, the shielding ring support rod, the lead and the shielding ring are at the same potential;
step 6: applying load, selecting a cross arm entity, applying zero potential on the cross arm, and loading the zero potential on an air boundary; selecting a shielding ring, a shielding ring support rod, a connecting plate of a lead and the shielding ring, and the outermost layer of the metal surface of the damper, wherein the potential of the outermost layer is high, and loading
A voltage;
and 7: and carrying out grid division on the shielding ring, the shielding ring support rod, the connecting plate of the wire and the shielding ring and the damper, applying a potential on each node, and then carrying out solution calculation to obtain the electric field intensity on the surface of the shielding ring.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711310584.2A CN108133089B (en) | 2017-12-11 | 2017-12-11 | Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711310584.2A CN108133089B (en) | 2017-12-11 | 2017-12-11 | Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108133089A CN108133089A (en) | 2018-06-08 |
| CN108133089B true CN108133089B (en) | 2021-07-20 |
Family
ID=62389206
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711310584.2A Active CN108133089B (en) | 2017-12-11 | 2017-12-11 | Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108133089B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108694296B (en) * | 2018-06-22 | 2021-09-24 | 武汉大学 | Finite element modeling method for electric field numerical calculation of large complex model |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8489372B2 (en) * | 2010-03-08 | 2013-07-16 | Livermore Software Technology Corp. | FEM-BEM coupling methods and systems for sliding contact interface |
| CN102201669A (en) * | 2010-04-30 | 2011-09-28 | 北京网联直流工程技术有限公司 | Method for selecting direct-current hardware fitting for +/- 660 kV valve hall |
| CN102368281B (en) * | 2011-11-14 | 2013-12-18 | 浙江中科电声研发中心 | Numerical value simulating method for loudspeaker magnetic circuit system |
| CN104217061B (en) * | 2014-07-08 | 2017-04-12 | 国家电网公司 | Temperature field simulation design method for low-voltage distribution cabinet |
| CN106446418A (en) * | 2016-09-27 | 2017-02-22 | 广东电网有限责任公司珠海供电局 | Cable terminal inflow water defect electric-field distribution analysis method and device |
| CN106528947A (en) * | 2016-10-17 | 2017-03-22 | 国网辽宁省电力有限公司电力科学研究院 | Method for obtaining potential distribution of 1,000kV metal oxide arrester |
| CN106919773A (en) * | 2017-04-25 | 2017-07-04 | 哈尔滨理工大学 | A kind of numerical computations APP for big stator winding conductor bars in electrical machines end model |
-
2017
- 2017-12-11 CN CN201711310584.2A patent/CN108133089B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108133089A (en) | 2018-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110287629B (en) | Computer device, equipment and method for determining structural parameters of high-voltage basin-type insulator | |
| Steven et al. | Evolutionary topology and shape design for general physical field problems | |
| Abdel-Salam | Electric fields | |
| CN108133089B (en) | Load application method for finite element numerical calculation of electrostatic field of shielding ring of power transmission line | |
| CN104899379A (en) | ANSYS-based crosslinking polyethylene insulated cable temperature field finite element calculation method | |
| CN110705172A (en) | Method for calculating distribution of instantaneous electric field of transformer under lightning impulse voltage | |
| CN114022641A (en) | Method for generating finite element grid of honeycomb sandwich panel with pit defects | |
| CN111985137B (en) | Multivariable cutting level set optimization method for porous structure topology optimization | |
| CN103984792A (en) | Simulation method for researching influence on potential and electric field distribution of rod insulator by suspended sand | |
| CN108694299B (en) | ICEM-CFD-based two-dimensional finite element neutronics steady-state calculation method | |
| CN112072687B (en) | Method for determining stray capacitance value of alternating current filter | |
| CN109117537A (en) | A kind of optimization method and device of high-voltage commutation transformer end square ring arrangement | |
| CN113378428A (en) | Electric field distribution estimation method for alternating current filter lightning arrester of extra-high voltage direct current converter station | |
| Desouky et al. | A new contribution in reducing electric field distribution within/around medium voltage underground cable terminations | |
| Manikandan et al. | Design and analysis of grounding systems for wind turbines using Finite element method | |
| Liu et al. | Finite element method and charge coupling method for distortion field of AC transmission line | |
| CN114722659A (en) | Electrode gap electric field non-uniform coefficient calculation method | |
| Kherd et al. | Generalized Biharmonic Wang-Ball Surface | |
| CN107045488A (en) | A kind of surface boundary member computational methods based on Bezier curved surface quadrangle parametric equations | |
| CN112668216A (en) | Method for predicting breakdown voltage of gap of live working of ground potential of power transmission line | |
| Vidya et al. | Gaussian process regression based model for prediction of discharge voltage of air gaps under positive polarity lightning impulse voltages | |
| Onyedikachi et al. | Modeling Performance Variation of Earthing Enhancement Materials in Different Soils and Conditions | |
| Chumak et al. | Estimation of the area of field emission of a carbon nanotube using modelling in COMSOL Multiphysics | |
| Li et al. | Research on optimal design of needle-shaped tower grounding device using thin-shell theory | |
| Subrahmanyam et al. | Image charge effect on metallic particle in single phase gas insulated bus duct (GIB) |
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 |