# Research on the Effects of Hydropneumatic Parameters on Tracked Vehicle Ride Safety Based on Cosimulation.

1. Introduction

Hydropneumatics is a technology typically used in vehicle suspensions. A hydropneumatic suspension has good properties such as nonlinear stiffness and damping, high power density, convenient tuning, and vertical position locking. Thus, it is widely used in tracked vehicles and improves ride comfort. In practical applications, common failures are the track separating from the road wheel or the sprocket. These failures may result in the vehicle losing control and bring very serious implications in terms of vehicle ride safety.

Previous studies have done some dynamic simulation analysis on the issue of track separating from the road wheel. It has been shown that improper tuning of the suspension parameters is the main reason for such failures [1, 2]. For a hydropneumatic suspension system, the primary parameters are the accumulator charge pressure and the damping. If these parameters are not tuned properly, or if they vary because of external gas leakage and internal oil leakage, then ride safety can be severely compromised. Currently, most research investigations focus on vehicle ride comfort, with little attention given to ride safety. Some dynamic analysis research has been conducted , but this work was based on a simplified mathematical model of the suspension and does not consider the dynamic influence of the track. Other research, which is usually conducted in a single simulation environment, is lacking systematic and quantitative analysis of how the hydropneumatic parameters affect the ride safety .

Determining the dynamic behavior of the track and its interactions with the hydropneumatic suspension system is very difficult to achieve mathematically. This is due to the complex nature of the nonlinear multibody system. Multi-body dynamics (MBD) simulation has long been recognized as an excellent method to predict the dynamic response of vehicles [5-9]. One of the key advantages of this approach is that it allows nonlinear elements such as the accumulator and the damping of the hydropneumatic suspension to be modeled readily.

The introduction of hydropneumatic suspension components that require extensive use of mathematical functions brings challenges for existing software that use the MBD approach and may restrict its use. MBD cosimulation gives a suitable framework for coupling software tools which specialize in different fields of mechanics without sacrificing overall accuracy, particularly if based on different mathematical methods .

The approach adopted here integrates the MBD software (RecurDynV8R3 ) used for modeling complex configurations of the tracked vehicle and the road model, with LMS Imagine.Lab AMESim14.0  used for the development of hydropneumatic systems. RecurDyn is a computer aided engineering tool for MBD analysis. It is capable of simulating vibrations, motions, stress, and rigid-flexible coupling analysis. Its Track (HM) module is excellent for modeling high mobility tracked vehicles. LMS Imagine.Lab AMESim is a software package used for modeling complex physical systems containing mechanical, hydraulic, electronic, thermal, and control components.

On the basis of these conditions, the detailed mechanical configurations and hydropneumatic elements can be modeled more readily using a cosimulation approach. The enormous and complicated calculations relating to the dynamic behavior of flexible objects such as tracks can be fulfilled. This leads to more accurate relationships between the parameter variations and the output variables. In this paper, the hydraulic-mechanical cosimulation on a virtual G road is conducted. Insights into how the hydropneumatic parameters affect the vehicle ride safety are achieved. The results can be readily used in the design, optimal control, and failure detection of a hydropneumatic suspension system.

This paper is organized as follows. A detailed description of the ride safety of tracked vehicles and the numerical model development is presented in Section 2; Section 3 shows the validation process and validity of the constructed model; Section 4 reports the cosimulation results and the analysis; finally, conclusions are presented in Section 5.

2. Model Development

2.1. Ride Safety and Its Evaluation Parameters. For a high mobility tracked vehicle, the suspension design and tuning should improve the ride comfort on the condition that the ride safety is guaranteed. From the point of view of surviving in the battlefield, the ride safety is generally more important than the ride comfort. The ride safety of wheeled vehicles is usually evaluated by the dynamic load of the wheel. In contrast to wheeled vehicles, the special structure of the track and its interactions with the running system brings new problems to the ride safety. The factors that affect ride safety and the evaluation parameters with their criteria are proposed as follows.

