CN112417624A - Shafting calculation method based on Q factor, calculation terminal and readable storage medium - Google Patents
Shafting calculation method based on Q factor, calculation terminal and readable storage medium Download PDFInfo
- Publication number
- CN112417624A CN112417624A CN202011441832.9A CN202011441832A CN112417624A CN 112417624 A CN112417624 A CN 112417624A CN 202011441832 A CN202011441832 A CN 202011441832A CN 112417624 A CN112417624 A CN 112417624A
- Authority
- CN
- China
- Prior art keywords
- shafting
- factor
- oil film
- rigidity
- diameter
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 20
- 238000013016 damping Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
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
-
- 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/20—Design reuse, reusability analysis or reusability optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Sliding-Contact Bearings (AREA)
Abstract
The invention provides a shafting calculation method based on Q factors, a calculation terminal and a readable storage medium. The method of the invention does not need to meet the avoidance rate, and even if the working rotating speed is the critical rotating speed, the method has no relation, thereby greatly expanding the design space of the shafting calculation method based on the Q factor.
Description
Technical Field
The invention relates to the field of shafting vibration, in particular to a shafting calculation method based on a Q factor, a calculation terminal and a readable storage medium.
Background
In the prior art, the shafting transverse vibration calculation standard is based on the avoidance rate of the critical rotating speed to the working rotating speed, and the avoidance rate of the critical rotating speed is required to be less than 15% and more than 20%, but practice shows that: some rotors do not meet these requirements and still perform well in the field. For example, the critical rotation speed is 97% of the operation rotation speed and the critical rotation speed is 92% of the operation rotation speed, both of which are easy to balance and operate smoothly on site. On the other hand, the opposite example is also provided, the calculation standard of the shafting transverse vibration is met, but the balance is difficult.
In the prior art, except for the rotor characteristic, only the bearing oil film rigidity is considered in the shafting critical rotating speed calculation, and the bearing oil film damping is not considered, so that the bearings have obvious influence on the shafting critical rotating speed calculation under most conditions. Typically, bearings are the primary source of damping, controlling the response of the rotor. The rigidity and the damping of the bearing have great influence on the calculation of the critical rotating speed of the shafting. Sometimes up to 40% or more.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a shafting calculation method based on a Q factor, which is characterized by comprising the following steps:
modeling the rigidity of the rotor;
calculating the rigidity inner and outer diameters of the rotor:
(1) obtaining the length of the stepped shaft section
L＝1/2(D2-D1)
Obtaining effective diameter on shaft section
Deff＝1/2(D2+D1)
D1，D2: is the diameter of the shaft section of two parts;
(2) obtaining the rigidity outer diameter of the wheel disc:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D+D0)
B: the thickness of the wheel disc;
d: the outer diameter of the contact surface of the wheel disc;
D0: the diameter of the intersection point of the lines in the 45-degree influence area;
(3) obtaining the rigidity inner diameter of the wheel disc:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D1+D0)
D1: inner diameter of contact surface of wheel disc
D2: outer diameter of wheel disc contact surface
D0：45°Diameter of intersection of lines in affected area
Calculating the rigidity and the damping of an oil film of the bearing, wherein the oil film between the shaft and the bearing has anisotropic rigidity and damping;
under the linear condition, the elasticity of the oil film is expressed by an oil film rigidity matrix [ K ] and an oil film damping matrix [ C ];
oil film rigidity coefficients Kxx, Kxy, Kyx and Kyy and oil film damping coefficients Cxx, Cxy, Cyx and Cyy are collectively called as oil film dynamic characteristic coefficients and all change along with the rotating speed.
And calculating the Q factor by adopting a transfer matrix method, and evaluating the safety of the calculation result.
The invention also provides a computing terminal for realizing the shafting computing method based on the Q factor, which comprises the following steps:
the memory is used for storing a computer program and a shafting calculation method based on a Q factor;
and the processor is used for executing the computer program and the shafting calculation method based on the Q factor so as to realize the steps of the shafting calculation method based on the Q factor.
The invention also provides a readable storage medium with the shafting calculation method based on the Q factor, and the readable storage medium stores a computer program which is executed by a processor to realize the steps of the shafting calculation method based on the Q factor.
According to the technical scheme, the invention has the following advantages:
the method provided by the invention is based on shafting calculation expansion of a Q factor and can be used for replacing the critical rotating speed avoidance rate. The method of the invention does not need to meet the avoidance rate, and even if the working rotating speed is the critical rotating speed, the method has no relation, thereby greatly expanding the design space of the shafting calculation method based on the Q factor.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings used in the description will be briefly introduced, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a Q factor graph;
FIG. 2 is a schematic of the stiffness outer diameter of the stepped shaft;
FIG. 3 is a schematic representation of the stiffness outer diameter at the wheel disc;
FIG. 4 is a schematic illustration of the stiffness inner diameter at the wheel disc.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides shafting calculation expansion based on a Q factor to replace the critical rotating speed avoidance rate. The method does not need to meet the avoidance rate, and even if the working rotating speed is the critical rotating speed, the shafting calculation method is greatly expanded.
The shafting calculation method based on the Q factor provided by the invention can be realized by electronic hardware, computer software or the combination of the two, and in order to clearly illustrate the interchangeability of hardware and software, the components and steps of each example have been generally described according to functions in the above description. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
Specifically, each rotor has its natural frequency of vibration called the critical rotational speed of stiffness, and the frequency of vibration of the rotor after the support system is taken into account is called the critical rotational speed of elasticity. When bearing oil films are considered, the safety criteria are: the critical rotating speed is kept away from the working rotating speed by below 15% and above 20%. When bearing oil film and damping are considered, the invention contents are as follows: and (3) shafting calculation expansion based on the Q factor.
The safety criterion of the shafting computing method is shown in figure 1 as a Q factor graph which is the relation between the Q factor and the critical rotating speed.
The portion formed by two oblique lines in fig. 1 is referred to as a river, and is qualified when the Q factor falls below the river and is unqualified when the Q factor falls in the river. It can be seen from fig. 1 that the river is decreasing around the operating speed, i.e. the Q factor is more demanding.
Wherein, Table 1 significance of Q factor
As the shafting calculation method based on the Q factor provided by the invention, firstly, the rigidity of the rotor is modeled, and the modeling method of the typical structure of the rotor comprises the following steps: the calculation mode of the inner and outer diameters of the rigidity can adopt a 45-degree influence area line method:
(1) for a stepped shaft segment, as shown in FIG. 2:
length of
L＝1/2(D2-D1)
Effective diameter on shaft section
Deff＝1/2(D2+D1)
D1，D2: two part shaft section diameter
(2) Stiffness outer diameter at the wheel disc, as shown in fig. 3:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D+D0)
B: thickness of wheel disc
D: outer diameter of wheel disc contact surface
D0: diameter of line intersection point of 45-degree influence area
(3) Stiffness bore at the wheel disc, as shown in fig. 4:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D1+D0)
D1: inner diameter of contact surface of wheel disc
D2: outer diameter of wheel disc contact surface
D0: diameter of line intersection point of 45-degree influence area
Bearing oil film stiffness and damping are then calculated, the oil film between the shaft and the bearing having anisotropic stiffness and damping. Under linearization, oil film elasticity can be expressed by an oil film stiffness matrix [ K ] and an oil film damping matrix [ C ].
Oil film stiffness coefficients Kxx, Kxy, Kyx and Kyy and oil film damping coefficients Cxx, Cxy, Cyx and Cyy are collectively called oil film dynamic characteristic coefficients and all of the coefficients are changed along with the rotating speed.
And then, calculating the Q factor by adopting an optimized transfer matrix method, and evaluating the safety of the calculation result. The specific evaluation mode can be set based on the actual shafting requirement environment, the shafting application equipment, the equipment use condition, the equipment implementation function and the like, and the specific evaluation mode is not limited.
The method provided by the invention is applied to different units, for example, the critical rotating speed 3032r/min of a certain steam turbine generator unit shafting is very close to the working rotating speed 3000r/min, the requirement of avoidance rate is not met, a Q factor method is adopted for calculation and examination, and the requirement of safety criterion is met.
TABLE 2 critical speed (r/min) of certain turbo generator set shafting
| mode | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| Critical speed of rotation | 795 | 1426 | 1440 | 1455 | 1699 | 1731 | 2082 | 2877 | 3032 | 3247 |
TABLE 3 certain turbo-generator set shafting Q factor
TABLE 4 shaft vibration data of bearings of certain steam turbine generator set shaft system
Practice proves that a certain steam turbine generator unit runs stably, has a very small vibration value, and is a safe and reliable unit.
Therefore, the shafting calculation expansion of the Q factor provided by the invention has higher universality, economy and safety. Shafting calculation expansion based on the Q factor can be used for replacing the critical rotating speed avoidance rate. In addition, the new method does not need to meet the previous critical rotating speed avoidance rate, and even if the working rotating speed is the critical rotating speed, the method has no relation, so that the shafting calculation space is greatly expanded.
Based on the method provided by the invention, the invention also provides a computing terminal for realizing the shafting computing method based on the Q factor, which comprises the following steps:
the memory is used for storing a computer program and a shafting calculation method based on a Q factor; and the processor is used for executing the computer program and the shafting calculation method based on the Q factor so as to realize the steps of the shafting calculation method based on the Q factor.
Based on the method provided by the invention, the invention also provides a readable storage medium with the shafting calculation method based on the Q factor, and the readable storage medium stores a computer program which is executed by a processor to realize the steps of the shafting calculation method based on the Q factor.
The computing terminal of the shafting computing method based on the Q factor may employ a mobile terminal such as a mobile phone, a smart phone, a notebook computer, a Digital broadcast receiver, a Personal Digital Assistant (PDA), a tablet computer (PAD), and the like, and a fixed terminal such as a Digital TV, a desktop computer, and the like.
The computing terminal of the shafting computing method based on the Q-factor may include a wireless communication unit, an audio/video (a/V) input unit, a user input unit, a sensing unit, an output unit, a memory, an interface unit, a controller, and a power supply unit, etc. It is to be understood that not all illustrated components are required to be implemented. More or fewer components may alternatively be implemented. Elements of the mobile terminal will be described in detail below.
The computing terminal of the shafting calculation method based on Q factor is the units and algorithm steps of the examples described in connection with the embodiments disclosed herein, and can be implemented in electronic hardware, computer software, or a combination of both, and in the above description the components and steps of the examples have been generally described in terms of functions in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The methods to which the present invention relates may write program code for performing the operations of the present disclosure in any combination of one or more programming languages, including an object oriented programming language such as Java, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computing device, partly on the user's device, as a stand-alone software package, partly on the user's computing device and partly on a remote computing device, or entirely on the remote computing device or server. In the case of a remote computing device, the remote computing device may be connected to the user computing device through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computing device (e.g., through the internet using an internet service provider).
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (3)
1. A shafting calculation method based on Q factors is characterized by comprising the following steps:
modeling the rigidity of the rotor;
calculating the rigidity inner and outer diameters of the rotor:
(1) obtaining the length of the stepped shaft section
L＝1/2(D2-D1)
Obtaining effective diameter on shaft section
Deff＝1/2(D2+D1)
D1，D2: is the diameter of the shaft section of two parts;
(2) obtaining the rigidity outer diameter of the wheel disc:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D+D0)
B: the thickness of the wheel disc;
d: the outer diameter of the contact surface of the wheel disc;
D0: the diameter of the intersection point of the lines in the 45-degree influence area;
(3) obtaining the rigidity inner diameter of the wheel disc:
length L is equal to effective diameter of B shaft section
Deff＝1/2(D1+D0)
D1: inner diameter of contact surface of wheel disc
D2: outer diameter of wheel disc contact surface
D0: diameter of line intersection point of 45-degree influence area
Calculating the rigidity and the damping of an oil film of the bearing, wherein the oil film between the shaft and the bearing has anisotropic rigidity and damping;
under the linear condition, the elasticity of the oil film is expressed by an oil film rigidity matrix [ K ] and an oil film damping matrix [ C ];
oil film rigidity coefficients Kxx, Kxy, Kyx and Kyy and oil film damping coefficients Cxx, Cxy, Cyx and Cyy are collectively called as oil film dynamic characteristic coefficients and all change along with the rotating speed;
and calculating the Q factor by adopting a transfer matrix method, and evaluating the safety of the calculation result.
2. A computing terminal for realizing a shafting computing method based on a Q factor is characterized by comprising the following steps:
the memory is used for storing a computer program and a shafting calculation method based on a Q factor;
a processor for executing the computer program and the method for shafting calculation based on Q-factors to realize the steps of the method for shafting calculation based on Q-factors as claimed in claim 1.
3. A readable storage medium having a Q-factor based shafting calculation method, characterized in that said readable storage medium has stored thereon a computer program for execution by a processor for carrying out the steps of the Q-factor based shafting calculation method according to claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011441832.9A CN112417624A (en) | 2020-12-11 | 2020-12-11 | Shafting calculation method based on Q factor, calculation terminal and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011441832.9A CN112417624A (en) | 2020-12-11 | 2020-12-11 | Shafting calculation method based on Q factor, calculation terminal and readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112417624A true CN112417624A (en) | 2021-02-26 |
Family
ID=74776704
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011441832.9A Pending CN112417624A (en) | 2020-12-11 | 2020-12-11 | Shafting calculation method based on Q factor, calculation terminal and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112417624A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104376157A (en) * | 2014-11-03 | 2015-02-25 | 哈尔滨汽轮机厂有限责任公司 | Shafting dynamic and static characteristic analysis calculating system of large steam turbine generator unit |
| US20190308608A1 (en) * | 2018-04-10 | 2019-10-10 | XL Hybrids | Dynamic hybrid vehicle system for adjusting motor rotary position |
| CN111730410A (en) * | 2020-07-21 | 2020-10-02 | 湖北文理学院 | Static pressure and dynamic and static pressure main shaft oil film rigidity damping real-time measurement method and device, detection device, storage medium and system |
-
2020
- 2020-12-11 CN CN202011441832.9A patent/CN112417624A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104376157A (en) * | 2014-11-03 | 2015-02-25 | 哈尔滨汽轮机厂有限责任公司 | Shafting dynamic and static characteristic analysis calculating system of large steam turbine generator unit |
| US20190308608A1 (en) * | 2018-04-10 | 2019-10-10 | XL Hybrids | Dynamic hybrid vehicle system for adjusting motor rotary position |
| CN111730410A (en) * | 2020-07-21 | 2020-10-02 | 湖北文理学院 | Static pressure and dynamic and static pressure main shaft oil film rigidity damping real-time measurement method and device, detection device, storage medium and system |
Non-Patent Citations (2)
| Title |
|---|
| 何江南，高进，陈丹: "1000MW核电汽轮机转子扭转刚度模化分析", 《东方汽轮机》 * |
| 申秀兰: "国产首台600MW超临界汽轮发电机组轴系振动特性的研究", 《动力工程》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Babu et al. | Vibration modeling of a rigid rotor supported on the lubricated angular contact ball bearings considering six degrees of freedom and waviness on balls and races | |
| Kim et al. | Bifurcation analysis of a rotor supported by five-pad tilting pad journal bearings using numerical continuation | |
| Brouwer et al. | Combined explicit finite and discrete element methods for rotor bearing dynamic modeling | |
| Chasalevris et al. | A novel semi-analytical method for the dynamics of nonlinear rotor-bearing systems | |
| Lee et al. | Unbalance response of a super-critical rotor supported by foil bearings—comparison with test results© | |
| Greenhill et al. | Iterative determination of squeeze film damper eccentricity for flexible rotor systems | |
| Han et al. | Dynamic analysis of an elastic ring squeeze film damper supported rotor using a semi-analytic method | |
| Imani et al. | Constrained nonlinear model predictive control for centrifugal compressor system surge including piping acoustic using closed coupled valve | |
| Kim et al. | Shooting with deflation algorithm-based nonlinear response and neimark-sacker bifurcation and chaos in floating ring bearing systems | |
| Untaroiu et al. | Balancing of flexible rotors using convex optimization techniques: optimum min-max LMI influence coefficient balancing | |
| Liu et al. | Stochastic dynamic analysis of the rotor–bearing system considering the randomness of the radial clearance | |
| Zhang et al. | Control of period-doubling and chaos in varying compliance resonances for a ball bearing | |
| Du et al. | Nonlinear dynamic analysis of a rigid rotor supported by a spiral-grooved opposed-hemisphere gas bearing | |
| Wang | Theoretical and nonlinear behavior analysis of a flexible rotor supported by a relative short herringbone-grooved gas journal-bearing system | |
| Silva Barbosa et al. | Kriging approach dedicated to represent hydrodynamic bearings | |
| CN112417624A (en) | Shafting calculation method based on Q factor, calculation terminal and readable storage medium | |
| Xu et al. | Application of multigrid for the solution of compressible reynolds equation of gas foil bearings | |
| Li et al. | Nonlinear analysis of stability and unbalanced response on spherical spiral-grooved gas bearings | |
| Zhou et al. | Nonlinear characteristics and stability optimization of rotor-seal-bearing system | |
| Yang et al. | Dynamic Characteristics of Elastic Ring Squeeze Film Damper Oscillated by Bearing Contact Force | |
| Grau et al. | An original definition of the profile of compliant foil journal gas bearings: static and dynamic analysis | |
| Li et al. | Dynamic responses of a wind turbine drivetrain under turbulent wind and voltage disturbance conditions | |
| Chasalevris | An Investigation on the Dynamics of High‐Speed Systems Using Nonlinear Analytical Floating Ring Bearing Models | |
| Liu et al. | Optimal vibration reduction of flexible rotor systems by the virtual bearing method | |
| Miranda et al. | Finite element method applied to the eigenvalue analysis of flexible rotors supported by journal bearings |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210226 |