CN108536906A - A kind of bushing temperature field distribution method of the calculating containing connected structure - Google Patents
A kind of bushing temperature field distribution method of the calculating containing connected structure Download PDFInfo
- Publication number
- CN108536906A CN108536906A CN201810170594.9A CN201810170594A CN108536906A CN 108536906 A CN108536906 A CN 108536906A CN 201810170594 A CN201810170594 A CN 201810170594A CN 108536906 A CN108536906 A CN 108536906A
- Authority
- CN
- China
- Prior art keywords
- connected structure
- calculating
- casing
- field distribution
- watchband
- 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
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/14—Pipes
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of bushing temperature field distribution method of the calculating containing connected structure, including S1 to establish the casing finite element model containing connected structure, divides volume mesh and surface grids；S2 calculates the contact resistance in connected structure, and is finally inversed by the volume resistivity of connected structure；S3 calculates the local pyrexia power of the whole heating power and connected structure of guide rod；S4 calculates the distribution of casing bulk temperature, according to result of calculation come the initial value of resistivity in amendment step S2, and iterates, and repeats step S2 S4, until the difference of measurement point is less than preset value in adjacent resistivity result of calculation twice.The contact resistance of existing different heat transfer mode, electric current at the kelvin effect and central conductor leader connecting portion in guide rod, between watchband-shaped contactor and guide rod in the casing that the present invention considers, it can be to carry out overall thermal linear accelerator with the casing with different connection structures, effectively reduce design cycle and design cost of the sleeve design in the process to thermal design.
Description
Technical field
The present invention relates to technical field of electric power, and in particular to a kind of bushing temperature field distribution side of the calculating containing connected structure
Method.
Background technology
Extra-high voltage casing is the critical component on many power equipments, such as transformer, converter station valve hall, the converting plant valve Room
Deng, play in the power system support conductor draw, electric insulation the effects that, safe and stable operation has been related to entirely
The reliability of electric system.Extra-high voltage casing is sent out other than needing by electric, mechanical force effect by central conductor leader ohm
Temperature gradient distribution caused by the joule heat of heat and dielectric is also an important factor for influencing safe operation.
At present for the calculating of extra-high voltage casing heat distribution mainly by numerical simulation means, it is contemplated that guide rod current-carrying feelings
The Joule heat of electric current under condition, ambient temperature situation will be under guide rod current-carrying capacities under the heat transfer types such as convection current, heat transfer and heat radiation
Guide rod is equivalent to a conductor, the Temperature Distribution of extra-high voltage casing is exported to calculate Joule heat, considers contact electricity
The presence of resistance.Numerical simulation result under this condition is that the highest position of bushing temperature is both present in central conductor leader and capacitance
The interface of fuse, and close to the centre position of conducting rod, the hottest spot temperature of casing, which should be less than, causes its electrical characteristic and reason
Change the limit value of characteristic variations.
However in recent years, the casing failure that a lot of defects due at conducting rod docking cause has occurred both at home and abroad, from event
Find that the position of hyperthermia and superheating trace is all often inside pipe casing event from the point of view of the disassembling analysis result for hindering casing, at conducting rod docking
The starting point of barrier.Thus, it will be seen that only considering that ohm fever of conducting rod is incomplete.In Practical Project, due to
The presence of connection structure in casing, to the contact resistance resistance value and conducting rod phase between the watchband-shaped contactor and conducting rod that connect
When, it is not possible to ignore.
Invention content
It is an object of the invention to overcome above-mentioned the deficiencies in the prior art, a kind of set tube temperature of the calculating containing connected structure is provided
Field distribution method is spent, to provide foundation for the design of casing, it is ensured that security reliability and stability when it runs.
To achieve the above object, the technical scheme is that：
A kind of bushing temperature field distribution method of the calculating containing connected structure, the method includes：
S1, the casing finite element model containing connected structure is established, divides volume mesh and surface grids；
Contact resistance in S2, calculating connected structure, and it is finally inversed by the volume resistivity of connected structure；
S3, calculate guide rod whole heating power and connected structure local pyrexia power；
S4, the distribution of casing bulk temperature is calculated, according to result of calculation come the initial value of resistivity in amendment step S2, and carried out
It iterates, repeats step S2-S4, until the difference of measurement point is less than preset in adjacent resistivity result of calculation twice
Value then calculates and terminates, exports result.
In step s 2, by establish watchband-shaped contactor equivalent submodel and consider contact surface area, surface roughness,
The case where electrode material, sensing finger pressure, gets off to calculate the size of the resistance value of the contact resistance in connected structure, and by contact resistance
On the equivalent volume resistivity for causing corresponding equivalent submodel of resistance value.
In step s3, for the conducting rod of cylindrical rule, directly using the AC resistance meter considered under kelvin effect
Formula is calculated to calculate calorific value；For watchband-shaped contactor, Coupled Electromagnetic-Thermal calculating is carried out with COMSOL Multiphysics,
In the hope of under the equivalent submodel of watchband-shaped contactor, calorific value when power current flows through the cannula center conducting rod with watchband-shaped contactor
Size.
In step s 4, using 3 D electromagnetic-coupled flow-thermal method, the temperature point of casing is calculated with finite element software
Cloth in calculating process and considers heat transfer, heat radiation and the influence of thermal convection current.
In step sl, set tube model is imported using APDL, obtains the finite element analysis model calculated, and pass through
APDL language establishes mesh generation program, the surface grids that the volume mesh in computational domain is coupled with the needs calculated in heat distribution into
Row divides.
The equivalent submodel of watchband-shaped contactor is established using COMSOL Multiphysics.
Compared with prior art, the present invention advantage is：
In the casing that the present invention considers kelvin effect in guide rod of existing different heat transfer mode, electric current and in
Contact resistance at heart conducting rod connecting portion, between watchband-shaped contactor and guide rod can be with different connection structures
Casing carries out overall thermal linear accelerator, effectively reduces the design cycle during sleeve design to thermal design and is designed to
This, has wide applicability and economy.
Description of the drawings
Fig. 1 is the flow chart provided in an embodiment of the present invention for calculating the bushing temperature field distribution method containing connected structure；
Fig. 2 is that the 500kV changes of current become dry type sleeve structure schematic diagram；
Fig. 3 is the enlarged structure schematic diagram at A in Fig. 2；
Fig. 4 is watchband-shaped contactor submodel schematic diagram；
Fig. 5 is electric current and volume heating power distribution schematic diagram；
Fig. 6 is based on 3 D electromagnetic-heat-flow coupling computational methods flow diagram；
Specific implementation mode
Present disclosure is described in further details with reference to the accompanying drawings and detailed description.
Embodiment：
As shown in fig.1, bushing temperature field distribution method of the calculating containing connected structure provided in this embodiment have including
Following steps：
S1, the casing finite element model containing connected structure is established, divides volume mesh and surface grids
The set tube model containing connected structure is built with 3 d modeling software, as Figure 2-3, which includes but not limited to lead
Bar, epoxy core, SF6Gas 5, flange arrangement and sheath, transformer oil 1, sheath 2, end metal fitting 3, internal shield 6, connection
Necessity structure such as component 4 imports set tube model with APDL (ANSYS Parametric Design Language), can be used
With the finite element analysis model of calculating, and mesh generation program is established by APDL language, in computational domain volume mesh with
The surface grids that the needs that heat distribution calculates couple are divided.
Contact resistance in S2, calculating connected structure, and it is finally inversed by the volume resistivity of connected structure
Establish the equivalent submodel of watchband-shaped contactor using COMSOL Multiphysics, and consider contact surface area,
The case where surface roughness, electrode material, sensing finger pressure, gets off the contact resistance resistance value size calculated in connected structure, and will
Its is equivalent to cause corresponding equivalent submodel (preferably section is the ring model of " bow " type, and watchband-shaped contactor is as shown in 8 in Fig. 4)
Volume resistivity on, by under complex situations at contact surface with kelvin effect electric current generate Joule heat be equivalent to submodel
Kelvin effect electric current generate Joule heat so that contact resistance calorific value more tallies with the actual situation, and facilitates later-mentioned step
Simplify and calculates.The equivalent submodel of watchband-shaped contactor is as shown in figure 4, including current－carrying conductor 7, watchband-shaped contactor 8, SF6Gas 9 and
Guide rod inner air 10.
S3, calculate guide rod whole heating power and connected structure local pyrexia power
Under conditions of considering skin effect, current field has certain skin depth, it will influences institute during design calculates
Equivalent resistance, to influence the Joule heat size that electric current flows through conductor.This step, which is intended to calculate, considers collection skin under power frequency
Conductor temp.-elevating amount after effect.For the conducting rod of cylindrical rule, directly using the AC resistance meter considered under kelvin effect
Formula is calculated to calculate calorific value；For watchband-shaped contactor, Coupled Electromagnetic-Thermal calculating is carried out with COMSOL Multiphysics,
In the hope of under the model, only a small range existing for watchband-shaped contactor calculate power current flow through watchband-shaped contactor when calorific value it is big
It is small.Two parts calculate output result be unit volume heating power, as heat distribution calculate heat source respectively to difference
Position is loaded；FEM calculation domain is reduced, and improves load mode, the effectively save calculating time.Electric current and volume
Heating power distribution is as shown in Figure 5.
S4, the distribution of casing bulk temperature is calculated；
Using 3 D electromagnetic-coupled flow-thermal method, the Temperature Distribution of casing is calculated with FEM software, is calculated
Need in the process consider heat transfer, heat radiation and thermal convection current influence, be used in combination result of calculation to actual temperature off field resistivity into
Row is corrected, and is iterated and (repeat step S2-S4) so that the calculated value approaching to reality value of resistivity, until temperature
When error between the adjacent iteration twice of field measurement point is less than 0.01K, then it is assumed that the thermo parameters method is approximately under actual conditions
Bushing temperature field distribution, stop calculate output result.3 D electromagnetic-heat-flow coupling computational methods flow is as shown in Figure 6.
It follows that the set provided in this embodiment for calculating the bushing temperature field distribution method containing connected structure and considering
Existing different heat transfer mode, electric current at the kelvin effect and central conductor leader connecting portion in guide rod, watch band contact in pipe
Refer to the contact resistance between guide rod, can be to carry out overall thermal linear accelerator with the casing with different connection structures,
Design cycle and design cost of the sleeve design in the process to thermal design are effectively reduced, there is wide applicability and economy
Property.
Above-described embodiment simply to illustrate that the present invention technical concepts and features, it is in the art the purpose is to be to allow
Those of ordinary skill cans understand the content of the present invention and implement it accordingly, and it is not intended to limit the scope of the present invention.It is all
It is the equivalent changes or modifications made according to the essence of the content of present invention, should all covers within the scope of the present invention.
Claims (6)
1. a kind of bushing temperature field distribution method of the calculating containing connected structure, which is characterized in that the method includes：
S1, the casing finite element model containing connected structure is established, divides volume mesh and surface grids；
Contact resistance in S2, calculating connected structure, and it is finally inversed by the volume resistivity of connected structure；
S3, calculate guide rod whole heating power and connected structure local pyrexia power；
S4, the distribution of casing bulk temperature is calculated, resistivity initial value is obtained according to result of calculation in amendment step S2, and carry out repeatedly
Iteration repeats step S2-S4, until the difference of measurement point is less than preset value in adjacent resistivity result of calculation twice, then
Calculating terminates, and exports result.
2. calculating the bushing temperature field distribution method containing connected structure as described in claim 1, which is characterized in that in step S2
In, by establishing the equivalent submodel of watchband-shaped contactor and considering contact surface area, surface roughness, electrode material, fingertip pressure
The case where power, gets off to calculate the size of the resistance value of the contact resistance in connected structure, and the equivalent cause of contact resistance resistance value is corresponding
On the volume resistivity of equivalent submodel.
3. calculating the bushing temperature field distribution method containing connected structure as claimed in claim 2, which is characterized in that in step S3
In, for the conducting rod of cylindrical rule, directly using the AC resistance calculation formula under consideration kelvin effect to calorific value into
Row calculates；For watchband-shaped contactor, Coupled Electromagnetic-Thermal calculating is carried out with COMSOL Multiphysics, in the hope of watchband-shaped contactor
Calorific value size when power current flows through the cannula center conducting rod with watchband-shaped contactor under equivalent submodel.
4. calculating the bushing temperature field distribution method containing connected structure as described in claim 1, which is characterized in that in step S4
In, using 3 D electromagnetic-coupled flow-thermal method, with the Temperature Distribution of finite element software calculating casing, in calculating process and examine
Consider heat transfer, heat radiation and thermal convection current influence.
5. calculating the bushing temperature field distribution method containing connected structure as described in claim 1, which is characterized in that in step S1
In, set tube model is imported using APDL, obtains the finite element analysis model calculated, and grid is established by APDL language and is drawn
Branch, the surface grids coupled with the needs calculated in heat distribution to the volume mesh in computational domain divide.
6. calculating the bushing temperature field distribution method containing connected structure as claimed in claim 2, which is characterized in that use
COMSOL Multiphysics establish the equivalent submodel of watchband-shaped contactor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810170594.9A CN108536906A (en) | 2018-03-01 | 2018-03-01 | A kind of bushing temperature field distribution method of the calculating containing connected structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810170594.9A CN108536906A (en) | 2018-03-01 | 2018-03-01 | A kind of bushing temperature field distribution method of the calculating containing connected structure |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108536906A true CN108536906A (en) | 2018-09-14 |
Family
ID=63486479
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810170594.9A Pending CN108536906A (en) | 2018-03-01 | 2018-03-01 | A kind of bushing temperature field distribution method of the calculating containing connected structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108536906A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110619138A (en) * | 2019-07-08 | 2019-12-27 | 国家电网有限公司 | Hardware fitting electric field distribution calculation method and device based on electric heating coupling field model |
| CN111368385A (en) * | 2018-12-25 | 2020-07-03 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for operating characteristics of sleeve in electric-thermal coupling field |
| CN112179496A (en) * | 2020-09-15 | 2021-01-05 | 西安交通大学 | Infrared diagnosis operation and maintenance method for ultra-high voltage and extra-high voltage wall bushing with contact finger structure |
| CN112798643A (en) * | 2020-12-09 | 2021-05-14 | 国网宁夏电力有限公司检修公司 | Maintenance method, medium and system for contact finger of wall bushing |
| CN113128059A (en) * | 2021-04-23 | 2021-07-16 | 西南交通大学 | Thermal equivalent analysis method for internal defects of high-voltage bushing |
| CN114297822A (en) * | 2021-11-30 | 2022-04-08 | 国网辽宁省电力有限公司锦州供电公司 | Method for calculating pressure and temperature of contact finger of plum blossom contact of high-voltage switch cabinet |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060145680A1 (en) * | 2004-12-31 | 2006-07-06 | Tatung Co., Ltd. | Method for analyzing the inner temperature field and flow field in power transformers |
| CN101825501A (en) * | 2010-05-13 | 2010-09-08 | 哈尔滨工业大学 | Method for measuring power connector contact temperature rise based on finite element analysis of three-dimensional thermal field |
| CN106326572A (en) * | 2016-08-26 | 2017-01-11 | 国网四川省电力公司电力科学研究院 | Calculation method of real-time current-carrying capacity of cable connector |
| CN106570289A (en) * | 2016-11-10 | 2017-04-19 | 河北工业大学 | Contact temperature rise measurement method of low-voltage electric appliance under large current on the basis of finite element analysis |
| CN106644163A (en) * | 2016-11-10 | 2017-05-10 | 河北工业大学 | Miniature circuit breaker internal temperature rise measuring method based on temperature field finite element analysis |
| CN106934099A (en) * | 2017-02-10 | 2017-07-07 | 华南理工大学 | A kind of method for determining strain clamp current-carrying capacity and distribution of three-dimensional temperature |
-
2018
- 2018-03-01 CN CN201810170594.9A patent/CN108536906A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060145680A1 (en) * | 2004-12-31 | 2006-07-06 | Tatung Co., Ltd. | Method for analyzing the inner temperature field and flow field in power transformers |
| CN101825501A (en) * | 2010-05-13 | 2010-09-08 | 哈尔滨工业大学 | Method for measuring power connector contact temperature rise based on finite element analysis of three-dimensional thermal field |
| CN106326572A (en) * | 2016-08-26 | 2017-01-11 | 国网四川省电力公司电力科学研究院 | Calculation method of real-time current-carrying capacity of cable connector |
| CN106570289A (en) * | 2016-11-10 | 2017-04-19 | 河北工业大学 | Contact temperature rise measurement method of low-voltage electric appliance under large current on the basis of finite element analysis |
| CN106644163A (en) * | 2016-11-10 | 2017-05-10 | 河北工业大学 | Miniature circuit breaker internal temperature rise measuring method based on temperature field finite element analysis |
| CN106934099A (en) * | 2017-02-10 | 2017-07-07 | 华南理工大学 | A kind of method for determining strain clamp current-carrying capacity and distribution of three-dimensional temperature |
Non-Patent Citations (3)
| Title |
|---|
| 吴晓文: "气体绝缘母线温升计算及监测方法的研究", 《中国博士学位论文全文数据库》 * |
| 王青于等: "应用三维电磁-热-流耦合场分析法计算换流变压器干式套管的温度场分布", 《中国电机工程学报》 * |
| 陈利花: "电力线路及设备的电热混合场分析", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368385A (en) * | 2018-12-25 | 2020-07-03 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for operating characteristics of sleeve in electric-thermal coupling field |
| CN111368385B (en) * | 2018-12-25 | 2023-06-27 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for sleeve operation characteristics in electrothermal coupling field |
| CN110619138A (en) * | 2019-07-08 | 2019-12-27 | 国家电网有限公司 | Hardware fitting electric field distribution calculation method and device based on electric heating coupling field model |
| CN112179496A (en) * | 2020-09-15 | 2021-01-05 | 西安交通大学 | Infrared diagnosis operation and maintenance method for ultra-high voltage and extra-high voltage wall bushing with contact finger structure |
| CN112798643A (en) * | 2020-12-09 | 2021-05-14 | 国网宁夏电力有限公司检修公司 | Maintenance method, medium and system for contact finger of wall bushing |
| CN112798643B (en) * | 2020-12-09 | 2022-07-05 | 国网宁夏电力有限公司检修公司 | Maintenance method, medium and system for contact finger of wall bushing |
| CN113128059A (en) * | 2021-04-23 | 2021-07-16 | 西南交通大学 | Thermal equivalent analysis method for internal defects of high-voltage bushing |
| CN114297822A (en) * | 2021-11-30 | 2022-04-08 | 国网辽宁省电力有限公司锦州供电公司 | Method for calculating pressure and temperature of contact finger of plum blossom contact of high-voltage switch cabinet |
| CN114297822B (en) * | 2021-11-30 | 2024-04-16 | 国网辽宁省电力有限公司锦州供电公司 | Method for calculating pressure and temperature of contact finger of plum blossom-shaped contact of high-voltage switch cabinet |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108536906A (en) | A kind of bushing temperature field distribution method of the calculating containing connected structure | |
| CN104217061B (en) | Temperature field simulation design method for low-voltage distribution cabinet | |
| CN105092050A (en) | Method and device for detecting overheat fault of GIS (Gas Insulated Switchgear) internal conductors | |
| CN204116460U (en) | A kind of measuring system being applicable to insulating material high field conduction characteristic | |
| CN108228990A (en) | The life prediction method and device of disc insulator | |
| CN106644163A (en) | Miniature circuit breaker internal temperature rise measuring method based on temperature field finite element analysis | |
| CN109000825B (en) | Cable containing harmonic current and terminal core temperature calculation method thereof | |
| Lin et al. | 3‐D thermal analysis for the valve‐side bushings of HVDC converter transformer and application in the evaluation of current‐carrying connections defects | |
| Jing et al. | Research and analysis of power transformer remaining life prediction based on digital twin technology | |
| Yuan et al. | Research on temperature field simulation of dry type air core reactor | |
| CN106294966A (en) | A kind of list telegram in reply cable core transient state temperature rise acquisition methods being independent of skin temperature | |
| Tang et al. | A homotopy continuation method for nonlinear electric field analysis | |
| CN207425373U (en) | A kind of non-homogeneous electrically heated nuclear fuel simulation stick | |
| Zhang et al. | Fault analysis and simulation of ground-fault accident on transformer in 500 kV AC/DC grids | |
| Huang et al. | Digital twin modeling and temperature field analysis of transformer windings | |
| Yu et al. | Research on finite element model of air-core reactor based on magnetic-thermal coupling simulation | |
| Chenying et al. | Simulation and analysis of power cable harmonic loss and temperature rise | |
| CN110083905B (en) | Hot spot temperature analysis method for dynamic capacity increase of isolating switch | |
| Zhengmin et al. | Temperature-rise simulation and research of a new 252kV GIS disconnector based on Ansys workbench | |
| Xu et al. | Magneto-thermal coupling simulation and experimental verification for a three-winding high-frequency transformer | |
| Miao et al. | Multi-physics field coupling calculation and test of temperature field for dry-type air-core reactor | |
| Cheng et al. | Temperature rise prediction of gas-insulated transmission lines based on thermal network model considering contact fingers | |
| CN206311230U (en) | A kind of temperature sensor | |
| Shiling et al. | Analytical Derivation and Thermal Stability Analysis of Electro-thermal Model Considering Temperature Non-linearity of SF 6 High Voltage Structure | |
| Hou et al. | Application of Transformer Heat Circuit Model in Power Substation Station Area |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180914 |
|
| RJ01 | Rejection of invention patent application after publication |