CN113378395B - Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor - Google Patents
Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor Download PDFInfo
- Publication number
- CN113378395B CN113378395B CN202110682418.5A CN202110682418A CN113378395B CN 113378395 B CN113378395 B CN 113378395B CN 202110682418 A CN202110682418 A CN 202110682418A CN 113378395 B CN113378395 B CN 113378395B
- Authority
- CN
- China
- Prior art keywords
- motor
- node
- ventilation
- temperature rise
- cooling gas
- 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.)
- Active
Links
- 238000009423 ventilation Methods 0.000 title claims abstract description 70
- 238000010168 coupling process Methods 0.000 title claims abstract description 26
- 238000001816 cooling Methods 0.000 title claims abstract description 23
- 230000002457 bidirectional effect Effects 0.000 title claims abstract description 11
- 239000012530 fluid Substances 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000008878 coupling Effects 0.000 claims abstract description 19
- 238000005859 coupling reaction Methods 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 16
- 239000000112 cooling gas Substances 0.000 claims description 55
- 238000007731 hot pressing Methods 0.000 claims description 12
- 239000007787 solid Substances 0.000 claims description 8
- 238000010438 heat treatment Methods 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 239000012809 cooling fluid Substances 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000004134 energy conservation Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000004804 winding Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 229910052802 copper Inorganic materials 0.000 description 4
- 239000010949 copper Substances 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Motor Or Generator Cooling System (AREA)
Abstract
The invention discloses a flow heat network modeling and bidirectional coupling method of a shaft radial mixed ventilation cooling motor, which is characterized in that: the motor is of a shaft radial mixed ventilation cooling structure, the coupling method is used for respectively establishing a fluid network model and a thermal network model, the mutual influence between temperature change and fluid flowing state is fully considered, and the motor global temperature rise prediction is realized through the bidirectional coupling of the fluid network model, so that the motor temperature rise prediction method has the advantage of short time required by the thermal network method to predict the temperature rise, and the accuracy of the motor temperature rise prediction can be improved, so that the motor temperature rise prediction method can be used for accurately predicting the temperature rise of the motor with a complex ventilation structure, and provides reference for the follow-up motor designer to predict the temperature rise of the motor.
Description
Technical Field
The invention relates to the technical field of motors, in particular to an analysis method for accurately predicting the temperature rise of a motor with radial mixed ventilation of a shaft by adopting a flow-heat network model for bidirectional coupling solution.
Background
The motor adopting the axial-radial mixed ventilation cooling is relatively uniform in temperature distribution along the axial direction, but the motor has a relatively complex structure, and modeling is relatively difficult when temperature rise prediction is carried out.
The motor temperature rise prediction method mainly comprises a numerical analysis method and a thermal network method: the former has higher calculation accuracy under the condition of accurate model and boundary conditions, but the global model has large modeling difficulty and longer calculation time; the latter calculation requires short time, but the calculation accuracy is limited by factors such as node division, parameter calculation, heat source distribution and the like.
There are few prior studies in which the interplay between temperature changes and fluid flow conditions is considered when predicting motor temperature rise using the present thermal network method. In practice, the temperature rise change affects the fluid flow characteristic parameter, and the change of the fluid flow characteristic parameter affects the convective heat transfer coefficient of the fluid-solid interface, thereby causing the temperature rise change of the motor.
Disclosure of Invention
Aiming at the current situation of a motor temperature rise prediction method, namely the current situation that the calculation accuracy of a numerical analysis method is high but the modeling difficulty of a global model is high, the calculation time of a thermal network method is short and the calculation accuracy is low, the invention provides a fluid thermal network modeling and bidirectional coupling method of a shaft radial hybrid ventilation cooling motor, a fluid network model and a thermal network model are respectively built for the motor, the mutual influence between the temperature change and the fluid flow state is considered, the fluid thermal network bidirectional coupling solution is carried out, and the motor temperature rise prediction is realized.
The method for modeling and bidirectionally coupling the flow heat network of the shaft radial mixed ventilation cooling motor is characterized by comprising the following steps of: the motor is a motor with a shaft radial mixed ventilation cooling structure, the coupling method is to establish a fluid network model and a thermal network model and perform bidirectional coupling solution to realize motor temperature rise prediction, and the specific implementation process of the method is as follows:
step one: based on the continuity equation and Bernoulli equation of the fluid, respectively constructing the fluid network model of the inner and outer wind paths of the vertical shaft radial mixed ventilation cooling structure motor. According to ventilation circuit A i The ventilation structure change of (i is more than or equal to 1 and less than or equal to n, i is a positive integer) divides the ventilation structure change into m ventilation branches (A) i,1 ,A i,2 ,...,A i,m M is a positive integer), will ventilate branch A i,j The pressure head element (centrifugal fan or rotor radial ventilation radial plate) of which (j is a positive integer) is equal to or more than 1 and less than m and is equivalent to a pressure source P i,j Will A i,j The pressure drop loss (along-path pressure drop loss or local pressure drop loss) of the ventilation branch is equivalent to the flow resistance Z i,j The pressure sources P of each branch are arranged according to the cooling gas flow path i,j And flow resistance Z i,j The connections form a motor fluid network model. In the model A i The ventilation circuit satisfies Σp i,j ＝∑Z i, j Q i,j 2 Wherein: q (Q) i,j For ventilation branch A i,j Is a flow rate of (a).
Step two: based on the motor structure and loss distribution, the motor is divided into x temperature rise regions (B 1 ,B 2 ,...,B x X is a positive integer), in the temperature rise region B s (s is not less than 1 and not more than x, s is a positive integer) y nodes (B) s,1 ,B s,2 ,...,B s,y Y is a positive integer), and the temperature rise region B s Is equally distributed to the nodes of the area, node B s,t (t is more than or equal to 1 and less than or equal to y, and t is a positive integer), and the loss obtained by distribution is equivalent to the heat source q of the node s,t The temperature rise change generated by the heat absorption of the cooling gas is equivalent to the equivalent hot-pressing source delta T of the node where the cooling gas is located s,t Node B s,t With node B s',t' Equivalent thermal resistance R s,t-s',t' Connection, connecting nodes according to the heat transfer path of the motor to form a thermal network model of the motor, wherein each node satisfies q s,t ＝∑(T s',t' -T s,t )/R s,t-s',t' Wherein: t (T) s,t Is node B s,t T is the temperature rise of s',t' Is in contact with node B s,t Adjacent node B s',t' Is a temperature rise of (c).
Step three: initial values are assigned to the parameters that change in the model of the streaming thermal network, respectively. Let the temperature rise T of the node where the cooling gas is located s,t =0, i.e. the cooling gas of each ventilation branch is at room temperature and the equivalent hot-pressing source Δt of the node where the cooling gas is located s,t Let flow rate Q of cooling gas of each ventilation branch i,j ＝1。
Step four: and solving the fluid network model. Temperature rise T of node where cooling gas is located s,t Carry over formula Z i,j ＝f(T s,t ,v,d e …) calculate the flow resistance of each ventilation branch, where f () is the calculation function of the flow resistance, v is the flow rate of the cooling gas, d e Equivalent diameter for the ventilation branch; to calculate the flow resistance Z i,j Substituting the values into the fluid network to solve to obtain new flow of cooling gas of each ventilation branchQuantity Q i,j '. Calculating the relative error of flow, if |Q i,j '-Q i,j |/Q i,j If the ratio is less than or equal to 0.5%, executing the seventh step, otherwise executing the fifth step and adding Q i,j ' assign to Q i,j 。
Step five: calculating heat convection resistance R in heat network model s,t-s',t' . Cooling gas flow Q of each ventilation branch calculated in the step four i,j Substituted into formula R s,t-s',t' ＝h(Q i,j ,λ,d e …) calculating heat convection resistance R of solid temperature rise node and fluid temperature rise node s,t-s',t' The coupling of the fluid network to the thermal network is realized, wherein h () is a calculation function of the thermal resistance, and lambda is the thermal conductivity of the cooling fluid.
Step six: and solving a thermal network model. The convection heat exchange thermal resistance R calculated in the fifth step is calculated s,t-s',t' Substituting into a thermal network model to solve to obtain the temperature rise T of each node of the motor s,t ' Heat flow q of the branch where the equivalent heat pressure source is located air-s,t The method comprises the steps of carrying out a first treatment on the surface of the The calculated heat flow q of the branch where the equivalent heat pressing source is located air-s,t And the cooling gas flow Q of each ventilation branch obtained by calculation in the step four i,j Substituting formula DeltaT s,t ＝q air-s,t /c v Q i,j Calculating to obtain new equivalent hot-pressing source delta T of the node where the cooling gas is located s,t ' wherein c v To cool the specific heat capacity of the gas. Calculating the relative error of the equivalent hot-pressing source if delta T s,t '-ΔT s,t |/ΔT s,t And less than or equal to 0.5 percent, executing the fourth step and calculating the T s,t ' assign to T s,t Coupling of the thermal network to the fluid network is achieved, otherwise step six is re-executed and deltat is set s,t ' assign to DeltaT s,t 。
Step seven: preserving the flow rate Q of cooling gas of each ventilation branch i,j Temperature rise T of each node of motor s,t 。
Compared with the prior art, the invention has the beneficial effects that:
1. the invention fully considers the mutual influence between the temperature change and the fluid flowing state, realizes the motor temperature rise prediction through the bidirectional coupling of the fluid heating network, has the advantage of short time required by the temperature rise prediction by a heating network method, can improve the accuracy of the motor temperature rise prediction, and provides reference for the subsequent motor designer to predict the motor temperature rise.
2. According to the motor temperature rise prediction method for shaft radial mixed ventilation, through the bidirectional coupling of the fluid field and the temperature field and the repeated iteration process, the motor temperature rise prediction method can be used for accurately predicting the temperature rise of the motor with a complex ventilation structure.
Drawings
For the purpose of illustrating the method, the invention is described in detail by the following detailed description and drawings.
FIG. 1 is a flow chart of a method of modeling and bi-directional coupling of a flow-thermal network of a shaft radial hybrid ventilated cooling motor;
FIG. 2 is a schematic diagram of a model of an electric motor with shaft radial mixing ventilation;
FIG. 3 is a schematic diagram of the cooling structure of an axial-radial hybrid ventilated motor;
FIG. 4 is a diagram of an internal and external air path fluid network model of an axial radial hybrid ventilated motor;
FIG. 5 is a diagram of the centrifugal fan and rotor web equivalent centrifugal fan p-Q characteristics;
FIG. 6 is a thermal network model diagram of an electric motor with shaft radial hybrid ventilation;
FIG. 7 is a graph of fluid network results after solution of the fluid thermal coupling;
FIG. 8 is a graph of thermal network results after solution of the fluid-thermal coupling;
reference numerals illustrate: 2-1 is a radiator; 2-2 is a cooling pipe; 2-3 is a motor main body; 2-4 are ventilation baffles; 2-5 are guide plates; 2-6 are outer wind path fans; 2-7 are inner air path fans; 4-1 is an external air path fluid network model; 4-2 is a cooler portion of the internal air path fluid network; 4-3 is an end portion of the internal wind path fluid network; 4-4 are motor body portions of the internal air path fluid network.
Detailed Description
For the purposes of clarity, technical solutions and advantages of the present patent, the method described in the present patent will be described in detail by the specific embodiments shown in the drawings, but it should be understood that these descriptions are exemplary and are not intended to limit the scope of the present patent.
Implementation example one:
the fluid network and the thermal network are mutually influenced, and the change of the parameters of a single network model not only can influence the parameters of another network model, but also can generate a trend of counteracting the change of the original parameters on the fluid network and the thermal network, so that the two-way coupling calculation can obtain more accurate prediction of the global temperature rise distribution of the motor.
This example operates in accordance with the flow chart shown in fig. 1, taking as an example an axial radial hybrid ventilation motor that is air-air cooled, as shown in fig. 2. The internal cooling mode of the motor is shaft radial mixed ventilation cooling, centrifugal fans are arranged on two sides of the interior of the motor, an axial ventilating duct is formed between a rotating shaft and the inner diameter of a rotor through supporting webs, a part of iron core of a motor main body is divided into 8 sections by 7 radial ventilating ducts, the iron core sections are numbered from a shaft extending end to a non-shaft extending end according to 1-8, and the radial ventilating ducts are numbered from the shaft extending end to the non-shaft extending end according to 1-7; the motor of this example is fitted with a cooler above it, which is divided into equally spaced 4 sections by 3 ventilation baffles, numbered 1-4 from the shaft extension to the non-shaft extension.
Step one: based on the continuity equation and Bernoulli equation of the fluid, respectively constructing the fluid network model of the inner and outer wind paths of the vertical shaft radial mixed ventilation cooling structure motor. For ease of calculation and making the following assumptions: considering that the inner surface of the stator and the outer surface of the rotor are smooth, and not considering the influence of the slotting on the flow of cooling gas; the cooling gas is considered as incompressible fluid, and an inner air path fluid network model and an outer air path fluid network model of the motor are respectively established according to the ventilation structure of the motor and the circulation path of the cooling gas.
Fig. 3 is a schematic diagram of a cooling structure of an electric motor with shaft radial mixed ventilation cooling, the cooling structure being composed of inner and outer air passages. The inner air path flows through the following loops: end cavity-shaft extension end (non-shaft extension end) -cooler-end cavity; end cavity-air gap-stator radial plenum-cooler-end cavity number 1 (2, 3, …, 7); end cavity-rotor support plenum-1 (2, 3, …, 7) rotor radial plenum-1 (2, 3, …, 7) stator radial plenum-cooler-end cavity. The outer air passage flow loop is as follows: cooler inlet-cooler deflector-cooling tube outlet.
A fluid network model as shown in fig. 4 is established based on the motor ventilation structure and the flow path of the cooling gas. The cooling gas in the inner air path and the cooling gas in the outer air path of the motor have only heat exchange relationship and have no direct relationship in flow, so that the inner air path fluid network and the outer air path fluid network are mutually independent.
The head elements of the fluid network are internal air path fans at both ends and centrifugal fans equivalent to rotor ventilating slot plates. Due to the non-linearity of the flow resistance, the head of the fan needs to determine the operating point on its p-Q characteristic. The p-Q characteristic of a prototype fan and rotor ventilated web equivalent centrifugal fan is shown in fig. 5. In the cooler area, the flow resistance is the along-path loss flow resistance generated by the cooling gas bypassing the radiating pipe and the pipeline bending flow resistance generated by the cooling gas turning; in the end cavity area of the motor, the flow resistance is the sudden shrinkage loss flow resistance of the cooling gas entering the interior of the motor and the sudden expansion loss flow resistance of the cooling gas entering the end cavity; in the motor end region, the flow resistance is the flow resistance of the loss along the path caused by the cooling gas bypassing the rotor end ring and the stator winding end; in the main body region of the motor, each ventilation channel comprises an along-way loss flow resistance, an inlet loss flow resistance and an outlet loss flow resistance are respectively arranged at the air inlet and the air outlet, and a sudden expansion and a sudden contraction loss flow resistance are arranged at the junction of the axial ventilation channel and the radial ventilation channel.
Step two: and establishing a motor thermal network model based on thermodynamic laws and energy conservation principles. To facilitate the calculation and make the following assumptions: the radial heat transfer and the axial heat transfer are not affected each other, and the circumferential heat transfer is ignored; the contact surfaces are well contacted, and the influence of contact thermal resistance is ignored; the loss of each area of the motor is uniformly distributed, and the loss is brought into the cooler through cooling gas, and a motor thermal network model is built according to the heat transfer path of the motor.
The motor is divided into a solid region and a fluid region. The solid region includes: stator yoke area, stator tooth area, stator winding area, rotor copper bar area, rotor tooth area, permanent magnet area, rotor yoke area and pivot area of No. 1-8 iron core section, stator winding area and rotor copper bar area of No. 1-7 radial air duct, stator winding area and rotor copper bar area of both sides tip. The fluid region includes: the air gap area of the No. 1-8 iron core section and the rotor support radial plate ventilating channel area, the No. 1-7 radial ventilating channel area, the two side end cavity areas, the inner area of the No. 1-4 cooler section and the cooling pipe area.
Each region is divided into 1 node except for 2 nodes of each stator winding region. In the solid area node, equivalent conduction thermal resistance is used for connecting between two adjacent temperature rising nodes along the axial direction in the same component and between two adjacent temperature rising nodes along the radial direction in the two contacted components; the adjacent solid component temperature rise nodes are connected with the fluid region temperature rise nodes by equivalent convective heat dissipation thermal resistance; all the active nodes are respectively connected with an independent heat source; the node satisfies q s,t ＝∑(T s',t' -T s,t )/R s,t-s',t' Wherein: t (T) s,t Is B s,t Temperature rise of node T s',t' Is equal to B s,t Node B with adjacent nodes s',t' Is a temperature rise of (c). Wherein, the nodes of the stator yoke region, the stator winding region, the stator tooth region, the rotor copper bar region, the rotor tooth region, the permanent magnet region and the rotor yoke region are active nodes, and the loss of the region where the nodes are distributed is equivalent to the heat source q of the node s,t The temperature rise change generated by the heat absorption of the cooling gas is equivalent to the equivalent hot-pressing source delta T of the node where the cooling gas is located s,t The nodes are connected according to the motor heat transfer path to form a thermal network model of the motor, as shown in fig. 6.
Step three: initial values are assigned to the parameters that change in the model of the streaming thermal network, respectively. Let the temperature rise T of the node where the cooling gas is located s,t =0, i.e. the cooling gas of each ventilation branch is at room temperature and the equivalent hot-pressing source Δt of the node where the cooling gas is located s,t Let flow rate Q of cooling gas of each ventilation branch i,j ＝1；
Step four: and solving the fluid network model. Temperature rise T of node where cooling gas is located s,t Carry over formula Z i,j ＝f(T s,t ,v,d e …) calculate the flow resistance of each ventilation branch, where v is the flow rate of the cooling gas, d e Equivalent diameter for the ventilation branch; to calculate the flow resistance Z i,j Substituting the values into the fluid network to solve to obtain new flow Q of cooling gas of each ventilation branch i,j '. Calculating the relative error of flow, if |Q i,j '-Q i,j |/Q i,j If the ratio is less than or equal to 0.5%, executing the seventh step, otherwise executing the fifth step and adding Q i,j ' assign to Q i,j 。
Step five: calculating heat convection resistance R in heat network model s,t-s',t' . Cooling gas flow Q of each ventilation branch calculated in the step four i,j Substituted into formula R s,t-s',t' ＝h(Q i,j ,λ,d e …) calculating heat convection resistance R of solid temperature rise node and fluid temperature rise node s,t-s',t' Where λ is the thermal conductivity of the cooling fluid.
Step six: and solving a thermal network model. The convection heat exchange thermal resistance R calculated in the fifth step is calculated s,t-s',t' Substituting into a thermal network model to solve to obtain the temperature rise T of each node of the motor s,t ' Heat flow q of the branch where the equivalent heat pressure source is located air-s,t The method comprises the steps of carrying out a first treatment on the surface of the The calculated heat flow q of the branch where the equivalent heat pressing source is located air-s,t And the cooling gas flow Q of each ventilation branch obtained by calculation in the step four i,j Substituting formula DeltaT s,t ＝q air-s,t /c v Q i,j Calculating to obtain new equivalent hot-pressing source delta T of the node where the cooling gas is located s,t ' wherein c v To cool the specific heat capacity of the gas. Calculating the relative error of the equivalent hot-pressing source, if meeting the requirement of delta T s,t '-ΔT s,t |/ΔT s,t And less than or equal to 0.5 percent, executing the fourth step and calculating the T s,t ' assign to T s,t Re-performing the calculation, otherwise re-performing step six and setting ΔT s,t ' assign to DeltaT s,t 。
Step seven: preserving the flow rate Q of cooling gas of each ventilation branch i,j Temperature rise T of each node of motor s,t The motor flow network coupling solution results are shown in fig. 7 and 8.
And (3) carrying out a temperature rise experiment on the prototype, wherein a dynamometer loading method is adopted during the experiment, so that the temperature rise experiment is carried out on the prototype under the working condition of rated load. After the motor stably operates for 30min, temperature rise data is recorded through a pre-embedded PT100 temperature sensor and a data acquisition processor until the motor is stopped. The relative error between the average temperature rise of the winding obtained by coupling of the motor flow heating network and the average temperature rise of the winding obtained by experiment is 3.15%, and the accuracy of the motor for predicting axial-radial mixed ventilation by the method is verified.
While the foregoing examples illustrate the basic logic and operational concepts of the present invention, those skilled in the art will appreciate that the present invention is not limited by the foregoing examples, and that the foregoing examples and description illustrate only the logic and concepts of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (1)
1. The method for modeling and bidirectionally coupling the flow heat network of the shaft radial mixed ventilation cooling motor is characterized by comprising the following steps of: the motor is a motor with a shaft radial mixed ventilation cooling structure, the coupling method is to establish a fluid network model and a thermal network model and perform bidirectional coupling solution to realize motor temperature rise prediction, and the specific implementation process of the method is as follows:
step one: based on a fluid continuity equation and a Bernoulli equation, respectively constructing an inner wind path fluid network model and an outer wind path fluid network model of the vertical shaft radial mixed ventilation cooling structure motor:
dividing the motor ventilation structure into n ventilation loops A according to the cooling gas flow path 1 ，A 2 ，...，A n According to ventilation circuit A i Dividing the ventilation structure change of (a) into m ventilation branches A i,1 ，A i,2 ，...，A i,m ；
Will ventilate branch A i,j Is equivalent to the pressure source P i,j Will ventilate branch A i,j Is equivalent to the flow resistance Z i,j The pressure sources P of each branch are arranged according to the cooling gas flow path i,j And flow resistance Z i,j The connections form a motor fluid network model in which the ventilation circuit A i Satisfy Sigma P i,j ＝∑Z i,j Q i,j 2 Wherein Q is i,j For ventilation branch A i,j I, j and m are positive integers, i is more than or equal to 1 and less than or equal to n, j is more than or equal to 1 and less than or equal to m, the pressure head element is a centrifugal fan or a radial ventilation radial plate of a rotor, and the pressure drop loss is a path pressure drop loss or a local pressure drop loss;
step two: establishing a motor thermal network model based on thermodynamic law and energy conservation principle: based on the motor structural type and loss distribution characteristics, dividing the motor into x temperature rise areas B according to a heat transfer path 1 ，B 2 ，...，B x In the temperature rise region B s Setting y node bs s,1 ，B s,2 ，...，B s,y Will raise the temperature in region B s Is equally distributed to the nodes of the area, node B s,t The loss obtained by distribution is equivalent to the heat source q of the node s,t The temperature rise change generated by the heat absorption of the cooling gas is equivalent to the equivalent hot-pressing source delta T of the node where the cooling gas is located s,t Node B s,t With node B s',t' Equivalent thermal resistance R s,t-s',t' Connection, connecting nodes according to the heat transfer path of the motor to form a thermal network model of the motor, wherein each node satisfies q s,t ＝∑(T s',t' -T s,t )/R s,t-s',t' Wherein: t (T) s,t Is node B s,t T is the temperature rise of s',t' Is in contact with node B s,t Adjacent node B s',t' Wherein x, s, t and y are positive integers, s is more than or equal to 1 and less than or equal to x, t is more than or equal to 1 and less than or equal to y;
step three: the parameters varied in the flow heating network model are respectively given initial values: let the temperature rise T of the node where the cooling gas is located s,t =0, i.e. the cooling gas of each ventilation branch is at room temperature and the equivalent hot-pressing source Δt of the node where the cooling gas is located s,t Let flow rate Q of cooling gas of each ventilation branch i,j ＝1；
Step four: solving a fluid network model: the node where the cooling gas is locatedTemperature rise T of (2) s,t Carry over formula Z i,j ＝f(T s,t ,v,d e …) calculate the flow resistance of each ventilation branch, where f () is the calculation function of the flow resistance, v is the flow rate of the cooling gas, d e Equivalent diameter for the ventilation branch; to calculate the flow resistance Z i,j Substituting the values into the fluid network to solve to obtain new flow Q of cooling gas of each ventilation branch i,j 'A'; calculating the relative error of flow, if |Q i,j '-Q i,j |/Q i,j If the ratio is less than or equal to 0.5%, executing the seventh step, otherwise executing the fifth step and adding Q i,j ' assign to Q i,j ；
Step five: calculating heat convection resistance R in heat network model s,t-s',t' : cooling gas flow Q of each ventilation branch calculated in the step four i,j Substituted into formula R s,t-s',t' ＝h(Q i,j ,λ,d e …) calculating heat convection resistance R of solid temperature rise node and fluid temperature rise node s,t-s',t' The coupling of the fluid network to the thermal network is realized, wherein h () is a calculation function of thermal resistance, and lambda is the thermal conductivity of cooling fluid;
step six: solving a thermal network model: the convection heat exchange thermal resistance R calculated in the fifth step is calculated s,t-s',t' Substituting into a thermal network model to solve to obtain the temperature rise T of each node of the motor s,t ' Heat flow q of the branch where the equivalent heat pressure source is located air-s,t The method comprises the steps of carrying out a first treatment on the surface of the The calculated heat flow q of the branch where the equivalent heat pressing source is located air-s,t And the cooling gas flow Q of each ventilation branch obtained by calculation in the step four i,j Substituting formula DeltaT s,t ＝q air-s,t /c v Q i,j Calculating to obtain new equivalent hot-pressing source delta T of the node where the cooling gas is located s,t ' wherein c v Specific heat capacity for cooling gas; calculating the relative error of the equivalent hot-pressing source if delta T s,t '-ΔT s,t |/ΔT s,t And less than or equal to 0.5 percent, executing the fourth step and calculating the T s,t ' assign to T s,t Coupling of the thermal network to the fluid network is achieved, otherwise step six is re-executed and deltat is set s,t ' assign to DeltaT s,t ；
Step seven: preserving the ventilation branchesFlow rate Q of cooling gas i,j Temperature rise T of each node of motor s,t 。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110682418.5A CN113378395B (en) | 2021-06-21 | 2021-06-21 | Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110682418.5A CN113378395B (en) | 2021-06-21 | 2021-06-21 | Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113378395A CN113378395A (en) | 2021-09-10 |
| CN113378395B true CN113378395B (en) | 2024-01-26 |
Family
ID=77578151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110682418.5A Active CN113378395B (en) | 2021-06-21 | 2021-06-21 | Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113378395B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104537154A (en) * | 2014-12-09 | 2015-04-22 | 天津大学 | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model |
| CN105243249A (en) * | 2015-11-18 | 2016-01-13 | 哈尔滨理工大学 | Three-phase induction motor transient temperature rise calculation model and measuring calculation method |
| CN107291970A (en) * | 2016-12-19 | 2017-10-24 | 沈阳工业大学 | A kind of computational methods in air to air cooling middle size motor temperature field |
| CN107391884A (en) * | 2017-08-28 | 2017-11-24 | 天津大学 | Double remaining permagnetic synchronous motor temperature rise calculation methods based on equivalent thermal network model |
| CN110414074A (en) * | 2019-07-05 | 2019-11-05 | 南京理工大学 | The equivalent thermal network modeling method of composite excitation axial magnetic field flux switch motor |
| CN111324974A (en) * | 2020-03-24 | 2020-06-23 | 北京交通大学 | Optimization method and device of air-cooled generator based on stator tooth internal cooling ventilation structure |
| WO2021119132A1 (en) * | 2019-12-09 | 2021-06-17 | Transport Phenomena Technologies, Llc | Devices, systems, and methods of hardware-in-the-loop simulation of heat exchange in an environment |
-
2021
- 2021-06-21 CN CN202110682418.5A patent/CN113378395B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104537154A (en) * | 2014-12-09 | 2015-04-22 | 天津大学 | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model |
| CN105243249A (en) * | 2015-11-18 | 2016-01-13 | 哈尔滨理工大学 | Three-phase induction motor transient temperature rise calculation model and measuring calculation method |
| CN107291970A (en) * | 2016-12-19 | 2017-10-24 | 沈阳工业大学 | A kind of computational methods in air to air cooling middle size motor temperature field |
| CN107391884A (en) * | 2017-08-28 | 2017-11-24 | 天津大学 | Double remaining permagnetic synchronous motor temperature rise calculation methods based on equivalent thermal network model |
| CN110414074A (en) * | 2019-07-05 | 2019-11-05 | 南京理工大学 | The equivalent thermal network modeling method of composite excitation axial magnetic field flux switch motor |
| WO2021119132A1 (en) * | 2019-12-09 | 2021-06-17 | Transport Phenomena Technologies, Llc | Devices, systems, and methods of hardware-in-the-loop simulation of heat exchange in an environment |
| CN111324974A (en) * | 2020-03-24 | 2020-06-23 | 北京交通大学 | Optimization method and device of air-cooled generator based on stator tooth internal cooling ventilation structure |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113378395A (en) | 2021-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Camilleri et al. | Predicting the temperature and flow distribution in a direct oil-cooled electrical machine with segmented stator | |
| Tong et al. | Totally enclosed self-circulation axial ventilation system design and thermal analysis of a 1.65-MW direct-drive PMSM | |
| CN109101693B (en) | Method for calculating rotor multi-physical field based on air path-fluid-heat transfer cooperative coupling | |
| CN104242549B (en) | Use the cooling of the motor of radially installed fan | |
| Zhou et al. | Ventilation cooling design for a novel 350-MW air-cooled turbo generator | |
| Weili et al. | Influence of the end ventilation structure change on the temperature distribution in the end region of large water–hydrogen–hydrogen cooled turbogenerator | |
| CN109359339A (en) | Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct | |
| CN109139396A (en) | Cooling system and cooling method of wind generating set and wind driven generator | |
| Yan et al. | Thermal analysis of a novel linear oscillating machine based on direct oil-cooling windings | |
| CN106096157B (en) | A kind of axial ventilation permanent magnet synchronous motor Calculation Method of Temperature Field | |
| CN113378395B (en) | Flow heat network modeling and bidirectional coupling method for shaft radial mixed ventilation cooling motor | |
| CN210246545U (en) | New energy automobile super high efficiency motor cooling system loop structure | |
| CN109067073B (en) | Method for optimizing heat dissipation structure of permanent magnet traction motor with totally-enclosed forced air cooling structure | |
| CN116720305A (en) | New energy automobile motor and electric control part internal temperature rise calculation method | |
| Akawung et al. | Thermal analysis of air cooling system for electric machines using lumped parameter and flow resistance network | |
| CN113890234B (en) | Closed motor cooling device with circulation convection between rotor holes | |
| CN115374645A (en) | Permanent magnet synchronous linear motor rotor thermal network modeling method | |
| CN115577652A (en) | Gas shell type cooling structure of turbine box and design method thereof | |
| Liu et al. | Study of three-dimensional temperature field with changed structures of air-cooled turbo generator stator | |
| CN114465388A (en) | Radial multi-parallel wind path inner-cooling type pumped storage generator motor rotor | |
| CN110671924B (en) | Air duct structure, dryer with air duct structure and dryer control method | |
| CN107291970A (en) | A kind of computational methods in air to air cooling middle size motor temperature field | |
| CN109869289B (en) | Ventilation cooling performance configuration method and device for high tower barrel and tower barrel | |
| CN114647260B (en) | Aircraft servo motor active heat control method based on variable density structure micro-channel heat sink performance regulation and control | |
| CN212627373U (en) | Cross internal circulation cooling structure of closed box type motor |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xu Yongming Inventor after: Meng Yang Inventor after: Cao Hengpei Inventor after: Chang Cuncun Inventor after: Pang Songyin Inventor before: Meng Yang Inventor before: Xu Yongming Inventor before: Cao Hengpei Inventor before: Chang Cuncun Inventor before: Pang Songyin |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230223 Address after: 213032 No. 666 Liaohe Road, Xinbei District, Changzhou City, Jiangsu Province Applicant after: CHANGZHOU INSTITUTE OF TECHNOLOGY Address before: 150080 Harbin University of Science and Technology, Harbin Road, Nangang District, Heilongjiang, 52 Applicant before: HARBIN University OF SCIENCE AND TECHNOLOGY |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |