Flywheel bearing assembly disturbance characteristic analysis method

Abstract

A flywheel bearing assembly disturbance characteristic analysis method relates to the technical field of bearing disturbance dynamics simulation analysis, and comprises the following steps: S1, obtaining initial conditions of a flywheel bearing assembly, and establishing a flywheel bearing assembly coordinate system; s2, establishing an elastic damping connection model between a flywheel frame coordinate system and a space fixed coordinate system; s3, establishing a disturbance calculation model of the flywheel bearing assembly; s4, analyzing a relative position and a motion relation between parts in the flywheel bearing assembly; s5, analyzing the interaction force between the parts in the flywheel bearing assembly, and establishing a nonlinear dynamic differential equation set of each part; s6, solving a nonlinear dynamic differential equation set according to a bearing dynamic theory; and S7, judging whether the dynamic equation set meets a convergence error or not. The method can accurately analyze the influence of different bearing rotating speeds, loads and size parameters on the overall disturbance characteristics of the flywheel assembly, and solves the technical problem that the disturbance characteristics of the flywheel bearing assembly cannot be accurately analyzed by the existing method.

Description

Technical Field

[0001] The present invention relates to the technical field of bearing disturbance dynamics simulation analysis, in particular to a method for analyzing disturbance characteristics of a flywheel bearing assembly. Background Art

[0002] With the recent development of aerospace, new-generation observation instruments have increasingly demanded higher pointing accuracy, which in turn places increasing demands on satellite stability. Microvibration disturbances generated by high-speed rotating equipment, such as flywheels, are a significant obstacle to improving satellite stability. Among the many factors that cause spacecraft microvibration, flywheel-generated disturbances are the largest source of high-frequency microvibrations throughout the satellite.

[0003] As the core component of the flywheel assembly that supports the rotation of the wheel body, the bearing has a complex relative motion state and mechanical environment among its internal parts. During the operation of the flywheel, the dynamic characteristics of the bearing are closely related to the overall disturbance output of the assembly. The design, processing and installation accuracy of the bearing directly affect the overall disturbance output and operation stability of the flywheel assembly.

[0004] Due to the small magnitude of micro-vibrations, experimental analysis of the effects of bearing preload, structural parameters, and surface defects on the overall disturbance of a flywheel assembly is difficult to achieve. Currently, there are few analytical methods for the disturbance characteristics of flywheel-bearing assemblies, and most of them focus on the dynamic characteristics of a single bearing, ignoring the impact of other components in the flywheel assembly on the overall disturbance characteristics. In particular, the inherent connection between the dynamic characteristics of the bearing and the overall disturbance characteristics of the flywheel assembly is not comprehensively considered, resulting in an inability to accurately analyze the disturbance characteristics of the flywheel-bearing assembly. Summary of the Invention

[0005] The purpose of the present invention is to provide a method for analyzing the disturbance characteristics of a flywheel bearing assembly. Based on the rolling bearing dynamics theory, the non-rigid installation of the flywheel bearing assembly is comprehensively considered, and a disturbance dynamics analysis model of the flywheel bearing assembly is established. The method can accurately analyze the influence of different bearing speeds, loads and dimensional parameters on the overall disturbance characteristics of the flywheel assembly, and can solve the technical problem that the existing analysis methods cannot accurately analyze the disturbance characteristics of the flywheel bearing assembly.

[0006] In order to achieve the above-mentioned purpose, the present invention adopts the following technical solutions.

[0007] A method for analyzing disturbance characteristics of a flywheel bearing assembly comprises the following steps: S1. Obtaining the initial conditions of the flywheel bearing assembly and establishing a flywheel bearing assembly coordinate system, including a space fixed coordinate system and a flywheel frame coordinate system; S2. Establishing an elastic damping connection model between the flywheel frame coordinate system and the space fixed coordinate system; S3. Establish a disturbance calculation model for the flywheel bearing assembly; S4. Analyze the relative position and motion relationship between the parts in the flywheel bearing assembly; S5. Analyze the interaction forces between the parts in the flywheel bearing assembly and establish a set of nonlinear dynamic differential equations for each part; S6. Based on the bearing dynamics theory, the Gstiff variable step size integration algorithm is used to solve the nonlinear dynamic differential equations; S7. Determine whether the dynamic equations satisfy the convergence error.

