CN116541969B - Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor - Google Patents
Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor Download PDFInfo
- Publication number
- CN116541969B CN116541969B CN202310778210.2A CN202310778210A CN116541969B CN 116541969 B CN116541969 B CN 116541969B CN 202310778210 A CN202310778210 A CN 202310778210A CN 116541969 B CN116541969 B CN 116541969B
- Authority
- CN
- China
- Prior art keywords
- temperature
- dimensional
- thermal resistance
- low
- temperature side
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000004088 simulation Methods 0.000 claims abstract description 27
- 238000009434 installation Methods 0.000 claims description 19
- 230000004907 flux Effects 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 7
- 238000009529 body temperature measurement Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect 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/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention provides a thermal resistance-based temperature field calculation method for a bolt connection part of a gas compressor, which comprises the following steps: acquiring a first high-temperature side-reference temperature and a first low-temperature side-reference temperature in a three-dimensional temperature field; constructing a two-dimensional simulation model of the bolt connection part; acquiring the application thermal resistance of the bolt connection part; evenly distributing application thermal resistance and calculating the initial value of a two-dimensional temperature field of the bolt connection part; calculating the side reference temperature of the second high temperature side and the side reference temperature of the second low temperature side in the two-dimensional temperature field; correcting the application thermal resistance; and adjusting the distribution proportion of the applied thermal resistance, and calculating the two-dimensional temperature field of the bolt connection part. The method can solve the problem that the two-dimensional temperature field and the three-dimensional temperature field are different, and has the advantages of simplicity, convenience, high calculation accuracy and high calculation efficiency of the two-dimensional temperature field.
Description
Technical Field
The invention belongs to the technical field of calculation of temperature fields of aeroengines, and relates to a thermal resistance-based calculation method of a temperature field of a bolt connection part of a gas compressor.
Background
In the field of aero-engines, the temperature field of the compressor has important significance for structural design and strength evaluation of the compressor. Currently, in order to ensure the calculation efficiency, the temperature field of the air compressor is mostly calculated by adopting a two-dimensional calculation method, and in the process of establishing a two-dimensional calculation model, a non-axisymmetric structure, such as a bolt, is generally subjected to simplified deletion treatment. However, as the bolt connection part of the compressor disk has obvious three-dimensional characteristics, after the calculation of the two-dimensional temperature field is simply carried out according to the addition boundary condition of the actual environment, a certain error exists between the calculation result and the three-dimensional temperature field. Especially when the boundary conditions at the two ends of the bolt are greatly different, the actual temperature difference of the installation planes at the two ends of the bolt connection is often larger than the two-dimensional calculation result.
Moreover, in the prior art, a larger temperature difference of a two-dimensional temperature field is usually realized by adding thermal contact resistance to each contact surface of a bolt connection part, but because errors are mainly caused by unique three-dimensional characteristics of the bolt connection part, the thermal resistance value is often inaccurate to calculate by using a thermal contact resistance calculation formula.
Therefore, it is necessary to design a method capable of reducing the difference between two-dimensional and three-dimensional temperature fields, improving the calculation accuracy of the two-dimensional temperature fields, and improving the calculation efficiency.
Disclosure of Invention
The invention discloses a thermal resistance-based temperature field calculation method for a bolt connection part of a gas compressor, which aims to reduce the difference between two-dimensional and three-dimensional temperature fields, improve the calculation precision of the two-dimensional temperature field and improve the calculation efficiency.
The technical scheme for realizing the aim of the invention is as follows: a thermal resistance-based calculation method for a temperature field of a bolt connection part of a gas compressor comprises the following steps:
s1, acquiring a three-dimensional temperature field of a bolt connection part of a compressor rotor, and calculating a first high-temperature side reference temperature and a first low-temperature side reference temperature under the three-dimensional temperature field;
s2, constructing a two-dimensional simulation model of the bolt connection part, wherein the two-dimensional simulation model of the bolt connection part comprises a low-temperature side installation edge, a middle disc and a high-temperature side installation edge, and the contact surface of the middle disc and the low-temperature side installation edge is defined as a first contact surface, and the contact surface of the middle disc and the high-temperature side installation edge is defined as a second contact surface;
s3, calculating the application thermal resistance of the bolt connection part based on the first high-temperature side reference temperature, the first low-temperature side reference temperature and the two-dimensional simulation model of the bolt connection part;
s4, uniformly distributing application thermal resistance to the first contact surface and the second contact surface, and calculating the initial value of the two-dimensional temperature field of the bolt connection part;
s5, obtaining a second high-temperature side reference temperature and a second low-temperature side reference temperature under the two-dimensional temperature field according to the initial value of the two-dimensional temperature field;
s6, calculating a first difference value between the second high-temperature side reference temperature and the first high-temperature side reference temperature and a second difference value between the second low-temperature side reference temperature and the first low-temperature side reference temperature, and acquiring a comparison result according to the first difference value, the second difference value and a preset comparison condition;
s7, correcting the application thermal resistance according to the comparison result, and repeating the steps S4-S6 until the comparison result meets the preset comparison condition;
and S8, adjusting the distribution proportion of the applied thermal resistance based on the corrected two-dimensional temperature field, and calculating the two-dimensional temperature field of the bolt connection part.
Further, in step S3, based on the first high-temperature side reference temperature, the first low-temperature side reference temperature, and the two-dimensional simulation model of the bolting part, an application thermal resistance of the bolting part is calculated, and the method includes:
s31, calculating the heat flux density of the mounting edge at the low temperature side according to the side temperature of the first low temperature side;
s32, calculating the total thermal resistance of the unit area by using the side temperature of the first low-temperature side, the side temperature of the first high-temperature side and the heat flux density;
s33, acquiring total heat conduction resistance of unit areas of a low-temperature side mounting edge, a middle disc and a high-temperature side mounting edge based on a two-dimensional simulation model of the bolt connection part;
s34, calculating the difference value of the total thermal resistance and the total heat conduction thermal resistance as the application thermal resistance of the unit area of the bolt connection part.
Further, in step S31, the calculation formula for calculating the heat flux density of the mounting side at the low temperature side according to the first low temperature side reference temperature is:wherein->Is the heat flux density; h is a heat exchange coefficient in boundary conditions of the simulation calculation model of the bolt connection part; />Is the temperature of the air flow; />For the wall temperature, a first low-temperature side reference temperature is selected；
In step S32, a calculation formula for calculating the total thermal resistance of the unit area by using the first low-temperature side reference temperature, the first high-temperature side reference temperature and the heat flux density is as follows;wherein->For total heat resistance->Is at a first high temperatureSide reference temperature.
Further, in step S5, the calculation formulas of the second high temperature side reference temperature and the second low temperature side reference temperature are:wherein L is a two-dimensional line segment of a high-temperature side bolt mounting surface or a low-temperature side bolt mounting surface in the two-dimensional simulation model of the bolt connection part; />The weighted average temperature of the line segment L, namely the second high-temperature side reference temperature or the second low-temperature side reference temperature; />The weighted average temperature of a line segment L formed by certain two adjacent nodes on the line segment L; />The area of the circular ring corresponding to the line segment l; />Is the area of the circular ring corresponding to the line segment L.
Further, in step S1, a three-dimensional temperature field of the bolt connection portion of the compressor rotor is obtained, and the method includes:
performing a temperature measurement test of the compressor rotor to obtain a three-dimensional temperature field of the bolt connection part;
or, based on the constructed three-dimensional simulation calculation model of the compressor rotor, setting boundary conditions, and calculating a three-dimensional temperature field of the bolt connection part.
The boundary conditions comprise heat exchange coefficient and heat exchange temperature.
Further, in step S1, a first high-temperature side reference temperature and a first low-temperature side reference temperature in the three-dimensional temperature field are calculated, and the method includes:
calculating the average temperature of the bolt mounting surface on the low temperature side of the bolt connection part as the first low temperature side base temperature and the average temperature of the bolt mounting surface on the high temperature side as the first high temperature side base temperature based on the three-dimensional temperature field;
or based on the three-dimensional temperature field, selecting the temperature of the concerned position on the bolt mounting surface of the bolt on the low temperature side as the first low temperature side basic temperature according to the requirement of a user, and selecting the temperature of the concerned position on the bolt mounting surface of the bolt on the high temperature side as the first high temperature side basic temperature.
Further, in step S6, the preset comparison conditions are: the first difference value is less than or equal to the high temperature side judgment threshold value, and the second difference value is less than or equal to the low temperature side judgment threshold value.
Further, in step S6, the comparison result includes:
s61, executing a step S8 when the first difference value and the second difference value meet comparison conditions;
and S62, when the first difference value and the second difference value do not meet the comparison condition, executing step S7.
Further, in step S7, the correction method of applying thermal resistance is as follows:
calculating a difference between the second high-temperature-side reference temperature and the second low-temperature-side reference temperature as a third difference; calculating a difference between the first high-temperature side reference temperature and the first low-temperature side reference temperature as a fourth difference;
when the third difference value is smaller than the fourth difference value, increasing the application thermal resistance; when the third difference value > the fourth difference value, the application heat resistance is reduced.
Further, in step S8, the method includes adjusting the applied thermal resistance distribution ratio based on the corrected two-dimensional temperature field, including:
s81, defining a middle disc reference temperature;
s82, acquiring a second middle disc reference temperature based on the corrected two-dimensional temperature field, and acquiring a first middle disc reference temperature based on the three-dimensional temperature field;
s83, setting a middle disc judgment threshold value, and calculating a difference value between the first middle disc reference temperature and the second middle disc reference temperature as a fifth difference value;
s84, when the fifth difference value is less than or equal to the middle disc judgment threshold value, the corrected two-dimensional temperature field is the two-dimensional temperature field of the bolt connection part; and when the fifth difference value is larger than the middle disc judgment threshold value, adjusting the distribution proportion of the application thermal resistance and recalculating the two-dimensional temperature field.
Compared with the prior art, the invention has the beneficial effects that: the method for calculating the temperature field of the bolt connection part of the gas compressor based on the thermal resistance can solve the problem that the two-dimensional temperature field and the three-dimensional temperature field of the bolt connection part of the gas compressor rotor are different, and has the advantages of simplicity, convenience, high calculation accuracy and high calculation efficiency.
Drawings
In order to more clearly illustrate the technical solution of the embodiments of the present invention, the drawings that are needed in the description of the embodiments will be briefly described.
FIG. 1 is a flow chart of a method for calculating a temperature field of a bolt connection part of a compressor based on thermal resistance;
FIG. 2 is a schematic diagram of a two-dimensional simulation model of a bolting site in an embodiment.
FIG. 3 is a schematic diagram of a computer device in an embodiment;
wherein, 1, the low temperature side installs the limit; 2. an intermediate plate; 3. a high temperature side mounting edge; 4. a first contact surface; 5. a second contact surface; 201. a memory; 202. a processor.
Detailed Description
The invention will be further described with reference to specific embodiments, and advantages and features of the invention will become apparent from the description. These examples are merely exemplary and do not limit the scope of the invention in any way. It will be understood by those skilled in the art that various changes and substitutions of details and forms of the technical solution of the present invention may be made without departing from the spirit and scope of the present invention, but these changes and substitutions fall within the scope of the present invention.
The embodiment discloses a thermal resistance-based method for calculating a temperature field of a bolt connection part of a gas compressor, which is shown in fig. 1, and comprises the following steps:
s1, acquiring a three-dimensional temperature field of a bolt connection part of a compressor rotor, and calculating a first high-temperature side reference temperature and a first low-temperature side reference temperature under the three-dimensional temperature field;
s2, constructing a two-dimensional simulation model of the bolt connection part, wherein the two-dimensional simulation model of the bolt connection part comprises a low-temperature side installation edge, a middle disc and a high-temperature side installation edge, a contact surface of the middle disc and the low-temperature side installation edge is defined as a first contact surface, and a contact surface of the middle disc and the high-temperature side installation edge is defined as a second contact surface.
What should be stated here is: the constructed two-dimensional simulation model of the bolt connection part can be simplified into a multi-layer flat wall one-dimensional heat conduction problem when the radial temperature difference is smaller than the transverse temperature difference, so that the heat conduction in the radial direction is ignored when the following application thermal resistance is calculated.
And S3, calculating the application thermal resistance of the bolting position based on the first high-temperature side reference temperature, the first low-temperature side reference temperature and the two-dimensional simulation model of the bolting position.
And S4, evenly distributing application thermal resistance to the first contact surface and the second contact surface, and calculating the initial value of the two-dimensional temperature field of the bolt connection part.
S5, obtaining the side temperature of the second high-temperature side and the side temperature of the second low-temperature side under the two-dimensional temperature field according to the initial value of the two-dimensional temperature field.
S6, calculating a first difference value between the second high-temperature side reference temperature and the first high-temperature side reference temperature and a second difference value between the second low-temperature side reference temperature and the first low-temperature side reference temperature, and obtaining a comparison result according to the first difference value, the second difference value and a preset comparison condition.
S7, correcting the application thermal resistance according to the comparison result, and repeating the steps S4-S6 until the comparison result meets the preset comparison condition.
And S8, adjusting the distribution proportion of the applied thermal resistance based on the corrected two-dimensional temperature field, and calculating the two-dimensional temperature field of the bolt connection part.
In an alternative embodiment, in the step S1, a three-dimensional temperature field of the bolt connection portion of the compressor rotor is obtained, and the three-dimensional temperature field of the bolt connection portion may be obtained by performing a temperature measurement test of the compressor rotor.
And the three-dimensional temperature field of the bolt connection part can be calculated through the constructed three-dimensional simulation calculation model of the compressor rotor by using FLUENT software or other software under given boundary conditions. The boundary conditions comprise heat exchange coefficient and heat exchange temperature.
Two methods are used for calculating the first high-temperature side basic temperature and the first low-temperature side basic temperature under the three-dimensional temperature field, wherein one method is based on the three-dimensional temperature field, the average temperature of the low-temperature side bolt mounting surface of the bolt connecting part is calculated to be used as the first low-temperature side basic temperature, and the average temperature of the high-temperature side bolt mounting surface is calculated to be used as the first high-temperature side basic temperature.
When the method is used for calculating the side temperature of the first low-temperature sideAnd a first high-temperature side-reference temperature +.>The formula of (2) is:wherein->Is the average temperature of the bolt mounting surface on the low temperature side (i.e., the first low temperature side-reference temperature +.>) Or the average temperature of the high temperature side bolt mounting surface (i.e., first high temperature side reference temperature +.>) A is a bolt mounting surface on the low temperature side or the high temperature side; />The average temperature of the grid unit a is mounted on the bolt at the low temperature side or the high temperature side; />Is the area of grid cell a; />Is a bolt at the low temperature side or the high temperature sideTotal area of the mounting surface.
The other method is based on a three-dimensional temperature field, the temperature of the concerned position on the bolt mounting surface of the bolt on the low temperature side is selected as the first low temperature side basic temperature according to the requirement of a user, and the temperature of the concerned position on the bolt mounting surface of the bolt on the high temperature side is selected as the first high temperature side basic temperature.
In the step S2, the constructed two-dimensional simulation model of the bolt connection part is shown in fig. 2, and the model includes a low-temperature side mounting edge 1, a middle plate 2, a high-temperature side mounting edge 3, a first contact surface 4, and a second contact surface 5;
in an alternative embodiment, in the step S3, assuming that the area of the first contact surface 4 is the same as that of the second contact surface 5 when the thermal resistance calculation is applied to the bolting part, the low temperature side bolt installation surface receives the heat conduction amount transferred from the high temperature side installation surface 3 to the low temperature side installation surface 1 after reaching the heat balance, and the heat exchanged between the low temperature side bolt installation surface and the air flow is equal, the heat flow density is also equal. Specifically, the process for calculating the application thermal resistance of the bolt connection part comprises the following steps:
s31, calculating the heat flux density of the low-temperature side mounting edge 1 according to the first low-temperature side surface temperature;
s32, calculating the total thermal resistance of the unit area by using the side temperature of the first low-temperature side, the side temperature of the first high-temperature side and the heat flux density;
s33, acquiring total heat conduction resistance of unit area of the low-temperature side mounting edge 1, the middle plate 2 and the high-temperature side mounting edge 3 based on a two-dimensional simulation model of the bolt connection part;
s34, calculating the difference value of the total thermal resistance and the total heat conduction thermal resistance as the application thermal resistance of the unit area of the bolt connection part.
In step S31, according to the first low-temperature side reference temperature, the calculation formula for calculating the heat flux density of the low-temperature side mounting edge 1 is:wherein->Is the heat flux density; h is the imitation of the bolt connection partTruly calculating heat exchange coefficients in boundary conditions of the model; />Is the temperature of the air flow; />For the wall temperature, a first low-temperature-side reference temperature is selected here>；
In step S32, the calculation formula for calculating the total thermal resistance per unit area by using the first low-temperature side reference temperature, the first high-temperature side reference temperature, and the heat flux density is as follows:wherein->For total heat resistance->Is the first high temperature side-side temperature.
In step S33, based on the two-dimensional simulation model of the bolting part, the calculation formula for obtaining the total heat conduction and resistance of the unit area of the low temperature side mounting edge 1, the middle plate 2 and the high temperature side mounting edge 3 is as follows:wherein->Is the total heat conduction resistance; />、/>、/>The wall thicknesses of the low temperature side mounting edge 1, the middle disc 2 and the high temperature side mounting edge 3 are respectively in sequence; />、/>、The thermal conductivity of the low temperature side mounting edge 1, the middle disc 2 and the high temperature side mounting edge 3 are sequentially shown, wherein the thermal conductivity is the thermal conductivity of each region at the average temperature.
In step S34, the thermal resistance per unit area is appliedThe calculation formula of (2) is +.>。
Further, in step S5, the calculation formulas of the second high temperature side reference temperature and the second low temperature side reference temperature are:wherein L is a two-dimensional line segment of a high-temperature side bolt mounting surface or a low-temperature side bolt mounting surface in the two-dimensional simulation model of the bolt connection part; />Is the weighted average temperature of line segment L, i.e. the second high-temperature side reference temperature +.>Or a second low-temperature side-based temperature +.>N is the correction times and is an integer more than or equal to 0; />The weighted average temperature of a line segment L formed by certain two adjacent nodes on the line segment L; />For line segment l pairThe area of the corresponding ring; />Is the area of the circular ring corresponding to the line segment L.
In an alternative embodiment, in step S5, the calculation process of the second high temperature side reference temperature and the second low temperature side reference temperature is:
s51, giving a two-dimensional simulation model of the bolt connection part, wherein boundary conditions of the two-dimensional finite element model comprise heat exchange coefficients and heat exchange temperatures;
s52, applying 0.5 times of application thermal resistance to the first contact surface 4 and the second contact surface 5 respectivelyCalculating an initial value of a two-dimensional temperature field;
s53, utilizing the formulaCalculating the second low temperature side mounting side wall surface temperature +.>And a second high temperature side mounting side wall surface temperature +.>。
In an alternative embodiment, in step S6, the preset comparison condition is: the first difference value is less than or equal to the high temperature side judgment threshold value, and the second difference value is less than or equal to the low temperature side judgment threshold value.
From the above, the second difference isThe first difference is->The low temperature side judgment threshold value is +.>The judgment threshold value at the high temperature side is +.>I.e. preset comparison condition +.>。
Further, in step S6, the method includes:
s61, when the first difference value and the second difference value meet comparison conditions, the fact that the application thermal resistance in the step S2 is reasonable and does not need correction is indicated, and the step S8 is executed to adjust the distribution proportion of the application thermal resistance;
and S62, when the first difference value and the second difference value do not meet the comparison condition, the fact that the application thermal resistance in the step S2 is unreasonably required to be corrected is indicated, and the step S7 is executed to correct the application thermal resistance.
In an alternative embodiment, in step S7, the correction method of applying thermal resistance is:
calculating a difference between the second high-temperature-side reference temperature and the second low-temperature-side reference temperature as a third difference; and calculating the difference between the first high-temperature side reference temperature and the first low-temperature side reference temperature as a fourth difference.
Wherein, when the third difference is less than the fourth difference, i.e.When the applied thermal resistance is smaller than the actual applied thermal resistance, when the third difference > the fourth difference, i.e. when +.>The applied thermal resistance is larger than the actual applied thermal resistance, and the applied thermal resistance is obtained by carrying out iterative correction on the heat flux density for n times>，/>For the heat flux density after the nth correction, the calculation formula is as follows:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>For the second high-temperature-side reference value in the n-1 th modified two-dimensional temperature field, is->A second low temperature side reference value in the modified two-dimensional temperature field for the n-1 th time; />Second low-temperature-side reference temperature +.A second low-temperature-side reference temperature obtained for using the n-1 th modified two-dimensional temperature field>The calculated heat flux density is calculated according to the following formula: />；
In an alternative embodiment, in the step S8, the method includes adjusting the applied thermal resistance distribution ratio based on the corrected two-dimensional temperature field, where:
s81, defining a middle disc reference temperature. The temperature of the concerned position on the middle disc can be selected as the middle disc reference temperature according to the user requirement;
s82, acquiring a second middle disk reference temperature based on the corrected two-dimensional temperature fieldAcquiring a first intermediate disk reference temperature +.>；
S83, setting a middle disk judgment thresholdCalculating the difference between the first intermediate disk reference temperature and the second intermediate disk reference temperature as a fifth differenceValue->；
S84, when the fifth difference value is less than or equal to the middle disc judgment threshold value, namelyThe corrected two-dimensional temperature field is the two-dimensional temperature field of the bolt connecting part; when the fifth difference > the middle disc judgment threshold, i.e., +.>Recalculating the two-dimensional temperature field by adjusting the applied thermal resistance distribution ratio so that +.>。
The method for calculating the temperature field of the bolt connection part of the gas compressor based on the thermal resistance can solve the problem that the two-dimensional temperature field and the three-dimensional temperature field of the bolt connection part of the gas compressor rotor are different, and has the advantages of simplicity, convenience, high calculation accuracy and high calculation efficiency.
In this embodiment, a computer device is further provided, as shown in fig. 3, including a memory 201, a processor 202, and a computer program stored in the memory and capable of running on the processor, where the processor implements the thermal resistance based method for calculating a temperature field of a bolt connection portion of a compressor when executing the computer program.
In particular, the computer device may be a computer terminal, a server or similar computing means.
In this embodiment, there is also provided a computer-readable storage medium storing a computer program for executing the above-described thermal resistance-based compressor bolting site temperature field calculation method.
In particular, computer-readable storage media, including both permanent and non-permanent, removable and non-removable media, may be used to implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer-readable storage media include, but are not limited to, phase-change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Disks (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable storage media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, alternatives, and improvements that fall within the spirit and scope of the invention.
Furthermore, it should be understood that although the present disclosure describes embodiments, not every embodiment is provided with a separate embodiment, and that this description is provided for clarity only, and that the disclosure is not limited to the embodiments described in detail below, and that the embodiments described in the examples may be combined as appropriate to form other embodiments that will be apparent to those skilled in the art.
Claims (9)
1. The method for calculating the temperature field of the bolt connection part of the gas compressor based on thermal resistance is characterized by comprising the following steps of:
s1, acquiring a three-dimensional temperature field of a bolt connection part of a compressor rotor, and calculating a first high-temperature side reference temperature and a first low-temperature side reference temperature under the three-dimensional temperature field;
s2, constructing a two-dimensional simulation model of the bolt connection part, wherein the two-dimensional simulation model of the bolt connection part comprises a low-temperature side installation edge, a middle disc and a high-temperature side installation edge, and the contact surface of the middle disc and the low-temperature side installation edge is defined as a first contact surface, and the contact surface of the middle disc and the high-temperature side installation edge is defined as a second contact surface;
s3, calculating the application thermal resistance of the bolting part based on the first high-temperature side reference temperature, the first low-temperature side reference temperature and the two-dimensional simulation model of the bolting part, wherein the application thermal resistance comprises the following components: s31, calculating the heat flux density of the mounting edge at the low temperature side according to the side temperature of the first low temperature side; s32, calculating the total thermal resistance of the unit area by using the side temperature of the first low-temperature side, the side temperature of the first high-temperature side and the heat flux density; s33, acquiring total heat conduction resistance of unit area of the low-temperature side installation edge, the middle disc and the high-temperature side installation edge based on a two-dimensional simulation model of the bolt connection part; s34, calculating the difference value of the total thermal resistance and the total heat conduction thermal resistance as the application thermal resistance of the unit area of the bolt connection part;
s4, uniformly distributing application thermal resistance to the first contact surface and the second contact surface, and calculating the initial value of the two-dimensional temperature field of the bolt connection part;
s5, obtaining a second high-temperature side reference temperature and a second low-temperature side reference temperature under the two-dimensional temperature field according to the initial value of the two-dimensional temperature field;
s6, calculating a first difference value between the second high-temperature side reference temperature and the first high-temperature side reference temperature and a second difference value between the second low-temperature side reference temperature and the first low-temperature side reference temperature, and acquiring a comparison result according to the first difference value, the second difference value and a preset comparison condition;
s7, correcting the application thermal resistance according to the comparison result, and repeating the steps S4-S6 until the comparison result meets the preset comparison condition;
and S8, adjusting the distribution proportion of the applied thermal resistance based on the corrected two-dimensional temperature field, and calculating the two-dimensional temperature field of the bolt connection part.
2. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1, wherein in step S31, a calculation formula for calculating a heat flux density of a mounting side at a low temperature side according to a first low temperature side reference temperature is:wherein->Is the heat flux density; h is a heat exchange coefficient in boundary conditions of the simulation calculation model of the bolt connection part; />Is the temperature of the air flow; />For the wall temperature, a first low-temperature-side reference temperature is selected here>；
In step S32, the calculation formula for calculating the total thermal resistance per unit area by using the first low-temperature side reference temperature, the first high-temperature side reference temperature, and the heat flux density is as follows:wherein->For total heat resistance->Is the first high temperature side-side temperature.
3. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1, wherein in step S5, the calculation formulas of the second high-temperature side reference temperature and the second low-temperature side reference temperature are:wherein L is a two-dimensional line segment of a high-temperature side bolt mounting surface or a low-temperature side bolt mounting surface in the two-dimensional simulation model of the bolt connection part;is the weighted average temperature of line segment L, alsoI.e. the second high temperature side-reference temperature or the second low temperature side-reference temperature; />The weighted average temperature of a line segment L formed by certain two adjacent nodes on the line segment L; />The area of the circular ring corresponding to the line segment l; />Is the area of the circular ring corresponding to the line segment L.
4. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1, wherein in step S1, a three-dimensional temperature field of a bolting part of a rotor of the compressor is obtained, the method comprising:
performing a temperature measurement test of the compressor rotor to obtain a three-dimensional temperature field of the bolt connection part;
or, based on the constructed three-dimensional simulation calculation model of the compressor rotor, setting boundary conditions, and calculating a three-dimensional temperature field of the bolt connection part.
5. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1 or 4, wherein in step S1, a first high-temperature side reference temperature and a first low-temperature side reference temperature in a three-dimensional temperature field are calculated, the method comprising:
calculating the average temperature of the bolt mounting surface on the low temperature side of the bolt connection part as the first low temperature side base temperature and the average temperature of the bolt mounting surface on the high temperature side as the first high temperature side base temperature based on the three-dimensional temperature field;
or based on the three-dimensional temperature field, selecting the temperature of the concerned position on the bolt mounting surface of the bolt on the low temperature side as the first low temperature side basic temperature according to the requirement of a user, and selecting the temperature of the concerned position on the bolt mounting surface of the bolt on the high temperature side as the first high temperature side basic temperature.
6. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1, wherein in step S6, the preset comparison condition is: the first difference value is less than or equal to the high temperature side judgment threshold value, and the second difference value is less than or equal to the low temperature side judgment threshold value.
7. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 6, wherein in step S6, the comparison result includes:
s61, executing a step S8 when the first difference value and the second difference value meet comparison conditions;
and S62, when the first difference value and the second difference value do not meet the comparison condition, executing step S7.
8. The method for calculating a temperature field of a bolting part of a compressor based on thermal resistance according to claim 1, wherein in step S7, the correction method for applying thermal resistance is as follows:
calculating a difference between the second high-temperature-side reference temperature and the second low-temperature-side reference temperature as a third difference; calculating a difference between the first high-temperature side reference temperature and the first low-temperature side reference temperature as a fourth difference;
when the third difference value is smaller than the fourth difference value, increasing the application thermal resistance; when the third difference value > the fourth difference value, the application heat resistance is reduced.
9. The method for calculating a thermal resistance-based temperature field at a bolting position of a compressor according to claim 1, wherein in step S8, the thermal resistance distribution ratio is adjusted based on the corrected two-dimensional temperature field, the method comprising:
s81, defining a middle disc reference temperature;
s82, acquiring a second middle disc reference temperature based on the corrected two-dimensional temperature field, and acquiring a first middle disc reference temperature based on the three-dimensional temperature field;
s83, setting a middle disc judgment threshold value, and calculating a difference value between the first middle disc reference temperature and the second middle disc reference temperature as a fifth difference value;
s84, when the fifth difference value is less than or equal to the middle disc judgment threshold value, the corrected two-dimensional temperature field is the two-dimensional temperature field of the bolt connection part; and when the fifth difference value is larger than the middle disc judgment threshold value, adjusting the distribution proportion of the application thermal resistance and recalculating the two-dimensional temperature field.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310778210.2A CN116541969B (en) | 2023-06-29 | 2023-06-29 | Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310778210.2A CN116541969B (en) | 2023-06-29 | 2023-06-29 | Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116541969A CN116541969A (en) | 2023-08-04 |
| CN116541969B true CN116541969B (en) | 2023-09-19 |
Family
ID=87454503
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310778210.2A Active CN116541969B (en) | 2023-06-29 | 2023-06-29 | Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116541969B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108828003A (en) * | 2018-06-09 | 2018-11-16 | 安徽华兴车辆有限公司 | A kind of solid material thermal contact resistance measuring device and measuring method |
| CN113268886A (en) * | 2021-06-18 | 2021-08-17 | 五凌电力有限公司株溪口水电厂 | Stator winding temperature correction method and device based on temperature field simulation model |
| CN113343532A (en) * | 2021-06-21 | 2021-09-03 | 哈尔滨工业大学 | Method for predicting evolution of transient high-temperature assembly parameters of countersunk head bolt hybrid connection structure |
| CN115563778A (en) * | 2022-10-09 | 2023-01-03 | 北京凌空天行科技有限责任公司 | Three-dimensional temperature field dynamic distribution calculation method in cold-hot wall conversion process |
| CN116297658A (en) * | 2023-03-21 | 2023-06-23 | 中国人民解放军国防科技大学 | Device and method for testing heat flux of transient temperature field based on heat shrinkage sheet |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101574304B1 (en) * | 2008-10-15 | 2015-12-03 | 엔티엔 가부시키가이샤 | Sensor-equipped bearing for wheel |
-
2023
- 2023-06-29 CN CN202310778210.2A patent/CN116541969B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108828003A (en) * | 2018-06-09 | 2018-11-16 | 安徽华兴车辆有限公司 | A kind of solid material thermal contact resistance measuring device and measuring method |
| CN113268886A (en) * | 2021-06-18 | 2021-08-17 | 五凌电力有限公司株溪口水电厂 | Stator winding temperature correction method and device based on temperature field simulation model |
| CN113343532A (en) * | 2021-06-21 | 2021-09-03 | 哈尔滨工业大学 | Method for predicting evolution of transient high-temperature assembly parameters of countersunk head bolt hybrid connection structure |
| CN115563778A (en) * | 2022-10-09 | 2023-01-03 | 北京凌空天行科技有限责任公司 | Three-dimensional temperature field dynamic distribution calculation method in cold-hot wall conversion process |
| CN116297658A (en) * | 2023-03-21 | 2023-06-23 | 中国人民解放军国防科技大学 | Device and method for testing heat flux of transient temperature field based on heat shrinkage sheet |
Non-Patent Citations (8)
| Title |
|---|
| Heat Source Forecast of Ball Screw Drive System Under Actual Working Conditions Based on On-Line Measurement of Temperature Sensors;Zhenjun Li 等;sensors;1-15 * |
| Modling and simulation of sparking in fastening assemblies;L. Chrmartin 等;HAL open science;1-7 * |
| 接触热阻对超声速飞行器结构响应分析的影响;张仡 等;科学技术与工程;20(12);4985-4991 * |
| 新型反应堆压力容器三维瞬态密封分析程序系统研制;杜雪松;中国优秀博士学位论文全文数据库 工程科技Ⅱ辑(第05期);C029-19 * |
| 某联合循环汽轮机高中压缸-螺栓系统有限元建模与多场耦合应力分析;黄超;中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(第02期);C042-1201 * |
| 汽轮机转子热应力在线监测的研究;王光定;中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(03);C039-207 * |
| 涡轮转子无螺栓挡板过渡态热固耦合响应分析;孙海鹤 等;燃气涡轮试验与研究;第33卷(第02期);44-47 * |
| 高温法兰连接系统温度场的有限元分析;喻九阳 等;武汉工程大学学报;第36卷(第10期);31-36 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116541969A (en) | 2023-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Réthoré et al. | A discrete force allocation algorithm for modelling wind turbines in computational fluid dynamics | |
| CN106021189A (en) | Multi-loop cable steady-state temperature rise acquisition method adapting to various boundary conditions | |
| CN108932392B (en) | Transient temperature calculation method based on improved triple reciprocity boundary element method | |
| CN106021676A (en) | Multi-circle cable steady-state temperature rise acquiring method based on transfer matrix | |
| Wang et al. | Solving of Two‐Dimensional Unsteady‐State Heat‐Transfer Inverse Problem Using Finite Difference Method and Model Prediction Control Method | |
| CN116541969B (en) | Thermal resistance-based temperature field calculation method for bolt connection part of gas compressor | |
| CN112345864A (en) | Method and device for detecting current-carrying capacity of overhead transmission line | |
| CN101650760B (en) | Heat dissipation simulation method and device of computer | |
| Chen et al. | Localized MFS for three‐dimensional acoustic inverse problems on complicated domains | |
| CN116384288B (en) | Compressible flow high-resolution numerical simulation method, medium and device | |
| CN116956414A (en) | Method for determining temperature effect of hyperstatic space structure based on indoor and outdoor temperature difference | |
| CN114295320B (en) | Wind measuring point determining method, system and readable storage medium | |
| CN113779742A (en) | Modeling method and device for simplified model of thermal resistance network and chip junction temperature prediction method | |
| CN116259428A (en) | Method and device for determining flow in primary loop of reactor and computer equipment | |
| CN116401974B (en) | Engine cooling performance evaluation method, device and equipment based on multidimensional | |
| CN116522683B (en) | Building energy consumption carbon emission determination method | |
| CN114924883B (en) | Method, device, equipment and readable medium for determining optimal process mapping | |
| CN107247835B (en) | Library capacity calculation precision verification method | |
| CN118168780B (en) | Pneumatic efficiency detection method, equipment, device and medium for rotary machine | |
| CN116187107B (en) | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium | |
| CN114219143B (en) | Global spectrum mode initial field vertical layering arbitrary interval smooth encryption method | |
| CN117079848B (en) | Nuclear power plant primary loop optimal temperature measurement point selection method | |
| CN114923809B (en) | Method and device for determining atmospheric density, storage medium and electronic equipment | |
| CN111815183B (en) | Gas consumption rationality analysis method and device for gas heating system | |
| CN117034672A (en) | Bridge arm reactor busbar parameter optimization method and device and computer equipment |
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 |