A rotor fault drop full process response simulation method, device and storage medium

Abstract

The present invention relates to the technical field of fault simulation, and discloses a method, device and storage medium for simulating the whole process of rotor fault drop. The present invention performs response simulation on the whole process of the fault drop of the rotor to be responded to in terms of tribology, dynamics and thermal multi-physics fields and their coupling only by combining the initial rotor parameter set of the rotor to be responded to, the initial operating condition data set of the drop and the coordinates and shapes of each component of the corresponding electrical system. The simulation process is simplified and the response simulation of the whole process of the fault drop of the rotor to be responded to is realized, which provides a theoretical basis for the research on the rotor drop behavior and the research on the failure process of the protective bearing.

Description

Technical Field

[0001] The present invention relates to the technical field of fault simulation, and in particular to a method, device and storage medium for simulating the full-process response of a rotor falling due to a fault. Background Art

[0002] As a typical representative of new types of suspension bearings, magnetic bearings have great advantages in terms of precision, efficiency, maintenance and application scenarios. They are often used in special environments such as high speed, vacuum, and ultra-clean. In order to avoid direct contact between the rotor and the stator when the electrical system fails, a protective bearing needs to be installed to provide temporary support for the rotor. When a rotor fails and falls, it will continuously collide, bounce and rub in the gap of the protective bearing, causing severe force and thermal effects that make the protective bearing prone to serious failure damage that is different from the fatigue failure of conventional rolling bearings, and even cause damage to the entire rotor system, resulting in a large amount of economic losses. At present, the reliability problem of the protective bearing has become one of the main bottlenecks restricting the current application and development of magnetic bearings.

[0003] In order to clarify the failure process of the protective bearing, it is first necessary to model the behavior of the protective bearing when the rotor fails and falls. At present, the rotor fall modeling method has made some progress, but there are still several shortcomings: (1) The modeling mainly includes dynamics and thermals, and lacks a description of the contact friction and wear behavior between the rotor and the protective bearing; (2) It only focuses on the state response at the beginning of the rotor fall, and lacks a study of the response of the entire fall process; (3) There is a lack of coupling between dynamic behavior and other physical fields. The above shortcomings limit the guiding significance of simulation calculations for the actual rotor fall and protective bearing failure process. Summary of the invention

[0004] In view of this, the present invention provides a method, device and storage medium for simulating the full process response of a rotor fault drop, so as to solve the problem that the shortcomings of the current rotor drop modeling method limit the guiding significance of simulation calculations for the actual rotor drop and protective bearing failure process.

[0005] In a first aspect, the present invention provides a method for simulating the whole process response of a rotor failure drop, the method comprising:

[0006] The initial rotor parameter set and the initial drop condition data set of the simulated rotor to be responded to are obtained, and the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded to are established according to the initial rotor parameter set and the initial drop condition data set; the whole process of the fault drop of the simulated rotor to be responded to is simulated according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, and the initial fault drop whole process response simulation is obtained, and the initial fault drop whole process response simulation result is obtained, and the initial fault drop whole process response simulation result includes multiple response simulation results of the simulated rotor to be responded to in terms of tribology, dynamics and thermal multi-physical fields and their coupling; it is judged whether the speed of the simulated rotor to be responded to is zero; when the speed of the simulated rotor to be responded to is zero, the initial fault drop whole process response simulation result is determined as the target fault drop whole process response simulation result of the simulated rotor to be responded to; when the speed of the simulated rotor to be responded to is not zero, the drop instantaneous condition data set of the simulated rotor to be responded to is obtained, and it is iterated repeatedly until the speed of the simulated rotor to be responded to is zero, and the target fault drop whole process response simulation result is obtained.

[0007] The method for simulating the whole process of rotor fault drop response provided by the present invention only combines the initial rotor parameter set of the simulated rotor to be responded to, the initial operating condition data set of the drop, and the coordinates and shapes of the corresponding electrical system components, and respectively performs response simulation on the whole process of the fault drop of the simulated rotor to be responded to in terms of tribology, dynamics and thermal multi-physical fields and their coupling. While simplifying the simulation process, the response simulation of the whole process of the fault drop of the simulated rotor to be responded to is realized, which provides a theoretical basis for the study of rotor drop behavior and the study of the failure process of protective bearings.

[0008] In an optional implementation, according to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, a response simulation is performed on the entire fault drop process of the rotor to be simulated, and an initial fault drop process response simulation result is obtained, including:

[0009] According to the initial rotor parameter set, the initial operating condition data set of the drop, and the coordinates and shapes of each component, the wear amount, total frictional heat generation, material thermal expansion effect, and lubricant viscosity-temperature effect between the simulated rotor to be responded to and the corresponding protective bearing are calculated respectively; according to the initial rotor parameter set, the initial operating condition data set of the drop, and the coordinates and shapes of each component, the rotor operation trajectory and the protective bearing operation trajectory are obtained through calculation using the preset Lagrangian equation; according to the wear amount, rotor operation trajectory, protective bearing operation trajectory, total frictional heat generation, material thermal expansion effect, and lubricant viscosity-temperature effect, the simulation results of the entire process response to the initial fault drop are determined.