[0008] Furthermore, in step S1, the initial conditions of the flywheel bearing assembly obtained include flywheel frame structural parameters, bearing center distance, bearing structural parameters, material properties, lubrication parameters and operating conditions.

[0009] Furthermore, in step S1, the established flywheel bearing assembly coordinate system also includes a bearing outer ring coordinate system, a bearing retainer center coordinate system, a retainer pocket center coordinate system, and a steel ball center coordinate system.

[0010] Furthermore, in step S2, when establishing the elastic damping connection model between the flywheel frame coordinate system and the space fixed coordinate system, the stiffness coefficient and damping coefficient of the connection between the flywheel frame coordinate system and the space fixed coordinate system are set.

[0011] Furthermore, in step S3, the established flywheel bearing assembly disturbance calculation model includes a calculation model of the reaction flywheel bearing assembly disturbance force and disturbance torque.

[0012] Furthermore, in step S4, the motion relationship between the steel ball and the cage relative to the outer ring, the motion relationship between the steel ball and the cage relative to the inner ring, and the motion relationship between the steel ball and the cage pocket center are determined based on the load, speed and bearing size parameters.

[0013] Furthermore, in step S5, based on the relative positions and motion parameters of the bearing parts in S4, the interaction forces between the rolling elements and raceways, the rolling elements and the cage, and the cage and the rings in the bearing are solved, and a group of nonlinear dynamic differential equations is established.

[0014] Furthermore, in step S7, the convergence error is preset to 1×10 -3 When the system residual meets the preset convergence threshold, the spatial position and motion state parameters of the bearing assembly at the current moment are recorded. When the system residual exceeds the convergence threshold, the current calculation result is used as the initial condition to update to the next iteration step, and the differential equation is solved again.

[0015] After adopting the above technical solution, the present invention has the following beneficial effects: 1. The analysis method of the present invention analyzes the flywheel assembly as a whole, models the components such as the bearing and the flywheel frame as a multi-body system, and associates the motion of each component through a coordinate system, which can more accurately analyze the disturbance characteristics of the flywheel assembly; 2. The present invention simulates the flexible connection characteristics of the elastic support by introducing stiffness and damping parameters between the flywheel frame coordinate system and the fixed coordinate system. This can realistically analyze the working environment of the component and the influence of the installation stiffness on the disturbance transmission. 3. The present invention establishes a fully parametric simulation model based on the rolling bearing dynamics theory. By replacing some experiments with numerical calculations, it can simulate the effects of different working conditions, bearing parameters, and installation conditions on disturbances, and predict disturbance characteristics without the need for a physical prototype. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 It is a schematic flow chart of the analysis method of the present invention.

[0017] Figure 2 Schematic diagram of the overall structure of the flywheel bearing assembly in the present invention (cross-section view).

[0018] Figure 3 It is a reference diagram of the bearing coordinate system in the present invention.

[0019] Figure 4 It is a schematic diagram of the non-rigid connection between the frame base and the support base.

[0020] Figure 5 It is a schematic diagram of the disturbance output analysis model of the flywheel bearing assembly.

[0021] Figure 6 It is a schematic diagram of the steel ball motion analysis.

[0022] Figure 7 It is a schematic diagram of the non-contact state in the motion relationship between the steel ball and the center of the cage pocket.

[0023] Figure 8 It is a schematic diagram of the state of the steel ball driving the cage in the motion relationship between the steel ball and the center of the cage pocket.

[0024] Figure 9 It is a schematic diagram of the state in which the retainer pushes the steel ball in the motion relationship between the steel ball and the center of the retainer pocket.

[0025] Figure 10 It is a schematic diagram of the force analysis of the steel ball.

[0026] Figure 11 It is a schematic diagram of the force analysis of the steel ball at a zero-one angle.

[0027] Figure 12 Is the cage in and Schematic diagram of plane force analysis.

[0028] Figure 13 Is the cage in and Schematic diagram of plane force analysis.

[0029] Figure 14 Is the cage in and Schematic diagram of plane force analysis.

[0030] Figure 15 It is a schematic diagram of the flywheel bearing dynamic analysis results.

[0031] Figure 16 It is a schematic diagram of the disturbance characteristics analysis results of the flywheel bearing assembly.

[0032] Description of the drawings: 1. Support seat, 2. Frame, 21. Frame base, 22. Frame shaft, 3. Wheel body, 4. Bearing, 41. Outer ring, 42. Inner ring, 43. Cage, 44. Steel ball. DETAILED DESCRIPTION

