CN105844069A - Temperature field calculation method and device of oil-immersed transformer - Google Patents

Temperature field calculation method and device of oil-immersed transformer Download PDF

Info

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
Application number
CN201610463081.8A
Other languages
Chinese (zh)
Other versions
CN105844069B (en
Inventor
李德波
许凯
钟俊
冯永新
杨贤
周丹
孟源源
柯春俊
刘石
林春耀
饶章权
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201610463081.8A priority Critical patent/CN105844069B/en
Publication of CN105844069A publication Critical patent/CN105844069A/en
Application granted granted Critical
Publication of CN105844069B publication Critical patent/CN105844069B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal 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

A kind of oil-filled transformer Calculation Method of Temperature Field and device
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.
CN201610463081.8A 2016-06-21 2016-06-21 A kind of oil-immersed transformer Calculation Method of Temperature Field and device Active CN105844069B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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