[0010] The rotor fault drop full process response simulation method provided by the present invention first calculates the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded to and the corresponding protective bearing by combining the initial rotor parameter set of the simulated rotor to be responded to, the initial operating condition data set of the drop, and the coordinates and shapes of the corresponding electrical system components. Further, the rotor running trajectory and the protective bearing running trajectory can be calculated in combination with the preset Lagrangian equation. Finally, the initial fault drop full process response simulation result can be obtained by combining the calculated wear amount, rotor running trajectory, protective bearing running trajectory, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect, thereby realizing the response simulation of the full process of the fault drop of the simulated rotor to be responded to.

[0011] In an optional implementation, the wear amount, total friction heat, material thermal expansion effect, and lubricant viscosity-temperature effect between the simulated rotor to be responded to and the corresponding protective bearing are calculated based on the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, including:

[0012] According to the coordinates and shapes of each component, determine whether there is overlap between the components; if there is overlap between the components, obtain multiple overlapping positions, multiple overlapping depths and the relative movement speed at the two interfaces corresponding to each overlapping position; calculate the wear amount according to the multiple overlapping positions, multiple overlapping depths and each relative movement speed; calculate the total friction heat generation according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component; establish the temperature nodes at the key positions of each component according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component; calculate the temperature change of each temperature node; calculate the material thermal expansion effect and the lubricant viscosity-temperature effect according to the temperature change of each temperature node.

[0013] The rotor fault drop whole process response simulation method provided by the present invention can calculate the wear amount between the simulated rotor to be responded to and the corresponding protective bearing through the overlap between the components. Furthermore, the material thermal expansion effect and the lubricant viscosity-temperature effect can be determined through the temperature change of each temperature node, thus realizing the response simulation of the whole process of the fault drop of the simulated rotor to be responded to.

[0014] In an optional embodiment, the wear amount is calculated according to a plurality of overlapping positions, a plurality of overlapping depths and respective relative movement speeds, including:

[0015] The multiple overlapping positions and multiple overlapping depths are calculated by the Lankarani-Nikravesh contact force model to obtain multiple normal contact forces; based on the multiple normal contact forces and the relative motion speeds, the multiple friction forces are calculated by the Coulomb friction model to obtain multiple friction forces; based on the multiple normal contact forces and the multiple friction forces, the wear amount is calculated by the Archard wear model.

[0016] In an optional implementation, calculating the temperature change of each temperature node includes:

[0017] The heat transfer impedance between each temperature node is calculated according to the transfer relationship of each temperature node; based on the heat transfer impedance between each temperature node, the temperature change of each temperature node is obtained through calculation of the thermal differential equation.

[0018] In a second aspect, the present invention provides a device for simulating the whole process response of a rotor failure drop, the device comprising:

[0019] The acquisition module is used to acquire the initial rotor parameter set and the initial drop condition data set of the simulated rotor to be responded to, and establish the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded to according to the initial rotor parameter set and the initial drop condition data set; the response simulation module is used to perform response simulation on the fault drop process of the simulated rotor to be responded to according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, and obtain the initial fault drop process response simulation result, which includes multiple response simulation results of the simulated rotor to be responded to in terms of tribology, dynamics and thermal multi-physical fields and their coupling; the judgment module is used to judge whether the speed of the simulated rotor to be responded to is zero; the determination module is used to determine the initial fault drop process response simulation result as the target fault drop process response simulation result of the simulated rotor to be responded to when the speed of the simulated rotor to be responded to is zero; the iteration module is used to obtain the fall instantaneous condition data set of the simulated rotor to be responded to when the speed of the simulated rotor to be responded to is not zero, and repeatedly iterate until the speed of the simulated rotor to be responded to is zero, and obtain the target fault drop process response simulation result.

[0020] In an optional implementation, the response simulation module includes:

[0021] The first calculation submodule is used to calculate the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded and the corresponding protection bearing according to the initial rotor parameter set, the initial operating condition data set of the drop and the coordinates and shapes of each component; the second calculation submodule is used to obtain the rotor operation trajectory and the protection bearing operation trajectory through the preset Lagrangian equation according to the initial rotor parameter set, the initial operating condition data set of the drop and the coordinates and shapes of each component; the determination submodule is used to determine the initial fault drop full process response simulation result according to the wear amount, rotor operation trajectory, protection bearing operation trajectory, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect.

[0022] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for simulating the full process response of a rotor fault drop according to the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0023] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for simulating the full process response of a rotor fault drop according to the first aspect or any corresponding embodiment thereof.

[0024] In a fifth aspect, the present invention provides a computer program product, comprising computer instructions for causing a computer to execute the method for simulating the full process response of a rotor fault drop according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0026] Figure 1 is a flow chart of a method for simulating the whole process of rotor failure drop response according to an embodiment of the present invention;