[0033] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the features and performance of a method for analyzing disturbance characteristics of a flywheel bearing assembly in the present invention are further described in detail below with reference to the accompanying drawings and embodiments.

[0034] Please see the attached Figures 1 to 16 , a method for analyzing the disturbance characteristics of a flywheel bearing assembly includes the following steps.

[0035] S1. Obtain the initial conditions of the flywheel bearing assembly, including the flywheel frame structural parameters, bearing center distance, bearing structural parameters, material properties, lubrication parameters, and operating conditions. Establish the flywheel bearing assembly coordinate system, including the fixed spatial coordinate system, the flywheel frame coordinate system, the bearing outer ring coordinate system, the bearing cage center coordinate system, the cage pocket center coordinate system, and the steel ball center coordinate system.

[0036] S2. Establish an elastic damping connection model between the flywheel frame coordinate system and the space fixed coordinate system, and set the stiffness coefficient and damping coefficient of the connection between the flywheel frame coordinate system and the space fixed coordinate system.

[0037] S3. Establish a disturbance calculation model for the flywheel bearing assembly, including a calculation model for the disturbance force and disturbance torque of the reaction flywheel bearing assembly.

[0038] S4. Analyze the relative positions and motion relationships between the parts in the flywheel bearing assembly. Based on the load, speed, and bearing size parameters, determine the motion relationship between the steel ball and retainer relative to the outer ring, the motion relationship between the steel ball and retainer relative to the inner ring, and the motion relationship between the steel ball and the retainer pocket center.

[0039] S5. Analyze the interaction forces between the parts in the flywheel bearing assembly. Based on the relative positions and motion parameters of the parts in the bearing in S4, solve the interaction forces between the rolling elements and raceways, the rolling elements and the cage, and the cage and the rings in the bearing, and establish a set of nonlinear dynamic differential equations.

[0040] S6. Based on the bearing dynamics theory, the Gstiff variable step size integration algorithm is used to solve the nonlinear dynamic differential equations.

[0041] S7, determine whether the dynamic equations meet the convergence error, the preset convergence error is When the system residual meets the preset convergence threshold, the spatial position and motion state parameters of the bearing assembly at the current moment are recorded. When the system residual exceeds the convergence threshold, the current calculation result is used as the initial condition to update to the next iteration step, and the differential equation is solved again.

[0042] In a specific implementation, a method for analyzing disturbance characteristics of a flywheel bearing assembly includes the following steps.

[0043] S1. Obtain the initial conditions of the flywheel bearing assembly and establish the flywheel bearing assembly coordinate system.

[0044] The initial conditions obtained in S1 include the flywheel frame structural parameters, bearing center distance, bearing structural parameters, material properties, lubrication parameters, load, and speed. The established coordinate systems include the fixed spatial coordinate system, the flywheel frame coordinate system, the bearing outer ring coordinate system, the bearing cage center coordinate system, the cage pocket center coordinate system, and the steel ball center coordinate system.

[0045] Specifically, if Figure 2 The figure shows the overall structure of the cannon bearing assembly, which includes a support base 1, a frame 2, and a wheel 3 rotatably mounted on the frame. The frame 2 includes a frame base 21 and a frame shaft 22. The frame base 21 is non-rigidly connected to the support base 1, and the wheel 3 is rotatably mounted on the frame shaft 22 via a pair of bearings 4.

[0046] In this embodiment, the flywheel bearing assembly used as an example adopts an angular contact ball bearing of model B7005, and its main parameters are shown in Table 1.

[0047] The material properties include the elastic modulus, Poisson's ratio, and density of each component in the bearing, as shown in Table 2.

[0048] Lubrication parameters include lubricant density, dynamic viscosity, viscosity-pressure coefficient, and thermal conductivity, etc. Working conditions include radial force, axial force, overturning moment, and rotational speed.