(1) Track separating from the road wheel or the sprocket: the failures of the track separating from the road wheel or the sprocket are usually caused by continuous running after the relative position of the track and the road wheel (sprocket) is skewed. The main causes may be (a) the gap between the track and the roadwheel which is too large (often over the teeth height) usually caused by variations in the hydropneumatic suspension parameters, (b) the track excessive jumping caused by insufficient track tensioning force , and (c) large lateral forces acting on the road wheel when the vehicle turns or runs on a side slope.

By comparing the above three causes, it is identified that the first two are more critical than the third. Thus, the track-wheel gap and the quantity of jumping are taken as evaluation parameters for determining whether track separation from the road wheel or sprocket has occurred.

(2) End-stop impact of suspension: at the end-stop of the suspension the arm hits into the bumper . This can cause enormous shock to the crew and may damage the arm and the hull. In a hydropneumatic suspension, the suspension stroke variation is categorized as a Gaussian random process. This means that if the RMS (root of mean square) value of the stroke of the actuation cylinder is higher than one-third of the designed stroke length, collision with the bumper is a possible event [15,16].

As the values of the proposed evaluation parameters listed above belong to a Gaussian random distribution, a statistical metric is more meaningful for evaluating the ride safety. Table 1 shows the ultimate evaluation parameters and their criteria used in this paper.

2.2. Vehicle and Road Modeling. As mentioned previously, a conventional calculation method based on the simplified mathematical model (high degree of linearization of the system stiffness and damping) is not suitable for a complicated nonlinear multibody system . The approach cannot effectively solve the interactions between the track and ground and the track and road wheels and the effect of track tightening force on the track. In this paper, the approach of MBD cosimulation is used to establish the vehicle dynamics model and the ground model.

(1) Vehicle model: the vehicle components were first modeled using a CAD toolkit. Then, the model is exported to the RecurDyn, where the mass, joint constraints, and the motions are assigned. The model, which is shown in Figure 1, consists of the hull, the suspension system, and the track assemblies . Each track assembly is composed of a sprocket, road wheels, an idler, rollers, and a track. According to the actual structure of the vehicle, the running system has 6 pairs of road wheels and arms with actuation cylinders; 3 pairs of supporting rollers; 1 pair of idlers; 1 pair of sprockets; and 101 track plates. The parameters values of these parts are in accordance with those of a real vehicle. The total mass of the structure is about 28.5 tons.

(2) Mechanical model of hydropneumatic suspension: the connection sketch of a single hydropneumatic suspension unit with the hull is shown in Figure 2. The actuation cylinder is connected to the hull by a revolution joint. The arm, which is fixed with a crank pin, is assembled in the hole of the bottom plate through a supporting bearing. The crank and the crankpin are fixed and turn simultaneously with the arm.

Based on the configuration of Figure 2, the equations of the vehicle's equivalent stiffness and damping are shown as follows:

k = [i.sup.2][k.sub.s] + [F.sub.s] di/df,

c = [i.sup.2][c.sub.d], (1)

where i is the suspension transmission ratio, [F.sub.s] is the gas spring force, [k.sub.s] is the gas spring stiffness, di/df is the derivative of i with the vertical displacement of the road wheel f, and [c.sub.d] is the damping coefficient of the hydropneumatic suspension. Their mathematical expressions are shown as follows:

i = [r.sub.s]/R, sin [gamma] / cos [delta] cos [alpha], (2)

[F.sub.s] = (p - [p.sub.a]) [A.sub.s], (3)

[k.sub.s] = n ([F.sub.s] + [p.sub.a][A.sub.s]) / w, (4)

[mathematical expression not reproducible], (5)

[c.sub.d] = d[F.sub.c]([??]) / d[??]. (6)

In (6), [F.sub.c] ([??]) is the damping force generated by the hydraulic restrictor of the suspension. It can be expressed as follows:

[F.sub.c] ([??]) = [DELTA]p ([??]) [A.sub.s]. (7)

For the above expressions, the definition of the other variables and their initial values are given in Table 2.

(3) Road model: the road model is derived based on the technique of solid modeling. A road roughness coefficient is an evaluation index which determines the road class. The following formula is used as the fitting expression of the power spectral density function of the international standard road:

