CN106815436A - The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field - Google Patents
The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field Download PDFInfo
- Publication number
- CN106815436A CN106815436A CN201710037842.8A CN201710037842A CN106815436A CN 106815436 A CN106815436 A CN 106815436A CN 201710037842 A CN201710037842 A CN 201710037842A CN 106815436 A CN106815436 A CN 106815436A
- Authority
- CN
- China
- Prior art keywords
- phase
- short circuit
- subdivision
- transformers
- phase transformers
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Abstract
The embodiment of the invention discloses the computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field, for solving because the model of current large-scale 220kV transformers is more complicated, the nonlinear characteristic of material is difficult to obtain, the amount of calculation of solution than larger, caused by the low technical problem of large-scale 220kV transformers 3 D electromagnetic solution efficiency.Present invention method includes:Set up 220kV three-phase transformer simplified models;It is several three-dimensional grid subdivision units by the 220kV three-phase transformers simplified model subdivision;The short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and driving source is set, carry out Electromagnetic Transient Analysis.
Description
Technical field
The present invention relates to technical field of electric power, more particularly to a kind of 220kV three-phase transformers B phases short circuit electromagnetic field calculating
Method and device.
Background technology
With continuing to develop for power construction, power equipment develops towards maximization direction.With transformer voltage grade
Be substantially improved, as the safe and stable operation of the transformer of one of power system main carriers, be related to that the entirety of power network resists
Risk ability.How to ensure transformer be chronically at good working order always power utility plant management weight in it
Weight.In recent years, repeatedly there is main transformer grave gas trip protection, pressure relief valve action caused by external short circuit impacts in network system
Oil spout etc..When because external circuit short circuit and caused by rush of current when winding conducting wire in by short circuit current numerical value can reach and
15~20 times of fixed number value.Transformer will produce larger electromagnetic force and electromagnetic consumable within the extremely short time.In electromagnetic force
In the presence of, winding is it may happen that vibration, so that the oil stream in influenceing fuel tank is moved, produces oil flow surge.Will accurately really
Oil stream motion in stand oil case is, it is necessary to the electromagnetism that the electromagnetic field intensity produced in short-circuit condition to transformer is distributed, winding is subject to
The temperature rise of power and oil stream has accurate understanding.
The stress of winding in transformer is accurately calculated, first should to transformer short-circuit when electromagnetic field carry out
It is accurate to calculate.Therefore, FInite Element is in the magnetic field analysis for being drawn into transformer very early.For the loss of transformer
Calculate, the calculating means often having are the time domain method with analytic formula and semiempirical formula, Orthogonal Decomposition synthetic method, but limitation
It is too big.Also generally all it is using 2D FInite Elements in terms of using computer-assisted analysis.
These methods can to a certain extent solve engineering problem.But the simplification that these methods contain in itself is too many,
The result applicability of calculating is smaller, there is provided data be difficult to support further transformer vibrate or oil flow surge analysis.
Therefore, the 3 D electromagnetic for this kind of large-scale 220kV transformers is solved, and what is be analyzed is rarely found.
To find out its cause, the model of its large-scale 220kV transformer is more complicated, the nonlinear characteristic of material is difficult to obtain, and asks
The amount of calculation of solution is than larger, low so as to result in large-scale 220kV transformers 3 D electromagnetic solution efficiency technical problem.
The content of the invention
The embodiment of the invention provides the computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field, solution
Determine because the model of current large-scale 220kV transformers is more complicated, the nonlinear characteristic of material is difficult to obtain, the meter of solution
Calculation amount than larger, caused by the low technical problem of large-scale 220kV transformers 3 D electromagnetic solution efficiency.
A kind of computational methods of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention, including:
Set up 220kV three-phase transformer simplified models;
It is several three-dimensional grid subdivision units by the 220kV three-phase transformers simplified model subdivision;
The short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and driving source is set,
Carry out Electromagnetic Transient Analysis.
Alternatively, it is that several three-dimensional grid subdivision units are specific by the 220kV three-phase transformers simplified model subdivision
Including:
It is several tetrahedron two-order units by the 220kV three-phase transformers simplified model subdivision.
Alternatively, the short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and is set
Driving source, carries out Electromagnetic Transient Analysis and specifically includes:
The short circuit in winding treatment corresponding with B is carried out to tetrahedron two-order unit several described, and it is short out to set dispatch from foreign news agency
Road transient state driving source;
The vector position free degree to the seamed edge of the tetrahedron two-order unit carries out the calculating of single order unit, to the tetrahedron two
The scalar potential free degree of the node of rank unit carries out second order unit and calculates, and solves Three dimensional transient magnetic field;
Result of calculation analysis is carried out to the Three dimensional transient magnetic field.
Alternatively, result of calculation analysis is carried out to the Three dimensional transient magnetic field to specifically include:
The Three dimensional transient magnetic field according to solving is analyzed the determination 220kV three-phase transformers simplified model pair
Internal magnetic field distribution character, each phase winding and the iron core ohmic loss answered, and each phase winding stressing conditions.
Alternatively, being analyzed according to the Three dimensional transient magnetic field for solving also includes:
Get the corresponding three-phase current transient wave of the 220kV three-phase transformers simplified model under B phase short-circuit conditions
The corresponding Magnetic Induction Density Distribution cloud atlas of the 220kV three-phase transformers simplified model and in specified path under shape, preset time instant
The Magnetic Induction Density Distribution curve map for obtaining is investigated under position.
A kind of computing device of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention, its feature exists
In, including:
Unit is set up, for setting up 220kV three-phase transformer simplified models;
Subdivision unit, for being several three-dimensional grid subdivision lists by the 220kV three-phase transformers simplified model subdivision
Unit;
Analytic unit, for carrying out the short circuit in winding treatment corresponding with B to three-dimensional grid subdivision unit several described,
And driving source is set, carry out Electromagnetic Transient Analysis.
Alternatively, subdivision unit, specifically for by the 220kV three-phase transformers simplified model subdivision for several on four sides
Body two-order unit.
Alternatively, analytic unit includes:
Short circuit treatment subelement, for carrying out the short circuit in winding corresponding with B to tetrahedron two-order unit several described
Treatment, and external circuit short circuit transient state driving source is set;
Computation subunit, by the seamed edge to the tetrahedron two-order unit vector position the free degree carry out single order unit based on
Calculate, the scalar potential free degree to the node of the tetrahedron two-order unit carries out the calculating of second order unit, solves Three dimensional transient magnetic field;
Analysis subelement, for carrying out result of calculation analysis to the Three dimensional transient magnetic field.
Alternatively, analysis subelement includes:
First analysis module, for being analyzed the determination transformer letter according to the Three dimensional transient magnetic field for solving
Change the corresponding internal magnetic field distribution character of model, each phase winding and iron core ohmic loss, and each phase winding stressing conditions.
Alternatively, analysis subelement also includes:
Second analysis module, for being analyzed according to the Three dimensional transient magnetic field for solving, gets in the short circuit of B phases
The 220kV tri- under the corresponding three-phase current transient waveform of the 220kV three-phase transformers simplified model, preset time instant under state
The corresponding Magnetic Induction Density Distribution cloud atlas of phase transformer simplified model and that the magnetic induction for obtaining is investigated under specified path position is strong
Degree scatter chart.
As can be seen from the above technical solutions, the embodiment of the present invention has advantages below:
The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention, bag
Include:Set up 220kV three-phase transformer simplified models;By the 220kV three-phase transformers simplified model subdivision for several are three-dimensional
Mesh generation unit;The short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and sets sharp
Encourage source, carry out Electromagnetic Transient Analysis, solve because the model of current large-scale 220kV transformers is more complicated, material it is non-
Linear characteristic is difficult to obtain, the amount of calculation of solution than larger, caused by large-scale 220kV transformers 3 D electromagnetic solution efficiency
Low technical problem.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also
Other accompanying drawings are obtained with according to these accompanying drawings.
Fig. 1 is the one of the computational methods of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention
The schematic flow sheet of individual embodiment;
Fig. 2 is a kind of the another of the computational methods of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention
The schematic flow sheet of one embodiment;
Fig. 3 is an iron core model graphics;
Fig. 4 is another iron core model graphics;
Fig. 5 is transformer model figure;
Fig. 6 is iron core BH curve maps;
Fig. 7 is iron core BP curve maps;
Fig. 8 is B phases short circuit transient current oscillogram;
Fig. 9 is iron core mesh generation schematic diagram;
Figure 10 is iron core, coil, arm-tie and folder mesh generation schematic diagram;
Figure 11 is integral grid subdivision front view;
Figure 12 is lattice statistical amount schematic diagram;
Figure 13 is second order tetrahedral grid schematic diagram;
High side voltage and testing current waveform diagram when Figure 14 is the short circuit of B phases;
External circuit schematic diagram when Figure 15 is the short circuit of B phases high pressure;
High voltage side current waveform diagram when Figure 16 is the short circuit of B phases;
Current voltage waveform diagram in when Figure 17 is the short circuit of B phases;
Figure 18 is certain moment Magnetic Induction Density Distribution cloud atlas of B phases short circuit;
Figure 19 is certain moment magnetic flux density vector figure of B phases short circuit;
Transformer core field strength distribution cloud atlas when Figure 20 is the short circuit of B phases;
Transformer iron core magnetic vector distribution map when Figure 21 is the short circuit of B phases;
Figure 22 is a kind of computing device of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention
The structural representation of one embodiment;
Figure 23 is a kind of computing device of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention
The structural representation of another embodiment.
Specific embodiment
The embodiment of the invention provides the computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field, solution
Determine because the model of current large-scale 220kV transformers is more complicated, the nonlinear characteristic of material is difficult to obtain, the meter of solution
Calculation amount than larger, caused by the low technical problem of large-scale 220kV transformers 3 D electromagnetic solution efficiency.
To enable that goal of the invention of the invention, feature, advantage are more obvious and understandable, below in conjunction with the present invention
Accompanying drawing in embodiment, is clearly and completely described, it is clear that disclosed below to the technical scheme in the embodiment of the present invention
Embodiment be only a part of embodiment of the invention, and not all embodiment.Based on the embodiment in the present invention, this area
All other embodiment that those of ordinary skill is obtained under the premise of creative work is not made, belongs to protection of the present invention
Scope.
This electromagnetism is calculated and emulated using ANSYS Maxwell R17.0 low frequency electromagnetic field finite element analysis softwares,
ANSYS Maxwell R17.0 are integrated in ANSYS Electronic Desktop R17.0 electronic desktops.ANSYS Maxwell
It is global foremost electromagnetic finite element analysis software, is widely used in motor, actuator, high-voltage electrical equipment, height and presses off
The Electromagnetic Calculation and scheme optimization of the electronic products such as pass and power product.
Refer to Fig. 1, a kind of calculating side of 220kV three-phase transformers B phases short circuit electromagnetic field provided in an embodiment of the present invention
One embodiment of method includes:
101st, 220kV three-phase transformer simplified models are set up;
In the present embodiment, when needing to be analyzed three-phase transformer 3 D electromagnetic field, get pass through first
The 220kV three-phase transformer simplified models that transformer basic-element model in Maxwell instruments is set up.
102nd, it is several three-dimensional grid subdivision units by 220kV three-phase transformer simplified models subdivision;
Mould is simplified by the 220kV three-phase transformers of the transformer basic-element model foundation in Maxwell instruments when getting
Type, it is necessary to pass through the static adaptive mesh of Maxwell instruments by 220kV three-phase transformer simplified models subdivision for several three
Dimension mesh generation unit.
103rd, the short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and sets sharp
Source is encouraged, Electromagnetic Transient Analysis are carried out.
When transformer simplified model subdivision is several three-dimensional grids by the static adaptive mesh by Maxwell instruments
After subdivision unit, it is necessary to pass through Maxwell instruments three-dimensional grid subdivision unit several described is carried out it is corresponding with B around
Group short circuit treatment, and driving source is set, carry out Electromagnetic Transient Analysis.
Refer to Fig. 2, a kind of computational methods of three-phase transformer 3 D electromagnetic provided in an embodiment of the present invention it is another
Individual embodiment includes:
201st, 220kV three-phase transformer simplified models are set up;
In the present embodiment, when needing to be analyzed three-phase transformer 3 D electromagnetic field, get pass through first
The 220kV three-phase transformer simplified models that transformer basic-element model in Maxwell instruments is set up.
202nd, it is several tetrahedron two-order units by 220kV three-phase transformer simplified models subdivision;
Mould is simplified by the 220kV three-phase transformers of the transformer basic-element model foundation in Maxwell instruments when getting
Type, it is necessary to pass through the static adaptive mesh of Maxwell instruments by 220kV three-phase transformer simplified models subdivision for several four
Face body two-order unit.
203rd, several tetrahedron two-order units are carried out with the short circuit in winding treatment corresponding with B, and it is short out to set dispatch from foreign news agency
Road transient state driving source;
When 220kV three-phase transformer simplified models subdivision is some by the static adaptive mesh by Maxwell instruments
After individual tetrahedron two-order unit, it is necessary to pass through Maxwell instruments several tetrahedron two-order units are carried out it is corresponding with B
Short circuit in winding treatment, and external circuit short circuit transient state driving source is set.
204th, the vector position free degree to the seamed edge of tetrahedron two-order unit carries out the calculating of single order unit, to tetrahedron second order list
The scalar potential free degree of the node of unit carries out the calculating of second order unit and solves Three dimensional transient magnetic field;
When several tetrahedron two-order units are carried out by Maxwell instruments correspondence short circuit in winding treatment after, it is necessary to
The vector position free degree to the seamed edge of tetrahedron two-order unit carries out the calculating of single order unit, to the mark of the node of tetrahedron two-order unit
The amount position free degree carries out the calculating of second order unit and solves Three dimensional transient magnetic field.
205th, it is analyzed according to the Three dimensional transient magnetic field for solving and determines that 220kV three-phase transformer simplified models are corresponding
Internal magnetic field distribution character, each phase winding and iron core ohmic loss, and each phase winding stressing conditions;
Solve after Three dimensional transient magnetic field, it is necessary to be analyzed determination 220kV according to the Three dimensional transient magnetic field for solving
The corresponding internal magnetic field distribution character of three-phase transformer simplified model, each phase winding and iron core ohmic loss, and each phase winding
Stressing conditions.
206th, the corresponding three-phase current transient wave of 220kV three-phase transformer simplified models under B phase short-circuit conditions is got
The corresponding Magnetic Induction Density Distribution cloud atlas of 220kV three-phase transformers simplified model and in specified path position under shape, preset time instant
It is lower to investigate the Magnetic Induction Density Distribution curve map for obtaining.
Solve after Three dimensional transient magnetic field, it is necessary to get the corresponding three-phase electricity of transformer simplified model under short-circuit condition
Stream transient waveform, the corresponding Magnetic Induction Density Distribution cloud atlas of transformer simplified model and under specified path position under preset time instant
The Magnetic Induction Density Distribution curve map that investigation is obtained.
It is described in detail with a concrete application scene below, application examples includes:
The non-shorting phase winding of transformer high, medium and low voltage is open circuit, and voltage regulation coil is not powered on stream, single-phase (taking B phases) senior middle school
Pressure short circuit in winding, carries out transient magnetic field analysis.
Power transformer grade 110KV, rated capacity 150MVA, nameplate and structural parameters are shown in Table 1.
Table 1
Inside transformer iron core is stacked using cold-reduced silicon sheet, calculates the modeling function using Maxwell to emulation
Object (i.e. 220kV three-phase transformers) completes the foundation of model.Coil uses cake formula structure, and structural parameters such as table 1 is calculated and used
Maxwell carries transformer coil basic-element model, and model as shown in Figure 3 and Figure 4, is built to transformer core and coil first
Mould.
Other structures, such as arm-tie, folder, Oiltank structure parameter are as shown in table 2.
Table 2
The transformer block mold set up based on Maxwell is as shown in Figure 5.
Material parameter is shown in Table 3:
Table 3
When transformer occurs B phase high pressure short circuit in windings, the transient current waveform of its driving source is as shown in Figure 8.
T- Ω algorithms are still used in Three dimensional transient, for low frequency transient magnetic field, maxwell equation group can
In the form of being written as formula (1)-formula (3).
ANSYS Maxwell softwares magnetic FEM calculates analysis and meets following Maxwell equations:
B=μ H (3)
On the basis of formula (1)-formula (3), two identities can be constructed, such as formula (4), (5) are shown.
When Three dimensional transient magnetic field is solved, the vector position free degree on its seamed edge employs single order unit and calculates Maxwell, and
The scalar potential free degree on node is calculated using second order unit.
In Three dimensional transient magnetic field, voltage source or current source can be called as model driving source, and winding is divided into two
Kind, one kind is twisted wire type winding, and one kind is entity winding.Wherein twisted wire winding does not consider Eddy Distribution, it is believed that the electricity in winding
Close is substantially uniformity, and entity winding then needs to calculate its Kelvin effect.At applied voltage source, the electric current on winding is not
Know, so the voltage for needing to carry out on winding loop is calculated.
For twisted wire type winding, because its resistance value is a lumped parameter, so when applied voltage source, Ke Yizhi
The numerical value that electric current is calculated by voltage source and D.C. resistance is connect, and for entity coil, because its resistance value and frequency, material
Etc. relevant, so the resistance is also referred to as AC resistance.Its alternating current resistance is not when to entity coil applied voltage source
Know, Maxwell is first to be calculated according to formula (6);
Vi=∫ ∫ ∫RiJ0i(E+v×B)dR (6)
Wherein, J0iIt is electric close on i-th loop, equally comprising the component on tri- directions of X, Y, Z.
If the magnetic field of change, then should also calculate winding back-emf, back-emf can be calculated according to formula (7).
Ei=∫ ∫ ∫RiHi·BidR (7)
In Three dimensional transient magnetic field, transient state include two implications, one be transient state electromagnetic process, one be transient state machine
Tool process.Because this calculating is not directed in mechanical transient process, therefore no longer processing mode is illustrated herein.
The analysis and calculating of electromagnetic field are generally summarized as non trivial solution of differentiating.For ODE, as long as by aiding in
After conditional decision arbitrary constant, its solution just turns into unique.For partial differential equation, unique subsidiary conditions can be turned into
Two kinds can be divided into, a kind of is the physical conditions residing for the border for expressing field, referred to as boundary condition;A kind is to determine the first of field
Beginning state, referred to as primary condition.Boundary condition and original condition are collectively referred to as definite condition.It is general to determine equation and definite condition as one
Individual entirety, referred to as Definite problem.Have in electromagnetic field solution procedure in various boundary conditions this solution procedurees, specifically should
With two categories below:
(1) Nuo Yiman boundary conditions
Wherein, τ is Nuo Yiman borders;N is the exterior normal vector of border Γ.
F (τ) and h (τ) be generic function, can be constant and zero, when for zero when claim homogeneous Nuo Yiman boundary conditions.Typically
Nuo Yiman boundary conditions are referred to as second kind boundary condition in electromagnetic problems, it defines the normal derivative point of boundary gesture
Cloth, the homogeneous Nuo Yiman boundary conditions mentioned in Maxwell calculating, i.e. normal derivative is zero, is default boundary condition, is not required to
User is wanted to specify.
(2) natural boundary conditions
Boundary condition on medium decomposition face, i.e., the tangential and Normal interface condition of different medium interface amounts belongs to certainly
Right boundary condition, is the default boundary condition of system in Maxwell calculating, it is not necessary to which user specifies.
In be considered in for actuating coil by whole Unified Model, the angle for project analysis is set out, and this analysis is to mould
Type is simplified, and simplified assumed condition has:
(1) local detail is simplified during whole model modeling;
(2) the indefinite place of model structure parameter is by empirical modeling;
(3) core material permeance BH curves are obtained with loss characteristic BP curves by picture information interpreting blueprints described point is provided
Arrive;
(4) the single-phase high pressure short circuit transient state excitation of transformer is provided by external circuit, and some parameters need to be adjusted according to calculating
It is whole;
(5) in addition to short-circuit Phase winding energization, remaining driving source is set to open circuit;
(6) because coil cake number is very more, the complexity of subsequent fluid calculating is considerably increased, so being extracted in result
When simplified summation process has been done to every mutually each grade winding.
It is to have influenceed that these are assumed or simplified on the accuracy analyzed, because the error that these factors are brought refuses meter
And.
The heavy difficult point solution of analysis process is as follows:
(1) electromagnetic fluid integration coupling
This project purpose be when there is the short circuit of single-phase high pressure in transformer, analysis on the influence of its interior oil stream, it is necessary to
The source condition that electromagnetic force suffered by each cake of each phase winding when being short-circuited is analysed as oily flow point.Each cake stress can be in fluid emulation
Indirect assignment in software Fluent, then carries out being calculated corresponding oil stream analysis result.But due to the structure of transformer
More complicated, actual condition during operation there may be than larger deviation with the electromagnetic force of direct addition, cause oily flow point to be analysed
Result it is inaccurate.
Now need using multiple physical field coupling function, using on electromagnetic simulation software Maxwell calculating transformer coils
Electric power, and mapped directly to Fluent softwares, calculate the influence to inside transformer Oil flow distribution.Needs are referred in particular to
Go out, when carrying out multiple physical field coupling based on ANSYS Workbench, the information such as all of model, simulation result is in difference
When being transmitted before software, be automatically obtained by the data mapping function of Workbench, and without again by other the 3rd
Square software, so as to ensure that the efficiency of work, and data transfer accuracy.
(2) model treatment
Because the appearance and size of transformer is larger, internal structure is complicated, and does not provide simulation model.Transformer whole model
Start from scratch foundation, during employ Maxwell and carry primitive to iron core and the modeling of high, medium and low and voltage regulation coil, significantly
Improve operating efficiency.
(3) pre-treatment such as material properties, working condition work
Involved material properties are calculated, such as iron core permeance, loss characteristic only provides picture, and interpreting blueprints work is not only
It is cumbersome and very high to required precision.The picture capture instrument that is carried using Maxwell of calculating, in directly reading electronic pictures
Curve, saves the manual described point of papery, enormously simplify intricate operation degree.But explanation is needed, this instrument is still operated with engineer, figure
Tablet quality, there is larger dependence, and interpreting blueprints work obtains the error that data are caused, and simulation calculation is not considered.
Short-circuit transient state excitation is provided by external circuit, and some parameters according to result of calculation depending on that need to be adjusted gained transient current
Winding is put on, with more real emulation actual conditions.
(4) mesh generation
Maxwell is that FInite Element is, by whole domain discretization, to be divided into many small based on finite element algorithm
Region, referred to as " unit " or " finite element ".Thus, the quality of mesh generation will directly affect the precision of result of calculation.And become
Depressor appearance and size is larger, and internal structure is complicated, and great challenge is brought to mesh generation.And the more difficult assurance of mesh generation work,
It is not very abundant engineer particularly with experience.
This is calculated the automatic adaptive mesh generation function of being possessed using Maxwell static fields and imports work(with dynamic grid
Can, static adaptive mesh is applied in transient magnetic field and is calculated.This method not only improves operating efficiency, and drops significantly
Low simulation accuracy improves operability credible with result of calculation to engineer and the dependence of software skilled operation degree
Degree.Fig. 9-Figure 12 cuts open result displaying for net.
(5) amount of calculation
Maxwell is that a quarter model facetization grid amount is in 930,000 or so, and Maxwell based on finite element algorithm
The three-dimensional grid subdivision unit for being used be the most stable tetrahedron two-order unit in finite element theory, i.e. each grid except
Using its 4 summits as outside node, while will also also serve as its calculate node per Tiao Leng centers, each grid has 10 sections
Point, as shown in figure 13, so computational accuracy is very high.
Further, since the main purpose that this electromagnetism is calculated is the analysis to each cake stress of each phase winding of transformer, so
On the basis of basic calculating amount, more than 15 Force Calculation matrix is with the addition of, it is clear that amount of calculation is very huge.
Using eHPC parallel computations, solved, calculating has used 10 cores to run simultaneously.Calculate and use 12 cores,
The server of 64G memory configurations is solved.
(6) result is extracted
Because transformer model is complicated, each mutually each grade winding cake number is more, is 140 multigroup per cake stress data, calculates
Amount is very more, and considerably increases the complexity during subsequent fluid is calculated.
Final determination result has done simplified summation process when extracting to every mutually each grade winding.
The single-phase high pressure short circuit transient state excitation of transformer is provided by external circuit, and some parameters need to be adjusted according to calculating,
Wherein, pressure regulation is played when middle pressure winding works, voltage regulation coil is not in working order.Operating mode is single-phase pressure short circuit in winding high, middle,
B phases are high, middle pressure short circuit in winding when high-tension coil voltage, electric current as shown in figure 14, it is 0 that B, C phase winding are set to open circuit i.e. electric current, outward
Circuit is built as shown in figure 15.
Emi analysis, the single-phase pressure high, middle of calculating transformer are carried out to the single-phase high pressure short circuit in winding transient process of transformer
During short circuit in winding, internal magnetic field distribution character, each phase winding and iron core ohmic loss, and each phase winding stressing conditions.
(1) excitation waveform
, when there is single-phase short circuit, high-pressure side A, B, C three-phase current transient waveform is as shown in figure 16, wherein waveform for transformer
It is short-circuit phase transient current.Short-circuit mutually high, medium voltage side current waveform is as shown in figure 17.
(2) field amount result
Whole result can show the Distribution of Magnetic Field of any moment transformer entire infrastructure, and Figure 18-Figure 21 is respectively different
During the short circuit of phase high pressure, synchronization inside transformer Magnetic Induction Density Distribution cloud atlas and polar plot.
Figure 22 is referred to, a kind of calculating of the 220kV three-phase transformers B phases short circuit electromagnetic field provided in the embodiment of the present invention
One embodiment of device includes:
Unit 301 is set up, for setting up 220kV three-phase transformer simplified models;
Subdivision unit 302, for being several three-dimensional grid subdivision lists by 220kV three-phase transformer simplified models subdivision
Unit;
Analytic unit 303, for being carried out at the short circuit in winding corresponding with B to three-dimensional grid subdivision unit several described
Reason, and driving source is set, carry out Electromagnetic Transient Analysis.
In the present embodiment, set up unit 301 and pass through to get the transformer basic-element model foundation by Maxwell instruments
220kV three-phase transformer simplified models, subdivision unit 302 is by the static adaptive mesh of Maxwell instruments by 220kV tri-
Phase transformer simplified model subdivision is several three-dimensional grid subdivision units, and analytic unit 303 is by Maxwell instruments to some
The individual three-dimensional grid subdivision unit carries out the short circuit in winding treatment corresponding with B, and sets driving source, carries out electromagnetic transient point
Analysis, solves because the model of current large-scale 220kV transformers is more complicated, and the nonlinear characteristic of material is difficult to obtain, and asks
The amount of calculation of solution than larger, caused by the low technical problem of large-scale 220kV transformers 3 D electromagnetic solution efficiency.
The above is to carry out detailed retouching to each unit of the computing device of electromagnetic field under three-phase transformer B phase short-circuit conditions
State, be below described in detail sub-unit, refer to Figure 23, a kind of 220kV three-phases provided in the embodiment of the present invention
Another embodiment of the computing device of transformer B phases short circuit electromagnetic field includes:
Unit 401 is set up, for setting up 220kV three-phase transformer simplified models;
Subdivision unit 402, for being several three-dimensional grid subdivision lists by 220kV three-phase transformer simplified models subdivision
Unit, subdivision unit 402, specifically for being several tetrahedron two-order units by 220kV three-phase transformer simplified models subdivision;
Analytic unit 403, for being carried out at the short circuit in winding corresponding with B to three-dimensional grid subdivision unit several described
Reason, and driving source is set, carry out Electromagnetic Transient Analysis.
Analytic unit 403 is specifically included:
Short circuit treatment subelement 4031, for carrying out the winding corresponding with B to tetrahedron two-order unit several described
Short circuit treatment, and external circuit short circuit transient state driving source is set;
Computation subunit 4032, by the seamed edge to tetrahedron two-order unit vector position the free degree carry out single order unit based on
Calculate, the scalar potential free degree to the node of tetrahedron two-order unit carries out the calculating of second order unit, solves Three dimensional transient magnetic field;
Analysis subelement 4033, for carrying out result of calculation analysis to Three dimensional transient magnetic field.
Analysis subelement 4033 is specifically included:
First analysis module 4033a, simplifies for being analyzed determination transformer according to the Three dimensional transient magnetic field for solving
The corresponding internal magnetic field distribution character of model, each phase winding and iron core ohmic loss, and each phase winding stressing conditions;
Second analysis module 4033b, for being analyzed according to the Three dimensional transient magnetic field for solving, gets short in B phases
The 220kV three-phases under the corresponding three-phase current transient waveform of 220kV three-phase transformers simplified model, preset time instant under line state
The corresponding Magnetic Induction Density Distribution cloud atlas of transformer simplified model and the magnetic induction intensity for obtaining is investigated under specified path position
Scatter chart.
In the present embodiment, set up unit 401 and pass through to get the transformer basic-element model foundation by Maxwell instruments
220kV three-phase transformer simplified models, subdivision unit 402 is by the static adaptive mesh of Maxwell instruments by 220kV tri-
Phase transformer simplified model subdivision is several three-dimensional grid subdivision units, and analytic unit 403 is by Maxwell instruments to some
The individual three-dimensional grid subdivision unit carries out the short circuit in winding treatment corresponding with B, and sets driving source, carries out electromagnetic transient point
Analysis, solves because the model of current large-scale 220kV transformers is more complicated, and the nonlinear characteristic of material is difficult to obtain, and asks
The amount of calculation of solution than larger, caused by the low technical problem of large-scale 220kV transformers 3 D electromagnetic solution efficiency.
It is apparent to those skilled in the art that, for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, may be referred to the corresponding process in preceding method embodiment, will not be repeated here.
In several embodiments provided herein, it should be understood that disclosed system, apparatus and method can be with
Realize by another way.For example, device embodiment described above is only schematical, for example, the unit
Divide, only a kind of division of logic function there can be other dividing mode when actually realizing, for example multiple units or component
Can combine or be desirably integrated into another system, or some features can be ignored, or do not perform.It is another, it is shown or
The coupling each other for discussing or direct-coupling or communication connection can be the indirect couplings of device or unit by some interfaces
Close or communicate to connect, can be electrical, mechanical or other forms.
The unit that is illustrated as separating component can be or may not be it is physically separate, it is aobvious as unit
The part for showing can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be according to the actual needs selected to realize the mesh of this embodiment scheme
's.
In addition, during each functional unit in each embodiment of the invention can be integrated in a processing unit, it is also possible to
It is that unit is individually physically present, it is also possible to which two or more units are integrated in a unit.Above-mentioned integrated list
Unit can both be realized in the form of hardware, it would however also be possible to employ the form of SFU software functional unit is realized.
If the integrated unit is to realize in the form of SFU software functional unit and as independent production marketing or use
When, can store in a computer read/write memory medium.Based on such understanding, technical scheme is substantially
The part for being contributed to prior art in other words or all or part of the technical scheme can be in the form of software products
Embody, the computer software product is stored in a storage medium, including some instructions are used to so that a computer
Equipment (can be personal computer, server, or network equipment etc.) performs the complete of each embodiment methods described of the invention
Portion or part steps.And foregoing storage medium includes:USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can store journey
The medium of sequence code.
The above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although with reference to preceding
Embodiment is stated to be described in detail the present invention, it will be understood by those within the art that:It still can be to preceding
State the technical scheme described in each embodiment to modify, or equivalent is carried out to which part technical characteristic;And these
Modification is replaced, and does not make the spirit and scope of the essence disengaging various embodiments of the present invention technical scheme of appropriate technical solution.
Claims (10)
1. computational methods of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field, it is characterised in that including:
Set up 220kV three-phase transformer simplified models;
It is several three-dimensional grid subdivision units by the 220kV three-phase transformers simplified model subdivision;
The short circuit in winding treatment corresponding with B is carried out to three-dimensional grid subdivision unit several described, and driving source is set, carried out
Electromagnetic Transient Analysis.
2. computational methods of 220kV three-phase transformers B phases according to claim 1 short circuit electromagnetic field, it is characterised in that will
The 220kV three-phase transformers simplified model subdivision is specifically included for several three-dimensional grid subdivision units:
It is several tetrahedron two-order units by the 220kV three-phase transformers simplified model subdivision.
3. computational methods of 220kV three-phase transformers B phases according to claim 2 short circuit electromagnetic field, it is characterised in that right
Several described three-dimensional grid subdivision units carry out the short circuit in winding treatment corresponding with B, and set driving source, carry out electromagnetism wink
State analysis is specifically included:
The short circuit in winding treatment corresponding with B is carried out to tetrahedron two-order unit several described, and external circuit short circuit wink is set
State driving source;
The vector position free degree to the seamed edge of the tetrahedron two-order unit carries out the calculating of single order unit, to the tetrahedron second order list
The scalar potential free degree of the node of unit carries out second order unit and calculates, and solves Three dimensional transient magnetic field;
Result of calculation analysis is carried out to the Three dimensional transient magnetic field.
4. computational methods of 220kV three-phase transformers B phases according to claim 3 short circuit electromagnetic field, it is characterised in that right
The Three dimensional transient magnetic field carries out result of calculation analysis and specifically includes:
It is corresponding that the Three dimensional transient magnetic field according to solving is analyzed the determination 220kV three-phase transformers simplified model
Internal magnetic field distribution character, each phase winding and iron core ohmic loss, and each phase winding stressing conditions.
5. computational methods of 220kV three-phase transformers B phases according to claim 4 short circuit electromagnetic field, it is characterised in that root
Being analyzed according to the Three dimensional transient magnetic field for solving also includes:
Get the corresponding three-phase current transient waveform of the 220kV three-phase transformers simplified model under B phase short-circuit conditions, pre-
The 220kV three-phase transformers simplified model corresponding Magnetic Induction Density Distribution cloud atlas is inscribed when putting and under specified path position
The Magnetic Induction Density Distribution curve map that investigation is obtained.
6. the computing device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field, it is characterised in that including:
Unit is set up, for setting up 220kV three-phase transformer simplified models;
Subdivision unit, for being several three-dimensional grid subdivision units by the 220kV three-phase transformers simplified model subdivision;
Analytic unit, for carrying out the short circuit in winding treatment corresponding with B to three-dimensional grid subdivision unit several described, and sets
Driving source is put, Electromagnetic Transient Analysis are carried out.
7. the computing device of 220kV three-phase transformers B phases according to claim 6 short circuit electromagnetic field, it is characterised in that cut open
Subdivision, specifically for being several tetrahedron two-order units by the 220kV three-phase transformers simplified model subdivision.
8. the computing device of 220kV three-phase transformers B phases according to claim 7 short circuit electromagnetic field, it is characterised in that point
Analysis unit includes:
Short circuit treatment subelement, for carrying out the short circuit in winding treatment corresponding with B to tetrahedron two-order unit several described,
And external circuit short circuit transient state driving source is set;
Computation subunit, the vector position free degree for the seamed edge to the tetrahedron two-order unit carries out the calculating of single order unit, right
The scalar potential free degree of the node of the tetrahedron two-order unit carries out second order unit and calculates, and solves Three dimensional transient magnetic field;
Analysis subelement, for carrying out result of calculation analysis to the Three dimensional transient magnetic field.
9. the computing device of 220kV three-phase transformers B phases according to claim 8 short circuit electromagnetic field, it is characterised in that point
Analysis subelement includes:
First analysis module, mould is simplified for being analyzed the determination transformer according to the Three dimensional transient magnetic field for solving
The corresponding internal magnetic field distribution character of type, each phase winding and iron core ohmic loss, and each phase winding stressing conditions.
10. according to the computing device of the 220kV three-phase transformer B phases short circuit electromagnetic field described in claim 9, it is characterised in that point
Analysis subelement also includes:
Second analysis module, for being analyzed according to the Three dimensional transient magnetic field for solving, gets in B phase short-circuit conditions
Under the phase transformations of the 220kV tri- under the corresponding three-phase current transient waveform of the 220kV three-phase transformers simplified model, preset time instant
The magnetic induction intensity point that the corresponding Magnetic Induction Density Distribution cloud atlas of depressor simplified model and being investigated under specified path position is obtained
Cloth curve map.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710037842.8A CN106815436A (en) | 2017-01-18 | 2017-01-18 | The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710037842.8A CN106815436A (en) | 2017-01-18 | 2017-01-18 | The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106815436A true CN106815436A (en) | 2017-06-09 |
Family
ID=59112400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710037842.8A Pending CN106815436A (en) | 2017-01-18 | 2017-01-18 | The computational methods and device of a kind of 220kV three-phase transformers B phases short circuit electromagnetic field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106815436A (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550435A (en) * | 2015-12-13 | 2016-05-04 | 东北电力大学 | Winding magnetic flux leakage and electrodynamic force calculation method based on transformer segmented model |
-
2017
- 2017-01-18 CN CN201710037842.8A patent/CN106815436A/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550435A (en) * | 2015-12-13 | 2016-05-04 | 东北电力大学 | Winding magnetic flux leakage and electrodynamic force calculation method based on transformer segmented model |
Non-Patent Citations (2)
Title |
---|
谢道伟: ""油浸式三相电力变压器电磁场及温度场的瞬态仿真分析"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
陈孝荣: ""自供电磁流变减振系统研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5692158A (en) | Methods for generating models of non-linear systems and components and for evaluating parameters in relation to such non-linear models | |
Annakkage et al. | Dynamic system equivalents: A survey of available techniques | |
CN106407499A (en) | Method and device for calculating three-dimensional electromagnetic fields of three-phase transformers | |
Liu et al. | A real-time nonlinear hysteretic power transformer transient model on FPGA | |
CN103018534B (en) | Determine the method and system of harmonic voltage | |
CN106055812A (en) | Three-dimensional transient electromagnetic field numerical simulation analysis method and device | |
CN101794321A (en) | Single-phase three-winding autotransformer model taking account of nonlinear influences of excitation impedance | |
CN104331544B (en) | A kind of three-phase full-bridge converter modeling method based on EIC principles | |
CN107423528B (en) | Method and device for analyzing short-circuit transient stress of power transformer winding | |
CN102609575A (en) | Power system transient stability simulating method based on implicit numerical integration | |
CN110705172A (en) | Method for calculating distribution of instantaneous electric field of transformer under lightning impulse voltage | |
CN110221237A (en) | A kind of capacitance type potential transformer harmonic voltage measurement method and apparatus | |
CN111724479B (en) | Three-dimensional modeling method and system for electrical primary equipment of transformer substation | |
CN106407576A (en) | An extracting method and system for electromagnetic radiation characteristic parameters of apparatuses | |
CN106650007A (en) | Computing method of leakage magnetic field and winding vibration displacement of transformer in steady state | |
CN108448631A (en) | The Dynamic Phasors modeling method that power distribution network is accessed containing distributed generation resource based on model reduction | |
CN106844950A (en) | The computational methods and device of a kind of 220kV three-phase transformers A phases short circuit electromagnetic field | |
Jurisic et al. | High frequency transformer model derived from limited information about the transformer geometry | |
Yao et al. | A highly efficient domain decomposition method applied to 3-D finite-element analysis of electromechanical and electric machine problems | |
CN109145447A (en) | A kind of controllable magnetic saturation reactor modeling method | |
CN106682362A (en) | Method and device for computing C-phase short-circuit electromagnetic field of 220kV three-phase transformer | |
CN106446372A (en) | Method and device for computing electromagnetic field of three-phase transformer under A-phase short-circuit state | |
CN104679936B (en) | Permanent magnetism saturation type fault current limiter magnetic current thermal coupling modeling method | |
CN103258124A (en) | Power frequency magnetic field assessment method for electromagnetic equipment of high voltage transformer substation | |
CN105243272A (en) | Transformer winding structure inductance calculation method |
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: 20170609 |
|
RJ01 | Rejection of invention patent application after publication |