Limitations To Suspension Performance in A Two-Degree-Of-Freedom Car Active Suspension
Limitations To Suspension Performance in A Two-Degree-Of-Freedom Car Active Suspension
Then Fig.2 stands for the following algebraic equations: It is obvious from (19) and (20) that the sprung mass position
transfer function HSP has an invariant “tire-hop” frequency at
𝑍(𝑠) = 𝐺11 (𝑠)𝑊(𝑠) + 𝐺12 (𝑠)𝑈(𝑠) (5)
𝜔1 = 𝑘𝑡 /𝑚𝑢 , where:
𝑌(𝑠) = 𝐺21 (𝑠)𝑊(𝑠) + 𝐺22 (𝑠)𝑈(𝑠) (6)
𝐻𝑆𝑃 (𝑠)/𝑠=𝑗 𝜔 1 = −𝑚𝑢 /𝑚𝑠 (22)
𝑈(𝑠) = 𝐾(𝑠)𝑌(𝑠) (7)
Similarly, from (16) and (18) the suspension deflection transfer
Manipulating the equations listed above, the following transfer function HSD has an invariant “rattle-space” frequency at
function Tzw(s) from w to z as a linear-fractional transformation of 𝜔2 = 𝑘𝑡 /(𝑚𝑢 + 𝑚𝑠 ) and:
K(s) can be derived:
−1 𝐺 𝐻𝑆𝐷 (𝑠)/𝑠=𝑗 𝜔 2 = −(1 + 𝑚𝑢 /𝑚𝑠 ) (23)
𝑇𝑧𝑤 = 𝐺11 + 𝐺12 𝐾 𝐼 − 𝐺22 𝐾 21 =
= 𝐺11 + 𝐺12 𝐼 − 𝐾𝐺22 −1 𝐾𝐺21 (8) Finally, from (20) and (21), the tire deflection transfer function
HTD does not have any invariant frequency point except 𝜔3 = 0,
It is shown in [1] that the set of all proper real-rational matrices
where:
K(s) stabilizing G(s) is parametrized by a free parameter 𝑄(𝑠) ∈
𝑅𝐻∞ as follows: 𝐻𝑇𝐷 (𝑠)/𝑠=𝑗 𝜔 3 = 0 (24)
−1
𝐾 = 𝑌 − 𝑀𝑄 𝑋 − 𝑁𝑄 =
−1
= 𝑋 − 𝑄𝑁 𝑌 − 𝑄𝑀 (9) 4. Transfer Functions Limitations
where: In the next, consider the standard block diagram shown in
M(𝑠), 𝑁(𝑠), 𝑋(𝑠), 𝑌(𝑠) ∈ 𝑅𝐻∞ , and Fig.2.
𝑀(𝑠), 𝑁(𝑠), 𝑋(𝑠), 𝑌(𝑠) ∈ 𝑅𝐻∞ can be found by comprime As an example, consider:
factorization approach of G22(s): w=zr , z=zs 𝑦 = 𝑧𝑢 , 𝑧𝑠 − 𝑧𝑢 , 𝑧𝑢 − 𝑧𝑟 𝑇 .
Then:
𝐺22 (𝑠) = 𝑁(𝑠)𝑀−1 (𝑠) = 𝑀 −1 (𝑠)𝑁(𝑠) (10)
𝐺11 (𝑠) 𝐺12 (𝑠)
𝑋(𝑠) −𝑌(𝑠) 𝑀(𝑠) 𝑌(𝑠) 𝐺(𝑠) = =
=𝐼 (11) 𝐺21 (𝑠) 𝐺22 (𝑠)
−𝑁(𝑠) 𝑀(𝑠) 𝑁(𝑠) 𝑋(𝑠) 0
1
𝑚 𝑠𝑠2
Substituting the equation (9) into (8) we obtain the transfer 𝑠𝑘 𝑡 −𝑠
matrix Tzw(s) from w to z in terms of the free parameter 𝑄(𝑠) ∈ 𝑚 𝑢 𝑠 2 +𝑘 𝑡 𝑚 𝑢 𝑠 2 +𝑘 𝑡
= −𝑘 𝑡 (𝑚 𝑢 +𝑚 𝑠 )𝑠 2 +𝑘 𝑡 (25)
𝑅𝐻∞ :
𝑚 𝑢 𝑠 2 +𝑘 𝑡 𝑚 𝑠 𝑠 2 (𝑚 𝑢 𝑠 2 +𝑘 𝑡 )
𝑇𝑧𝑤 (𝑠) = −𝑚 𝑢 𝑠 2 −1
𝑚 𝑢 𝑠 2 +𝑘 𝑡 𝑚 𝑢 𝑠 2 +𝑘 𝑡
= 𝐺11 (𝑠) + 𝐺12 (𝑠)𝑀(𝑠) 𝑌(𝑠) − 𝑄(𝑠)𝑀(𝑠) 𝐺21 (𝑠) =
= 𝐺11 (𝑠) + 𝐺12 (𝑠) 𝑌(𝑠) − 𝑀(𝑠)𝑄(𝑠) 𝑀(𝑠)𝐺21 (𝑠) The limitations of all achievable closed-loop transfer functions
(12) 𝑇𝑧𝑤 (𝑠) = 𝐻𝑆𝑃 (𝑠) are derived from the right and left comprime
factorization of G22(s), i.e.:
As the parameter Q(s) varies over the set of all stable proper
functions, the equation (12) parametrizes all achievable transfer −𝑚𝑠 𝑠 3
functions Tzw(s). 𝑝4 (𝑠)
−1
If it is assumed that the tire does not leave the ground, for the (𝑚𝑢 + 𝑚𝑠 )𝑠 2 + 𝑘𝑡 𝑚𝑠 𝑠 2 (𝑚𝑢 𝑠 2 + 𝑘𝑡 )
one-quarter car model (Fig.1) the linear differential equations of 𝐺22 = =
𝑝4 (𝑠) 𝑝4 (𝑠)
motion are:
−𝑚𝑠 𝑠 2 𝑀 −1 (𝑠)
𝑚𝑠 𝑧𝑠 = 𝑢 − 𝑓 (13) 𝑝4 (𝑠)
𝑁(𝑠)
𝑚𝑢 𝑧𝑢 = −𝑢 + 𝑘𝑡 (𝑧𝑟 − 𝑧𝑢 ) (14) −1
𝑚 𝑠𝑠 𝑚 𝑠𝑠2
0 1
where zu and zs are measured from the static equilibrium position. 𝑝 2 (𝑠) 𝑝 2 (𝑠)
𝑝 2 (𝑠)
𝑚 𝑢 𝑚 𝑠𝑠2 −𝑚 𝑠 𝑘 𝑡
First, let the load disturbance is absent (f=0). Adding equations = 0 𝑚 𝑢 +𝑚 𝑠 (26)
𝑝 2 (𝑠) 𝑝 2 (𝑠) 𝑝 2 (𝑠)
(13) and (14) we obtain the invariant equation of: −1 𝑠
0 0
𝑝 2 (𝑠) 𝑝 2 (𝑠)
𝑚𝑠 𝑧𝑠 + 𝑚𝑢 𝑧𝑢 = 𝑘𝑡 (𝑧𝑟 − 𝑧𝑢 ) (15) 𝑀 −1
𝑁 (𝑠)
that is independent on the suspension force u. The following where p2(s) and p4(s) are Hurwitz polynomials of degree 2 and
transfer functions will be investigated: 4, respectively. Then:
𝐻𝑆𝑃 (𝑠) = 𝑍𝑠 (𝑠)/𝑍𝑟 (𝑠) (16) 𝐻𝑆𝑃 (𝑠) =
𝐻𝑆𝐷 (𝑠) = 𝑍𝑠 (𝑠) − 𝑍𝑢 (𝑠) /𝑍𝑟 (𝑠) (17) = 𝐺11 (𝑠) + 𝐺12 (𝑠)𝑀(𝑠) 𝑌(𝑠) − 𝑄(𝑠)𝑀(𝑠) 𝐺21 (𝑠) =
−𝑠𝑘𝑡
𝐻𝑇𝐷 (𝑠) = 𝑍𝑢 (𝑠) − 𝑍𝑟 (𝑠) /𝑍𝑟 (𝑠) (18) 1 𝑠(𝑚 𝑢 𝑠 2 +𝑘 𝑡 ) ∗
=− 𝑌(𝑠) 𝑘𝑡 + 𝑄 (𝑠) (27)
𝑝 4 (𝑠) 𝑝 4 (𝑠)
𝑚𝑢 𝑠 2
3. Invariant Properties
𝑄3 𝑠
Manipulating the equation (15) we can derive the following where𝑄 ∗ 𝑠 = , 𝑄(𝑠) = 𝑄1 𝑠 𝑄2 𝑠 𝑄3 (𝑠) .
𝑝2 𝑠
invariant identities:
It follows from (27) that thanks to the term 𝑌(𝑠) − 𝑄(𝑠)𝑀(𝑠)
(𝑚𝑢 + 𝑚𝑠 )𝑠 2 + 𝑘𝑡 𝐻𝑆𝑃 (𝑠) − 𝑚𝑢 𝑠 2 + 𝑘𝑡 𝐻𝑆𝐷 (𝑠) = 𝑘𝑡 (19) bounded for large s:
𝑚𝑠 𝑠 2 𝐻𝑆𝑃 (𝑠) + 𝑚𝑢 𝑠 2 + 𝑘𝑡 𝐻𝑇𝐷 (𝑠) = −𝑚 𝑢 𝑠 2 (20) 𝑙𝑖𝑚𝑠 2 𝐻𝑆𝑃 (𝑠) < ∞ (28)
𝑠→∞
(𝑚𝑢 + 𝑚𝑠 )𝑠 2
+ 𝑘𝑡 𝐻𝑇𝐷 (𝑠) + 𝑚𝑠 𝑠 2 𝐻𝑆𝐷 (𝑠) = (21)
i.e., the resulting rate of decay is of second degree:
= −(𝑚𝑢 + 𝑚𝑠 )𝑠 2
−2
𝐻𝑆𝑃 (𝑠)/𝑠→∞ = 𝑂(𝑠 ) (29) and suspension deflection velocity are measured, is as follows
[2],[3]:
𝑠(𝑚 𝑢 𝑠 2 +𝑘 𝑡)
It is obvious from (27) that the member has two 1
𝑝 4 (𝑠) 𝐻𝑧𝑤 (𝑠)/𝑠→∞ = 𝑂(𝑠 −3 ) (infinite frequency constraint)
imaginary axis zeros - at s=0 and 𝑠 = 𝑗𝜔1 = 𝑗 𝑘𝑡 /𝑚𝑢 - which can
not be canceled by the denominator of 𝑄 ∗ (𝑠) ∈ 𝑅𝐻∞ . With respect
1
H zw ( s ) / s 0 1 O( s 2 )
(zero frequency constraint)
to (27) and the Bezout identity (11) it follows, that:
1 𝑚𝑢 𝑘𝑡
−𝑠𝑘𝑡 𝐻𝑧𝑤 (𝑗𝜔1 ) = − for 𝜔1 =
1 𝑚𝑠 𝑚𝑢
𝐻𝑆𝑃 (𝑠)/𝑠=0 = − 𝑌(𝑠) 𝑘𝑡 /𝑠=0 =
𝑝4 (𝑠) and analogically:
𝑚𝑢 𝑠 2
𝑘𝑡 𝑘𝑡
=− 𝑌 (𝑠)/𝑠=0 = 1 (30) 2 (𝑠)/
𝐻𝑧𝑤 𝑠→∞ = − 𝑠 −2 + 𝑂(𝑠 −3 ) (infinity freq. constraint)
𝑝 4 (𝑠) 2 𝑚𝑢
adjusting the ratio of the unsprung and sprung masses. From the 6. Results
result:
Using a simple linear two-degree of freedom car suspension model
3 𝑘𝑡
𝐻𝑧𝑤 (𝑠)/𝑠→∞ = −1 + 𝑠 −2 + 𝑂(𝑠 −3 ) (38) in Fig.1, it was shown that there are still some limitations to
𝑚𝑢
suspension performance even in the fully active state. It has been
3
it is evident, that 𝐻𝑧𝑤 (𝑗𝜔) must tend to one from above as 𝜔 shown in the paper that there are some frequency points and
3
tends to ∞ and it turns out that 𝐻𝑧𝑤 (𝑗𝜔) cannot be made less than frequency ranges where the transfer functions have modulus strictly
or equal to one at all frequencies. Since the right-hand side of (37) greater than one i.e. where road and load disturbance amplification
3 (𝑗𝜔) to be less than occurs. On the base of the Bode integral theorems, it has been
is non-negative then it is not possible for 𝐻𝑧𝑤
or equal to one at all frequencies since that would make the left- shown that the transfer functions must be greater in modulus to at
hand side of the equation (37) negative. It has been shown in [3], least the same extent that it is less than one when measured in terms
that no matter what signals are used for feedback, the tire deflection of the area on a Bode magnitude plot. In such a case there is a
transfer function must amplify road disturbances at some possibility to shift frequency ranges where disturbance
frequencies. This fact is valid even for full state feedback used in amplification occurs to a ”more advantageous place” or to make
the control loop. magnitudes lower spreading the frequency range.
A similar theorem is valid for transfer functions where Acknowledgment
𝐻(𝑠)/𝑠→0 = 1 + 𝑂(𝑠 2 )
This research has been supported by MSMT project INTER-
Theorem 2. Let 𝐻(𝑠) belongs to 𝑅𝐻∞ and satisfies𝐻(𝑠)/𝑠→0 = VECTOR 17019.
1 + 𝑂(𝑠 2 ).Let 𝑧𝑖 , 𝑖 = 1, . . , 𝑛 𝑛 ∈ 𝑁 are zeros of 𝐻(𝑠) with
𝑅𝑒( 𝑧𝑖 ) > 0. Then:
∞
𝑙𝑜𝑔 𝐻(𝑗𝜔)
𝑑𝜔
=𝜋 𝑛 𝑅𝑒 (𝑧 𝑖 )
(39)
References
0 𝜔2 𝑖=1 𝑧 2
𝑖
[1] Freudenberg, J. S., Looze, D. P., Right-half Plane Poles and
Similarly to the consequences of Theorem 3, the result Zeros and Design Tradeoffs in Feedback Systems, IEEE
1 (𝑠)/ 2
𝐻𝑧𝑤 𝑠→0 = 1 + 𝑂(𝑠 ) from Section 1 is the case when Transactions on Automatic Control, 30/6, 1985
1
𝐻𝑧𝑤 (𝑗𝜔) cannot be less than or equal to one at all frequencies
since that would make the left-hand side of (39) negative. This fact [2] Hyniova, K., Frequency Response Limitations for Active
is again valid even when full state feedback is introduced in the Suspension of Vehicles, Micro-CAD 99 – International Computer
control loop [3]. Science Conference, Miskolc, pp. 209-217, 2006,
In such cases that were mentioned above designers have the [3] Hyniova, K., Achievable Dynamic Responses for Active
only possibility ”to shift” the frequencies where the amplifications Suspension System of Vehicles, Habilitation thesis, Prague, 2000
occur to ”more advantageous places” or ”to spread” the ranges [4] Zhou, K., Doyle, J.C., Essentials of Robust Control, Prentice
where amplifications occur making the amplification lower the Hall, London, 1998
positive area of the Bode magnitude plot, i.e. the area where
𝐻(𝑗𝜔) is greater than one (is grater to at least the same extent than
the negative area where 𝐻(𝑗𝜔) is less than one) by choosing a
proper feedback controller. Analogically, similar results can be
derived for arbitrarily chosen measurements and load disturbances.