Parameterized modeling flow distribution pair leakage rapid calculation method and device

Abstract

The invention discloses a parameterized modeling flow distribution pair leakage rapid calculation method and device, and the method comprises the steps: firstly calculating the residual pressing force and residual overturning moment of a flow distribution pair according to the structural parameters of a plunger pump, and then solving through vector synthesis and a force balance equation set, and determining the thickness of a central oil film of the flow distribution pair, the overturning angle of the cylinder body and the overturning direction of the cylinder body, and finally substituting the parameters into a leakage flow calculation formula considering the overturning of the cylinder body to obtain the leakage flow of the flow distribution pair. According to the quick calculation method for leakage of the flow distribution pair based on parametric modeling, the applicability of a traditional analytical model under the simple working condition is reserved, the influence of cylinder overturning under the high-pressure and high-speed working condition is considered, the application range of the traditional model is expanded, and therefore quick and accurate estimation of the leakage flow of the flow distribution pair under the wide-area working condition is achieved, and the calculation efficiency is improved. And a rapid and accurate mathematical model is provided for high-precision and high-energy-efficiency control of the EHA system.

Description

Technical Field

[0001] The present invention belongs to the technical field of hydraulic pumps, and particularly relates to a rapid calculation method and device for the leakage of a port plate pair with parametric modeling. Background Art

[0002] As a new generation of servo actuator system, the Electro-Hydrostatic Actuator (EHA) has broad application prospects in the fields of aerospace, construction machinery, etc. due to its advantages of high efficiency, compactness, and plug-and-play. The EHA system adjusts the output flow by regulating the speed and direction of the motor pump, thereby realizing the position servo control of the actuator. This requires the piston pump to accurately perform the feed-forward compensation of the leakage flow under the working conditions in all four quadrants to ensure the dynamic performance and control accuracy of the EHA system. Among them, the port plate pair has the largest sealing size among all friction pairs and is the main source of the leakage flow of the piston pump, accounting for more than 60% of the total leakage. Moreover, under the working conditions of high pressure and high speed, the cylinder block is extremely prone to tipping, resulting in a sharp increase in leakage. Therefore, it is crucial to quickly and accurately calculate the leakage flow of the port plate pair.

[0003] The traditional analytical model assumes the port plate pair as a parallel gap with a constant height and derives a clear leakage flow expression. However, this requires an experienced engineer to set a reasonable gap height h to ensure the prediction accuracy. Especially under the working conditions of high pressure and high speed, the tipping of the cylinder block will cause a significant increase in the leakage flow, and the simple parallel slit flow assumption will no longer be applicable. Another method for calculating leakage is to establish a numerical calculation model based on the friction lubrication theory. This method takes into account the dynamic behavior of the cylinder block tipping under the working conditions of high pressure and high speed, as well as the complex thermo-elastohydrodynamic lubrication effect on the friction interface, and can more accurately reflect the non-linear leakage characteristics of the port plate pair under a wide range of working conditions. However, the calculation of the numerical model often involves multiple loop iterations and coupled solutions, which are time-consuming and occupy a large amount of computing resources, and are difficult to be directly applied to the control system of the EHA.

[0004] Therefore, there is an urgent need for a calculation method that can quickly and accurately estimate the leakage flow of the port plate pair under a wide range of working conditions. Summary of the Invention

[0005] The purpose of the present invention is to provide a rapid calculation method and device for the leakage of a port plate pair with parametric modeling in view of the deficiencies of the prior art, so as to realize the rapid and accurate estimation of the leakage flow of the port plate pair under a wide range of working conditions. This helps the control system to accurately and timely compensate for the flow loss, and further realize the high-precision and high-energy efficiency control of the EHA system.

[0006] The purpose of the present invention is achieved by the following technical solutions: In the first aspect, the present invention provides a rapid calculation method for the leakage of a port plate pair with parametric modeling, including the following steps:

[0007] (1) Obtain the structural parameters of the plunger pump;

[0008] (2) For the valve plate pair, adopt the design method of residual clamping force, and calculate the residual compressive force and the residual overturning moments in the horizontal and vertical directions on the valve plate pair;

[0009] (3) Through vector synthesis, determine the tilting orientation of the cylinder block based on the residual overturning moments in the horizontal and vertical directions; (4) Based on the residual compressive force and the residual overturning moments in the horizontal and vertical directions, construct the force balance equations when the valve plate pair overturns, and solve to obtain the central oil film thickness and the tilting angle

[0010] (5) Substitute the central oil film thickness, the tilting angle of the cylinder block, and the tilting orientation into the leakage flow rate calculation formula considering the tilting of the cylinder block to obtain the leakage flow rate of the valve plate pair.

[0011] Furthermore, the structural parameters of the plunger pump in step (1) include: the number of plungers z, the swash plate angle β, the plunger diameter d p , the plunger mass m p , the mass of the slipper m s , the high-pressure zone wrap angle Φ of the valve plate H , the radius R of the plunger distribution circle p , the inner diameter R1 of the inner seal band in the valve plate pair, the outer diameter R2 of the inner seal band in the valve plate pair, the inner diameter R3 of the outer seal band in the valve plate pair, the outer diameter R4 of the outer seal band in the valve plate pair, the distance L from the origin to the centroid of the plunger-slipper assembly at the inner dead center m0 , the distance L from the origin to the center of the plunger ball head at the inner dead center ss0 and the distance L from the origin to the center of the spline pair sc0 .

[0012] Furthermore, take the tilting orientation in step (3) as the y-axis, pointing from the origin to the direction of the minimum clearance height, take the cylinder block rotation center as the z-axis, pointing to the swash plate direction, and determine the x-axis according to the right-hand rule to establish a coordinate system. Calculate the reaction force and the reaction moment on the seal band in this coordinate system, and then construct the force balance equation between the reaction force and the residual compressive force on the valve plate pair and the moment balance equation between the reaction moment and the residual overturning moment on the valve plate pair.

[0013] Furthermore, the specific force balance equations when the valve plate pair overturns are constructed as follows:

[0014]

[0015]

[0016] In the formula, F rc and M rcThe reaction force and moment provided by the sealing belt when the cylinder block overturns, respectively. Solve the above force balance equations, and determine the central oil film thickness h0 and the overturning angle γ.

[0017] Further, the reaction force F rc and the reaction moment M rc provided by the sealing belt in step (4) are calculated as follows:

[0018]

[0019]

[0020] In the formula, R i and R o are the inner diameter and outer diameter of the sealing belt respectively, h max and h min are the maximum clearance height and minimum clearance height on the sealing belt respectively; among them,

[0021]

[0022]

[0023]

[0024]

[0025] In the formula, h0 is the central oil film thickness, γ is the cylinder block overturning angle, K is the reaction force stiffness, δ is the reaction force critical threshold. When the clearance height is lower than this value, it indicates that the oil film thickness is very small and the reaction force begins to be generated.

[0026] Further, the leakage flow rate of the flow distribution pair in step (5) includes the leakage flow rate Q Ls (R1, R2) of the inner seal and the leakage flow rate Q Ls (R3, R4) of the outer seal, where R1 is the inner diameter of the inner sealing belt of the flow distribution pair, R2 is the outer diameter of the inner sealing belt of the flow distribution pair, R3 is the inner diameter of the outer sealing belt of the flow distribution pair, and R4 is the outer diameter of the outer sealing belt of the flow distribution pair.

[0027] Further, the leakage flow rate Q Ls of the sealing belt considering the cylinder block overturning in step (5) is calculated as follows:

[0028]

[0029] In the formula, μ is the viscosity of the oil.

[0030] In a second aspect, the present invention further provides a rapid calculation device for the leakage of a flow distribution pair with parametric modeling, including a memory and one or more processors. Executable code is stored in the memory, and when the processor executes the executable code, the rapid calculation method for the leakage of a flow distribution pair with parametric modeling as described above is implemented.

[0031] In a third aspect, the present invention further provides a computer-readable storage medium, on which a program is stored. It is characterized in that when the program is executed by a processor, the rapid calculation method for the leakage of a flow distribution pair with parametric modeling as described above is implemented.

[0032] In a fourth aspect, the present invention further provides a computer program product, including a computer program. When the computer program is executed by a processor, the rapid calculation method for the leakage of a flow distribution pair with parametric modeling as described above is implemented.

[0033] The beneficial effects of the present invention are as follows:

[0034] 1. This method can quickly and accurately solve the leakage flow rate of the flow distribution pair under different working conditions, and the calculation efficiency is greatly improved compared with the numerical model, providing an accurate mathematical model for the real-time compensation of the flow loss in the EHA system.