[0027] Figure 2 is a flow chart of another method for simulating the whole process of rotor failure drop response according to an embodiment of the present invention;

[0028] Figure 3 It is a flow chart of a simulation calculation method for predicting the whole process response of a magnetic bearing rotor failure drop according to an embodiment of the present invention;

[0029] Figure 4 is a structural block diagram of a rotor fault drop full process response simulation device according to an embodiment of the present invention;

[0030] Figure 5 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0031] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0032] An embodiment of the present invention provides a method for simulating the response of a rotor during the entire process of a rotor failure drop. By performing response simulation on the entire process of a rotor failure drop to be simulated in terms of tribology, dynamics, thermal multi-physical fields and their coupling, the simulation process is simplified while the response simulation of the entire process of a rotor failure drop to be simulated is achieved, providing a theoretical basis for the study of rotor drop behavior and the study of the failure process of protective bearings.

[0033] According to an embodiment of the present invention, an embodiment of a method for simulating the full process response of a rotor fault drop is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0034] In this embodiment, a method for simulating the whole process response of a rotor failure drop is provided, which can be used in electronic devices such as computers, mobile phones, tablet computers, etc. Figure 1 FIG. 4 is a flow chart of a method for simulating the whole process of rotor failure drop response according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0035] Step S101, obtaining an initial rotor parameter set and a drop initial operating condition data set of a simulated rotor to be responded to, and establishing coordinates and shapes of various components of an electrical system corresponding to the simulated rotor to be responded to according to the initial rotor parameter set and the drop initial operating condition data set.

[0036] The initial rotor parameter set may include parameters such as rotor geometry parameters and protective bearing geometry parameters.

[0037] Furthermore, when a rotor fails and falls, collisions, bounces, and frictions will continuously occur in the gap of the protective bearing. Therefore, the initial operating condition data set of the fall may include: rotor speed, gap between the rotor and the protective bearing, collision speed, relative motion speed at the contact interface between the rotor and the protective bearing, contact relative sliding distance, rotor kinetic energy, rotor potential energy and other data.

[0038] Specifically, the coordinates and shapes of various components of the electrical system in three-dimensional space can be further established based on the obtained initial rotor parameter set of the simulated rotor to be responded to and the initial drop condition data set.

[0039] Step S102, based on the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, a response simulation is performed on the entire fault drop process of the rotor to be responded to, to obtain the initial fault drop process response simulation result.

[0040] The initial fault drop whole process response simulation results may include multiple response simulation results of the rotor to be responded to in terms of tribology, dynamics, thermal multi-physics fields and their coupling.

[0041] Specifically, the whole process of the fault drop of the simulated rotor to be responded to is simulated in terms of tribology, dynamics, thermal multi-physics fields and their coupling, realizing the response simulation of the whole process of the fault drop of the simulated rotor to be responded to, and solving the problem that the shortcomings of the current rotor drop modeling method limit the guiding significance of simulation calculations on the actual rotor drop and protective bearing failure process.

[0042] Step S103, determining whether the rotation speed of the simulated rotor to be responded to is zero.

[0043] Specifically, after the entire response simulation of the fault drop of the simulated rotor to be responded to is completed, it is determined whether the rotation speed of the simulated rotor to be responded to is zero at this time.

[0044] Step S104: when the rotation speed of the simulated rotor to be responded is zero, the initial fault drop full process response simulation result is determined as the target fault drop full process response simulation result of the simulated rotor to be responded.

[0045] Specifically, if the speed of the simulated rotor to be responded to is zero at this time, the response simulation ends, and the obtained initial fault drop full process response simulation result is used as the final target fault drop full process response simulation result.

[0046] Step S105, when the rotation speed of the simulated rotor to be responded is not zero, obtain the fall real-time operating condition data set of the simulated rotor to be responded, and iterate repeatedly until the rotation speed of the simulated rotor to be responded is zero, and obtain the full process response simulation result of the target fault fall.

[0047] Specifically, if the rotation speed of the simulated rotor to be responded to is not zero at this time, the real-time working condition data set of the simulated rotor to be responded to in the current falling state is continuously acquired.

[0048] Furthermore, based on the obtained real-time operating condition data set and combined with the initial rotor parameter set, the response simulation of the entire fault drop process of the simulated rotor to be responded to is repeated until the simulation is stopped when the speed of the simulated rotor to be responded to is zero, and the corresponding target fault drop entire process response simulation result is output.

[0049] The method for simulating the whole process of rotor fault drop response provided in the present embodiment only combines the initial rotor parameter set of the simulated rotor to be responded to, the initial operating condition data set of the drop, and the coordinates and shapes of the corresponding electrical system components to perform response simulation on the whole process of the fault drop of the simulated rotor to be responded to in terms of tribology, dynamics and thermal multi-physics fields and their coupling. While simplifying the simulation process, it realizes the response simulation of the whole process of the fault drop of the simulated rotor to be responded to, providing a theoretical basis for the study of rotor drop behavior and the study of protection bearing failure process.

