CN112257288B - Method for calculating resistance of axial center type centrifugal ventilator - Google Patents
Method for calculating resistance of axial center type centrifugal ventilator Download PDFInfo
- Publication number
- CN112257288B CN112257288B CN202011224016.2A CN202011224016A CN112257288B CN 112257288 B CN112257288 B CN 112257288B CN 202011224016 A CN202011224016 A CN 202011224016A CN 112257288 B CN112257288 B CN 112257288B
- Authority
- CN
- China
- Prior art keywords
- resistance
- aperture
- centrifugal ventilator
- unit
- aperture unit
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- 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
-
- 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
Abstract
The application provides a method for calculating the resistance of an axial center type centrifugal ventilator, which comprises the following steps: determining the arrangement of aperture units of the axial centrifugal ventilator according to the structural form of the axial centrifugal ventilator and the flow direction of airflow flowing through the axial centrifugal ventilator, and constructing a resistance characteristic unit network frame; determining the parameter value of the corresponding resistance characteristic unit according to the structural parameter and the operation parameter of the resistance characteristic unit; determining the working environment of the axial center type centrifugal ventilator; and calculating the resistance of the aperture unit and the integral resistance of the axial centrifugal ventilator with a plurality of aperture units according to the parameter values in the working environment. The method provided by the application can be established through a network architecture according to the structural characteristics and the size of the axle center type centrifugal ventilator, the theoretical expression of the resistance of the axle center type centrifugal ventilator is directly obtained, and the calculation method is simple and efficient and is suitable for resistance calculation under all working conditions.
Description
Technical Field
The application belongs to the technical field of an axial center type centrifugal ventilator, and particularly relates to a resistance calculation method of the axial center type centrifugal ventilator.
Background
The axial centrifugal ventilator is the ventilator which is most widely applied in the aircraft engine, and the resistance characteristic of the axial centrifugal ventilator is the most important design index of a lubricating and ventilating system. At present, the method for acquiring the resistance of the axial center type centrifugal ventilator mainly comprises a simulation analysis method and a test fitting method. The simulation analysis method mainly simulates the operation condition of the axial centrifugal ventilator through computer software to obtain the resistance characteristic of the axial centrifugal ventilator, and needs CAD software to carry out modeling, grid classification software to carry out grid classification and CFD software to carry out simulation analysis, so that the resistance of the axial centrifugal ventilator under specific rotating speed and environmental conditions can be obtained; the resistance characteristic of the test fitting method is generally obtained by simulating the use working condition on a tester through an axial center type centrifugal ventilator test piece. The resistance value of the axial centrifugal ventilator under the conditions of fixed rotating speed, flow, pressure and temperature can be obtained by the two methods.
However, the simulation analysis method needs to perform steps such as model establishment, gridding division, simulation analysis and the like, the calculation process is complex and long in period, only discrete points under fixed conditions can be obtained, a large amount of manpower and time are required, and the efficiency is low. In the test fitting method, the resistance value of the ventilator is measured on a tester through a test piece, the test piece, test equipment and test personnel are required to be put into, only discrete resistance points under a fixed condition can be obtained, the resistance characteristic of the ventilator under a fixed working condition is obtained through fitting, the fitting precision of the test piece is closely related to the precision of the test equipment, the skill of the test personnel and the selection of the test points, the cost is high, and the efficiency is low; if the ventilator structure is changed, the test result loses the reference value;
at the scheme design stage of the existing aircraft engine lubricating system, dozens of state points in a typical section of an engine need to be calculated and evaluated, and the scheme of the system needs to be continuously adjusted and improved.
Disclosure of Invention
The object of the present application is to provide a method for calculating the resistance of an axial centrifugal ventilator, which solves or reduces at least one of the problems of the background art.
The technical scheme of the application is as follows: a method for calculating the resistance of an axial centrifugal ventilator comprises the following steps:
determining the arrangement of aperture units of the axial centrifugal ventilator according to the structural form of the axial centrifugal ventilator and the flow direction of airflow flowing through the axial centrifugal ventilator, and constructing a resistance characteristic unit network frame;
determining the parameter value of the corresponding resistance characteristic unit according to the structural parameter and the operation parameter of the resistance characteristic unit;
determining the working environment of the axial center type centrifugal ventilator;
calculating the aperture unit resistance and the integral resistance of the axial centrifugal ventilator with a plurality of aperture units according to the parameter values under the working environment, wherein the aperture unit resistance satisfies the following conditions:
wherein j represents the aperture unit number, Δ P j To the pore unit resistance under the corresponding number, P j To the ambient pressure of the aperture unit under the corresponding number, λ j For the aperture cell velocity coefficient, T, at the corresponding number j For the pore diameter unit static temperature, T, under the corresponding numbering j * Is the total temperature of the aperture unit under the corresponding number, R is the gas constant, k is the gas constant, χ j The aperture unit rotation correction factor f under the corresponding number 1j For the correction factor f for the entrance to the aperture unit under the corresponding number 2j Correcting coefficients for the aperture unit outlets under the corresponding numbers;
integral resistance delta P of axial centrifugal ventilator z Satisfy the requirement of
n is the number of the aperture units, and j is more than or equal to 1 and less than or equal to n.
In the present application, the rotation correction coefficient χ corresponding to the numbered aperture unit j Satisfies the following conditions:
if the aperture unit is at the inlet of the axial centrifugal ventilator, the speed v of the fluid entering the aperture unit is vertical to the direction of the tangential speed u, then the rotation correction coefficient chi j According to the velocity coefficient beta uv Curve of the rotation correction coefficient χTaking values, wherein the velocity coefficient beta uv Is the ratio of the fluid inlet flow velocity v to the unit rotational tangential velocity u;
otherwise, rotating the correction coefficient χ j 1 is taken.
In the present application, the inlet correction factor f of the corresponding numbered aperture unit 1j The algorithm of (1) is as follows:
in the formula, S 1j Fluid inlet area of the aperture unit under the corresponding number, S 2j To the area of the aperture unit under the corresponding number, S 3j Let τ be the inlet fill factor and η the inlet relaxation factor for the fluid outlet area of the pore size unit under the corresponding numbering.
In the present application, the inlet correction factor f of the corresponding numbered aperture unit 2j The algorithm is as follows:
in the formula, S 1j Fluid inlet area of the aperture unit under the corresponding number, S 3j The fluid exit area of the aperture unit under the corresponding number.
In the present application, the value of the gas constant k is 1.4.
On the other hand, the technical scheme provided by the application is as follows: an axial centrifugal ventilator, the aperture unit resistance of the axial centrifugal ventilator and the overall resistance of the axial centrifugal ventilator with a plurality of aperture units are obtained by any one of the methods described above.
According to the method provided by the application, the network architecture can be built according to the structural characteristics and the size of the axial centrifugal ventilator, the resistance theoretical expression of the axial centrifugal ventilator is directly obtained, and the calculation method is simple and efficient; the method can be suitable for resistance calculation under all working conditions, the resistance value of the axial ventilator under any condition can be obtained, and when the structure changes, related structure or characteristic parameters can be directly modified to obtain new resistance characteristics; finally, the method can adjust and improve the resistance characteristic according to the performance requirement in the design process of the lubricating system, is suitable for multi-point simulation calculation and evaluation of the system, has strong applicability, and greatly shortens the calculation cost and the calculation period.
Drawings
In order to more clearly illustrate the technical solutions provided in the present application, the drawings will be briefly described below. It is to be understood that the drawings described below are merely exemplary of some embodiments of the application.
FIG. 1 is a flow chart of a method for calculating the resistance of an axial centrifugal ventilator according to the present application.
FIG. 2 is a schematic view of a hub-type centrifugal ventilator aperture unit in the present application.
FIG. 3 is a graph of rotation correction factor versus speed factor according to an embodiment of the present application.
FIG. 4 is a graph of inlet fill factor versus aspect ratio for an embodiment of the present application.
FIG. 5 is a graph of inlet relaxation factor versus fillet ratio for an embodiment of the present application.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present application clearer, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the drawings in the embodiments of the present application.
The method for calculating the resistance of the axial centrifugal ventilator mainly comprises the steps of calculating the resistance of the aperture unit and calculating the total resistance (after the resistance of a plurality of aperture units) of the axial centrifugal ventilator.
As shown in FIG. 1, the method for calculating the resistance of the axial centrifugal ventilator comprises the following steps:
step 1: building resistance feature unit network architecture
According to the structural form of the axial centrifugal ventilator, the arrangement and connection modes of the aperture units are determined from bottom to top or from top to bottom according to the airflow direction, and a resistance characteristic unit network framework is built.
Referring to fig. 1, which shows a typical axial centrifugal ventilator structure according to an embodiment of the present application, an aperture unit 1 is mainly generated by an aperture channel having a circular shape or a certain regular shape, and a flow channel of a fluid 2 is subjected to a sudden expansion and a sudden contraction when passing through the aperture unit 1. The inlet flow passage area 3 of the aperture unit is S 1 The aperture unit flow passage area 4 is S 2 The exit area 5 of the aperture unit is S 3 The aperture unit length 6 is l.
And 2, step: resistance feature cell structure size parameter
Assigning values to the resistance unit models according to the structural parameters and the operating parameters of each resistance unit, wherein the structural parameters comprise the area S of the inlet flow passage 1 Outlet flow passage area S 3 The area of the aperture unit is S 2 Length l, the operating parameters include angular velocity ω.
And 3, step 3: defining the working environment of a centrifugal ventilator
According to the operating environment of the axial centrifugal ventilator, the environmental pressure P, the temperature T and the fluid mass flow Q of the axial centrifugal ventilator are defined m 。
And 4, step 4: calculating the aperture unit resistance
The resistance calculation method of the aperture unit is as follows:
in the formula: j represents an aperture cell number; delta P i The resistance of the aperture unit under the corresponding aperture number; p j The environmental pressure of the aperture unit under the corresponding aperture number is obtained; lambda [ alpha ] j The aperture unit velocity coefficient under the corresponding aperture number; t is j The pore diameter unit static temperature under the corresponding pore diameter number; t is j * The total temperature of the aperture unit under the corresponding aperture number; r is a gas constant; k is a gas constant, which typically takes a value of 1.4; f. of 1j The correction coefficient of the inlet of the aperture unit under the corresponding aperture number is obtained; f. of 2j And the correction coefficient of the outlet of the aperture unit under the corresponding aperture number is shown. Chi-type food processing machine j The aperture under the corresponding aperture numberA unit rotation correction factor chi if the aperture unit is at the inlet of the axial centrifugal ventilator and the velocity v of the fluid entering the aperture unit is perpendicular to the tangential velocity u thereof j Is shown in figure 2 as the rotational correction coefficient χ and velocity coefficient β uv Relation curve (speed coefficient curve beta) uv The ratio of the fluid inlet flow velocity to the unit rotational tangential velocity), otherwise the rotational correction factor χ j 1 is taken.
The calculation method of the inlet and outlet correction coefficients comprises the following steps:
in the formula, S 1j Fluid inlet areas of correspondingly numbered aperture units, S 2j Numbering the area of the aperture unit correspondingly, S 3j Corresponding values for the inlet filling coefficient τ and the inlet relaxation coefficient η for the fluid outlet area of the corresponding numbered pore size unit are given in the handbook of hydromechanical resistances and friction, n.e. ikeliqike (1985), which can be seen in the curves given in fig. 4 and 5 for the inlet filling coefficient τ and the length to diameter ratio and for the inlet relaxation coefficient η and the fillet ratio.
And 5: calculating the whole resistance (total resistance) of the axle center type centrifugal ventilator
The integral resistance of the axial centrifugal ventilator is the sum of the resistances of a plurality of aperture units, and the expression is as follows:
in the formula,. DELTA.P z N is the number of aperture units, and j is more than or equal to 1 and less than or equal to n.
According to the calculation method of the axial centrifugal ventilator, a network framework can be constructed according to the structural characteristics and the size of the axial centrifugal ventilator, the resistance theoretical expression of the axial centrifugal ventilator is directly obtained, and the calculation method is simple and efficient; the method is a universal resistance algorithm of the axial ventilator, is suitable for resistance calculation under all working conditions, can obtain the resistance value of the axial ventilator under any condition, and can directly modify related structures or characteristic parameters when the structures change to obtain new resistance characteristics; finally, in the design process of the lubricating system, the method can adjust and improve the resistance characteristic according to the performance requirement, is suitable for multi-point simulation calculation and evaluation of the system, has strong applicability, and greatly shortens the calculation cost and the calculation period.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present application should be covered within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (6)
1. A method for calculating the resistance of an axial centrifugal ventilator is characterized by comprising the following steps:
determining the arrangement of aperture units of the axial centrifugal ventilator according to the structural form of the axial centrifugal ventilator and the flow direction of airflow flowing through the axial centrifugal ventilator, and constructing a network frame of a resistance characteristic unit;
determining the parameter value of the corresponding resistance characteristic unit according to the structural parameter and the operating parameter of the resistance characteristic unit;
determining the working environment of the axial center type centrifugal ventilator;
calculating the aperture unit resistance and the integral resistance of the axial centrifugal ventilator with a plurality of aperture units according to the parameter values under the working environment, wherein the aperture unit resistance satisfies the following conditions:
wherein j represents an aperture cell number, Δ P j To the pore unit resistance under the corresponding number, P j To the ambient pressure of the aperture unit under the corresponding number, λ j For the aperture cell velocity coefficient, T, at the corresponding number j Pore diameter Unit static temperature, T, at corresponding number j * The total temperature of the pore diameter unit under the corresponding number, R is a gas constant, k is a gas constant, χ j For the aperture unit rotation correction factor under the corresponding number, f 1j For the correction factor f for the entrance to the aperture unit under the corresponding number 2j Correcting coefficients for the aperture unit outlets under the corresponding numbers;
integral resistance delta P of axial centrifugal ventilator z Satisfy the requirements of
n is the number of the aperture units, and j is more than or equal to 1 and less than or equal to n.
2. The method of claim 1, wherein the rotation correction factor χ for the numbered aperture unit is calculated j Satisfies the following conditions:
if the aperture unit is at the inlet of the axial centrifugal ventilator, the speed v of the fluid entering the aperture unit is vertical to the direction of the tangential speed u, then the rotation correction coefficient chi j According to the velocity coefficient beta uv Taking a value with the curve of the rotation correction coefficient chi, wherein the speed coefficient beta uv Is the ratio of the fluid inlet flow velocity v to the unit rotational tangential velocity u;
otherwise, rotating the correction coefficient χ j 1 is taken.
3. The method of claim 1 wherein the inlet correction factor f for a numbered aperture unit is calculated 1j The algorithm of (1) is as follows:
in the formula, S 1j Fluid inlet area of the aperture unit under the corresponding number, S 2j To the area of the aperture unit under the corresponding number, S 3j For the fluid outlet area of the corresponding numbered aperture unit, τ is the inlet fill factor and η is the inlet relaxation factor.
4. A method as claimed in claim 3 wherein the inlet correction factor f for a numbered aperture unit is calculated 2j The algorithm is as follows:
in the formula, S 1j Fluid inlet area of the aperture unit under the corresponding number, S 3j The fluid outlet area of the aperture unit under the corresponding number.
5. The method of claim 1, wherein the gas constant k is 1.4.
6. An axial centrifugal ventilator, wherein the aperture unit resistance of the axial centrifugal ventilator and the overall resistance of the axial centrifugal ventilator having a plurality of aperture units are obtained by the calculation method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011224016.2A CN112257288B (en) | 2020-11-05 | 2020-11-05 | Method for calculating resistance of axial center type centrifugal ventilator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011224016.2A CN112257288B (en) | 2020-11-05 | 2020-11-05 | Method for calculating resistance of axial center type centrifugal ventilator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112257288A CN112257288A (en) | 2021-01-22 |
| CN112257288B true CN112257288B (en) | 2022-08-19 |
Family
ID=74268200
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011224016.2A Active CN112257288B (en) | 2020-11-05 | 2020-11-05 | Method for calculating resistance of axial center type centrifugal ventilator |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112257288B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106446316A (en) * | 2016-05-13 | 2017-02-22 | 中国航空工业集团公司沈阳发动机设计研究所 | Design method for centrifugal ventilator adopting honeycomb structure |
| CN110489829A (en) * | 2019-07-31 | 2019-11-22 | 中国航发沈阳发动机研究所 | A kind of air system Design Method of Fixture Elements based on discharge characteristic |
-
2020
- 2020-11-05 CN CN202011224016.2A patent/CN112257288B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106446316A (en) * | 2016-05-13 | 2017-02-22 | 中国航空工业集团公司沈阳发动机设计研究所 | Design method for centrifugal ventilator adopting honeycomb structure |
| CN110489829A (en) * | 2019-07-31 | 2019-11-22 | 中国航发沈阳发动机研究所 | A kind of air system Design Method of Fixture Elements based on discharge characteristic |
Non-Patent Citations (3)
| Title |
|---|
| 叶轮式通风器分离效率计算方法;张小彬等;《航空动力学报》;20180423(第04期);全文 * |
| 多翼离心风机叶轮的正交设计与试验分析;何立博等;《流体机械》;20160930(第09期);全文 * |
| 金属泡沫离心通风器内流体流动阻力特性研究;陈会敏;《现代制造技术与装备》;20190415(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112257288A (en) | 2021-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102938003B (en) | Method for predicting aeroelasticity stability numerical value of turbomachinery with error frequency included | |
| CN110442833A (en) | One kind assessing Wind turbines health state evaluation method based on various dimensions SCADA data | |
| CN104657596B (en) | A kind of large-scale new compressor performance prediction fast modeling method migrated based on model | |
| CN111859746B (en) | Method for predicting variable working condition performance of turbomachinery based on flow field reconstruction | |
| CN113048086B (en) | Low-noise unequal distance heart fan optimization design method based on radial basis function neural network model | |
| CN109753716A (en) | Core based on flow field simulation/thermal power steam turbine group fluid excitation numerical computation method and system | |
| CN107545110A (en) | A kind of dynamic stress accelerated life test section preparation method | |
| CN106650977A (en) | Short-term power prediction method used for newly-built wind farm | |
| CN113627101A (en) | Wind turbine wake flow simulation method based on improved AD/RSM model | |
| CN116227357A (en) | Wind power plant wake prediction method, model and system thereof | |
| Bühler et al. | Actuator line method simulations for the analysis of wind turbine wakes acting on helicopters | |
| CN112257288B (en) | Method for calculating resistance of axial center type centrifugal ventilator | |
| CN111159941A (en) | Method for simulating transient numerical value of flow field in automobile hydraulic torque converter | |
| Schmitz et al. | Novel performance prediction of a transonic 4.5 stage compressor | |
| CN112326226B (en) | Resistance calculation method for impeller type centrifugal ventilator | |
| JP2019532215A (en) | How to determine the load on a wind turbine | |
| CN114065430B (en) | State data processing method and system for planetary gear box | |
| CN109766611B (en) | Wind farm power simplified prediction method considering terrain gradient | |
| CN116541950A (en) | Method for estimating air performance of mixed exhaust engine by using test data | |
| CN116341407A (en) | Aerodynamic design method for test bed | |
| CN112199907B (en) | Analysis method for multi-index ventilation effect of hydropower station underground powerhouse based on CFD | |
| CN115329415A (en) | Wind vibration coefficient determination method and system for large-span power transmission tower | |
| Pérot et al. | Axial fan noise aeroacoustics predictions and inflow effect on tonal noise using LBM | |
| CN113361058A (en) | Method and device for determining a wind parameter representative of a wind farm | |
| CN117454721B (en) | Wind power plant wake superposition effect evaluation method and medium based on digital simulation experiment |
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 |