[0035] 2. This method retains the applicability of the traditional analytical model under simple working conditions, and considers the influence of the cylinder block overturning under high-pressure and high-speed working conditions, expanding the application range of the traditional model, enabling it to more accurately reflect the non-linear leakage characteristics of the flow distribution pair under wide-range working conditions, and having a high prediction accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 is a flow chart of a rapid calculation method for the leakage of a flow distribution pair with parametric modeling according to the present invention;

[0037] Figure 2 is a force analysis diagram of the flow distribution pair;

[0038] Figure 3 is a schematic diagram of the cylinder block overturning orientation;

[0039] Figure 4 is the clearance height of the flow distribution pair under different overturning degrees;

[0040] Figure 5 is a comparison between the method proposed by the present invention and the numerical solution model.

[0041] Figure 6 is a schematic structural diagram of a rapid calculation device for the leakage of a flow distribution pair with parametric modeling provided by the present invention.

[0042] Reference numerals:

[0043] 1. Distribution plate, 2. Cylinder block, 3. Plunger, 4. Return spring plate, 5. Slipper, 6. Swash plate, 7. Central spring, 8. Discharge port of plunger pump, 9. Suction port of plunger pump. Detailed implementation manners

[0044] To more clearly illustrate the embodiments of the present invention, the present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.

[0045] As Figure 1 shown, the present invention provides a fast calculation method for the leakage of a distribution pair with parametric modeling, including the following steps.

[0046] Step 1: Input the structural parameters of the plunger pump. Taking a certain 7 mL / r plunger pump as an example, the main structural parameters input are shown in the following table:

[0047] Table 1 Main structural parameters of the plunger pump

[0048]

[0049] Step 2: Conduct a force analysis on the cylinder block to calculate the remaining clamping force F rz and the remaining overturning moment M rx and M ry .

[0050] The cylinder block rotates and slides on the surface of the distribution plate, forming a distribution pair interface. On the one hand, the high-pressure oil in the plunger chamber acts on the step surface at the bottom of the cylinder hole, causing the cylinder block to be axially pressed by the clamping force F cz and tightly pressed on the distribution plate. On the other hand, an oil film supporting force F f is generated between the high-speed rotating cylinder block and the distribution plate, pushing the cylinder block away from the distribution plate, so that while ensuring the close sealing between the cylinder block and the distribution plate, it also has good lubrication and load-bearing performance. However, once the action lines of F cz and F f do not coincide, a moment that causes the cylinder block to overturn will inevitably be formed. To ensure the normal and stable operation of the plunger pump, a design method of remaining clamping force needs to be adopted for the distribution pair. The axial force balance coefficient B b of the distribution pair is

[0051] (1)

[0052] where the balance coefficient B b takes 0.9 - 1, the oil film supporting force F f is slightly less than the clamping force F cz , generating a remaining clamping force F rz =F cz -F f and forming a counter-moment to balance its overturning moment, so that the cylinder block is in a balanced state.

[0053] As Figure 2 shown, a rectangular coordinate system O-X-Y-Z is established for force analysis. The origin O is located at the center of the port plate, the Z-axis coincides with the rotation center of the cylinder block and points to the swash plate direction, the Y-axis is vertically upward, and the X-axis points from the oil suction port of the port plate to the oil discharge port side. The central spring force F s and the plunger chamber pressure F pn act together on the cylinder block to form an axial compression force, and act on the port plate through the cylinder block, as shown in the following formula:

[0054] (2)

[0055] (3)

[0056] In the formula, d p is the plunger diameter, p n is the pressure of the nth plunger chamber. p n can be determined by the following formula:

[0057] (4)

[0058] In the formula, p H and p L are the oil discharge port and oil suction port pressures of the plunger pump respectively, θ n is the angular displacement of the nth plunger chamber. The compression force generated by the plunger chamber on the oil discharge side is greater than that on the oil suction side, resulting in the cylinder block tipping over to the high-pressure side of the port plate. The calculation of the corresponding hydraulic tipping moment is shown in the following formula:

[0059] (5)

[0060] (6)

