CN117454685A - Quantitative analysis method and system for influence of adjustment process parameters on rotor system - Google Patents
Quantitative analysis method and system for influence of adjustment process parameters on rotor system Download PDFInfo
- Publication number
- CN117454685A CN117454685A CN202311219350.2A CN202311219350A CN117454685A CN 117454685 A CN117454685 A CN 117454685A CN 202311219350 A CN202311219350 A CN 202311219350A CN 117454685 A CN117454685 A CN 117454685A
- Authority
- CN
- China
- Prior art keywords
- information
- rotor system
- rolling bearing
- bearing
- calculation function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 86
- 230000008569 process Effects 0.000 title claims abstract description 38
- 238000004445 quantitative analysis Methods 0.000 title claims description 15
- 238000005096 rolling process Methods 0.000 claims abstract description 175
- 238000012360 testing method Methods 0.000 claims abstract description 44
- 238000004458 analytical method Methods 0.000 claims abstract description 42
- 230000007246 mechanism Effects 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 94
- 238000004364 calculation method Methods 0.000 claims description 64
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000004088 simulation Methods 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000008278 dynamic mechanism Effects 0.000 claims description 4
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 4
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 4
- 238000009434 installation Methods 0.000 claims 1
- 230000000750 progressive effect Effects 0.000 claims 1
- 230000004044 response Effects 0.000 abstract description 6
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 15
- 230000000694 effects Effects 0.000 description 12
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000013139 quantization Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000008093 supporting effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005111 flow chemistry technique Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The method comprises the steps of obtaining basic parameter set information of a rolling bearing, constructing a rotor system association model of the rotor system according to the basic parameter set information, and generating analysis data set information of the rotor system according to a test data set and the rotor system association model. The method and the device can fully consider the influence of the bearing adjustment process parameters on the dynamic characteristics of the rotor system, are favorable for accurately determining the influence of the bearing adjustment process parameters such as the adjustment load and the connection load on the critical rotation speed, the dynamic stability and the unbalanced response of the rotor system, can more accurately predict the dynamic characteristics of the rotor system, and provide basis for analyzing the influence mechanism of the bearing adjustment process parameters such as the adjustment load and the connection load on the service performance of the rotor system.
Description
Technical Field
The application relates to the technical field of rotor dynamics, in particular to a quantitative analysis method and a quantitative analysis system for influence of an adjustment process parameter on a rotor system.
Background
Rolling bearing-rotor systems are widely used in various mechanical devices, such as wind power spindles, machine tool spindles, aeroengines; the performance of the rolling bearing-rotor system directly affects the operating efficiency and stability of the overall machine. In practical application, the adjustment process of the rolling bearing can directly influence the adjustment parameters and further directly influence the bearing performance, so that the dynamic characteristics of the rotor system can be predicted better by considering the influence of the adjustment parameters on the dynamic characteristics of the rotor system.
At present, many researches are carried out on the dynamic characteristics of a rolling bearing-rotor system by using a finite element method, but the influence of a bearing adjustment process parameter on the dynamic characteristics of the rotor system is not considered, the influence of the bearing adjustment process parameter on the dynamic characteristics of the rotor system is not beneficial to analysis, and a new rotor system analysis method needs to be provided.
Disclosure of Invention
Based on the above, the embodiment of the application provides a method and a system for quantitatively analyzing the influence of the adjustment process parameters on the rotor system, so as to solve the problem that the influence of the adjustment process parameters of the bearing on the dynamic characteristics of the rotor system is not beneficial to analysis in the prior art.
In a first aspect, an embodiment of the present application provides a method for quantitatively analyzing an influence of an adjustment process parameter on a rotor system, where the rotor system includes a bearing seat, a rolling bearing, and a rotating shaft, where the rolling bearing is mounted on the bearing seat, and the rotating shaft is mounted on the rolling bearing, the method includes:
acquiring basic parameter set information of the rolling bearing, wherein the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first groove bottom diameter information of the inner ring of the rolling bearing, second groove bottom diameter information of the outer ring of the rolling bearing, first poisson ratio information of the bearing seat, second poisson ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing;
constructing a rotor system association model of the rotor system according to the basic parameter set information, wherein the rotor system association model is used for describing a mathematical model for associating the adjustment process parameters with the rotor system;
And generating analysis data set information of the rotor system according to the test data set and the rotor system association model.
Compared with the prior art, the beneficial effects that exist are: according to the quantitative analysis method for the influence of the load adjustment process parameters on the rotor system, terminal equipment can firstly acquire basic parameter set information of the rolling bearing, then construct a rotor system association model of the rotor system according to the basic parameter set information, and then generate analysis data set information of the rotor system according to the test data set and the rotor system association model, so that the influence of the load adjustment process parameters on the dynamic characteristics of the rotor system is fully considered, the dynamic characteristics of the rotor system can be predicted more accurately, a basis is provided for analyzing the influence mechanism of the load adjustment process parameters such as the load adjustment load and the connection load on the service performance of the rotor system, and the problem that the influence of the load adjustment process parameters on the dynamic characteristics of the rotor system is not beneficial to analysis is solved to a certain extent.
In a second aspect, an embodiment of the present application provides a quantitative analysis system for an influence of an adjustment process parameter on a rotor system, the quantitative analysis system being applied to the rotor system, the rotor system including a bearing housing, a rolling bearing and a rotating shaft, the rolling bearing being mounted on the bearing housing, the rotating shaft being mounted on the rolling bearing, the system including:
Basic parameter set information acquisition module: the method comprises the steps of obtaining basic parameter set information of the rolling bearing, wherein the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first groove bottom diameter information of the inner ring of the rolling bearing, second groove bottom diameter information of the outer ring of the rolling bearing, first poisson ratio information of the bearing seat, second poisson ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing;
the rotor system association model building module: the rotor system association model is used for constructing a rotor system association model of the rotor system according to the basic parameter set information, wherein the rotor system association model is used for describing a mathematical model for associating the adjustment process parameters with the rotor system;
an analysis data set information generation module: and the analysis data set information is used for generating the rotor system according to the test data set and the rotor system association model.
In a third aspect, embodiments of the present application provide a terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to the first aspect as described above when the computer program is executed.
In a fourth aspect, embodiments of the present application provide a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method of the first aspect described above.
It will be appreciated that the advantages of the second to fourth aspects may be found in the relevant description of the first aspect and are not repeated here.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings used in the description of the embodiments or the prior art will be briefly described below.
FIG. 1 is a flow chart of a quantitative analysis method according to an embodiment of the present application;
FIG. 2 is a first schematic illustration of a rotor system provided in an embodiment of the present application;
FIG. 3 is a second schematic view of a rotor system provided in an embodiment of the present application;
FIG. 4 is a flowchart of step S200 in a quantitative analysis method according to an embodiment of the present application;
FIG. 5 is a first schematic diagram of a finite element model provided in an embodiment of the present application;
FIG. 6 is a second schematic diagram of a finite element model provided in an embodiment of the present application;
FIG. 7 is a flowchart of step S300 in a quantitative analysis method according to an embodiment of the present application;
FIG. 8 is a vibration mode diagram of a rotor system provided in an embodiment of the present application;
FIG. 9 is a first schematic diagram of a threshold rotational speed result provided by an embodiment of the present application;
FIG. 10 is a second schematic diagram of a threshold rotational speed result provided by an embodiment of the present application;
FIG. 11 is a first schematic diagram of an imbalance response provided by an embodiment of the present disclosure;
FIG. 12 is a second schematic diagram of an imbalance response provided by an embodiment of the present disclosure;
fig. 13 is a schematic flow chart after step S320 in the quantization analysis method according to an embodiment of the present application;
FIG. 14 is a first schematic illustration of stability provided by an embodiment of the present application;
FIG. 15 is a second schematic illustration of stability provided by an embodiment of the present application;
FIG. 16 is a block diagram of a quantization analysis system provided in an embodiment of the present application;
fig. 17 is a schematic diagram of a terminal device according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system configurations, techniques, etc. in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
In the description of this application and the claims that follow, the terms "first," "second," "third," etc. are used merely to distinguish between descriptions and should not be construed to indicate or imply relative importance.
Reference in the specification to "one embodiment" or "some embodiments" or the like means that a particular feature, structure, or characteristic described in connection with the embodiment is included in one or more embodiments of the application. Thus, appearances of the phrases "in one embodiment," "in some embodiments," "in other embodiments," and the like in the specification are not necessarily all referring to the same embodiment, but mean "one or more but not all embodiments" unless expressly specified otherwise. The terms "comprising," "including," "having," and variations thereof mean "including but not limited to," unless expressly specified otherwise.
In order to illustrate the technical solutions described in the present application, the following description is made by specific examples.
Referring to fig. 1, fig. 1 is a flow chart of a method for quantitatively analyzing an influence of a tuning process parameter on a rotor system according to an embodiment of the present disclosure. In this embodiment, the execution subject of the quantization analysis method is a terminal device. It will be appreciated that the types of terminal devices include, but are not limited to, cell phones, tablet computers, notebook computers, ultra-Mobile Personal Computer (UMPC), netbooks, personal digital assistants (Personal Digital Assistant, PDA), etc., and embodiments of the present application do not impose any limitation on the specific type of terminal device.
Referring to fig. 1, the quantitative analysis method provided in the embodiment of the present application includes, but is not limited to, the following steps:
in S100, basic parameter set information of the rolling bearing is acquired.
Without loss of generality, the quantitative analysis method can be applied to rotor systems, i.e. rolling bearing-rotor systems; the rotor system comprises a bearing seat, a rolling bearing and a rotating shaft, wherein the rolling bearing can be arranged on the bearing seat, and the rotating shaft can be arranged on the rolling bearing. For ease of understanding the embodiments of the present application, reference is made to fig. 2 and 3, fig. 2 being a schematic diagram of a model of a rolling bearing-rotor system, and fig. 3 being a schematic diagram of a kinetic model of a rolling bearing-rotor system.
Specifically, the basic parameter set information is used for describing a data set of a plurality of basic parameters of the rolling bearing, the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first bottom diameter information of an inner ring of the rolling bearing, second bottom diameter information of an outer ring of the rolling bearing, first Poisson's ratio information of the bearing seat, second Poisson's ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing, wherein the inner diameter information is used for describing the inner diameter of the rolling bearing, the outer diameter information is used for describing the outer diameter of the inner ring of the rolling bearing and the bearing seat hole of the bearing seat in the diameter direction, the first bottom diameter information of the groove is used for describing the bottom diameter of the inner ring of the rolling bearing seat, the first Poisson's ratio information is used for describing the first elastic modulus of the rolling bearing, and the second Poisson's ratio information is used for describing the elastic modulus of the rolling bearing.
Illustratively, when the basic parameters of the rolling bearing are shown in table 1 below, the inner diameter information of the rolling bearing is 20 mm, the outer diameter information of the rolling bearing is 42 mm, the first groove bottom diameter information is 5.94 mm, and the second groove bottom diameter information is 6.27 mm.
Table 1 basic parameter table of rolling bearing
Bearing inner diameter d/mm | 20 |
Bearing outer diameter D/mm | 42 |
Diameter D of rolling element b /mm | 5.5 |
Diameter of center circle D m /mm | 31 |
Radius r of inner race groove i /mm | 2.970 |
Outer race groove radius r e /mm | 3.135 |
Original contact angle alpha o /° | 15 |
Number of rolling elements Z | 13 |
Specifically, the terminal device may first obtain basic parameter set information of the rolling bearing, where the basic parameter set information may be customized by a user, or may be manually input after the user measures the real rolling bearing.
In S200, a rotor system correlation model of the rotor system is constructed from the basic parameter set information.
Without loss of generality, when parameters such as the structure, the quality and the like of the rotor system are unchanged, main parameters affecting the dynamic characteristics of the rotor system are the rotating speed of the rotating shaft and the rigidity of the rolling bearing, and the rolling bearing is affected by the rotating speed, radial load (such as self gravity of the system), load adjustment, connecting load and other adjustment parameters of the rotating shaft in actual operation, so when the load adjustment and the connecting load are changed, the contact angle of the rolling bearing is changed, thereby changing the rigidity of the rolling bearing and affecting the dynamic performance of the rotor system. The rotor system association model provided by the embodiment of the application is used for describing a mathematical model for associating the adjustment process parameters with the rotor system, and can associate the adjustment process parameters such as the adjustment load, the connection load and the like with the dynamic characteristics of the rotor system by taking the rigidity matrix of the rolling bearing as a bridge.
Specifically, the rotor system association model comprises an inner ring groove bottom diameter increase value calculation function, an outer ring groove bottom diameter decrease value calculation function, a bearing original initial contact angle calculation function, a corrected initial contact angle calculation function, a contact angle relation calculation function, an original equivalent stiffness coefficient information calculation function, a bearing stiffness information calculation function, a rotor system stiffness matrix calculation function and a rotor system dynamics mechanism function.
The basic parameter set information further comprises, for example, first point contact parameter information of an inner ring of the rolling bearing, second point contact parameter information of an outer ring of the rolling bearing, first curvature function information of the inner ring of the rolling bearing, second curvature function information of the outer ring of the rolling bearing, radial play information of the rolling bearing, first raceway curvature radius coefficient information of the inner ring of the rolling bearing, second raceway curvature radius coefficient information of the outer ring of the rolling bearing, rolling element diameter information of the rolling bearing, rolling element number information of the rolling bearing and load-adjusting information of the rotor system, wherein the first point contact parameter information is used for describing a point contact parameter of the inner ring of the rolling bearing, the second point contact parameter information is used for describing a point contact parameter of the outer ring of the rolling bearing, the first curvature function information is used for describing a curvature function of the inner ring of the rolling bearing, the second curvature function information is used for describing a curvature function of the outer ring of the rolling bearing, the radial play information is used for describing a radial play coefficient of the rolling bearing, the first raceway curvature radius coefficient information is used for describing a raceway curvature radius coefficient of the inner ring of the rolling bearing, the second raceway curvature radius coefficient information is used for describing a raceway curvature coefficient of the outer ring of the rolling bearing, the rolling element diameter information is used for describing a raceway curvature coefficient of the rolling element number of the rolling element, and the load-adjusting information of the rotor system is used for describing a number of rolling element number of rolling elements.
It should be noted that, each data in the basic parameter set information may be user-defined, manually input, and/or calculated by using the existing technology, so that details are omitted.
Specifically, the terminal device may construct a rotor system association model of the rotor system according to the basic parameter set information.
In some possible implementations, to facilitate analysis of the effect of bearing adjustment process parameters on rotor system dynamics, referring to fig. 4, step S200 includes, but is not limited to, the steps of:
in S210, inner ring groove bottom diameter increase value information and outer ring groove bottom diameter decrease value information are generated from the basic parameter set information, the inner ring groove bottom diameter increase value calculation function, and the outer ring groove bottom diameter decrease value calculation function.
Specifically, the inner ring groove bottom diameter increase value information is used to describe the groove bottom diameter increase value of the inner ring of the rolling bearing, and the outer ring groove bottom diameter decrease value information is used to describe the groove bottom diameter decrease value of the outer ring of the rolling bearing. The terminal equipment can input basic parameter set information into a preset inner ring groove bottom diameter increase value calculation function to generate inner ring groove bottom diameter increase value information; and inputting the basic parameter set information into a preset outer ring groove bottom diameter reduction value calculation function to generate outer ring groove bottom diameter reduction value information.
In some possible implementations, the inner ring groove bottom diameter increase value calculation function and the outer ring groove bottom diameter decrease value calculation function may both combine with elastic mechanics, where the inner ring groove bottom diameter increase value calculation function may be:
in delta F Information indicating an increase value of the diameter of the bottom of the inner ring groove; d represents inner diameter information of the rolling bearing; Δf 1 Representing first interference information; d, d F Representing first groove bottom diameter information;
in some possible implementations, the outer ring groove bottom diameter reduction value calculation function may be:
in delta E Information representing the diameter reduction value of the outer ring groove bottom; d represents the outer diameter information of the rolling bearing; Δf 2 Representing second interference information; d (D) E Representing second groove bottom diameter information; mu (mu) h Representing first poisson's ratio information; mu (mu) b Representing second poisson's ratio information; e (E) h Representing first elastic modulus information; e (E) b Representing second elastic modulus information.
In S220, the original initial contact angle information of the rolling bearing is generated according to the basic parameter set information and the bearing original initial contact angle calculation function.
Specifically, the raw initial contact angle information is used to describe the raw initial contact angle of the rolling bearing. The terminal equipment can input basic parameter set information into a preset bearing original initial contact angle calculation function to generate original initial contact angle information of the rolling bearing.
In some possible implementations, the bearing raw initial contact angle calculation function may be:
wherein alpha is o Representing raw initial contact angle information; u (u) r Representing radial play information; f (f) i Representing first channel radius of curvature coefficient information; f (f) e Representing second channel radius of curvature coefficient information; d (D) b Representing rolling element diameter information.
In S230, corrected initial contact angle information of the rolling bearing is generated from the inner ring groove bottom diameter increase value information, the outer ring groove bottom diameter decrease value information, the bearing original initial contact angle information, and the corrected initial contact angle calculation function.
Specifically, the corrected initial contact angle information is used to describe the corrected initial contact angle of the rolling bearing. The terminal equipment can input the inner ring groove bottom diameter increasing value information, the outer ring groove bottom diameter decreasing value information and the bearing original initial contact angle information into a preset correction initial contact angle calculation function to generate the correction initial contact angle information of the rolling bearing.
In some possible implementations, since the connection load is directly proportional to the interference, the interference may directly replace the connection load, and in order to fully consider the effect of the connection load on the gap, the above-mentioned modified initial contact angle calculation function may be:
Wherein, alpha' o Representing corrected initial contact angle information; u (u) r Representing radial play information; f (f) i Representing first channel radius of curvature coefficient information; f (f) e Representing second channel radius of curvature coefficient information; d (D) b Representing rolling element diameter information; delta E Information representing the diameter reduction value of the outer ring groove bottom; delta F And information indicating the increase value of the inner ring groove bottom diameter.
In S240, the original equivalent stiffness coefficient information of the rolling bearing is generated from the basic parameter set information and the original equivalent stiffness coefficient information calculation function.
Specifically, the original equivalent stiffness coefficient information is used to describe the original equivalent stiffness coefficient of the rolling bearing. The terminal equipment can input the basic parameter set information into a preset original equivalent stiffness coefficient information calculation function to generate original equivalent stiffness coefficient information of the rolling bearing.
In some possible implementations, the original equivalent stiffness coefficient information calculation function may be:
wherein K is n Representing original equivalent stiffness coefficient information; delta i Representing first point contact parameter information; delta o Representing the secondPoint contact parameter information; Σρ i Representing first curvature function information; Σρ o Representing second curvature function information.
In S250, real-time contact angle information is generated according to the corrected initial contact angle information, the original equivalent stiffness coefficient information and the contact angle relation calculation function.
In particular, real-time contact angle information is used to describe a contact angle of a rolling bearing with high realism and real-time. The terminal equipment can input the corrected initial contact angle information and the original equivalent stiffness coefficient information into a preset contact angle relation calculation function to generate real-time contact angle information.
In some possible implementations, in order to fully take into account the rotor system in modulating the load F α Under the action of the above, the contact angle relation calculation function can be:
wherein, alpha represents real-time contact angle information; f (F) α Representing load information of adjustment; k (K) n Representing original equivalent stiffness coefficient information; z represents rolling element number information; g represents a calculation factor, g=f e +f i -1;
α′ o Indicating corrected initial contact angle information.
In S260, bearing stiffness information of the rolling bearing is generated according to the real-time contact angle information, the original equivalent stiffness coefficient information and the bearing stiffness information calculation function.
In particular, the bearing stiffness information of the rolling bearing is used to describe the bearing stiffness of the rolling bearing. The terminal equipment can input the real-time contact angle information and the original equivalent stiffness coefficient information into a preset bearing stiffness information calculation function to generate bearing stiffness information of the rolling bearing.
In some possible implementations, the bearing stiffness information calculation function may be combined with the Hertz contact theory, and the bearing stiffness information calculation function may be:
K bz =1.5K n ·cos 2 α(δ α sinα)·0.5Z,
wherein K is bz Representing bearing stiffness information; k (K) n Representing original equivalent stiffness coefficient information; alpha represents real-time contact angle information; z represents rolling element number information; delta α And the axial displacement information after the rolling bearing is preloaded.
In S270, bearing stiffness matrix information of the rolling bearing is generated from the bearing stiffness information, the rotor system stiffness matrix calculation function, and the rotor system dynamics mechanism function.
In particular, the bearing stiffness matrix information is used to describe the bearing stiffness matrix of the rolling bearing. The terminal equipment can combine the bearing rigidity information input, the rotor system rigidity matrix calculation function and the rotor system dynamics mechanism function to generate the bearing rigidity matrix information of the rolling bearing.
In some possible implementations, in order to fully consider the situation that the rotor system is a rotating vibration system, there is a gyroscopic effect, and fully consider the situation that the supporting mode of the rolling bearing-rotor system has nonlinear influencing factors such as rigidity, damping and the like, the dynamic mechanism function of the rotor system may be combined with the characteristics of the rotor system and the finite element theory, and the dynamic mechanism function of the rotor system may be:
Wherein { F } represents external force vector information of the rotor system; [ M ]]Quality matrix information representing the rotor system; [ C]Damping matrix information representing the rotor system; [ G]Gyro matrix information representing the rotor system; [ K ]]Stiffness matrix information representing the rotor system; { q } represents generalized coordinate vector information of the rotor system;a first derivative representing generalized coordinate vector information of the rotor system; />A second derivative representing generalized coordinate vector information of the rotor system;
the rotor system stiffness matrix calculation function may be:
[K]=[K bz ]+[K S ]-ω 2 [M C ],
in [ K ] bz ]Bearing stiffness matrix information representing the rolling bearing; [ K ]]Stiffness matrix information representing the rotor system; [ K ] S ]Rotating shaft rigidity matrix information representing the rotating shaft; omega represents angular velocity information of the rotating shaft; [ M ] C ]Matrix information is added to the centrifugal force representing the rotation axis.
In some possible implementations, referring to fig. 5 and 6, the terminal device may build a finite element model of the rotor system based on Ansys finite element software, the finite element model of the rotor system including the spindle, the turntable, and the bearing stiffness characteristics, and may simulate the bearing stiffness characteristics with bearing units (e.g., COMBI 214), and may simulate the bearing stiffness characteristics with SOLID units in the Ansys Workbench (e.g., SOLID 186), while the rotor system may be divided into hexahedral grids with a Multizone grid division method, while material properties of the rotor system, such as modulus of elasticity, poisson ratio, and density of the rotor system, may be input by a user.
In some possible implementations, the basic parameters of the rotor system finite element model may be the data in table 2 below.
Table 2 basic parameters of finite element model of rotor system
Span of rotating shaft | 1000mm |
Diameter of spindle | 90mm |
Outer diameter of turntable | 360mm |
Inner diameter of turntable | 90mm |
Thickness of turntable | 80mm |
Density of material | 7.850×10 3 kg/m 3 |
Modulus of elasticity | 2×10 11 Pa |
Poisson's ratio | 0.3 |
Network element | Hexahedron |
In some possible implementations, the connection relationships and imposed boundary conditions of the components in the rotor system finite element model may be based on the actual condition of the rotor system; for example, the terminal device may simulate the connection relationship between the turntable and the rotating shaft by using a bound contact type, and then apply a rotation speed to the rotating shaft and set a constraint type of Displacement at the bearing to limit the Z-direction (i.e. axial) movement, and not limit the X, Y-direction (i.e. radial) movement, so as to enable the X, Y-direction movement to be free. Therefore, the restriction of adding Displacement at the bearing is realized by combining the characteristic of simulating the rigidity of the bearing, the finite element model of the rotor system is more similar to the actual working condition, the restriction and supporting effect of the rolling bearing on the rotor system are simulated, and the method is beneficial to the implantation of the loading and unloading load, the connecting load and other loading and unloading parameters.
In S300, analysis dataset information for the rotor system is generated from the test dataset and the rotor system correlation model.
Specifically, the terminal device may input a test data set into a preset rotor system association model, and generate analysis data set information of the rotor system, where the analysis data set information is used to describe analysis data related to an effect of the bearing adjustment process parameter on the dynamic characteristics of the rotor system.
In some possible implementations, to further facilitate analysis of the effect of bearing adjustment process parameters on rotor system dynamics, referring to fig. 7, step S300 includes, but is not limited to, the steps of:
in S310, a test dataset is input into the rotor system association model based on the Lanczos algorithm, and a simulation result set is generated.
Without loss of generality, the test data set comprises a plurality of test assembly load information and test interference information, and the simulation result set comprises simulation bearing stiffness matrix information output by a plurality of rotor system association models.
Specifically, the terminal device may input the test data set into the rotor system association model based on the Lanczos algorithm, and generate a simulation result set, where the simulation result set is used to describe high-validity result data output by the rotor system association model.
Illustratively, referring to table 3 below, the load-on-load in table 3 is test load-on-load information; when the interference between the bearing inner ring and the rotating shaft is 4 mm and the rolling bearing rotates 4000 revolutions per minute, the terminal equipment can determine the bearing rigidity information according to the test assembly load information.
TABLE 3 load F adjustment α With bearing stiffness K bz Relation table of (2)
For example, referring to Table 4 below, the magnitude of the interference in Table 4 is the test interferenceQuantity information; when the load F is adjusted α When the rolling bearing rotates at 15000 rpm for 30N, the terminal equipment can determine the bearing rigidity information according to a plurality of pieces of interference magnitude information.
TABLE 4 interference X and bearing stiffness K bz Relation table of (2)
In S320, a critical rotational speed dataset of the rotor system is determined from the simulation result set.
Specifically, the terminal device may determine a critical rotational speed dataset of the rotor system according to the simulation result set.
In some possible implementation manners, in order to facilitate intuitively determining the influence of the bearing adjustment process parameter on the dynamic characteristics of the rotor system, referring to fig. 8, the terminal device may open the coiol is gyro effect based on the Workbench, and then select the Lanczos algorithm to perform multi-load step modal analysis on the rotor system at different rotation speeds, so as to directly obtain a campbel diagram and a vibration pattern diagram. Referring to fig. 9 and fig. 10, after the terminal device generates a campbel l diagram, the terminal device may determine the critical rotation speed through the campbel l diagram, and analyze to obtain the result of the influence of the load of adjustment and the load of connection on the critical rotation speed of the rotor system.
In some possible implementations, referring to fig. 11 and 12, the terminal device may perform harmonic response analysis on the rotor system based on the complete method of Workbench and opening the coiol is gyro effect, and determine an unbalanced response of the rotor system, so as to obtain an influence result of the load adjustment and the load connection on the unbalanced response of the rotor system.
In S330, analysis dataset information of the rotor system is generated from the test dataset and the critical rotational speed dataset.
In particular, the terminal device may generate analysis dataset information of the rotor system from the test dataset and the critical rotation speed dataset
In some possible implementations, to further facilitate analysis of the effect of the bearing adjustment process parameters on the rotor system dynamics, referring to fig. 13, after step S320, the method further includes, but is not limited to, the steps of:
in S321, a dynamic stability feature value of the rotor system is determined according to the test data set and a preset feature value calculation function.
In particular, the test data set may also include amplitude information of the rotor system; the terminal device can input the test data set into a preset characteristic value calculation function to determine the dynamic stability characteristic value of the rotor system.
In some possible implementations, in order to fully consider the damping effect when analyzing the influence of the tuning parameters such as the tuning load and the connection load on the dynamic stability of the rotor system, and extract the characteristics of the dynamic mechanism function of the rotor system by combining the speed method, the characteristic value calculation function may be:
Y=λ+iω,
wherein lambda represents a dynamic stability characteristic value, and a smaller dynamic stability characteristic value represents a more stable rotor system; y represents amplitude information of the rotor system; i represents an imaginary unit; ω is modal frequency information of the rotor system.
In S322, state information of the rotor system is determined according to the dynamic stability feature value and a preset steady state determination function.
Specifically, the terminal device may input the state stability characteristic value to a preset stable state determining function, to determine state information of the rotor system.
In some possible implementations, the steady state determination function may be:
in the formula, ment RotorSystem Status information representing the rotor system; when lambda is<When 0, the state information indicates a steady state; when λ=0, the state information indicates a gradual steady state, indicating that the system damping is 0; when lambda > 0, the shape The state information represents a destabilized state;
accordingly, the step S330 includes, but is not limited to, the following steps:
in S331, analysis dataset information of the rotor system is generated from the test dataset, the critical rotation speed dataset, the dynamic stability feature value and the state information.
In particular, the terminal device may generate analysis dataset information for the rotor system based on the test dataset, the critical rotational speed dataset, the dynamic stability feature value, and the status information.
In some possible implementation manners, referring to fig. 14 and 15, the terminal device may open a pool is gyroscopic effect based on Workbench, then perform multi-load step modal analysis on the rotor system at different rotation speeds by combining with a speed method, directly obtain a value of λ, and then analyze dynamic stability by the value of λ to obtain an effect result of adjusting load and connecting load on the dynamic stability of the rotor system.
The implementation principle of the quantitative analysis method for the influence of the adjustment process parameters on the rotor system in the embodiment of the application is as follows: the terminal equipment can firstly acquire basic parameter set information of the rolling bearing, then construct a rotor system association model of the rotor system according to the basic parameter set information, then generate analysis data set information of the rotor system according to the test data set and the rotor system association model, fully consider the influence of the bearing adjustment process parameters on the dynamic characteristics of the rotor system, more accurately predict the dynamic characteristics of the rotor system, and be favorable for analyzing the influence of the bearing adjustment process parameters on the dynamic characteristics of the rotor system from multiple dimensions.
It should be noted that, the sequence number of each step in the above embodiment does not mean the sequence of execution sequence, and the execution sequence of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiment of the present application.
The embodiment of the present application also provides a quantitative analysis system for the influence of an adjustment process parameter on a rotor system, which is applied to the rotor system, wherein the rotor system comprises a bearing seat, a rolling bearing and a rotating shaft, the rolling bearing is mounted on the bearing seat, the rotating shaft is mounted on the rolling bearing, for convenience of description, only the part relevant to the present application is shown, as shown in fig. 16, the system 160 comprises:
the basic parameter set information acquisition module 161: the method comprises the steps of acquiring basic parameter set information of a rolling bearing, wherein the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first groove bottom diameter information of the inner ring of the rolling bearing, second groove bottom diameter information of the outer ring of the rolling bearing, first poisson ratio information of the bearing seat, second poisson ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing;
Rotor system correlation model building module 162: the method comprises the steps of constructing a rotor system association model of a rotor system according to basic parameter set information, wherein the rotor system association model is used for describing a mathematical model for associating an adjustment process parameter with the rotor system;
the analysis data set information generation module 163: for generating analysis dataset information of a rotor system from a test dataset and a rotor system correlation model
Optionally, the analysis data set information generating module 163 further includes:
simulation result set generation submodule: the method comprises the steps of inputting a test data set into a rotor system association model based on a Lanczos algorithm to generate a simulation result set, wherein the test data set comprises a plurality of pieces of test assembly load information and test interference magnitude information, and the simulation result set comprises simulation bearing rigidity matrix information output by the rotor system association model;
a critical rotation speed dataset determination submodule: the method comprises the steps of determining a critical rotation speed data set of a rotor system according to a simulation result set;
an analysis dataset information generation sub-module: and the analysis data set information is used for generating the analysis data set information of the rotor system according to the test data set and the critical rotation speed data set.
It should be noted that, because the content of information interaction and execution process between the modules is based on the same concept as the method embodiment of the present application, specific functions and technical effects thereof may be referred to in the method embodiment section, and details are not repeated herein.
The embodiment of the present application further provides a terminal device, as shown in fig. 17, where the terminal device 170 of this embodiment includes: a processor 170, a memory 172 and a computer program 173 stored in the memory 172 and executable on the processor 170. The processor 170, when executing the computer program 173, implements the steps in the above-described flow processing method embodiment, such as steps S100 to S300 shown in fig. 1; alternatively, the processor 170, when executing the computer program 173, performs the functions of the respective modules in the above-described apparatus, for example, the functions of the modules 161 to 163 shown in fig. 16.
The terminal device 170 may be a desktop computer, a notebook computer, a palm top computer, a cloud server, etc., and the terminal device 170 includes, but is not limited to, a processor 170 and a memory 172. It will be appreciated by those skilled in the art that fig. 17 is merely an example of the terminal device 170 and is not meant to be limiting as the terminal device 170, may include more or fewer components than shown, or may combine certain components, or different components, e.g., the terminal device 170 may also include an input-output device, a network access device, a bus, etc.
The processor 170 may be a central processing unit (Central Processing Unit, CPU), other general purpose processor, digital signal processor (Digital Signal Processor, DSP), application specific integrated circuit (Application Specific Integrated Circuit, ASIC), field-programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, etc.; a general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 172 may be an internal storage unit of the terminal device 170, such as a hard disk or a memory of the terminal device 170, or the memory 172 may be an external storage device of the terminal device 170, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card) or the like, which are provided on the terminal device 170; further, the memory 172 may also include both an internal storage unit and an external storage device of the terminal device 170, the memory 172 may also store the computer program 173 and other programs and data required by the terminal device 170, and the memory 172 may also be used to temporarily store data that has been output or is to be output.
An embodiment of the present application also provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the steps of the various method embodiments described above. Wherein the computer program comprises computer program code, the computer program code can be in the form of source code, object code, executable file or some intermediate form, etc.; the computer readable medium may include: any entity or device capable of carrying computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth.
The foregoing are all preferred embodiments of the present application, and are not intended to limit the scope of the present application in any way, therefore: all equivalent changes of the method, principle and structure of the present application should be covered in the protection scope of the present application.
Claims (9)
1. The method for quantitatively analyzing the influence of the installation and adjustment process parameters on the rotor system is applied to the rotor system, and the rotor system comprises a bearing seat, a rolling bearing and a rotating shaft, wherein the rolling bearing is arranged on the bearing seat, and the rotating shaft is arranged on the rolling bearing, and is characterized in that the method comprises the following steps:
acquiring basic parameter set information of the rolling bearing, wherein the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first groove bottom diameter information of the inner ring of the rolling bearing, second groove bottom diameter information of the outer ring of the rolling bearing, first poisson ratio information of the bearing seat, second poisson ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing;
Constructing a rotor system association model of the rotor system according to the basic parameter set information, wherein the rotor system association model is used for describing a mathematical model for associating the adjustment process parameters with the rotor system;
and generating analysis data set information of the rotor system according to the test data set and the rotor system association model.
2. The method of claim 1, wherein the base parameter set information further comprises first point contact parameter information of an inner ring of the rolling bearing, second point contact parameter information of an outer ring of the rolling bearing, first curvature function information of the inner ring of the rolling bearing, second curvature function information of the outer ring of the rolling bearing, radial play information of the rolling bearing, first channel radius of curvature coefficient information of the inner ring of the rolling bearing, second channel radius of curvature coefficient information of the outer ring of the rolling bearing, rolling body diameter information of the rolling bearing, rolling body number information of the rolling bearing, and load adjustment information of the rotor system; the rotor system association model comprises an inner ring groove bottom diameter increase value calculation function, an outer ring groove bottom diameter decrease value calculation function, a bearing original initial contact angle calculation function, a corrected initial contact angle calculation function, a contact angle relation calculation function, an original equivalent stiffness coefficient information calculation function, a bearing stiffness information calculation function, a rotor system stiffness matrix calculation function and a rotor system dynamics mechanism function; the constructing a rotor system association model of the rotor system according to the basic parameter set information comprises the following steps:
Generating inner ring groove bottom diameter increase value information and outer ring groove bottom diameter decrease value information according to the basic parameter set information, the inner ring groove bottom diameter increase value calculation function and the outer ring groove bottom diameter decrease value calculation function;
generating original initial contact angle information of the rolling bearing according to the basic parameter set information and the original initial contact angle calculation function of the bearing;
generating corrected initial contact angle information of the rolling bearing according to the inner ring groove bottom diameter increase value information, the outer ring groove bottom diameter decrease value information, the original initial contact angle information of the bearing and the corrected initial contact angle calculation function;
according to the basic parameter set information and the original equivalent stiffness coefficient information calculation function, original equivalent stiffness coefficient information of the rolling bearing is generated;
generating real-time contact angle information according to the corrected initial contact angle information, the original equivalent stiffness coefficient information and the contact angle relation calculation function;
generating bearing stiffness information of the rolling bearing according to the real-time contact angle information, the original equivalent stiffness coefficient information and the bearing stiffness information calculation function;
And generating bearing rigidity matrix information of the rolling bearing according to the bearing rigidity information, the rotor system rigidity matrix calculation function and the rotor system dynamics mechanism function.
3. The method of claim 2, wherein the inner race groove bottom diameter increase value calculation function is:
in delta F Increasing value information for the diameter of the bottom of the inner ring groove; d is the inner diameter information of the rolling bearing; Δf 1 Is the first interference information; d, d F The first groove bottom diameter information;
the outer ring groove bottom diameter reduction value calculation function is as follows:
in delta E The diameter of the outer ring groove bottom is reduced to obtain information; d is the outer diameter information of the rolling bearing; Δf 2 Is the second interference information; d (D) E The second groove bottom diameter information; mu (mu) h For the first poisson ratio information; mu (mu) b For the second poisson ratio information; e (E) h Information about the first elastic modulus; e (E) b Information about the second elastic modulus;
the original initial contact angle calculation function of the bearing is as follows:
wherein alpha is o Is the original initial contact angle information; u (u) r Is the radial play information; f (f) i Information on a radius of curvature coefficient of the first channel; f (f) e Radius of curvature coefficient information for the second channel; d (D) b Diameter information of the rolling bodies;
the modified initial contact angle calculation function is:
wherein, alpha' o For the corrected initial contact angle information; u (u) r Is the radial play information; f (f) i Information on a radius of curvature coefficient of the first channel; f (f) e Radius of curvature coefficient information for the second channel; d (D) b Diameter information of the rolling bodies; delta E The diameter of the outer ring groove bottom is reduced to obtain information; delta F Increasing value information for the diameter of the bottom of the inner ring groove;
the original equivalent stiffness coefficient information calculating function is as follows:
wherein K is n The original equivalent stiffness coefficient information is obtained; delta i The first point contact parameter information; delta o The second point contact parameter information is the second point contact parameter information; Σρ i -providing the first curvature function information; Σρ o -providing said second curvature function information;
the contact angle relation calculation function is as follows:
wherein alpha is the real-time contact angle information; f (F) α Load information is adjusted for the load; k (K) n The original equivalent stiffness coefficient information is obtained; z is the number information of the rolling bodies; g is a calculation factor, g=f e +f i -1;α′ o For the corrected initial contact angle information;
the bearing rigidity information calculating function is as follows:
bz 2
K=1.5K n ·cosα(δ α sinα)·0.5Z,
wherein K is bz -providing said bearing stiffness information; k (K) n The original equivalent stiffness coefficient information is obtained; alpha is the real-time contact angle information; z is the number information of the rolling bodies; delta α Axial displacement information after the rolling bearing is pre-tightened;
the dynamic mechanism function of the rotor system is as follows:
wherein { F } is the external force vector information of the rotor system; [ M ]]Is saidMass matrix information of the rotor system; [ C]Damping matrix information for the rotor system; [ G]Gyro matrix information for the rotor system; [ K ]]Stiffness matrix information for the rotor system; { q } is the generalized coordinate vector information of the rotor system;a first derivative of generalized coordinate vector information for the rotor system; />Second derivative of generalized coordinate vector information for the rotor system;
the rigidity matrix calculation function of the rotor system is as follows:
4. a method according to claim 3, wherein said generating analytical dataset information for the rotor system from a test dataset and the rotor system correlation model comprises:
inputting a test data set into the rotor system association model based on Lanczos algorithm to generate a simulation result set, wherein the test data set comprises a plurality of pieces of test assembly load information and test interference magnitude information, and the simulation result set comprises a plurality of pieces of simulation bearing rigidity matrix information output by the rotor system association model;
Determining a critical rotation speed data set of the rotor system according to the simulation result set;
and generating analysis data set information of the rotor system according to the test data set and the critical rotation speed data set.
5. The method of claim 4, wherein after said determining a critical speed dataset for the rotor system from the simulation result set, the method further comprises:
determining a dynamic stability characteristic value of the rotor system according to the test data set and a preset characteristic value calculation function;
determining state information of the rotor system according to the dynamic stability characteristic value and a preset stable state determining function;
wherein, the eigenvalue calculation function is:
Y=λ+iω,
wherein lambda is the dynamic stability characteristic value; y is the amplitude information of the rotor system; i is an imaginary unit; omega is modal frequency information of the rotor system;
the steady state determination function is:
in the formula, ment RotorSystem Status information for the rotor system; when lambda is<0, the state information is in a stable state; when λ=0, the state information is a progressive steady state; when lambda > 0, the state information is in an unstable state;
Accordingly, the generating analysis dataset information of the rotor system according to the test dataset and the critical rotation speed dataset comprises:
and generating analysis data set information of the rotor system according to the test data set, the critical rotation speed data set, the dynamic stability characteristic value and the state information.
6. The utility model provides a quantitative analysis system of dress adjustment technological parameter to rotor system influence is applied to the rotor system, the rotor system includes bearing frame, antifriction bearing and pivot, antifriction bearing install in the bearing frame, the pivot install in antifriction bearing, its characterized in that, the system includes:
basic parameter set information acquisition module: the method comprises the steps of obtaining basic parameter set information of the rolling bearing, wherein the basic parameter set information comprises inner diameter information of the rolling bearing, outer diameter information of the rolling bearing, first interference information of an inner ring of the rolling bearing and a rotating shaft in the diameter direction, second interference information of an outer ring of the rolling bearing and a bearing seat hole of a bearing seat in the diameter direction, first groove bottom diameter information of the inner ring of the rolling bearing, second groove bottom diameter information of the outer ring of the rolling bearing, first poisson ratio information of the bearing seat, second poisson ratio information of the rolling bearing, first elastic modulus information of the bearing seat and second elastic modulus information of the rolling bearing;
The rotor system association model building module: the rotor system association model is used for constructing a rotor system association model of the rotor system according to the basic parameter set information, wherein the rotor system association model is used for describing a mathematical model for associating the adjustment process parameters with the rotor system;
an analysis data set information generation module: and the analysis data set information is used for generating the rotor system according to the test data set and the rotor system association model.
7. The system of claim 6, wherein the analysis dataset information generation module comprises:
simulation result set generation submodule: the simulation method comprises the steps of inputting a test data set into a rotor system association model based on a Lanczos algorithm to generate a simulation result set, wherein the test data set comprises a plurality of pieces of test assembly load information and test interference magnitude information, and the simulation result set comprises a plurality of pieces of simulation bearing rigidity matrix information output by the rotor system association model;
a critical rotation speed dataset determination submodule: determining a critical rotation speed data set of the rotor system according to the simulation result set;
an analysis dataset information generation sub-module: and the analysis data set information of the rotor system is generated according to the test data set and the critical rotation speed data set.
8. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any of claims 1 to 5 when the computer program is executed.
9. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311219350.2A CN117454685A (en) | 2023-09-20 | 2023-09-20 | Quantitative analysis method and system for influence of adjustment process parameters on rotor system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311219350.2A CN117454685A (en) | 2023-09-20 | 2023-09-20 | Quantitative analysis method and system for influence of adjustment process parameters on rotor system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117454685A true CN117454685A (en) | 2024-01-26 |
Family
ID=89578826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311219350.2A Pending CN117454685A (en) | 2023-09-20 | 2023-09-20 | Quantitative analysis method and system for influence of adjustment process parameters on rotor system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117454685A (en) |
-
2023
- 2023-09-20 CN CN202311219350.2A patent/CN117454685A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111141518B (en) | Model-based asymmetric rotor bearing system unbalance identification method | |
CN109556765A (en) | A kind of blade non-contacting vibration strain measurements conversion method | |
Zhao et al. | A novel transient balancing technology of the rotor system based on multi modal analysis and feature points selection | |
Jia et al. | Dynamic analysis of Jeffcott rotor under uncertainty based on Chebyshev convex method | |
Liu et al. | Stochastic dynamic analysis of the rotor–bearing system considering the randomness of the radial clearance | |
Zhang et al. | A comprehensive deep learning library benchmark and optimal library selection | |
Cheng et al. | Research on the effect of structural and material parameters on vibrations based on quasi-static model of bearings | |
CN117454685A (en) | Quantitative analysis method and system for influence of adjustment process parameters on rotor system | |
CN116839560B (en) | Hemispherical resonator gyroscope and hemispherical resonator quality leveling method, equipment and medium thereof | |
CN113268908A (en) | Response solving method and device of rotor system | |
CN114936483B (en) | Finite element equivalent modeling method for hydrostatic simulation analysis of large-scale bearing | |
CN112556931B (en) | Particle swarm algorithm-based modal dynamic balance method for high-speed bearing rotor system | |
CN112231950B (en) | Noise simulation prediction method, device and equipment of single-stage centrifugal blower | |
CN115310327A (en) | Simulation analysis method and device for multi-working-condition combination | |
Kang et al. | Field simultaneous estimation of residual unbalance and bearing dynamic coefficients for double-disk rotor-bearing system using dual augmented Kalman filter | |
JP6238796B2 (en) | Calculation method, calculation apparatus, and calculation system | |
Liu et al. | An efficient dynamic modeling technique for a central tie rod rotor | |
CN110309615B (en) | Method for predicting natural frequency of rotating blade | |
CN113378322B (en) | Method, device and equipment for optimizing structural parameters of rotating piece and storage medium | |
CN115310235B (en) | Simulation analysis coupling method for acquiring key parameters of machine tool spindle | |
Jauhari et al. | Evaluation of dynamic behavior a machine tool spindle system through modal and unbalance response analysis | |
CN114544124B (en) | Triaxial vibration stress spectrum processing method, device, equipment and storage medium | |
Bernhauser et al. | Efficient Transient Run-up Simulations for the Investigation of Acoustically Relevant Vibrations of Turbocharger Structures with Floating-Ring Bearings | |
da Silva Tuckmantel et al. | An analysis on the supporting structure representative model in rotating systems | |
Xia et al. | Based on ANSYS crown tapered rolling bearings dynamics simulation analysis |
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 |
Country or region after: China Address after: 528000 No. 18, Jiangwan Road, Chancheng District, Guangdong, Foshan Applicant after: Foshan University Address before: 528000 No. 18, Jiangwan Road, Chancheng District, Guangdong, Foshan Applicant before: FOSHAN University Country or region before: China |
|
CB02 | Change of applicant information |