[0050] In this embodiment, a method for simulating the whole process response of a rotor failure drop is provided, which can be used in electronic devices such as computers, mobile phones, tablet computers, etc. Figure 2 FIG. 4 is a flow chart of a method for simulating the whole process of rotor failure drop response according to an embodiment of the present invention. Figure 2 As shown, the process includes the following steps:

[0051] Step S201, obtain the initial rotor parameter set and the initial drop condition data set of the simulated rotor to be responded, and establish the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded according to the initial rotor parameter set and the initial drop condition data set. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0052] Step S202, based on the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, a response simulation is performed on the entire fault drop process of the rotor to be responded to, to obtain the initial fault drop process response simulation result.

[0053] Specifically, the above step S202 includes:

[0054] Step S2021, based on the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, respectively calculate the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded to and the corresponding protective bearing.

[0055] Among them, the material thermal expansion effect refers to the changes in the spatial position and gap caused by the thermal expansion of the material of the simulated rotor to be responded to; the lubricant viscosity-temperature effect refers to the changes in the lubrication performance of the lubricant caused by the viscosity-temperature effect.

[0056] Furthermore, in order to avoid direct contact between the rotor and the stator when the electrical system fails, a protective bearing needs to be installed to provide temporary support for the rotor. When the rotor fails and falls, it will continuously collide, bounce and rub in the gap of the protective bearing. Therefore, this embodiment further calculates the wear amount and total friction heat between the simulated rotor to be responded and the corresponding protective bearing.

[0057] In some optional implementations, the above step S2021 includes:

[0058] Step a1, judging whether there is overlap between the components according to the coordinates and shapes of the components.

[0059] Specifically, the geometric relationship between the components can be obtained according to the coordinates and shapes of the components, and then it is determined whether there is overlap between the components according to the geometric relationship, and the overlapping parts are determined to be in contact.

[0060] Step a2: if there is overlap between the components, obtain multiple overlapping positions, multiple overlapping depths and the relative movement speed at two interfaces corresponding to each overlapping position.

[0061] Specifically, if there is overlap between components, the contact positions and depths between the components are acquired, that is, multiple overlap positions and multiple overlap depths are acquired.

[0062] Furthermore, based on the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, the relative motion velocity v at the two interfaces of each basic position between the components can be calculated. s .

[0063] Step a3, calculating the wear amount according to the multiple overlapping positions, the multiple overlapping depths and the relative movement speeds.

[0064] In some optional implementations, the above step a3 includes:

[0065] Step a31, calculating the multiple overlapping positions and the multiple overlapping depths through the Lankarani-Nikravesh contact force model to obtain multiple normal contact forces.

[0066] Step a32, based on multiple normal contact forces and each relative motion speed, multiple friction forces are calculated by Coulomb friction model.

[0067] Step a32, based on multiple normal contact forces and multiple friction forces, the wear amount is calculated by using the Archard wear model.

[0068] First, based on multiple overlapping positions and multiple overlapping depths, the Lankarani-Nikravesh contact force model can be used to calculate the normal contact force F at each contact position. n , as shown in the following relation (1):

[0069]

[0070] Where: K represents the contact stiffness coefficient; δ represents the penetration depth; n represents the contact coefficient; c e represents the material restitution coefficient; represents the initial collision velocity. Among them, K, δ, and n can be determined by the overlapping position and overlapping depth.

[0071] Secondly, the Coulomb friction model is used to calculate the friction force F for contacts with non-zero relative motion speed, that is, with tangential relative sliding. t , as shown in the following relation (2):

[0072]

[0073] Where: μ represents the friction coefficient.

[0074] Finally, the Archard wear model can be used to calculate the material contact wear depth caused by contact relative sliding and friction, that is, the wear amount h w , as shown in the following relation (3):

[0075]

[0076] Where: k represents the friction coefficient; s represents the relative sliding distance of contact; H represents the Brinell hardness value of the material; S c Represents the contact area.

[0077] Step a4, calculating the total friction heat according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component.

[0078] Specifically, the total friction heat generated may be composed of radial contact heat generated by the rotor and the protective bearing, axial contact heat generated by the rotor and the protective bearing, and internal heat generated by the protective bearing.

[0079] Step a5, establishing the temperature nodes at the key positions of each component according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component.

[0080] Specifically, the temperature nodes at key positions of various components may include: a rotor radial contact surface, a rotor axial contact surface, a protective bearing inner ring, a protective bearing rolling element, a protective bearing outer ring and a bearing seat.

[0081] Step a6, calculating the temperature change of each temperature node.

[0082] In some optional implementations, the above step a6 includes:

[0083] Step a61, calculating the heat transfer impedance between each temperature node according to the transfer relationship of each temperature node.

[0084] Step a62, based on the heat transfer impedance between each temperature node, the temperature change of each temperature node is obtained through calculation using a thermal differential equation.

[0085] Specifically, the heat transfer resistance between each temperature node can be calculated according to the transfer relationship of each temperature node.

[0086] Furthermore, a heat transfer grid model is established, and the temperature change ΔQ of each node is calculated using the first-order thermal differential equation.

[0087] The thermal differential equation is shown in the following equation (4):

[0088]

[0089] In the formula: m represents mass; C p represents specific heat capacity; T represents the node temperature of each temperature node; ΔQ represents the heat flux of each temperature node, which is equal to the difference between the incoming heat and the outgoing heat, and is the ratio of the node temperature difference to the heat transfer impedance between the nodes.

[0090] Step a7, calculating the material thermal expansion effect and the lubricant viscosity-temperature effect according to the temperature change of each temperature node.

[0091] Specifically, according to the temperature change of each temperature node, the influence of the temperature change of each component on the motion state can be calculated, which may include:

[0092] (1) Material thermal expansion effect: The thermal expansion of the material causes changes in the spatial position and clearance, affecting the rotor drop dynamics behavior and the friction and wear in the protective bearing. The material thermal expansion ε can be calculated using the following relationship (5):

[0093]

[0094] In the formula: ε represents the thermal expansion; ξ represents the thermal expansion coefficient; γ represents the Poisson's ratio of the material; d represents the diameter; the subscripts i, o, and b represent the inner ring, outer ring, and rolling element of the protective bearing, respectively.

[0095] Among them, ε i This will directly cause the change in the clearance between the rotor and the protective bearing, and there is also an induced load F caused by thermal expansion in the protective bearing. th , as shown in the following relation (6):

[0096] F th =K[ε b -0.5(ε i -ε o )cosα b ] 3 / 2 (6)

[0098] (2) Lubricant viscosity-temperature effect: The lubricant is affected by the viscosity-temperature effect, which causes changes in lubrication performance. That is, the lubricant viscosity-temperature effect refers to the fact that the viscosity of the lubricant decreases sharply as the temperature rises. Furthermore, the relationship between temperature and viscosity can be described by the Reynolds viscosity-temperature equation, as shown in the following equation (7):

[0099]

[0100] Where: η represents the viscosity at temperature T; η 0 Indicates temperature as T 0 The viscosity at β 0 Represents the viscosity-temperature coefficient.

[0101] Step S2022, based on the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, the rotor running trajectory and the protective bearing running trajectory are calculated through the preset Lagrangian equation.

[0102] Specifically, the Lagrangian equation of the rotor is established and solved, and the coordinates and speed parameters of each degree of freedom of the rotor are calculated to form the running trajectory and state of the rotor. The Lagrangian equation is shown in the following relation (8):

[0103]

[0104] Where: T 1 represents the rotor kinetic energy; V represents the rotor potential energy; q i represents generalized coordinates; Q i Represents generalized force.

[0105] Further, according to the stress state of the inner ring of the protection bearing, the friction torque M of the inner ring of the protection bearing can be calculated. The friction torque M is used to calculate the rotation state and rotation speed of the protection bearing.

[0106] Furthermore, the Lagrangian equation of the protective bearing is established and solved, as shown in the above relationship (8), and the coordinates and speed parameters of each degree of freedom of the inner ring and outer ring of the protective bearing are calculated to form the running trajectory and state of the protective bearing.

[0107] Step S2023, determining the simulation result of the entire process response of the initial fault drop according to the wear amount, rotor running track, protection bearing running track, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect.

[0108] Specifically, once the wear amount, rotor running trajectory, protective bearing running trajectory, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect are all calculated, the first response simulation of the simulated rotor to be responded to can be completed and the corresponding initial fault drop full process response simulation results can be formed.

[0109] Step S203, determine whether the speed of the simulated rotor to be responded to is zero. Figure 1 Step S103 of the illustrated embodiment will not be described in detail here.

[0110] Step S204: when the speed of the simulated rotor to be responded is zero, the initial fault drop full process response simulation result is determined as the target fault drop full process response simulation result of the simulated rotor to be responded. Figure 1 Step S104 of the illustrated embodiment will not be described in detail here.

[0111] Step S205: when the speed of the simulated rotor to be responded is not zero, obtain the fall real-time working condition data set of the simulated rotor to be responded, and iterate repeatedly until the speed of the simulated rotor to be responded is zero, and obtain the full process response simulation result of the target fault fall. Figure 1 Step S105 of the illustrated embodiment will not be described in detail here.

[0112] The rotor fault drop response simulation method provided in this embodiment first calculates the wear amount and total friction heat between the simulated rotor to be responded and the corresponding protective bearing by combining the initial rotor parameter set of the simulated rotor to be responded, the initial operating condition data set of the drop, and the coordinates and shapes of the corresponding electrical system components. At the same time, the wear amount between the simulated rotor to be responded and the corresponding protective bearing can be calculated by the overlap between the components. Further, the material thermal expansion effect and the lubricant viscosity-temperature effect can be determined by the temperature change of each temperature node. Further, the rotor running trajectory and the protective bearing running trajectory can be calculated in combination with the preset Lagrangian equation. Finally, the initial fault drop response simulation result can be obtained by combining the calculated wear amount, rotor running trajectory, protective bearing running trajectory, total friction heat, material thermal expansion effect and lubricant viscosity-temperature effect. Therefore, by implementing the present invention, while simplifying the simulation process, the response simulation of the fault drop of the simulated rotor to be responded to is realized, which provides a theoretical basis for the study of rotor drop behavior and the study of protective bearing failure process.

[0113] In one example, a simulation calculation method for predicting the full-process response of a magnetic bearing rotor failure drop is provided, such as Figure 3 As shown, including:

[0114] Step 1: Input the initial rotor system parameters and initial drop conditions, including rotor geometry parameters, protection bearing geometry parameters, rotor and protection bearing clearance, rotor speed and other parameters, and set the simulation single-step calculation step length to 1×10-6s;

[0115] Step 2: Establish the coordinates and shapes of each component in three-dimensional space, determine whether there is overlap between components based on geometric relationships, determine whether there is contact, and calculate the normal contact force F at each contact position based on the contact position and depth. n , where F n The Lankarani-Nikravesh contact force model shown in the above relationship (1) is used for calculation;

[0116] Step 3: Calculate the relative velocity of the two interfaces at each contact position, and calculate the friction force F for contacts where the relative velocity is not zero, that is, there is tangential relative sliding. t , where F t The Coulomb friction model shown in the above relationship (2) is used for calculation;

[0117] Step 4: Calculate the material contact wear depth h caused by contact relative sliding and friction w , where h w Calculated using the Archard wear model shown in the above relationship (3);

[0118] Step 5: Establish and solve the Lagrangian equation of the rotor, calculate the coordinates and speed parameters of each degree of freedom of the rotor, and form the running trajectory and state of the rotor; wherein the Lagrangian equation is shown in the above relationship (8);

[0119] Step 6: Calculate the friction torque on the inner ring of the protective bearing according to the stress state of the inner ring of the protective bearing, wherein the friction torque is used to calculate the rotation state and rotation speed of the protective bearing;

[0120] Step 7: Establish and solve the Lagrangian equation of the protective bearing, which is the same as step 5, and calculate the coordinates and velocity parameters of each degree of freedom of the inner and outer rings of the protective bearing;

[0121] Step 8: Calculate the heat generated by the friction process between the components of the system: the total heat generated is composed of the radial contact heat between the rotor and the protective bearing, the axial contact heat between the rotor and the protective bearing, and the internal heat of the protective bearing;

[0122] Step 9: Establish the temperature nodes at the key positions of each component, including: the rotor radial contact surface, the rotor axial contact surface, the protective bearing inner ring, the protective bearing rolling element, the protective bearing outer ring and the bearing seat, and calculate the heat transfer impedance between each temperature node according to the transfer relationship of each temperature node;

[0123] Step 10: Establish a system heat transfer grid model and use the first-order thermal differential equation to calculate the temperature change of each node, where the thermal differential equation is shown in the above relationship (4);

[0124] Step 11: Calculate the effect of temperature changes of each component of the system on the motion state, which can be divided into two aspects:

[0125] (1) Thermal expansion of materials causes changes in spatial position and gap;

[0126] (2) The lubricant is affected by the viscosity-temperature effect, causing changes in lubrication performance;

[0127] Step 12: Determine whether the rotor speed drops to zero at this time. If true, end the calculation and output all the results of the current calculation. If not, return to step 2 to continue the calculation and solution.

[0128] Among them, the calculated coupling effects of material thermal expansion effect and lubricant viscosity-temperature effect in the dynamic model and contact model are as follows:

[0129] The coupling between system thermals and dynamics is reflected in the gap change caused by thermal expansion of the material (dynamics → thermal). This gap change will cause changes in the motion state and contact relationship of each component, thereby causing changes in dynamics (thermal → dynamics). The coupling between thermals and contact friction is reflected in the influence of the viscosity-temperature effect of the lubricant on the friction force (thermal → tribology). The change in friction force caused by the viscosity-temperature effect brings about changes in contact heating (tribology → thermal).

[0130] This example provides a simulation calculation method for predicting the response of the entire process of a magnetic bearing rotor failure falling. It can simulate the axial trajectory of each degree of freedom of the rotor during the fall, the contact force of the protective bearing, the friction and wear of the rotor and the protective bearing, and the heat and temperature rise between the components. By reasonably simplifying the calculation process, the simulation capability of the entire process from the occurrence of the fall to the rotor shutdown is achieved, providing characteristic parameters for the failure research of the protective bearing.