[0061] In the formula, z is the number of plungers, R p is the radius of the plunger distribution circle. The plunger shoe assembly generates a rotational centrifugal force F cpsn while rotating with the cylinder block, and generates a reciprocating inertia force F ipsn while reciprocating in the cylinder bore. These two inertia forces act directly on the cylinder block through the plunger pair:

[0062] (7)

[0063] (8)

[0064] In the formula, m p and m sare the masses of the plunger and the slider respectively, β is the swash plate angle, and ω is the rotational angular velocity of the cylinder block. The plane where the center of mass of the plunger and slipper assembly lies is not parallel to the surface of the valve plate, resulting in the cylinder block tipping towards the outer dead center of the valve plate. The calculation of the corresponding inertial tipping moment is shown in the following formula:

[0065] (9)

[0066] (10)

[0067] In the formula, L m0 is the distance from the origin to the center of mass of the plunger and slipper assembly at the inner dead center. The high-pressure oil in the plunger cavity and the compressed spring press the plunger and slider assembly against the swash plate. The lateral force F ssyn exerted by the reaction of the swash plate on each plunger and slipper assembly is expressed as:

[0068] (11)

[0069] The lateral force acts on the cylinder block, and the resulting tipping moment is:

[0070] (12)

[0071] (13)

[0072] The lateral forces of all the plunger and slipper assemblies are offset by the spline reaction force F sc in the Y-axis direction:

[0073] (14)

[0074] In the formula, L ss0 is the distance from the origin to the center of the plunger ball head at the inner dead center. The tipping moment generated by the spline reaction force on the cylinder block is:

[0075] (15)

[0076] (16)

[0077] In the formula, L sc0 is the distance from the origin to the center of the spline. In summary, the pressing force F c and the tipping moment M t acting on the cylinder block are:

[0078] (17)

[0079] (18)

[0080] Consider the following mathematical approximation to simplify the above formula:

[0081] (19)

[0082] (20) (21)

[0083] (22)

[0084] And written in component form as:

[0085] (23)

[0086] (24)

[0087] (25)

[0088] Among them,

[0089] (26)

[0090] (27)

[0091] (28)

[0092] The static pressure oil film supporting force of the port plate on the cylinder block consists of three parts, including the static pressure supporting force of the outer seal band, the static pressure supporting force of the inner seal band, and the hydraulic supporting force of the high-pressure oil in the oil discharge waist-shaped groove on the cylinder block. On the oil discharge side, the calculation formula for the oil film pressure distribution of the port plate pair is:

[0093] (29)

[0094] In the formula, R1, R2, R3, and R4 are the inner diameter of the inner seal band, the outer diameter of the inner seal band, the inner diameter of the outer seal band, and the outer diameter of the port plate pair of the port plate pair respectively. Integrating the above formula, the supporting force F of the static pressure oil film of the port plate pair on the cylinder block is obtained fH and the anti-overturning moment M fHx and M fHy :

[0095] (30)

[0096] (31)

[0097] (32)

[0098] After arranging the above formula, it can be obtained:

[0099] (33)

[0100] (34)

[0101] (35)

[0102] In the formula, Φ H is the wrap angle of the high-pressure area of the port plate, and Φ L is the wrap angle of the low-pressure area of the port plate, and Φ L = 2π - Φ H . Among them,

[0103] (36)

[0104] (37)

[0105] To sum up, the remaining clamping force F rz received by the port plate pair is:

[0106] (38)

[0107] The remaining overturning moment M rx and M ry are respectively:

[0108] (39)

[0109] (40)

[0110] Step 3: Determine the overturning orientation ζ of the cylinder block through vector synthesis.

[0111] As Figure 3 shown, the combined action of the remaining overturning moments M rx and M ry determines the overturning orientation ζ of the cylinder block:

[0112] (41)

[0113] Step 4: Derive the analytical expression of the reaction force of the seal band of the port plate pair, construct a force balance equation set, and solve the central oil film thickness h0 and the overturning angle γ.

[0114] The unbalanced moment causes the cylinder block to overturn relative to the port plate, and solid contact occurs at the places where the oil film of the seal band is very thin, thus providing an additional reaction force to maintain the force balance of the cylinder block. Assume that there is a linear relationship between the clearance height of the port plate pair and the reaction force p r of the solid contact on the seal band, that is:

[0115] (42)

