Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a gear pump maintenance method based on digital twins.
The purpose of the invention is realized by the following technical scheme.
In one aspect of the invention, a predictive maintenance method based on a digital twin gear pump comprises the following steps:
measuring and receiving working condition parameters of the gear pump in real time, wherein the working condition parameters comprise working pressure p, flow Q, rotating speed n and torque MpAnd a temperature T;
calculating volumetric efficiency and limit speed of the gear pump based on the operating condition parameters to assess the operating condition of the gear pump and make a prediction of failure, wherein,
in the formula: omega
minIs the limit rotation speed; p is the working pressure; mu is the dynamic viscosity of the oil liquid; c is the clearance between the bearing and the bearing bush; r
aIs the bearing radius; eta
VIs volume efficiency; v
2,eIs the displacement at the working pressure; v
2,iDisplacement at no load pressure; q. q.s
V2,eIs the flow at the working pressure; q. q.s
V2,IIs the flow at no load pressure; n is
eIs the rotation speed under the working pressure; n is
iIs the rotation speed at the no-load pressure,
calculating the central positions of a driving gear and a driven gear of the gear pump and the abrasion loss of a friction pair, wherein the central positions are as follows:
in the formula: m is1Is the mass of the gear 1; f. ofb1Is the bearing counter force borne by the gear 1; f. ofp1Is the pressure to which the gear 1 is subjected; f. ofmgThe meshing force applied to the gear; j. the design is a square1Is the moment of inertia of gear 1; theta1Is the angle of rotation of the gear 1; mp1Torque generated by the pressure to which the gear 1 is subjected; mmTorque generated for the motor; r isb1Is the radius of the meshing point of the gear 1; m is2Is the mass of the gear 2; f. ofb2Is the bearing counter force borne by the gear 2; f. ofp2Is the pressure to which the gear 2 is subjected; j. the design is a square2Is the moment of inertia of gear 2; theta2Is the angle of rotation of the gear 2; mp2Torque generated by the pressure to which the gear 2 is subjected; r isb2Is the radius of the meshing point of the gear 2,
the wear amount of the friction pair comprises the end face of the gear and the floatingThe abrasion loss among the movable shaft sleeve, the tooth top and the inner wall of the pump shell, the abrasion loss between the end surface of the gear and the floating shaft sleeve determine the size of a gap hf between the end surface of the gear and the floating shaft sleeve, and the abrasion loss between the tooth top and the inner wall of the pump shell determine the size of a gap hi between the pump shell of the gear pump and the tooth top of the gear; the abrasion loss of the friction pair is as follows:
in the formula: q is the abrasion loss of the friction pair; e is the energy consumed by wear per stroke; n is the number of strokes; s is the sliding distance of each stroke; c is the wear system constant; tau is
maxThe maximum shear force to which the part is subjected;
obtaining a flow rate through the gap based on the pressure differential and the size of the gap, and obtaining a pressure for each chamber based on a relationship between the gap and the flow rate, wherein,
in the formula: b is the tooth width; omega is the gear pump rotating speed; r isextIs the addendum circle radius; bfIs the tooth height; r ismIs the average tooth height; mu is dynamic viscosity of oil liquid; ltThe tooth crest thickness; lfThe tooth thickness of the pitch circle; p is a radical ofiIs the ith chamber pressure; p is a radical ofi-1Is the i-1 th chamber pressure; qh,iThe tooth crest clearance flow is obtained; qf,iIs the tooth flank clearance flow;
under adiabatic, isentropic conditions, and pi+1>piThe pressure in the ith control body can be obtained according to the continuity equation:
ΔQi=Qi+1-Qi,
in the formula: b isoilIs the bulk modulus of elasticity; omega is the gear pump rotating speed; viIs the ith control volume; theta is a gear rotation angle; delta QiIs the amount of change in the flow in the ith control fluid; qi+1The flow rate of the ith +1 th control body flowing into the ith control body through the gap is set; qiThe flow rate of the ith control body flowing out through the gap is measured;
modeling an inlet control body, an equal-volume control body, an outlet control body and a closed cavity control body in sequence by using the above formula, and obtaining the pressure of the whole internal flow field by using M, wherein the model is as follows:
in the formula: p is a radical ofn(θ)|1The pressure of the nth chamber of the gear 1; p is a radical ofm(θ)|2Pressure of the mth cavity of the gear 2; p is a radical oft(theta) is the pressure of the gear pump sealed cavity; p is a radical ofout(θ) is gear pump outlet pressure; p is a radical ofin(θ) is gear pump inlet pressure;
in the gear pump external flow field, the Fourier expansion of the pressure pulsation signal at the measuring point is as follows:
in the formula: a. theoLIs a direct current component; a. the1Lcos(ωoLt+ψ1) Is the fundamental component of the frequency conversion; n is not less than 2, AnLcos(nωoLt+ψn) Is a harmonic component of the frequency conversion; n-1 … 6; a. theoHIs a constant; a. themHGenerating a pressure pulsation signal amplitude for the tooth frequency; k is a radical ofamIs an amplitude modulation coefficient; k is a radical offmIs a frequency modulation coefficient; zetamIs a phase angle; omegaoHIs the tooth frequency; m is 1 … 6; p'outL(t) is an amplitude modulated signal; p ″)outL(t) is a frequency-modulated signal,
and establishing a digital twin model of the gear pump, and taking the outlet pressure pulsation of the gear pump as a monitoring signal of the predictive maintenance model.
In the method, during the expansion of the pressure pulsation signal at the measuring point of the gear pump external flow field, time domain analysis, frequency domain analysis, wavelet analysis, functional analysis or matrix analysis is carried out on the measured data and the simulation data.
In the method, the time domain analysis comprises time domain statistical analysis or correlation analysis, and the frequency domain analysis comprises FFT, coherent analysis, cepstrum analysis or autoregressive spectrum analysis.
In the method, the central positions of the driving gear and the driven gear and the abrasion loss of the friction pair are calculated by solving a differential equation system.
In the method, a differential equation set of pressure pulsation is solved through a four-order and five-order Runge-Kutta single-step algorithm in an internal flow field of the gear pump.
In the method, in the external flow field of the gear pump, a homogeneous turbulence spectrum analysis method is adopted, energy distribution of vortexes of various scales in turbulence at a measuring point of the gear pump under different working conditions is learned through a neural network, and a basis function coefficient of a heuristic function space is corrected, so that parameters of a pressure pulsation mathematical model are optimized.
In the method, the flow between the gear pump end surface gap and the gear crest gap is equivalent to the Couett-Poiseuille flow between the flat plates, and the flow passing through the gap is as follows:
in the formula: b is the tooth width; omega is the gear pump rotating speed; r isextIs the addendum circle radius; bfIs the tooth height; r ismIs the average tooth height; mu is dynamic viscosity of oil liquid; ltThe tooth crest thickness; lfThe tooth thickness of the pitch circle; p is a radical ofiIs the ith chamber pressure; p is a radical ofi-1Is the i-1 th chamber pressure; qh,iThe tooth crest clearance flow is obtained; qf,iIs the flank clearance flow.
The above description is only an overview of the technical solutions of the present invention, and in order to make the technical means of the present invention more clearly apparent, and to make the implementation of the content of the description possible for those skilled in the art, and to make the above and other objects, features and advantages of the present invention more obvious, the following description is given by way of example of the specific embodiments of the present invention.
Detailed Description
Specific embodiments of the present invention will be described in more detail below with reference to the accompanying drawings fig. 1 to 4. While specific embodiments of the invention are shown in the drawings, it should be understood that the invention may be embodied in various forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
It should be noted that certain terms are used throughout the description and claims to refer to particular components. As one skilled in the art will appreciate, various names may be used to refer to a component. This specification and claims do not intend to distinguish between components that differ in name but not function. In the following description and in the claims, the terms "include" and "comprise" are used in an open-ended fashion, and thus should be interpreted to mean "include, but not limited to. The description which follows is a preferred embodiment of the invention, but is made for the purpose of illustrating the general principles of the invention and not for the purpose of limiting the scope of the invention. The scope of the present invention is defined by the appended claims.
For the purpose of facilitating understanding of the embodiments of the present invention, the following description will be made by taking specific embodiments as examples with reference to the accompanying drawings, and the drawings are not to be construed as limiting the embodiments of the present invention.
For better understanding, fig. 1 is a schematic diagram of a digital twin based gear pump predictive maintenance, and as shown in fig. 1, a digital twin based gear pump maintenance method comprises the following steps:
measuring and receiving working condition parameters of the gear pump in real time, wherein the working condition parameters comprise working pressure p, flow Q, rotating speed n and torque MpThe temperature T;
calculating volumetric efficiency and limit speed of the gear pump based on the operating condition parameters to assess the operating condition of the gear pump and make a prediction of failure, wherein,
in the formula: omega
minIs the limit rotation speed; p is the working pressure; mu is the dynamic viscosity of the oil liquid; c is the clearance between the bearing and the bearing bush; r
aIs the bearing radius; eta
VIs volume efficiency; v
2,eIs the displacement at the working pressure; v
2,iDisplacement at no load pressure; q. q.s
V2,eIs the flow at the working pressure; q. q.s
V2,IIs the flow at no load pressure; n is
eIs the rotation speed under the working pressure; n is
iIs the rotation speed at the no-load pressure,
calculating the central positions of a driving gear and a driven gear of the gear pump and the abrasion loss of a friction pair, wherein the central positions are as follows:
in the formula: m is1Is the mass of the gear 1; f. ofb1Is the bearing counter force borne by the gear 1; f. ofp1Is the pressure to which the gear 1 is subjected; f. ofmgThe meshing force applied to the gear; j. the design is a square1Is the moment of inertia of gear 1; theta1Is the angle of rotation of the gear 1; mp1Torque generated by the pressure to which the gear 1 is subjected; mmTorque generated for the motor; r isb1Is the radius of the meshing point of the gear 1; m is2Is the mass of the gear 2; f. ofb2Is the bearing counter force borne by the gear 2; f. ofp2Is the pressure to which the gear 2 is subjected; j. the design is a square2Is the moment of inertia of gear 2; theta2Is the angle of rotation of the gear 2; mp2Torque generated by the pressure to which the gear 2 is subjected; r isb2Is the radius of the meshing point of the gear 2,
the abrasion loss of the friction pair comprises abrasion among the gear end, the floating shaft sleeve, the tooth top and the inner wall of the pump shell, and the abrasion loss of the friction pair is as follows:
in the formula: q is the abrasion loss of the friction pair; e is the energy consumed by wear per stroke; n is the number of strokes; s is the sliding distance of each stroke; c is the wear system constant; tau is
maxThe maximum shear force to which the part is subjected;
in the internal flow field of the gear pump, a gear pump pressure pulsation model is established based on a centralized parameter method, an inlet control body, an isometric control body, an outlet control body and a closed cavity control body are respectively modeled,
in the formula: p is a radical ofn(θ)|1The pressure of the nth chamber of the gear 1; p is a radical ofm(θ)|2Pressure of the mth cavity of the gear 2; p is a radical oft(theta) is the pressure of the gear pump sealed cavity; p is a radical ofout(θ) is gear pump outlet pressure; p is a radical ofin(θ) is gear pump inlet pressure;
in the gear pump external flow field, the Fourier expansion of the pressure pulsation signal at the measuring point is as follows:
in the formula: a. theoLIs a direct current component; a. the1Lcos(ωoLt+ψ1) Is the fundamental component of the frequency conversion; n is not less than 2, AnL cos(nωoLt+ψn) Is a harmonic component of the frequency conversion; n-1 … 6; a. theoHIs a constant; a. themHGenerating a pressure pulsation signal amplitude for the tooth frequency; k is a radical ofamIs an amplitude modulation coefficient; k is a radical offmIs a frequency modulation coefficient; zetamIs a phase angle; omegaoHIs the tooth frequency; m is 1 … 6; p'outL(t) is an amplitude modulated signal; p ″)outLAnd (t) is a frequency modulated signal.
In a preferred embodiment of the method, time domain analysis, frequency domain analysis, wavelet analysis, functional analysis or matrix analysis is performed on the measured data and the simulation data.
In a preferred embodiment of the method, the time domain analysis comprises time domain statistical analysis or correlation analysis, and the frequency domain analysis comprises FFT, coherence analysis, cepstral analysis or autoregressive spectral analysis.
In the preferred embodiment of the method, the central positions of the driving gear and the driven gear and the abrasion amount of the friction pair are calculated by solving a differential equation system.
In the preferred embodiment of the method, the differential equation system of the pressure pulsation is solved by ODE45 in the gear pump internal flow field.
In the preferred embodiment of the method, in the gear pump external flow field, a homogeneous turbulence spectrum analysis method is adopted, and energy distribution of vortexes of various scales in turbulence at a measuring point of the gear pump under different working conditions is learned through a neural network, so that a basis function coefficient of a tentative function space is corrected, and parameters of a pressure pulsation mathematical model are optimized.
In one embodiment, a digital twin based gear pump maintenance method comprises:
constructing a working condition parameter updating module, wherein the pressure sensor, the flow sensor and the rotating speed-torque sensor are arranged on the physical experiment table, and the temperature sensor respectively obtains the pressure p, the flow Q, the rotating speed n and the torque MpTemperature T, etc.
Constructing a state evaluation module, wherein the volumetric efficiency eta of the gear pump is calculated according to the actually measured flow of the flow sensorVAnd obtaining the limit rotation speed omega by an expression between the pressure and the rotation speedminAnd the gear pump operation can be guided.
Constructing a model dynamic parameter calculation module, wherein a dynamic model of the gear pump is established, and the central positions of the driving gear and the driven gear are calculated by solving a differential equation set; the gear pump tooth crest and end surface abrasion mainly refers to abrasive wear, tooth crest and end surface abrasion amount is obtained according to an abrasion amount calculation formula and through combination of experimental measurement, and basic parameters of the gear pump are shown in table 1.
TABLE 1
And (2) constructing a gear pump internal flow field module, wherein a gear pump control body pressure pulsation equation set is established by adopting a centralized parameter method, the flow of tooth tops and end surfaces is assumed to be laminar flow between plates in the model, the influence of turbulent flow is considered by an inlet, an outlet and a closed cavity control body, and then the differential equation set of pressure pulsation is solved through ODE45 to obtain the pressure pulsation value in each control body.
And constructing a gear pump outflowing flow field module, wherein a homogeneous turbulence spectrum analysis method is adopted, pressure pulsation at a measuring point is approximately considered to be a stable random process, and the pressure pulsation is considered to be the superposition of vortexes with various different scales, and the energy distribution of the vortexes with various different scales in turbulence at the measuring point of the gear pump under different working conditions is learned through a neural network to correct a basis function coefficient of a heuristic function space, so that parameters of a pressure pulsation mathematical model are optimized.
And constructing a physical entity module, wherein a gear pump performance test bench is constructed according to the requirements of the gear pump performance test with reference to national standards.
FIG. 2 is a schematic diagram of a human-computer interaction interface based on the digital twin gear pump prediction maintenance method, as shown in FIG. 2, a working condition parameter module on the interface constantly reflects parameters such as pressure p, flow Q, rotating speed n, temperature T and the like of the operation of the gear pump; the performance evaluation module mainly displays the volumetric efficiency eta of the gear pumpVLimit rotation speed omegaminThe total working time h and the volumetric efficiency make a prediction evaluation on the failure of the gear pump, and the minimum rotating speed mainly prompts an operator to the minimum rotating speed at which the lubricating oil film of the gear pump bearing is damaged and abraded under a certain load, for example, when the load is 14MPa, the minimum rotating speed is 2000r/min, so that the abrasion of the gear pump bearing is reduced by increasing the rotating speed or reducing the load; the characteristic signal module mainly displays the information of the selected characteristic signals, such as pressure pulsation signals, flow signals, vibration signals and the like, and the change information of the characteristic signals can reflect the performance degradation and the occurrence of faults of the gear pump.
FIG. 3 is a time domain graph of characteristic signal pressure pulsation based on the digital twin gear pump predictive maintenance method according to the present invention, in which the gear pump speed was 2466.6r/min, the load was 10MPa, the temperature was 22 ℃, and the flow rate was 9866.4 mL/min. The figure shows that the obvious periodic fluctuation of the pressure pulsation signal is mainly caused by the periodic motion of the gear pump, and 12 small fluctuations are included in a large fluctuation, wherein the small fluctuation represents outlet pressure pulsation caused by oil removal when one tooth is rotated, and the large fluctuation represents pulsation caused by periodic change of leakage when the gear pump rotates for one circle. The amplitude of the pulsation is influenced by the rotating speed and the outlet pressure, and if the pulsation is worn and broken, the leakage amount is increased, so that the influence factor of large fluctuation is increased.
FIG. 4 is a frequency domain plot of the characteristic signal pressure pulsation of the predictive maintenance method based on the digital twin gear pump according to the present invention, the gear pump speed was 2466.6r/min, the load was 10MPa, the temperature was 22 deg.C, and the flow rate was 9866.4mL/min during the experiment. The figure shows that obvious periodic fluctuation components of the pressure pulsation signal mainly comprise gear rotating frequency, gear pump meshing frequency (tooth frequency) and frequency multiplication of the tooth frequency. The frequency conversion component is influenced most obviously by the rotating speed, the proportion of the frequency conversion component is gradually smaller along with the increase of the rotating speed, the tooth frequency gradually becomes the main pulse frequency, but the absolute amplitude of the frequency conversion basically keeps unchanged, the relative amplitude changes obviously, and the frequency conversion is reflected in a frequency domain graph to modulate the tooth frequency.
The method establishes a digital twin model of the gear pump, and takes the outlet pressure pulsation of the gear pump as a monitoring signal of a predictive maintenance model. Firstly, establishing a pressure pulsation model of an internal flow field of the gear pump by adopting a centralized parameter method to obtain the pressure in each control body inside the gear pump; secondly, a local homogeneous turbulence pulsation model is established for the position of a pressure sensor measuring point on an outlet pipeline of the gear pump, a mathematical model of pressure pulsation at the measuring point is obtained by adopting a homogeneous turbulence spectrum analysis method, energy distribution of various scale vortexes in turbulence at the measuring point of the gear pump under different working conditions is learned through a neural network, a basis function coefficient of a heuristic function space is corrected, high-fidelity simulation of the virtual model on a physical entity model is realized, and meanwhile, the virtual model makes an instructive strategy on the operation condition of the physical entity by calculating the volumetric efficiency and the limit rotating speed of the gear pump. By adopting the method, the real-time interaction between the gear pump virtual model and the physical entity can be realized, so that the composition of turbulence at the outlet measuring point of the gear pump is disclosed, and a new thought is provided for the performance monitoring and maintenance of the gear pump.
Although the embodiments of the present disclosure are described above with reference to the drawings, the technical solutions of the present disclosure are not limited to two operating conditions of different rotation speeds and different loads, and include other kinds of operating conditions. The particular embodiments disclosed above are illustrative and explanatory only and are not restrictive. Those skilled in the art, having the benefit of this disclosure, may effect numerous modifications thereto and changes may be made without departing from the scope of the disclosure as set forth in the claims that follow.