[0049] The coordinate system of each part of the flywheel bearing assembly is established, such as Figure 3 As shown, where: is a space-fixed coordinate system, which is a space-fixed coordinate system; is the flywheel frame coordinate system, The axis is along the axis of the shaft system, and the motion state of each part in the bearing can be analyzed through the reference coordinate system; The outer ring coordinate system, also known as the outer ring overall coordinate system, is the overall coordinate system of the outer rings of the bearings on both sides, the outer sleeve and the wheel body. The origin of the outer ring coordinate system coincides with the center of the shaft system. The axis of rotation of the shaft and the outer ring coincide; is the center coordinate system of the cage, where the subscript for Indicates the left side, for Indicates that on the right side, the origin of the cage center coordinate system coincides with the cage center of mass. The axis of rotation of the shaft and the cage coincide; is the center coordinate system of the cage pocket, the origin of the coordinate system of the center coordinate system of the cage pocket is Side cage The centers of the pockets coincide, Along the radial direction of the cage, Along the circumferential direction of the cage, Along the bearing axis; is the steel ball center coordinate system, the coordinate origin of the steel ball center coordinate system is Side bearing The centers of the steel balls coincide, Along the radial direction of the bearing, Along the bearing circumference, Along the axial direction of the bearing.

[0050] S2. Establish a non-rigid connection between the flywheel frame coordinate system and the space fixed coordinate system.

[0051] When establishing a non-rigid connection between the flywheel frame coordinate system and the space fixed coordinate system, set the stiffness coefficient and damping coefficient of the connection between the flywheel frame coordinate system and the space fixed coordinate system, such as Figure 4 As shown. Among them, the support base is bound to the fixed coordinate system in space, and the support base marking point 、 、 and , by setting 、 、 and middle 、 、 The displacement stiffness and damping in the direction, as well as the rotational stiffness and damping, are related to the flywheel frame 、 、 and Perform non-rigid connection. The calculation formula is as follows: , (1) Where, 、 and for 、 、 The force component value in the axial direction; 、 and for 、 、 The value of the moment component in the axial direction; 、 and for 、 、 Axis direction 、 Relative displacement component value between marker points; 、 and for 、 、 Axis direction 、 The relative rotation displacement component value between the marked points; 、 are the stiffness coefficient and the damping coefficient; 、 and for 、 、 Axis direction 、 The relative moving speed component value between the marking points; , and for 、 、 Axis direction 、 The relative rotational velocity component value between the marked points; 、 and for 、 、 The initial force component value in the axial direction; 、 and for 、 、 The initial moment component value in the axial direction.

[0052] S3. Establish a disturbance calculation model for the flywheel bearing assembly.

[0053] The disturbance calculation model of the flywheel bearing assembly established in S3 includes the calculation of the disturbance force and disturbance torque of the reaction flywheel bearing assembly. The disturbance calculation model of the flywheel bearing assembly is as follows: Figure 5 The calculation formula is as follows: , (2) Where, is the overall mass of the reaction flywheel assembly; 、 is the lever arm; is the coordinate point number; 、 、 The coordinate point of the frame base is 、 、 acceleration in direction; 、 and For the framework 、 、 Directional disturbance force; 、 and For the framework 、 、 Disturbance torque in direction; , and is the torque compensation coefficient.

[0054] S4. Analyze the relative position and motion relationship between the various parts of the flywheel bearing assembly.

[0055] When analyzing the relative position and motion relationship between the various parts in the bearing, the motion relationship of the steel ball and retainer relative to the outer ring, the motion relationship of the steel ball and retainer relative to the inner ring, and the motion relationship of the steel ball relative to the center of the retainer pocket are determined based on boundary conditions such as load, speed and bearing size parameters.

[0056] like Figure 6 As shown, the steel ball and the groove surface are in contact The linear speeds relative to the cage are: , (3) , (4) Where, is the rotation speed of the steel ball; is the diameter of the steel ball; is the pitch diameter; is the attitude angle; is the external contact angle; is the speed of the outer ring relative to the cage; Is a dimensionless parameter, and the calculation formula is: , (5) . (6) At the contact point Department: , (7) Where, is the speed of the inner ring relative to the cage, and its relationship expression is: , (8) Where, is the angular velocity of the inner circle; is the angular velocity of the outer ring.

[0057] From equations (5) and (6), we can get: , (9) . (10) According to equations (9) and (10), the angular velocity of the steel ball is: , (11) In the formula, if the inner and outer ring speed directions are different, take “ ”, otherwise, take “ ”.

[0058] , (12) , (13) . (14) Substituting equations (13) and (14) into equation (12), the angular velocity of the steel ball can be obtained as: , (15) In the formula, if the inner and outer ring speed directions are different, take “ ”, otherwise, take “ ”.