[0116] where h is the height of the clearance between the flow distribution pairs, K is the stiffness of the reaction force, and δ is the critical threshold of the reaction force. When the clearance height is lower than this value, it indicates that the oil film thickness is very small and the reaction force begins to be generated. The clearance height at any point (r, θ) on the sealing band is expressed as:

[0117] (43)

[0118] where h0 is the central oil film thickness and γ is the tilting angle of the cylinder block. A coordinate system o-x-y-z is defined to describe the reaction force generated on the sealing band, as Figure 3 shown. The z-axis is consistent with the Z-axis, the y-axis coincides with the tilting azimuth angle, points from the origin to the direction of the minimum clearance height, and the x-axis is determined according to the right-hand rule. The clearance height of the sealing band in the circumferential direction can be rewritten as:

[0119] (44)

[0120] where θ r is the angular displacement in the o-x-y-z coordinate system.

[0121] As Figure 4 shown, according to the severity of the cylinder block tilting, it can be divided into mild tilting, severe tilting, and complete tilting. For one of the sealing bands of the flow distribution pairs, when the cylinder block is mildly tilted, the tilting angle is small, and the clearance height on the sealing band is lower than the height threshold δ, that is:

[0122] (45)

[0123] where R m is the average radius of the sealing band. At this time, the reaction force F rc and the reaction moment M rc can be calculated as:

[0124] (46)

[0125] (47)

[0126] where R i and R o are the inner diameter and outer diameter of the sealing band respectively. When the cylinder block is severely tilted, the tilting angle is large, and the flow distribution pairs present a significant wedge-shaped clearance, and part of the clearance height will exceed the critical threshold, that is:

[0127] (48)

[0128] At this time, the reaction force F rc and the reaction moment M rc can be calculated as:

[0129] (49)

[0130] (50)

[0131] Where θ r1 and θ r2 are the angular displacements when the clearance height is equal to the critical threshold, and the calculation formula is:

[0132] (51)

[0133] When the cylinder block is completely overturned, the cylinder block will be completely separated from the port plate, and the clearance height is all above the critical threshold, that is:

[0134] (52)

[0135] At this time, the reaction force F rc and the reaction torque M rc are both 0. Therefore, the reaction forces and torques of the sealing strip under different degrees of cylinder overturning are:

[0136] (53)

[0137] (54)

[0138] Where h max is the maximum clearance height h max =h0 + R m γ, h min is the minimum clearance height h max =h0 - R m γ, where

[0139] (55)

[0140] (56)

[0141] (57)

[0142] (58)

[0143] To sum up, the force balance equation of the port pair includes force balance and torque balance, as shown in the following formula:

[0144] (59)

[0145] (60)

[0146] By solving the force balance equations, the central oil film thickness h0 and the tilting angle γ can be obtained.

[0147] Step 5: Derive the calculation formula for the leakage flow rate of the flow distribution pair considering the tilting of the cylinder block, and calculate the leakage flow rate of the flow distribution pair.

[0148] Considering that the leakage flow of the flow distribution pair is laminar flow in the gap, and assuming steady state, continuity, incompressibility, and neglecting the body force, according to the knowledge of fluid mechanics, the calculation formula for the leakage of the flow distribution pair along the circumferential direction of the sealing band is:

[0149] (61)

[0150] In the formula, μ is the viscosity of the oil. The total leakage flow rate passing through a single sealing band can be obtained by integrating the above formula:

[0151] (62)

[0152] After arrangement:

[0153] (63)

[0154] Therefore, the total leakage flow rate of the flow distribution pair is:

[0155] (64)

[0156] It is worth mentioning that the calculation formula for the leakage flow rate of the flow distribution pair considering the tilting of the cylinder block derived in the present invention is a more extensive development on the traditional classical analytical model. When the tilting angle is 0, that is, γ = 0, the above leakage flow rate calculation formula will degenerate into the traditional classical analytical model:

[0157] (65)

[0158] To further illustrate the superiority of the present invention, for the examples given in the present invention, the leakage flow rates of the flow distribution pair at different pressures and rotational speeds are calculated respectively by using the numerical model and the analytical model proposed in the present invention. The calculation results are as Figure 5 shown. It can be seen that the calculation results of the analytical model are in good agreement with the numerical model, and can more accurately reflect the non-linear leakage characteristics of the flow distribution pair under a wide range of working conditions. More importantly, the calculation time required for the numerical model to solve the leakage flow rate under one working condition is about 2.5 h. For the 12 rotational speed working conditions (1000 r / min to 12000 r / min) and 5 pressure working conditions (7 MPa to 35 MPa) given in the figure, the total calculation time reaches 150 h. However, by using the analytical method proposed in the present invention, the time required to solve the leakage flow rates of all working conditions is only 0.6 s, with extremely high calculation efficiency.

