CN109902443A - A kind of composite insulator internal flaw research method based on COMSOL - Google Patents

A kind of composite insulator internal flaw research method based on COMSOL Download PDF

Info

Publication number
CN109902443A
CN109902443A CN201910250440.5A CN201910250440A CN109902443A CN 109902443 A CN109902443 A CN 109902443A CN 201910250440 A CN201910250440 A CN 201910250440A CN 109902443 A CN109902443 A CN 109902443A
Authority
CN
China
Prior art keywords
composite insulator
electric field
comsol
internal flaw
research method
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
Application number
CN201910250440.5A
Other languages
Chinese (zh)
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.)
Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd
North China Electric Power University
State Grid Hebei Energy Technology Service Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd
North China Electric Power University
State Grid Hebei Energy Technology Service 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 Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd, North China Electric Power University, State Grid Hebei Energy Technology Service Co Ltd filed Critical Electric Power Research Institute of State Grid Hebei Electric Power Co Ltd
Priority to CN201910250440.5A priority Critical patent/CN109902443A/en
Publication of CN109902443A publication Critical patent/CN109902443A/en
Pending legal-status Critical Current

Links

Landscapes

  • Insulators (AREA)

Abstract

The composite insulator internal flaw research method based on COMSOL that the invention discloses a kind of, a kind of composite insulator internal flaw research method based on COMSOL, the composite insulator threedimensional model of 1:1 is established including step A, using the geometry module in COMSOL simulation software;Step B, different degrees of conductive defect is set in composite insulator zero defect and different location, analyzes its space axial electric field.The present invention establishes composite insulator zero defect and the asynchronous threedimensional model of position difference defect level respectively, have studied the distribution character of its space axial electric field, the present invention can improve the deficiency of existing method, by analyzing axial electric field curve, it can be found that influence of the different defects to Insulators Used.

Description