[0059] S5. Analyze the interaction forces between the parts in the bearing and establish a nonlinear equation group for each part, including the nonlinear differential equation group for the steel ball, cage and ring.

[0060] A nonlinear equation group for each component in the bearing is established. Based on the relative positions and motion parameters of each component in the bearing, the interaction forces between the rolling elements and raceways, the rolling elements and cages, and the cages and rings in the bearing are solved, and a nonlinear differential equation group is established.

[0061] The normal contact force between the steel ball and the channel is: , (16) Where, is the elastic approach; is the load-deformation constant.

[0062] like Figures 7-9 As shown, according to the center of the pocket With the center of the steel ball The motion state is judged based on the contact situation. is the center of the cage pocket, is the center of the steel ball, and When overlapping, the steel ball has no contact with the cage; Ahead When the ball drives the cage, the pocket center moves ; Hysteresis When the cage pushes the steel ball, the center of the pocket moves .

[0063] The contact force between the steel ball and the cage is: , (17) Where, is the linear approximation constant, ; is the load-deformation constant; is the cage pocket clearance, , is the cage pocket diameter; For the The displacement of the center of a pocket relative to the center of the steel ball in the hole.

[0064] According to the motion analysis of the steel ball in the bearing and the contact relationship between the steel ball and the inner and outer rings, the force analysis of the steel ball in the bearing is carried out, such as Figures 10-11 shown.

[0065] Subscript in the figure 、 Indicates inner and outer channels; Indicates the steel balls; 、 are the major and minor axes of the contact ellipse between the steel ball and the channel. 、 For the Normal contact load between a steel ball and the inner and outer grooves; 、 、 and For the The drag force between the steel ball and the channel in the contact area; 、 For the Normal contact load between each steel ball and the cage in the tangential and axial directions; 、 For the The inertial force component of the steel ball; 、 For the The rolling friction resistance on the surface of a steel ball; 、 For the Sliding friction resistance on the surface of a steel ball; 、 and For the The component of the inertial force acting on each steel ball; 、 and For the The Coriolis inertial force on the steel ball.

[0066] According to the motion and force analysis, the steel ball is in the inertial coordinate system The nonlinear dynamic differential equations under are as follows: , (18) Where, 、 and The steel ball in the steel ball coordinate system The moment of inertia along the three coordinate axes; 、 For the steel ball 、 moment of inertia in the direction; 、 and is the angular velocity of the steel ball , , Directional component; 、 、 is the angular velocity of the steel ball , , The angular acceleration component in the direction.

[0067] Outer ring in inertial coordinate system The nonlinear dynamic differential equations under are as follows: , (19) Where, is the overall mass of the outer ring; 、 and is the acceleration of the center of mass of the outer ring; 、 and is the principal moment of inertia of the outer ring as a whole; 、 and is the angular velocity of the outer ring as a whole; 、 is the overall angular acceleration of the outer ring; is the outer ring groove curvature radius coefficient.

[0068] Figures 12-14 middle, 、 The amount of movement of the steel ball relative to the center of the pocket in the tangential and axial directions; is the pocket width; is the radius of the steel ball; is the eccentricity of the cage; 、 are the two components of the offset; is the cage azimuth angle; For the A steel ball azimuth; is the minimum oil film thickness; 、 and To guide the normal, tangential forces and moments acting on the cage by the rib; It is the axial friction force of the cage and the guide ring on the cage when the cage contacts the guide ring; is the normal force between the steel ball and the cage pocket in the direction of movement; is the axial force between the steel ball and the cage pocket; 、 The friction force between the steel ball and the cage pocket in the axial and radial directions of the cage; is the cage gravity; For the cage in the reference coordinate system The angular velocity below.

[0069] The nonlinear dynamic differential equation of the cage is as follows: , (20) Where, For the quality of the cage; 、 and is the acceleration of the cage center of mass; 、 and is the moment of inertia of the cage; 、 The gyroscopic moment on the cage is 、 Directional component; 、 、 、 、 and are the components of the Coriolis inertia force and the drag inertia force acting on the cage in the non-inertial coordinate system; 、 and is the cage angular acceleration; is the number of steel balls; is the pitch diameter of the steel ball.

[0070] S6. According to the bearing dynamics theory, the GSTIFF variable step size integration algorithm is used to solve the nonlinear dynamic differential equations.