[0159] Corresponding to the embodiment of the rapid calculation method for the leakage of the port plate pair with parametric modeling described above, the present invention also provides an embodiment of a rapid calculation device for the leakage of the port plate pair with parametric modeling.

[0160] Referring to Figure 6 , an embodiment of a rapid calculation device for the leakage of the port plate pair with parametric modeling provided by an embodiment of the present invention includes a memory and one or more processors. Executable code is stored in the memory. When the processor executes the executable code, it is used to implement the rapid calculation method for the leakage of the port plate pair with parametric modeling in the above embodiment.

[0161] The embodiment of the rapid calculation device for the leakage of the port plate pair with parametric modeling provided by the present invention can be applied to any device with data processing capabilities. The device with data processing capabilities can be a device or apparatus such as a computer. The device embodiment can be implemented by software, or by hardware, or by a combination of software and hardware. Taking software implementation as an example, as a logically meaningful device, it is formed by the processor of any device with data processing capabilities reading the corresponding computer program instructions in the non-volatile memory into the memory for operation. From the hardware level, as Figure 6 shown, it is a hardware structure diagram of any device with data processing capabilities where the rapid calculation device for the leakage of the port plate pair with parametric modeling provided by the present invention is located. Except for Figure 6 the processor, memory, network interface, and non-volatile memory shown, the device with data processing capabilities where the device in the embodiment is located usually also includes other hardware according to the actual functions of the device with data processing capabilities, which will not be elaborated here.

[0162] The specific implementation process of the functions and roles of each unit in the above device can be specifically seen in the implementation process of the corresponding steps in the above method, which will not be elaborated here.

[0163] For the device embodiment, since it basically corresponds to the method embodiment, the relevant parts can refer to the partial description of the method embodiment. The device embodiments described above are only illustrative. The units described as separate components may or may not be physically separated. The components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of the present invention. Those of ordinary skill in the art can understand and implement it without creative work.

[0164] An embodiment of the present invention further provides a computer-readable storage medium, on which a program is stored. When the program is executed by a processor, a rapid calculation method for leakage of a flow distribution pair with parametric modeling in the above embodiment is implemented.

[0165] The computer-readable storage medium may be an internal storage unit of any device with data processing capabilities described in any of the foregoing embodiments, such as a hard disk or memory. The computer-readable storage medium may also be an external storage device of any device with data processing capabilities, such as a plug-in hard disk, a Smart Media Card (SMC), an SD card, a Flash Card, etc. equipped on the device. Further, the computer-readable storage medium may also include both an internal storage unit and an external storage device of any device with data processing capabilities. The computer-readable storage medium is used to store the computer program and other programs and data required by any device with data processing capabilities, and may also be used to temporarily store data that has been output or will be output.

[0166] The present invention also provides a computer program product, including a computer program. When the computer program is executed by a processor, a rapid calculation method for leakage of a flow distribution pair with parametric modeling is implemented.

[0167] The above-described embodiments are only descriptions of the preferred embodiments of the present invention, and do not limit the scope of the present invention. Without departing from the design spirit of the present invention, various deformations and improvements made by those of ordinary skill in the art to the technical solution of the present invention shall fall within the protection scope determined by the claims of the present invention.

Claims

1. A rapid calculation method for the leakage of a port plate pair with parametric modeling, characterized in that, Including the following steps: (1) Obtain the structural parameters of the plunger pump; (2) For the valve plate pair, adopt the design method of residual clamping force to calculate the residual compressive force and the residual overturning moments in the horizontal and vertical directions suffered by the valve plate pair; (3) Through vector synthesis, determine the tilting orientation of the cylinder block based on the residual overturning moments in the horizontal and vertical directions; (4) Based on the residual compressive force and the residual overturning moments in the horizontal and vertical directions, construct the force balance equations when the valve plate pair overturns, and solve to obtain the central oil film thickness and the tilting angle (5) Substitute the central oil film thickness, the tilting angle of the cylinder block and the tilting orientation into the leakage flow rate calculation formula considering the tilting of the cylinder block to obtain the leakage flow rate of the valve plate pair.

2. The rapid calculation method for the leakage of the port plate pair with parametric modeling according to claim 1, characterized in that: The structural parameters of the plunger pump in step (1) include: the number of plungers z, the swash plate angle β, the plunger diameter d p , the mass of the plunger m p , the mass of the slipper m s , the high-pressure zone wrap angle Φ of the valve plate H , the radius R of the plunger distribution circle p , the inner diameter R1 of the inner seal band of the valve plate pair, the outer diameter R2 of the inner seal band of the valve plate pair, the inner diameter R3 of the outer seal band of the valve plate pair, the outer diameter R4 of the outer seal band of the valve plate pair, the distance L from the origin to the centroid of the plunger-slipper assembly at the inner dead center m0 , the distance L from the origin to the center of the plunger ball head at the inner dead center ss0 and the distance L from the origin to the center of the spline pair sc0 .

3. The rapid calculation method for leakage of the port plate pair with parametric modeling according to claim 1, characterized in that: Take the tilting orientation in step (3) as the y-axis, pointing from the origin to the direction of the minimum clearance height, take the rotation center of the cylinder block as the z-axis, pointing to the swash plate direction, and determine the x-axis according to the right-hand rule to establish a coordinate system. Calculate the reaction force and the reaction moment on the sealing belt in this coordinate system, and then construct the force balance equation between the reaction force and the residual compressive force suffered by the valve plate pair and the moment balance equation between the reaction moment and the residual overturning moment suffered by the valve plate pair.

4. The rapid calculation method for the leakage of the flow distribution pair with parametric modeling according to claim 2, characterized in that: The specific construction of the force balance equations when the valve plate pair overturns is as follows: ； ； where, F rc and M rc are respectively the reaction force and moment provided by the sealing belt when the cylinder block overturns; solve the above force balance equations, and determine the central oil film thickness h0 and the overturning angle γ.

5. The rapid calculation method for the leakage of the flow distribution pair with parametric modeling according to claim 4, characterized in that: The reaction force F provided by the sealing tape in step (4) rc and the reaction moment M rc The calculation formulas are as shown in the following formula: ； ； wherein, R i and R o are the inner diameter and outer diameter of the sealing strip respectively, h max and h min are the maximum clearance height and minimum clearance height on the sealing strip respectively; wherein, ； ； ； ； In the formula, h0 is the central oil film thickness, γ is the tilting angle of the cylinder block, K is the reaction force stiffness, and δ is the critical threshold of the reaction force. When the clearance height is lower than this value, it indicates that the oil film thickness is very small and the reaction force begins to be generated.

6. The rapid calculation method for the leakage of the port plate pair by parametric modeling according to claim 1, characterized in that: The leakage flow rate of the flow distribution pair in step (5) includes the leakage flow rate Q Ls (R1, R2) of the inner seal and the leakage flow rate Q Ls (R3, R4) of the outer seal, where R1 is the inner diameter of the inner seal band of the flow distribution pair, R2 is the outer diameter of the inner seal band of the flow distribution pair, R3 is the inner diameter of the outer seal band of the flow distribution pair, and R4 is the outer diameter of the outer seal band of the flow distribution pair.

7. The rapid calculation method for the leakage of the flow distribution pair with parametric modeling according to claim 5, characterized in that: In step (5), the leakage flow rate Q of the sealing strip considering the overturning of the cylinder block is considered Ls The calculation formula is as follows: ； In the formula, μ is the viscosity of the oil.

8. A rapid calculation device for the leakage of a flow distribution pair with parametric modeling, comprising a memory and one or more processors, wherein executable code is stored in the memory, characterized in that, When the processor executes the executable code, it realizes a rapid calculation method for the leakage of a parametrically modeled valve plate pair as described in any one of claims 1-7.

9. A computer-readable storage medium having a program stored thereon, characterized in that, When the program is executed by the processor, it realizes a rapid calculation method for the leakage of a parametrically modeled valve plate pair as described in any one of claims 1-7.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it realizes a rapid calculation method for the leakage of a parametrically modeled valve plate pair as described in any one of claims 1-7.

Images