A kind of composite insulator internal flaw research method based on COMSOL
Technical field
The present invention relates to composite insulator axial electric field research field more particularly to a kind of compound inslations based on COMSOL Sub- internal flaw research method.
Background technique
Composite insulating material gradually substitutes traditional glass and pottery since weight is small and the functional feature of anti-pollution Porcelain insulating material.However, the limitation widely applied main problem of composite insulator is that existing live monitoring method cannot The state of insulation of composite insulator is diagnosed well.Live detection method both domestic and external mainly uses visual examination method, infrared ray Thermometry, ultraviolet imaging method.Faulty insulator electric discharge dynamic characteristic is detected based on ultraviolet imagery technology, when deterioration is exhausted When edge is located in the middle part of insulator chain, sensitivity is lower.With infrared temperature-test technology, based on insulator chain temperature curve feature Faulty insulator diagnostic method, and the insulator in 220kV substation is detected to improve correctness.But degradation compared with When low, temperature rise difference is not obvious.
Summary of the invention
The object of the present invention is to provide a kind of composite insulator internal flaw research method based on COMSOL, energy of the present invention The deficiency for enough improving existing method, to insulator, there are the variations of the axial electric field of internal flaw comprehensive understanding, and judgement is fallen vacant Fall into position.
In order to solve the above technical problems, the present invention adopts the following technical scheme:
A kind of composite insulator internal flaw research method based on COMSOL of the present invention, it is a kind of based on the compound of COMSOL Interior insulator Study of Defects method, including
Step A, the composite insulator threedimensional model of 1:1 is established using the geometry module in COMSOL simulation software;
Step B, different degrees of conductive defect is set in composite insulator zero defect and different location, analyzes its space Axial electric field.
Further, in step B, include the following steps,
It step B1, is that corresponding material is arranged in each part of model according to the composite insulator string of actual motion;
Step B2, the boundary condition of emulation is combined to apply corresponding electricity to compound inslation submodel according to the actual operation Gesture;
Step B3, the subdivision of different scale is carried out for different piece;
Step B4, space axial electric field curve graph is calculated to obtain.
Further, in step B1, the material of full skirt and sheath is silicon rubber, and head and tail portion connection gold utensil material are Stainless steel.
Further, in step B2, comprising the following steps:
In the case where ignoring space charge effects, solved using corresponding electrostatic methods, therefore insulator chain surface And the electric field in surrounding medium solves and meets:
In formula: parameter E is the electric field strength of dielectric area, unit N/C;Parameter ε is the dielectric constant of medium;Parameter ρ For space charge volume density, unit c/m3;Potential functionMeet formula (3)
When calculating the field distribution in composite insulator umbrella skirt surface and surrounding medium, the relative dielectric constant of air-shed is 1;The relative dielectric constant of full skirt is 6;The relative dielectric constant of sheath is 4;The dielectric constant of steel feet and chapeau de fer is 1000;
The calculating of composite insulator string space electric field is an open domain computational problem, but open domain is can not directly to carry out accurately The problem of calculating, converting finite field for infinitely great air-shed problem using finite element modeling technology;When calculating, using 1 class boundary condition are as follows:
In above-mentioned formula: Γ0For air-shed boundary and insulator low-pressure end;ΓkFor insulator high-voltage end;Formula (4)-(5) group At complete boundary value problem.
Further, in step B3, by local refinement subdivision method, high-precision is carried out for insulation subregion and is cutd open Point, minimum grid selection of dimension 5mm is to guarantee computational accuracy;Air-shed big to size and regular simultaneously carries out rough subdivision To reduce unnecessary calculation amount in simulation process, therefore the maximum mesh size of threedimensional model air-shed is taken as 50mm, Ensure under conditions of certain precision and can improve to a certain extent calculating speed.
Further, in step B4, space axial direction curve is the axis that is calculated at insulator case surface 3mm To electric field.
Compared with prior art, advantageous effects of the invention are as follows:
The present invention establishes composite insulator zero defect and the asynchronous threedimensional model of position difference defect level respectively, grinds The distribution character of its space axial electric field is studied carefully, the present invention can improve the deficiency of existing method, bent by analysis axial electric field Line, it can be found that influence of the different defects to Insulators Used.
Detailed description of the invention
The invention will be further described for explanation with reference to the accompanying drawing.
Fig. 1 is electric field component measurement point distribution figure in mesohigh end of the present invention;
Fig. 2 is intermediate voltage terminal electric field component measurement point distribution figure in the present invention;
Fig. 3 is electric field component measurement point distribution figure in mesolow end of the present invention;
Fig. 4 is the line chart that axial electric field is calculated in high-voltage end in different measurement points for defect;
Fig. 5 is the line chart that axial electric field is calculated in intermediate ends in different measurement points for defect;
Fig. 6 is the line chart that axial electric field is calculated in ground terminal in different measurement points for defect;
Specific embodiment
Referring to Fig.1-3, a specific embodiment of the invention the following steps are included:
Step A, the 110kV composite insulator threedimensional model of 1:1 is established using the geometry module in COMSOL simulation software.
Step B, different degrees of conductive defect is set in composite insulator zero defect and different location, analyzes its space Axial electric field.Specifically comprise the following steps:
It B1, is that corresponding material, this specific implementation is arranged in each part of model according to the composite insulator string of actual motion In the step B1 of example, the material of full skirt and sheath is silicon rubber, and head and tail portion connection gold utensil material are stainless steel.
B2, the boundary condition of emulation is combined to apply corresponding potential to compound inslation submodel according to the actual operation. In the case where ignoring space charge effects, can be solved using corresponding electrostatic methods.Therefore insulator chain surface and week It encloses the electric field in medium and solves and meet:
In formula: parameter E is the electric field strength of dielectric area, unit N/C;Parameter ε is the dielectric constant of medium;Parameter ρ For space charge volume density, unit c/m3.Potential functionMeet formula (3)
When calculating the field distribution in composite insulator umbrella skirt surface and surrounding medium, the relative dielectric constant of air-shed is 1;The relative dielectric constant of full skirt is 6;The relative dielectric constant of sheath is 4;The dielectric constant of steel feet and chapeau de fer is 1000.
The calculating of composite insulator string space electric field is an open domain computational problem, but open domain is can not directly to carry out accurately It calculates, the problem of present invention converts finite field for infinitely great air-shed problem using finite element modeling technology, to solve It has determined the thorny problem that open domain is directly calculated.When calculating, using the 1st class boundary condition are as follows:
In above-mentioned formula: Γ0For air-shed boundary and insulator low-pressure end;ΓkFor insulator high-voltage end.Formula (4)-(5) group At complete boundary value problem.
B3, comprehensively consider operation time and precision, the subdivision of different scale is carried out for different piece.The corresponding boundary values of electric field Problem is solved on COMSOL Multiphysics software using finite element numerical computational methods.Using tetrahedron conduct Basic unit carries out mesh generation to entire model.In view of the size of insulator is much smaller than the size of entire air-shed, in order to Operation efficiency is improved, software runing time is reduced, selects local refinement subdivision method.It, can be with by local refinement subdivision method The region of the tiny complexity such as full skirt, plug for insulator carries out high-precision subdivision, and minimum grid selection of dimension 5mm is to guarantee Computational accuracy.Air-shed big to size and regular simultaneously carries out rough subdivision to reduce unnecessary meter in simulation process Calculation amount, therefore the maximum mesh size of threedimensional model air-shed can be taken as 50mm, under conditions of ensuring certain precision again Calculating speed can be improved to a certain extent.Model quantity after subdivision there are about 2,500,000 units, lay by Simulation Calculation 36 core processors, 128GB memory work station on, computation model single calculation needs the used time about 1.5 hours.
B4, space axial electric field curve graph is calculated to obtain, space axial direction curve is to calculate at insulator case surface 3mm Obtained axial electric field.
The line chart of axial electric field is calculated in different measurement points in high-voltage end, intermediate ends, ground terminal respectively in defect And analysis result.
One, internal flaw is in high-voltage end
From fig. 4, it can be seen that " U " font is presented in axial electric field distribution curve when insulator chain is intact, maximum field is strong Degree is 9.47kV/cm, and field minimum intensity is 0.21kV/cm.Due to the presence of internal flaw, so that axial electric field strength is in phase The position answered can decline.Meanwhile near the terminal position of existing defects, axial electric field strength is significantly increased, and with scarce It falls into the increase of severity and increases.This is because caused by the internal field near defect tip increases.When defect cylindrical body When length is 35mm, the field strength of the 2nd measurement point is 7.78kV/cm, and the field strength relative to adjacent measurement points is big.Work as defect When cylinder length is 105mm, the field strength of the 4th measurement point is 6.85kV/m, and the field strength relative to adjacent measurement points is big. When defect cylinder length is 210mm, the field strength of the 4th measurement point is 5.14kV/m, the field relative to adjacent measurement points Qiang Yao great.
Two, internal flaw is in intermediate ends
As can be known from Fig. 5, when internal flaw is in intermediate ends, the variable quantity of axial electric field strength and radial electric field intensity All increase with the increase of defect severity.From figure 3, it can be seen that the electric field strength near defect terminal position becomes Change amount is big.
When defect cylinder length is 35mm, the field strength of the 12nd and the 13rd measurement point is relative to adjacent measurement points Field strength it is big.When defect cylinder length is 105mm, the field strength of the 11st and the 14th measurement point is relative to adjacent survey The field strength of amount point is big.When defect cylinder length is 210mm, the field strength of the 9th and the 16th measurement point is relative to phase The field strength of adjacent measurement point is big.There are the different of internal flaw with high-voltage end is the electric field of defect both-side ends attachment for this Intensity all changes, therefore this feature helps to detect the position of internal flaw and severity.
Three, internal flaw is in ground terminal
As can be known from Fig. 6, due to the presence of internal flaw, the electric field strength near defect terminal position is changed significantly increasing Greatly.When defect cylinder length is 35mm, the field strength of the 23rd measurement point is big relative to the field strength of adjacent measurement points.When When defect cylinder length is 105mm, the field strength of the 21st measured place is big relative to the field strength of adjacent measurement points.When defect circle When cylinder length is 210mm, the field strength of the 18th measured place is big relative to the field strength of adjacent measurement points.
From the present invention:
1) conductive defect of composite insulator inside different location causes field distribution variation to be characterized in different, Electric field strength near defect terminal position is changed significantly increase.This facilitates the position for detecting internal flaw.
2) there are when conductive defect inside composite insulator, the radial electric field component of measurement point is distributed between sub-pieces Variation to axial electric field component distribution variation to become apparent.
3) there are the conductive defects of different severity to cause the amount of distortion of field distribution can be with inside composite insulator The increase of the severity of conductive defect and increase, detect the severity of internal flaw.
Embodiment described above is only that preferred embodiment of the invention is described, and is not carried out to the scope of the present invention It limits, without departing from the spirit of the design of the present invention, those of ordinary skill in the art make technical solution of the present invention Various changes and improvements, should all fall into claims of the present invention determine protection scope in.

