6.1 Simulink Model 84(a)(b)Figure 6.4: (a) Block scheme <strong>for</strong> the longitudinal reaction <strong>for</strong>ce; (b) Block scheme <strong>for</strong> the contact point dynamics.
6.1 Simulink Model 85as explained in Section 4.1. Such a dynamics <strong>of</strong> the contact point when the tire is releasingis so complex that it would require an ¨ad hoc¨ study, which is not provided in this thesis.In order to overcome this issue, we consider N ic = m t g and we assume that the conditionson F icx given in (4.13) still hold. However, when the value <strong>of</strong> F icx switches from F icx 0 toF icx = 0, the value <strong>for</strong> v icx could not be zero as it should be. In such a case, as the contactpoint dynamics switches from m t˙v icx 0 to m t˙v icx = 0, the velocity v icx will stay stable tothe value it had be<strong>for</strong>e switching, so that also the tire reaction <strong>for</strong>ce will keep its stable valueinstead <strong>of</strong> per<strong>for</strong>ming the expected stretching/releasing behavior. The solution to this issuecan be found by simply rewriting (4.13) as follows:f ix = −D x (∆v ix − v icx ) − K x∫(∆v ix − v icx )dtm t˙v icx + D x v icx + µ kx N ci sign(v icx ) = F icx′⎧⎪⎨ 0, <strong>for</strong> |FF icx ′ ix ′=| < µ sxN i ∧ f i (t − ∆t) = 0⎪⎩ Fix ′ , <strong>for</strong> |F′ ix | > µ kxN i ∧ f i (t − ∆t) 0∫F ix ′ = m t˙v ix + D x v ix + K x (v ix − v icx )dt(6.4)where f i = 1, 0 is a flag indicating whether the contact point is moving or not.We notice that, in this case, when the value <strong>of</strong> F icx switches from F ′ icx 0 to F′ icx= 0, thesecond equation in (6.4) represents an asymptotically stable system instead <strong>of</strong> a simply stablesystem, there<strong>for</strong>e v icx will tend to zero independently from the value it had be<strong>for</strong>e switching.Briefly, as we do not know the dynamics <strong>of</strong> N ci allowing the contact point to stop, we areimposing v icx = 0 when F ′ icx= 0. However, it must be noticed that the overall dynamicsdoes not change from the theoretical one, except <strong>for</strong> the fact that the contact point does notstop when Ficx ′ = 0 as <strong>for</strong> the theoretical case, but it stops some time later. Moreover, as weit will be shown in the next Section, the system represented by the second equation in (6.4)is over-damped, i.e. D x > 2 √ m t K x , and the corresponding time constant is relatively small,i.e.√mtK x≪ 0.1. As a consequence, the contact point stops immediately after F ′ icx= 0. Itis worth noticing also that such approximation <strong>of</strong> the contact point dynamics does not takeinto account the case in which v ix is relatively high, or N ci relatively low, so that v icx couldnever reach zero, as explained in Section 4.1. However, we are considering only the caseω z < 64 degs≈ 1 rad,that is v s ix < 0.15 m , and the acquired data from the <strong>for</strong>ce sensor alwaysspresents the saw-tooth shape when pulling the robot up to 0.1 m , there<strong>for</strong>e we can assumes