Tractor-implement vibration model is developed for simulation as shown in Fig. 2 and mass of tractor with implements and moment of inertia relative to central axis is used to represent tractor body structure. Due to surface roughness front and rear tires of tractor are subjected to displacement excitations qfz(t) and qrz(t) respectively. Because there is a difference in amplitudes of qfz(t) and qrz(t) therefore angular displacement with respect to centre of mass arises. Due to this angular displacement with respect to mass centre roll, yaw and pitch movements of tractor came into the picture.
The ground excitation relation between front and rear axles of the tractor-implement system at any instant i by considering the time lag between qfz(t)and qrz(t) can be expressed as
$${q}_{fzi}\left(t\right)={q}_{rzi}\left(t+\tau \right) \left(1\right)$$
Where qfzi(t) and qrzi(t) are displacement excitations at front and rear axles respectively, Ʈ is time lag and can be calculated as
$$\tau =\frac{{l}_{bf}+{l}_{br}}{v} \left(2\right)$$
Where lbf is distance of mass centre between chassis and front axle, lbr is the distance of mass centre between chassis and rear axle and v is the velocity of tractor-implement system.
The equations of motion can be expressed as
$$\left({M}_{T}+{M}_{I}\right)\ddot{z}={f}_{1}+{f}_{2}-\left({M}_{T}+{M}_{I}\right)g \left(3\right)$$
The displacement of wheels can be calculated as
$${z}_{1}=z+\left({l}_{1}+{l}_{3}\frac{{M}_{I}}{{M}_{T}+{M}_{I}}\right) \left(4\right)$$
\({z}_{2}=z+\left({l}_{2}-{l}_{3}\frac{{M}_{I}}{{M}_{T}+{M}_{I}}\right) \left(5\right)\) The above equations are used to calculate the effect of tillage implement on displacement of wheels.
The mathematical model of tractor-implement system is shown in Fig. 3. The model represents a real system with parameters to passenger seat mass (ms), force transmitted to the front left and right corner of tractor chassis (Ff1, Ff2 ),Tractor body mass (mb) tyre stiffness (kt), damping coefficients (ct) vertical axle displacements (zf1, zf2, zr1, zr2). Road excitations have been represented by Zr. Pitch and roll of tractor-implement system are represented by Ɵ and ϕ respectively.
The vibratory forces transmitted by each wheel at the point of contact between the axle and the wheel can be formulated as follows:
$${k}_{t}\left({z}_{a}-{z}_{R}\right)+{c}_{t}\left({\dot{z}}_{a}-{\dot{z}}_{R}\right) \left(6\right)$$
The governing dynamical equations for the forces transmitted to chassis are given as:
$${F}_{f1}={k}_{tf}\left({z}_{af1}-{z}_{Rf1}\right)+{c}_{tf}\left({\dot{z}}_{af1}-{\dot{z}}_{Rf1}\right) \left(7\right)$$
$${F}_{f2}={k}_{tf}\left({z}_{af2}-{z}_{Rf2}\right)+{c}_{tf}\left({\dot{z}}_{af2}-{\dot{z}}_{Rf2}\right) \left(8\right)$$
$${F}_{r1}={k}_{tr}\left({z}_{ar1}-{z}_{Rr1}\right)+{c}_{tr}\left({\dot{z}}_{ar1}-{\dot{z}}_{Rr1}\right) \left(9\right)$$
$${F}_{r2}={k}_{tr}\left({z}_{ar2}-{z}_{Rr2}\right)+{c}_{tr}\left({\dot{z}}_{ar2}-{\dot{z}}_{Rr2}\right) \left(10\right)$$
The linear accelerations of the points of force application to the tractor body may be computed as following:
$${\ddot{z}}_{afc}={\ddot{z}}_{cg}-{l}_{f}\ddot{\theta } \left(11\right)$$
$${\ddot{z}}_{a{r}_{1}}={\ddot{z}}_{cg}-{l}_{r}\ddot{\theta }+{t}_{1}\ddot{\varphi } \left(12\right)$$
$${\ddot{z}}_{ar2}={\ddot{z}}_{cg}-{l}_{r}\ddot{\theta }-{t}_{1}\ddot{\varphi } \left(13\right)$$
The linear accelerations of the joints between front axle and front wheels are indicated as follows:
$$\ddot{\alpha }=\frac{1}{{I}_{xxa}}\left({F}_{f1}{t}_{1}+{F}_{f2}{t}_{2}\right) \left(14\right)$$
Finally, the mathematical equations for the linear accelerations of the axles of front wheel joints can be computed as follows:
$${\ddot{z}}_{a{f}_{1}}={\ddot{z}}_{afc}+{t}_{1}\ddot{\alpha } \left(15\right)$$
$${\ddot{z}}_{af2}={\ddot{z}}_{afc}+{t}_{2}\ddot{\alpha } \left(16\right)$$