Claims (6)

1. a kind of composite insulator internal flaw research method based on COMSOL, it is characterised in that: including
Step A, the composite insulator threedimensional model of 1:1 is established using the geometry module in COMSOL simulation software;
Step B, different degrees of conductive defect is set in composite insulator zero defect and different location, analyzes its space axial direction Electric field.
2. the composite insulator internal flaw research method according to claim 1 based on COMSOL, it is characterised in that: step In rapid B, include the following steps,
It step B1, is that corresponding material is arranged in each part of model according to the composite insulator string of actual motion;
Step B2, the boundary condition of emulation is combined to apply corresponding potential to compound inslation submodel according to the actual operation;
Step B3, the subdivision of different scale is carried out for different piece;
Step B4, space axial electric field curve graph is calculated to obtain.
3. the composite insulator internal flaw research method according to claim 2 based on COMSOL, it is characterised in that: step In rapid B1, the material of full skirt and sheath is silicon rubber, and head and tail portion connection gold utensil material are stainless steel.
4. the composite insulator internal flaw research method according to claim 2 based on COMSOL, it is characterised in that: step In rapid B2, comprising the following steps:
In the case where ignoring space charge effects, solved using corresponding electrostatic methods, therefore insulator chain surface and week It encloses the electric field in medium and solves and meet:
In formula: parameter E is the electric field strength of dielectric area, unit N/C;Parameter ε is the dielectric constant of medium;Parameter ρ is sky Between electric charge volume density, unit c/m3;Potential functionMeet formula (3)
When calculating the field distribution in composite insulator umbrella skirt surface and surrounding medium, the relative dielectric constant of air-shed is 1;Umbrella The relative dielectric constant of skirt is 6;The relative dielectric constant of sheath is 4;The dielectric constant of steel feet and chapeau de fer is 1000;
The calculating of composite insulator string space electric field is an open domain computational problem, but open domain is directly to be accurately calculated , the problem of converting finite field for infinitely great air-shed problem using finite element modeling technology;When calculating, using the 1st class Boundary condition are as follows:
In above-mentioned formula: Γ0For air-shed boundary and insulator low-pressure end;ΓkFor insulator high-voltage end;Formula (4)-(5) constitute Complete boundary value problem.
5. the composite insulator internal flaw research method according to claim 2 based on COMSOL, it is characterised in that: step In rapid B3, by local refinement subdivision method, high-precision subdivision, minimum grid selection of dimension 5mm are carried out for insulation subregion To guarantee computational accuracy;Air-shed big to size and regular simultaneously carries out rough subdivision need not in simulation process to reduce The calculation amount wanted, therefore the maximum mesh size of threedimensional model air-shed is taken as 50mm, under conditions of ensuring certain precision Calculating speed can be improved to a certain extent again.
6. the composite insulator internal flaw research method according to claim 2 based on COMSOL, it is characterised in that: step In rapid B4, space axial direction curve is the axial electric field that is calculated at insulator case surface 3mm.
CN201910250440.5A 2019-03-29 2019-03-29 A kind of composite insulator internal flaw research method based on COMSOL Pending CN109902443A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910250440.5A CN109902443A (en) 2019-03-29 2019-03-29 A kind of composite insulator internal flaw research method based on COMSOL

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910250440.5A CN109902443A (en) 2019-03-29 2019-03-29 A kind of composite insulator internal flaw research method based on COMSOL