[G.sub.q] (N) = [G.sub.q] ([N.sub.0]) [(N/[N.sub.0]).sup.-[omega]], (8)

where N is the spatial frequency ([m.sup.-1]), [N.sub.0] is the reference spatial frequency (0.1 [m.sup.-1]), [G.sub.q] ([N.sub.0]) is the road power spectrum value at the reference spatial frequency, also called the road roughness coefficient ([m.sup.2]/[m.sup.-1] = [m.sup.3]), and [omega] is the frequency index which is the diagonal slope in a double logarithmic coordinate axis [20, 21].

In this paper, G class road is used to conduct the virtual simulation. The road file is first written in Matlab and introduced to the RecurDyn as shown in Figure 1.

2.3. Hydraulic Model of the Hydropneumatic Suspension. The hydraulic system of the hydropneumatic suspension and the interface module are modeled using AMESim, as shown in Figure 3. The model consists of 12 hydropneumatic suspension units, one pair of tension cylinders, and one interface module. Each unit is composed of an actuation cylinder, a damping tube, an accumulator, and one pair of relief valves. The hydraulic port of the actuation cylinder is connected with the accumulator and the mechanical port is linked with the interface module. The disturbance from the road causes the actuation cylinder to extend and retract. Thus, the oil flows into or out of the accumulator through the damping tube. The relief valves are set to be closed in order to study the effect of the damping tube. The tensioning cylinders are used for tightening the two tracks.

(1) Mathematical model of the actuation cylinder: the output force of the actuation cylinder contains the gas spring force [F.sub.s], the damping force generated by oil flow [F.sub.c], and the friction force [F.sub.f] [22,23]. The total output force can be written as

F = [F.sub.s] + [F.sub.c] + [F.sub.f], (9)

where [F.sub.s], [F.sub.c] are defined in Section 2.1 and [F.sub.f] is expressed by

[mathematical expression not reproducible]. (10)

The static friction force and the dynamic friction force are calculated from the following, respectively.

[F.sub.sta] = [F.sub.f piston] + [F.sub.f rod], (11)

[mathematical expression not reproducible], (12)

where [F.sub.f piston] and [F.sub.f rod] are the static friction forces of the piston and the rod, [v.sub.M] is the relative velocity between the piston and the cylinder barrel, [v.sub.tr1] and [v.sub.tr2] are the slip velocities for fully developed dynamic friction of the piston and the rod, and [[mu].sub.D1] and [[mu].sub.D2] are the dynamic friction coefficients of the piston and the rod.

(2) Mathematical model of the damping tube: a long and thin damping tube is arranged between the actuation cylinder and the accumulator. According to the theory of hydraulic fluid dynamics, the pressure loss through the tube can be written as 

[DELTA][p.sub.l] = [lambda] l/d [rho][v.sup.2.sub.F]/2 sign ([v.sub.F]), (13)

where [lambda] is the flow restriction coefficient along the tube,

[mathematical expression not reproducible], (14)

the Reynolds number is defined as

Re = [v.sub.F]d/v (15)

and the sign function is given as

[mathematical expression not reproducible], (16)

where [v.sub.F] is the fluid flowing velocity, l and d are the length and flow diameter of the damping tube, and v and [rho] are the kinematic viscosity and the density of the fluid.

(3) Mathematical model of the accumulator: the gas chamber of the accumulator is filled with pressurized nitrogen. The state changing process of the gas can be described by the following equation.

p[V.sup.n] = [p.sub.0][V.sub.0.sup.n], (17)

where n is the gas polytrophic exponent, V is the actual gas volume, [V.sub.0] is the initial gas volume, p is the actual gas pressure, and [p.sub.0] is the initial gas charge pressure.

