This invention relates to a system and method of protecting against rollover in a motor vehicle.
Dynamic control systems have been recently introduced in automotive vehicles for measuring the body characteristics of the vehicle and controlling the dynamics of the vehicle based on the measured body characteristics. For example, certain systems measure vehicle characteristics to prevent vehicle rollover and for tilt control (or body roll). Tilt control maintains the vehicle body on a plane or nearly on a plane parallel to the road surface, and rollover control maintains the vehicle wheels on the road surface. Certain systems use a combination of yaw control and tilt control to maintain the vehicle body horizontal while turning. Commercial examples of these systems are known as active rollover prevention (ARP) and rollover stability control (RSC) systems.
The spectrum of conditions that may occur during the operation of the vehicle is too large to be practically evaluated during the development and production of the vehicle. As a result, the tuning of the rollover stability control system for the vehicle is typically performed with an extreme roof load to provide sufficient confidence that the system will perform suitably over road conditions that the vehicle will experience when being driven.
However, when the rollover stability control system is tuned in the roof-loaded condition, the gains are higher than those that would result from tuning in the normal-loaded condition. Thus, the system becomes very sensitive to small disturbances.
Moreover, conventional systems consider the longitudinal vehicle dynamics to estimate the mass of the system. Hence, these systems do not provide an indication about the way the mass is distributed with respect to the roll axis (i.e., the roll moment of inertia).
In satisfying the above need, as well as overcoming the enumerated drawbacks and other limitations of the related art, the present invention provides a system and method that estimates a parameter related to the mass of the loaded vehicle as well as the mass distribution.
In a general aspect of the invention, the system detects a lateral acceleration and roll rate of the vehicle and estimates a mass distribution parameter. The system then generates a tuned mass distribution parameter that is based on the the lateral acceleration, the roll rate, and the mass distribution parameter and introduces the tuned mass distribution parameter to a rollover stability control system.
Further features and advantages will become apparent from the following description, and from the claims.
Referring now to
With reference to
ΣM=O Eqn. 1
which yields
I{umlaut over (θ)}=mhu−K θ−c{dot over (θ)}, Eqn. 2
where
I is the moment of inertia in the roll direction,
{umlaut over (θ)} is the roll acceleration,
{dot over (θ)} is the roll rate,
θ is the roll angle,
m is the total mass,
K is the roll stiffness,
c is the roll damping coefficient,
u is the lateral acceleration, and
h is the height of the center of gravity from the roll axis.
Rearranging Eqn. 2 provides
Thus, in state space, the continuous time system roll model is
x(t)=Ax(t)+Bu(t), x(0)=0 if the initial time is set to zero,
The C matrix is chosen depending on the type of sensor employed. In this case, a roll rate sensor is being employed, hence C=[0 1].
Note that the above discussion is directed to obtaining a solution for the state vector x(t) in continuous time. Therefore, the system described in Eqn. 4 is typically discretized according to the expression
{dot over (x)}(k)=Adx(k−1)+Bdu(k) Eqn. 5
y(k)=Cx(k)
where k identifies the kth time step and
Ad=In+AT,
Bd=BT,
and where
In is the nth order identity matrix, which in this case is a second order identity matrix, and
T is the time step.
Converting the discretized state space equation (Eqn. 5) to transfer function space identified here as z yields:
Hence, U(z) is the lateral acceleration in z space and Y(z) is the corresponding roll rate.
Expanding Eqn. 7 yields:
Y(z)+d1Y(z)z−1+d2Y(z)z−2=n1U(z)z−1−n2U(z)z−2, or
Y(z)=n1U(z)z−1−n2U(z)z−2−d1Y(z)z−1−d2Y(z)z−2 Eqn. 8
which can be generalized as
Y(·)=UT(·)*{circumflex over (θ)} Eqn. 9
where {circumflex over (θ)}=[n1, n2, d1, d2] is the parameter vector and Y(·), UT(·) are known (i.e., measured).
Since d1, and d2 are not functions of m and h, d1 and d2 can be calculated in advance so that only n1 and n2 need to be estimated. Observing that n1=n2, the inverse z transform of the transfer function is
Y(k)+d1Y(k−1)+d2Y(k−2)=n1[U(k−1)−U(k−2)], Eqn. 10
In this way, the variables of Eqn. 9 are scalar. The estimated parameter n1 is a function of the vehicle mass and moment of inertia of the body about the roll axis.
Turning again to
The estimate for the mass distribution parameter n1 can be performed by a recursive least squares (RLS) method or any other suitable method. An example of an RLS algorithm used in conjunction with the system shown in
phi=P*u;
gamma(k)=phi/(u*phi+lambda);
Y_hat(k)=u*n_hat(:,k−2);
n_hat(:,k−1)=n_hat(:,k−2)+gamma(k)*(Y(k)−Y_hat(k));
if (abs(U(k))>3
n_hat(:,k−1)=n_hat(:,k−2);
Other embodiments are within the scope of the following claims.