[0131] In this embodiment, a rotor fault drop whole process response simulation device is also provided, which is used to implement the above embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0132] This embodiment provides a rotor failure drop whole process response simulation device, such as Figure 4 As shown, including:

[0133] The acquisition module 401 is used to acquire the initial rotor parameter set and the initial drop condition data set of the simulated rotor to be responded to, and establish the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded to according to the initial rotor parameter set and the initial drop condition data set.

[0134] The response simulation module 402 is used to perform response simulation on the entire fault drop process of the rotor to be responded to according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, and obtain the response simulation results of the entire initial fault drop process. The response simulation results of the entire initial fault drop process include multiple response simulation results of the rotor to be responded to in terms of tribology, dynamics and thermal multi-physical fields and their coupling.

[0135] The judgment module 403 is used to judge whether the rotation speed of the simulated rotor to be responded is zero.

[0136] The determination module 404 is used to determine the initial fault drop full process response simulation result as the target fault drop full process response simulation result of the simulated rotor to be responded to when the speed of the simulated rotor to be responded to is zero.

[0137] The iteration module 405 is used to obtain the real-time operating condition data set of the simulated rotor to be responded to when the speed of the simulated rotor to be responded to is not zero, and repeatedly iterate until the speed of the simulated rotor to be responded to is zero, so as to obtain the simulation result of the whole process of the target fault drop response.

[0138] In some optional implementations, the response simulation module 402 includes:

[0139] The first calculation submodule is used to calculate the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded and the corresponding protective bearing according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component.

[0140] The second calculation submodule is used to obtain the rotor running trajectory and the protection bearing running trajectory through calculation using a preset Lagrangian equation based on the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component.

[0141] The submodule is used to determine the simulation results of the whole process response of the initial fault drop according to the wear amount, rotor running track, protection bearing running track, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect.

[0142] In some optional implementations, the first calculation submodule includes:

[0143] The judging unit is used to judge whether there is overlap between the components according to the coordinates and shapes of the components.

[0144] The acquisition unit is used to acquire a plurality of overlapping positions, a plurality of overlapping depths and a relative movement speed at two interfaces corresponding to each overlapping position if there is overlap between the components.

[0145] The first calculation unit is used to calculate the wear amount according to a plurality of overlapping positions, a plurality of overlapping depths and each relative movement speed.

[0146] The second calculation unit is used to calculate the total friction heat according to the initial rotor parameter set, the initial operating condition data set of the drop, and the coordinates and shapes of each component.

[0147] A unit is established to establish temperature nodes at key positions of each component according to an initial rotor parameter set, an initial drop condition data set, and coordinates and shapes of each component.

[0148] The third calculation unit is used to calculate the temperature change of each temperature node.

[0149] The fourth calculation unit is used to calculate the material thermal expansion effect and the lubricant viscosity-temperature effect according to the temperature change of each temperature node.

[0150] In some optional implementations, the first computing unit includes:

[0151] The first calculation subunit is used to calculate the multiple overlapping positions and the multiple overlapping depths through the Lankarani-Nikravesh contact force model to obtain multiple normal contact forces.

[0152] The second calculation subunit is used to obtain multiple friction forces based on multiple normal contact forces and various relative motion speeds through Coulomb friction model calculation.

[0153] The third calculation subunit is used to obtain the wear amount based on multiple normal contact forces and multiple friction forces through Archard wear model calculation.

[0154] In some optional implementations, the third computing unit includes:

[0155] The fourth calculation subunit is used to calculate the heat transfer impedance between each temperature node according to the transfer relationship of each temperature node.

[0156] The fifth calculation subunit is used to obtain the temperature change of each temperature node through calculation based on the heat transfer impedance between each temperature node through thermal differential equation.

[0157] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0158] The rotor fault drop full process response simulation device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0159] The embodiment of the present invention also provides a computer device having the above Figure 4 The rotor failure drop full process response simulation device shown.

[0160] See also Figure 5 , Figure 5 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Figure 5As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 5 A processor 10 is taken as an example.

[0161] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0162] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0163] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0164] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive; the memory 20 may also include a combination of the above types of memory.

[0165] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0166] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0167] A part of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the existence of the computer program instruction in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc., and accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium accessible to the computer.

[0168] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for simulating the whole process response of a rotor failure drop, characterized in that: The method comprises: Acquire an initial rotor parameter set and a drop initial operating condition data set of the simulated rotor to be responded to, and establish the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded to according to the initial rotor parameter set and the drop initial operating condition data set; According to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, a response simulation is performed on the entire fault drop process of the simulated rotor to be responded to, to obtain an initial fault drop process response simulation result, wherein the initial fault drop process response simulation result includes multiple response simulation results of the simulated rotor to be responded to in terms of tribology, dynamics, thermal multi-physics fields and their coupling; Determining whether the rotation speed of the simulated rotor to be responded to is zero; When the rotation speed of the simulated rotor to be responded is zero, the initial fault drop full process response simulation result is determined as the target fault drop full process response simulation result of the simulated rotor to be responded; When the rotation speed of the simulated rotor to be responded is not zero, the fall instantaneous operating condition data set of the simulated rotor to be responded is obtained, and the process is iterated repeatedly until the rotation speed of the simulated rotor to be responded is zero, thereby obtaining the full process response simulation result of the target fault fall.

2. The method according to claim 1, characterized in that According to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, a response simulation is performed on the fault drop process of the simulated rotor to be responded to, and a response simulation result of the initial fault drop process is obtained, including: According to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, respectively calculate the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded and the corresponding protection bearing; According to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component, a rotor running trajectory and a protective bearing running trajectory are obtained through calculation using a preset Lagrangian equation; The simulation result of the full process response of the initial fault drop is determined according to the wear amount, the rotor running track, the protective bearing running track, the total friction heat generation, the material thermal expansion effect and the lubricant viscosity-temperature effect.

3. The method according to claim 2, characterized in that According to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, the wear amount, total friction heat, material thermal expansion effect, and lubricant viscosity-temperature effect between the simulated rotor to be responded to and the corresponding protective bearing are calculated, including: Determine whether there is overlap between the components according to the coordinates and shapes of the components; If there is overlap between the components, a plurality of overlapping positions, a plurality of overlapping depths and a relative movement speed at two interfaces corresponding to each overlapping position are obtained; calculating the wear amount according to the plurality of overlapping positions, the plurality of overlapping depths and the respective relative movement speeds; Calculating the total friction heat according to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component; Establishing temperature nodes at key positions of each component according to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component; Calculate the temperature change of each temperature node; According to the temperature change of each temperature node, the thermal expansion effect of the material and the viscosity-temperature effect of the lubricant are calculated.

4. The method according to claim 3, characterized in that Calculating the wear amount according to the plurality of overlapping positions, the plurality of overlapping depths and the relative movement speeds includes: Calculate the multiple overlapping positions and the multiple overlapping depths using a Lankarani-Nikravesh contact force model to obtain multiple normal contact forces; Based on the multiple normal contact forces and the relative motion speeds, multiple friction forces are obtained through Coulomb friction model calculation; The wear amount is obtained by calculating the Archard wear model based on the multiple normal contact forces and the multiple friction forces.

5. The method according to claim 3, characterized in that: Calculate the temperature changes of each temperature node, including: Calculate the heat transfer impedance between each temperature node according to the transfer relationship of each temperature node; According to the heat transfer impedance between each temperature node, the temperature change of each temperature node is obtained through calculation of thermal differential equation.

6. A rotor fault drop full process response simulation device, characterized in that: The device comprises: An acquisition module is used to acquire an initial rotor parameter set and a falling initial working condition data set of the simulated rotor to be responded to, and establish the coordinates and shapes of each component of the electrical system corresponding to the simulated rotor to be responded to according to the initial rotor parameter set and the falling initial working condition data set; a response simulation module, for performing response simulation on the entire fault drop process of the simulated rotor to be responded to according to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component, to obtain a response simulation result of the entire initial fault drop process, wherein the response simulation result of the entire initial fault drop process includes multiple response simulation results of the simulated rotor to be responded to in terms of tribology, dynamics, thermal multi-physics fields and their coupling; A judgment module, used for judging whether the speed of the simulated rotor to be responded is zero; A determination module, used for determining the initial fault drop full process response simulation result as the target fault drop full process response simulation result of the simulated rotor to be responded to when the speed of the simulated rotor to be responded to is zero; The iteration module is used to obtain the actual operating condition data set of the fall of the simulated rotor to be responded to when the rotation speed of the simulated rotor to be responded to is not zero, and repeatedly iterate until the rotation speed of the simulated rotor to be responded to is zero, so as to obtain the simulation result of the whole process response of the target fault fall.

7. The device according to claim 6, characterized in that The response simulation module comprises: A first calculation submodule is used to calculate the wear amount, total friction heat generation, material thermal expansion effect and lubricant viscosity-temperature effect between the simulated rotor to be responded and the corresponding protection bearing according to the initial rotor parameter set, the initial drop condition data set and the coordinates and shapes of each component; A second calculation submodule is used to obtain a rotor running trajectory and a protection bearing running trajectory through calculation using a preset Lagrangian equation according to the initial rotor parameter set, the initial drop condition data set, and the coordinates and shapes of each component; A determination submodule is used to determine the simulation result of the whole process response of the initial fault drop according to the wear amount, the rotor running track, the protection bearing running track, the total friction heat generation, the material thermal expansion effect and the lubricant viscosity-temperature effect.

8. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the rotor fault drop full process response simulation method according to any one of claims 1 to 5 by executing the computer instructions.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the rotor fault drop full process response simulation method according to any one of claims 1 to 5.

10. A computer program product, characterized in that It comprises computer instructions, and the computer instructions are used to enable a computer to execute the rotor fault drop full process response simulation method according to any one of claims 1 to 5.

Images