a is a graphical illustration of aerodynamic side force versus time of a vehicle during straight driving at 120 kph with side winds of 40 kph occurring at 2 seconds;
b is a graphical illustration of aerodynamic moment versus time of the vehicle of
c is a graphical illustration of lateral acceleration versus time of the vehicle of
d is a graphical illustration of yaw rate versus time of the vehicle of
a is a graphical illustration of aerodynamic side force versus time of a vehicle during straight driving at 120 kph with side winds of 80 kph occurring at 2 seconds;
b is a graphical illustration of aerodynamic moment versus time of the vehicle of
c is a graphical illustration of lateral acceleration versus time of the vehicle of
d is a graphical illustration of yaw rate versus time of the vehicle of
a is a graphical illustration of actual and estimated lateral velocities versus time of a simulated vehicle traveling at 120 kph in response to a 80 kph crosswind occuring after 2 seconds with the estimate obtained by applying the disclosed system and method for estimating lateral velocity;
b is a graphical illustration of actual and estimated aerodynamic side force versus time of the vehicle of
a is a graphical illustration of lateral acceleration versus time of a simulated vehicle traveling at 120 kph in response to a 80 kph crosswind occuring after 2 seconds for a conventional vehicle and a vehicle applying the disclosed system and method for estimating and compensating the lateral disturbances using the front steering input;
b is a graphical illustration of yaw rate versus time of the vehicle of
c is a graphical illustration of lateral deviation versus time of the vehicle of
d is a graphical illustration of front steering angle versus time of the vehicle of
Referring to
Mzs=Fyse (Eq. 1)
wherein e is the distance between the center of mass 16 of the vehicle 12 and the center of aerodynamic pressure 14. In one aspect, for a given vehicle, the distance e may be known, at least approximately, thereby leaving only Fys as the unknown variable.
Under steady-sate conditions, the side force Fys imposed upon the vehicle 12 due to the side winds 10 may be modeled as follows:
and the yaw moment Mzs may be modeled as follows:
wherein ρ is the air density, A is the frontal area of the vehicle 12, νw is the wind velocity relative to the vehicle, L is vehicle wheelbase, CS and CM, are the side force and the yawing moment coefficients, respectively, both of which may depend upon the relative wind angle βw. Applying Eq. 1 to Eqs. 2 and 3 it can be seen that the side force and the yaw moment coefficients may be proportional to each other as follows:
C
M
=C
Se/L (Eq. 4)
The aerodynamic side force Fys and yaw moment Mzs may directly affect the equations for lateral and yaw motions. The magnitudes of the aerodynamic force and moment may depend upon the square of the wind velocity relative to the vehicle 12. Under most driving conditions they may be relatively small and may become significant only when both the vehicle forward speed and the side wind velocity are relatively large.
Those skilled in the art will appreciate that the magnitudes of the lateral forces acting on the vehicle seldom exceed 10 percent of the vehicle weight and, therefore, the resulting lateral acceleration response seldom exceeds 0.1 g and the resulting yaw rate is usually below 4 deg/s. For example, referring to
Referring to
Fyg=mg sin γ (Eq. 5)
wherein m is total mass of the vehicle 12′ and g is the acceleration of gravity.
Those skilled in the art will appreciate that there is no yaw moment with respect to the center of gravity. However, the lateral force Fyg may directly affect the equation of lateral motion of the vehicle, but only indirectly the yaw motion equation. In addition, the gravity component may directly affect the measured lateral acceleration, since the gravity force may have the same effect on the accelerometer as the inertial force. Thus the measured lateral acceleration aym may be modeled as follows:
a
ym
=a
y
−g sin γ={dot over (ν)}y+νxΩ−g sin γ (Eq. 6)
wherein ay is the actual lateral acceleration, νx and νy are the longitudinal and lateral velocities, respectively, and Ω is the vehicle yaw rate.
In one aspect, the linear bicycle model may be applied, wherein the estimates are more accurate when the vehicle remains within the linear range of handling. Therefore, in addition to the tire forces, the external forces and moment acting on the vehicle 12 may include the lateral force due to side wind Fys, the aerodynamic yawing moment Mzs and the lateral force due to the road bank angle Fyg. In the linear range of handling, the vehicle dynamics in the yaw plane may be modeled as follows:
wherein cf and cr denote cornering stiffness values of both front and both rear tires, respectively, a and b are the distances of the front and rear axles to the center of mass of vehicle, Izz is the vehicle yaw moment of inertia and δf is the front wheel steering angle.
In one aspect, Eq. 7 may be obtained from the balance of forces in the lateral direction and Eq. 8 may be obtained from the balance of moments about the vertical axis.
Substituting the derivative of lateral velocity {dot over (ν)}y from Eq. 6 into Eq. 7 and denoting the aerodynamic side force disturbance per unit mass w follows:
yields the following system of equations:
Since vehicle parameters may be known (at least approximately), the lateral acceleration aym and the yaw rate Ω may be measured and the vehicle speed νx may be estimated, Eqs. 10 and 11 are a system of two equations with only two unknown values, namely lateral velocity νy and the disturbance w due to the aerodynamic force. Solving Eqs. 10 and 11 for the unknown variables yields:
From Eqs. 12 and 13 an estimate of the lateral force disturbance due to side wind Fys can be obtained using only known signals. The time derivative of yaw rate (occurring in Eq. 12) may in practice be approximated by passing the measured yaw rate Ω through a high pass filter in order to reduce the effect of measurement noise. After the lateral force disturbance Fys is determined, the estimate of the yawing moment due to side wind can be determined, if desired, from Eq. 1.
The estimate of bank angle disturbance, g sin γ, can then be obtained from Eq. 6 to yield:
g sin γ={dot over (ν)}y+νxΩ−aym (Eq. 14)
and the estimate may then be passed through a low pass filter in order to reduce the effect of noise. The derivative of lateral velocity in Eq. 14 may be obtained by differentiating the lateral velocity obtained from Eq. 12. The total lateral force disturbance per unit mass of the vehicle is the sum of the disturbances resulting from the side wind w and due to the bank angle g sin γ.
Those skilled in the art will appreciate that some of the equations described above may be modified to accommodate vehicles equipped with an active rear steer system. For example, when the rear wheels may be steered with an angle δr, Eqs. 10 and 11 may take the following form:
Consequently the term cf(a−e)δf in Eq. 12 may be replaced with cf(a−e)δf−cr(b+e)δr and the term
in Eq. 13 may be replaced with
to yield the following equations:
When the vehicle is subjected to a sufficiently large total disturbance and the driver does not provide any significant steering correction, an automatic steering or brake correction may be determined. Since the yaw rate may have a large influence on the vehicle deviation from the desired path in the linear handling range of the vehicle, the correction may be selected to eliminate the steady state value of yaw rate caused by the disturbance. When the vehicle is equipped with either active front steer (AFS) or active rear steer (ARS) system, a front or rear steering correction, respectively, may be applied. If the vehicle posses a brake-based electronic stability control (ESC) system, an asymmetric brake intervention may be initiated.
In one aspect, the objective may be to determine the front steering correction. Therefore, Eq. 9 may be substituted into Eqs. 7 and 8 to yield:
At steady-state, the time derivatives of lateral velocity and yaw rate may generally equal zero (i.e., {dot over (ν)}y=0 and {dot over (Ω)}=0) such that Eqs. 19 and 20 become:
wherein the symbols ss refer to the steady-state values.
Solving Eqs. 21 and 22 for the unknown values νyss and Ωss yields the following steady state value of yaw rate resulting from the disturbances and the front steering correction:
wherein the symbol L denotes vehicle wheelbase (i.e., L=a+b) and Ku is the understeer gradient, which is given by:
It follows from Eq. 23 that in order to cancel the effect of disturbances on the vehicle yaw rate, the front steering angle correction must be:
This nominal front steering correction may then be passed through a low pass filter in order to smooth out the command signal and make it more compatible with the dynamics of the steering system. Subsequently, the command may be passed through a high pass filter in order to gradually phase out the steering correction when the disturbance approaches steady-state.
If the vehicle is equipped with an active rear steer system (instead of the active front steer), then following the same approach, the steady-state value of vehicle yaw rate may be:
wherein δf is the front steering angle (due to driver steering) and δr is the rear wheel steering angle correction.
In order to cancel the effect of the disturbances, the rear wheel steering correction may be given by the following equation:
If the vehicle is not equipped with an active steering system, but features a brake-based ESC system, then an automatic brake system intervention can be used instead of steering angle correction. In this case the corrective yaw moment may be generated by applying the difference in braking forces between the left and right sides of the vehicle (i.e., ΔFxLR). During driving in the linear handling range, this brake intervention may impart the yaw moment to vehicle, which may be given by the following equation:
wherein tw denotes the track width.
The vehicle may then be subjected to two yawing moments: Mzs=Fyse=mew (Eq. 1) due to side wind and the corrective moment Mzcor (Eq. 28). Therefore, Eq. 23 for the steady-state value of yaw rate may be modified by replacing the yaw moment mew in the last term by the sum of moments mew−0.5ΔFxLRtw. This yields:
Requiring that the braking correction ΔFxLR cancels the effect of disturbances on the steady state value of yaw rate yields the following value of the brake force difference:
Using Eq. 24, the above equation can be written as:
The difference in braking forces may be expressed in terms of the difference in the wheel circumferential speeds, since the latter is often more convenient to use. In the linear range of tire operation, the longitudinal tire force may be a linear function of tire longitudinal slip:
wherein cx denotes the longitudinal tire stiffness and ΔνLR is the difference in circumferential speeds of left and right wheels. The minus sign suggests that during braking of the left wheels a negative difference between the speed of the left and right wheels may be generated.
It follows from Eq. 32 that the difference in longitudinal tire forces given by Eq. 31 corresponds to the commanded difference between the (circumferential) speeds of left and right wheel ΔνLR:
The relationships described above may form the basis for estimating and rejecting disturbances associated with side wind and changes in road bank angle.
Referring to
The system may determine whether the vehicle is in the linear handling range by evaluating the magnitude of the product of yaw rate and speed, νxΩ, and possibly the measured lateral acceleration, aym. If the values are small enough (e.g., less than about 2 m/s2), the vehicle may be in the linear range of handling. Additional conditions may be considered, for example, the magnitude of the difference between the desired and measured yaw rates should be small enough (e.g., less than about 4 deg/s).
At block 506, the system may determine whether the driver has provided a significant steering correction. This may be accomplished by evaluating the magnitude and the rate of change of the steering angle inputted by the driver. If both are below their respective thresholds, the system may conclude that the driver has not provided a sufficient steering correction. In one aspect, the thresholds may be speed dependent. In another aspect, the magnitude and the rate of change of the desired yaw rate may be used, which may be derived from the steering angle and speed.
Thus, if the driver provides correction, there may be no need for additional system correction and, therefore, the process may end at block 508. However, if the driver does not provide correction, the system may continue to block 510 and may begin the estimation process.
At block 510 the process may determine the estimate of lateral velocity νy using, for example, the measured lateral acceleration ay, yaw rate Ω, front wheel steering angle δf, the rear wheel steering angle δr, estimated vehicle speed νx and/or any other known parameters of vehicle. In one aspect, the lateral velocity may be computed from Eqs. 12 and/or 17.
At block 512, the process may determine the estimate of the lateral force disturbance due to the side wind using, for example, measured lateral acceleration aym, yaw rate Ω, and estimated longitudinal and lateral velocities νx and νy.
The lateral force disturbance due to the side wind per unit mass of vehicle w may be determined from Eq. 13 and, if desired, the yaw moment disturbance due to side wind may be computed from Eq. 1. If the vehicle is equipped with an active rear steer system, then Eq. 13 may be replaced with Eq. 18.
At block 514, the process may determine the estimate of the lateral force disturbance in the form of a gravity component due to the bank angle of the road g sin γ. In one aspect, the process may use the following signals and Eq. 14 to estimate the lateral force disturbance associated with the road bank angle: measured yaw rate Ω, lateral acceleration aym and the estimated vehicle longitudinal and lateral velocities νx and νy. The estimate may be passed through a low pass filter to reduce the effect of noise.
At block 516, the process may determine the total lateral force disturbance estimate (per unit mass) as the sum of the lateral force due to crosswind w and the gravity component due to the bank angle g sin γ.
At block 518, the process may determine whether the magnitude and the rate of change of the lateral disturbances are sufficiently large to warrant the automatic intervention of the control system.
For example, the magnitude of total disturbance estimate, |w+g sin γ|, and/or its rate of change must exceed threshold values. A combination of the magnitude of total disturbance and its rate of change may be considered. If the condition is not satisfied, the steering correction or a brake intervention may not be applied and the process may end at block 520. If the magnitude is sufficient to warrant intervention, the process may proceed to block 522.
At block 522, the process may determine the magnitude and/or the direction of the steering correction or brake intervention necessary to counter the effect of the lateral disturbance.
In one aspect, the nominal front steering correction may be determined from Eq. 25. This nominal front steering correction may then be passed through a low pass filter in order to smooth out the command signal and make it more compatible with the dynamics of the steering system. Subsequently, the command may be passed through a high pass filter in order to gradually phase out the steering correction when the disturbance approaches steady-state. An example of such filter is the one with a transfer function of s/(s+p) where s is the Laplace operand and p=0.3 rad/s, though those skilled in the art will appreciate that this may be achieved in numerous different ways.
If an automatic brake system intervention is used instead of steering angle correction, then the corrective yaw moment may be generated based upon the difference in braking forces between the left and right side of vehicle ΔFxLR. In one aspect, this difference, if expressed in terms of braking forces, may be given by Eq. 31. In another aspect, Eq. 33 may be used with the difference in wheel speeds ΔvLR.
a, 6b, 7a, 7b, 7c and 7d provide example results of the application of the above algorithm to a vehicle in high-fidelity simulation. The vehicle was driven straight at about 120 kph when, at 2 seconds into the simulation, the vehicle entered a crosswind having a speed of about 80 kph. In
Thus, the disclosed system may provide temporary support to the driver when a sudden gust of wind and/or a change in the bank angle of the road causes a significant lateral deviation from the desired path. After a predetermined amount of time, the disclosed system may return control to the driver.
In one aspect, the disclosed system does not require any expensive vision/optical systems. Therefore, the disclosed system may not add any additional hardware beyond what is already available on a vehicle equipped with an ESC system and/or a controlled steering system (e.g., active front steer, active rear steer or electric power steer). The only addition may be software.
Although various aspects of the disclosed estimation and compensation system and method have been shown and described, modifications may occur to those skilled in the art upon reading the specification. The present disclosure includes such modifications and is limited only by the scope of the claims.