Publications (1)

Publication Number Publication Date
CN109902443A true CN109902443A (en) 2019-06-18

Family

ID=66955009

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910250440.5A Pending CN109902443A (en) 2019-03-29 2019-03-29 A kind of composite insulator internal flaw research method based on COMSOL

Country Status (1)

Country Link
CN (1) CN109902443A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111829663A (en) * 2020-06-19 2020-10-27 南方电网科学研究院有限责任公司 Composite insulator defect classification diagnosis method based on surface temperature distribution
CN112182940A (en) * 2020-10-23 2021-01-05 国网天津市电力公司 High-voltage switch cabinet insulation partition plate arrangement optimization method considering microenvironment
CN113076672A (en) * 2021-04-08 2021-07-06 江苏省送变电有限公司 Composite insulator bird pecking defect risk grade evaluation method based on improved SSD
CN113325280A (en) * 2021-05-26 2021-08-31 国网陕西省电力公司电力科学研究院 GIS insulator defect electric field simulation method based on finite element analysis
CN113486542A (en) * 2021-02-02 2021-10-08 上海大学 Slag applied electric field modeling method based on COMSOL
CN114155756A (en) * 2021-11-11 2022-03-08 国网江苏省电力有限公司技能培训中心 Porcelain insulator crack defect simulation device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771924A (en) * 2016-12-28 2017-05-31 国网内蒙古东部电力有限公司检修分公司 A kind of utilization photoelectricity field sensor detects the detecting system and method for defects of insulator
CN108717071A (en) * 2018-08-31 2018-10-30 国网河北能源技术服务有限公司 A kind of composite insulator umbrella skirt breakage diagnostic method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771924A (en) * 2016-12-28 2017-05-31 国网内蒙古东部电力有限公司检修分公司 A kind of utilization photoelectricity field sensor detects the detecting system and method for defects of insulator
CN108717071A (en) * 2018-08-31 2018-10-30 国网河北能源技术服务有限公司 A kind of composite insulator umbrella skirt breakage diagnostic method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GBAH KONE: "Numerical investigation of electric field distortion induced by internal defects in composite insulators", 《HIGH VOLT》 *
王永强等: "GIS盆式绝缘子气隙缺陷下电场变化规律仿真研究", 《绝缘材料》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111829663A (en) * 2020-06-19 2020-10-27 南方电网科学研究院有限责任公司 Composite insulator defect classification diagnosis method based on surface temperature distribution
CN112182940A (en) * 2020-10-23 2021-01-05 国网天津市电力公司 High-voltage switch cabinet insulation partition plate arrangement optimization method considering microenvironment
CN112182940B (en) * 2020-10-23 2022-03-11 国网天津市电力公司 High-voltage switch cabinet insulation partition plate arrangement optimization method considering microenvironment
CN113486542A (en) * 2021-02-02 2021-10-08 上海大学 Slag applied electric field modeling method based on COMSOL
CN113076672A (en) * 2021-04-08 2021-07-06 江苏省送变电有限公司 Composite insulator bird pecking defect risk grade evaluation method based on improved SSD
CN113076672B (en) * 2021-04-08 2024-02-02 江苏省送变电有限公司 Composite insulator bird pecking defect risk level assessment method based on improved SSD
CN113325280A (en) * 2021-05-26 2021-08-31 国网陕西省电力公司电力科学研究院 GIS insulator defect electric field simulation method based on finite element analysis
CN114155756A (en) * 2021-11-11 2022-03-08 国网江苏省电力有限公司技能培训中心 Porcelain insulator crack defect simulation device
CN114155756B (en) * 2021-11-11 2024-03-26 国网江苏省电力有限公司技能培训中心 Porcelain bottle crack defect simulation device