If the accumulator is loaded or unloaded rapidly, the thermodynamic process of the gas state change belongs to an adiabatic process and the value of the polytrophic exponent is n = 1.4. Otherwise, the thermodynamic process is the isothermal process and n = 1 [25, 26]. When considering its application in a real vehicle, which is under typical road disturbances, the reciprocating impact on the accumulator is usually instantaneous. However, certain heat exchange processes exist, so the state change can be considered as a process that is between adiabatic and isothermal. In this paper, the polytrophic exponent n is set as 1.3 [27, 28].

(4) Interface module: the function of the interface module is to carry out the exchange of the simulation data between the two simulation platforms, namely, RecurDyn and AMESim. Here, each program executes its respective simulation simultaneously. At each time step (0.001 sec), both codes update one another with new state values before advancing to the next step. Simulations begin when RecurDyn calculates the stroke and velocity of the actuation cylinder. The hydropneumatic system in AMESim calculates the force of the actuation cylinder and feed it back to the MBD vehicle model.

The interface module has 14 input ports and 28 output ports. The input ports f 1~f 12 are designated to the actuation cylinders' output force and f 13~f 14 denote the left and right track tensioning cylinders' output force. The output ports 1~ 24 represent the stroke and velocity of the actuation cylinders and 25~28 are for the track tensioning cylinders. The interface module and its connection with the suspension system are depicted in Figure 3.

3. Model Verification

In order to validate the model, a single hydropneumatic unit is tested using the suspension test rig shown in Figure 4. The parameters in Table 3 are used in the simulation models from (13) to (17) and are in accordance with those in the test rig. In the AMESim parameter setting mode, the values are assigned to each hydraulic model. Relative to all the forces in the actuation cylinder, the friction force is comparatively tiny. Thus, a simplification is to set the value of the fiction parameter to the system default, which is suitable for most cases. The vehicle is designated to traverse a triangle barrier at the speed of 5 kmph as illustrated in Figure 5. Based on the above conditions, the disturbance on the road wheel is clear. On the test rig, the excitation of the vibration table is also set the same as that in the simulation.

Figure 6 shows the test and simulation results of the actuation cylinder (L1) stroke. In the initial position, the track assemblies are not in contact with the ground; otherwise, the initial impact of the track on the ground would result in the simulations failing. Thus, the vehicle descends to the ground firstly and the actuation cylinders are compressed from the initial extended position to the static balance position. This is shown as an initial stroke compression in the simulation curve in Figure 6. After approximately one second, the vibrations are attenuated almost to zero and the vehicle reaches its static balance position. The change in stroke after 1.5 s is caused by the vehicle acceleration and the small subsequent oscillations are from the impact of road wheel with the gap between the shoe plates. The large compression in stroke after 4 s is the result of the road wheel bouncing when traversing the barrier. The 2nd~6th road wheel passing over the barrier resulted in the signal oscillations after 5 s. In general, the simulation results agree well with the test data.

The test and simulation results of the road wheel's (L1) acceleration are illustrated in Figure 7. In the simulation curve, except for the large acceleration caused by the barrier impact, there exist some small fluctuations induced for the same reason as discussed in Figure 6. In general, the test and simulation results matched very well and the constructed simulation model is validated.

4. Cosimulation Analysis

Using the modeling and simulation environment, the virtual tests on the G road are conducted under different conditions such as vehicle speed, flow diameter of the damping tube; accumulator charge pressure, and track tensioning pressure. The influences of the above parameters on the vehicle ride safety are investigated.

4.1. Effect of the Damping. The damping of the hydropneumatic suspension is directly set by the flow diameter (d) of the tube installed between the cylinder and the accumulator. The simulations which explore the effect of different tube diameters are conducted at different vehicle speeds. Figure 8 shows that there is always an optimal diameter which renders the RMS value of the track-wheel gap ([D.sub.rms]) to be minimum for any speeds. The optimal diameter for minimum [D.sub.rms] decreases as the speed increases. The decrease of diameter on the left side of the optimal value has a stronger influence on [D.sub.rms] than the increase of the diameter on the other side.

