CN114372321B - Temperature design method for water interlayer shell of double-screw compressor - Google Patents

Temperature design method for water interlayer shell of double-screw compressor Download PDF

Info

Publication number
CN114372321B
CN114372321B CN202111523459.6A CN202111523459A CN114372321B CN 114372321 B CN114372321 B CN 114372321B CN 202111523459 A CN202111523459 A CN 202111523459A CN 114372321 B CN114372321 B CN 114372321B
Authority
CN
China
Prior art keywords
temperature
rotor
angle
male rotor
shell
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
Application number
CN202111523459.6A
Other languages
Chinese (zh)
Other versions
CN114372321A (en
Inventor
徐明照
袁玮玮
李宁
王辰
张泉明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Qiyao Screw Machinery Co ltd
711th Research Institute of CSIC
Original Assignee
Shanghai Qiyao Screw Machinery Co ltd
711th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Qiyao Screw Machinery Co ltd, 711th Research Institute of CSIC filed Critical Shanghai Qiyao Screw Machinery Co ltd
Priority to CN202111523459.6A priority Critical patent/CN114372321B/en
Publication of CN114372321A publication Critical patent/CN114372321A/en
Application granted granted Critical
Publication of CN114372321B publication Critical patent/CN114372321B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Applications Or Details Of Rotary Compressors (AREA)

Abstract

The invention aims to solve the technical problem of providing a method for designing the temperature of a water interlayer shell of a double-screw compressor, which comprises the steps of calculating the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled and a relation curve equation between the temperature of process medium gas and the temperature drop of the inner wall of the water interlayer shell according to given process medium gas parameters, machine type parameters of the double-screw compressor and cooling water parameters, and subtracting the temperature drop of the inner wall of the water interlayer shell from the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled to obtain the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled. By utilizing the design method, the temperature distribution of the inner wall of the shell of the double-screw compressor before and after the water interlayer is cooled can be quickly calculated according to the given cooling water parameters.

Description

Temperature design method for water interlayer shell of double-screw compressor
Technical Field
The invention relates to a screw compressor technology.
Background
The screw compressor belongs to a positive displacement compressor, and realizes gas compression along with the continuous reduction of the closed volume in the rotating process by forming the closed volume between a male rotor and a female rotor and a shell. In order to achieve high operation efficiency, it is necessary to control the leakage clearances at the positions of the teeth gaps between the male and female rotors, the tip gaps between the rotors and the casing, the gaps between the rotors and the casing exhaust end surface, and the like, within a small range. However, under some working conditions, the intake and exhaust temperatures of the screw compressor generate huge differences, which causes large thermal deformation of the rotor and the casing, and thus the leakage gaps at various positions change irregularly, which causes the reduction of the operation reliability and the operation efficiency of the main engine, and meanwhile, the uncertainty of the gaps at various positions caused by the thermal deformation causes the difficulty in performance prediction during the design of the main engine, and increases the theoretical and actual design deviation.
The discharge temperature of a dry screw compressor depends mainly on the physical properties of the medium and the operating pressure ratio, which can generally be calculated by the following formula:
Figure BDA0003408628570000011
in the formula:
T d discharge temperature of the compressor, K
T s Suction temperature of the compressor, K
ε o External pressure ratio of compressor
m-multiple square process index
When the external pressure is high or the multi-aspect index of gas is high, the exhaust temperature of the screw compressor is high, and the inner wall surface of the shell is of a double-hole structure, the temperatures of the inner wall surface and the inner wall surface are different, so that the thermal expansion amount of the whole surface is uneven, the gap between the rotor and the inner wall surface of the shell is changed, the risk of the rotor rubbing the cylinder at the position with the smaller gap is increased, and the leakage loss at the position with the larger gap is increased, so that the efficiency of the main engine is reduced.
At present, in order to control the thermal deformation of the screw compressor, cooling media such as atomized cooling oil, atomized cooling water and the like are mostly sprayed at the inlet of the screw compressor, so that a good cooling effect is obtained. However, for some processes, the incorporation of any impurity gas is strictly prohibited and the above method is not feasible. For such working conditions, forced convection heat transfer is often used to achieve temperature control, such as a water jacket shell. A layer of shell is added on the outer side of the original cylinder, and circulating cooling water is introduced between the two layers of shells, so that the temperature of the cylinder is reduced. Because the volumes between the teeth of the screw compressor are spirally and axially pushed, the pressure and the temperature are continuously changed in the pushing process, so that the temperature distribution of the wall surface of the cylinder is more complex, and the cooling design of the water interlayer is also difficult.
For the calculation of the cooling effect of the water jacket, there are currently two main methods:
(1) In the experimental method, the temperature distribution condition of the shell of the testing machine is obtained by introducing cooling water with different flow rates to the testing machine under the specified working condition. The main disadvantage is that the test cost is very high, and the time for reaching the steady state balance of the temperature in the experiment is long, so the process needs to be repeated once each time the parameters are adjusted, and the labor, material and time costs are considerable.
(2) And the CFD simulation method is used for simulating the internal flow field of the screw compressor through CFD software and then obtaining the temperature distribution of the shell of the screw compressor through hot fluid-solid coupling. The main defects are long time consumption, a movable grid technology is needed for obtaining a more accurate result due to unsteady flow of gas in the screw compressor, the quality requirement on the grid is high due to the tiny leakage gaps among the volumes among all the teeth, and the calculation result is very easy to disperse in the iterative process. And the calculation process needs to be repeated once when the parameters are adjusted once, so that the whole calculation period is very long, and the project implementation progress is seriously delayed.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for designing the temperature of a water interlayer shell of a double-screw compressor, and by utilizing the method, the temperature distribution of the inner wall of the shell of the double-screw compressor before and after the water interlayer is cooled can be quickly calculated according to given cooling water parameters.
The temperature design method of the water interlayer shell of the double-screw compressor comprises the following steps:
establishing a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the male rotor and the male rotor angle of the air inlet end face of the male rotor and a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the female rotor and the female rotor angle of the air inlet end face of the female rotor according to given process medium gas parameters and model parameters of the double-screw compressor;
obtaining a relational expression between the temperature of any point on the tooth profile of the male rotor section and the male rotor angle of the male rotor section according to the relational expression between the temperature of any point on the tooth profile of the male rotor air inlet end face and the male rotor angle of the male rotor air inlet end face, and obtaining a relational expression between the temperature of any point on the tooth profile of the female rotor section and the female rotor angle of the female rotor section according to the relational expression between the temperature of any point on the tooth profile of the female rotor air inlet end face and the female rotor angle of the female rotor air inlet end face;
obtaining the temperature distribution of the inner wall of the water interlayer shell before cooling the water interlayer according to the rotor angle corresponding to each point on the longitudinal section of the shell vertical to the axis of the rotor, the relational expression between the temperature of any point on the tooth profile of the cross section of the male rotor and the male rotor angle of the cross section of the male rotor and the relational expression between the temperature of any point on the tooth profile of the cross section of the female rotor and the female rotor angle of the cross section of the female rotor;
based on a cooling water convective heat transfer model, calculating the temperature drop of the inner wall of the water interlayer shell when the process medium gas is at different temperatures according to given cooling water parameters and process medium gas parameters to obtain a plurality of discrete point calculation results, and performing curve fitting on the plurality of discrete point calculation results to obtain a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell;
according to a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell, subtracting the temperature drop of the inner wall of the water interlayer shell from the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled to obtain the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled;
and comparing the maximum temperature value of the inner wall of the water interlayer shell after the water interlayer is cooled with a preset temperature threshold, and if the maximum temperature value is lower than the temperature threshold, indicating that the given cooling water parameters are reasonable.
In the method for designing the temperature of the water jacket shell of the double-screw compressor, if the highest value of the temperature of the inner wall of the water jacket shell after the water jacket is cooled is greater than or equal to the temperature threshold value, the flow rate of the cooling water is changed, the temperature distribution of the inner wall of the water jacket shell after the water jacket is cooled at different flow rates is calculated, and a predicted relation curve between the highest value of the temperature of the inner wall of the water jacket shell and the flow rate of the cooling water is obtained.
The invention has at least the following advantages and characteristics:
(1) According to the embodiment of the invention, the temperature distribution of the inner wall surface of the water interlayer shell of the double-screw compressor can be quickly calculated according to the given cooling water parameters, the process medium gas parameters and the double-screw compressor type parameters;
(2) The embodiment of the invention can provide a predicted change curve of the highest value of the temperature of the inner wall of the shell and the amount of cooling water, can greatly reduce the iteration times, and provides guidance for adjusting the parameters of the cooling water, thereby quickly obtaining a target value;
(3) The temperature distribution result of the inner wall of the shell, which is obtained by calculation in the embodiment of the invention, can be used for temperature control design of the double-screw compressor, can also be used for research on thermal deformation of the screw compressor, and provides reference for clearance design and structural design of the screw compressor.
Drawings
Fig. 1 is a schematic flow chart showing a temperature design method of a water jacket shell of a twin-screw compressor according to an embodiment of the present invention.
Fig. 2 shows a schematic representation of the change in tooth space volume of a twin-screw compressor.
Fig. 3 shows a schematic diagram of a cooling water convection heat exchange model according to an embodiment of the present invention.
Detailed Description
The temperature design method of the water interlayer shell of the double-screw compressor comprises the following steps:
establishing a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the male rotor and the male rotor angle of the air inlet end face of the male rotor and a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the female rotor and the female rotor angle of the air inlet end face of the female rotor according to given process medium gas parameters and model parameters of the double-screw compressor;
obtaining a relational expression between the temperature of any point on the tooth profile of the male rotor section and the male rotor angle of the male rotor section according to the relational expression between the temperature of any point on the tooth profile of the male rotor air inlet end face and the male rotor angle of the male rotor air inlet end face, and obtaining a relational expression between the temperature of any point on the tooth profile of the female rotor section and the female rotor angle of the female rotor section according to the relational expression between the temperature of any point on the tooth profile of the female rotor air inlet end face and the female rotor angle of the female rotor air inlet end face;
obtaining the temperature distribution of the inner wall of the water interlayer shell before cooling the water interlayer according to a rotor angle corresponding to each point on the inner wall of the water interlayer shell on a longitudinal section of the shell perpendicular to the axis of the rotor, a relational expression between the temperature of any point on the tooth profile of the section of the male rotor and the angle of the male rotor on the section of the male rotor, and a relational expression between the temperature of any point on the tooth profile of the section of the female rotor and the angle of the female rotor on the section of the female rotor;
based on a cooling water convective heat transfer model, calculating the temperature drop of the inner wall of the water interlayer shell when the process medium gas is at different temperatures according to given cooling water parameters and process medium gas parameters to obtain a plurality of discrete point calculation results, and performing curve fitting on the plurality of discrete point calculation results to obtain a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell;
according to a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell, subtracting the temperature drop of the inner wall of the water interlayer shell from the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled to obtain the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled;
and comparing the maximum temperature value of the inner wall of the water interlayer shell after the water interlayer is cooled with a preset temperature threshold, and if the maximum temperature value is lower than the temperature threshold, indicating that the given cooling water parameters are reasonable.
Further, if the highest value of the temperature of the inner wall of the water jacket shell after the water jacket is cooled is larger than or equal to the temperature threshold, the flow rate of the cooling water is changed, the temperature distribution of the inner wall of the water jacket shell after the water jacket is cooled at different flow rates is calculated, and a prediction relation curve of the highest value of the temperature of the inner wall of the water jacket shell and the flow rate of the cooling water is obtained.
The method for designing the temperature of the water jacket shell of the twin-screw compressor according to the present invention will be described in detail with reference to a specific embodiment. Referring to fig. 1, the method for designing the temperature of the water jacket shell of the twin-screw compressor of this embodiment includes the following steps:
a. given process media gas parameters and screw compressor model parameters were received as shown in table 1.
TABLE 1 Process Medium gas parameters and screw compressor model parameters
Figure BDA0003408628570000051
The process medium gas parameters also comprise medium physical parameters: specific heat at constant pressure C p Specific heat at constant volume C v The kinematic viscosity v, the prandtl number Pr and the Reynolds number Re.
b. Given cooling water parameters are received as shown in table 2.
TABLE 2 Cooling Water parameters
Figure BDA0003408628570000052
c. According to the given process medium gas parameters and the model parameters of the double-screw compressor, a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the male rotor and the male rotor angle of the air inlet end face of the male rotor and a relational expression between the temperature of any point on the tooth profile of the air inlet end face of the female rotor and the female rotor angle of the air inlet end face of the female rotor are established.
Temperature T (alpha) of any point on the tooth profile of the air inlet end surface of the male rotor 1 ) Male rotor angle alpha with male rotor inlet end face 1 The relation of (A) is as follows:
Figure BDA0003408628570000061
temperature T (alpha) of any point on tooth form of air inlet end surface of female rotor 2 ) Female rotor angle alpha with female rotor air inlet end face 2 The relation of (A) is as follows:
Figure BDA0003408628570000062
wherein: t is s Is the suction temperature of the compressor; t is d Is the discharge temperature of the compressor; alpha is alpha 1s Is the axial suction angle of the male rotor; alpha (alpha) ("alpha") 2s Is the axial suction angle of the female rotor;
Figure BDA0003408628570000063
is the torsion angle coefficient; tau is 1z Is the male rotor twist angle; i is the gear ratio of the male rotor to the female rotor; k is an isentropic index, k is,
Figure BDA0003408628570000064
C p is the constant pressure specific heat of the process medium gas, C v Is the constant specific heat of the process medium gas;
Figure BDA0003408628570000065
the angle of rotation of the male rotor at the end of the first stage of the compressor compression process;
Figure BDA0003408628570000066
is the internal compression corner of the male rotor.
The compression process of a twin-screw compressor is generally divided into two stages: the first stage is that the rotor tooth space volume is reduced until the contact line with unchanged shape and length is formed between a pair of tooth space volumes of the female rotor and the male rotor (only when the pressure of the screw compressor is relatively small, the compression process is terminated at the stage); the second phase ends from the first phase until the male rotor has rotated through a male rotor twist angle tau 1z Until now.
The change of the tooth space volume of the twin-screw compressor is shown in FIG. 2, when the female rotor teeth 1 'rotate to the tooth space area A immediately after the female rotor teeth 1' invade the male rotor teeth 1 01 At the position, i.e. the compression start point, and the starting point of the reduction of the tooth space volume, the tooth crest radial line O of the corresponding male rotor tooth 1 1 W and O 1 O 2 The angle between the connecting lines is beta. O is 1 Is the axis of the male rotor, O 2 Is the axis of the female rotor. The angle of rotation of the male rotor at this position being specified to be zero, i.e.
Figure BDA0003408628570000067
(the male rotor corner as referred to herein)
Figure BDA0003408628570000068
I.e. the male rotor angle alpha of the inlet end face of the male rotor as described above 1 )。
Figure BDA0003408628570000071
Figure BDA0003408628570000072
Figure BDA0003408628570000073
In the formula:
r 1 -outer circle radius of male rotor
A-rotor center distance
r 2 -outer circle radius of female rotor
Figure BDA0003408628570000074
The angle of rotation of the male rotor at the end of the first stage
z 1 Male rotor tooth number
i-gear ratio
z 2 Number of female rotor teeth
Then in the first stage
Figure BDA0003408628570000075
Internal volume ratio epsilon of screw compressor v Angle of rotation with male rotor
Figure BDA0003408628570000076
Is close toIs made of
Figure BDA0003408628570000077
In the formula:
Figure BDA0003408628570000078
coefficient of torsion angle
τ 1z Male rotor twist angle
In the second stage (
Figure BDA00034086285700000718
) Internal volume ratio ε v Angle of rotation with male rotor
Figure BDA0003408628570000079
The relationship between is
Figure BDA00034086285700000710
Make the male rotor corner
Figure BDA00034086285700000711
Angle of rotation of female rotor
Figure BDA00034086285700000712
(female rotor angle as used herein)
Figure BDA00034086285700000713
I.e. the female rotor angle alpha of the female rotor inlet end face 2 ) The internal volume ratio epsilon can be obtained by the meshing relation of the male rotor and the female rotor v Angle of rotation with female rotor
Figure BDA00034086285700000714
The relationship between them.
In the first stage
Figure BDA00034086285700000715
Internal volume ratio epsilon v Angle of rotation with female rotor
Figure BDA00034086285700000716
The relationship between is
Figure BDA00034086285700000717
Second stage
Figure BDA0003408628570000081
Internal volume ratio epsilon v Angle of rotation with female rotor
Figure BDA0003408628570000082
The relationship between is
Figure BDA0003408628570000083
The temperature distribution of the inner wall surface of the shell before cooling is calculated, the temperature values are the same along the spiral line according to the spiral characteristics of the rotor, so that the three-dimensional problem can be converted into the two-dimensional problem, and the temperature distribution of the inner wall surface of the shell can be obtained according to the temperature distribution on the air inlet end surface.
d. Obtaining a relational expression between the temperature of any point on the tooth profile of the male rotor section and the male rotor angle of the male rotor section according to the relational expression between the temperature of any point on the tooth profile of the male rotor air inlet end face and the male rotor angle of the male rotor air inlet end face, and obtaining a relational expression between the temperature of any point on the tooth profile of the female rotor section and the female rotor angle of the female rotor section according to the relational expression between the temperature of any point on the tooth profile of the female rotor air inlet end face and the female rotor angle of the female rotor air inlet end face;
suppose that a point on the tooth profile of the male rotor section is F 1 ，F 1 Point with axis of male rotor O 1 Connecting line F between 1 O 1 And O 1 The angle between W (point W is shown in FIG. 2) is the male rotor angle β of the male rotor section 1 On the assumption of tooth profile of the female rotor sectionOne point of is F 2 ，F 2 And the axis O of the female rotor 2 Connecting line F between 2 O 2 And O 2 The included angle between W is the angle beta of the female rotor 2 。
Temperature T (beta) at any point on the tooth profile of the male rotor section 1 ) Male rotor angle beta with male rotor section 1 The relationship between them is:
Figure BDA0003408628570000084
temperature T (beta) of any point on tooth profile of cross section of female rotor 2 ) Angle beta of the female rotor with respect to the cross-section of the female rotor 2 The relationship between them is:
Figure BDA0003408628570000091
in the formula: lambda 1 Is the length-diameter ratio, lambda, of the male rotor 2 Is the aspect ratio of the female rotor; d 1 The diameter of the outer circle of the male rotor, d 2 The diameter of the outer circle of the female rotor; delta l is the axial distance between the cross section of the rotor and the air inlet end face of the rotor; t is s Is the suction temperature of the compressor; t is d Is the discharge temperature of the compressor; alpha is alpha 1s Is the axial suction angle of the male rotor; alpha is alpha 2s Is the axial suction angle of the female rotor;
Figure BDA0003408628570000092
is the torsion angle coefficient; tau. 1z Is the male rotor twist angle; i is the gear ratio of the male rotor to the female rotor; k is an isentropic index, and k is,
Figure BDA0003408628570000093
C p is the constant pressure specific heat of the process medium gas, C v Is the constant specific heat of the process medium gas;
Figure BDA0003408628570000094
the angle of rotation of the male rotor at the end of the first stage of the compressor compression process;
Figure BDA0003408628570000095
is the internal compression corner of the male rotor.
T (. Beta.) described above 1 ) And T (. Beta.) 2 ) Has taken the helical characteristics of the rotor into full account, by introducing the axial distance Δ l between the cross section of the rotor and the inlet end face of the rotor, and integrating it into the independent variable β 1 And beta 2 The temperature of any point on the three-dimensional tooth form of the rotor can be directly calculated without considering further conversion of the spiral characteristics of the rotor (namely, the point on the air inlet end surface is not required to be converted to the point on the three-dimensional tooth form along a spiral line equation), so that the calculation process is greatly simplified, and the calculation speed is improved.
e. And obtaining the temperature distribution of the inner wall of the water interlayer shell before cooling the water interlayer according to the rotor angle corresponding to each point on the inner wall of the water interlayer shell on the longitudinal section of the shell perpendicular to the axis of the rotor, the relational expression between the temperature of any point on the tooth profile of the section of the male rotor and the angle of the male rotor on the section of the male rotor, and the relational expression between the temperature of any point on the tooth profile of the section of the female rotor and the angle of the female rotor on the section of the female rotor.
During the rotation of the rotor, the temperature of any point on the inner wall of the housing is periodically changed, and the position of the rotor shown in fig. 2 is still exemplified as the position of 0 °, and if the point on the inner wall of the housing is on the male rotor side, the angle theta of the male rotor corresponding to the point on the longitudinal section of the housing perpendicular to the rotor axis can be obtained (assuming that the point is K) 1 ，K 1 Point and male rotor axis O 1 A connecting line K between 1 O 1 And O 1 The angle between W is theta), and theta is substituted into the above-mentioned T (beta) 1 ) Among the formulas (in this case beta) 1 Equal to theta), the temperature of the point can be obtained, and the period change interval of the temperature of the point is
Figure BDA0003408628570000101
Taking periodic variation interval
Figure BDA0003408628570000102
OfSubstituting discrete points into T (beta) 1 ) In the formula, a periodic temperature change curve of the point on the inner wall of the shell can be further obtained; if a point on the inner wall of the housing is on the female rotor side, it can be obtained that the point corresponds to a female rotor angle θ (assuming that the point is K) on a longitudinal cross-section of the housing perpendicular to the rotor axis 2 ，K 2 And point female rotor axis O 2 Connecting line K between 2 O 2 And O 2 The angle between W is theta), and theta is substituted into the above-mentioned T (beta) 2 ) Among the formulas (in this case beta) 2 Equal to theta), the temperature of the point can be obtained, and the periodic variation interval of the point is
Figure BDA0003408628570000103
Taking the interval of periodic variation
Figure BDA0003408628570000104
Substituting a plurality of discrete points in (b) into the above-mentioned T (β) 2 ) From the formula, the periodic temperature variation curve of the point on the inner wall of the shell can be further obtained.
f. Based on a cooling water convective heat transfer model, calculating the temperature drop of the inner wall of the water interlayer shell when the process medium gas is at different temperatures according to given cooling water parameters and process medium gas parameters to obtain a plurality of discrete point calculation results, and performing curve fitting on the plurality of discrete point calculation results to obtain a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell.
The process medium gas contacts with the inner wall surface of the shell to carry out convective heat transfer, and the cooling water contacts with the outer wall surface of the shell to carry out convective heat transfer. Due to the irregular structure of the screw compressor housing, the flow of the fluid in the compressor is also complicated, and it is difficult to calculate and obtain an analytic solution to the problem of heat transfer. The convective heat transfer model of the cooling water in this embodiment is shown in fig. 3, and a shell infinitesimal with an angle d γ and an axial length l is taken, as shown by a shaded part in fig. 3, because the angle is very small, the infinitesimal can be approximated as a flat plate, and it is assumed that the temperature of the fluid is kept constant when the fluid flows through the flat plate, the boundary layer effect is ignored, and the heat radiation is not considered.
According to the principles related to heat transfer [3] The correlation formula of the fluid sweep-out full-plate long isothermal flat plate laminar heat exchange is as follows:
Figure BDA0003408628570000105
the correlation formula of the turbulent heat exchange of the fluid sweep full-plate long isothermal flat plate is as follows:
Figure BDA0003408628570000106
in the formula:
Nu l -nussel number
The surface heat transfer coefficient h and the heat exchange amount phi of the flat plate are
Figure BDA0003408628570000107
φ＝h·S·Δt
In the formula:
λ l -thermal conductivity
length of l-plate
S-fluid flow area
Delta t-temperature difference of heat exchange
When the steady state is reached, the temperature of the shell is kept unchanged, namely the heat exchange quantity phi of the medium gas and the inner wall surface of the cylinder g Equal to the heat exchange amount phi between the outer wall surface of the cylinder and the cooling water w Therefore, a plurality of groups of calculation results of the gas temperature and the temperature drop (namely the delta t) of the inner wall of the water interlayer shell are obtained through calculation. And performing curve fitting on the plurality of groups of calculation results of the gas temperature and the shell temperature drop discrete point to obtain a relation curve equation between the gas temperature and the temperature drop of the inner wall of the water interlayer shell.
g. According to a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell, subtracting the temperature drop of the inner wall of the water interlayer shell from the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled to obtain the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled.
h. Comparing the highest temperature value of the inner wall of the water interlayer shell after the water interlayer is cooled with a preset temperature threshold, and if the highest temperature value is lower than the temperature threshold, indicating that the given cooling water parameter is reasonable and meets the design requirement, wherein the current cooling water parameter is a target value; otherwise, judging that the given cooling water parameter is unqualified, changing the flow of the cooling water, calculating the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled when the flow of the cooling water is different, obtaining a prediction relation curve of the highest temperature value of the inner wall of the water interlayer shell and the flow of the cooling water, referring to the prediction relation curve, returning to the step b to adjust the cooling water parameter, and repeating the calculation process until the calculated highest temperature value of the inner wall of the water interlayer shell after the water interlayer is cooled is lower than the temperature threshold.