Similar Documents

Publication Publication Date Title
CN109902443A (en) A kind of composite insulator internal flaw research method based on COMSOL
CN104776938B (en) Cable connector core temperature inversion method and system based on cable surface temperature
CN105676081B (en) Direct current cables based on ultralow frequency dielectric loss squeezes out insulation ag(e)ing state evaluating method
CN102520323B (en) Hyperspectrum-based composite insulator aging operating state detection method
CN103630803B (en) A kind of cable partial discharge live correction method
CN101639502A (en) System for automatically measuring charge distribution on surface of solid medium
CN109239539A (en) A kind of device and experimental method for studying epoxy insulation material Characteristics of Partial Discharge
CN102519846B (en) Hyperspectrum-based composite insulator hydrophobicity detection method
CN108387826B (en) Ultra-high voltage split conductor corona onset field strength prediction method considering altitude correction
CN1719269A (en) Insulator charged detection instrument and its implement method
CN102707131B (en) Very fast transient overvoltage (VFTO) automatic measuring system for high-sensitivity medium window
Ma et al. Measurement and simulation of charge accumulation on a disc spacer with electro-thermal stress in SF 6 gas
CN111679103A (en) High-voltage cable water-blocking tape electrical performance test system based on three-electrode method
CN210775723U (en) Simulation test device for heating characteristics of zero-value insulator string
CN113158383B (en) Method for evaluating actual field intensity of direct-current cable joint by using coaxial simplified model
Bian et al. The effect of surface roughness on corona-generated electromagnetic interference for long-term operating conductors
CN103901329B (en) Method for determining corona discharge ultraviolet detection sensitivity of high-voltage electrical equipment
CN105572510B (en) A kind of DC power transmission line corona loss method of testing based on corona cage
CN112924763A (en) Method for evaluating transient electric shock below high-voltage alternating-current transmission line
CN114400118A (en) Aging method of composite insulator core rod
Wei et al. Inverse application of charge simulation method in detecting faulty ceramic insulators and processing influence from tower
Zheng et al. Three-dimensional modeling and optimization of a fully-covered interdigital capacitive sensor for power cable insulation detection
CN116430176A (en) Diagnosis method for transformer oil paper insulation state under lightning impulse voltage
CN115128447A (en) Method for detecting state quantity of middle casing pipe in GIS of thermal power energy storage system
CN109341926A (en) Contraction Type Electric Cable Attachment in Crane Ore Group interfacial pressure detection method is pressed in one kind

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: 20190618

RJ01 Rejection of invention patent application after publication