4.1 Tire Lateral Force 56Figure 4.1: A scheme <strong>of</strong> a spring-mass-damper system with friction between the mass and the ground.∫ t∆Y = Y ic − Y i = y ic = − v iy (τ)dτ (4.1)0where y ic denotes the position <strong>of</strong> the contact point with respect to the local frame.The mass does not move until the resultant <strong>of</strong> all the <strong>for</strong>ces applied to the contact point (F ic )does not overcome the maximum static friction <strong>for</strong>ce F smax= µ sy N i , i.e. until the condition|F ic | ≤ F smax is satisfied, with N i representing the normal reaction to the ground. Conversely,when such a condition is not anymore satisfied, the mass starts moving and the contact pointposition is described by the following equation:and relation (4.1) becomes:m t ÿ ic + D y ẏ ic + K y y ic = F ky = −µ ky N i sign(v ic ) (4.2)∫ t∆Y = y ic = (v ic (τ) − v iy (τ))dτ (4.3)t kwhere v ic = Ẏ ic represents the velocity <strong>of</strong> the contact point with respect to the fixed frame,expressed in the local frame.It is worth noticing that relation (4.3) still holds if from t k the mass is not moving, i.e. v ic =0. By considering v iy as input <strong>of</strong> the system and substituting (4.3) in (4.2), we obtain thefollowing equation:m t (˙v ic (t) − ˙v iy (t)) + D y (v ic (t) − v iy (t)) + K y∫(v ic (t) − v iy (t))dt = −µ ky N i sign(v ic ) (4.4)Let us suppose that the contact point starts moving with v ic > 0, so that sign(v ic ) = 1. Bysolving equation (4.4) in the frequency domain, we obtain:
4.1 Tire Lateral Force 57µ ky NV ic (s) = V iy (s) −(4.5)m t s 2 + D y s + K ywhere V ic (s) and V iy (s) are the Laplace trans<strong>for</strong>m <strong>of</strong> v ic and v iy respectively.By calculating the Laplace inverse trans<strong>for</strong>m <strong>of</strong> equation (4.5), we obtain the following relationin the time domain:where h(t) = L −1 1m t s 2 +D y s+K y.v ic (t) = v iy (t) − µ ky N i h(t) (4.6)We notice that, if the parameters m t , D y , K y are defined such that the system is stable, thecontribution <strong>of</strong> the friction in (4.6) initially increases, reaching its maximum, and then goesexponentially to zero, either oscillating or directly, depending on whether the system is underor over-damped. Thereby, if we suppose that v ic (t) > 0 ∀t, i.e. v iy (t) > µ ky N i h(t) ∀t,the contact point will keep moving with a velocity v ic (t) → v iy (t). Conversely, if we havev ic (t) < 0 <strong>for</strong> some values <strong>of</strong> t, equation (4.6) does not hold anymore, as the friction <strong>for</strong>cechanges its sign, and the nonlinearity <strong>of</strong> the sign() function can affect the system stability.However, it can be proved, by substituting the function sign(v ic ) with arctan(k s v ic ), k s ¨big¨,that the system is at least simply stable [21]. In particular, if |v ic |, |˙v ic | become relativelysmall, i.e. |v ic |, |˙v ic | < µ kyN ik s, we can consider as if the contact point stops. Thereby, the systemsuddenly switches to the static configuration, in which the friction coefficient switches fromthe kinetic to the static one, so that the contact point does not move until the elastic <strong>for</strong>cedoes not overcome the static friction <strong>for</strong>ce, i.e. until the condition |F ic | ≤ F smax is satisfied.The system describing the dynamics <strong>of</strong> the contact point presents, there<strong>for</strong>e, a sort <strong>of</strong> hysteresis.More precisely, when the condition |F ic | ≤ F smaxis not anymore satisfied the dynamicsswitches from v ic = 0 to (4.4), and when the condition |v ic |, |˙v ic | < µ kyN ik sfrom (4.4) to v ic = 0, and so on.is satisfied it switchesBy considering the local frame in Figure 4.1 as the frame fixed to the i th wheel and takinginto account the a<strong>for</strong>ementioned considerations, the lateral reaction <strong>for</strong>ce perceived from thei th wheel can be written as:f iy = −D y (v iy (t) − v ic (t)) − K y∫(v iy (t) − v ic (t))dtwhere the dynamics <strong>of</strong> v ic is either v ic = 0 or (4.4), depending whether the i th contact pointis moving or not, and the wheel center velocity v iy is related to the robot linear and angular