Similar to Figure 8, it can be seen from Figure 9 that there always exists an optimal diameter which minimizes the RMS value of the road wheel dynamic force ([F.sub.rms]), for any vehicle speeds. On the left side of the optimal value, the decrease of the diameter also has a stronger influence on [F.sub.rms]. As the vehicle speed increases, the optimal diameter decreases. The higher the vehicle speed is, the stronger the effect the diameter has. Lower speed generally results in less dynamic force, but an improperly tuned diameter may cause higher dynamic force for a lower speed than that for a higher speed.

The RMS value of the stroke ([S.sub.rms]) of the actuation cylinder for different damping diameters and vehicle speeds on a G road are shown in Figure 10. It can be seen that [S.sub.rms] increases as the diameter and vehicle speed increase. For 47 kmph, as the diameter is more than 10 mm, [S.sub.rms] exceeds one-third of the designed stroke length; thus an end-stop impact of the suspension will be a possible event.

4.2. Effect of the Accumulator Charge Pressure. The accumulator charge pressure determines the basic stiffness of the suspension. Through simulations and analysis of various charge pressures, the effects of the charge pressure on the RMS value of the track-wheel gap, the road wheel dynamic force, and the stroke are obtained.

Figure 11 shows how the RMS value of the wheel-track gap ([D.sub.rms]) changes as the charge pressure ([p.sub.c]) changes. As the pressure decreases down to 2 MPa, [D.sub.rms] increases slowly. Once the pressure drops below 2 Mpa, [D.sub.rms] increases rapidly. Thus, if [D.sub.rms] is over one-third of the grouser height, the failure of the track separating from the wheel will be possible. It is known that the height of the grouser is about 58 mm. Therefore, if the pressure is less than about 2 MPa, the probability of occurrence of the failure increases significantly.

Variations of the RMS value of the road wheel dynamic force ([F.sub.rms]) with charge pressure ([p.sub.c]) are illustrated in Figure 12 for various vehicle speeds. It can be seen that [F.sub.rms] increases with an increase in charge pressure at all speeds.

The RMS value of the stroke (Srms) for different vehicle speeds and charge pressures are given in Figure 13. Notice that when the pressure is less than about 1.5 MPa, [S.sub.rms] increases sharply as [p.sub.c] decreases. The designed stroke length of the actuator cylinder is 200 mm; if [p.sub.c] drops below about 1.5 MPa at the speed of 47 kmph, the end-stop impact of the suspension has a high probability of occurrence.

4.3. Effect of the Track Tensioning Force. The track tensioning force and its fluctuations play an important role in the ride safety of the vehicle. A decline of tensioning force will cause a rise in fluctuations and more frequent track jumping from the sprocket. Thus, the failure of the track separating from the sprocket will occur. In the hydropneumatic suspension system, a pair of hydraulic tensioning cylinders is used to adjust the track. The tensioning force is directly related to the actuating pressure [p.sub.T] and the track tensioning force fluctuations can be evaluated by its standard deviation ([F.sub.D]).

The RMS value of jumping quantity ([W.sub.rms]) and the standard deviation of the tensioning force ([F.sub.D]) of the track for various pressures are given in Figure 14. Both [W.sub.rms] and [F.sub.D] increase as the pressure decreases. When the jumping quantity is over the length of the effective mesh area on the sprocket teeth, the failure of the track separating from the sprocket may happen. A tighter track provides better ride safety, but, on the other hand, the force acting on the sprocket, the idler and the roller by the track will be higher. Thus, the friction force between the track and the above components increases . The loss of the chassis output power on the running system will also be greater.

5. Conclusions

Factors which influence the ride safety of a tracked vehicle are analyzed. Systematic evaluation parameters of ride safety and their criteria are proposed. Using a cosimulation technique, the vehicle MBD model and the hydropneumatic suspension model are built and verified. Through simulations at various vehicle speeds, accumulator charge pressures, and damping diameters, the effect they have on ride safety are investigated and quantified. The key conclusions are as follows.

(1) Different optimal damping parameters exist for the minimum wheel-track gap and the wheel dynamic force at different vehicle speeds. The value of the optimal diameter decreases as the vehicle speed increases. If the diameter is not tuned well, the dynamic wheel force at a low speed may be higher than that at a high speed.

