A parameterized modeling flow distribution pair leakage rapid calculation method and device

Abstract

The application discloses a parameterized modeling flow distribution pair leakage quick calculation method and device, first, according to the structural parameters of the plunger pump, the residual compression force and the residual overturning moment of the flow distribution pair are calculated, then the center oil film thickness, the cylinder overturning angle and the cylinder overturning direction of the flow distribution pair are determined through vector composition and force balance equation set solving, finally, the above parameters are brought into the leakage flow calculation formula considering the cylinder overturning, and the leakage flow of the flow distribution pair is obtained. The parameterized modeling flow distribution pair leakage quick calculation method provided by the application retains the applicability of the traditional analytical model under simple working conditions, considers the influence of the cylinder overturning under high pressure and high speed working conditions, expands the application range of the traditional model, so that the quick and accurate estimation of the leakage flow of the flow distribution pair under wide range working conditions is realized, and a quick and accurate mathematical model is provided for the high-precision and high-energy-efficiency control of the EHA system.

Description

Technical Field

[0001] This invention belongs to the field of hydraulic pump technology, specifically relating to a method and device for rapid calculation of leakage in a flow distribution pair using parametric modeling. Background Technology

[0002] Electro-hydraulic actuators (EHAs), as a new generation of servo actuation systems, have broad application prospects in aerospace, engineering machinery, and other fields due to their advantages of high efficiency, compactness, and plug-and-play functionality. EHA systems regulate output flow by adjusting the speed and direction of the motor pump, thereby achieving position servo control of the actuator. This requires the plunger pump to accurately perform feedforward compensation for leakage flow under four-quadrant operating conditions to ensure the dynamic performance and control accuracy of the EHA system. Among all friction pairs, the flow distribution pair has the largest sealing size and is the main source of leakage flow in the plunger pump, accounting for more than 60% of the total leakage. Moreover, under high-pressure and high-speed conditions, cylinder overturning is highly likely, leading to a surge in leakage. Therefore, quickly and accurately calculating the leakage flow of the flow distribution pair is crucial.

[0003] Traditional analytical models assume the distribution pair to be a parallel gap with a constant height, deriving a clear expression for the leakage flow rate. However, this requires experienced engineers to set a reasonable gap height *h* to ensure prediction accuracy. Especially under high-pressure, high-speed conditions, cylinder overturning significantly exacerbates leakage, rendering the simple parallel gap flow assumption inapplicable. Another method for calculating leakage is to establish a numerical model based on friction lubrication theory. This method considers the dynamic behavior of cylinder overturning under high-pressure, high-speed conditions, as well as the complex thermo-elasto-hydraulic lubrication effects at the friction interface. It can more accurately reflect the nonlinear leakage characteristics of the distribution pair under wide-range conditions. However, numerical model calculations often involve multiple iterative cycles and coupled solutions, which are time-consuming and consume significant computational resources, making them difficult to directly apply to EHA control systems.

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

[0005] The purpose of this invention is to address the shortcomings of existing technologies by providing a method and apparatus for rapid calculation of leakage in a flow distribution pair using parametric modeling. This enables rapid and accurate estimation of leakage flow in the flow distribution pair under a wide operating range, which helps the control system to accurately and in real-time compensate for flow loss, thereby achieving high-precision and high-efficiency control of the EHA system.

[0006] The objective of this invention is achieved through the following technical solution: Firstly, this invention provides a method for rapid calculation of leakage in a distribution sub-system using parametric modeling, comprising the following steps:

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

[0008] (2) The design method of residual clamping force is adopted for the distribution pair. The residual clamping force and the residual overturning moment in the horizontal and vertical directions are calculated.

[0009] (3) The cylinder overturning orientation is determined by vector synthesis based on the remaining overturning moment in the horizontal and vertical directions;

[0010] (4) Based on the remaining clamping force and the remaining overturning moment in the horizontal and vertical directions, a set of force balance equations is constructed when the distribution pair overturns, and the center oil film thickness and overturning angle are obtained by solving the equations.

[0011] (5) Substitute the center oil film thickness, cylinder overturning angle and overturning orientation into the leakage flow calculation formula that takes into account cylinder overturning to obtain the leakage flow of the distribution pair.

[0012] Furthermore, the structural parameters of the plunger pump in step (1) include: the number of plungers z, the swashplate angle β, and the plunger diameter d. p plunger mass m p Slipper mass m s Φ, high-voltage zone wrap angle of the distribution plate H , Piston distribution circle radius R p The inner diameter of the inner sealing strip of the distribution pair is R1, the outer diameter of the inner sealing strip of the distribution pair is R2, the inner diameter of the outer sealing strip of the distribution pair is R3, the outer diameter of the outer sealing strip of the distribution pair is R4, and the distance L from the origin to the center of mass of the plunger slipper assembly at the inner dead point is L. m0 The distance L from the origin to the center of the plunger ball head at the inner dead point ss0 and the distance L from the origin to the center of the spline subcenter sc0 .

[0013] Furthermore, taking the overturning orientation in step (3) as the y-axis, pointing from the origin to the direction of the minimum clearance height, and taking the cylinder rotation center as the z-axis, pointing to the swashplate direction, and the x-axis as determined by the right-hand rule, a coordinate system is established. Under this coordinate system, the reaction force and reaction torque on the sealing strip are calculated, and then the force balance equation between the reaction force and the remaining clamping force on the distribution pair, as well as the torque balance equation between the reaction torque and the remaining overturning torque on the distribution pair, are constructed.

[0014] Furthermore, the force balance equations for the overturning of the distribution pair are constructed as follows:

[0015]

[0016]

[0017] In the formula, F rc and M rcThese represent the reaction force and torque provided by the sealing strip when the cylinder block overturns. Solve the above force balance equations, taking the central oil film thickness h0 and the overturning angle γ as fixed.

[0018] Furthermore, the reaction force F provided by the sealing strip in step (4) rc and reaction torque M rc The calculation formula is shown below:

[0019]

[0020]

[0021] In the formula, R i and R o These are the inner and outer diameters of the sealing strip, h, respectively. max and h min These are the maximum and minimum gap heights on the sealing strip, respectively; where,

[0022]

[0023]

[0024]

[0025]

[0026] In the formula, h0 is the center oil film thickness, γ is the cylinder overturning angle, K is the reaction force stiffness, and δ 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.

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

[0028] Furthermore, step (5) takes into account the leakage flow rate Q of the sealing strip due to cylinder overturning. Ls The calculation formula is as follows:

[0029]

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

[0031] Secondly, the present invention also provides a fast calculation device for leakage of a distribution sub-component based on parametric modeling, comprising a memory and one or more processors, wherein the memory stores executable code, and when the processor executes the executable code, it implements the fast calculation method for leakage of a distribution sub-component based on parametric modeling.

[0032] Thirdly, the present invention also provides a computer-readable storage medium having a program stored thereon, characterized in that, when the program is executed by a processor, it implements the aforementioned method for rapid calculation of leakage in a parametrically modeled distribution sub-channel.

[0033] Fourthly, the present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements the aforementioned method for rapid calculation of leakage in a parametrically modeled distribution sub-unit.

[0034] The beneficial effects of this invention are as follows:

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

[0036] 2. This method retains the applicability of the traditional analytical model under simple working conditions and considers the impact of cylinder overturning under high pressure and high speed conditions, thus expanding the application range of the traditional model and enabling it to more accurately reflect the nonlinear leakage characteristics of the distribution pair under wide-range working conditions, with high prediction accuracy. Attached Figure Description

[0037] Figure 1 This is a flowchart of a rapid calculation method for leakage of a distribution sub-unit using parametric modeling, as described in this invention.

[0038] Figure 2 This is a force analysis diagram of the distribution pair;

[0039] Figure 3 This is a schematic diagram showing the tilting position of the cylinder block;

[0040] Figure 4 The clearance height of the distribution pair under different overturning degrees;

[0041] Figure 5 This is a comparison between the method proposed in this invention and the numerical solution model.

[0042] Figure 6 This is a schematic diagram of a rapid calculation device for leakage of a distribution sub-unit using parametric modeling, provided by the present invention.

[0043] Figure label:

[0044] 1. Distributor plate, 2. Cylinder block, 3. Piston, 4. Return plate, 5. Slipper, 6. Swashplate, 7. Center spring, 8. Piston pump drain port, 9. Piston pump suction port. Detailed Implementation

[0045] 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.

[0046] like Figure 1 As shown, the present invention provides a method for rapid calculation of leakage in a distribution sub-system using parametric modeling, comprising the following steps.

[0047] Step 1: Input the plunger pump structural parameters. Taking a 7mL / r plunger pump as an example, the main structural parameters to be input are shown in the table below:

[0048] Table 1 Main structural parameters of plunger pump

[0049]

[0050] Step 2: Perform a stress analysis on the cylinder block and calculate the remaining clamping force F on the distribution pair. rz With the remaining overturning moment M rx and M ry .

[0051] The cylinder block rotates and slides on the surface of the distributor plate, forming the distributor interface. On one hand, the high-pressure oil in the plunger cavity acts on the stepped surface at the bottom of the cylinder bore, causing the cylinder block to be subjected to axial clamping force F. cz The force of the oil film presses tightly against the distributor plate, while the high-speed rotating cylinder block generates an oil film supporting force F between itself and the distributor plate. f This pushes the cylinder block away from the distributor plate, ensuring a close seal between the cylinder block and the distributor plate while also providing good lubrication and load-bearing performance. However, once F cz and F f The lines of action do not coincide, inevitably creating a torque that could overturn the cylinder block. To ensure the normal and stable operation of the plunger pump, the distribution pair needs to employ a design method using residual clamping force. The axial force balance coefficient B of the distribution pair. b for

[0052] (1)

[0053] Where the balance coefficient B b Take 0.9~1, oil film support force F f Slightly less than the clamping force F cz This generates residual clamping force F. rz =F cz -F f This generates a counter-torque to balance the overturning torque, keeping the cylinder in a balanced state.

[0054] like Figure 2 As shown, an OXYZ rectangular coordinate system is established for force analysis. The origin O is located at the center of the distributor plate. The Z-axis coincides with the cylinder block rotation center and points towards the swashplate. The Y-axis is vertically upward, and the X-axis points from the distributor plate's oil inlet to the oil outlet side. The central spring force F... s and plunger chamber pressure F pn The combined forces create an axial clamping force on the cylinder block, which is then transmitted to the distributor plate via the cylinder block, as shown in the following formula:

[0055] (2)

[0056] (3)

[0057] In the formula d p p is the plunger diameter. n Let p be the pressure in the nth plunger chamber. n It can be determined by the following formula:

[0058] (4)

[0059] In the formula p H and p L These represent the pressures at the discharge and suction ports of the plunger pump, respectively. n Let be the angular displacement of the nth plunger chamber. The clamping force generated by the pressure in the discharge-side plunger chamber is greater than that in the suction-side, causing the cylinder block to tilt towards the high-pressure side of the distributor plate. The corresponding hydraulic tilting moment is calculated as follows:

[0060] (5)

[0061] (6)

[0062] In the formula, z represents the number of plungers, and R p Let F be the radius of the plunger distribution circle. The plunger slipper assembly generates a rotational centrifugal force F as the cylinder rotates. cpsn On the other hand, the reciprocating motion within the cylinder bore generates a reciprocating inertial force F. ipsn These two inertial forces act directly on the cylinder block through the plunger assembly:

[0063] (7)

[0064] (8)

[0065] In the formula, m p and m sLet denoted by , and β be the swashplate tilt angle, and ω be the cylinder rotational angular velocity. The plane containing the center of mass of the plunger-slipper assembly is not parallel to the surface of the distributor plate, causing the cylinder block to tilt towards the outer dead center of the distributor plate. The calculation of the corresponding inertial overturning moment is shown in the following formula:

[0066] (9)

[0067] (10)

[0068] In the formula L m0 Let F be the distance from the origin to the center of mass of the plunger slide assembly at the inner dead point. The high-pressure oil and compressed spring within the plunger chamber press the plunger slide assembly against the swashplate, and the reaction force F generated by the swashplate on each plunger slide assembly... ssyn Represented as:

[0069] (11)

[0070] The lateral force acting on the cylinder block results in the following overturning moment:

[0071] (12)

[0072] (13)

[0073] The lateral forces of all plunger and slipper assemblies are counteracted by the spline reaction force F in the Y-axis direction. sc Cancel each other out:

[0074] (14)

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

[0076] (15)

[0077] (16)

[0078] In the formula L sc0 This is the distance from the origin to the center of the spline. In summary, the clamping force F acting on the cylinder block... c and overturning moment M t for:

[0079] (17)

[0080] (18)

[0081] Consider the following mathematical approximation to simplify the above equation:

[0082] (19)

[0083] (20) (twenty one)

[0084] (twenty two)

[0085] And write it in fractional form as follows:

[0086] (twenty three)

[0087] (twenty four)

[0088] (25)

[0089] in,

[0090] (26)

[0091] (27)

[0092] (28)

[0093] The hydrostatic oil film support force of the distribution system on the cylinder block consists of three parts: the hydrostatic support force of the outer sealing strip, the hydrostatic support force of the inner sealing strip, and the hydraulic support force of the high-pressure oil in the drain groove on the cylinder block. On the drain side, the calculation formula for the distribution system's auxiliary oil film pressure distribution is:

[0094] (29)

[0095] In the formula, R1, R2, R3, and R4 are the inner diameter, outer diameter, inner diameter, and outer diameter of the inner sealing strip of the distribution pair, respectively. Integrating the above formula, we obtain the supporting force F of the hydrostatic oil film of the distribution pair on the cylinder block. fH and overturning moment M fHx and M fHy :

[0096] (30)

[0097] (31)

[0098] (32)

[0099] After rearranging the above formula, we can obtain:

[0100] (33)

[0101] (34)

[0102] (35)

[0103] In the formula Φ H For the high-voltage zone wrap angle of the distribution plate, Φ L The low-pressure area wrap angle Φ of the distribution plate L = 2π-Φ H .in,

[0104] (36)

[0105] (37)

[0106] In summary, the remaining clamping force F on the distribution pair rz for:

[0107] (38)

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

[0109] (39)

[0110] (40)

[0111] Step 3: Determine the cylinder overturning orientation ζ through vector synthesis.

[0112] like Figure 3 As shown, the remaining overturning moment M rx and M ry The combined effect determines the overturning orientation of the cylinder block:

[0113] (41)

[0114] Step 4: Derive the analytical expression of the reaction force of the flow distribution pair sealing strip, construct a set of force balance equations, and solve for the central oil film thickness h0 and the overturning angle γ.

[0115] The unbalanced torque causes the cylinder block to overturn relative to the distributor plate, resulting in solid contact where the oil film on the sealing strip is very thin. This provides additional reaction force to maintain the force balance of the cylinder block. Assume the reaction force p of the distributor clearance height and the solid contact on the sealing strip... r The relationship between them is linear, that is:

[0116] (42)

[0117] In the formula, h is the height of the distribution pair clearance, 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 a reaction force begins to be generated. The clearance height at any point (r, θ) of the sealing strip is expressed as:

[0118] (43)

[0119] In the formula, h0 is the thickness of the central oil film, and γ is the cylinder overturning angle. The coordinate system oxyz is defined to describe the reaction force generated on the sealing strip, such as... Figure 3 As shown, the z-axis coincides with the Z-axis, the y-axis coincides with the overturning azimuth angle and 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 circumferential sealing strip can be rewritten as:

[0120] (44)

[0121] In the formula θ r denoted as angular displacement in the oxyz coordinate system.

[0122] like Figure 4 As shown, based on the severity of cylinder rollover, it can be divided into mild rollover, severe rollover, and complete rollover. For one of the sealing strips of the distribution pair, when the cylinder rollover is mild, the rollover angle is small, and the gap height on the sealing strip is below the height threshold δ, that is:

[0123] (45)

[0124] In the formula R m Let F be the average radius of the sealing strip. At this point, the reaction force F on the sealing strip... rc and reaction torque M rc It can be calculated as follows:

[0125] (46)

[0126] (47)

[0127] In the formula R i and R o These are the inner and outer diameters of the sealing strip, respectively. When the cylinder body tilts severely, the tilt angle is large, and the flow distribution pair exhibits a significant wedge-shaped gap. Part of the gap height will exceed the critical threshold, i.e.:

[0128] (48)

[0129] At this moment, the reaction force F on the sealing strip rc and reaction torque M rc It can be calculated as follows:

[0130] (49)

[0131] (50)

[0132] In the formula θ r1 and θ r2 The angular displacement is calculated when the gap height equals the critical threshold. The formula is:

[0133] (51)

[0134] When the cylinder block completely tilts over, it will detach completely from the distributor plate, and the clearance height will be above the critical threshold, i.e.:

[0135] (52)

[0136] At this moment, the reaction force F on the sealing strip rc and reaction torque M rc All are 0. Therefore, the reaction force and torque of the sealing strip under different cylinder tilting degrees are:

[0137] (53)

[0138] (54)

[0139] In the formula, h max h is the maximum gap height on the sealing strip. max =h0+R m γ, h min The minimum gap height h on the sealing strip max =h0-R m γ, where,

[0140] (55)

[0141] (56)

[0142] (57)

[0143] (58)

[0144] In summary, the force balance equations for the flow distribution pair include force balance and torque balance, as shown in the following equation:

[0145] (59)

[0146] (60)

[0147] The thickness of the central oil film h0 and the overturning angle γ can be obtained by solving the force balance equations.

[0148] Step 5: Derive the formula for calculating the leakage flow rate of the distribution pair considering cylinder block overturning, and calculate the leakage flow rate of the distribution pair.

[0149] Considering that the leakage flow of the distribution pair is laminar flow in the gap, and assuming it to be steady, continuous, incompressible, and with negligible body forces, the leakage calculation formula for the distribution pair along the circumferential direction of the sealing strip, based on fluid mechanics, is as follows:

[0150] (61)

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

[0152] (62)

[0153] Summarized as follows:

[0154] (63)

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

[0156] (64)

[0157] It is worth mentioning that the formula for calculating the leakage flow rate of the distribution pair considering cylinder overturning, derived in this invention, is a more extensive development of the traditional classical analytical model. When the overturning angle is 0, i.e., γ=0, the above leakage flow rate calculation formula will degenerate into the traditional classical analytical model:

[0158] (65)

[0159] To further illustrate the advantages of this invention, for the examples given in this invention, the leakage flow rate of the distribution pair under different pressures and rotational speeds was calculated using both the numerical model and the analytical model proposed in this invention. The calculation results are as follows: Figure 5 As shown in the figure, the analytical model's calculation results are in good agreement with the numerical model, accurately reflecting the nonlinear leakage characteristics of the distribution pair under a wide range of operating conditions. More importantly, the numerical model requires approximately 2.5 hours of computation to solve for the leakage flow rate under one operating condition. For the 12 speed conditions (1000 r / min to 12000 r / min) and 5 pressure conditions (7 MPa to 35 MPa) given in the figure, the total computation time reaches 150 hours. In contrast, the analytical method proposed in this invention solves for the leakage flow rate under all operating conditions in just 0.6 seconds, demonstrating extremely high computational efficiency.

[0160] Corresponding to the aforementioned embodiment of a method for rapid calculation of leakage in a parametrically modeled distribution pair, the present invention also provides an embodiment of a device for rapid calculation of leakage in a parametrically modeled distribution pair.

[0161] See Figure 6 The present invention provides a fast calculation device for leakage of a distribution sub-type based on parametric modeling, comprising a memory and one or more processors. The memory stores executable code, and when the processor executes the executable code, it is used to implement a fast calculation method for leakage of a distribution sub-type based on parametric modeling as described in the above embodiment.

[0162] The embodiment of the parametric modeling-based fast calculation device for distribution sub-leakage provided by this invention can be applied to any device with data processing capabilities, such as a computer. The device embodiment can be implemented in software, hardware, or a combination of both. Taking software implementation as an example, as a logical device, it is formed by the processor of any data processing device loading the corresponding computer program instructions from non-volatile memory into memory for execution. From a hardware perspective, such as... Figure 6 The diagram shown is a hardware structure diagram of any device with data processing capabilities, where the parametric modeling flow distribution sub-leakage fast calculation device provided by this invention is located. (Except for...) Figure 6 In addition to the processor, memory, network interface, and non-volatile memory shown, any data processing device in the embodiment may also include other hardware depending on the actual function of the data processing device, which will not be described in detail here.

[0163] The specific implementation process of the functions and roles of each unit in the above device can be found in the implementation process of the corresponding steps in the above method, and will not be repeated here.

[0164] For the device embodiments, since they basically correspond to the method embodiments, the relevant parts can be referred to in the description of the method embodiments. The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of the present invention according to actual needs. Those skilled in the art can understand and implement this without creative effort.

[0165] This invention also provides a computer-readable storage medium storing a program thereon, which, when executed by a processor, implements a method for rapid calculation of leakage in a distribution sub-unit using parametric modeling as described in the above embodiments.

[0166] The computer-readable storage medium can be an internal storage unit of any data processing device described in any of the foregoing embodiments, such as a hard disk or memory. The computer-readable storage medium can also be an external storage device of any data processing device, such as a plug-in hard disk, smart media card (SMC), SD card, flash card, etc., equipped on the device. Furthermore, the computer-readable storage medium can include both internal storage units and external storage devices of any data processing device. The computer-readable storage medium is used to store the computer program and other programs and data required by the data processing device, and can also be used to temporarily store data that has been output or will be output.

[0167] The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements the aforementioned method for rapid calculation of leakage in a parametrically modeled distribution sub-unit.

[0168] The embodiments described above are merely preferred embodiments of the present invention and are not intended to limit the scope of the present invention. Various modifications and improvements made by those skilled in the art to the technical solutions of the present invention without departing from the spirit of the present invention should fall within the protection scope defined by the claims of the present invention.

Claims

1. A method for rapid calculation of leakage in a distribution pair using parametric modeling, characterized in that, Includes the following steps: (1) Obtain the structural parameters of the plunger pump; (2) The design method of residual clamping force is adopted for the distribution pair. The residual clamping force and the residual overturning moment in the horizontal and vertical directions are calculated. (3) The cylinder overturning orientation is determined by vector synthesis based on the remaining overturning moment in the horizontal and vertical directions; (4) Based on the remaining clamping force and the remaining overturning moment in the horizontal and vertical directions, a set of force balance equations is constructed when the distribution pair overturns, and the center oil film thickness and overturning angle are obtained by solving the equations. (5) Substitute the center oil film thickness, cylinder overturning angle and overturning orientation into the leakage flow calculation formula that takes into account cylinder overturning to obtain the leakage flow of the distribution pair.

2. The method for rapid calculation of leakage in a distribution pair using 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 swashplate angle β, and the plunger diameter d. p plunger mass m p Slipper mass m s Φ, high-voltage zone wrap angle of the distribution plate H , Piston distribution circle radius R p The inner diameter of the inner sealing strip of the distribution pair is R1, the outer diameter of the inner sealing strip of the distribution pair is R2, the inner diameter of the outer sealing strip of the distribution pair is R3, the outer diameter of the outer sealing strip of the distribution pair is R4, and the distance L from the origin to the center of mass of the plunger slipper assembly at the inner dead point is L. m0 The distance L from the origin to the center of the plunger ball head at the inner dead point ss0 and the distance L from the origin to the center of the spline subcenter sc0 .

3. The method for rapid calculation of leakage in a distribution pair using parametric modeling according to claim 1, characterized in that: Using the overturning orientation in step (3) as the y-axis, pointing from the origin to the minimum clearance height direction, and the cylinder rotation center z-axis pointing to the swashplate direction, the x-axis is determined according to the right-hand rule to establish a coordinate system. Under this coordinate system, the reaction force and reaction torque on the sealing strip are calculated, and then the force balance equation between the reaction force and the remaining clamping force on the distribution pair and the torque balance equation between the reaction torque and the remaining overturning torque on the distribution pair are constructed.

4. The method for rapid calculation of leakage in a distribution pair using parametric modeling according to claim 2, characterized in that: The force balance equations for the overturning of the distribution pair are constructed as follows: In the formula, F rc and M rc Let F be the reaction force and torque provided by the sealing strip when the cylinder block overturns; solve the above force balance equations, taking the center oil film thickness h0 and the overturning angle γ as constants. rz M is the remaining clamping force on the distribution pair. rx and M ry The remaining overturning moments are in the x and y directions.

5. The method for rapid calculation of leakage in a distribution pair using parametric modeling according to claim 4, characterized in that: The reaction force F provided by the sealing strip in step (4) rc and reaction torque M rc The calculation formula is shown below: In the formula, R i and R o These are the inner and outer diameters of the sealing strip, h, respectively. max and h min These are the maximum and minimum gap heights on the sealing strip, respectively; where, In the formula, h0 is the center oil film thickness, γ is the cylinder overturning angle, K is the reaction force stiffness, and δ 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.

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

7. The method for rapid calculation of leakage in a distribution pair using parametric modeling according to claim 5, characterized in that: Step (5) takes into account the leakage flow rate Q of the sealing strip due to cylinder overturning. Ls The calculation formula is as follows: In the formula, μ is the viscosity of the oil, and p H ζ represents the discharge port pressure of the plunger pump, and ζ represents the cylinder overturning position.

8. A fast calculation device for leakage of a distribution sub-unit using parametric modeling, comprising a memory and one or more processors, wherein the memory stores executable code, characterized in that, When the processor executes the executable code, it implements a method for rapid calculation of leakage in a distribution sub-system using parametric modeling 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 implements a method for rapid calculation of leakage of distribution sub-suppliers using parametric modeling 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 implements a method for rapid calculation of leakage in a distribution sub-system using parametric modeling as described in any one of claims 1-7.

Images