Claims (7)

1. The temperature design method of the water interlayer shell of the double-screw compressor is characterized by comprising the following steps:
establishing a relational expression between the temperature of any point on the tooth profile of the air inlet end surface of the male rotor and the angle of the male rotor of the air inlet end surface of the male rotor and a relational expression between the temperature of any point on the tooth profile of the air inlet end surface of the female rotor and the angle of the female rotor of the air inlet end surface of the female rotor according to given process medium gas parameters and model parameters of the double-screw compressor;
obtaining a relational expression between the temperature of any point on the tooth profile of the male rotor section and the male rotor angle of the male rotor section according to the relational expression between the temperature of any point on the tooth profile of the male rotor air inlet end face and the male rotor angle of the male rotor air inlet end face, and obtaining a relational expression between the temperature of any point on the tooth profile of the female rotor section and the female rotor angle of the female rotor section according to the relational expression between the temperature of any point on the tooth profile of the female rotor air inlet end face and the female rotor angle of the female rotor air inlet end face;
obtaining the temperature distribution of the inner wall of the water interlayer shell before cooling the water interlayer according to a rotor angle corresponding to each point on the inner wall of the water interlayer shell on a longitudinal section of the shell perpendicular to the axis of the rotor, a relational expression between the temperature of any point on the tooth profile of the section of the male rotor and the angle of the male rotor on the section of the male rotor, and a relational expression between the temperature of any point on the tooth profile of the section of the female rotor and the angle of the female rotor on the section of the female rotor;
based on a cooling water convective heat transfer model, calculating the temperature drop of the inner wall of the water interlayer shell when the process medium gas is at different temperatures according to given cooling water parameters and process medium gas parameters to obtain a plurality of discrete point calculation results, and performing curve fitting on the plurality of discrete point calculation results to obtain a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell;
according to a relation curve equation between the temperature of the process medium gas and the temperature drop of the inner wall of the water interlayer shell, subtracting the temperature drop of the inner wall of the water interlayer shell from the temperature distribution of the inner wall of the water interlayer shell before the water interlayer is cooled to obtain the temperature distribution of the inner wall of the water interlayer shell after the water interlayer is cooled;
and comparing the highest temperature value of the inner wall of the water interlayer shell after the water interlayer is cooled with a preset temperature threshold, wherein if the highest temperature value is lower than the temperature threshold, the given cooling water parameter is reasonable.
2. The method for designing the temperature of the water jacket shell of the twin-screw compressor as recited in claim 1, wherein if the maximum temperature value of the inner wall of the water jacket shell after the cooling of the water jacket is greater than or equal to the temperature threshold value, the flow rate of the cooling water is changed, the temperature distribution of the inner wall of the water jacket shell after the cooling of the water jacket is calculated at different flow rates of the cooling water, and a predicted relationship curve between the maximum temperature value of the inner wall of the water jacket shell and the flow rate of the cooling water is obtained.
3. The method for designing temperature of water jacket shell of twin-screw compressor according to claim 1, wherein temperature T (β) of any point on the profile of the male rotor cross section 1 ) Male rotor angle beta with male rotor section 1 The relationship between them is:
Figure FDA0003408628560000021
negative rotationTemperature T (beta) at any point on tooth profile of sub-section 2 ) Angle beta of the female rotor with respect to the cross-section of the female rotor 2 The relationship between them is:
Figure FDA0003408628560000022
in the formula: lambda 1 Is the length-diameter ratio, lambda, of the male rotor 2 Is the aspect ratio of the female rotor; d 1 The diameter of the outer circle of the male rotor, d 2 The diameter of the outer circle of the female rotor; delta l is the axial distance between the cross section of the rotor and the air inlet end face of the rotor; t is a unit of s Is the suction temperature of the compressor; t is a unit of d Is the discharge temperature of the compressor; alpha is alpha 1s Is the axial suction angle of the male rotor; alpha is alpha 2s Is the axial suction angle of the female rotor;
Figure FDA0003408628560000023
is the torsion angle coefficient; tau. 1z Is the male rotor twist angle; i is the gear ratio of the male rotor to the female rotor; k is an isentropic index, k is,
Figure FDA0003408628560000024
C p is the constant pressure specific heat of the process medium gas, C v The specific heat of the process medium gas with constant volume is used;
Figure FDA0003408628560000025
the angle of the male rotor at the end of the first stage of the compressor compression process;
Figure FDA0003408628560000026
is the internal compression corner of the male rotor.
4. The method for designing temperature of water jacket shell of twin-screw compressor according to claim 1 or 3, wherein temperature T (α) at any point on the profile of the inlet end face of the male rotor 1 ) Male rotor angle alpha with male rotor inlet end face 1 The relation of (A) is as follows:
Figure FDA0003408628560000031
temperature T (alpha) of any point on tooth form of air inlet end surface of female rotor 2 ) Female rotor angle alpha with female rotor inlet end face 2 The relation of (A) is as follows:
Figure FDA0003408628560000032
wherein: t is s Is the suction temperature of the compressor; t is a unit of d Is the discharge temperature of the compressor; alpha is alpha 1s Is the axial suction angle of the male rotor; alpha is alpha 2s Is the axial suction angle of the female rotor;
Figure FDA0003408628560000033
is the torsion angle coefficient; tau. 1z Is the male rotor twist angle; i is the gear ratio of the male rotor to the female rotor; k is an isentropic index, and k is,
Figure FDA0003408628560000034
C p is the constant pressure specific heat of the process medium gas, C v The specific heat of the process medium gas with constant volume is used;
Figure FDA0003408628560000035
the angle of rotation of the male rotor at the end of the first stage of the compressor compression process;
Figure FDA0003408628560000036
is the internal compression corner of the male rotor.
5. The method for designing the temperature of the water jacket shell of the twin-screw compressor according to claim 1, wherein the cooling water parameters include a cooling water temperature and a cooling water flow rate.
6. The twin screw compression of claim 1The temperature design method of the water interlayer shell of the machine is characterized in that the given process medium gas parameters comprise constant pressure specific heat C p Constant volume specific heat C v The system comprises a compressor, a process medium gas, a dynamic viscosity upsilon, a Prandtl number Pr, a Reynolds number Re, a compressor inlet temperature, a compressor inlet pressure, a process medium gas mass flow and an industrial medium gas component.
7. The method for designing the temperature of the water jacket shell of the twin-screw compressor according to claim 1, wherein the model parameters of the twin-screw compressor comprise a model number, an aspect ratio of a male rotor, an aspect ratio of a female rotor and an internal volume ratio.
CN202111523459.6A 2021-12-13 2021-12-13 Temperature design method for water interlayer shell of double-screw compressor Active CN114372321B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111523459.6A CN114372321B (en) 2021-12-13 2021-12-13 Temperature design method for water interlayer shell of double-screw compressor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111523459.6A CN114372321B (en) 2021-12-13 2021-12-13 Temperature design method for water interlayer shell of double-screw compressor