(2) As the accumulator charge pressure decreases, the RMS value of the wheel-track gap and the stroke increases but the RMS value of road wheel dynamic force decreases. For the hydropneumatic suspension studied in this paper, when the pressure drops below about 2 MPa, the probability of the end-stop impact of the suspension increases significantly.

(3) As the actuating pressure of the track tensioning cylinder decreases, the track tensioning force decreases. An insufficient tensioning force leads to a looser track, a rise in fluctuations, and more frequent track jumping from the sprocket. Thus, the ride safety deteriorates. In order to improve the ride safety, the track tensioning force should be increased sufficiently.

https://doi.org/10.1155/2017/1256536

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this manuscript.

References

 X. Dong, The dynamic analysis and simulation of tracked vehicles [Ph.D. thesis], Hunan University, Changsha, China, 2007.

 W. Wang, H. Huang, L. Gu et al., "Dynamic simulation on suspension vibration characteristic of tracked vehicle," Journal of Vibration and Shock, vol. 25, no. S, pp. 1016-1018, 2006.

 J. Y. Wong, "Dynamics of tracked vehicles," Vehicle System Dynamics, vol. 28, no. 2-3, pp. 197-219, 1997.

 Y. Zhu, Y. Di, and X. Guo, "Simulation of running system of tracked vehicle's performances based on ADAMS," Drive System Technique, vol. 2, no. 2, pp. 29-31, 2008.

 M. C. Levesley, R. Ramli, N. Stembridge, and D. A. Crolla, "Multi-body co-simulation of semi-active suspension systems," Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, vol. 221, no. 1, pp. 99-115, 2007.

 M. Haiba, D. C. Barton, P. C. Brooks, and M. C. Levesley, "Using a quarter-vehicle multi-body model to estimate the service loads of a suspension arm for durability calculations," Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, vol. 217, no. 2, pp. 121-133, 2003.

 M. C. Levesley, S. A. Kember, D. C. Barton, P. C. Brooks, and O. M. Querin, "Dynamic Simulation of Vehicle Suspension Systems for Durability Analysis," Materials Science Forum, vol. 440-441, pp. 103-110, 2003.

 Y. G. Liao and H. I. du, "Cosimulation of multi-body-based vehicle dynamics and an electric power steering control system," Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, vol. 215, no. 3, pp. 141-151, 2001.

 J. Maiorana, B. P. Minaker, D. Zhang et al., "Cosimulation of active suspension," in Proceedings of the SAE 2005 World Congress & Exhibition, pp. 363-371, 2005.

 S. Dietz, G. Hippmann, and G. Schupp, "Interaction of vehicles and flexible tracks by co-simulation of multibody vehicle systems and finite element track models," Vehicle System Dynamics, vol. 37, pp. 372-384, 2003.

 FunctionBay Inc., "RecurdynV8R3," 2014, http://www.recu.rdyn.com.

 LMS Imagine Lab, "AMESim14.0," 2015, http://www.amesim.com.

 Y. Zhang, L. Chen, and W. Luo, "Simulation of off wheel problem of high speed tracked vehicle based on Recurdyn," Coal Mine Machinery, vol. 35, no. 2, pp. 49-51, 2014.

 Y. Defu, Theory and Practice of Vehicle Suspension Shock Absorber, Weapon Industry Press, 2003.

 F. Ding, Dynamics of Suspension System of Tracked Armored Vehicle, National Defence Industry Press, 2004.

 U. Solomon and C. Padmanabhan, "Semi-active hydro-gas suspension system for a tracked vehicle," Journal of Terramechanics, vol. 48, no. 3, pp. 225-239, 2011.

 C. Zhou, Vehicle Suspension Design and Theory, Peking University press, 2011.

 Z. Tao, G. Zhiqiang, and Z. Zhili, "Simulation of off-wheel problem of high speed tracked vehicle's running system," Journal of Henan University of Science & Technology, vol. 27, no. 6, pp. 12-15, 2006.

 J. Guo, L. Gu, and X.-L. Li, "Effect of the damping of suspension system on the riding comfort of a tracked vehicle," Beijing Ligong Daxue Xuebao/Transaction of Beijing Institute of Technology, vol. 21, no. 4, pp. 440-444, 2001.

 H.-B. Jiang, J.-T. Geng, X.-L. Zhang, and T. Wang, "Study on co-simulation of vehicle suspension system employing inerter based on virtual protoype model," Zhendongyu Chongji/Journal of Vibration and Shock, vol. 29, no. 12, pp. 221-223, 2010.

 E. Guglielmino, T. Sireteanu, C. W. Stammers, G. Ghita, and M. Giuclea, Semi-Active Suspension Control: Improved Vehicle Ride and Road Friendliness, Springer, London, UK, 2008.

 W. Bauer, Hydropneumatic Suspension Systems, Springer Heidelberg, London, UK, 2010.

 Y.-C. Wang, X.-H. Gao, J.-H. Guo, and C.-Q. Zhang, "Research on virtual prototype of hydro-pneumatic suspension," Xitong Fangzhen Xuebao / Journal of System Simulation, vol. 18, no. 8, pp. 2183-2186, 2006.

 J. Wei, H. Du, X. Fang, and K. Guo, "Road-friendliness of inter-connected hydro-pneumatic suspension based on ADAMS/ Simulink/AMESim," Nongye Jixie Xuebao/Transactions of the Chinese Society of Agricultural Machinery, vol. 41, no. 10, pp. 11-17, 2010.

 G. Lin, D. Ma, and Z.-L. Zhu, "Vibration of a missilelauncher with a multi-spindled interconnected hydro-pneumatic suspension," Journal of Vibration and Shock, vol. 32, no. 12, pp. 144-148, 2013.

 D. Zhou, G. Wang, Y. Liu, and X. Guo, "Effects of parameters of hydro-pneumatic suspension cylinder on vehicle ride comfort," Nongye Jixie Xuebao/Transactions of the Chinese Society of Agricultural Machinery, vol. 35, no. 5, pp. 25-28, 2004.

 H. Wang, The vibration dynamics of the heavy vehicle with composite interconnecting oil-air suspension [Ph.D. thesis], Beijing Institute of Technology, Beijing, China, 1996.

 J. Liu and H. Wang, "Analysis on stiffness characteristic of hydro-pneumatic suspension system of a rubber-wheel vehicle," Mining & Processing Equipment, vol. 37, no. 13, pp. 44-47, 2009.

 B. Liu, Study on dynamics of driving device of high-speed tracked vehicle [Ph.D. thesis], Beijing Institute of Technology, Beijing, China, 2015.

Shousong Han, Zhiqiang Chao, and Xiangbo Liu

Department of Mechanical Engineering, Academy of Armored Forces Engineering,

Beijing 100072, China

Correspondence should be addressed to Shousong Han; hssbeijing@163.com

Received 14 February 2017; Revised 24 May 2017; Accepted 15 June 2017; Published 20 September 2017

Caption: FIGURE 1: Vehicle dynamic model and G road model.

Caption: FIGURE 2: Connection sketch of hydropneumatic suspension.

Caption: FIGURE 3: Hydraulic system model of the hydropneumatic suspension.

Caption: FIGURE 4: Test rig of the hydropneumatic suspension.

Caption: FIGURE 5: Modeling and simulation of the vehicle traversing a barrier.

Caption: FIGURE 6: Test and simulation results of the stroke.

Caption: FIGURE 7: Test and simulation results of the road wheel acceleration.

Caption: FIGURE 8: The RMS value of the track-wheel gap versus the flow diameter.

Caption: FIGURE 9: The RMS value of the road wheel dynamic force versus the flow diameter for various speeds.

Caption: FIGURE 10: The RMS value of the stroke versus the flow diameter for various speeds.

Caption: FIGURE 11: The RMS value of the track-wheel gap versus the accumulator charge pressure.

Caption: FIGURE 12: The RMS value of the road wheel dynamic force versus the accumulator charge pressure for various speeds.

Caption: FIGURE 13: The RMS value of the stroke versus the accumulator charge pressure for various speeds.

Caption: FIGURE 14: The standard deviation of track tension force and the RMS value of jumping quantity versus the actuating pressure of the track tensioning cylinder.
```TABLE 1: Evaluation parameters and criteria of the ride safety.

Parameter                Symbol      Unit

RMS value of the       [D.sub.rms]    mm
track-wheel gap
RMS value of the       [F.sub.rms]    kN
wheel dynamic
force
RMS value of the       [C.sub.rms]    mm
stroke of the
actuation cylinder
RMS value of the       [W.sub.rms]    mm
jumping quantity
of the track

Parameter                             Criteria

RMS value of the          <(1/3) height of the track teeth
track-wheel gap
wheel dynamic
force
RMS value of the            <(1/3) designed stroke length
stroke of the
actuation cylinder
RMS value of the       <(1/3) effective engage length between
jumping quantity        the track and the sprocket teeth
of the track

TABLE 2: Suspension initial value of the design variables.

Variable                     Symbol        Unit         Value

Length of axle arm              R           m           0.40
Road wheel diameter             d           m           0.60
Angle of the axle        [[alpha].sub.j]   deg          -14.0
arm from x at the
static position
Initial angle of the     [[alpha].sub.0]   deg          -28.0
axle arm from x
Angle of the axle        [[alpha].sub.d]   deg          30.0
arm from x at
the top end
Designed                    [D.sub.w]       m           0.39
displacement of the
Designed stroke of              S           m           0.20
the actuation
cylinder
Section area of the         [A.sub.s]       m2    4.4 x [10.sup.-3]
actuation cylinder
Working pressure of             p           Pa        Variable
the actuation
cylinder
Standard atmosphere         [p.sub.a]       Pa    1.01 x [10.sup.5]
pressure
Displacement of the             f           m         Variable
The x axis position             w           m         Variable
of the piston
Velocity of the               [??]         m/s        Variable
actuation cylinder
Differential             [DELTA]p ([??])    Pa        Variable
pressure of the
hydraulic restrictor
Angle of the axle            [alpha]       deg        Variable
arm from x
Angle of the crank           [beta]        deg        Variable
from y
The angle between            [gamma]       deg        Variable
the crank and the
actuation cylinder
The gas polytrophic             n           --           1.3
exponent

TABLE 3: Initial value of hydraulic variables.

Hydraulic variables           Symbol         Unit

Initial charge pressure      [p.sub.0]        Pa
of the accumulator
Volume of the accumulator    [V.sub.0]    [m.sup.3]
Length of the damping tube       l            m
Flow diameter of the             d            m
damping tube
Section area of the          [A.sub.s]    [m.sup.2]
actuation cylinder
Initial track                [p.sub.T]        Pa
tensioning pressure
Section area of the          [A.sub.s]    [m.sup.2]
track tensioning
cylinder
Density of the oil             [rho]     kg/[m.sup.3]
Kinematic viscosity              v           cSt
of the oil

Hydraulic variables                Value

Initial charge pressure       6.8 x [10.sup.6]
of the accumulator
Volume of the accumulator    0.85 x [10.sup.-3]
Length of the damping tube          1.2
Flow diameter of the          6 x [10.sup.-3]
damping tube
Section area of the          4.4 x [10.sup.-3]
actuation cylinder
Initial track                 6.2 x [10.sup.6]
tensioning pressure
Section area of the          2.0 x [10.sup.-3]
track tensioning
cylinder
Density of the oil           0.87 x [10.sup.3]
Kinematic viscosity                  10
of the oil
```
Title Annotation: Printer friendly Cite/link Email Feedback Research Article Han, Shousong; Chao, Zhiqiang; Liu, Xiangbo Shock and Vibration Report Jan 1, 2017 5761 Probabilistic Entropy EMD Thresholding for Periodic Fault Signal Enhancement in Rotating Machine. An Experimental Research on the Vibration of Surface-Textured Journal Bearings. Automotive safety Dynamics Dynamics (Mechanics) Mathematical models Suspension systems