[0071] S7. Determine whether the dynamic equations have converged, that is, whether they meet the convergence error.

[0072] The default convergence error is When the system residual meets the preset convergence criteria, the spatial position and motion state parameters of the bearing assembly at the current moment are recorded; if the error exceeds the convergence threshold, the current calculation result is used as the initial condition to update to the next iteration step, and the differential equation is solved again. The calculation process is as follows: Figure 1 shown.

[0073] The simulation results are as follows: According to the above parameter settings and simulation steps in this embodiment, the simulation results refer to Figure 15 and 16 . Figure 15 The dynamic characteristics of the flywheel bearing assembly during operation include the collision force between the steel ball and the cage and the change of the cage speed over time. Figure 16 For the flywheel bearing assembly as a whole 、 、 The disturbance force and disturbance torque of the axis.

[0074] The analysis results show the intrinsic connection between the dynamic characteristics of the bearing and the disturbance characteristics of the flywheel assembly, and the stability of the bearing operation affects the stability of the disturbance output of the flywheel assembly.

[0075] It should be noted that the present invention demonstrates its main features, basic principles and advantages, but is not limited to the above-mentioned model. The present invention is also applicable to flywheel bearing assemblies with the same structure but different dimensional parameters and can be analyzed and calculated. It should be clear to those skilled in the art that the present invention is not limited to the above-mentioned embodiments. The description of the embodiments and the specification is only used to illustrate the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention will produce various changes and improvements according to actual conditions, and these changes and improvements all fall within the scope of protection claimed by the present invention. The scope of protection of the present invention is determined by the attached claims and their equivalents.

Claims

1. A method for analyzing disturbance characteristics of a flywheel bearing assembly, characterized by: The following steps are included: S1. Obtaining the initial conditions of the flywheel bearing assembly and establishing a flywheel bearing assembly coordinate system, including a space fixed coordinate system and a flywheel frame coordinate system; S2. Establishing an elastic damping connection model between the flywheel frame coordinate system and the space fixed coordinate system; S3. Establish a disturbance calculation model for the flywheel bearing assembly; S4. Analyze the relative position and motion relationship between the parts in the flywheel bearing assembly; S5. Analyze the interaction forces between the parts in the flywheel bearing assembly and establish a set of nonlinear dynamic differential equations for each part; S6. Based on the bearing dynamics theory, the Gstiff variable step size integration algorithm is used to solve the nonlinear dynamic differential equations; S7. Determine whether the dynamic equations satisfy the convergence error.

2. A method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, characterized in that: In step S1, the initial conditions of the flywheel bearing assembly obtained include flywheel frame structural parameters, bearing center distance, bearing structural parameters, material properties, lubrication parameters and operating conditions.

3. The method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, wherein: In step S1, the established flywheel bearing assembly coordinate system also includes the bearing outer ring coordinate system, the bearing cage center coordinate system, the cage pocket center coordinate system and the steel ball center coordinate system.

4. A method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, characterized in that: In step S2, when establishing the elastic damping connection model between the flywheel frame coordinate system and the space fixed coordinate system, the stiffness coefficient and damping coefficient of the connection between the flywheel frame coordinate system and the space fixed coordinate system are set.

5. The method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, wherein: In step S3, the established flywheel bearing assembly disturbance calculation model includes a calculation model of the reaction flywheel bearing assembly disturbance force and disturbance torque.

6. A method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, characterized in that: In step S4, the motion relationship between the steel ball and the cage relative to the outer ring, the motion relationship between the steel ball and the cage relative to the inner ring, and the motion relationship between the steel ball and the cage pocket center are determined based on the load, speed and bearing size parameters.

7. A method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, characterized in that: In step S5, based on the relative positions and motion parameters of the bearing components in S4, the interaction forces between the rolling elements and raceways, the rolling elements and the cage, and the cage and the rings in the bearing are solved, and a group of nonlinear dynamic differential equations is established.

8. The method for analyzing disturbance characteristics of a flywheel bearing assembly according to claim 1, wherein: In step S7, the convergence error is preset to When the system residual meets the preset convergence threshold, the spatial position and motion state parameters of the bearing assembly at the current moment are recorded. When the system residual exceeds the convergence threshold, the current calculation result is used as the initial condition to update to the next iteration step, and the differential equation is solved again.

Images