Abstract
In order to conveniently analyze the dynamic performance of tracked vehicles, mathematic models are established based on the actual structure of vehicles and terrain mechanics when they are moving on the soft random terrain. A discrete method is adopted to solve the coupled equations to calculate the acceleration of the vehicle’s mass center and tractive force of driving sprocket. Computation results output by the model presented in this paper are compared with results given by the model, which has the same parameters, built in the multi-body dynamic software. It shows that the steady state calculation results are basically consistent, while the model presented in this paper is more convenient to be used in the optimization of structure parameters of tracked vehicles.
1. Introduction
Tracked vehicles are widely used in the field of military, agriculture, and construction industries. They are applied on the ground that has a low carrying capacity or wheeled vehicles cannot go through (swamp, snowfield, etc.) or are used for obstacle-crossing and climbing, and so forth because the tracks can be treated as self-carried road to increase the vehicle’s passing ability. Nowadays, there are lots of modeling methods about tracked vehicles which primarily treat the tracks as flexible [1–5] or rigid ones [6–9].
Ma and Perkins employed a nonlinear finite element method to represent the track segment contacting with the terrain [1, 2]; Park et al. divided the track system contacting with terrain into a finite number of small elements and used structural analysis and the numerical iterative method to calculate the drawbar pull force [3]; Ma treated the tracks as flexible tracks and exerted the excitation directly on the road wheels to reduce the complexity of the model [4]; Han et al. regarded tracks as flexible ones and used springs to replace tracks between road wheels in order to build the ride comfort model of the tracked vehicle [5].
The assumption of flexible tracks cannot fully reflect the dynamic performance of tracks. Rigid models are adopted by more researchers. Merhof and Hackbarth illustrated the actual structure of the tracks and set up their dynamic equations [6]; Rubinstein and Hitron described each track’s positions, movements, and the interaction with terrain by coordinate transformation in three-dimensional space [7]; Ma et al. regarded the tracks and the road wheels as a whole and analyzed each part in the track system [8]; Hu took heavy tracked vehicle as a research object and divided the track system into 4 parts, namely, the tracks in contact with terrain, the tracks in contact with front idler, the tracks in contact with driving sprocket, and the tracks between front idler and driving sprocket [9].
Although the assumption of rigid tracks is approaching the truth, the effects of structure parameters on the performance are difficult to analyze, and generally the model cannot be used in the structural optimization. Therefore, this paper presents a simplified dynamic model based on the assumption of rigid tracks. A discrete method is used to solve the coupled equations. The computation results, namely, the acceleration of the vehicle’s mass center and the tractive force of driving sprocket, are compared with those obtained from RecurDyn to verify the model’s validity.
2. Dynamic Model of Tracked Vehicles
2.1. Random Road Excitations
In this paper grade road surface is taken as the random excitations for the tracked vehicles. The road model is built by making use of power spectral density function in time frequency domain [10–12] and the iteration of trigonometric series method [13]; the result is illustrated in Figure 1. The tracked vehicle’s highest speed is 60 km/h, the length of a track is 0.153 m, and so the time that a wheel passes by a track is 0.009 s. Thus the time interval for the discrete sampling of terrain excitation is 0.009 s in order to simulate the vehicle run on the road surface.
2.2. Track-Terrain Contact
Currently the most widely used model for analyzing soil sinkage is Bekker’s equation [14, 15]:where, is the pressure, kpa; is the pressure-sinkage exponent; is the soil cohesive coefficient, kN/mn+1; is the soil internal friction coefficient, kN/mn+2; is the width of a track, m; is the sinkage, m.
Actually the tracked vehicle will have repeated loads on the same terrain when different wheels of the vehicle pass by it. This phenomenon can be more precisely described by the loading-unloading curve of soil model introduced by Wong [16] than Bekker’s equation, as shown in Figure 2.
As it can be seen from Figure 2, OA is the continuous loading curve which can be described by (1); BA is the unloading and repeated loading curve in which the slope can be seen as a constant value , and the curve can be described by [17]where is the maximum sinkage, m; is the sinkage, m.
In order to build and solve the dynamic models of tracked vehicles, some simplifications, which have little effects on the final calculation results, are needed: every is sequentially exerted on the two pins of a track so that the terrain under the track would be complete and continuous; this will lay the foundation for the calculation of shear forces acting on tracks by terrain; in reality the angle of track changes when a boogie passes by it. However, since the time is very short, it is reasonable to assume that the track keeps the same angle during this period so that the coupled relationships among wheel, track, and terrain can be decoupled; the average sinkage of both pins of a track after the wheel leaves it is regarded as the sinkage of the track when the wheel is on the middle of it, as shown in Figure 3.
At the first sampling time, terrain excitations on both pins of the track under the first road wheel are original excitations and . At the second sampling time, the excitations on the track under the first wheel are and that has sinkage and springback after the first sampling time and so are the following excitations on the track. For the tracks under the second to the sixth road wheels, the terrain excitations on both pins of them are excitations that have sinkage and springback after the former sampling time. Besides, there are 2 tracks at the right and left sides of one terrain excitation and at least one road wheel has passed by before. Therefore, all of the former sinkage and springback should be taken into consideration; namely, the original excitation should subtract the maximum average sinkage among all the ones and then add the springback to form the new excitations acting on the later tracks.
2.3. Mechanical Analysis of Tracks
In this model the inertial force and inertial moment of a track are neglected due to the fact that the mass of a track is very small compared with the whole vehicle. Thus, the forces and moments acting on a track can be illustrated in Figure 4.
In order to compute the shear force acting on a track, the first step is to calculate , and , . Here, it is supposed that the time period that a road wheel passes a track is . During this period the angle of a track, vertical and angular acceleration of the vehicle’s mass center, the vertical force acting on a track from road wheel , and the sinkage of excitation are calculated every , other physical variables are calculated every . Therefore, there are 4 calculation points in except the midpoint. Figure 5 illustrates this concept.
Firstly, let us analyze the interaction between the first road wheel and the track under it, and the analysis result also applies to other wheels and tracks. It is supposed that the tracked vehicle is initially static on a level terrain. After a the angle of the track contacting with the road wheel (, ) changes to because of the excitation from terrain. , acting on the track can be obtained from Figure 6:
Two force components of the pretensioning force, namely, and , are acting on the right pin. Figure 7 illustrates all the forces acting on the track at the first calculation point in the time period .
The pressure distribution acting on the track from terrain is obtained by Figure 7 and the shear force (mentioned in Section 2.4) at the first calculation point can be calculated.
In order to solve the coupled equations, , , , , , and are kept unchanged at the next calculation point. At the second calculation point, , , and can be obtained according to equilibrium conditions. For instance, in Figure 8, we can have
In Figure 9, moves forward, the pressure distribution acting on the track from terrain and the shear force can be calculated by , , and (5) and the unchanged values of , , , , , and . Based on this method, the similar values can be obtained when the road wheel is at the midpoint, , and of the track.
In the second time period , road wheel () moves on to track (). At this junction point the tractive force acting on tracks () and () keeps unchanged. According to Figure 10, after the pressure distribution and shear force are obtained at the fourth calculation point in the first time period , , , and can be calculated as follows:
As illustrated in Figure 11, after the road wheel moves onto track (), , can be obtained from , which is the sum of force components of , in the direction along with track . equals (obtained by (8)):
To obtain the values at the following calculation points, the method is the same as the calculation points in the first . When the sixth time period comes, random terrain excitation begins to be exerted on the track under the second road wheel because there are 5 tracks between two wheels, as illustrated in Figure 12.
In Figure 12, , and act on the same pin as , , and . So they are active and reactive forces (moments). Then , , and can be computed by the equilibrium condition of track (); namely,
The method to calculate , , and is the same as the one used to calculate the force or moment values of the track under the first road wheel in the first . In this way, forces or moments acting on the tracks under and between the road wheels can be obtained.
As illustrated in Figure 13, suppose track contacts with the sixth road wheel. After calculating , , and , can be obtained bywhere is the driving force given by sprocket, kN; is the rotation angle of the vehicle’s mass center, rad; is the angle between the tracks connected with the driving sprocket and the level terrain when the vehicle is static, rad.
2.4. Shear Forces Acting on Tracks
When tracks are pulled from terrain by the driving sprocket, tractive forces will be exerted on them by terrain. On the soft terrain, the tractive force is mainly due to the shear force. Therefore, the tractive force given by the driving sprocket can be obtained by calculating the shear force acting on tracks from terrain [18]. When a road wheel passes by a track in a time period , the pressure distribution under the track changes according to the wheels’ positions, and so do forces and moments on both pins of the track. Figure 14 illustrates the simplification of forces and moments acting on a track at a specific calculation point.
After the simplification, the forces and moments exerted on a track are equal to a force and a moment . For a rigid and small track, the pressure distribution under the track can be regarded as a linear distribution [19]. There are 3 different situations according to the resultant force and moment , as illustrated in Figure 15:where , are the maximum and minimum values of the distribution, kpa; is the length of a track, m; is the width of a track, m; is the eccentricity, , m.
(a)
(b)
(c)
A local coordinate can be placed on a track so that the pressure distribution function can be described aswhere is the pressure distribution function under a track, kpa; is the slope; is the intercept, m.
Substitute (13) into Janosi shear stress equation [20] (14), shear force acting on tracks from terrain according to different pressure distribution functions can be calculated:
Equation (14) can be reduced towhere is the distance that road wheels pass by a track, m; is the width of a track, m; is the soil cohesion coefficient, kN/m2; is the internal friction angle of soil, rad; is the shear modulus deformation of soil, m; is the sliding rate of vehicle.
2.5. Dynamic Equations of Tracked Vehicles
The following assumptions are made to set up dynamic equations of tracked vehicles: the vehicle is bilaterally symmetrical; excitations exerted on both left and right tracks are exactly the same; road wheels and tracks are rigid. Based on these assumptions, a half vehicle model can be used to build the dynamic model of the vehicle, as illustrated in Figure 16.
Initially the vehicle is static and in equilibrium. So the dynamic equations in direction are as follows:where is the mass of half vehicle, kg; is the displacement of mass center in direction, m; is the forces acting on vehicle’ s body from the suspension system, kN; , is the forces acting on front idler and driving sprocket of tracks, kN; , is the force components in direction acting on front idler and driving sprocket of tracks, kN; , is the angle between the tracks connecting the driving sprocket and the level road when the vehicle is static, rad; is the rotation angle of the mass center of the vehicle, rad.
Pitching movement equations for the vehicle body,where is the inertial moment when the vehicle rotates around the mass center, kg·m2; is the force components in direction acting on front idler and driving sprocket, kN; , is the horizontal and vertical distances of road wheels’ mass center from vehicle’s mass center, are positive, are negative; , is the horizontal and vertical distances of front idler’s and driving sprocket’s mass center from vehicle’s mass center, m; is positive, is negative.
Forces exerted by suspension system,
Since the mass center of the tracked vehicle rotates slightly when it is moving on the grade E random road surface, sowhere is the damping coefficient, kN/(m/s); is the spring coefficient, kN/m; is the excitations exerted on a track’s right pin, m; is the stroke of the suspension system, m.
The vertical force acting on the track from road wheel, ,where is the mass of a road wheel, kg; is the acceleration of gravity, m/s2; is the road wheel’s acceleration, m/s2; is the angle of track , rad.
3. Decoupling of Dynamic Model
3.1. Coupled Dynamic Model
There exist couplings in the dynamic model of tracked vehicles: if the shear force under a track is to be calculated, it is necessary to know the forces and moments acting on both pins of the track and the vertical force acting on the track by road wheels. But in order to get these forces and moments, the shear force must be first calculated, as mentioned in Section 2.3; in order to solve (16)–(18), the shear force and displacements of springs must be known. But, as described in , the shear force should be obtained through the vertical force acting on the track from road wheel and so are the displacements of springs. Obviously, the vertical force is related to the vibration and pitching movement of the tracked vehicle’s mass center, which means (16)–(18) are required to calculate the force. The coupling relationships among them are illustrated in Figure 17.
3.2. Model Decoupling Process
The process is shown in Figure 18. Discretization is adopted to solve the coupled dynamic model. A very short sampling time is taken to precisely calculate the derivation and integral of equations or to obtain all the parameters that are needed to calculate other ones, and then substitute the parameters obtained at the former calculation point into the equations to get the parameters for the next calculation point. In the decoupling process, the pressure distribution and the shear force are calculated every when a road wheel passes by a track. The vertical and angular velocities of the mass center and so forth are calculated every , in which , calculated at the fourth calculation point in a , must be known beforehand. When the road wheel is at the midpoint of a track, the method introduced in Section 2.3, namely, the calculation equations for the pressure distribution at every calculation point, can be adopted to get the pressure distribution values, and thus the average sinkage of the track can be obtained.
4. Case Study
The main parameters of a tracked vehicle are shown in Table 1. The slipping rate of the vehicle , the sampling time period s. Parameters of the soil, that is, sandy loam, are kN/mn+1, kN/mn+2, kN/m2, rad, m, and , and the maximum sinkage = 0.2 m.
A virtual prototype model of the tracked vehicle was built using the module “Track-HM” in RecurDyn to test the accuracy of calculation, as illustrated in Figure 19. The parameters used to build the virtual prototype are the same as the ones used in this paper, as shown in Table 1. Considering it needs some time for the virtual prototype model to accelerate to the required speed, a Step function is used to drive the model: where is the time, s; is the initial time, 0; is the initial function value, 0; is the end time, 6 s; is the maximum angular velocity of the driving sprocket. The maximum speed of vehicle is 60 km/h, the pitch radius of driving sprocket is 0.25 m, so its value is 67 rad/s.
The vehicle model is initially static and the friction between components is neglected. When the vehicle is in a uniform rectilinear motion after 6 seconds, the vertical acceleration of the mass center of the vehicle and the tractive force of driving sprocket can be compared with the results of this paper.
It can be seen from Figure 20 that when the tracked vehicle reaches stable state, the vertical acceleration of the vehicle body’s mass center calculated by the model presented in this paper is basically consistent with the result from RecurDyn, and the acceleration values are in the range of −6 m/s2~6 m/s2, which meets the requirement that the maximum vertical acceleration of tracked vehicle’s mass center should be less than 1 g. The maximum values of the model given in this paper and RecurDyn are 5.808 m/s2 and 5.938 m/s2, respectively. The error between them is 2.19%.
It can be seen from Figure 21 that when the tracked vehicle reaches stable state, the tractive force calculated by our model is basically consistent with the result from RecurDyn if the friction among components of the vehicle and centrifugal force of the tracks are not taken into consideration. The error between our model and RecurDyn is 8.77%.
The results of the model presented in this paper have some deviations from the results output by RecurDyn. This is because our model is built without consideration of inertia forces of tracks and slippage of the vehicle. When a road wheel passes by a track, the track will certainly rotate a little angle. Although the time period that a road wheel on a track is very short, neglecting this angular change may affect the result. Discretization of the computation process is used to decouple the coupled problems, computation errors may be inevitably induced into the results. However, generally the computation errors are in an acceptable range, which indicates the model introduced in this paper is useful for the future optimization of design variables.
5. Conclusion
This paper introduces a dynamic model of tracked vehicle when the vehicle moves on the soft random terrain. The contributions of the paper include the soft random terrain excitation is introduced into the model, which makes it much closer to the reality; a detailed track modeling method is introduced; a decoupling method is introduced to solve the coupled dynamic model. Future work would be directed towards improving the modeling of tracks by taking the inertia forces, frictions, and so forth into consideration, so that the computation errors can be further reduced.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research work is funded by the National Natural Science Foundation of China with the Grant nos. 51175457 and 51475416. The first author is also supported by the Fundamental Research Funds for the Central Universities of China.