Publications (2)

Publication Number Publication Date
CN114372321A CN114372321A (en) 2022-04-19
CN114372321B true CN114372321B (en) 2023-03-17

Family

ID=81140523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111523459.6A Active CN114372321B (en) 2021-12-13 2021-12-13 Temperature design method for water interlayer shell of double-screw compressor

Country Status (1)

Country Link
CN (1) CN114372321B (en)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1328513C (en) * 2005-11-21 2007-07-25 西安交通大学 Method for designing molded lines of rotor of helical lobe compressor
CN106484999A (en) * 2016-10-11 2017-03-08 江苏辛普森新能源有限公司 A kind of heat source tower heat pump helical-lobe compressor interior volume specific ratio optimization method
CN113407894A (en) * 2021-06-25 2021-09-17 中国船舶重工集团公司第七一一研究所 Screw rotor tooth space gap thermal state variation calculation method

Also Published As

Publication number Publication date
CN114372321A (en) 2022-04-19

Similar Documents

Publication Publication Date Title
Morini et al. Analysis of a scroll machine for micro ORC applications by means of a RE/CFD methodology
Wu et al. Effects of lubricating oil on the performance of a semi-hermetic twin screw refrigeration compressor
Stosic On heat transfer in screw compressors
EP3018349B1 (en) Rotary screw compressor
Shen et al. A study of clearance height on the performance of single-screw expanders in small-scale organic Rankine cycles
CN110610037A (en) Thermodynamic simulation method for counter-flow heat regenerator
Li et al. Simulation of dry screw vacuum pumps based on chamber model and thermal resistance network
Bianchi et al. Numerical CFD simulations on a small-scale ORC expander using a customized grid generation methodology
Wang et al. Study on a novel screw rotor with variable cross-section profiles for twin-screw vacuum pumps
Mustafin et al. Analysis of influence of screw compressor construction parameters and working condition on rotor temperature fields
CN114372321B (en) Temperature design method for water interlayer shell of double-screw compressor
Stošić et al. The design of a twin-screw compressor based on a new rotor profile
CN110374870B (en) Method for arranging screw rotor gaps of double-screw compressor
Zhang et al. Geometric design and performance analysis of conical-rotor screw vacuum pump with adjustable flow field clearances
Lu et al. Performance improvement of single screw compressor by meshing clearance adjustment used in refrigeration system
Andres et al. Cfd simulation of a two stage twin screw compressor including leakage flows and comparison with experimental data
Chen et al. Numerical simulation and design methodology of a novel asymmetric cylinder profile for sliding vane vacuum pumps
Wang et al. Clearance distribution design and thermal deformation analysis of variable-pitch screw rotors for twin-screw vacuum pumps
Hsieh et al. Numerical analysis of performance, rotor temperature distributions, and rotor thermal deformation of an R134a screw compressor
CN110067600B (en) Vortex expander
Spille-Kohoff et al. CFD simulation of a dry scroll vacuum pump with clearances, solid heating and thermal deformation
Zhang et al. Three-dimensional CFD simulation of a Roots blower for the hydrogen circulating pump
Ma et al. CFD simulation and experimental studying of a dry screw vacuum pump
Xin et al. Innovative Claw Rotor Design for Claw Vacuum Pump and its Power Consumption Analysis
Lu et al. Quasi 1D modelling of conical rotary compressors

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
CB02 Change of applicant information

Address after: 201108 Shanghai city Minhang District Huaning Road No. 3111

Applicant after: The 711 Research Institute of China Shipbuilding Corp.

Applicant after: SHANGHAI QIYAO SCREW MACHINERY Co.,Ltd.

Address before: 201108 Shanghai city Minhang District Huaning Road No. 3111

Applicant before: Shanghai Marine Diesel Engine Research Institute

Applicant before: SHANGHAI QIYAO SCREW MACHINERY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant