1. Introduction
Axial piston pumps, integral as the dynamic components within hydraulic systems, significantly influence the overall performance and longevity of these systems. Presently, axial piston pumps are categorized into two distinct types: swash-plate axial piston pumps and inclined-axis axial piston pumps. Both variants operate based on the principle of piston reciprocation within the cylinder block, which varies the sealed working volume to facilitate the processes of oil intake and discharge. However, they are distinguished by their internal configurations; a swash-plate axial piston pump induces piston reciprocation through the inclination of the swash plate, while an inclined-axis axial piston pump achieves this motion by setting a specific rotational angle for the cylinder block.
Volumetric efficiency, a pivotal parameter indicative of an axial piston pump’s performance, mirrors the internal fluid dynamics and potential energy losses. It serves as a critical reference for pump selection and application. The efficiency’s scale is predominantly affected by the pump’s geometric design, operational conditions, and the characteristics of the fluid in use.
South Korean scholars, including Beak, have utilized geometric approaches to investigate the performance characteristics of inclined-axis pumps, elucidating the kinematic principles of conical piston-driven cylinder motion and delineating the operational scope of conical piston drives. Experimental validations have confirmed the significant impact of the piston inclination angle and the phase angle between the piston assembly and the cylinder on the performance of inclined-axis pumps [
1,
2]. Researchers such as Hong have concentrated on the frictional losses within inclined-axis pumps through mathematical theory, with the goal of pinpointing the primary factors affecting output torque. Experimental assessments of friction in mechanical components, including pistons, spherical heads, bearings, and distributor disks, were conducted. These empirical data were integrated into theoretical models of piston kinematics and pump friction, enabling the calculation of the frictional torque exerted by the piston. A comparison with the actual input torque substantiated the accuracy of their model. Their findings suggest that the viscous frictional forces from the distributor disk and bearings are the predominant sources of frictional losses, whereas the frictional force and torque of the piston are negligible in friction loss studies [
3].
Kim and associates explored the influence of the distributor disk dimensions on the pressure pulsation within the piston chamber of inclined-axis pumps, as well as on the overall output pressure pulsation and torque pulsation of the pump. Utilizing Amesim for model establishment, they juxtaposed simulation analyses with experimental research outcomes. Their findings indicated that the dimensions of the distributor disk have a substantial effect on the pressure pulsation within the piston chamber, with the overall output pressure pulsation being correlated with the basal pressure of the piston chamber yet exerting a minor impact on the pump’s torque [
4,
5,
6,
7]. Furthermore, the size of the distributor disk significantly influences the pump housing’s vibration [
8].
Professor Jung, J.Y. from Chonbuk National University in South Korea conducted profound theoretical analyses of the principle of rotation in conical piston-driven cylinders, followed by experimental validation. Addressing the angular discrepancy in inclined-axis axial piston pumps, he proposed a design scheme involving transmission shaft inconsistencies [
9,
10]. Professor Jung, J.Y. also delved into the lubrication properties of a conical piston ring, observing that the thickness of the oil film between the piston ring and the cylinder is increased with an increasing rotational speed but remains largely unaffected by the pump’s outlet pressure. Conversely, the frictional force between the piston ring and the cylinder escalates with both rotational speed and pump outlet pressure [
1,
7]. Simulation studies and experimental validations on the flow pulsation characteristics of inclined-axis axial piston pumps were also undertaken by Professor Jung J.Y. [
6].
Roccatello, A. and Nervegna, N. of the Polytechnic University of Milan, Italy, engineered an innovative cylinder block driving mechanism that employs gear engagement between the main shaft flange and the cylinder block. This advancement has ameliorated the stress conditions on the pistons and resolved the issue of angular discrepancy commonly found in conventional inclined-axis pumps [
11]. Furthermore, Professor Abuhaiba, M. from the Islamic University of Gaza and Professor Olson, W. W. from the University of Toledo conducted an exhaustive kinematic analysis of an articulated inclined-axis pump, which, like the aforementioned design, is devoid of angular discrepancy. Their comprehensive study provides substantial guidance for the design of articulated inclined-axis pumps [
12].
Scholars from Zhejiang University, including Zhang Junhui, investigated the effects of the clearance between the piston and the cylinder on the efficiency of pumps. They formulated a mathematical expression for the leakage phenomenon considering the centrifugal and reciprocating inertial forces acting on the piston, which holds significant implications for the structural design of inclined-axis pumps in the future [
13]. Lin Jing et al. utilized mathematical spatial coordinate transformation to analyze the motion mechanisms of inclined-axis pumps, identifying three potential interference scenarios. They established criteria for detecting interference and proposed methods for resolution, thus providing direction for the design and manufacturing processes for inclined-axis axial piston pumps [
14]. Quan Lingxiao constructed a hydraulic and mechanical integrated model of an inclined-axis pump using Amesim and ADAMS, conducting simulations to assess the influence of the valve core spring and damping on the pump’s variable control performance and output flow characteristics [
15]. Shu Wang optimized the design of axial piston pumps using the concept of pressure carryover, formulated the mathematical relationship between the valve-plate geometry and volumetric efficiency, and introduced new terminologies, such as Discrepancy of Pressure Carryover (DPC) and Carryover Cross-porting (CoCp). This study offers guidance for the design of high-efficiency valve plates, with the results validated through simulations and experimental testing [
16]. N. D. Maring compared the valve-plate designs in axial-piston pumps, demonstrating that a trapped-volume design without slots has improved efficiency compared to traditional designs. The study provides theoretical support for trapped-volume optimization in valve-plate design [
17].
In summary, scholarly research on inclined-axis axial piston pumps has predominantly concentrated on aspects such as the vibration characteristics of the inclined-axis axial piston pump, the mechanism underlying the formation of the oil film, the influence of the leakage between frictional pairs on the pump’s volumetric efficiency, and the variable characteristics of the inclined-axis axial piston pump. However, the kinematic properties of hingeless inclined-axis axial piston pumps have not been articulated. In the hingeless design of an axial piston pump, the axis of the piston is inclined at an angle τ to the axis of the cylinder bore. This angular relationship not only induces a phase angle difference ϕ between the cylinder block and the transmission shaft but also introduces a layer of complexity to the piston’s motion.
3. Simulation
The finite element analysis presented in this paper is performed utilizing the Simerics+4.6.0 (Simerics Inc., Seattle, WA, USA), which applies the finite volume method to the simulation of an inclined-axis axial piston pump. The equations associated with the finite volume method are presented as follows:
The equation of mass conservation [
18]:
Equation for the conservation of momentum [
18]:
Equation for the conservation of energy [
18]:
The turbulence model utilized for the inclined-axis piston pump is the standard k-ε model, which exhibits enhanced stability, computational precision, and economic suitability for the simulation of inclined-axis axial piston pumps.
Boundary Condition
In order to achieve simulation outcomes that closely reflect actual working conditions, the setting of appropriate boundary conditions is of significant importance. Considering the impact of pressure reflection, the simulation of the inclined-axis piston pump presented in this paper incorporates the extension of the lengths of both the inlet and outlet pipelines.
A throttling orifice is employed at the terminus to simulate authentic operating scenarios. The dimensions of the throttling orifice are ascertained based on the subsequent equation:
The fidelity of numerical simulations of fluid dynamics is intricately linked to the quantity of the model’s grid cells. In general terms, an increased number of grid cells typically results in enhanced simulation precision. Nonetheless, the constraints of computational resources necessitate a balanced approach to determining the optimal number of grid cells. Within the preprocessing phase, the meshing functionality inherent to the software is utilized to combine various body meshes, surface meshes, and other controlled dimensions, thereby generating a multitude of grid models. Following this, preliminary simulations are executed for each model, and a comparative analysis is undertaken to ascertain the most appropriate grid model.
Table 1 provides the control parameter information for the mesh of the fluid domain.
In order to ascertain the independence of the fluid domain from the grid, the pressure within the piston cavity is chosen for surveillance, as illustrated in
Figure 2. The pressure at the designated monitoring point and the resultant output flow rate are utilized as the criteria for assessment, with the outcomes of the simulation depicted in
Figure 3 and
Figure 4.
Referring to
Figure 2, it is evident that prior to the 0.096 s mark in the operation of the piston pump, the pressures at an identical monitoring point across the models exhibit a fundamentally similar trajectory of variation. Upon the piston’s transit into the high-pressure region, there occurs a significant oscillation in the hydraulic pressure within the piston cavity, with a positive overshoot in the hydraulic pressure observed at the 0.096 s mark. The simulated peak hydraulic pressure within the piston cavity during its operation in the high-pressure region, as ascertained by Model A, is identified at 16.99 MPa. Model B’s simulation predicts a peak hydraulic pressure of 14.2 MPa in the high-pressure region. Model C’s simulation delineates a peak hydraulic pressure of 14.6 MPa in the high-pressure region, while Model D’s simulation indicates a peak hydraulic pressure of 14 MPa within the same region.
The inclined-axis piston pump, as a result of its operational mechanism, encounters cyclical pressure fluctuations. As depicted in
Figure 3, the trend in these fluctuations is uniform across the models examined. Model A records a peak pressure of 11.62 MPa and a minimum pressure of 6.4 MPa; it should be noted that the repetition of values for Model A likely indicates a typographical error in the source text. Model B exhibits a peak pressure of 11.35 MPa and a minimum pressure of 6.4 MPa. Model C’s peak pressure is documented at 11.8 MPa, with a minimum pressure of 6.8 MPa. Meanwhile, Model D has a peak pressure of 11.3 MPa and a minimum pressure of 6.43 MPa.
The inclined-axis piston pump, as a consequence of its operational principles, is subject to periodic flow pulsations. As depicted in
Figure 4, the trend in these pulsations is uniform across the models analyzed. Model A achieves a peak flow rate of 239 L/min, with a nadir of 171.3 L/min; Model B registers a peak flow rate of 239.4 L/min and a minimum of 172 L/min; Model C exhibits a peak flow rate of 239.13 L/min and a minimum of 175.72 L/min; whereas Model D demonstrates a peak flow rate of 239.7 L/min, with its minimum flow rate also recorded at 175.72 L/min.
Upon comparative assessment of the pressures within the piston cavities, the output pressures, and the output flow rates among the various models, Model C was determined to be the most suitable and was therefore selected to proceed with in the subsequent stages of the work. The methodology for mesh division is detailed, and the specifics can be observed in
Figure 5.
The piston pump parameters are as shown in
Table 2. The parameters of the oil medium and some working conditions are shown in
Table 3.
The rest of the simulation parameters are shown in
Table 3.
The distribution mechanism of this pump is equipped with a spherical oil film, where the termini of the oil film are established as pressure output interfaces. Damping holes are situated at the piston’s extremities to emulate the leakage that transpires between the piston and the flange during its operational phase, with the resultant pressures equated to the discharge oil pressure at 0.2 MPa. The simulations are executed under load pressures of 10 MPa and 20 MPa, respectively. Simulated pressure cloud diagrams of the plunger are shown in
Figure 6 and
Figure 7.
In summary, the correction factor for the inclined-axis pump, when operated at the maximum cylinder swing angle of 25°, has been determined through calculation as follows:
Upon securing the simulation outcomes, they are subsequently amplified by the correction factor to yield the refined results. The simulation results are shown in
Table 4 and
Table 5. The revised simulation results are shown in
Table 6 and
Table 7.
These observations indicate that at an identical rotational speed, there is a negative correlation between the load pressure and the output flow rate; as the load pressure escalates, the output flow rate diminishes. This phenomenon can be attributed to the increased leakage at the distributor and the piston ball joint, which becomes more pronounced with higher pressures, consequently leading to a decrease in the output flow rate at a constant rotational speed.
The output flow rate of the pump is ascertained through experimental documentation, utilizing the following depicted experimental apparatus for the measurement of the pump’s output flow rate.
4. Verification
The output flow rate of the pump was measured using a test recording of the pump’s output flow rate with the experimental setup shown in
Figure 8 and
Figure 9.
A stable oil supply to the test pump was ensured by employing two makeup pumps. Thereafter, the load pressure of the pump was measured using a pressure gauge. Concurrently, the flow rate of the test pump was monitored with a flowmeter. Additionally, a proportional relief valve was utilized to regulate the load pressure of the pump.
In order to compare the accuracy of the simulation outcomes with the actual performance, experimental results were procured under conditions identical to those of the simulation. The parameters set in the experiment are shown in
Table 8.
A comparison between the experimental and simulation results is presented, as shown in
Figure 10 and
Figure 11.
These observations indicate that at an identical rotational speed, the output flow rate of the pump diminishes with an increase in load pressure.
Examination of the graph reveals that for a constant load pressure across different rotational speeds, there exists a variance between the output flow rates derived from the simulation and those garnered from the experimental data. These discrepancies amplify at escalating rotational speeds, a consequence of the additional displacement ∆s, which stems from the angular variation between the cylinder block and the piston—an element not factored into the simulation. This oversight results in an error in the calculated displacement. The correlation between the output flow rate
q and the displacement
V is articulated as follows:
The discrepancies in the flow rates for the various operating conditions mentioned above are detailed in
Table 9 and
Table 10.
The discrepancies in the output flow rates illustrated in the preceding figure, once translated into the variances in displacement, are articulated in
Table 11 and
Table 12.
Upon reviewing
Table 10, it becomes apparent that the displacement values ascertained via simulation are not congruent with those obtained from the experimental measurements. When the rotational speed is held constant, the variance between the displacements from simulation and experimentation enlarges with an elevation in the load pressure. This is primarily due to the increase in load, which results in greater leakage through the distributor pair oil film and the piston oil film during actual testing. The finite element simulations, however, do not account for the leakage from the piston pair oil film, thereby assuming a leakage value of zero for this component. As a result, with an increase in load pressure, the disparity between the experimental measurements and the simulation results becomes more significant.
Upon reviewing
Figure 12 and
Figure 13, it is evident that the discrepancies between the outcomes post-adjustment with the correction factor and the results obtained from experimentation are significantly diminished. For each operation scenario, a comparative analysis of the error between the finite element simulation results and the experimental data, as well as a comparative analysis of the error between the adjusted simulation results and the experimental results, is shown in
Figure 14 and
Figure 15.
Inspection of
Figure 14 and
Figure 15 reveals that across all operating conditions, the outcomes post-adjustment with the correction factor
k are more congruent with the results garnered from experimentation. At a load pressure of 10 MPa, the mean discrepancy between the finite element simulation results and the experimental data is 8.95%, while the mean discrepancy between the adjusted simulation results and the experimental data is reduced to 0.23%. Similarly, at a load pressure of 20 MPa, the average error between the unmodified simulation results and the experimental results is 9.15%, and this is minimized to an average error of 0.36% with the application of the correction factor
k. This adjustment notably diminishes the error rate relative to the experimental outcomes, thereby augmenting the precision of the simulation results.
To encapsulate, in the context of finite element simulation pertaining to an inclined-axis axial piston pump, simulation outcomes that solely account for the forced displacement s engendered by the transmission shaft flange are not adequately precise. It is imperative to incorporate the additional displacement Δs, which arises from the fluctuation in the angularity τ between the piston axis and the cylinder bore axis, to attain accurate results. The application of the correction factor k serves to align the finite element simulation outcomes more closely with the results obtained from experimentation.