ap-6-10.dvi Acta Polytechnica Vol. 50 No. 6/2010 A Model of Active Roll Vehicle Suspension I. Čech Abstract This paper describes active suspension with active roll for four-wheel vehicle (bus) by means of an in-series pump actuator with doubled hydropneumatic springs. It also gives full control law with no sky-craping. Lateral stiffness and solid axle geometry in the mechanical model are not neglected. Responses to lateral input as well as responses to statistical unevennesses show considerable improvement of passengers comfort and safety when cornering. Keywords: active suspension, roll-yaw model of a four-wheel vehicle, cross control, solid axle, hydraulic control, pump actuator. 1 Introduction A mathematical description of a model of active roll vehicle suspension is made in steady state of the har- monic input in symbolic form. The complex am- plitudes and effective values are denoted by capi- tal letters and the instantaneous values, and also the constant values, are denoted in lower-case let- ters. Differentiations are substituted by the operator s = i2πf , second differentiations by s2, and so on. As our model includes two more stages of differenti- ations than in usual mechanical models, a complex symbolic form is necessary to handle the model. Only linear relations are used. A complete 4-wheel vehicle is dealt with here, but no heave input is presumed. The parameters relate to a tall bus. The computed answers are compared with passive suspension. The results are given for both vertical and hori- zontal input, in spectral form and in impulse form. Unlike other proposed solutions, e.g. [1, 2, 3] our project includes the static load control, and should be low-powered. This is because of the in-series char- acter of the control, so that static displacement com- pensation takes no power. In addition, the controller uses no throttling and works against no static load. The control scheme is simple, so that there are no problems with stability. There is also the possibility of achieving active suspension using a source of force parallel to air sus- pension (which includes the static load displacement control). This source of force can be implemented by a linear electric motor. However, the power pro- duced by the linear motor is small in comparison with a rotational electric motor with gearing. 2 Mechanical model Dealing only with linear relations, we first presume the unevennesses decomposition diagram in Fig. 1. It shows an analysis of the vertical input (unevennesses) z1–z4 of the four wheels 1–4 of the model into input zks1 , zks2 of the heave-pitch model [4], which will not be dealt with here, and input zkp1 , zkp2 of the four- wheel roll model. This input will be marked zk1 , zk2 . Indices 1, 2 relate to front and rear axles. The me- chanical schema of this roll-yaw model is shown in Fig. 2. Fig. 1: Analysis of vertical inputs This figure shows a model of a suspension with a solid rear axle. We will consider the radii of gyra- tion of the body rbx, of the seat rsx, and of the axle rwx. (The mass of the axle is assumed to be concen- trated in the wheels.) Then there are the heights of the mass centers of the body ztb, of the seat zts and of the wheel ztw. The track is denoted by yw, the lateral distance of the seats by ys, and the distance of the spring settings by yc. The height of the joint of the solid axle is marked by zq, and the lateral stiff- ness is marked by ky. The lateral displacements of the wheels at stiffnesses ky1 , ky2 are denoted by Yky1 , Yky2 . Vertical antiphased displacements of the body (above the wheel), seat and wheel are noted by Zb, Zs, and Zw1 , Zw2 . 16 Acta Polytechnica Vol. 50 No. 6/2010 Fig. 2: Mechanical scheme The most important input of this model is the lateral acceleration Ac. We will take into account the lateral sliding of the rolling tyre. Sliding acts similarly to damper. Its damping is proportional to the weight of the vehicle and reciprocally proportional to the travelling veloc- ity and the constant of the tyre tgα. The stiffness of a damper is its damping multiplied by operator s. So kypk = s · ag(mb +ms +mw)/vx/ tgα The lateral stiffness ky consists of this sliding stiff- ness kypk and the lateral stiffness of the tyre kwy act- ing in series. So 1/ky = 1/kwy + 1/kypk, We will use kwyre = 0.6kw, tgα = 0.1. The lateral input of the model is due to the lateral acceleration Ac of the body. The solid axle rotates around the center of the anti-phased unevennesses C, the body rotates around joint Q and translates in lateral direction with the lateral displacement of the joint. Only small angular displacements are assumed. The lateral stiffness in front and rear, and also the wheel masses front and rear correspond to the position of the mass centre, i.e. ky1/ky2 = mw1/mw2 = xw2/xw1 The roll model rotates around the on-axis, but its angle α α = arctg(zq/xw) is small and approximately taken as cosα = 1. Nev- ertheless, the efficient value of the axle joint height under the mass centre of the body is zqt = xw1/xw · zq 17 Acta Polytechnica Vol. 50 No. 6/2010 3 The equations of the mechanical model Two auxiliary constants are used to enable the equa- tions to be used also for independent suspension, namely tn = 1, nz = 0 when with solid axle, tn = 0, nz = 1 when no solid axle is used. (With indepen- dent suspension the inertia forces of the wheel are supported by the body.) If Y z1, Y z2 are the lateral displacements of the yaw motion, then the lateral force originating from the body mass mb is due to the acceleration of cor- nering, the roll movement of the body, the roll of the joint, and the lateral shift of the mass center, so Sqb=4mbAc+ 4mbs2[(ztb − zqt)Zb + zqtZw2 ] · 2/yw− 4mb · s2 · (xw2Yky1 + xw1Yky2)/xw (3.1) and the lateral force of the seat Sqs=4msAc+ 4mss2[(zts − zqt)Zb + zqtZw2 ] · 2/yw− 4ms · s2 · (xw2Yky1 + xw1Yky2)/xw (3.2) The vertical forces in the suspension between body and wheel, front and rear Sbw1 = kb1(Zb − Zw!− Zr1) + (3.3) kst1(Zb − Zw!), Sbw2 = kb2(yc/yw · (Zb − Zw2)− Zr2) + (3.4) kst2yc/yw · (Zb − Zw2) Now the equation of moments to the body around the on-axis: moments of the weight of the displaced mass cen- tres 0 = −8[(ztb − zqt)mb + (zts − zqt)ms]ag/ywZb+ moment of the inertia force from the roll motion +8s2mbr 2 bx +msr 2 sx]/yw · Zb moments of the forces in the suspension +ywSbw1 + ycSbw2 moment from the suspension of the seat −2ysks(Zs − ysZb/yw) moments from the lateral forces +(ztb − zqt)Sqb + (zts − zqt)Sqs moments of the roll motion of the wheels with inde- pendent suspension +2 · (mw1 + nzmw2)[2s 2(z2tw + r2w)/yw · Zb− 2ztwag/yw · Zb] moments from the lateral motion of the wheels +2ztw[(mw1 + nzmw2)Ac− s2(mw1Yky1 + nzmw2Yky2)] (3.5) Equation of the forces vertical to the seat 0 = s2msZs + ks(Zs − ys/yw · Zb) (3.6) Equation of vertical forces to the front wheels 0 = s2mw1Zw1 + kw1(Zw1 − Zk1)− Sbw1 (3.7) Equation of moments to the solid rear axle: From the inertia of the axle 0 = 2s2mw2(y 2 w/4 + tn · z2tw + tn · r2wx)· 2/yw · Zw2 from the unevennesses +ywkw2(Zw2 − Zk2) from the suspension −ycSbw2 from the lateral forces in the joint Sqb, Sqs +zqxw1/xw · (Sqb + Sqs) from the lateral motion of the axle −2tn · mw2(ztw − zq)(Ac − s2 · Yky2) from the weight of the displaced mass centre −2tn · mw2ztw · 2agZw2/yw from the weight of the body and seat to the displaced axle + 8(mb +ms)agzqZw2/yw · xw1/xw (3.8) Equation of lateral forces Lateral forces between the wheels and the road 0 = −2ky1Yky1 − 2ky2Yky2 inertia forces of the body and seat +Sqb + Sqs from solid axle roll +s2 · mw2 · ztw · 2/yw · (tn · Zw2 + nz · Zb) 18 Acta Polytechnica Vol. 50 No. 6/2010 from the lateral shift of the wheels +2(mw1 +mw2)Ac− s2 · (mw1Yky1 +mw2Yky2) (3.9) Finally, the equation of moments to the body around the z-axis passing through the center of grav- ity of the body: inertia moments of the body and seat 0= s2(4mbr 2 bz + 4msr 2 sz)· (Yky1 − Yky2 − zq · 2/yw · Zw2)/xw plus the inertia moments of the wheels +s2[2mw1(r 2 wz + x2w1 + y2w/4)+ 2mw2(r 2 wz + x2w2 + y2w/4)]· (Yky1 − Yky2 − zq · 2/yw · Zw2)/xw and the moment of the forces between the wheels and the road +2xw1ky1Yky1−2xw2ky2(Yky2+zq ·2/yw ·Zw2) (3.10) 4 The control Fig. 3 shows the control scheme of the suspension. In it, sensors are marked by their sensitivity constants. We will assume that active suspension will be used only on the rear axle. On the left side is the scheme of the front axle, with its spring kb1 and the damping bb1 and with static level control consisting of sensor fz, very low pass VLP, hydraulic valve V , central accumulator CA and drain DR. The dynamics of this control will be not dealt with here, i.e. Zr1 = 0. Some explanation at first. Elements marked SP are pressure sensors. A ring with minus in it produces the difference of two sig- nals. A ring withN in it means a negator — an element converting the phase by 180 deg. A ring with P in it means an amplifier, with am- plification proportional to the control input marked by an arrow. Hydraulic linkage is shown by double lines, and the signal paths are shown by single lines. The direc- tions of the signals are marked by triangular arrows. No feed-back through them is assumed. Signals com- ing to a common point add together. Crossing signal paths are not connected. On the left side, the equipment of front wheels 1, 2 is shown. On the right side, the equipment of rear wheels 3, 4 is shown. For reasons of stability, active suspension can be used only on one axle — let it be the rear axle. (Cars with front drive can use the front axle.) The lateral acceleration meter (with sensitivity constant tar) is mounted on the body at a height of zta. Fig. 3: Control scheme 19 Acta Polytechnica Vol. 50 No. 6/2010 The active suspension on the rear axle (right part of the figure) consists of the sensors of the displace- ment wheel-body marked by fz, producing the mean value of the front and rear displacements, and of the sensor of the vertical acceleration, marked by 1/fa, and of the actuator embodied by a pump driven by an electric motor. The actuator is marked A, and its control amplification is ζ2, the common amplication of signals is ζ. There is also a control velocity sen- sor, marked by its constant ζr, to develop a source of velocity by means of a negative feedback. Thus, the force of the actuator (proportional to the mass mb +ms) is Sa= ζ2sζr · Zr2 + ζζ2ζd[ζh · s2/(2πfa) · Zb+ 2πfzζb(Zb − Zw1/2− Zw2/2)− (Ac − s2Yky2 + 2s 2((zta − zqt)Zb+ zqtZw2)/yw − 2agZb/yw)tar]b0· (mb +ms)/m0 (4.1) where b0 is a formal constant valued 1N/ms and m0 is a constant with dimension kg. and where suit- able damping with damping coefficient βb is intro- duced by the element, marked by its transmission ζb = 1 + s · 2βb/(2πfz) Two filters are also included: a low pass with transmission ζd and cut-off frequency fd, which helps to prevent undamping of the wheel, and a high-pass with transmission ζh and cut-off frequency fh, which eliminates the false signal of the acceleration sensor of the tilted body. Namely ζd = 1/(1 + s/(2πfd)), ζh = s/(1 + s/(2πfh)) To the negative feedback: This control equa- tion makes the actuator approximately a source of hydraulic current, i.e. the control displacement is not dependent on the force Sbw2 . (This can also be achieved by very high inner stiffness of the actuator, but with great losses of power.) The force of the actuator is at the same time Sr = krZr2 + k0Zr2 − kb2(Zb2 − Zw2 − Zr2) (4.2) where Zr2 is the control displacement and k0 is the stiffness of the balancing spring in which the same static pressure is maintained as that in spring kb by means of sensors SP , differentiator D, very low pass V LP and valve V . (In the computed examples k0 = kbre.) The inner stiffness of the actuator kr consists of the damping of pump br and the inertia of the pump and the electric motor, so kr = sbr + s2mr Combining equations 4.1 and 4.2, a single control equation can be written, as follows 0= ζ2sζrZr2 + ζζ2ζd[ζh · s2/(2πfa) · Zb+ 2πfzζb(Zb − Zw1/2− Zw2/2)− (Ac − s2Yky + 2s2((zta − zqt)Zb+ zqtZw2)/yw − 2agZb/yw)tar]− [(sbr + s2 · mr)Zr + k0Zr− kb · (Zb − Zw − Zr)]m0/(mb +ms)/b0 (4.3) Cross control is provided on the front axle. This cross control uses the measurements of the pressure transducers, marked SP (in Fig. 3). The sum of the loads in the rear (after filtering through a very low pass V LP to get static values sst1 , sst2) controls the amplification of the proportional amplifier P of the load difference of the front wheel load. The outputs of the two proportional amplifiers are compared, and their difference controls the pumps of the front wheels with the appropriate phase. An equation to fulfil this aim can be written as follows 0= sZr1 + ζdpvcc[kbre2 (Zb+ yw/yc · Zr2 − Zw2)/(xw1/xw)/sst1− kbre1(Zb + Zr1 − Zw1)/(xw2/xw)/sst2 ] (4.4) where vcc is the control constant and where the trans- mission of the low path ζdp = 2πfdp/(2πfd + s) with the characteristic frequency fdp delays the an- swer to antiphased unevennesses. The influence of the inner stiffness of the actuator is not included here (a source of displacement is assumed). Equations 3.1 to 3.10, 4.3, 4.4 make a set of equa- tions for unknown quantities Sbw1 , Sbw2 , Sqb, Sqs, Zb, Zs, Zw1 , Zw2 , Yky1 , Yky2 , Zr1 , Zr2 . 5 Input and output The spectrum of a surface of medium roughness is used for the statistical form of the road unevennesses (ride velocity vx = 60 km/h): Zfk = √ (vx)/f · √(1/2− 1/2/(1 + (2π · f/4.5 · yw/v4x))/ (1 + (rpn · f/vx)2) · 0.64 mm. (5.1) where the track yw = 1.8 m. A correction to the tire diameter rpn is added. The effective value of this spectrum is 7.7 mm. The impulse input of the shape 1-cos (also “trans- lated sinusoidal”) was used in its spectral form τf =sin(2πfti)/(2πf)+ ti/2 · [sin(τ1)/τ1 + sin(τ2)/τ2], (5.2) 20 Acta Polytechnica Vol. 50 No. 6/2010 where τ1 = π − 2πfti, τ2 = π + 2πfti and where ti is the duration of the impulse at its half-height. The maximum value of the impulse unevenness zki is assumed to be proportional to its duration (pro- gressive impulses). These inputs are shown in the graphs of Zfbw or Zbw. The lateral input (lateral acceleration) is a sort of trapezoid-shaped impulse, but with rounded edges. It has been obtained by adding the two above- mentioned curves 1− cos with ti/2 (double cosinus), the second delayed by ti from the first. This input will be attached to the graphs of lat- eral seat acceleration. The graphs of the frequency characteristics Hf have the effective values attached. These effective values follow the formula√ 2 ∫ H2f df within limits 0 : fmax The pulse-effective values, given with time histories ht, are the effective values of the whole answer di- vided by the pulse duration, ie.√∫ h2tdt/ti) within limits t = 0− inf The vertical statistical input is due to the central path between the wheels. The time interval between the front and the rear of the vertical input depends on the axle base xw and riding speed vx tx = xw/vx. So the relation of the vertical input front and rear is Zk2 = Zk1(cos(stx) + i sin(stx) 6 Criteria Lateral acceleration in the rear seat Asy2 =Ac − s2 · Yky2+ s2[(ztb − zqt)Zb + zqZw2 ] · 2/yw− 2agZb/yw Vertical acceleration in the seat As = s2Zs Vertical acceleration of the body at radius yw/2 yw/2 Ab = s2Zb Body-wheel displacement front and rear Zbw1 =Zb − Zw1 , Zbw2 =Zb − Zw2 . Load transfer ratio front and rear Sw1/sst1 = kw1(Zw1 − Zk1)/sst1 , Sw2/sst2 = kw2(Zw2 − Zk2)/sst2 . where the static forces are sst1 = 2(mb +ms +mw)agxw2/xw, sst2 = 2(mb +ms +mw)agxw1/xw Seat-body displacements front and rear Zsb1 =Zs1 − Zb, Zsb2 =Zs2 − Zb Dynamic-to-static ratio of lateral forces between wheels and road, front and rear Sy1/sst1 = ky1Yky1/sst1 , Sy2/sst2 = ky2Yky2/sst2 . Spring displacements, front and rear Zp1 = Zb − Zw1 − Zr1 , Zp2 = (Zb − Zw2)yw/yc − Zr2 The forces in the suspension, front and rear, are Sbw1 = kb1(Zb1 − Zw1), Sbw2 = kb2(Zb2 − Zw2 − Zr2) 7 Universal parameters and used parameters The universal parameters (natural frequencies, dam- ping coefficients) defined for the heave-pitch model are also used in the roll model, but we must also consider the anti-roll bar. If the anti-roll bar stiffness-to-spring-stiffness kb ratio κstab is given, then kb= kbre + sbb = 4π2f2b (mb +ms)(1 + κstab)+ s · 4πβbfb(mb +ms) Parameters used with a high bus V variant X − A, i.e. cross control front and active control rear 21 Acta Polytechnica Vol. 50 No. 6/2010 fb = 1 Hz βb = 0.4 mb = 1 000 kg ztb = 0.9 m rbx = 0.6 m yw = 1.8 m fw = 10 Hz βw = 0.05 mw = 150 kg ztw = 0.5 m rwx = 0.3 m yc = 1.2 m fs = 3 Hz βs = 0.28 ms/mb = 0.5 zts = 2.7 m rsx = 0.4 m ys = 1.2 m ζ = 3 fz = 0.4 Hz βa = 1 fa = 0.3 Hz fd = 1 Hz fh = 0.1 Hz tar = 0.075 s ζ2 = 100 ζr = 0.9 br = bb/1 000mr = mb/1 000 b0 = 1 Ns/m m0 = 1 kg zta = 0.9 m vcc = 1 m/s fcc = 3 Hz κstab1 = 0 κstab2 = 0 zq = 0.3 m xw = 5.5 m xw1/xw2 = 2 rbz = 2 m rsz = 2 m rwz = 0.3 m Fig. 4: Time histories due to the impulse of lateral ac- celeration. A thick line means an active variant (act), a thin line means a passive variant (pas). The attached numbers are imp-eff. values. A semicolon is used to mark an index. In the graph of power, the attached value refers to energy consumption [Ws/1 000 kg] Stiffnesses kb, kw are assumed to be proportional to the static load front and rear: kb1 = kb · 2xw2/xw, kb2 = kb · 2xw1/xw, kw1 = kw · 2xw2/xw, kw2 = kw · 2xw1/xw For variant P − P , i.e. for a vehicle with passive suspension front and rear, the parameters differ as follows: κstab1 =1.2 κstab2 =1.4 zq =0.9 m ζ2=100 The parameters of actuators mr and br are only estimated. 8 Example of modelling results In the graphs in Figs. 4, 5, the active variant act (passive suspension front, active rear) is indicated by more prominent lines (lines with dots) than those of the passive variant (pas–pas). In some suitable graphs, the appropriate input is shown by interrupted lines. Effective values (with frequency characteristics) or impulse-effective values (effective values of the an- swers divided by the impulse duration, with time his- tories) are attached to the labels. Fig. 5: Frequency characteristics due to unevennesses. A thick line means an active variant (act), a thin line means a passive variant (pas). The attached numbers are imp- eff. values. A semicolon is used to mark an index 22 Acta Polytechnica Vol. 50 No. 6/2010 active passive Steady-state values due to ac = 1 m/s2 Lateral acceleration in the seat asy [mm/s2] 0.77 1.35 Body-wheel travel zbw2 [mm] 25 25 Vertical dynamic force to static force ratio, front sw1/sst1 12 20 Vertical dynamic force to static force ratio, rear sw2/sst2 14 20 Steady-state values due to in-phase unevenness zk1 = zk2 = 0.01 m Lateral acceleration in the seat Asy [m/s2] 0.11 0.15 Vertical dynamic force to static force ratio, front Sw1/sst1 0.02 0.022 Vertical dynamic force to static force ratio, rear Sw2/sst2 0.017 0.019 Steady-state values due to cross unevenness zk1 = 0.01 m, zk2 = −0.01 m Lateral acceleration in the seat Asy [m/s2] 0.028 0.023 Vertical dynamic force to static force ratio, front Sw1/sst1 0.000 42 −0.07 Vertical dynamic force to static force ratio, rear Sw2/sst2 0.009 0.03 9 Stability Stability against self-exciting oscillation is an impor- tant criterion of the system. It is most important with the roll model, where there is also the possi- bility of roll-over. Asymptotic stability is used, i.e. the real part of the eigenvalues of the equation ma- trix must be negative. Stability was checked for all parameters of a simplified model. The parameters of nominal value hnom were varied between 1/100 to 100 multiples of its given value hjm. When no insta- bility was found, output 100 was given. If instability was found at parameter value hcrit, the stability rate hcrit/hnom was put out. A minimal stability rate of about 2.2 of the passive model increases to 5 with active suspension. 10 Conclusions based on the example With vertical statistical input, the lateral acceler- ation in the seat is slightly increased and the lat- eral force is substantially decreased. The lateral load transfer is slightly deteriorated, but the effective val- ues incorporating the values from the heave-pitch model are little influenced. It has been shown that it is possible by means of cross control to distribute the load changes due to cornering proportional to the static load of the axle. The lateral input is the main advantage of active suspension, thanks to the active roll: From the viewpoint of ride comfort, active roll substantially suppresses the lateral force during cor- nering. From the viewpoint of ride safety, active roll also substantially suppresses the load transfer to the outside wheels in a curve. Thus the vehicle can take curves at a higher speed or with more safety. With active roll suspension there is no need to use a solid axle for anti-roll purposes. This means more comfort for the passenger and there is reduced dam- age to the road, due to a big reduction in the lateral forces between wheel and road. A solid axle was used in our examples, but with reduced height of the roll axis. With active suspension there is also no need for anti-roll bars. Active roll technology also enables the design of a slender passenger car for two passengers sitting in tandem (we can call this a tandemo) with, e.g., an 0.8 m track. References [1] Sampson, D. J. M., Cebon, D.: Active Roll Control of Single Unit Heavy Road Vehicles, http://www.cvdc.org/recent papers/ SampsonCebon VSD02.pdf [2] Vacuĺın, O., Valášek, M., Svoboda, J.: Influence of Vehicle Tilting on Its Performance, Conf. In- teraction and Feedbacks, Praha, 2005. [3] Masashi, Yamashita: Application of H-Y Con- trol to Active Suspension Systems, Automatica, Vol. 30, 1994. [4] Čech, I.: A Pitch-Plane Model of a Vehicle with Controlled Suspension. Vehicle System Dynam- ics, Vol. 23, 1998, pp. 133–148. Ing. Ilja Čech E-mail: cech10600jes11@seznam.cz Jesenická 11, 106 00 Prague 10, Czech Republic 23 Acta Polytechnica Vol. 50 No. 6/2010 List of symbols Complex amplitudes, spectral values and instantaneous values Ac Afc ac lateral (centripetal) acceleration [m/s2] As Afs as vertical acceleration of the seat [m/s2] Asy Afty asy lateral acceleration in the seat above the center of gravity [m/s2] Ax Afx ax longitudinal acceleration input [m/s2] Sa Sfa sa force developed by the actuator [N] Sbw Sfbw sbw force between body and wheel [N] Sq Sfq sq lateral inertia force due to yaw motion [N] Sp Sfp sp force on the spring [N] Sqb Sfqb sqb lateral force in the joint due to the body [N] Sqs Sfqs sqs lateral force in the joint due to the seat [N] Ssb Sfsb ssb vertical force between seat and body [N] Sw Sfw sw vertical force between wheel and road surface [N] Sy Sfy sy lateral force between wheel and road surface [N] Vr Vfr vr velocity of movement of the actuator [m/s] Yar Yfar yar lateral displacement of the body at lateral acceleration meter height [m] Yky Yfky yky lateral displacement in lateral stiffnesses ky1 , ky2 Zb Zfb zb vertical displacement (for the roll model at radius yw/2) [m] Yz Yfz yz lateral displacement of the body due to yaw motion [m] Zbw Zfbw zbw body-wheel displacement [m] Zk Zfk zk vertical input displacement [m] Zr Zfr zr control displacement [m] Zs Zfs zs vertical seat displacement [m] Zw Zfw zw vertical wheel displacement [m] Quantities without physical dimension βa relative damping rate of the control βb relative damping rate of the body βs relative damping rate of the seat ζ common amplification of signals ζ2 amplification of the actuator ζb corrective damping element ζd transmission of low pass filter ζh transmission of high pass filter ζr sensitivity constant of feedback sensor κstab ani-roll bar stiffness κyw lateral stiffness to radial stiffness ratio nz, tn nz = 0, tn = 1 for solid axle, nz = 1, tn = 0 for no solid axle Other quantities ag gravity acceleration [m/s2] bb damping of passive suspension [Ns/m] bpk slip damping of the tyre [Ns/m] b0 formal constant of the actuator b0 = 1 [Ns/m] br damping of the actuator [Ns/m] bs seat damping [Ns/m] bw wheel damping [Ns/m] f frequency [Hz] fa 1/(2πfa) sensitivity constant of transducer of vertical acceleration [Hz] fb natural frequency of passive suspension [Hz] fd, fh characteristic frequency of the filter [Hz] 24 Acta Polytechnica Vol. 50 No. 6/2010 fr sensitivity constant of the control velocity sensor [Hz] fz 2πfz constant of the displacement transducer [Hz] fε 1/(2πfε) constant of the roll acceleration transducer [Hz] hj nominal value of the parameter hkrit critical value of the parameter k0 spring rate of the balancing spring [N/m] kb stiffness of the body-wheel suspension [N/m] kbre spring rate of the body-wheel spring [N/m] kr inner stiffness of the actuator [N/m] ks stiffness of the seat-body spring [N/m] kst stiffness of the roll-bar [N/m] kw radial stiffness of the tyre [N/m] kwyre real lateral stiffness of the tyre [N/m] ky lateral stiffness [N/m] kypk sliding stiffness [N/m] m0 formal constant m0 = 1 kg mb quarter body mass [kg] ms quarter seat mass [kg] mc total mass of vehicle mb +ms +mw [kg] mr actuator mass transferred to the pump radius [kg] mw wheel mass [kg] p power [W] rpn tyre radius [m] rbx gyration radius of the body to the longitudinal axis [m] rbz gyration radius of the body to the vertical axis [m] rsx gyration radius of the seat to the longitudinal axis [m] rsz gyration radius of the seat to the vertical axis [m] rwx gyration radius of the wheel to the longitudinal axis [m] rwz gyration radius of the wheel to the vertical axis [m] s operator i2πf t time [s] tar sensitivity constant of the lateral acceleration sensor [Hz] tgal constant of inverse proportionality between slip rotation and load ti duration of input impulse at half-height [s] vx travelling speed [m/s] yc spring distance [m] ys seat distance [m] yw wheel track [m] za height of the lateral acceleration meter [m] zki maximum value of the vertical input impulse [m] zq height of the joint [m] zqt height of the axis on under the center of gravity [m] zta height of the lateral acceleration transducer [m] xw axle base [m] xw1 distance between the mass centre of the body and the front axle [m] xw2 distance between the mass centre of the body and the rear axle [m] ztb height of the body mass centre [m] zts height of the seat mass centre [m] ztw height of the wheel mass centre [m] Indices 1 front axle 2 rear axle b body f spectral values re real part s seat x longitudinal direction, axis y lateral direction, axis w wheel 25