CN105844069A - Temperature field calculation method and device of oil-immersed transformer - Google Patents
Temperature field calculation method and device of oil-immersed transformer Download PDFInfo
- Publication number
- CN105844069A CN105844069A CN201610463081.8A CN201610463081A CN105844069A CN 105844069 A CN105844069 A CN 105844069A CN 201610463081 A CN201610463081 A CN 201610463081A CN 105844069 A CN105844069 A CN 105844069A
- Authority
- CN
- China
- Prior art keywords
- model
- temperature
- transformator
- oil
- interface
- 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.)
- Granted
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
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)
- Housings And Mounting Of Transformers (AREA)
Abstract
The embodiments of the invention disclose a temperature field calculation method and device of an oil-immersed transformer. The method and the device are used for solving the technical problem of the existing technical schemes that the calculation accuracy is low due to the fact that merely overall loss of a coil or core is calculated in a calculation process instead of loss distribution. The temperature field calculation method of the oil-immersed transformer, disclosed by the embodiments of the invention, comprises the steps: simplifying a transformer model modeled through a first mode; carrying out grid division on the simplified transformer model, and defining corresponding boundaries and interfaces; merging divided grids through a second mode, and generating corresponding interfaces; acquiring a turbulence model, attribute parameters and conditional parameters determined through the second mode; coupling electromagnetic fields and temperature fields so as to determine temperature based total loss; carrying out simulation processing according to the turbulence model, the attribute parameters, the conditional parameters and the total loss as well as a determined solving mode and determined convergence criteria, thereby acquiring a corresponding temperature field distribution result.
Description
Technical field
The present invention relates to technical field of electric power, particularly relate to a kind of oil-filled transformer Calculation Method of Temperature Field
And device.
Background technology
Oil-filled transformer, as very important electrical equipment a kind of in power system, plays transmission electric energy
With the effect changing system voltage grade.It plays huge in reasonable distribution and the economic use of electric energy
Big effect, is therefore widely used in power engineering.The safety of power system is mainly by electric power
The safe operation of equipment determines, and transformator is one of most important high pressure equipment, so should ensure that it
Safe and reliable operation.During transformer production manufacture and product analysis, people are the most continuous
Explore and test the operational reliability improving transformator.
Transformator, during normal work, inevitably produces heat.If there is hot-spot,
Then may change electromagnetism and the hot property of transformer part, jeopardize the properly functioning of transformator, therefore transformation
Design and the exploitation of transformator are had very important significance by the control of device temperature rise.
The transformator in past, owing to capacity is little, thermoelectricity temperature rise is inconspicuous, does not cause too many attention, leads
Cause to its research seldom.Add complexity and the test difficulty of its problem self, hardly result in accurately number
Value, so the research to transformator heating problem is made slow progress always.The true eddy current of transformator eddy-current loss goods
Loss density is excessive is likely to cause hot-spot problem, and then causes some ill effects, and impact becomes
The safe operation of depressor.
In order to the heating of transformator is controlled, in transformator, typically there is employing cooling structure.According to
Statistics, the oil-filled transformer of current more than 80% uses the Natural Oil Circulation Power type of cooling.This cold
But the direct result of mode is that copper oil temperature difference strengthens, and is that the heat dissipation problem of transformator is more prominent.Tradition side
In method, the computing formula of transformator temperature rise is that the method by hot road and substantial amounts of experimental data are summed up and obtained
's.But due to the change of insulation system, it is no longer can meet mix insulation structure liquid soaking transformer
Design requirement.
Along with the development of existing CAE (computer-aided engineering) technology, transformator is directly emulated,
The temperature rise of transformator is obtained by the way of numerical computations.The problem of temperature rise of transformator relate to electromagnetism,
Multiple subject such as calorifics and hydrodynamics, is an extremely complex problem.Therefore in existing computing hardware
With under software, generally it is required for calculating temperature rise by certain simplification means.These simplification means obtain
The temperature rise of transformator, in certain scope of application, the acceptable during result obtained.But,
Simplify the calculating error caused, be sometimes difficult to estimate, result even in the result mistake of calculating, finally
Threaten the safety of transformator.
At present when carrying out the calculating of three-dimensional temperature field of transformator, it is broadly divided into two steps: first pass through reality
Test, experience or numerical computations, directly obtain the overall losses on coil, iron core;Then transformation is being calculated
During the three-dimensional temperature field of device, will loss as thermal source load, thus solve and obtain transformer temperature distribution.
Although existing technical scheme can improve the design efficiency of design of transformer teacher, but in the mistake calculated
Calculated in journey is only the overall losses on coil or iron core rather than loss distribution, thus leads
Cause the technical problem that computational accuracy is low.
Summary of the invention
A kind of oil-filled transformer Calculation Method of Temperature Field of embodiment of the present invention offer and device, solve
Existing technical scheme is calculated during calculating is only the overall damage on coil or iron core
Consumption rather than loss distribution, and the technical problem that the computational accuracy that causes is low.
A kind of oil-filled transformer Calculation Method of Temperature Field that the embodiment of the present invention provides, including:
The transformer model set up by first mode is carried out simplification process;
Transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding border and boundary
Face;
Merged the described grid divided by the second pattern, and generate the described interface of correspondence;
Get turbulence model, property parameters and the conditional parameter determined by described second pattern;
Carry out the coupling processing in electromagnetic field and temperature field, determine total losses based on temperature;
According to described turbulence model, described property parameters, described conditional parameter and described total losses, and tie
Close the Solution model determined and convergence criterion carries out simulation process, obtain corresponding thermo parameters method result.
Preferably, the transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding sides
Boundary and interface specifically include:
According to the attribute of computational fields or geometric properties or geometry, described transformator simplified model is divided yardstick enter
Row piecemeal processes, and generates a plurality of sub-block;
A plurality of described sub-blocks are carried out stress and strain model process, and defines corresponding border and interface.
Preferably, get turbulence model, property parameters and the condition determined by described second pattern to join
Also include before number:
Calculate the Inlet Reynolds Number of described transformator simplified model, determine stream in described transformator simplified model
Move as turbulent flow.
Preferably, described property parameters is material properties parameter, and described conditional parameter is boundary condition parameter.
Preferably, carry out the coupling processing in electromagnetic field and temperature field, determine that total losses based on temperature are concrete
Including:
Derive the geometric coordinate of thermal source in the fluid model corresponding with described transformator simplified model;
Importing described geometric coordinate in electromagnetic model, interpolation calculation goes out the power density values of respective coordinates;
Described power density values is led back in described fluid model, and is applied to the described geometric coordinate of correspondence
On, determine total losses based on temperature.
Preferably, according to described turbulence model, described property parameters, described conditional parameter and described total damage
Consumption, and combine the Solution model that determines and convergence criterion carries out simulation process, obtain corresponding temperature field and divide
Cloth result specifically includes:
The described convergence criterion of correspondence is determined according to described turbulence model;
Described transformator simplified model is carried out initialization process;
According to described turbulence model, described property parameters, described conditional parameter and described total losses, and tie
Close the Solution model determined and convergence criterion the described transformator simplified model after initializing is carried out at emulation
Reason, obtains corresponding thermo parameters method result.
A kind of oil-filled transformer Temperature calculating device that the embodiment of the present invention provides, including:
Simplify processing unit, for the transformer model set up by first mode is carried out simplification process;
Stress and strain model unit, the transformator simplified model after processing simplification carries out stress and strain model, and
The corresponding border of definition and interface;
Interface signal generating unit, for being merged the described grid divided by the second pattern, and generates correspondence
Described interface;
Acquiring unit, with getting turbulence model, property parameters and the bar determined by described second pattern
Part parameter;
Coupling processing unit, for carrying out the coupling processing in electromagnetic field and temperature field, determines based on temperature
Total losses;
Temperature calculating unit, for according to described turbulence model, described property parameters, described condition ginseng
Number and described total losses, and combine the Solution model that determines and convergence criterion carries out simulation process, it is right to obtain
The thermo parameters method result answered.
Preferably, stress and strain model unit specifically includes:
Partition subelement, for described transformator simplified model according to the attribute of computational fields or geometry
Feature or geometry divide yardstick and carry out piecemeal process, generate a plurality of sub-block;
Stress and strain model subelement, for a plurality of described sub-blocks being carried out stress and strain model process, and define right
Answer border and interface;
Coupling processing unit specifically includes:
Derive subelement, for deriving thermal source in the fluid model corresponding with described transformator simplified model
Geometric coordinate;
Power density computation subelement, for importing described geometric coordinate, interpolation calculation in electromagnetic model
Go out the power density values of respective coordinates;
Apply subelement, for described power density values is led back in described fluid model, and be applied to right
On the described geometric coordinate answered, determine total losses based on temperature.
Preferably, oil-filled transformer Temperature calculating device also includes:
Entrance computing unit, for calculating the Inlet Reynolds Number of described transformator simplified model, determines described
In transformator simplified model, flowing is turbulent flow.
Preferably, Temperature calculating unit specifically includes:
Convergence determines subelement, for determining the described convergence criterion of correspondence according to described turbulence model;
Initialize subelement, for described transformator simplified model is carried out initialization process;
Temperature calculating subelement, for according to described turbulence model, described property parameters, described condition
Parameter and described total losses, and combine the Solution model determined and convergence criterion to the described change after initializing
Depressor simplified model carries out simulation process, obtains corresponding thermo parameters method result;
Wherein, described property parameters is material properties parameter, and described conditional parameter is boundary condition parameter.
As can be seen from the above technical solutions, the embodiment of the present invention has the advantage that
A kind of oil-filled transformer Calculation Method of Temperature Field of embodiment of the present invention offer and device, wherein,
Oil-filled transformer Calculation Method of Temperature Field includes: carry out the transformer model set up by first mode
Simplification processes;Transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding border
And interface;Merged the grid divided by the second pattern, and generate the interface of correspondence;Get logical
Turbulence model, property parameters and the conditional parameter that second pattern of crossing determines;Carry out electromagnetic field and temperature field
Coupling processing, determines total losses based on temperature;According to turbulence model, property parameters, conditional parameter and
Total losses, and combine the Solution model that determines and convergence criterion carries out simulation process, obtain corresponding temperature
Field distribution result.In the present embodiment, carry out grid by the transformator simplified model after simplification is processed and draw
Point, and define corresponding border and interface;Merged the grid divided by the second pattern, and generate correspondence
Interface;Get turbulence model, property parameters and the conditional parameter determined by the second pattern;Enter
Row electromagnetic field and the coupling processing in temperature field, determine total losses based on temperature;According to turbulence model, genus
Property parameter, conditional parameter and total losses, and combine the Solution model determined and convergence criterion is carried out at emulation
Reason, obtains corresponding thermo parameters method result, solves existing technical scheme and counts during calculating
Obtain is only the overall losses on coil or iron core rather than loss distribution, and the calculating caused
The technical problem that precision is low.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to reality
Execute the required accompanying drawing used in example or description of the prior art to be briefly described, it should be apparent that below,
Accompanying drawing in description is only some embodiments of the present invention, for those of ordinary skill in the art,
On the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
One reality of a kind of oil-filled transformer Calculation Method of Temperature Field that Fig. 1 provides for the embodiment of the present invention
Execute the schematic flow sheet of example;
Another of a kind of oil-filled transformer Calculation Method of Temperature Field that Fig. 2 provides for the embodiment of the present invention
The schematic flow sheet of embodiment;
One reality of a kind of oil-filled transformer Temperature calculating device that Fig. 3 provides for the embodiment of the present invention
Execute the structural representation of example;
Another of a kind of oil-filled transformer Temperature calculating device that Fig. 4 provides for the embodiment of the present invention
The structural representation of embodiment;
Fig. 5 is the application examples schematic diagram of Fig. 2;
Fig. 6 is section speed vector figure;
Fig. 7 is section Temperature Distribution cloud atlas.
Detailed description of the invention
A kind of oil-filled transformer Calculation Method of Temperature Field of embodiment of the present invention offer and device, solve
Existing technical scheme is calculated during calculating is only the overall damage on coil or iron core
Consumption rather than loss distribution, and the technical problem that the computational accuracy that causes is low.
For making the goal of the invention of the present invention, feature, the advantage can be the most obvious and understandable, below will
In conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground describes, it is clear that the embodiments described below are only a part of embodiment of the present invention, and not all
Embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art are not making creativeness
The all other embodiments obtained under work premise, broadly fall into the scope of protection of the invention.
Refer to Fig. 1, the one of a kind of oil-filled transformer Calculation Method of Temperature Field that the embodiment of the present invention provides
Individual embodiment includes:
101, the transformer model set up by first mode is carried out simplification process;
In the present embodiment, when needs transformer temperature Flow Field Numerical simulation calculation, it is necessary first to by
The transformer model of one Model Establishment carries out simplification process.
102, the transformator simplified model after processing simplification carries out stress and strain model, and define corresponding border and
Interface;
After the transformer model set up by first mode is carried out simplification process, need to simplification
Transformator simplified model after reason carries out stress and strain model, and defines corresponding border and interface.
103, merged the grid divided by the second pattern, and generate the interface of correspondence;
Transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding border and friendship
After interface, need to merge, by the second pattern, the grid divided, and generate the interface of correspondence.
104, turbulence model, property parameters and the conditional parameter determined by the second pattern is got;
After being merged the grid divided by the second pattern, and generating the interface of correspondence, need to obtain
To the turbulence model determined by the second pattern, property parameters and conditional parameter.
105, carry out the coupling processing in electromagnetic field and temperature field, determine total losses based on temperature;
After getting turbulence model, property parameters and the conditional parameter determined by the second pattern, need
Electromagnetic field to be carried out and the coupling processing in temperature field, determine total losses based on temperature.
106, according to turbulence model, property parameters, conditional parameter and total losses, and combine determine solve
Pattern and convergence criterion carry out simulation process, obtain corresponding thermo parameters method result.
When carrying out the coupling processing of electromagnetic field and temperature field, after determining total losses based on temperature, need
According to turbulence model, property parameters, conditional parameter and total losses, and combine the Solution model and receipts determined
Hold back criterion and carry out simulation process, obtain corresponding thermo parameters method result.
In the present embodiment, carry out stress and strain model by the transformator simplified model after simplification is processed, and fixed
The corresponding border of justice and interface;Merged the grid divided by the second pattern, and generate the interface of correspondence;
Get turbulence model, property parameters and the conditional parameter determined by the second pattern;Carry out electromagnetic field and
The coupling processing in temperature field, determines total losses based on temperature;According to turbulence model, property parameters, bar
Part parameter and total losses, and combine the Solution model that determines and convergence criterion carries out simulation process, it is right to obtain
The thermo parameters method result answered, solve existing technical scheme calculate during calculated only
It is only the overall losses on coil or iron core rather than loss distribution, and the skill that the computational accuracy that causes is low
Art problem.
The above is that the process to oil-filled transformer Calculation Method of Temperature Field is described in detail, and below will
Being described in detail detailed process, refer to Fig. 2, a kind of oil immersed type that the embodiment of the present invention provides becomes
Another embodiment of depressor Calculation Method of Temperature Field includes:
201, the transformer model set up by first mode is carried out simplification process;
In the present embodiment, when needs transformer temperature Flow Field Numerical simulation calculation, it is necessary first to by
The transformer model of one Model Establishment carries out simplification process.
202, according to the attribute of computational fields or geometric properties or geometry, transformator simplified model is divided yardstick to enter
Row piecemeal processes, and generates a plurality of sub-block;
After the transformer model set up by first mode is carried out simplification process, need transformator
Simplified model divides yardstick according to the attribute of computational fields or geometric properties or geometry and carries out piecemeal process, generates
A plurality of sub-blocks.
203, a plurality of sub-blocks are carried out stress and strain model process, and define corresponding border and interface;
Carry out when transformator simplified model is divided yardstick according to the attribute of computational fields or geometric properties or geometry
Piecemeal processes, and after generating a plurality of sub-block, needs a plurality of sub-blocks are carried out stress and strain model process, and
The corresponding border of definition and interface.
204, merged the grid divided by the second pattern, and generate the interface of correspondence;
When a plurality of sub-blocks being carried out stress and strain model process, and after defining corresponding border and interface, need
To merge, by the second pattern, the grid divided, and to generate the interface of correspondence.
205, the Inlet Reynolds Number of calculating transformer simplified model, in determining transformator simplified model, flowing is
Turbulent flow;
After being merged the grid divided by the second pattern, and generating the interface of correspondence, need to calculate
The Inlet Reynolds Number of transformator simplified model, flows for turbulent flow in determining transformator simplified model.
206, turbulence model, property parameters and the conditional parameter determined by the second pattern is got;
When the Inlet Reynolds Number of calculating transformer simplified model, in determining transformator simplified model, flowing is for rapid
After stream, need to get turbulence model, property parameters and the conditional parameter determined by the second pattern.
Property parameters is material properties parameter, and conditional parameter is boundary condition parameter.
207, the geometric coordinate of thermal source in the fluid model corresponding with transformator simplified model is derived;
After getting turbulence model, property parameters and the conditional parameter determined by the second pattern, need
Derive the geometric coordinate of thermal source in the fluid model corresponding with transformator simplified model.
208, importing geometric coordinate in electromagnetic model, interpolation calculation goes out the power density values of respective coordinates;
When deriving in the fluid model corresponding with transformator simplified model after the geometric coordinate of thermal source, need
Importing geometric coordinate in electromagnetic model, interpolation calculation goes out the power density values of respective coordinates.
209, power density values is led back in fluid model, and be applied on the geometric coordinate of correspondence, determine
Total losses based on temperature;
When importing geometric coordinate in electromagnetic model, after interpolation calculation goes out the power density values of respective coordinates,
Need to lead back in fluid model power density values, and be applied to correspondence geometric coordinate on, determine based on
The total losses of temperature.
210, the convergence criterion of correspondence is determined according to turbulence model;
When power density values is led back in fluid model, and it is applied on the geometric coordinate of correspondence, determines base
After the total losses of temperature, need to determine the convergence criterion of correspondence according to turbulence model.
211, transformator simplified model is carried out initialization process;
After determine the convergence criterion of correspondence according to turbulence model, need transformator simplified model is carried out
Initialization processes.
212, according to turbulence model, property parameters, conditional parameter and total losses, and combine determine solve
Pattern and convergence criterion carry out simulation process to the transformator simplified model after initializing, and obtain corresponding temperature
Degree field distribution result.
After transformator simplified model is carried out initialization process, need according to turbulence model, attribute ginseng
Number, conditional parameter and total losses, and combine the Solution model determined and convergence criterion to the change after initializing
Depressor simplified model carries out simulation process, obtains corresponding thermo parameters method result.
It is described in detail, as it is shown in figure 5, application examples includes with a concrete application scenarios below:
(1) model simplification and process
The Temperature calculating of transformator, had both related to the flowing of transformer oil, also relate to heat at coil and
Transmission between other structure members, Gu belong to stream-Conjugate Heat Transfer analysis.Calculate at transformator flow field simulation
In, need, on the premise of ensureing computational accuracy, model is done some and reasonably simplify, thus control emulation
The calculating scale analyzed, simplify object chooses the influence degree being based on concrete structure to Temperature calculating
Determine.
Oil-filled transformer mainly includes copper cash, stay, insulating cylinder, pad, oil baffle, iron core, metal
Folder, connector and tank envelope etc., main pyrotoxin be three copper coils, metal clip, iron core and
Tank envelope.Three coils are rule part, and its heating power accounts for more than the 95% of total heating power, therefore
The radiating effect of coiler part is most important to the heat radiation of whole transformator, models by truth as far as possible.
Simultaneously in order to ensure the symmetry of model, in the outlet of fuel tank two ends, face need to move to plane of symmetry center.
The simplification of table 1 physical model
For reducing calculating scale, this calculates and only considers the 1/2 of real transformer model, the most only selects one
The model of side carries out simulation analysis.
When block mold is processed, it is necessary first to according to the difference of material and purposes, all solids is carried out
Piecemeal, the solid of different materials attribute and different purposes wants independent adult, such as copper cash and insulating cylinder, one
It is to have thermal source, a not heating, therefore want independent adult;Secondly, material and heating attribute basic
The region caused, can merge, such as tank envelope and folder;Finally, fluid domain and domain complementarity,
After the computational fields defining entirety and solid domain, fluid domain also determines.
In order to control the stress and strain model scale of fluid domain, on the premise of not affecting computational accuracy, need to be to solid
Body territory carries out part process.Including by material piecemeal, remove chamfering, local alignment, shape simplification etc..
(2) stress and strain model
Right existing model is the 1/2 of real transformer block mold, but this 1/2 model is carried out integral net
Lattice divide, and remain extremely difficult, so the method that multi-block technique divides should be used.First to existing mould
Type carries out piecemeal, is respectively divided grid, and while grid division, on the one hand defines corresponding border,
Ready for next arranging boundary condition;On the one hand define the interface of piecemeal, i.e. interface,
To ensure after importing in FLUENT, interface can be set, complete interface between irregular net
Definition.
The piecemeal principle of block mold mainly divides three classes, and the first kind is that the attribute according to computational fields carrys out piecemeal, as
Computational fields exists stagnant zone and rotary area, it is possible to according to the kinetic characteristic in territory simultaneously, be divided into quiet
Territory and two pieces of dynamic territory, be respectively divided grid;Equations of The Second Kind is that the geometry complexity according to computational fields carrys out piecemeal,
If a region exists multiple irregular geometry, these geometry can be divided into different blocks,
It is respectively divided in grid, or region and there is the irregular geometry of class, and remaining part is all regular,
Just therefore region can be divided into Rule section and irregularities, process respectively;3rd class is basis
The division yardstick of the geometry of computational fields carrys out piecemeal, as physical dimension compares the situation that geometry around is much smaller
Under, typically can process by independent piecemeal.
The piecemeal of this model processes based on second principle, i.e. model is come by its rule degree
Piecemeal.Wherein, though the quantity such as copper cash, stay, insulating cylinder and pad are various, but regular shape, contact
Relation is simple and clear, and has periodic characteristics in circumferential direction, is therefore attributed to a class;And tank envelope,
The component shape such as folder, core insulation is loaded down with trivial details and contact relation is complicated, is therefore attributed to another kind of.
The instrument of this stress and strain model uses ICME CFD, based on to the control of number of grid and Flow Field Calculation
The raising of precision, uses hexahedral mesh to divide for relatively simple part region as far as possible, and for
Generate the region that hexahedral mesh is relatively difficult, then use tetrahedral grid.In ICEM CFD, can
To use the division methods of BLOCK that model is carried out cutting whole into sections, generate hexahedral mesh.
When carrying out the operation of concrete stress and strain model, should be noted following Railway Project:
1) mutually internally inconsistent due to the geometry of three coils, if using BLOCK method to it
Carry out universal formulation, the most scrappy and bad control of BLOCK block, therefore this grid can be caused
Divide and use piecemeal to divide, coil is divided into three different sub-blocks divide respectively.Respectively
Connect three coils is the insulating cylinder that two thickness is less, owing to insulating cylinder surface contact relation is multiple
Miscellaneous, therefore select the middle face interface as piecemeal of insulating cylinder, loop construction is just by two insulation
The median plane of cylinder has been divided into three sub-blocks.Additionally, due to the contacting of coil two ends and other structural member
Relation is the most complicated, therefore respectively clips a part of region at the upper and lower two ends of coil and uses tetrahedral grid
Dividing mode, remaining part uses BLOCK method to carry out hexahedral mesh division.
2) copper cash in coil and insulating element about are except regular geometry and contact relation
Outside Jian Dan, the most also there is cyclophysis, therefore coiler part is being carried out hexahedron
During stress and strain model, this characteristic can be made full use of.First with the BLOCK side in ICEM CFD
Method is syncopated as the block in one of them cycle, is defining each piecemeal title and grid point distribution
After, use swivel replication method, obtain the BLOCK structure in whole coil period, then repair
Positive association relation, obtains the hexahedral mesh of coiler part, then generates mesh.
3) when selecting interface, the convenience of stress and strain model to be taken into full account and economy, to four
On the premise of face volume mesh division work is without impact, hexahedral mesh is farthest facilitated to divide.
When dividing interface, the thing phase of two side areas should be considered, if interface both sides is arbitrary
While be solid, then need on interface, arrange temperature coupling (coupled wall);If both sides
It is fluid, the most directly bonds;If the existing solid domain in interface both sides, there is again fluid domain,
It is accomplished by interface, by different region bursts, bonding respectively.At this computation model
In, the solid of interface both sides is all consistent with fluid mass physical property, therefore has only to according to material
Material attribute difference carries out burst to interface, does adhesion process the most again.
4) method of interface can make the model energy piecemeal of complexity carry out stress and strain model, but to interface
On grid have certain requirement, owing to overall region is divided into multiple subdomain by interface, and sub
Grid between territory is by realizing connecting at the enterprising row interpolation of interface.Therefore, both sides are having a common boundary
Mesh scale on face cannot differ too big, and magnitude is relatively.
5) also need to define region and border while grid division, so import at grid
In FLUENT, these regions and border can automatically generate.
(3) boundary condition
3.1 turbulence model
Calculating for this, be single fluid in transformator, transformer oil enters into little fuel tank from entrance, after
Enter rectification hole rectification through little tank bottoms, then flow into the runner of coil, subsequently from coil other end stream
Going out to enter large tank, then the outlet from two ends is flowed out.Transformer model only one of which owing to calculating enters
Mouthful, Inlet Reynolds Number can be calculated according to its parameter.
May determine that according to Inlet Reynolds Number, in transformator, flowing is turbulent flow.In FLUENT, it is available for choosing
The turbulence model selected has a lot, is most commonly used that k-ε model.Standard k-ε model has higher stability,
Economy and computational accuracy, have wide range of applications, and is suitable for the turbulent flow of high reynolds number, but is not suitable for eddy flow etc.
The flowing that anisotropy is stronger;RNG k-ε model is applicable to calculate low reynolds number turbulent flow, than master pattern more
Consider rotation effect well, strong eddy flow is moved computational accuracy and increases;Realizable k-ε model is relatively
The advantage of first two k-ε model is that Realizable k-ε model has continuity constraint for turbulent viscosity, can
To keep eddy stress consistent with true turbulent flow, it is possible to simulation plane and the diffusion of Circular Jet more accurately
Speed, calculates at eddy flow, to the CALCULATION OF BOUNDARY-LAYER of band direction pressure gradient with separate in the problems such as stream calculation,
Result of calculation more meets truth, simultaneously in the Complex Flows separating stream calculation and band Secondary Flow calculates
Also do well.Calculate for this, owing to flowing itself belongs in low speeds flow, and coil runner
Flowing velocity is less, the diffusion of entrance flowing Circular Jet to be considered, k-ε model and other model phase in addition
Ratio, the consumption to calculating resource is less, is also easier to convergence.Therefore considering, this calculates selection
Realizable k-ε model is more rational.
Additionally, due in coil bulk flow road, fluid-flow rate is less, directly use turbulent model not
Too suitable, therefore, need the when that this subregion being calculated to process as laminar flow zone,
Laminar zone option is selected when the cell zone condition of coil inside runner is defined.
The generation of 3.2 interfaces
Based on before the piecemeal of transformer model being processed, it is now desired to merge the net of piecemeal in FLUENT
Lattice, the multi-block technique integrating documents that will individually generate is a file.After importing all grids,
FLUENT additionally provides the function of grid optimization, and tetrahedral grid can be converted into polyhedron grid,
While ensureing mesh quality, reduce number of grid as far as possible, interface both sides net can also be reduced simultaneously
Lattice different scale.
After mesh update and conversion complete, it is possible to carry out interface setting in FLUENT.Friendship is set
Interface can ensure that the grid of different piecemeal is connected with each other, due to the grid on interface both sides in computational fields
Will not be just the same, therefore calculation procedure can contact the grid on both sides by the way of interpolation, makes flux believe
Breath can mutually transmit in interface both sides.When arranging interface, if having side is solid, then
Must be provided with interface is Coupled Wall, to ensure both sides temperature coupling;If both sides are fluid,
Then need not this arrange.
3.3 initial conditions
In the calculating in temperature field, operation to be performed is needed to include opening energy equation, material behavior, thermal resistance
Loss, the inlet flow rate of fluid.
1) owing to computational fields relating to the calculating in temperature field, so needing to open energy equation.
2) all material related in model is set up, as the common used materials such as water and air can be from FLUENT
Reading in the material depot carried, other self-defining material then needs defined in FLUENT.Temperature field
The material properties that relates to of calculating include density, specific heat capacity and pyroconductivity, for fluid, except this
Outside three parameters, go back toughness, and actually material properties and temperature is the most relevant, does not the most consider material
The temperature characterisitic of material, is because in the range of temperature of this computation model, the material properties of fluid
Change less, be therefore set to constant.The attribute of solid then takes the attribute under room temperature.
This material property used in calculating is as shown in the table:
Table 2 computational fields material property table
3) for oil tank system, needing to be respectively provided with entrance and exit boundary condition, porch here needs
Speed entrance to be arranged, exit then is respectively necessary for arranging pressure export according to service condition difference.
Meanwhile, this calculates the impact in view of gravity, therefore must add action of gravity option.This calculates
In, gravity direction, in Z-direction, therefore needs to apply the acceleration of gravity of a-Z direction.
4), while oil is cold, model have also contemplated that the cross-ventilation effect with surrounding, and it neutralizes the external world
Air contact is tank envelope, and ambient temperature is set to 18 DEG C, the convection transfer rate root of tank envelope
It is taken as 7W/m2K according to experience.
5) after completing all settings, grid can be checked, determining such as boundary condition and
It is the most reasonable that interface etc. are arranged.
In simulation calculation, the boundary condition that be there is a need to arrange is as shown in the table.
Table 3 boundary condition
3.4 electromagnetic fields couple with temperature field
Electromagnetism and the coupling of fluid, be mainly reflected in transformer temperature field calculate in thermal source applying mode on.
The effect in magnetic field in work, generally, the loss of thermal source is come from owing to the thermal source of transformator is lost
Distribution is not the most uniform, but changes according to the change of geometric position.In order to enable more preferably topotype
Intend distribution is lost really, it is considered to directly the loss calculation result of electromagnetism be applied on the thermal source of fluid,
Ensure that the calculating of transformer temperature field is close to truth by both couplings.
Owing to the calculating of electromagnetism is to carry out in MAXWELL, and the calculating in temperature field is at FLUENT
In carry out, therefore the coupling of electromagnetism and fluid mean that between two software data transmission.And it is this
There are two problems when the data of type are transmitted, on the one hand the geometry of two softwares is the most consistent, because electric
The calculating of magnetic and fluid is different to the requirement of model, and different software can for same geometric model
Different simplification requirements can be had, it is contemplated that the core of electromagnetism and fluid coupling is thermal source part,
Therefore coupling successful premise is that in two softwares, the geometry of thermal source part should keep consistent;On the other hand,
MAXWELL uses Finite Element Method, and FLUENT uses finite volume method, both
Method is required for carrying out computational fields the division of grid, and electromagnetism and fluid calculate the requirement to grid also
Differing, therefore in the case of geometry is completely the same, the grid that both divide still has bigger difference,
Here the method using interpolation realizes the data transmission on different grid node, be reduced as far as because
The inconsistent error brought of grid.
For the geometry requirement of coupling, this calculates the model of application in FLUENT is physical model
1/2, the model used in electromagnetism calculates is then full mould, although keeping consistent on geometric position, but
It is that model scope size is not consistent, it is difficult to carry out direct-coupling calculating, it is therefore desirable to calculated by electromagnetism
To thermal source interpolation be distributed in FLUENT.For this calculate, the basic ideas of coupling be divided into
Under several steps:
1) geometric coordinate of thermal source in fluid model is derived;
2) importing above-mentioned coordinate in electromagnetic model, interpolation calculation goes out the power density values of respective coordinates;
3) power density values is led back in fluid model, and be applied in respective coordinates
Owing to two kinds of softwares import inconsistent with data form everywhere, also need to lead in concrete operations engineering
Cross the mutual conversion that coding carries out between two kinds of data.
The concrete operation step of data coupling is as follows:
1) in FLUENT, derive the geometric coordinate of thermal source
In FLUENT, it is selected from the menu File-> Interpolate..., can open Interpolate data
Dialog box, selects to need the heat source region of output in Cell Zones, clicks the output of Write... button i.e.
Can, owing to program needs, export file name should be fixed as test1.ip.
2) the ip file that FLUENT generates is converted to MAXWELL discernible pts file,
Converse routine is realized by C language, and function is to change the form of data file, but does not revises
Concrete data, the filename of Program Generating is points.pts.
3) in MAXWELL, import point.pts, derive the loss value of respective nodes.Menu selects
Select Maxwell 3D-> Fields-> Calculator, open Fields Calculator dialog box,
Quantity drop-down list under Input selects Total-Loss (or OhmicLoss), point
Hitting Export, the filename of output is fixed as maxwell_res.fld.
4) the loss file that MAXWELL generates is converted back FLUENT discernible ip file,
Owing to ip file needs the geometric data of thermal source, file and the test1.ip that before generate therefore are lost
File should be placed in same file folder, and the c program of conversion is also placed in this document folder, program
Can automatically read the value of loss file, and combine the geometric data in test1.ip, generate one newly
The file of fixing entitled source.ip, this document contains the final loss value needing and applying.
5) in FLUENT model, add UDS, these UDS and can accommodate the damage imported by outside
Consumption value, clicks Define-> User-Defined-> Scalars... in a menu, opens User-Defined
Scalars dialog box, such as the most undefined UDS, then needing to increase UDS quantity is 1,
Solution Zones selects all zones or all solid zones, divides for corresponding heat source region
Join UDS.
6) in FLUENT, import source.ip file, open Interpolate data dialog box equally,
In Cell Zones, select the heat source region of correspondence, click Read button and select respective file
?.
7) in FLUENT, compile my_source.c, i.e. apply the UDF of thermal source.Menu selects
Define-> User-Defined-> Functions-> Compiled, opens Compiled UDFs dialogue
Frame, adds my_source.c in Source Files, clicks Build button, compile by rear,
Click Load button.
8) UDF of applying thermal source is added in corresponding heat source region, after UDF has compiled,
The Source Terms of corresponding thermal source, selects Energy, will have more in drop-down list on the right
One udf my_source::libudf, chooses this, clicks OK and completes this setting.
After coupling, it is impossible to re-start the initialization of model, because initialization can empty existing interior
Data relevant in depositing, it is therefore proposed that after initialization is complete, then carry out the coupling behaviour of electromagnetism and fluid
Make.It addition, for compiling the UDF passed through under monokaryon state, cannot directly use under many nuclear states,
Need to recompilate.
After completing coupling, additionally write UDF in order to verify the total losses value being applied on thermal source.Result
As shown in the table.Wherein fuel tank, folder and silicon steel are consistent with average operating mode, are uniform thermals source.And three
Individual coil is to apply thermal source by the method for coupling, and concrete loss value is as shown in the table.
Table 4 each computational fields heating power is distributed
Computational fields title | Computational fields thermal power (kw) | Percentage ratio (%) |
High-tension coil (Xianquan1) | 95.181 | 63.78 |
Intermediate voltage winding (Xianquan2) | 43.427 | 29.10 |
Low-voltage coil (Xiangquan3) | 5.336 | 3.57 |
Silicon steel | 1.300 | 0.87 |
Folder and fuel tank | 3.999 | 2.68 |
3.5 solve setting
After completing the arranging of boundary condition, need solver is configured, below figure institute is specifically set
Show.Wherein most solves parameter and keeps default value, when arranging solver, should be noted discrete scheme
Select.General flow problem will not directly use second order accuracy when calculating, but first uses single order precision
Calculate the field at an initial stage, restrain because single order form is easier than Second-order Up-wind form, but precision is poor.
After checking the result of single order precision, second order form can be used instead on the premise of guaranteeing the convergence of single order precision
Improve the precision of calculating.For the calculating in temperature field, typically after model has been set up, first calculate
Information of flow, and check flow field result, if Flow Field Calculation convergence, when result also compares reasonable simultaneously,
Just temperature field is incorporated to model calculate together.
3.6 convergence criterions and initialization
The amount of convergence of transformer model includes equation of continuity, the equation of momentum and energy equation respectively, due to this
In have employed k-ε model, therefore convergence criterion includes the convergence of k and ε, and k represents Turbulent Kinetic, ε
Represent turbulence dissipation rate.The selection of general convergency value rule of thumb determines, for general considerations, physics
The convergence criterion of amount is 0.001.During Practical Calculation, can suitably turn down residual error convergency value, with
Ask and obtain the higher simulation result of precision.The higher model of required precision is calculated, cannot be complete
The when of meeting convergence criterion, can select to monitor the temperature of concern part, if temperature is in the mistake of iteration
Journey varies less, then it is believed that calculate and restrained.When example convergence encounters difficulties, it is also possible to
Its convergence is improved by adjusting relaxation factor.
Before starting to calculate, need model is initialized.The convergence initializing stream field in flow field
Process can have a certain impact, and different initialization modes therefore can be selected to calculate respectively, and by receiving
Hold back curve to determine the selection of initialization mode.For this model, the initialization side of selection standard here
Method.After completing to initialize, the maximum step number of setup algorithm, can start to calculate.
3.7 simulation result
Owing in transformator, computational fields is more, and flow behavior complexity differs, and therefore checks result
Time, by checking that the feature profile information chosen flows and profiling temperatures in analyzing system.
3.8 oil Flow Velocity distributions
Before fluid enters each coil runner, fluid carries out rectification, so in first rectified hole after entering fuel tank
After flow through lower yoke insulation after enter in coil runner, therefore fluid divides in the flowing of lower yoke insulation position
Cloth directly affects flow and the flow velocity of each coil runner.
From fig. 6 it can be seen that flow velocity is the fastest in high-tension coil, and flow direction and geometric direction
Unanimously;In intermediate voltage winding, secondly, this coil inner flow passage is affected flow velocity by oil baffle, is distributed in " Z " font,
In runner direction with to flow consistent region flow velocity higher;Low-voltage coil always flows end stream owing to not having fluid
Entering, only contact with other regional fluids at end, near insulating cylinder position of opening, fluid velocity is not 0,
There is certain flow velocity, so can speed up and take away the heat distributed in low-voltage coil.
3.9 Temperature Distribution
It can be seen that high-tension coil temperature is the highest from section Temperature Distribution Fig. 7, intermediate voltage winding secondly,
Low-voltage coil is minimum.Owing to this thermal source loaded is ansoft coupling thermal source loss, the unit in coil
Volume heating power is spatially distributed inconsistent, and therefore, in each coil, coil temperature also differs
Causing, wherein Temperature Distribution difference most notably intermediate voltage winding, this also reflects coil in reality from side
The loading condition of regional in work.
3.10 with comparison of test results
The simulation value of the simulation value and temperature of choosing inlet pressure respectively contrasts with test, result such as following table
Shown in, its error, all within 20%, meets computational accuracy requirement.
Table 5 inlet pressure simulation value and experiment value comparing result
Numbering | Experiment value (Pa) | Simulation value (Pa) | Error (%) |
1 | 3.70E+04 | 3.38E+04 | -8.65 |
2 | 3.80E+04 | 3.38E+04 | -11.05 |
3 | 3.75E+04 | 3.38E+04 | -9.87 |
All coils temperature value takes volume mean temperature, for top-oil temperature, selects more to connect from fuel tank top
Near plane, as the plane of reference, takes the liquid mass flow mean temperature in this plane and carries out data analysis ratio
Right.
Table 6 transformer temperature simulation value contrasts with experiment value
Position | Experiment value (DEG C) | Simulation value (DEG C) | Error (%) |
Top-oil temperature | 44.5 | 46.3 | 3.96 |
High-tension coil | 69.9 | 74.5 | +6.58 |
Intermediate voltage winding | 70.3 | 63.3 | -9.95 |
Low-voltage coil | 49.0 | 57.7 | +17.76 |
In the present embodiment, carry out stress and strain model by the transformator simplified model after simplification is processed, and fixed
The corresponding border of justice and interface;Merged the grid divided by the second pattern, and generate the interface of correspondence;
Get turbulence model, property parameters and the conditional parameter determined by the second pattern;Carry out electromagnetic field and
The coupling processing in temperature field, determines total losses based on temperature;According to turbulence model, property parameters, bar
Part parameter and total losses, and combine the Solution model that determines and convergence criterion carries out simulation process, it is right to obtain
The thermo parameters method result answered, solve existing technical scheme calculate during calculated only
It is only the overall losses on coil or iron core rather than loss distribution, and the skill that the computational accuracy that causes is low
Art problem.
Further, by using electromagnetic field and the coupling analysis in temperature field, transformator work can precisely be obtained
Make the Temperature Distribution under state, go the question meanings such as hot-spot great for prevention transformation;The electricity provided
Magnetic field and the coupling process in temperature field, the problem that can effectively solve to be similar to, it is not limited to transformator
Analysis of heat transmission.
Refer to Fig. 3, a kind of oil-filled transformer Temperature calculating device provided in the embodiment of the present invention
One embodiment includes:
Simplify processing unit 301, for the transformer model set up by first mode is carried out at simplification
Reason;
Stress and strain model unit 302, the transformator simplified model after processing simplification carries out stress and strain model,
And define corresponding border and interface;
Interface signal generating unit 303, for being merged the grid divided by the second pattern, and generates correspondence
Interface;
Acquiring unit 304, with getting turbulence model, property parameters and the condition determined by the second pattern
Parameter;
Coupling processing unit 305, for carrying out the coupling processing in electromagnetic field and temperature field, determines based on temperature
Total losses;
Temperature calculating unit 306, is used for according to turbulence model, property parameters, conditional parameter and total losses,
And combine the Solution model that determines and convergence criterion carries out simulation process, obtain corresponding thermo parameters method knot
Really.
In the present embodiment, the transformator simplified model after being processed simplification by stress and strain model unit 302 is entered
Row stress and strain model, and define corresponding border and interface;Interface signal generating unit 303 is by the second pattern
Merge the grid divided, and generate the interface of correspondence;Acquiring unit 304 gets by the second pattern
Turbulence model, property parameters and the conditional parameter determined;Coupling processing unit 305 carries out electromagnetic field and temperature
The coupling processing of degree field, determines total losses based on temperature;Temperature calculating unit 306 is according to turbulent flow mould
Type, property parameters, conditional parameter and total losses, and combine the Solution model that determines and convergence criterion is carried out
Simulation process, obtains corresponding thermo parameters method result, solves existing technical scheme in the mistake calculated
Calculated in journey is only the overall losses on coil or iron core rather than loss distribution, and causes
The low technical problem of computational accuracy.
The above is that each unit to oil-filled transformer Temperature calculating device is described in detail, below
Each subelement will be described in detail, refer to Fig. 4, a kind of oil immersion provided in the embodiment of the present invention
Formula transformer temperature field calculates another embodiment of device and includes:
Simplify processing unit 401, for the transformer model set up by first mode is carried out at simplification
Reason;
Stress and strain model unit 402, the transformator simplified model after processing simplification carries out stress and strain model,
And define corresponding border and interface;
Stress and strain model unit 402 specifically includes:
Partition subelement 4021, for transformator simplified model according to the attribute of computational fields or geometry
Feature or geometry divide yardstick and carry out piecemeal process, generate a plurality of sub-block;
Stress and strain model subelement 4022, for a plurality of sub-blocks being carried out stress and strain model process, and define right
Answer border and interface;
Interface signal generating unit 403, for being merged the grid divided by the second pattern, and generates correspondence
Interface;
Entrance computing unit 404, for the Inlet Reynolds Number of calculating transformer simplified model, determines transformator
In simplified model, flowing is turbulent flow.
Acquiring unit 405, with getting turbulence model, property parameters and the condition determined by the second pattern
Parameter;
Coupling processing unit 406, for carrying out the coupling processing in electromagnetic field and temperature field, determines based on temperature
Total losses;
Coupling processing unit 406 specifically includes:
Derive subelement 4061, for deriving thermal source in the fluid model corresponding with transformator simplified model
Geometric coordinate;
Power density computation subelement 4062, for importing geometric coordinate, interpolation calculation in electromagnetic model
Go out the power density values of respective coordinates;
Apply subelement 4063, for power density values being led back in fluid model, and be applied to correspondence
On geometric coordinate, determine total losses based on temperature.
Temperature calculating unit 407, is used for according to turbulence model, property parameters, conditional parameter and total losses,
And combine the Solution model that determines and convergence criterion carries out simulation process, obtain corresponding thermo parameters method knot
Really.
Temperature calculating unit 407 specifically includes:
Convergence determines subelement 4071, for determining the convergence criterion of correspondence according to turbulence model;
Initialize subelement 4072, for transformator simplified model is carried out initialization process;
Temperature calculating subelement 4073, for according to turbulence model, property parameters, conditional parameter and total
Loss, and combine the Solution model that determines and convergence criterion the transformator simplified model after initializing is carried out
Simulation process, obtains corresponding thermo parameters method result;
Wherein, property parameters is material properties parameter, and conditional parameter is boundary condition parameter.
In the present embodiment, the transformator simplified model after being processed simplification by stress and strain model unit 402 is entered
Row stress and strain model, and define corresponding border and interface;Interface signal generating unit 403 is by the second pattern
Merge the grid divided, and generate the interface of correspondence;Acquiring unit 405 gets by the second pattern
Turbulence model, property parameters and the conditional parameter determined;Coupling processing unit 406 carries out electromagnetic field and temperature
The coupling processing of degree field, determines total losses based on temperature;Temperature calculating unit 407 is according to turbulent flow mould
Type, property parameters, conditional parameter and total losses, and combine the Solution model that determines and convergence criterion is carried out
Simulation process, obtains corresponding thermo parameters method result, solves existing technical scheme in the mistake calculated
Calculated in journey is only the overall losses on coil or iron core rather than loss distribution, and causes
The low technical problem of computational accuracy.
Further, by using electromagnetic field and the coupling analysis in temperature field, transformator work can precisely be obtained
Make the Temperature Distribution under state, go the question meanings such as hot-spot great for prevention transformation;The electricity provided
Magnetic field and the coupling process in temperature field, the problem that can effectively solve to be similar to, it is not limited to transformator
Analysis of heat transmission.
Those skilled in the art is it can be understood that arrive, and for convenience and simplicity of description, above-mentioned retouches
The specific works process of the system stated, device and unit, is referred to the correspondence in preceding method embodiment
Process, does not repeats them here.
In several embodiments provided herein, it should be understood that disclosed system, device and
Method, can realize by another way.Such as, device embodiment described above is only shown
Meaning property, such as, the division of described unit, be only a kind of logic function and divide, actual can when realizing
There to be other dividing mode, the most multiple unit or assembly can in conjunction with or be desirably integrated into another
System, or some features can ignore, or do not perform.Another point, shown or discussed each other
Coupling direct-coupling or communication connection can be the INDIRECT COUPLING by some interfaces, device or unit
Or communication connection, can be electrical, machinery or other form.
The described unit illustrated as separating component can be or may not be physically separate, makees
The parts shown for unit can be or may not be physical location, i.e. may be located at a place,
Or can also be distributed on multiple NE.Can select according to the actual needs part therein or
The whole unit of person realizes the purpose of the present embodiment scheme.
It addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit,
Can also be that unit is individually physically present, it is also possible to two or more unit are integrated in a list
In unit.Above-mentioned integrated unit both can realize to use the form of hardware, it would however also be possible to employ software function list
The form of unit realizes.
If described integrated unit realizes and as independent production marketing using the form of SFU software functional unit
Or when using, can be stored in a computer read/write memory medium.Based on such understanding, this
The part that the most in other words prior art contributed of technical scheme of invention or this technical scheme
Completely or partially can embody with the form of software product, this computer software product is stored in one
In storage medium, including some instructions with so that computer equipment (can be personal computer,
Server, or the network equipment etc.) perform completely or partially walking of method described in each embodiment of the present invention
Suddenly.And aforesaid storage medium includes: USB flash disk, portable hard drive, read only memory (ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD
Etc. the various media that can store program code.
The above, above example only in order to technical scheme to be described, is not intended to limit;
Although being described in detail the present invention with reference to previous embodiment, those of ordinary skill in the art should
Work as understanding: the technical scheme described in foregoing embodiments still can be modified by it, or to it
Middle part technical characteristic carries out equivalent;And these amendments or replacement, do not make appropriate technical solution
Essence depart from various embodiments of the present invention technical scheme spirit and scope.
Claims (10)
1. an oil-filled transformer Calculation Method of Temperature Field, it is characterised in that including:
The transformer model set up by first mode is carried out simplification process;
Transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding border and boundary
Face;
Merged the described grid divided by the second pattern, and generate the described interface of correspondence;
Get turbulence model, property parameters and the conditional parameter determined by described second pattern;
Carry out the coupling processing in electromagnetic field and temperature field, determine total losses based on temperature;
According to described turbulence model, described property parameters, described conditional parameter and described total losses, and tie
Close the Solution model determined and convergence criterion carries out simulation process, obtain corresponding thermo parameters method result.
Oil-filled transformer Calculation Method of Temperature Field the most according to claim 1, it is characterised in that
Transformator simplified model after processing simplification carries out stress and strain model, and defines corresponding border and boundary mask
Body includes:
According to the attribute of computational fields or geometric properties or geometry, described transformator simplified model is divided yardstick enter
Row piecemeal processes, and generates a plurality of sub-block;
A plurality of described sub-blocks are carried out stress and strain model process, and defines corresponding border and interface.
Oil-filled transformer Calculation Method of Temperature Field the most according to claim 1, it is characterised in that
Also include before getting turbulence model, property parameters and the conditional parameter determined by described second pattern:
Calculate the Inlet Reynolds Number of described transformator simplified model, determine stream in described transformator simplified model
Move as turbulent flow.
Oil-filled transformer Calculation Method of Temperature Field the most according to claim 1, it is characterised in that
Described property parameters is material properties parameter, and described conditional parameter is boundary condition parameter.
Oil-filled transformer Calculation Method of Temperature Field the most as claimed in any of claims 1 to 4,
It is characterized in that, carry out the coupling processing in electromagnetic field and temperature field, determine that total losses based on temperature are concrete
Including:
Derive the geometric coordinate of thermal source in the fluid model corresponding with described transformator simplified model;
Importing described geometric coordinate in electromagnetic model, interpolation calculation goes out the power density values of respective coordinates;
Described power density values is led back in described fluid model, and is applied to the described geometric coordinate of correspondence
On, determine total losses based on temperature.
Oil-filled transformer Calculation Method of Temperature Field the most according to claim 5, it is characterised in that
According to described turbulence model, described property parameters, described conditional parameter and described total losses, and combine really
Fixed Solution model and convergence criterion carry out simulation process, obtain corresponding thermo parameters method result and specifically wrap
Include:
The described convergence criterion of correspondence is determined according to described turbulence model;
Described transformator simplified model is carried out initialization process;
According to described turbulence model, described property parameters, described conditional parameter and described total losses, and tie
Close the Solution model determined and convergence criterion the described transformator simplified model after initializing is carried out at emulation
Reason, obtains corresponding thermo parameters method result.
7. an oil-filled transformer Temperature calculating device, it is characterised in that including:
Simplify processing unit, for the transformer model set up by first mode is carried out simplification process;
Stress and strain model unit, the transformator simplified model after processing simplification carries out stress and strain model, and
The corresponding border of definition and interface;
Interface signal generating unit, for being merged the described grid divided by the second pattern, and generates correspondence
Described interface;
Acquiring unit, with getting turbulence model, property parameters and the bar determined by described second pattern
Part parameter;
Coupling processing unit, for carrying out the coupling processing in electromagnetic field and temperature field, determines based on temperature
Total losses;
Temperature calculating unit, for according to described turbulence model, described property parameters, described condition ginseng
Number and described total losses, and combine the Solution model that determines and convergence criterion carries out simulation process, it is right to obtain
The thermo parameters method result answered.
Oil-filled transformer Temperature calculating device the most according to claim 7, it is characterised in that
Stress and strain model unit specifically includes:
Partition subelement, for described transformator simplified model according to the attribute of computational fields or geometry
Feature or geometry divide yardstick and carry out piecemeal process, generate a plurality of sub-block;
Stress and strain model subelement, for a plurality of described sub-blocks being carried out stress and strain model process, and define right
Answer border and interface;
Coupling processing unit specifically includes:
Derive subelement, for deriving thermal source in the fluid model corresponding with described transformator simplified model
Geometric coordinate;
Power density computation subelement, for importing described geometric coordinate, interpolation calculation in electromagnetic model
Go out the power density values of respective coordinates;
Apply subelement, for described power density values is led back in described fluid model, and be applied to right
On the described geometric coordinate answered, determine total losses based on temperature.
Oil-filled transformer Temperature calculating device the most according to claim 7, it is characterised in that
Oil-filled transformer Temperature calculating device also includes:
Entrance computing unit, for calculating the Inlet Reynolds Number of described transformator simplified model, determines described
In transformator simplified model, flowing is turbulent flow.
10. fill according to the oil-filled transformer Temperature calculating described in any one in claim 7 to 9
Put, it is characterised in that Temperature calculating unit specifically includes:
Convergence determines subelement, for determining the described convergence criterion of correspondence according to described turbulence model;
Initialize subelement, for described transformator simplified model is carried out initialization process;
Temperature calculating subelement, for according to described turbulence model, described property parameters, described condition
Parameter and described total losses, and combine the Solution model determined and convergence criterion to the described change after initializing
Depressor simplified model carries out simulation process, obtains corresponding thermo parameters method result;
Wherein, described property parameters is material properties parameter, and described conditional parameter is boundary condition parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610463081.8A CN105844069B (en) | 2016-06-21 | 2016-06-21 | A kind of oil-immersed transformer Calculation Method of Temperature Field and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610463081.8A CN105844069B (en) | 2016-06-21 | 2016-06-21 | A kind of oil-immersed transformer Calculation Method of Temperature Field and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105844069A true CN105844069A (en) | 2016-08-10 |
CN105844069B CN105844069B (en) | 2019-02-12 |
Family
ID=56576300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610463081.8A Active CN105844069B (en) | 2016-06-21 | 2016-06-21 | A kind of oil-immersed transformer Calculation Method of Temperature Field and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105844069B (en) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106295018A (en) * | 2016-08-15 | 2017-01-04 | 广东电网有限责任公司电力科学研究院 | A kind of transformator phantom construction method and system |
CN106354971A (en) * | 2016-09-14 | 2017-01-25 | 中国石油大学(华东) | Novel transformer state estimation method based on field-circuit coupling analysis |
CN106570274A (en) * | 2016-11-07 | 2017-04-19 | 广东电网有限责任公司电力科学研究院 | High overload transformer temperature rise calculation method and apparatus |
CN107291969A (en) * | 2016-12-19 | 2017-10-24 | 沈阳工业大学 | The body net coupling analysis method that a kind of Transformer Winding Temperature Rise is calculated |
CN107368638A (en) * | 2017-07-10 | 2017-11-21 | 广东电网有限责任公司电力科学研究院 | A kind of traction electric machine transformer multiple physical field method for numerical simulation and device |
CN108388744A (en) * | 2018-03-08 | 2018-08-10 | 云南电网有限责任公司电力科学研究院 | A kind of emulation mode and system of stable state transformer fuel factor |
CN108491675A (en) * | 2018-05-29 | 2018-09-04 | 中车青岛四方机车车辆股份有限公司 | The computational methods and device of braking system maximum temperature value in a kind of braking process |
CN109085198A (en) * | 2018-07-27 | 2018-12-25 | 昆明理工大学 | A kind of experimental provision and application method measuring transformer oil convection transfer rate |
CN109543246A (en) * | 2018-10-31 | 2019-03-29 | 深圳供电局有限公司 | Method and system for determining distribution of direct-current magnetic bias temperature field of transformer |
CN109858074A (en) * | 2018-12-13 | 2019-06-07 | 云南电网有限责任公司电力科学研究院 | A kind of simulation method of the transformer oil flow surge based on finite volume method |
CN109977549A (en) * | 2019-03-27 | 2019-07-05 | 成都安世亚太科技有限公司 | A kind of novel vertoro thermal model processing method |
CN110427693A (en) * | 2019-07-30 | 2019-11-08 | 中国恩菲工程技术有限公司 | Fluid emulation calculation method in industrial reactor with spray gun |
CN111611753A (en) * | 2020-05-13 | 2020-09-01 | 广东省智能制造研究所 | Design method of blanket with uniform temperature rise |
CN111881611A (en) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | Simulation processing method and device for physical field of motor |
CN112364465A (en) * | 2020-12-07 | 2021-02-12 | 中国市政工程华北设计研究总院有限公司 | Method for calculating height of supporting leg of air-temperature gasifier |
CN113032948A (en) * | 2019-12-24 | 2021-06-25 | 上海电力大学 | Method for obtaining temperature field of oil-immersed transformer |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110068572A (en) * | 2009-12-16 | 2011-06-22 | (주)명진종합건설 | Direct cooling device for oil filled transformer using thermoelectric element |
CN104036125A (en) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | Method for accurately calculating temperature field in oil-immersed transformer |
CN104198068A (en) * | 2013-12-05 | 2014-12-10 | 国家电网公司 | Temperature monitoring device and temperature monitoring method for winding of oil immersed transformer |
CN105045966A (en) * | 2015-06-29 | 2015-11-11 | 上海电机学院 | Hybrid calculation method for inner temperature field of transformer |
-
2016
- 2016-06-21 CN CN201610463081.8A patent/CN105844069B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110068572A (en) * | 2009-12-16 | 2011-06-22 | (주)명진종합건설 | Direct cooling device for oil filled transformer using thermoelectric element |
CN104198068A (en) * | 2013-12-05 | 2014-12-10 | 国家电网公司 | Temperature monitoring device and temperature monitoring method for winding of oil immersed transformer |
CN104036125A (en) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | Method for accurately calculating temperature field in oil-immersed transformer |
CN105045966A (en) * | 2015-06-29 | 2015-11-11 | 上海电机学院 | Hybrid calculation method for inner temperature field of transformer |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106295018B (en) * | 2016-08-15 | 2019-10-22 | 广东电网有限责任公司电力科学研究院 | A kind of transformer simulation model construction method and system |
CN106295018A (en) * | 2016-08-15 | 2017-01-04 | 广东电网有限责任公司电力科学研究院 | A kind of transformator phantom construction method and system |
CN106354971A (en) * | 2016-09-14 | 2017-01-25 | 中国石油大学(华东) | Novel transformer state estimation method based on field-circuit coupling analysis |
CN106570274A (en) * | 2016-11-07 | 2017-04-19 | 广东电网有限责任公司电力科学研究院 | High overload transformer temperature rise calculation method and apparatus |
CN106570274B (en) * | 2016-11-07 | 2020-11-17 | 广东电网有限责任公司电力科学研究院 | Temperature rise calculation method and device for high-overload transformer |
CN107291969A (en) * | 2016-12-19 | 2017-10-24 | 沈阳工业大学 | The body net coupling analysis method that a kind of Transformer Winding Temperature Rise is calculated |
CN107368638A (en) * | 2017-07-10 | 2017-11-21 | 广东电网有限责任公司电力科学研究院 | A kind of traction electric machine transformer multiple physical field method for numerical simulation and device |
CN108388744A (en) * | 2018-03-08 | 2018-08-10 | 云南电网有限责任公司电力科学研究院 | A kind of emulation mode and system of stable state transformer fuel factor |
CN108491675A (en) * | 2018-05-29 | 2018-09-04 | 中车青岛四方机车车辆股份有限公司 | The computational methods and device of braking system maximum temperature value in a kind of braking process |
CN109085198B (en) * | 2018-07-27 | 2022-04-12 | 昆明理工大学 | Experimental measurement device for measuring convection heat transfer coefficient of transformer oil and use method |
CN109085198A (en) * | 2018-07-27 | 2018-12-25 | 昆明理工大学 | A kind of experimental provision and application method measuring transformer oil convection transfer rate |
CN109543246B (en) * | 2018-10-31 | 2023-08-18 | 深圳供电局有限公司 | Method and system for determining distribution of DC magnetic bias temperature field of transformer |
CN109543246A (en) * | 2018-10-31 | 2019-03-29 | 深圳供电局有限公司 | Method and system for determining distribution of direct-current magnetic bias temperature field of transformer |
CN109858074A (en) * | 2018-12-13 | 2019-06-07 | 云南电网有限责任公司电力科学研究院 | A kind of simulation method of the transformer oil flow surge based on finite volume method |
CN109977549B (en) * | 2019-03-27 | 2021-10-22 | 成都安世亚太科技有限公司 | Novel thermal analysis model processing method for transformer rectifier |
CN109977549A (en) * | 2019-03-27 | 2019-07-05 | 成都安世亚太科技有限公司 | A kind of novel vertoro thermal model processing method |
CN110427693A (en) * | 2019-07-30 | 2019-11-08 | 中国恩菲工程技术有限公司 | Fluid emulation calculation method in industrial reactor with spray gun |
CN110427693B (en) * | 2019-07-30 | 2023-03-24 | 中国恩菲工程技术有限公司 | Fluid simulation calculation method in industrial reactor with spray gun |
CN113032948A (en) * | 2019-12-24 | 2021-06-25 | 上海电力大学 | Method for obtaining temperature field of oil-immersed transformer |
CN113032948B (en) * | 2019-12-24 | 2022-04-12 | 上海电力大学 | Method for obtaining temperature field of oil-immersed transformer |
CN111611753A (en) * | 2020-05-13 | 2020-09-01 | 广东省智能制造研究所 | Design method of blanket with uniform temperature rise |
CN111611753B (en) * | 2020-05-13 | 2024-01-30 | 广东省智能制造研究所 | Design method of blanket with uniform temperature rise |
CN111881611A (en) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | Simulation processing method and device for physical field of motor |
CN111881611B (en) * | 2020-07-31 | 2023-07-14 | 珠海格力电器股份有限公司 | Simulation processing method and device for motor physical field |
CN112364465A (en) * | 2020-12-07 | 2021-02-12 | 中国市政工程华北设计研究总院有限公司 | Method for calculating height of supporting leg of air-temperature gasifier |
CN112364465B (en) * | 2020-12-07 | 2022-04-15 | 中国市政工程华北设计研究总院有限公司 | Method for calculating height of supporting leg of air-temperature gasifier |
Also Published As
Publication number | Publication date |
---|---|
CN105844069B (en) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105844069A (en) | Temperature field calculation method and device of oil-immersed transformer | |
Chen et al. | CFD investigation on thermal-hydraulic behaviors of a wire-wrapped fuel subassembly for sodium-cooled fast reactor | |
Guo et al. | Theory of the lattice Boltzmann equation: lattice Boltzmann model for axisymmetric flows | |
CN107368638A (en) | A kind of traction electric machine transformer multiple physical field method for numerical simulation and device | |
CN104978456B (en) | A kind of universal valve design optimization method | |
CN110147559B (en) | Converter multidisciplinary optimization design method based on multi-physics coupling | |
CN107885955A (en) | The emulation mode and device of electric motor temperature field | |
Chen et al. | Thermal–hydraulic performance of a 5× 5 rod bundle with spacer grid in a nuclear reactor | |
CN111027112B (en) | Porous medium simulation method for fast reactor rod bundle assembly coupling heat transfer model | |
Gajbhiye et al. | Validation and verification of a robust 3-D MHD code | |
Končar et al. | Simulation of turbulent flow in MATIS-H rod bundle with split-type mixing vanes | |
CN110175377A (en) | A kind of Numerical Prediction Method about hydrofoil cloud cavitation | |
Sun et al. | Heat transfer performance prediction of Taylor–Couette flow with longitudinal slits using artificial neural networks | |
CN106202834A (en) | The computational methods of a kind of three-phase transformer oil stream motion and device | |
Li et al. | Cooling water jacket design of motorized spindle system using multi-objective topology optimization | |
Lim et al. | 2-D lumped-parameter thermal modelling of axial flux permanent magnet generators | |
Karaman et al. | Numerical investigation of flow characteristics through simple support grids in a 1 x 3 rod bundle | |
Sergio et al. | Optimization methodology assessment for the inlet velocity profile of a hydraulic turbine draft tube: part II—performance evaluation of draft tube model | |
Ye et al. | Numerical analysis of the flow behavior in a helically coiled once through steam generator | |
Chen et al. | Challenge analysis and schemes design for the CFD simulation of PWR | |
CN115169265B (en) | Method, system, equipment and medium for analyzing mixing coefficient based on numerical analysis | |
Huang et al. | Numerical simulations of fluid–structure interaction based on Cartesian grids with two boundary velocities | |
Wang et al. | Development and application of CFD and subchannel coupling analysis code for lead‐cooled fast reactor | |
CN115422808A (en) | Transformer temperature field model order reduction method based on Krylov subspace | |
Tang et al. | Multi-physics topology optimization for thermal-flow problems applied to additively manufactured heat exchangers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |