CN107992664A - A kind of change of current power transformation-thermal coupling field computation method based on improvement gridless routing - Google Patents
A kind of change of current power transformation-thermal coupling field computation method based on improvement gridless routing Download PDFInfo
- Publication number
- CN107992664A CN107992664A CN201711207015.5A CN201711207015A CN107992664A CN 107992664 A CN107992664 A CN 107992664A CN 201711207015 A CN201711207015 A CN 201711207015A CN 107992664 A CN107992664 A CN 107992664A
- Authority
- CN
- China
- Prior art keywords
- electric field
- field
- temperature
- calculation
- calculating
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The present invention is claimed a kind of based on the change of current change electro thermal coupling field computation method for improving gridless routing, including definite extra-high voltage converter change empire paper material mathematics nonlinear model, and the non-linear electro thermal coupling field computation model of two-dimensional axial symmetric；Change of current change electro thermal coupling field computation is carried out using non-mesh method is improved, model stationing mode is adaptively adjusted according to size, exchange rheology electric field and temperature field are respectively calculated, the local subdomain size of borderline region adjustment is calculated closing on, itself and global boundary-intersected are avoided, is calculated on global border using point collocation.In electric field and temperature field coupling data transmit, when needing to calculate temperature field thermal force, part subdomain Gauss integration point local subdomain in electric Field Calculation model in temperature field is searched according to geometric coordinate, and according to the part subdomain internal segment count out interpolation obtain thermic load value, electric field material attribute is updated, realizes that the change of current becomes the iterative calculation of electro thermal coupling field and the regularity of distribution obtains.The present invention, which improves, calculates time and precision.
Description
Technical Field
The invention belongs to the technical field of evaluation of the insulation level of extra-high voltage equipment, and particularly relates to a calculation method of an extra-high voltage converter transformer-thermal coupling field.
Background
Because of the imbalance of domestic resource distribution, extra-high voltage is being vigorously constructed, which is an important strategy for thirteen-five planning, and especially with the proposal of an energy internet construction strategy, the safety and reliability of an extra-high voltage system are greatly concerned.
The reliability of the extra-high voltage converter transformer serving as key equipment of an extra-high voltage system directly determines the safe operation of the system, but the current converter transformer key technology still depends on foreign technologies mainly due to insulation problems. Because the converter transformer insulation component has an alternating current-direct current composite action working condition and the conductivity of oil and oil paper is greatly influenced by temperature, the influence of an electric field and a temperature field of the converter transformer insulation component must be analyzed when the converter transformer insulation problem is analyzed, namely an electric-thermal coupling field calculation model of the converter transformer is constructed. The Harbin university of mechanics and the Sigan university of transportation analyze the electric field distribution rule under the AC/DC working condition, the North China university of electric power considers the nonlinear problem of the electric field, and the influence of actual temperature distribution on the electric field is not considered during the existing design, but the temperature and the corresponding oil and oil paper electric parameters under the most severe insulation condition are selected, so that the problem of overlarge insulation margin design exists.
Therefore, the inventor finds that the converter transformer-thermal coupling problem is rarely researched at home and abroad, the grid analysis method adopted due to the large size difference of the converter transformer insulating parts has the problems of large grid number, low calculation efficiency and the like, and the coupling data transmission of the electric field and the temperature field also has the problem of mapping precision, so that the improvement on the aspects is needed.
With the rapid development of the ultrahigh voltage system in China, further intensive research is needed on the insulation design and heat dissipation design problems of the ultrahigh voltage converter transformer. At this time, research on a calculation method of the coupling characteristic of the extra-high voltage converter transformer electric field-thermal field is necessary.
Disclosure of Invention
The present invention is directed to solving the above problems of the prior art. The method for calculating the converter transformer-thermal coupling field based on the improved meshless method is effectively selected for improving the insulation margin of the converter transformer. The technical scheme of the invention is as follows:
a method for calculating a converter transformer-thermal coupling field based on an improved gridless method comprises the following steps:
101. determining the structural size and the properties of insulating oil and insulating paperboard materials of the extra-high voltage converter transformer calculation model, determining the coupling characteristics of an electric field and a temperature field of the extra-high voltage converter transformer according to the composite action operation condition of alternating current and direct current borne by the converter transformer, and establishing an extra-high voltage electricity-thermal coupling field calculation model;
102. setting initial values of an electric field and temperature of the whole model when the converter transformer is initially calculated, respectively solving the electric field and the temperature field by adopting an improved gridless method, adjusting the size of a node local sub-domain close to a boundary in the calculation by adopting the improved gridless method, avoiding intersection with a global boundary, and calculating at the boundary node by adopting a point matching method;
103. searching a local sub-domain of a Gaussian integration point of a local sub-domain of the temperature field in the electric field calculation model according to the geometric coordinates, interpolating according to the number of nodes in the local sub-domain in the electric field calculation model to obtain a heat load value, updating the material attribute of the electric field by adopting a method of searching the Gaussian integration point of the local sub-domain, and realizing electric-thermal coupling iterative calculation;
104. and (4) repeating the steps 101-103 to obtain the distribution rule of the electric-thermal field of the key insulating component of the converter transformer.
Further, in step 102, the electric field and the temperature field are respectively solved by using an improved gridless method, and the solved field and temperature field control equations are as follows:
electric field control equation: considering that the alternating current and direct current composite operation condition exists in the converter transformer, an electric field control equation is as follows:
wherein gamma is the conductivity and epsilon is the relative dielectric constant; x, y and z respectively represent the coordinates of the calculation area;
temperature field control equation: the control equation of the internal heat conduction of the converter transformer at the steady state is as follows:
wherein λ is thermal conductivity, W/(m.K); t is the medium temperature, DEG C; q v Is the heat generation rate per unit volume, W/m 3 。
Wherein the temperature field calculates the required Q v Can be calculated by the following formula.
Wherein J is a total current density including a source current density and an eddy current density, A/m 2 And σ represents the conductivity, S/m
Further, in step 102, the electric field and the temperature field are respectively solved by using an improved gridless method, and the solved electric field boundary conditions and temperature field boundary conditions are as follows:
electric field boundary conditions: when considering the electro-thermal coupling effect of the converter transformer insulating part, the direct current condition of the valve side winding needs to be considered, high voltage is loaded on the valve side winding and the electrostatic ring part respectively during calculation, and the earth potential is loaded on the grid side winding and the electrostatic ring connected with the grid side winding, wherein the specific mathematical expressions are respectively as follows:
in addition, the mathematical relation of the conductivity and the temperature of the insulating oil and the insulating paper needs to be considered when a conversion current transformation-thermal coupling model is constructed, the expression of which is shown as (5),
γ＝κ 0 exp(-F/k.T)(5)
wherein k is 0 F is a correlation coefficient and can be obtained by fitting according to the conductivity change rule of the insulating oil and the insulating paper at different temperatures;
temperature field boundary: the surface of the oiled paper radiates heat to the outside through natural convection, and the boundary conditions can be expressed as:
wherein h is the surface convection heat transfer coefficient, W/(m) 2 ·K)；T f The surface temperature of the heating element is DEG C; t is amb Ambient temperature, deg.C;
and the other boundary of the converter transformer insulating member is set to be thermally insulated as shown in formula (7).
Further, the step 103 specifically includes the steps of: setting a temperature initial value, calculating material conductivity according to the temperature initial value, searching a local subdomain of a node in an electric field calculation model according to the coordinate position of a Gaussian integral point in the local subdomain of the node when calculating the thermal load of a temperature field, searching the number of surrounding nodes according to the size of the local subdomain of the node, interpolating by adopting a point matching method of an RPIMp-shaped function in a non-grid method to obtain a thermal load value, searching a mapping point of the node material attribute in the electric field calculation model in the temperature field by adopting the point matching method based on the RPIMp-shaped function, interpolating to obtain an electric field conductivity updating parameter according to the number of nodes in the mapped local subdomain, calculating electric field distribution, judging whether the difference value of two adjacent calculation results meets the calculation error requirement, and if not, performing the next iteration until the difference value of the two adjacent iteration calculation results meets the control accuracy requirement.
Further, when the improved meshless method is used for calculating the electric field and the temperature field, the method comprises the following steps:
calculating a current transformation electric field to obtain electric field distribution of the current transformation electric field, and calculating a model according to a current transformation temperature field; searching local sub-domains to which Gaussian integration points in each local sub-domain belong in an electric field model according to the geometric coordinates, calculating the electric field distribution of the local sub-domains by adopting an RPIMp shape function, and calculating the heat q in the unit volume of the local sub-domains according to a loss calculation formula;
and calculating the temperature distribution of the heat q by adopting an improved non-grid method according to the obtained heat q and combining a temperature field control equation and boundary conditions, calculating the temperature t of each node in an electric field model, and updating the material property according to the relation between the conductivity and the temperature.
The invention has the following advantages and beneficial effects:
the invention adopts the non-grid calculation of the current conversion transformation electric field and the temperature field, thereby avoiding grid subdivision and reducing the calculation amount. Meanwhile, aiming at the problem of irregular boundary, an improved non-grid method for self-adaptive adjustment of a local sub-domain is provided, and aiming at the problem of data transmission of coupling of an electric field and a temperature field, a coupling data transmission method utilizing Gaussian integral point mapping and node interpolation in the local sub-domain is provided on the basis of the improved non-grid method, so that rapid calculation of a converter transformer-thermal coupling field can be realized, and guidance is provided for insulation design.
Drawings
FIG. 1 is a flow chart of the calculation of the converter transformer-thermal coupling field by the improved gridless method of the preferred embodiment of the present invention;
fig. 2 is a modified meshless method local sub-domain adaptive adjustment method for calculating a converter transformer electric field and a temperature field and a conventional meshless calculation region determination method provided by an embodiment of the present invention;
FIG. 3 is a flow chart of a current converting transformer/temperature field calculation based on an improved MLPG meshless method according to an embodiment of the present invention;
fig. 4 is a schematic diagram of electric field and temperature field coupling data transmission by using gaussian integration points and node interpolation in a mapping local sub-domain in the improved meshless method according to the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described in detail and clearly with reference to the accompanying drawings. The described embodiments are only some of the embodiments of the present invention.
The technical scheme for solving the technical problems is as follows:
referring to fig. 1, determining structural parameters and material properties of an extra-high voltage converter transformer calculation model, determining electric field and temperature field coupling characteristics of the extra-high voltage converter transformer calculation model according to an operation condition, and further establishing a two-dimensional axisymmetric electric-thermal coupling field calculation model of the extra-high voltage converter transformer calculation model;
referring to fig. 2, according to the shape of the calculation model, the calculation model is subjected to point distribution discretization, meanwhile, in order to enhance the calculation efficiency and accuracy, a dense point distribution mode is adopted in a region with a small size, a plurality of nodes are also designed on the boundary, different point distribution modes are adopted for the electric field and the temperature field, so that the requirements of the calculation accuracy of the electric field and the temperature field on the discretization mode are met, and corresponding material parameters and boundary conditions are set.
Referring to fig. 3, for the condition that the converter transformer bears the alternating current-direct current combined action, the non-linear problem of the electrical-thermal coupling is considered, an improved non-grid method is adopted to respectively perform solving operation on the electric field and temperature field calculation models in the step 1, a local sub-domain close to a global boundary node is adjusted in the calculation to avoid intersection with the global boundary, and a point matching method is adopted to directly perform calculation on the boundary.
Referring to fig. 4, on the basis of the above improved non-grid method, it is necessary to transmit electric field and temperature field coupling data, first setting a temperature initial value, calculating material conductivity according to the value, searching a local sub-domain of a node where the node is located in an electric field calculation model according to a coordinate position of a gaussian integral point in a node local sub-domain when calculating temperature field thermal load, searching a number of surrounding nodes according to the size of the node local sub-domain, obtaining a thermal load value thereof by using an RPIMp shape function interpolation in the non-grid method, searching a mapping point of a node material property in the electric field calculation model in the same method in the temperature field, obtaining a conductivity update parameter thereof by using the interpolation of the number of nodes in the mapped local sub-domain, calculating electric field distribution, and judging whether a difference value between two adjacent calculation results meets a calculation error requirement, if not, entering next iteration until the difference value between the two adjacent iteration calculation results meets a control accuracy requirement. And 5, calculating the change rules of the electric field and the temperature field under different operation conditions of the converter transformer according to the calculation method, and comparing and analyzing the efficiency and the precision of the proposed calculation method by utilizing the existing simulation tool so as to verify the effectiveness of the proposed method.
The electric field and temperature field control equation and the applied boundary conditions in the second step are as follows:
1) Electric field control equation
Considering that the alternating current and direct current composite operation condition exists in the converter transformer, an electric field control equation is as follows:
wherein gamma is the conductivity and epsilon is the relative dielectric constant
2) Temperature field control equation
The control equation of the internal heat conduction of the converter transformer at the steady state is as follows:
wherein λ is thermal conductivity, W/(m.K); t is the medium temperature, DEG C; q v Is the heat generation rate per unit volume, W/m 3 。
Wherein the temperature field calculates the required Q v Can be calculated by the following formula.
Wherein J is a total current density including a source current density and an eddy current density, A/m 2 。
3) Boundary condition of electric field
When considering the electro-thermal coupling effect of the converter transformer insulating part, the condition that direct current exists in a valve side winding needs to be considered, high voltage is loaded on the valve side winding and an electrostatic ring part respectively according to an experimental scheme during calculation, and the ground potential is added on a grid side winding and an electrostatic ring connected with the grid side winding, wherein the specific mathematical expressions are respectively as follows:
in addition, the mathematical relation of the conductivity and the temperature of the insulating oil and the insulating paper needs to be considered when a conversion current transformation-thermal coupling model is constructed, and the expression is shown as (19).
γ＝κ 0 exp(-F/k.T) (19)
Wherein k is 0 And F is a correlation coefficient and can be obtained by fitting according to the conductivity change rule of the insulating oil and the insulating paper at different temperatures.
4) Temperature field boundary
The surface of the oiled paper radiates heat to the outside through natural convection, and the boundary conditions can be expressed as:
wherein h is the surface convection heat transfer coefficient, W/(m) 2 ·K)；T f The surface temperature of the heating element is DEG C; t is amb Is at ambient temperature, DEG C.
And the other boundary of the converter transformer insulating member is set to be thermally insulated as shown in equation (21).
The specific implementation conditions of the improved non-grid method in the third step and the fourth step are as follows:
1) A shape function construction method and an equation are adopted, a RPIMp method is adopted to construct a shape function without a grid algorithm, and a Galerkin discrete equation calculated by an electric field and a temperature field is obtained.
2) In the conventional MLPG method, calculation is performed according to a node local sub-domain, the shape of the local sub-domain is generally circular or rectangular, as shown in FIG. 1, a node which is not intersected with a boundary in a calculation domain is directly subjected to Gaussian calculation, but whether the node local sub-domain is intersected with the boundary or not needs to be judged when an irregular boundary condition exists, if the node local sub-domain is not intersected with the boundary, calculation is directly performed according to an internal node method, if the node local sub-domain is intersected with the boundary, the intersection region is determined, the shape of the local sub-domain and the intersection region are re-determined, the position relation between a Gaussian integral point and the boundary needs to be considered when Gaussian integration is adopted, the node is ensured to be in the calculation domain range, and formulas (4), (6) to (7) are calculated in an integral mode for the boundary node.
3) And for the nodes close to the global boundary, calculating the minimum distance d between the node and the boundary close to the node, multiplying the minimum distance d by a scaling coefficient k, wherein k <1, and taking the minimum distance as the radius of the local sub-domain, so that no intersection region exists between the local sub-domain of the node close to the global boundary and the global boundary. Further, nodes are arranged on the boundary according to the degree of shape irregularity, and since the RPIMp shape function has a functional property, the function is adopted as a shape function of the boundary node calculation, and a lattice-fitting method is adopted, so that the boundary node electric field and temperature field calculations satisfy formulas (4), (6) to (7).
4) When the improved meshless method is used for calculating the electric field and the temperature field, the method comprises the following steps:
1. the temperature is started and initialized.
2. And calculating the current conversion transformation electric field by adopting an improved non-grid method to obtain the electric field distribution of the current conversion transformation electric field. According to the converter transformer temperature field calculation model, local sub-domains of Gaussian integration points in each local sub-domain in the electric field model are searched according to geometric coordinates, the electric field distribution of the local sub-domains is calculated by adopting an RPIMp shape function, and the heat q in the unit volume of the local sub-domains is calculated according to a loss calculation formula.
3. And (3) calculating the temperature distribution of the material by adopting an improved non-grid method according to the obtained heat q and by combining a temperature field control equation and boundary conditions, calculating the temperature t of each node in the electric field model by adopting the same method in the step (2), and updating the material property according to the relation between the conductivity and the temperature.
4. And updating and calculating the distribution of the converter transformer field and the distribution of the temperature field again according to the updated material properties, repeating the steps 2) and 3) until a stable state is reached, and calculating the distribution rule of the converter transformer field and the temperature field.
It should be noted that, in the foregoing system embodiment, each included system unit is only divided according to functional logic, but is not limited to the above division as long as the corresponding function can be implemented; in addition, specific names of the functional units are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present invention.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by relevant hardware instructed by a program, and the program may be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc.
The above examples are to be construed as merely illustrative and not limitative of the remainder of the disclosure. After reading the description of the invention, the skilled person can make various changes or modifications to the invention, and these equivalent changes and modifications also fall into the scope of the invention defined by the claims.
Claims (5)
1. A method for calculating a converter transformer-thermal coupling field based on an improved gridless method is characterized by comprising the following steps of:
101. determining the structural size and the material properties of insulating oil and insulating paper boards of the ultra-high voltage converter transformer calculation model, determining the coupling characteristics of an electric field and a temperature field of the ultra-high voltage converter transformer according to the operating condition of the converter transformer under the combined action of alternating current and direct current, and establishing an ultra-high voltage electricity-thermal coupling field calculation model;
102. setting initial values of an electric field and temperature of the whole model when the converter transformer is initially calculated, respectively solving the electric field and the temperature field by adopting an improved gridless method, adjusting the size of a node local sub-domain close to a boundary in the calculation by adopting the improved gridless method, avoiding intersection with a global boundary, and calculating at the boundary node by adopting a point matching method;
103. searching a local sub-domain of a Gaussian integral point of a local sub-domain of the temperature field in the electric field calculation model according to the geometric coordinates, interpolating according to the number of nodes in the local sub-domain in the electric field calculation model to obtain a heat load value, updating the material attribute of the electric field in a mode of searching the Gaussian integral point of the local sub-domain, and realizing electric-thermal coupling iterative calculation;
104. and (4) repeating the steps 101-103 to obtain the distribution rule of the electric-thermal field of the key insulating component of the converter transformer.
2. The method for calculating the conversion current transformer-thermal coupling field based on the improved gridless method according to claim 1, wherein the step 102 is implemented by solving the electric field and the temperature field respectively by using the improved gridless method, and the solved field and temperature field control equations are as follows:
electric field control equation: considering that the alternating current and direct current composite operation condition exists in the converter transformer, an electric field control equation is as follows:
wherein gamma is the conductivity, and epsilon is the relative dielectric constant; x, y and z respectively represent the coordinates of the calculation area;
temperature field control equation: the control equation of the internal heat conduction of the converter transformer at the steady state is as follows:
wherein λ is thermal conductivity, W/(m.K); t is the medium temperature, DEG C; q v Is the heat generation rate per unit volume, W/m 3 ；
Wherein the temperature field calculates the required Q v Can be calculated by the following formula
Wherein σ represents the conductivity, S/m, J is the total current density including the source current density and the eddy current density, A/m 2 。
3. The method for calculating the converter transformer-thermal coupling field based on the improved meshless method according to the claim 1 or 2, wherein the step 102 adopts the improved meshless method to respectively solve the electric field and the temperature field, and the boundary conditions of the electric field and the temperature field are as follows:
electric field boundary conditions: when considering the electro-thermal coupling effect of the converter transformer insulating part, the direct current condition of the valve side winding needs to be considered, high voltage is respectively loaded on the valve side winding and the electrostatic ring part during calculation, and the ground potential is respectively added on the grid side winding and the electrostatic ring connected with the grid side winding, wherein the specific mathematical expressions are respectively as follows:
in addition, the mathematical relation of the conductivity and the temperature of the insulating oil and the insulating paper needs to be considered when a conversion current transformation-thermal coupling model is constructed, the expression of which is shown as (5),
γ＝κ 0 exp(-F/k.T) (5)
wherein k is 0 F is a correlation coefficient and can be obtained by fitting according to the conductivity change rule of the insulating oil and the insulating paper at different temperatures;
temperature field boundary: the surface of the oiled paper radiates heat to the outside through natural convection, and the boundary conditions can be expressed as:
wherein h is the surface convection heat transfer coefficient, W/(m) 2 ·K)；T f The surface temperature of the heating element is DEG C; t is amb Ambient temperature, deg.C;
and the other boundary of the converter transformer insulating member is set to be thermally insulated as shown in formula (7).
4. The improved meshless method based current conversion power transformation-thermal coupling field calculation method according to claim 3, wherein the step 103 specifically comprises the steps of: setting a temperature initial value, calculating material conductivity according to the temperature initial value, searching a local subdomain of a node in an electric field calculation model according to a coordinate position of a Gaussian integration point in a node local subdomain when calculating the temperature field heat load, searching the number of surrounding nodes according to the size of the node local subdomain, interpolating by adopting a point distribution method of an RPIMp shape function in a grid-free method to obtain a heat load value, searching a mapping point of a node material attribute in the electric field calculation model in the temperature field by adopting a point distribution method based on the RPIMp shape function, solving an electric field conductivity updating parameter according to the number interpolation of the nodes in the mapped local subdomain, calculating electric field distribution, judging whether a difference value of two adjacent calculation results meets a calculation error requirement, and if not, entering next iteration until the difference value of the two adjacent iteration calculation results meets a control accuracy requirement.
5. The improved meshless method based converter transformer-thermal coupling field calculation method is characterized in that when the improved meshless method is used for calculating the electric field and the temperature field, the method comprises the following steps:
calculating a current transformation electric field to obtain electric field distribution of the current transformation electric field, and calculating a model according to a current transformation temperature field; searching local sub-domains to which Gaussian integration points in each local sub-domain belong in an electric field model according to the geometric coordinates, calculating the electric field distribution of the local sub-domains by adopting an RPIMp-shaped function, and calculating the heat q in the unit volume of the local sub-domains according to a loss calculation formula;
and calculating the temperature distribution of the heat q by adopting an improved non-grid method according to the obtained heat q and combining a temperature field control equation and boundary conditions, calculating the temperature t of each node in an electric field model, and updating the material property according to the relation between the conductivity and the temperature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711207015.5A CN107992664B (en) | 2017-11-27 | 2017-11-27 | Improved gridless method-based current conversion transformer-thermal coupling field calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711207015.5A CN107992664B (en) | 2017-11-27 | 2017-11-27 | Improved gridless method-based current conversion transformer-thermal coupling field calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107992664A true CN107992664A (en) | 2018-05-04 |
| CN107992664B CN107992664B (en) | 2021-04-13 |
Family
ID=62033356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711207015.5A Active CN107992664B (en) | 2017-11-27 | 2017-11-27 | Improved gridless method-based current conversion transformer-thermal coupling field calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107992664B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763688A (en) * | 2018-05-17 | 2018-11-06 | 西北工业大学 | A kind of numerical computation method of electrostatic suspension liquid alloy heat transfer |
| CN109117458A (en) * | 2018-06-28 | 2019-01-01 | 浙江省电力有限公司电力科学研究院 | A kind of DC equipment contamination accumulation characteristics calculation method based on improvement gridless routing |
| CN109783855A (en) * | 2018-12-11 | 2019-05-21 | 重庆大学 | A kind of change of current based on upper non-mesh method becomes the calculation method of space charge |
| CN109858141A (en) * | 2019-01-28 | 2019-06-07 | 天津大学 | 220kVGIL insulator method for equalizing voltage based on nonlinear conductance epoxy resin |
| CN110765685A (en) * | 2019-10-18 | 2020-02-07 | 南方电网科学研究院有限责任公司 | Simulation method and device for multi-physical-field coupling of reactor and storage medium |
| CN111368385A (en) * | 2018-12-25 | 2020-07-03 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for operating characteristics of sleeve in electric-thermal coupling field |
| CN111814380A (en) * | 2020-07-16 | 2020-10-23 | 重庆大学 | Method for rapidly calculating temperature distribution of oil-immersed transformer winding by adopting two-dimensional axisymmetric model |
| CN111985075A (en) * | 2020-07-03 | 2020-11-24 | 国网山东省电力公司电力科学研究院 | Temperature distribution calculation method and system suitable for zinc oxide arrester |
| CN112084724A (en) * | 2020-09-08 | 2020-12-15 | 北京航空航天大学 | CFD local rapid calculation method |
| CN113723017A (en) * | 2021-07-14 | 2021-11-30 | 中国航发沈阳发动机研究所 | Method and device for determining temperature field of wall surface of rotating disk of aircraft engine |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106294914A (en) * | 2015-06-04 | 2017-01-04 | 长沙理工大学 | A kind of converter transformer valve-side three-dimensional electric field emulation mode considering Anisotropic Nonlinear |
| CN106442124A (en) * | 2016-11-30 | 2017-02-22 | 华南理工大学 | Shape-memory alloy material fatigue performance test system in electricity-thermal coupling field |
-
2017
- 2017-11-27 CN CN201711207015.5A patent/CN107992664B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106294914A (en) * | 2015-06-04 | 2017-01-04 | 长沙理工大学 | A kind of converter transformer valve-side three-dimensional electric field emulation mode considering Anisotropic Nonlinear |
| CN106442124A (en) * | 2016-11-30 | 2017-02-22 | 华南理工大学 | Shape-memory alloy material fatigue performance test system in electricity-thermal coupling field |
Non-Patent Citations (1)
| Title |
|---|
| 梁永春: "高压电力电缆温度场和载流量评估研究动态", 《高压电技术》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763688B (en) * | 2018-05-17 | 2022-08-05 | 西北工业大学 | Numerical calculation method for heat transfer of electrostatic suspension liquid alloy |
| CN108763688A (en) * | 2018-05-17 | 2018-11-06 | 西北工业大学 | A kind of numerical computation method of electrostatic suspension liquid alloy heat transfer |
| CN109117458A (en) * | 2018-06-28 | 2019-01-01 | 浙江省电力有限公司电力科学研究院 | A kind of DC equipment contamination accumulation characteristics calculation method based on improvement gridless routing |
| CN109783855A (en) * | 2018-12-11 | 2019-05-21 | 重庆大学 | A kind of change of current based on upper non-mesh method becomes the calculation method of space charge |
| CN111368385A (en) * | 2018-12-25 | 2020-07-03 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for operating characteristics of sleeve in electric-thermal coupling field |
| CN111368385B (en) * | 2018-12-25 | 2023-06-27 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Simulation calculation method for sleeve operation characteristics in electrothermal coupling field |
| CN109858141A (en) * | 2019-01-28 | 2019-06-07 | 天津大学 | 220kVGIL insulator method for equalizing voltage based on nonlinear conductance epoxy resin |
| CN110765685A (en) * | 2019-10-18 | 2020-02-07 | 南方电网科学研究院有限责任公司 | Simulation method and device for multi-physical-field coupling of reactor and storage medium |
| CN111985075A (en) * | 2020-07-03 | 2020-11-24 | 国网山东省电力公司电力科学研究院 | Temperature distribution calculation method and system suitable for zinc oxide arrester |
| CN111985075B (en) * | 2020-07-03 | 2024-03-26 | 国网山东省电力公司电力科学研究院 | Temperature distribution calculation method and system suitable for zinc oxide lightning arrester |
| CN111814380A (en) * | 2020-07-16 | 2020-10-23 | 重庆大学 | Method for rapidly calculating temperature distribution of oil-immersed transformer winding by adopting two-dimensional axisymmetric model |
| CN111814380B (en) * | 2020-07-16 | 2023-10-24 | 重庆大学 | Method for rapidly calculating temperature distribution of winding of oil-immersed transformer by adopting two-dimensional axisymmetric model |
| CN112084724A (en) * | 2020-09-08 | 2020-12-15 | 北京航空航天大学 | CFD local rapid calculation method |
| CN113723017A (en) * | 2021-07-14 | 2021-11-30 | 中国航发沈阳发动机研究所 | Method and device for determining temperature field of wall surface of rotating disk of aircraft engine |
| CN113723017B (en) * | 2021-07-14 | 2024-01-02 | 中国航发沈阳发动机研究所 | Method and device for determining temperature field of wall surface of rotary disk of aero-engine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107992664B (en) | 2021-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107992664B (en) | Improved gridless method-based current conversion transformer-thermal coupling field calculation method | |
| CN108052697B (en) | Gallery pipe GI L three-dimensional temperature field and expansion deformation calculation method based on workbench | |
| CN109086514B (en) | Reactor temperature field calculation method based on multi-physical field coupling | |
| CN105045966A (en) | Hybrid calculation method for inner temperature field of transformer | |
| Minaei et al. | Thermal resistance capacity model for short-term borehole heat exchanger simulation with non-stiff ordinary differential equations | |
| CN106021676A (en) | Multi-circle cable steady-state temperature rise acquiring method based on transfer matrix | |
| CN104050307B (en) | Microwave passive circuit 9 Cr 2 steel using electromagnetic heating Method In The Whole-process Analysis based on time domain spectral element method | |
| Kim et al. | A hybrid reduced model for borehole heat exchangers over different time-scales and regions | |
| Sowmya et al. | Analysis of heat transfer through different profiled longitudinal porous fin by differential transformation method | |
| Javed et al. | Vertical ground heat exchangers: A review of heat flow models | |
| Yun et al. | A dual-reciprocity boundary element approach for axisymmetric nonlinear time-dependent heat conduction in a nonhomogeneous solid | |
| Ranjan et al. | Closed form solution for a conductive–convective–radiative annular fin with multiple nonlinearities and its inverse analysis | |
| Kim et al. | Three-dimensional numerical modeling of vertical ground heat exchangers: Domain decomposition and state model reduction | |
| Molavi et al. | A novel methodology for combined parameter and function estimation problems | |
| Valadbeygi et al. | Prediction of Heating Energy Consumption in Houses via Deep Learning Neural Network | |
| Gaba et al. | Thermal performance of functionally graded parabolic annular fins having constant weight | |
| Bornschlegell et al. | Thermal optimization of a single inlet T-junction | |
| Azizian et al. | A dynamic thermal based reliability model of cast-resin dry-type transformers | |
| Moguen et al. | Pressure–velocity coupling for unsteady low Mach number flow simulations: An improvement of the AUSM+-up scheme | |
| CN115730491A (en) | Method for calculating temperature field distribution in transformer and related equipment | |
| CN115859720A (en) | Method for calculating temperature field distribution in transformer and related equipment | |
| Ye et al. | Numerical research on hydraulic and thermal performance of the motor water jackets based on the orthogonal experiment | |
| Zhou et al. | Numerical study on characteristics of flow and thermal fields of tandem cylinders | |
| Cisilino | Topology optimization of 2D potential problems using boundary elements | |
| Zia Zahedi et al. | Nonlinear Adaptive Magneto-Thermal Analysis at Bushing Regions of a Transformers Cover Using Finite Difference Method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |