Open access peer-reviewed chapter

Fuzzy Control Strategy for Cooperative Non-holonomic Motion of Cybercars with Passengers Vibration Analysis

By Francesco Maria Raimondi and Maurizio Melluso

Published: January 1st 2010

DOI: 10.5772/6973

Downloaded: 1500

1. Introduction

The cybercars are electric road wheeled non-holonomic vehicles with fully automated driving capabilities. They contribute to sustainable mobility and are employed as passenger vehicles. Non-holonomic mechanics describes the motion of the cybercar constrained by non-integrable constraints, i.e. constraints on the system velocities that do not arise from constraints on the configuration alone. First of all there are thus with dynamic non-holonomic constraints, i.e. constraints preserved by the basic Euler-Lagrange equations (Bloch, 2000, Melluso, 2007, Raimondi & Melluso, 2006a). Of course, these constraints are not externally imposed on the system but rather are consequences of the equations of motion of the cybercar, and so it sometimes convenient to treat them as conservation laws rather than constraints per se. On the other hand, kinematic non-holonomic constraints are those imposed by kinematics, such as rolling constraints. The goal of the motion control of cybercars is to allow the automated vehicle to go from one terminal to another while staying on a defined trajectory and maintaining a set of performance criteria in terms of speeds, accelerations and jerks. There are many results concerning the issue of kinematic motion control for single car (Fierro & Lewis, 1997). The main idea behind the kinematic control algorithms is to define the velocity control inputs which stabilize the closed loop system. These works are based only on the steering kinematics and assume that there exists perfect velocity tracking, i.e. the control signal instantaneously affects the car velocities and this is not true. Other control researchers have target the problems of time varying trajectories tracking, regulating a single car to a desired position/orientation and incorporating the effects of the dynamical model to enhance the overall performance of the closed loop system. The works above are based on a backstepping approach, where the merging of kinematic and dynamic effects leads to the control torques applied to the motors of the wheels. A Fuzzy dynamic closed loop motion control for a single non-holonomic car based on backstepping approach and oriented to stability analysis of the motion errors has been developed by Raimondi & Melluso (2005). In Raimondi & Melluso (2006b) and Raimondi & Melluso (2007a) adaptive fuzzy motion control systems for single non-holonomic automated vehicles with unknown dynamic and kinematic parameters and Kalman’s filter to localize the car have been presented. With regards to the problems of cooperative control of multiple cybercars, a number of techniques have been developed for omni-directional (holonomic) wheeled cars (see Gerkey & Mataric, 2002, La Valle & Hutchinson, 1998). Decentralized algorithms have been dealt with for holonomic cars by Lumelsky & Harinarayan (1997). With regards to the cooperation of multiple non-holonomic cybercars, few results have been published. On this subject, an approach based on the definition of suitable functions of inverse kinematics to control the motion of a platoon of autonomous vehicles has been presented by Antonelli & Chiaverini (2006). The problem of controlling multiple non-holonomic vehicles by using fuzzy control so that they converge to a source has been studied by Driessen et al, (1999). However, since the cars do not have passengers on board, all the studies above do not consider the problem of the acceleration and jerk. For fully automated operation with passengers, a trajectory planning method that produces smooth trajectories with low acceleration is required. The jerk, i.e. the derivative of the acceleration, adversely affects the efficiency of the control algorithms and passengers comfort, so that it has to be reduced. Not many results have been published on this subject (Labakhua et al., 2006, Panfeng et al., 2007).

In this chapter a new closed loop fuzzy control system for non-holonomic motion of multiple cybercars in presence of passengers is proposed. The control strategy merges an innovative decentralized planning trajectory algorithm and a new fuzzy motion control law. About the cooperation, if the target position is fixed, then a number of cybercars has to reach the target one, without to come into collision with the other closest vehicles. The trajectories are planned as the desired time evolution for the position and orientation of some representative point of each cybercar. Forward trajectories are planned only, i.e. trajectories without manoeuvres. In other words all the cooperative cybercars should not stop, except, of course, at the initial and final position. Therefore circular trajectories with continuous curvature have been chosen. Since, for example, in airport the cybercars move in preferential roads without obstacles, the environment in which they move is considered free of obstacle. To ensure the trajectory tracking of all the cooperative cybercars, a new control strategy based on fuzzy inference system is proposed and developed. The fuzzy system generates the control torques for all the cybercars. The cybercars are still employed to transport passengers which are inevitable exposed to vibrations (Birlik & Sezgin, 2007). The acceleration is adopted as preferred measurement of the human vibration exposure. Therefore, with respect to other control theories, the parameters of the fuzzy controller developed in this chapter may be tuned with respect of the ISO 2631-1 standard, which proposes a comfort scale using a mean acceleration index.

This chapter is organized as follows. Section 2 presents the dynamical model of multiple cybercars which has to be employed to project the dynamic fuzzy control system. Also the acceleration model is formulated to develop a control strategy where the passenger comfort is ensured. Section 3 presents a new decentralized cooperative trajectory planner, where the aim is that all the cybercars must reach a target position without collisions between them. Section 4 presents the closed loop fuzzy motion control system, where the asymptotical stability of the motion errors given by the difference between the reference trajectory planned in Section 3 and the actual trajectory of each cybercar is proved by using the Lyapunov’s theorem and the Barbalat’s Lemma (Slotine & Li, 1991). The parameters of the fuzzy control law are investigated to ensure a good level comfort of the passengers. In this sense, the adjustment of the saturation values of the fuzzy dynamic control surfaces guarantees low values of the longitudinal, lateral accelerations and jerks. In Section 5 experimental tests in a Matlab environment are employed to confirm the effectiveness of the proposed motion control strategy. Some conclusions are drawn in Section 6.


2. Model formulation of multiple non-holonomic cybercars

Consider a system made up of rnon-holonomic cybercars. A schematization of the system in open chain configuration is shown in Fig. 1. Now indicating withqi(t)Rnthe time varying coordinates of the position and orientation of the i-cybercar,

Figure 1.

System ofrvehicles in open chain configuration

the complete system is subject tor×mnon-holonomic constraints given by:



ρ(q1(t)q2(t)qr(t)q˙1(t)q˙2(t)q˙r(t))=[A1(q1(t))q˙1(t)A2(q2(t))q˙2(t)      ...Ar(qr(t))q˙r(t)]E2
andAi(qi(t))Rm×nare the matrices of the non-holonomic constraints of the motion of each vehicle. Consider the following well known dynamic model (Fierro & Lewis, 1997) :
Mi(qi(t))q¨i(t)+Ci(qi(t)q˙i(t))q˙i(t)=Ei(qi(t))τi(t)AiT(qi(t))λii=1... rE3

whereMi(qi(t))Rn×nis a positive definite inertial matrix,Ci(qi(t)q˙i(t))Rn×nis the centripetal Coriolis matrix,τi(t)Rpis a vector of the torques applied to the right and left wheels of each vehicle,Ei(qi(t))Rn×pis a mapping between the torques above and the forces applied along the components ofqi(t)Ai(qi(t))are the matrices given by (2), whileλiRmis a vector of Lagrange multipliers. LetSi(qi(t))Rn×(nm)be a full rank matrix made up of a set of smooth and linearly independent vectors spanning the null space ofAi(qi(t)), i.e.:


Due to the non-holonomic constraints (1), it is possible to find velocity vectorsvi(t)Rnm, such that:

q˙i(t)=Si(q i(t))vi(t)i=1... rE5

Referring to the i-vehicle shown in Fig. 2, letP0(x0iy0i)be the reference point of the motion and indicate the orientation withθi. Indicate with(XY)and with(XcYc)the ground and the body references respectively. Letribe the radius of the wheels. Indicate withbithe distance from the center of the wheel to the longitudinal axis and withdithe distance from the reference point to the mass centerPc. For the later description indicate withmiandIithe mass and the inertia of the i-vehicle respectively. In view of the previous considerations, the following model can be written:

q˙i(t)=Si(qi(t))vi(t)Si(qi(t))=(cosθi(t)0sinθi(t)001) viT(t)=[ui(t)  ωi(t)] qiT(t)=[x0i(t)  y0i(t)   θi(t)]  i=1... rE6

whereui(t)andωi(t)are the linear and angular velocities of the i-cybercar of the cooperative system.

Figure 2.

Cybercar with reference systems

To design a control law which consider the ISO 2631-1, and, therefore, to analyze the vibration of the passengers during the motion, it is necessary to obtain the acceleration model. Let the acceleration vector of each cooperative vehicle be:

awiT(t)=[a^xi(t)  a^yi(t)]i=1,...,rE7

The accelerationsa^xi(t)anda^yi(t)are the longitudinal and lateral accelerationsrespectively of the i-cybercar of the cooperative system. The accelerations (7) could be obtained through the following steps. By considering the kinematical model of a non-holonomic vehicle (cf. eq. 6), it is possible to calculate the accelerations along the axes of the ground reference (X,Y) (see Fig. 2) in function of the linear and angular velocities. It results:

axi(t)=d2x0i(t)dt2=x¨0i(t)    ayi(t)=d2y0i(t)dt2=y¨0i(t)E8

Differentiating (6) leads to:

[x¨0i(t)y¨0i(t)θ¨i(t)]=[θ˙i(t)sinθi(t)0θ˙i(t)cosθi(t)000][ui(t)ωi(t)]+               +[cosθi(t)0sinθi(t)001][u˙i(t)ω˙i(t)]  i=1... rE9

Therefore the human body is subjected to forces along the X and Y axes. Now it is necessary to project the forces above along the axes of the body reference, i.e.XcandYc, so that the lateral and longitudinal accelerations (7) can be calculated as:


After some calculations it results:

a^xi(t)= ui(t)θ˙i(t) sinθi(t)cosθi(t)+u˙i(t) cos2θi(t)+            +ui(t)θ˙i(t) cosθi(t) sinθi(t)+u˙isin2θi(t)=u˙i;a^yi=ui(t)θ˙i(t)sinθi(t)cosθi(t)+u˙i(t)cos2θi(t)+        +ui(t)θ˙i(t)cosθi(t)sinθi(t)+u˙i(t)sinθi(t) =uiθ˙iE11

Due to non-holonomic constraints given by (2), the components ofqi(i=1…r) vector are not independent. For this reason the dynamic model (3) cannot be employed directly for the motion control. A dynamic model in body fixed coordinates has to be used. Substituting the equation (6) into model (3) leads to:



M¯i=SiT(qi(t))Mi(qi(t))Si(qi(t))=(mi00Ii) V¯mi(ωi(t))=SiT(qi(t))(Mi(q(t)i)S˙i(qi(t))+Vmi(qi(t))Si(qi(t)))=(0diωimidiωimi0) i=1...rE13

3. Decentralized cooperative trajectory planning algorithm

In this paragraph a new decentralized trajectory planning algorithm is developed for the non-holonomic cooperative motion of the cybercars. The problem is the following. All the cybercars must reach a target position without collisions between them. Once the distances between the initial positions of the cybercars and the target are known, a decentralized algorithm permits to plane circular trajectories intersection-free. After communication to each cybercar of the initial position and target coordinates, the trajectory planner of each cybercar provide to plane a circular trajectory independently of each other. This means we have a decentralized planner. More precisely, indicate with[x0i(0)    y0i(0)   θi(0)]the initial position and orientation of each cybercar. An automated vehicle moving with constant linear and angular velocities, tracks a circular trajectory. Therefore a reference motion may be planned along a circumference that includes the initial coordinates above and the position of the target. With reference to Fig. 3, letCi(xoi(0)y0i(0)θi(0)be the initial position of the i-cybercar. LetB(xTyTθT)the position and orientation of the target. LetΔxiandΔyibe the shiftings along the tangent and radial directions respectively of the i-vehicle. The i-cybercar moves from Ci to Balong a circular trajectory of which the radius is equal toλi, while the center is located in D.

Figure 3.

Motion planning

It yields:

BCi=d˜i=(xTx0i(0))2+(yTy0i(0))2γi=arctan [(yTy0i(0))/(xTx0i(0))] i=1...rE14

whereγiis the angle betweenBCiand the xaxis. Consider the following angular relation (see Fig.3):

δi=γiθi(0)==arctan [(yTy0i(0))/(xTx0i(0))]θi(0)E15

The length of the line BA is equal to the distanceΔxi. Therefore, if we consider the triangleCiAB, then:

Δxi=d˜icosδi;  Δyi=d˜isinδiE16

The angular shiftingαibetweenCiand B results as it follows:


From observation of the triangle DAB, it results:


The solution of the equation (13) with respect toλiis:


Now the values of the reference angular (ωri) and linear (uri) velocities of each cooperative vehicle may be calculated as follows:

ωri=αiΔT;  uri=ωiλi;i=1...rE20

whereΔTis a fixed look-ahead time interval chosen by the designer. Let us consider multiple automated cybercars in an initial open chain configuration (cf. Fig. 1), i.e. collinear and with the same orientations given byΔxi(i=1,2). The algorithm above allows circular trajectories without intersections to be planned, so that the vehicles will avoid collisions while moving. Each trajectory is planned independently of the others. This means we have a decentralized cooperation of the vehicles. Fig. 4 shows an example, where two vehicles are considered in open chain configurationC1C2. One observes that the first vehicle of the open chain follows a circular trajectory fromC1to the target BalongΔx1, while the second vehicle follows a circular trajectory fromC2to the target one alongΔx2. The distance betweenC1and Bis smaller than the distance betweenC2and B, so that, based on the equation (9), it isd˜1d˜2. Consequently, based on the equations (16) and (19), the radius of the circumference tracked by the first vehicle (i.e.λ1) is smaller than the radius of the circumference tracked by the second vehicle (i.e.λ2). Since the cybercars are initially collinear and have the same orientations, and the circumferences must be include both the initial positions of the vehicles and the target position, the trajectories are without intersections and the vehicles can reach the target without collisions. The method can be used for r cooperative vehicles in initial open chainC1C2...Cr, so that each vehicle can reach the target without coming into collision with other vehicles.

Note that the vehicles have to be in open chain configuration initially, i.e. collinear. If there is a vehicle which is not mutually collinear, it must reach a collinear position. On this subject, some studies have focused on modelling formations of non-holonomic vehicles (Bicho & Monteiro, 2003).

Figure 4.

Trajectory planning without collisions for multiple vehicles in initial open chain configuration

4. Fuzzy dynamic closed loop motion control for cooperative cybercars with passengers comfort

Consider the i-cybercar of the cooperative system (cf. Fig. 2).The kinematical model is given by (6), while the dynamical model is given by (12) and (13). Employing the values of linear and angular velocities given by (20) and using the kinematical model (6) lead to the following equations for the circular reference motion of each cybercar:

x˙ri(t)=uricosθri(t);y˙ri(t)=urisinθri(t);θ˙ri(t)=ωri   i=1...rE21

Let the following vectors:

qiT(t)=[x0i(t)  y0i(t)  θi(t)] i=1...rE22

be the position and orientation of each cybercar. One defines the following motion errors between the planned circular reference trajectories the state variables given by (22) as it follows:

ei(t)=[cosθi(t)sinθi(t)0sinθi(t)cosθi(t)0001]×[xri(t)x0i(t)yri(t)y0i(t)θri(t)θi(t)]=[exi(t)eyi(t)eθ i(t)]      i=1...rE23

The errorsexi(t)andeyi(t)are said to be the longitudinal and lateral motion errors respectively, whileeθi(t)is the orientation error. The following velocity control laws are employed for each cybercar:

vc i(t)=[uci(t)ωci(t)]=[uricoseθi(t)+K1iexi(t)ωri+uri(K2ieyi(t)+K3isineθ i(t))] K1iK2iK3i0i=1... rE24

Replacing (24) into (6) leads to the following closed loop mathematical model:

e˙i(t)=[(ωri+uri(K2ieyi(t)+K3isineθ i(t)))eyi(t)K1iexi(t)(ωri+uri(K2ieyi(t)+K3isineθ i(t)))exi(t)+uisineθ i(t)uri(K2ieyi(t)+K3isineθ i(t))] i=1... rE25

It is possible to formulate the following theorem.

Theorem 1. Consider the cooperative system (6), in closed loop with the velocity control laws (24). If the linear and angular reference velocities given by (21) are limited functions, then the equilibrium state of the closed loop model (25) is the origin of the state space and it is asymptotically stable.

Proof. From model (25), it is evident that, ife˙i(t)=0(i=1…r), thenei(t)=0(i=1…r),so that the equilibrium point is the origin of the state space. Consider the following extended vectore(t)R(r×n)×1which contains the motion errors of all the cooperative vehicles:

eT(t)=[e1(t)  e2(t)    .    .    .   er(t)] E26

The following Lyapunov’s function is chosen:

V(e(t))=12i=1rK1i(exi2(t)+eyi2(t))+2i=1rK1iK2i(1coseθ i(t)) E27

The function (27) is definite positive. By calculating the time derivative of the function (27) and substituting the equations (25) into result, it yields:

V˙(e(t))=i=1rK1iexi2(t)i=1ruri2K1iK3iK2isin2eθ i(t) E28

The function (28) does not depend on lateral motion errorseyi(t)(i=1…r), so that it is equal to zero in correspondence of the inputs[0   eyi(t)  0]. Therefore the function (28) is semi-definite negative. The conclusion is that the components of the vector (26) are stable and bounded. Since all the motion errors and the reference velocities are bounded, the second time derivative of the function (28) is bounded, therefore Barbalat’s lemma implies that the function (28) converges to zero when t diverges, so that the longitudinal motion errors and the orientation errors of all the cooperative vehicles converge to zero. From the second and third equations of system (20) it follows that:

limte˙yi(t)=0i=1... rE29


eyi()=e¯yii=1... rE30

wheree¯yiis a constant value. Since the orientation errors converge to zero, from the third equation of the system (25) it results:

uiK2ie¯yi=0i=1... rE31

It can be concluded that the lateral motion errors converge asymptotically to zero. Q.E.D.

By employing the kinematical control strategy (24), it is difficult to control directly the lateral and longitudinal accelerations which are responsible of harmful effects on the passengers. For this reason a fuzzy dynamical control strategy is developed below, where the properties of the fuzzy maps assures the Lyapunov’s stability of the motion errors given by (25), while the saturation properties of the maps ones permit to control directly the maximum acceleration of each vehicle of the cooperative system during the motion. Letη˜i(t)the time varying error between the velocity control laws given by (24) and the physical velocityvi(t)of each cybercar (i.e. the solution of the differential equations (12)):

η˜i(t)=[η˜1i(t) η˜2i(t)]=[vci(t)vi(t)] i=1....rE32

The fuzzy inference mechanism is explained below. The inputs of the fuzzy system are the errors (32). The fuzzy rules forη˜ji(j=1,2and i=1,2…r) are the following:

1) ifη˜jiis negative andη˜˙jiis negative then the outputΣjihas a negative big value;

2) ifη˜jiis negative andη˜˙jiis positive, then the outputΣjihas a negative small value;

3) ifη˜jiis positive andη˜˙jiis negative, then the outputΣjihas a positive small value;

4) ifη˜jiis positive andη˜˙jiis positive, then the outputΣjihas a positive big value.

Now we assume the following dynamical control laws (Raimondi & Melluso, 2007b):

SiT(qi(t))Ei(qi(t))τi(t)==M¯i(v˙ci(t)+s˙i(t))+V¯mi(ωi(t))vi(t)=[Fi(t)Ti(t)] i=i... r  E33

whereSi(qi(t))M¯iandV¯mi(ωi(t))are given by (6) and (13) respectively,Fi(t)andTi(t)are the linear force and the angular momentum applied to the i-cybercar, whiles˙i(t)is the output of the fuzzy inference system, so that it is:

 s˙i(t)=Σi(η˜ i(t)η˜˙ i(t))=[Σ1i(η˜1i(t) η˜˙1i(t)) Σ2i(η˜2i(t) η˜˙2i(t))] i=1...rE34

whereΣji(j=1,2i=1…r) are the nonlinear input-output maps of the fuzzy inference system. The maps above depends on the choice of the input and output memberships. Figs 5 and 6 show the membership functions for obtainingΣ1i(i=1…r).

Figure 5.

Input membership functions

Figure 6.

Output membership function

The memberships employed for obtaining the control surfacesΣ2i(i=1…r)have the same symmetry. The input values of the memberships of Fig. 5 are chosen in such a way to be in the range of the physical velocity errors obtained in our application. The outputsΣji(j=1,2i=1,2,…r)given by (34) represent the variations of the torques applied to the right and left wheels of each cooperative cybercar for positive or negative changes of the speed errors given by (32). Due to the choice of the form of the memberships functions, the outputs above are nonlinear uniformly continuous and bounded functions. They contribute to evaluating the acceleration or deceleration of the cybercars, to track the planned reference motion with good dynamic perfomances.

Assumption 1.The following properties for the nonlinear fuzzy maps given by (33) are assumed:


The assumption above is a necessary condition to ensure that each cooperative cybercar can follow the planned trajectory; also by varying the valuesδij, the acceleration saturation values can be changed, so that the motion of each cybercar can be controlled with low accelerations. Note that, if we choice different number or form of the membership functions, then the properties given by (35) cannot be guaranteed.

Substituting (33) into dynamical model (12) leads to:

η˜˙i(t)+Σi(η˜i(t)η˜˙i(t))=0i=1... rE36

It can be written:

Σ¯i(η˜i(t)η˜˙i(t))=0i=1... rE37

so that:

η˜˙i(t)=fi(η˜i(t))i=1... rE38

Considering equations (25) and (38) leads to the following closed loop model of the fuzzy dynamic control system:

e˙i(t)=[e˙xi(t)e˙yi(t)e˙θ i(t)]=[(ωri+uri(K2ieyi(t)+K3isineθi(t)))eyi(t)K1iex i(t)(ωri+uri(K2ieyi(t)+K3isin(eθ i)))exi(t)+urisineθ i(t)uri(K2ieyi(t)+K3isineθ i(t))] η˜˙ i(t)=fi(η˜ i(t)) i=1... rE39

Now the following theorem may be formulated.

Theorem 2.Consider the closed loop system given by (39). Under the assumption 1, the equilibrium state of the modeled system is the origin of the state space and it is asymptotically stable.

Proof.Consider the following extended state vector:

e¯(t)=[e¯1(t)e¯2(t)e¯r(t)] E40


e¯i(t)=[exi(t)eyi(t)eθi(t)η˜1i(t)η˜2i(t)]=[ei(t)η˜i(t)] i=1... rE41

From equations (39) it results:

[e˙xi(t) e˙yi(t) e˙θi(t)] T=0   t;   i=1...rE42

so that:

[exi(t) eyi (t)eθi(t)] T=0t; i=1...rE43

Also, ifη˜˙i(t)=0, then:

Σi(η˜ i(t)0)=0i=1...rE44

so that, based on the property 4 given in the Assumption 1, it can be said that:

η˜ i(t)=0t,   i=1...rE45

This implies that:


Therefore the equilibrium point of model (39) is the origin of the state space. The following Lyapunov’s function is chosen:

V(e¯(t))=12i=1rK1i(exi2(t)+eyi2(t))+2i=1rK1iK2i(1coseθi(t))            +12i=1rη˜iT(t)η˜i(t)++j=12i=1r0η˜jiΣji(ξji,0)dξji K1iK2i0E47

The first, second and third terms of function (47) are always positive. Now, from the property 6 of the fuzzy maps it results:

0η˜ji(t)(Σji(ηji(t)η˜˙ji(t))Σji(0,η˜˙ji(t)))j=1,2 ;   i=1,2,...rE48

Consequently, ifη˙ji(t)=0, then it yields:


so that, based on the property 1, one obtains:

0η˜ji(t)Σji(ηji(t),0)j=1,2;   i=1,2...rE50

It can be concluded that:

j=12i=1r0η˜jiΣji(ξji,0)dξji0    η˜jij=1,2;   i=1,2...rE51

Therefore the function (47) is positive definite. Calculating the time derivative of the function (47) and replacing (39) into it lead to:

V˙(e¯)=i=1rK1iexi2(t)i=1ruri2K1iK3iK2isin2(eθi(t))i=1rη˜iT(t)Σi(η˜i(t)η˜˙i(t))+          i=1rΣiT(η˜i(t)η˜˙i(t))Σi(η˜i(t)0)  E52

The first and second terms of (52) are negative. Based on the property 5 (see Assumption 1), the elements of the summation of the third term of (52) are positive numbers, so that the term above is negative. From property 5 and inequality (50) it yields:

η˜ji2(t)Σji(η˜ji(t),0)Σji(η˜ji(t)η˜˙ji(t))0j=1,2;   i=1,2,...,rE53

so that, ifΣji(η˜ji(t),0)is positive, thenΣji(η˜ji(t)η˜˙ji(t))is positive. Therefore the elements of the summation of the fourth term of function (52) are positive. Note that function (52) does not depend on the lateral error, therefore it can be concluded that the function above is negative semi-definite. Therefore the components of the vector (40) are bounded. Now it is possible to calculate the second time derivative of function (47). Based on the previous considerations, it is a bounded function, so that, by applying the Barbalat’s lemma, it follows that:


From (52) it can be concluded that the errorsexieθiη˜1iη˜2i(i=1,…,r) converge to zero. Replacing (24) into (32) leads to:

η˜i(t)=[uci(t)ui(t)ωci(t)ωi(t)]=[η˜1i(t)η˜2i(t)]==[uricos(eθi(t))+K1iexi(t)ui(t)ωri+K2iurieyi(t)+K3iurisineθ i(t)ωi(t)]  i=1...rE55

Since the errorsexieθ iη˜1iη˜2iconverge asymptotically to zero, the lateral errorseyiof each cooperative cybercar converges asymptotically to zero. Q.E.D.

Fig. 7 illustrates the block scheme of the Fuzzy dynamical motion closed loop control system for a single cybercar.

With regards to the passengers comfort, several factor influence vibration discomfort in relation to passenger activities, e.g. seated posture, use of backrest. Passengers usually adopt their posture to attenuate the intensity of vibrations and jerks in order to perform their activities satisfactorily. However the transmission of vibrations on the human body is higher if a passenger uses armrest, backrest or places boot feet on the floor. Therefore attenuation of vibration exposure is a very important requirement of a motion control system for cybercars. There are various means by which the vibration may be expressed, such as

Figure 7.

Block scheme of the dynamical fuzzy control system for cybercar motion

displacement, velocity and acceleration. Of these physical quantities acceleration is generally adopted as preferred measured of quantifying the severity of human vibration exposure (Suzuki, 1998). From (10) it appears that the accelerations which cause vibrations on the human body depend on the curvature of the trajectory. Indicate withawiR2the following vector:

awiT=[a^xi(t)  a^yi(t)]     i=1... rE56

wherea^xi(t)anda^yi(t)are the accelerations given by (11). Indicate witha¯xianda¯yithe r.m.s. values of the accelerations given by (56). Indicate the jerks with the following vector:

j iT(t)=[jxi(t)  jyi(t)]i=1...r E57

The componentsjxi(t)andjyi(t)are said to be the lateral and longitudinal jerksof each cybercar, i.e. the rate of change of the accelerations as it follows:

jxi(t)=da^xi(t)dt  jyi(t)=da^yi(t)dt i=1...r E58

Indicate withj¯xiandj¯yithe R.M.S. values of the jerks given by (58). The ISO 2631-1 Standard relates comfort of the passengers with the r.m.s. overall acceleration which causes vibrations acting on the human body defined as:


wherea¯xianda¯yiare given by (56)a¯ziis the acceleration on the z axis perpendicular to plane Xc,Yc (see Fig. 2), whileαβandγare multiplying factors. Since each vehicle moves on the plane, the accelerationa¯ziis equal to zero. Based on the ISO 2631-1, the relations between the values of the overall acceleration given by (58) and the passenger comfort of the i-cybercar are given by the Table 1.

R.M.S. overall accelerationPassenger comfort level
a˜wi0.315m/s2Not uncomfortable
0.315a˜wi0.63m/s2A little uncomfortable
0.5a˜wi1m/s2Fairy uncomfortable
1.25a˜wi2.5m/s2Very uncomfortable
a˜wi2.5m/s2Extremely uncomfortable

Table 1.

ISO 2631-1 Standard

By using our fuzzy approach, it is possible to obtain low values of the lateral and longitudinal accelerations in easy way. In fact the saturation values of the outputs of the fuzzy maps, i.e. the valuesδjigiven by the third property of (35), represent a saturation of the linear and angular accelerations of the cooperative cybercars, so that they are degree of freedom and can be chosen by the designer. In particular, from equation (9) it can be seen that the accelerations along the axes of the ground reference depend on the linear accelerationu˙i(t)and on the angular accelerationω˙i(t). Now, during the motion, the values ofu˙i(t)andω˙i(t)depend on the outputs of the dynamical fuzzy control law given by (33) and (34). Fig. 8 shows the typical fuzzy control surfaces obtained by using the fuzzy control law (33), where the acceleration given by (34) satisfies the properties of the assumption 1.

Figure 8.

Fuzzy mapping

The saturation property of the fuzzy map causes the accelerations above to be bounded, and so, after a few attempts, the designer can choose the ranges of the crisp output values of the fuzzy system in order to satisfy the ISO 2631-1 standard and to reduce the vibrations and the jerks. It is evident that the the output values of Fig. 8 fall within a small range, since a motion with low acceleration and jerk is desired. The simulation experiments described in the next section confirm the efficiency of our algorithm in terms of cooperation, stability of the fuzzy motion control system and passenger comfort.

5. Simulation experiments

In this performed simulations the efficiency of the cooperative fuzzy motion control law proposed and developed in this chapter and the good level comfort of the passengers during the motion of the cybercar is illustrated. The parameters of the cybercars have been chosen based on existing cybercars (McDonald & Voge, 2003). The weight of a cybercar is 300kg, the width is 1.45m, the height is 1.6m, while the length is 3.7m. Referring to Fig. 2, the kinematical parameters are chosen as:

ri=0.40m;  bi=0.725m;  i=1,2,3E60

The dynamical parameters are:

Mi=300kg;  di=1.5m;  i=1,2,3E61

The parameters of the speed control law (24) are given by:


The reference trajectories of each cybercar were generated using the decentralized algorithm developed in Section 3, so that the initial motion error values are equal to zero. In fact the circumferences have been generated based on the distance between the initial position of the cybercars and the position of the target. Initially the vehicles are in open chain configuration along y-direction. The initial positions of the three cooperative cybercars are the following:

x01(t=0)=2m;  x02(t=0)=2m ;  x03(t=0)=2m;y01(t=0)=2m;y02(t=0)=3m;y03(t=0)=4m;θ01(t=0)=1,74rad;θ02(t=0)=1.74rad; θ03(t=0)=1.74radE63

All the generalized coordinates given by (63) refer to a ground reference whose origin is shown in Fig. 9. The position coordinates of the target with respect to the ground reference are:

[xTyT]=[1.7m2.8m]  E64

To analyze the performances in terms of passenger comforts, we compare for cases with reference to the parametersδij given by the third property of (35) :

a. low values of the parametersδji  (i=1,2; j=1,2,3);

b. high values of the parametersδji  (i=1,2; j=1,2,3)

Case a- The absolute saturation values of the crisp outputs of the fuzzy control surfaces are:

δ1i=|Σ1i_sat|= 0.1m/s2;  δ2i=|Σ2i_sat|=0.1rad/s2;i=1,2,3E65

Fig. 9 shows the planar trajectories of the cybercars as planned by using the algorithm given in Section 3.

Initially the cybercars are in open chain configuration. The trajectories are intersections-free and therefore there are not collisions during the motion.

The graphs of Fig. 10 and 11 show the time evolutions of the velocity errors given by (32). Due to the dynamics of the cybercars, there are not perfect velocity tracking, i.e. the speed control laws (24) do not affect instantaneously the linear and angular velocities, but the errors converge to zero after some times.

The most significant graphs which illustrate the stability performances of the motion errors of each cybercar are drawn in Figs. 122 and 13, where the time evolutions of the longitudinal, lateral and orientation errors given by (23) are shown.

Figure 9.

Motion of three cooperative cybercars

Figure 10.

Velocity errors of the cybercars 1 and 2

Figure 11.

Velocity errors of the cybercars 3

Figure 12.

Tracking errors of cybercars 1 and 2.

Figure 13.

Tracking errors of cybercar 3

Now we investigate on the passenger comforts with the saturation values given by (65). Figs. 144 and 15 shows the time evolution of the accelerations given by (56) which are responsible of vibrations on the human body, while Table 2 illustrates the r.m.s values of the same accelerations and the overall acceleration given by the mean index (59).

Figure 14.

Lateral and longitudinal accelerations of the cybercars 1 and 2.

Figure 15.

Lateral and longitudinal accelerations of the cybercars 3.

Longitudinal accelerationLateral accelerationOverall Acceleration

Table 2.

r.m.s. and mean accelerations of all the cybercars

Note that the values of the r.m.s. overall accelerations are between “Not uncomfortable” and “A little uncomfortable” (see Table 1 and 2), so that the passengers comfort level is very good.

Case b- With respect to (65), the saturation values of the fuzzy control surfaces are increased as follows:

δ1i=|Σ1i_sat|= 0.5m/s2;  δ2i=|Σ2i_sat|=0.5rad/s2;i=1,2,3E66

The comfort of the passengers are also studied in this case. On this subject Figs. 16 and 17 illustrate the lateral and longitudinal accelerations of the cybercars. The r.m.s values of the accelerations above and the mean acceleration given by the index (59) are listed in Table 3, while the values of the jerks given by (58) in cases of low and high saturation values of the fuzzy control surfaces are illustrated in table 4.

Figures 16 and 17 and the results of the table 3 show that the overall accelerations are between “fairy uncomfortable” and “uncomfortable”, so that the comfort of the passengers during the motion is bad. By the results shown in Tables 2-4 it is evident that, in case of low values saturation of the fuzzy control surfaces, the accelerations and the jerks are reduced, which means ride passengers comfort enhancement. Therefore the designer can be choice the parameters of the fuzzy controller to optimize the vibrations acting on the human body of the passengers.

Figure 16.

Lateral and longitudinal accelerations of the cybercars 1 and 2

Figure 17.

Lateral and longitudinal accelerations of the cybercar 3

Longitudinal accelerationLateral accelerationOverall Acceleration

Table 3.

r.m.s. and mean accelerations of all the cybercars

Jerks with low saturation values of the fuzzy control surfaces (cfr. (65)).Jerks with high saturation values of the fuzzy control surfaces (cfr. (66)).

Table 4.

r.m.s. lateral and longitudinal jerks of all the cybercars


6. Conclusion

In this chapter a new fuzzy cooperative control algorithm for multiple fully automated cybercars, where the parameters of fuzzy controller may be tuned to obtain low vibrations on the body of the passengers, has been developed. A generalized mathematical model for multiple cybercars to project the fuzzy control system and an acceleration model to ensure the comfort of the passengers have been formulated. A new decentralized trajectory planner which guarantee the absence of collisions between the closest vehicles has been presented. A new fuzzy control strategy which stabilizes all the cooperative vehicles in the planned trajectories has been developed, where the asymptotical stability of the motion errors has been proved by using Lyapunov’s theorem and Barbalat’s lemma. Good passengers comfort levels during the motion has been ensured by tuning of the saturation of the fuzzy maps. In the simulation tests an example in case of motion control of three automated vehicles has been developed. Trajectories without intersections have been generated and, by choosing suitable input-output values of the fuzzy maps, the stability of the motion errors and very good passengers comfort levels based on ISO 2631-1 Standard have been obtained.


This work was realized with the contribution of the MIUR ex-60%.

All sections have been equally and jointly developed by the authors.

© 2010 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Francesco Maria Raimondi and Maurizio Melluso (January 1st 2010). Fuzzy Control Strategy for Cooperative Non-holonomic Motion of Cybercars with Passengers Vibration Analysis, Motion Control, Federico Casolo, IntechOpen, DOI: 10.5772/6973. Available from:

chapter statistics

1500total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Characteristics of Mechanical Noise during Motion Control Applications

By Mehmet Emin Yüksekkaya

Related Book

First chapter

Introduction to Infrared Spectroscopy

By Theophile Theophanides

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us