# Patent

## Title
Dynamic road gradient estimation

## Abstract

A vehicle and vehicle system are provided with a controller that is configured to generate output indicative of a vehicle mass estimation. The vehicle mass estimation is based on a longitudinal acceleration and a wheel torque when at least one of the longitudinal acceleration, a vehicle speed and a yaw rate indicate an occurrence of a qualified event. The controller is further configured to generate output indicative of a dynamic road gradient estimation based on the vehicle speed, the wheel torque and the vehicle mass estimation.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of U.S. provisional Application No. 61/695,886 filed Aug. 31, 2012, the disclosure of which is incorporated in its entirety by reference herein.


TECHNICAL FIELD
One or more embodiments relate to a vehicle system and method for estimating road gradient and vehicle mass.
BACKGROUND
Road grade and vehicle mass are two parameters that affect a number of vehicle control functions. Such functions include powertrain control, energy management control, transmission control, brake control and stability control. For example, conventional vehicles include an internal combustion engine which provides torque for propelling the vehicle. The engine provides torque through the transmission to the drive wheels, even when the engine is idling. Such torque is referred to as creep-torque. When a vehicle is stopped on an inclined surface, this creep-torque supplements torque applied by a brake system (brake torque) to maintain the vehicle in a stationary position. If the engine were to shutdown, while the vehicle was stopped on the incline, the engine creep-torque would be eliminated and additional brake torque would be needed to maintain the vehicle position and prevent the vehicle from rolling backwards. The amount of additional brake torque needed is based on road grade and vehicle mass.
A micro-hybrid vehicle may enable a stop/start strategy for starting and stopping the vehicle engine during a driving event. The engine is shutdown if no power is required (e.g. while waiting at a traffic light). As soon as power is requested, the engine is automatically restarted. By avoiding unnecessary engine idling, the vehicle's fuel economy will be improved. A micro-hybrid having start/stop functionality may adjust brake pressure to maintain a vehicle position when the engine is shutdown while the vehicle is stopped on an incline. Such a micro-hybrid estimates road grade and vehicle mass, and adjusts brake pressure based on these estimates.
SUMMARY
In one embodiment, a vehicle is provided with a controller that is configured to generate output indicative of a vehicle mass estimation. The vehicle mass estimation is based on a longitudinal acceleration and a wheel torque when at least one of the longitudinal acceleration, a vehicle speed and a yaw rate indicate an occurrence of a qualified event. The vehicle mass estimation is independent of a road gradient estimation.
In another embodiment, a vehicle system is provided with a controller that is configured to generate output indicative of a vehicle mass estimation based on a longitudinal acceleration and a wheel torque in response to a qualified event. The controller is further configured to generate output indicative of a dynamic road gradient estimation based on a vehicle speed, a wheel torque and the vehicle mass estimation.
In yet another embodiment, a method for estimating a dynamic road gradient is provided. Output is generated that is indicative of a vehicle mass estimation based on a longitudinal acceleration and a wheel torque in response to a qualified event. Output is also generated that is indicative of a dynamic road gradient estimation based on a vehicle speed, a wheel torque and the vehicle mass estimation.
As such the vehicle system provides advantages by estimating the vehicle mass independent of the road gradient estimation using an event-seeking strategy. The mass is estimated at vehicle conditions, or qualified events, where the signal to noise ratio of certain inputs is high. The vehicle system then estimates vehicle mass using a recursive least square (RLS) algorithm when the sought event occurs. The vehicle system re-initializes certain parameters in the mass estimation algorithm at vehicle conditions corresponding to potential mass changes (e.g., when the vehicle stops) in order to increase the sensitivity of the algorithm. The vehicle system estimates a dynamic road gradient (RGEdyn) based on the estimated vehicle mass. Generally the dynamic RGE has lower quality conditions as compared to static and kinematic RGE because it involves more signals in its estimation, and the noise and error in each signal will be propagated in the final RGE output. However, the dynamic RGE does not rely on the longitudinal acceleration signal and therefore it has different signal coverage from that of the static and the kinematic estimations. Additionally, the vehicle system estimates the vehicle mass independent of the road gradient estimation, thereby isolating noise in RGEdyn from the vehicle mass estimate.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic diagram of a vehicle having a vehicle system for estimating road gradient and vehicle mass according to one or more embodiments;
FIG. 2 is a side view of the vehicle of FIG. 1 illustrated on a road inclined at a gradient;
FIG. 3 is a bottom view of the vehicle of FIG. 1;
FIG. 4 is a schematic block diagram illustrating a method for arbitrating road gradient estimation, as performed by the vehicle system of FIG. 1, according to one or more embodiments;
FIGS. 5A-5E are time plots illustrating road gradient estimations controlled by the method of FIG. 4;
FIG. 6 is a flow chart illustrating a method for arbitrating road gradient estimations according to one or more embodiments;
FIG. 7 is a time plot illustrating vehicle speed and acceleration;
FIG. 8 is a flow chart illustrating a method for evaluating accelerometer convergence status according to one or more embodiments;
FIG. 9 is a flow chart illustrating a method for locking the static road gradient output according to one or more embodiments;
FIGS. 10A-10C are time plots illustrating vehicle speed, acceleration and road gradient estimations as controlled by the method of FIG. 9;
FIG. 11 a schematic block diagram illustrating a method for kinematic road gradient estimation according to one or more embodiments;
FIG. 12 is a flow chart illustrating a method for estimating kinematic road gradient according to one or more embodiments;
FIG. 13A-13D are time plots illustrating vehicle speed, acceleration and road gradient estimations as controlled by the method of FIG. 12;
FIG. 14 is a flow chart illustrating a method for determining a vehicle driving situation according to one or more embodiments;
FIG. 15 is a flow chart illustrating a method for estimating a vehicle lateral velocity according to one or more embodiments;
FIG. 16 is a flow chart illustrating a method for estimating a body roll angle according to one or more embodiments;
FIG. 17 is a diagram illustrating a method for selecting suspension parameters, according to one or more embodiments;
FIG. 18 is a schematic block diagram illustrating a method for vehicle mass estimation according to one or more embodiments;
FIG. 19 is a flow chart illustrating a method for estimating the vehicle mass according to one or more embodiments;
FIG. 20 is a flow chart illustrating a method for resetting vehicle mass estimation parameters according to one or more embodiments;
FIG. 21 is a flow chart illustrating a method for updating vehicle mass estimation parameters according to one or more embodiments;
FIG. 22 is a time plot illustrating vehicle mass estimations;
FIG. 23 is another time plot illustrating vehicle mass estimations;
FIG. 24 is yet another a time plot illustrating vehicle mass estimations;
FIG. 25 is a flow chart illustrating a method for estimating the quality of an estimator according to one or more embodiments;
FIGS. 26A-26C are time plots illustrating vehicle speed, acceleration and road gradient estimations as controlled by the method of FIG. 25; and
FIGS. 27A-27G are time plots illustrating vehicle speed, and quality factors as controlled by the method of FIG. 25.



DETAILED DESCRIPTION
As required, detailed embodiments of the present invention are disclosed herein; however, it is to be understood that the disclosed embodiments are merely exemplary of the invention that may be embodied in various and alternative forms. The figures are not necessarily to scale; some features may be exaggerated or minimized to show details of particular components. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a representative basis for teaching one skilled in the art to variously employ the present invention.
With reference to FIG. 1, a vehicle system for estimating road gradient and vehicle mass is illustrated in accordance with one or more embodiments and is generally referenced by numeral 10. The vehicle system 10 is depicted within a vehicle 12. The vehicle system 10 includes a controller, such as an engine control module (ECM) 14, an internal combustion engine (ICE) 16 and a vehicle system controller (VSC) 18 that are in communication with each other. The VSC 18 receives input that corresponds to vehicle speed (Vx), acceleration (a), yaw rate (r) and torque (Tpwt, Tbrk), and then provides output that corresponds to a road gradient estimation (RGE) and estimated vehicle mass (Mv) in response to the input.
The illustrated embodiment depicts the vehicle 12 as a micro-hybrid vehicle, which is a vehicle that is propelled by the engine 16, and the engine 16 is repeatedly started and stopped to conserve fuel. An enhanced starter motor 20 is coupled to an engine crankshaft. The starter motor 20 receives electrical power and provides output torque to the crankshaft for starting the engine 16. Although illustrated and described in the context of a micro-hybrid vehicle 12, it is understood that embodiments of the vehicle system 10 may be implemented on other types of vehicles, such as a conventional powertrain vehicle, a battery electric vehicle (BEV), which is powered by an electric motor without assistance of an engine, and a hybrid electric vehicle (HEV), which is powered by an electric motor and an engine.
The vehicle 12 includes a transmission 22 for adjusting the output torque of the engine 16. Torque from the engine 16 is transferred through the transmission 22 to a differential 24 by a transmission output shaft 26. Axle half shafts 28 extend from the differential 24 to a pair of driven wheels 30 to provide drive torque for propelling the vehicle 12.
The vehicle 12 includes a shifter 32 for manually selecting a transmission gear. The shifter 32 includes a sensor (not shown) for providing an output signal that corresponds to a selected transmission gear (e.g., PRNDL). A transmission control module (TCM) 34 communicates with the shifter 32 and the transmission 22 for adjusting the transmission gear ratio based on the shifter selection. Alternatively the shifter 32 may be mechanically connected to the transmission 22 for adjusting the transmission gear ratio.
The vehicle 12 includes a braking system which includes a brake pedal 36, and a booster and a master cylinder which are generally referenced by brake activation block 38 in FIG. 1. The braking system also includes an ABS brake control module 40 that is connected to wheel brake assemblies 42 and the brake activation block 38 by a series of hydraulic lines 44 to effect friction braking. The wheel brake assemblies 42 are oriented at each wheel 30 and may be configured as caliper, or drum brake assemblies. Each wheel brake assembly 42 includes a sensor (not shown) for providing a wheel speed signal (ws). In one or more embodiments, the ABS brake control module 40 provides a calculated vehicle speed signal (Vx) to the VSC 18 based on the wheel speed signals.
The braking system also includes sensors for providing information that corresponds to current brake characteristics. The braking system includes a position switch for providing a brake pedal state (Sbp) signal that corresponds to a brake pedal position (e.g., applied or released). In other embodiments, the braking system includes a position sensor (not shown) for measuring pedal position. The braking system also includes one or more pressure sensors for providing a brake pressure (Pbrk) signal that corresponds to an actual brake pressure value within the brake system (e.g., brake line pressure, master cylinder pressure, or pressure at the wheel assemblies).
The vehicle 12 includes an accelerator pedal 48 with a position sensor for providing an accelerator pedal position (APP) signal that corresponds to a driver request for propulsion. The ECM 14 controls the throttle of the engine 16 based on the APP signal.
The vehicle 12 includes an energy storage device, such as a battery 50. The battery 50 supplies electrical energy to the vehicle controllers and the starter motor 20, as generally indicated by dashed lines in FIG. 1. The vehicle 12 may include a single battery 50, such as a conventional low voltage battery, or multiple batteries, including a high voltage battery. Additionally, the vehicle 12 may include other types of energy storage devices, such as capacitors or fuel cells.
The vehicle 12 also includes one or more inertial sensors 52 which provide output that corresponds to vehicle acceleration (a). In one embodiment, the inertial sensor 52 provides output corresponding to longitudinal acceleration. In another embodiment, the inertial sensors 52 provide output corresponding to longitudinal acceleration, lateral acceleration and vertical acceleration, all of which are generally referenced by the variable “a” in FIG. 1. In one or more embodiments, the inertial sensors 52 also provide output (r) corresponding to a yaw rate. In other embodiments, the sensor 52 is an inclinometer (not shown) which provides output corresponding to a road gradient combined with a vehicle pitch angle.
The vehicle 12 includes a navigation system 54 for receiving destination information from a driver. The navigation system 54 also communicates with a plurality of systems/networks that are external to the vehicle. For example the navigation system 54 may include transceivers for communicating with satellites (e.g. GPS), and computers (via cellular networks, hardwire connections or radio frequency (RF) waves) (not shown). Such networks are described in detail in U.S. application Ser. No. 12/744,092 to Pleet et al., which is hereby incorporated by reference. The navigation system 54 determines a distance between a present location of the vehicle 12 and the destination selected by the driver, and provides this distance, along with speed, acceleration and elevation information associated with the vehicle 12, its target destinations, or other relevant GPS waypoints to the VSC 18.
The VSC 18 communicates with other vehicle systems, sensors and controllers for coordinating their function. As shown in the illustrated embodiment, the VSC 18 receives a plurality of input signals (e.g., Sbp, Pbrk, engine speed (Ne), vehicle speed (Vx), powertrain torque (Tpt), brake torque (Tbt), a, and r) from various vehicle sensors and controllers. Although it is shown as a single controller, the VSC 18 may include multiple controllers that may be used to control multiple vehicle systems according to an overall vehicle control logic, or software. The vehicle controllers, including the VSC 18 generally include any number of microprocessors, ASICs, ICs, memory (e.g., FLASH, ROM, RAM, EPROM and/or EEPROM) and software code to co-act with one another to perform a series of operations. The controllers also include predetermined data, or “look up tables” that are based on calculations and test data and stored within the memory. The VSC 18 communicates with other vehicle systems and controllers (e.g., the ECM 14, the TCM 34, etc.) over one or more wired or wireless vehicle connections using common bus protocols (e.g., CAN and LIN).
With reference to FIGS. 2 and 3, the vehicle 12 is illustrated with a number of forces, accelerations and moments acting thereon. The vehicle 12 is shown with reference to a coordinate system centered about a center of gravity of the vehicle (CG). The coordinate system includes a longitudinal axis (x), a lateral axis (y) and a vertical axis (z). Referring to FIG. 2, the vehicle 12 is positioned upon an inclined surface. The surface has a gradient (road gradient) that is represented by the symbol αr. Additionally, the vehicle pitch angle relative to the road gradient is represented by symbol θ. The weight of the vehicle is represented by Mg, which is illustrated acting upon the CG of the vehicle 12, where M is estimated vehicle mass, and g represents the acceleration due to gravity (9.8 m/s2). The longitudinal velocity of the vehicle 12 is represented by symbol Vx, and is determined based on the wheel speed signals (Ws shown in FIG. 1), according to one or more embodiments. The longitudinal acceleration and the vertical acceleration of the vehicle 12 are represented by variables axs, azs respectively, where the superscript letter(s) indicates that the signal is provided by a sensor, such as the inertial sensors 52 shown in FIG. 1. Referring to FIG. 3 the yaw angular rate acts about the vertical axis (z) and is represented by letter r. The inertial sensor(s) 52 are mounted on the vehicle 12 offset from the CG, and the vehicle system 10 calculates acceleration values at the CG based on the offset. The longitudinal offset between the sensors 52 and the CG is represented by ld, and the lateral offset is represented by lc in FIG. 3.
FIG. 4 depicts a simplified block diagram of a method 410 for road gradient estimation arbitration in accordance with one or more embodiments. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 410 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The vehicle system includes a number of algorithms for estimating road gradient, where each algorithm is suited for different vehicle conditions. The method 410 selects a road gradient estimation algorithm based on current vehicle conditions.
The VSC receives inputs at signal processing block 412. The inputs include vehicle states and sensor signals, (e.g., axs, ays, azs, r) derived or calculated values (e.g., Vx, Tpwt, Tbrk, Fpwt, Fbrk), and quality factors which are all generally represented by input 414. In some embodiments, the vehicle system receives inputs from the navigation system that correspond to vehicle speed and acceleration, which are generally represented by input 418.
At block 420 the VSC estimates road gradient (RGE) using static, kinematic and dynamic estimation algorithms, which are represented by sub-blocks 422, 424 and 426 respectively. The RGE values are provided to an RGE arbitration and diagnostic block 428 which selects which estimation (RGEst, RGEkin, or RGEdyn) to provide as an output (RGEout) for other vehicle systems based on a number of vehicle conditions.
The static RGE algorithm 422 estimates road gradient based on the longitudinal acceleration input (axs). Static RGE is suited for vehicle standstill conditions, such as engine start/stop functionality when a micro-hybrid is stopped on an incline. When the vehicle is stopped on an incline (as shown in FIG. 2) the longitudinal acceleration value (axs) is biased due to deflection within the vehicle suspension system, which is represented by the vehicle pitch angle (θ). The static RGE algorithm compensates for the vehicle pitch angle, as described in detail with respect to FIGS. 7-10C. Generally, the static RGE is the most accurate algorithm among the three RGE estimations since its estimation involves the minimal estimation input and it is the most directly derived estimation results. However, the static RGE is only available at vehicle standstill conditions.
The kinematic RGE algorithm 424 estimates RGE based on input corresponding to the vehicle's kinematic properties including the relationship between vehicle speed, accelerations and yaw rate. The algorithm processes the input using a Kalman Filter as described in detail with respect to FIGS. 11-13D. The kinematic RGE algorithm 424 is suited for normal vehicle motion conditions (e.g., Vx>5 kph) because some of the underlying calculations have low accuracy at low vehicle speeds. Additionally, a Kalman filter is a dynamic filter, a good estimation is usually not available until the estimation state has converged. The transient time depends on the difference between the initial state and the true road grade state as well as the estimation conditions like: signal frequency component richness and signal to noise ratio, etc.
The dynamic RGE algorithm 426 estimates RGE based on input corresponding to vehicle speed, wheel torques (brake torque and powertrain output torque), and other drag forces acting on the vehicle (e.g., aerodynamic drag, rolling resistance and the road gradient load). The vehicle system includes an algorithm for estimating vehicle mass, which is included in block 426, and described in detail with respect to FIGS. 18-24. Generally the dynamic RGE has lower quality conditions as compared to static and kinematic RGE because it involves more signals in its estimation, and the noise and error in each signal will be propagated in the final RGE output. However, the dynamic RGE does not rely on the longitudinal acceleration signal (axs), and therefore it has different signal coverage from that of the static and the kinematic estimations. For example, in the case of sensor error conditions the longitudinal acceleration signal (axs) may not be available, and then the dynamic RGE is the only available estimate.
At block 430 the vehicle system performs driving state monitoring and estimation algorithm adjustment. The estimation algorithm adjustment includes feedforward input compensation for the RGEkin algorithm of block 424, which is described in detail with respect to FIGS. 14-17.
At block 432 the vehicle system evaluates vehicle operation conditions, such as vehicle stability, oscillation and speed, and provides a corresponding judgment to block 428. At block 434 the vehicle system determines a quality factor for each RGE algorithm (QFst, QFkin, QFdyn), along with an overall RGE quality factor (QFRGE). In one or more embodiments the QFRGE is provided with the RGEout to other vehicle systems. Such quality evaluations are described in detail with respect to FIGS. 25-27G. Generally, a quality factor evaluates both the availability and the accuracy of the estimated information. The quality factor is selected from a range of 0-3, where a quality value of three (3) indicates full quality, two (2) indicates some degradation, (1) indicates poor quality, and a value of zero indicates that estimation is not currently available.
FIGS. 5A-5E illustrate the impact of the method 410 for arbitrating road gradient estimations of FIG. 4. FIGS. 5A-5E include 5 graphs of measurement data taken over a common period of time on a road having a known actual road gradient. FIG. 5A illustrates vehicle speed (Vx) and the static road gradient estimation (RGEst). FIG. 5B illustrates vehicle speed (Vx) and the kinematic road gradient estimation (RGEkin). FIG. 5C illustrates vehicle speed (Vx) and the road gradient estimation output (RGEout). FIG. 5D illustrates the road gradient estimation output (RGEout) and the actual road gradient (αr). FIG. 5E illustrates the percent error between RGEout and αr.
FIG. 5A-5C illustrate how the method 410 arbitrates between the static and kinematic road gradient estimations during a drive cycle where the vehicle makes a number of stops. Vehicle speed (Vx) is illustrated on each of FIGS. 5A-5C. The vehicle makes three stops in the illustrated embodiment, which are referenced by numerals 512, 514 and 516.
As depicted in FIG. 5A, the static road gradient estimate (RGEst) is available when the vehicle stops, or is at standstill. At the first stop 512, the RGEst indicates that the road gradient is approximately 20%, as represented by numeral 518. The RGEst indicates that the road gradient is approximately 0% at the second stop, and −12% at the third stop, as indicated by numerals 520 and 522 respectively. When the vehicle is moving, static road gradient estimation is not available, and the RGEst waveform is equal to zero, according to one or more embodiments.
As depicted in FIG. 5B, the kinematic road gradient estimate (RGEkin) is available when the vehicle is moving or during normal vehicle motion. Prior to the first stop 512, the RGEkin indicates that the road gradient increases from approximately 0% to 20%, as represented by numeral 524. At the first stop 512, RGEkin is not available, and RGEkin latches/locks to the last RGEkin value as represented by numeral 526. Between the first stop 512 and the second stop 514, the RGEkin decreases from approximately 20% to 0% as represented by numeral 528. Between the second stop 514 and the third stop 516, the RGEkin decreases from approximately 0% to −12% as represented by numeral 530. After the third stop 516, the RGEkin increases from −12% to 0%, as represented by numeral 532.
As depicted in FIG. 5C, the output road gradient estimate (RGEout) as described with reference to FIG. 4, generally selects RGEst at standstill conditions and RGEkin when the vehicle is moving. FIG. 5D illustrates RGEout along with the actual road gradient (αr). FIG. 5E depicts the percent error (Error) between RGEout and αr. For this example, the error is approximately +/−4%.
Table A depicts a comparison of the inputs and vehicle conditions associated with the algorithms for estimating road gradient, and is shown below:














Estimation Input
Estimation Availability

















Wheel

Very




RGE
Vehicle

Torque

Low
Normal
Convergence


Algorithms
Speed
Acceleration
& Force
Standstill
Speeds
Motion
Phase





Static

X

X


Yes*


Kinematic
X
X



X
Yes


Dynamic
X

X
X

X
Yes









Column 1 of Table A illustrates the three road gradient estimation algorithms: static, kinematic and dynamic. Columns 2-4 depict the primary inputs that each algorithm uses for its estimate. For example, the vehicle system estimates a static road gradient based on vehicle acceleration. The vehicle system estimates a kinematic road gradient based on vehicle speed and acceleration, and the vehicle system estimates a dynamic road gradient based on vehicle speed and wheel torque/force. These inputs may be measured values, or derived values as described with reference to FIG. 1. Additionally, in some embodiments, the vehicle system uses inputs in addition to those depicted in FIG. 5, for the corresponding road gradient estimates.
Columns 5-8 depict the vehicle conditions at which each road grade estimate is available. For example, the static road gradient estimation is available when the vehicle is at standstill. The kinematic road gradient estimation is available when the vehicle is in normal motion, (e.g., vehicle speed is greater than 5 mph). The dynamic road gradient estimation is available when the vehicle is at standstill, or during normal vehicle motion. Although column 6 of the table indicates that none of the estimates are available at very low vehicle speeds, (e.g., less than 5 mph), the vehicle system still provides an accurate road gradient estimate during such conditions, as described below with reference to FIGS. 25-26C.
Column 7 depicts the convergence phase of each estimate. The convergence phase relates to the amount of time it takes the vehicle system to estimate a road gradient using each algorithm during a transition between vehicle conditions (e.g., from standstill to very low speed, to normal motion).
With reference to FIG. 6, a method for arbitrating road gradient estimations is illustrated in accordance with one or more embodiments and is generally referenced by numeral 610. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 610 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. In one or more embodiments, the operations of the method 610 correspond to the operations of method 410 described with reference to FIG. 4. The vehicle system includes a number of algorithms for estimating road gradient, where each algorithm is suited for different vehicle conditions. The method 610 selects a road gradient estimation algorithm based on estimation status, estimation readiness and estimation quality factors.
The estimation status factor for each RGE algorithm (static, kinematic, and dynamic) is defined as a Boolean variable (TRUE/FALSE) and represented by variables (STATUSst, STATUSkin, and STATUSdyn). The estimation status factor indicates whether the estimation condition is satisfied and if the estimation function is updating its state at a corresponding RGE channel (RGEst, RGEkin, and RGEdyn). A status value of TRUE indicates that the estimation method is actively monitoring the vehicle dynamic states (e.g., Vx, axs, ays, azs, and r) and the road gradient information and its estimation channel are updating.
The estimation readiness factor for each RGE algorithm (static, kinematic, and dynamic) is also defined as a Boolean variable (TRUE/FALSE) and represented by variables (READYst, READYkin, and READYdyn). The estimation readiness generally indicates the present usefulness of the RGE function. More specifically, the static estimation readiness (READYst,) indicates whether RGE has already provided a road grade estimate for the present vehicle location. The kinematic and dynamic estimation readiness values (READYkin, and READYdyn) indicate that their estimation state is validated with quality factor greater than 2 (degraded, or full quality) for more than a threshold period of time. In the presence of an unqualified estimation condition, the accumulated time under such condition has not exceeded another time threshold such that the estimation state has not deviated out of the degraded boundary.
The estimation quality factor for the RGE output as well as that for each RGE algorithm (static, kinematic, and dynamic) corresponds to a value between 0-3, and is represented by variables (QFRGE, QFst, QFkin, and QFdyn). The estimation quality factor evaluates both the availability and the accuracy of the estimated road gradient information. QFRGE is the general RGE output signal quality evaluation, which is based on the quality factor of the selected RGE algorithm, and further processed according to the vehicle's dynamic situation after the update event. Such quality evaluation is described in detail with respect to FIGS. 25-26C. For the evaluated quality factor from individual channels, QFst, QFkin, and QFdyn, a value of 3 indicates full quality and a value of 2 indicates that estimation is still available but it is possible that the current value is outside of the specified precision range.
At operation 612, the vehicle system evaluates the current estimation status of the static RGE function to determine if STATUSst is equal to TRUE. If the determination at operation 612 is positive, the vehicle system proceeds to operation 614. At operation 614, the vehicle system evaluates the current estimation readiness of the static RGE function to determine if READYst is equal to TRUE. If the determination at operation 614 is positive, the vehicle system proceeds to operation 616. At operation 616, the vehicle system evaluates the current estimation quality of the static RGE function to determine if QFst is equal to three (full quality). If the determination at operation 616 is positive, the vehicle system proceeds to operation 618, and selects the static road gradient estimation by setting the RGE output equal to RGEst. If the determination at operation 612 is negative, the vehicle system proceeds to operation 620.
At operation 620, the vehicle system evaluates the current estimation status of the kinematic RGE function to determine if STATUSkin is equal to TRUE. If the determination at operation 620 is positive, the vehicle system proceeds to operation 622. At operation 622, the vehicle system evaluates the current estimation readiness of the kinematic RGE function to determine if READYkin is equal to TRUE. If the determination at operation 622 is positive, the vehicle system proceeds to operation 624. At operation 624, the vehicle system evaluates the current estimation quality of the kinetic RGE function to determine if QF1 is equal to three (full quality). If the determination at operation 616 is positive, the vehicle system proceeds to operation 626 and selects the kinematic road gradient estimation by setting the RGE output equal to RGEkin.
The static and kinematic road gradient estimations are preferred over RGEdyn because RGEdyn has lower quality conditions due to the additional signals in its estimation, and the noise in each signal will be propagated in the final RGE output. Therefore even if the static and kinematic quality factors are less than three, the vehicle system may still select RGE or RGEkin over RGEdyn.
If the determination at operation 616 is negative, then QFst is less than 3, and the vehicle system proceeds to operation 628. At operation 628 the vehicle system sets a temporary quality factor value (QFtmp) equal to QFst, and a temporary road gradient estimate (RGEtmp) equal to RGEst. Similarly, if the determination at operation 624 is negative, then QFkin is less than 3, and the vehicle system proceeds to operation 630. At operation 630 the vehicle system sets a temporary quality factor value (QFtmp) equal to QFkin, and a temporary road gradient estimate (RGEtmp) equal to RGEkin. After both operations 628 and 630, the vehicle system proceeds to operation 632 to evaluate the overall RGE quality factor (QFRGE).
At operation 632 the vehicle system determines if QFRGE is less than 3. If the determination is negative, the vehicle system proceeds to operation 634 and does not update RGE. If the overall RGE has full quality, the vehicle system is not going to select a different RGE of lesser quality. If the determination at operation 632 is negative, then the vehicle system proceeds to operation 636.
At operation 636, the vehicle system evaluates the current estimation status and the current estimation readiness of the dynamic RGE function to determine if both STATUSdyn and READYdr, are equal to TRUE. If the determination at operation 636 is negative, the vehicle system proceeds to operation 638. At operation 638 the vehicle system evaluates the temporary quality factor QFtmp to determine if it is greater than or equal to 2. If the determination at operation 638 is negative, the vehicle system proceeds to operation 634 and does not update RGE. If the determination at operation 638 is positive, the vehicle system proceeds to operation 640 and selects the temporary road gradient (static or kinematic) by setting RGE equal to RGEtmp. If the determination at operation 636 is positive, the vehicle system proceeds to operation 642.
At operation 642, the vehicle system evaluates the temporary quality factor (QFtmp) and the dynamic quality factor (QFdyn), to determine if the minimum value of QFtmp and QFdyn is greater than or equal to 2. If the determination at operation 642 is positive, the vehicle system proceeds to operation 644. At operation 644 the vehicle system determines if QFtmp is greater than QFdyn. If the determination at operation 644 is positive, the vehicle system proceeds to operation 640 and selects the temporary road gradient (static or kinematic) by setting RGE equal to RGEtmp. If the determination at operation 644 is negative (QFtmp is not greater than QFdyn), then the vehicle system proceeds to operation 646 and selects the dynamic road gradient by setting RGE equal to RGEdyn.
If the determination at operations 614, 620 or 622 is negative, the vehicle system proceeds to operation 648. At operation 648, the vehicle system evaluates the current estimation status and the current estimation readiness of the dynamic RGE function to determine if both STATUSdyn and READYdyn are equal to TRUE.
If the determination at operation 648 is positive, the vehicle system proceeds to operation 650. At operation 650, the vehicle system evaluates the current estimation quality of the dynamic RGE function to determine if QFdyn is greater than or equal to two. If the determination at operation 650 is positive, the vehicle system proceeds to operation 646 and selects the dynamic road gradient estimation by setting the RGE output equal to RGEdyn. If the determination at operations 648 or 650 are negative, the vehicle system proceeds to operation 652 and does not update the RGE output, by latching the RGE output value to the previous RGE output value. After operations 618, 626, 640, 646, and 652 the vehicle system proceeds to operation 654. At operation 654 the vehicle system provides the RGE output and QF, and then returns to operation 612 for another iteration of the method 610.
As such the vehicle system provides advantages over existing systems by arbitrating or selecting a road gradient estimation from a number of different estimations (static, kinematic, and dynamic). Each road gradient estimation is suited for different vehicle conditions. By arbitrating between these different estimations, the vehicle system provides a road gradient estimation over a wide range of vehicle conditions.
Referring back to FIG. 2, the static and kinematic road grade estimates are based on vehicle kinematic equations. Kinematics generally refers to description of how a body moves in space and time. The vehicle kinematic equations are based on the vehicle acceleration values (axs, ays, azs) that are provided by the inertial sensors (52). The coordinate system (x, y, z) is oriented relative to the road and about the center of gravity (CG) of the vehicle. However, the inertial sensors 52 are mounted offset from the coordinate system (x, y, z), where lc is a lateral offset and ld is a longitudinal offset. Additionally, the inertial sensors 52 provide output signals that are offset from the coordinate system due to the vehicle suspension. The vehicle pitch angle (θ) represents an angular offset from the longitudinal (x) axis. The vehicle roll angle (φ) (not shown) represents an angular offset from the lateral (y) axis, and the vehicle yaw rate (r) represents an angular offset from the vertical (z) axis. The vehicle system estimates the road gradient (αr), and therefore separates the vehicle offsets (e.g., vehicle pitch angle (θ)) from the road gradient (αr).
The governing kinematic vehicle frame equations are represented below:

axs+lc{dot over (r)}=({dot over (V)}x−r*Vy)cos θ+g sin(αr+0)  (1)

ays+ld{dot over (r)}=({dot over (V)}y+r*Vx)cos φv−g sin φ cos(αr+θ)  (2)

azs={dot over (V)}x sin θ+g cos(αr+θ)+naz  (3)

where axs, ays, azs are the vehicle longitudinal, lateral and vertical accelerations. lc is the sensor lateral offset displacement form the vehicle rotational center (CG) in a turning maneuver, and ld is the sensor longitudinal offset displacement. Vx, and Vy, are longitudinal speed and lateral speed; and {dot over (V)}x, {dot over (V)}y, and {dot over (V)}z, are the calculated longitudinal, lateral and vertical accelerations; θ is the vehicle pitch angle and r is the yaw rate. φ is the total vehicle roll angle, which corresponds to the sum of the vehicle body roll angle φv and the road bank angle φb. Additionally, n represents all of the vertical acceleration disturbances caused by vehicle suspension oscillation and road surface fluctuation.

The vehicle system estimates a static road gradient based on the vehicle longitudinal acceleration and equation 1. Since the vehicle is static, or at standstill, the longitudinal velocity (Vy), the longitudinal acceleration ({dot over (V)}x), yaw rate, and the derivative of the yaw rate are set to zero, and the kinematic equations are simplified. Equation 4 illustrates an equation for calculating longitudinal acceleration at static conditions, and is shown below:

axs=g(αr+θ)≈g sin αr+gθ  (4)

The vehicle pitch angle (θ) due to suspension deflection is proportional to the road gradient angle as θ≈Kp (g sin αr) where Kp represents the stiffness of the vehicle body at steady state. Kp is a predetermined value based on the magnitude of axs. In one or more embodiments, the VSC includes stored data corresponding to a lookup table of pre-calibrated data with the value of Kp interpolated with respect to the vector of axs values across its positive and negative operating range. Equation 5 illustrates an intermediate derivation step, and equation 6 provides an equation for calculating the road gradient (αr) at static conditions, as shown below:




  
    
      
        
          a
          x
          s
        
        ≈
        
          
            g
            ⁡
            
              (
              
                1
                +
                
                  
                    K
                    p
                  
                  ·
                  g
                
              
              )
            
          
          ⁢
          sin
          ⁢
          
            
          
          ⁢
          
            α
            r
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          α
          r
        
        =
        
          a
          ⁢
          
            
          
          ⁢
          
            sin
            (
            
              
                a
                x
                s
              
              
                g
                ⁡
                
                  (
                  
                    1
                    +
                    
                      
                        K
                        p
                      
                      ⁢
                      g
                    
                  
                  )
                
              
            
            )
          
        
      
    
    
      
        (
        6
        )
      
    
  




Generally, the static RGE is the most accurate algorithm among the three RGE estimations since it involves the minimal estimation input and it is the most directly derived estimation results. However, the static RGE is only available at vehicle standstill conditions.
FIG. 7 illustrates measurement data taken as the vehicle stops on a grade, according to one or more embodiments. FIG. 7 illustrates input signals vehicle speed (Vx), and longitudinal acceleration (axs), along with road grade (αr) data calculated according to equation 6. At time t1, the vehicle is decelerating as indicated by the decreasing slope of Vx. At time t2, the vehicle stops. After the vehicle stops, there is noise present on the longitudinal acceleration signal (axs) due to vehicle body oscillations, as referenced by numeral 712. This noise is amplified by the road gradient calculation, as illustrated by the road gradient αr waveform, and referenced by numeral 714. At time t3, the longitudinal acceleration signal (axs) has converged enough to provide quality information for road grade estimation, as referenced by numeral 716. The time difference between t2 and t3 is approximately 3 s in the illustrated embodiment. Therefore a vehicle system relying on the road gradient (αr) calculation for RGE, may wait approximately 3 s before using the value. This time delay results in decreased fuel efficiency for a micro-hybrid utilizing engine start/stop functionality.
With reference to FIG. 8, a method for evaluating accelerometer convergence status is illustrated in accordance with one or more embodiments and is generally referenced by numeral 810. The method is implemented using software code contained within the VSC, according to one or more embodiments. While the vehicle system is estimating the static road gradient, its estimation state and quality may be continuously monitored. High quality estimation occurs when the accelerometer output (axs) is converged and (axs) is kept in such a steady state for a minimal amount of time. The accelerometer output (axs) is assigned a convergence status factor (CONst) that indicates whether the raw signal is in steady state or not in three level of label code. Where CONst that is equal to 0, indicates that the convergence is not stable; a CONst equal to 1 indicates that the convergence is somewhat stable; and a CONst equal to 2 indicates that the convergences is very stable.
Additionally, the convergence status signal (CONst) is also used to evaluate the static RGE output quality. For example, in one or more embodiments, a CONst of 2 (very stable), corresponds to a highest or full quality estimation (QFst=3). A CONst of 1 (somewhat stable) corresponds to a degraded quality estimation (QFst=2). A CONst of 0 (not stable) corresponds to a non-qualified estimation (QFst=1), or no data (QFst=0). The convergence status signal (CONst) is based on the derivative of longitudinal acceleration.
At operation 812 the vehicle system evaluates the vehicle speed (Vx) and the convergence status factor (CONst) to determine if both Vx is less than a low speed threshold; and if the accelerometer signal quality is ok (e.g., if CONst=2). If the determination at operation 812 is negative, the vehicle system proceeds to operation 814. At operation 814 the vehicle system evaluates the convergence status factor (CONst) to determine if CONst indicates that the signal is stable (e.g., if CONst=2). If the determination at operation 814 is negative, the vehicle system proceeds to operation 816 and resets the timers and flags by setting a stable timer, a converging timer, a nonstable timer, and the convergence status factor (CONst) all to zero (CONst=Not Stable). The vehicle system also sets a Gradefound flag equal to FALSE, when there is an RGE value available. Gradefound indicates whether a valid estimation (RGEst) has already been provided for a present vehicle stop event. A Gradefound that is equal to TRUE indicates that a valid RGEst has already been provided, and Gradefound that is equal to FALSE indicates that a valid RGEst has not already been provided or is no longer valid. If the determination at operation 812 is positive, the vehicle system proceeds to operation 818.
At operation 818 the vehicle system computes a longitudinal jerk ({dot over (a)}xs), or a first order derivative of the longitudinal acceleration using a low pass filter. At operation 820, the vehicle system compares the jerk ({dot over (a)}xs) to a converging threshold to determine if {dot over (a)}xs is less than the converging threshold. If the determination at operation 820 is negative, then the vehicle system proceeds to operation 822 and the nonstable timer accrues with cycle time Dt. The vehicle system also proceeds to operation 822 if the determination at operation 814 is positive. At operation 824 the vehicle system evaluates the non-stable timer, to determine if the non-stable timer is greater than a non-stable timer threshold. If the determination at operation 824 is positive, the vehicle system proceeds to operation 826 and sets the converging timer, the stable timer and the converging status factor (CONst) all to zero (CONst=Not Stable). If the determination at operation 820 is positive, the vehicle system proceeds to operation 828.
At operation 828 the vehicle system evaluates the longitudinal jerk ({dot over (a)}xs), to determine if {dot over (a)}xs is less than a stabilizing threshold. If the determination at operation 828 is negative, the vehicle system proceeds to operation 830 and the converging timer accrues with cycle time Dt. At operation 832 the vehicle system evaluates the converging timer to determine if the converging timer is greater than the converging time threshold. If the determination at operation 832 is positive, the vehicle system proceeds to operation 834 and resets the nonstable and stable timers by setting them equal to zero; sets the converging status factor to somewhat stable (CONst=1); and sets the Gradefound to TRUE. If the determination at operation 828 is positive, the vehicle system proceeds to operation 836.
At operation 836 the stable timer accrues with cycle time Dt. At operation 838 the vehicle system evaluates the stable timer to determine if the stable timer is greater than a stabilizing time threshold. If the determination at operation 838 is positive, the vehicle system proceeds to operation 840 and resets the nonstable and converging timers by setting them equal to zero; sets the converging status factor to very stable (CONst=2); and sets the Gradefound to TRUE.
With reference to FIG. 9, a method for locking the static road gradient output is illustrated in accordance with one or more embodiments and is generally referenced by numeral 910. The method is implemented using software code contained within the VSC, according to one or more embodiments. Generally once the vehicle system determines that the convergence status is very stable (e.g., CONst=2), then RGEst is locked to avoid any impact from future vehicle disturbances. The lock will be released when the vehicle system detects that the vehicle is moving. A STATIClock signal indicates whether RGEst is locked, unlocked or quarantined. A lock status of quarantined, indicates an intermediate state between lock and unlock where the vehicle system is waiting for additional information.
At operation 912 the vehicle system enters an initialization state that sets the STATIClock signal equal to unlocked. At operation 914, the vehicle system evaluates the STATIClock signal to determine if the static road gradient (RGEst) is locked. If the determination at operation 914 is negative, the vehicle system proceeds to operation 916. At operation 916 the vehicle system determines if the STATIClock signal is quarantined. If the determination at 916 is negative, the vehicle system proceeds to operation 918. At operation 918 the vehicle system determines if the GRADEfound is TRUE. The GRADEfound signal is based on the current RGEout signal. If the determination at operation 918 is positive, the vehicle system proceeds to operation 920. At operation 920 the vehicle system evaluates the convergence status signal to determine if CONst is “very” stable, (e.g., if CONst=2). If the determination at operation 920 is positive, the vehicle system proceeds to operation 922 and sets STATIClock to locked. Then in operation 924 the static road gradient estimate (RGEst) is locked. If the determination at operation 914 is positive (e.g., STATIClock=locked) then the vehicle system proceeds to operation 926.
At operation 926 the vehicle system evaluates the grade found (GRADEfound) to determine if GRADEfound=TRUE. If the determination at operation 926 is positive, the vehicle system proceeds to operation 928 and evaluates the convergence status signal to determine if CONst is not stable, (e.g., if CONst=0). If the determination at operation 928 is positive, the vehicle system proceeds to operation 930 and sets STATIClock to quarantined. After operation 930, and if the determination at operations 928 or 934 are negative, the vehicle system proceeds to operation 924 and locks the static road gradient estimate (RGEst). If the determination at operation 916 is positive, the vehicle system proceeds to operation 932.
At operation 932 the vehicle system evaluates the grade found signal to determine if GRADEfound=TRUE. If the determination at operation 932 is positive, the vehicle system proceeds to operation 934 and evaluates the convergence status signal to determine if CONst is very stable, (e.g., if CONst=2). If the determination at operation 934 is positive, the vehicle system proceeds to operation 922 and sets STATIClock to locked. If the determination at operation 932, 918, 920 or 926 is negative, the vehicle system proceeds to operation 936 sets STATIClock to locked. After operations 936 and 924, the vehicle system proceeds to operation 938, and then returns to operation 914.
With reference to FIGS. 10A-10C, the vehicle system estimates the static road gradient (RGEst) based on the methods 810 and 910 of FIGS. 8 and 9. FIGS. 10A-10C include 3 graphs of measurement data taken over a common period of time and are generally referenced by numeral 1010. Additionally, the vehicle system estimates an initial value for RGEst based on either the kinematic or dynamic road gradient estimates, and then filters the longitudinal acceleration signal (axs) using a low pass filter having a variable bandwidth to determine RGEst, according to one or more embodiments.
FIG. 10A illustrates vehicle speed (Vx) and a signal representing the variable bandwidth of the low pass filter (flt) where the filter has a narrow bandwidth when flt is equal to 1, and a wide bandwidth when flt is equal to zero. FIG. 10B illustrate the measured longitudinal acceleration value (axs). FIG. 10C illustrates the static road grade estimation (RGEst) along with the kinematic road grade estimation (RGEkin). In one or more embodiments the bandwidths are distributed along a frequency range between 2.5 Hz and 0.04 Hz; where the wide bandwidth corresponds to the high end of the frequency range and the narrow bandwidth corresponds to the low end of the frequency range. In other embodiments, the wide bandwidth corresponds to approximately 1 Hz and the narrow bandwidth corresponds to approximately 0.05 Hz.
The vehicle system estimates an initial value for RGEst based on either the kinematic or dynamic road gradient estimates. If no estimates are available, the vehicle system will use a default value. Since the static RGE function only starts after the vehicle speed is lower than a predefined speed threshold, (e.g., 1 kph), the initial value of the estimation algorithm is a factor in the convergence time, or the time delay until RGEst converges to the actual road gradient. Fortunately, the vehicle system includes multiple road gradient algorithms. The other kinematic and dynamic estimation algorithms primarily work in the condition while the vehicle is in motion. Thus, before the start of the static RGE algorithm, information about the road gradient several seconds before the current time or several meters after the vehicle (assuming the vehicle is driving forward) shall be available from the other RGE algorithms and they are the valid reference for the static RGE algorithm to start with. After that, the static RGE algorithm takes over the main estimation task to continue monitoring the road gradient variation while the other two algorithms are paused in absence of qualified estimation conditions. More specifically, the static RGE algorithm, when estimation conditions satisfy, will take the output of either the kinematic or the dynamic estimation algorithm as its initial value depending on which one has the highest quality evaluation. When neither of them has satisfied quality, a default value will be used instead. The qualified static estimation condition is that the vehicle speed will be lower than a speed threshold and such speed threshold will be smaller than the speed thresholds that used to determine on-hold of the other two RGE algorithms.
At time t1, the vehicle starts decelerating as indicated by the decreasing slope of Vx. At time t2, the vehicle stops and the RGE becomes active. Rather than starting from 0%, and converging to the actual road gradient, the RGE starts from the last RGEkin value, as shown at point 1012 where RGEst is equal to RGEkin, and both of which are approximately equal to −22%. Then as RGE accumulates data, it converges to actual road gradient of approximately −20% at time t3 as referenced by numeral 1014.
After the vehicle stops, there is noise present on the longitudinal acceleration signal (axs) due to vehicle body oscillations, as referenced by numeral 1016. The vehicle system filters the longitudinal acceleration signal (axs) using a low pass filter having a variable bandwidth, as referenced by the fit waveform at 1018. The vehicle system computes a filtered longitudinal acceleration value by dynamically filtering the axs data using a time constant selected based on first and higher order derivatives of axs. In one embodiment, the vehicle system compares a longitudinal jerk ({dot over (a)}xs), or a first order derivative of the longitudinal acceleration, to a threshold jerk. If the jerk ({dot over (a)}xs), is greater than the threshold jerk, then a time constant corresponding to a lower bandwidth is selected for the low pass filter. Conversely, if the jerk is less than the threshold jerk, a time constant corresponding to a higher bandwidth is selected for the low pass filter. For example, the threshold jerk may be either the converging threshold 820 or the stabilizing threshold 828 described above with reference to FIG. 8. In the illustrated embodiment, the fit waveform corresponds to the bandwidth, where fit equals 1 when a low bandwidth is applied and flt equals 0 when a high bandwidth is selected.
In another embodiment, a bandwidth is selected from predetermined data (e.g., a lookup table) based on a second order derivative of the longitudinal acceleration äxs. By filtering axs, the road gradient estimation (RGEst) converges relatively quickly. For example, RGEst converges at time t3, as referenced by numeral 1014, and the time delay between t2 and t3 is approximately 1 s.
Referring back to FIG. 7, existing systems for evaluating road gradient may wait for the longitudinal acceleration signal to converge before making a road gradient estimate. Thus such a system would wait approximately 3 seconds before estimating road gradient for the vehicle conditions illustrated in FIG. 7. However, as shown in FIG. 10A-10C, the vehicle system reduces this time delay by using an initial value for RGEst that is based on other estimations, and by dynamically filtering the longitudinal acceleration.
As such the vehicle system provides advantages over existing systems by estimating static road gradient using kinematic equations and dynamic filtering of the longitudinal acceleration. The vehicle system also sets an initial value of the static road gradient estimation equal to a road gradient provided by a kinematic or dynamic road gradient estimation. By dynamically filtering the longitudinal acceleration, and by using another estimation for the initial value, the vehicle system reduces the time delay before a quality static road gradient estimated value may be provided, which in turn allows for the engine to be shutdown earlier in a start/stop application which results in improved fuel economy.
With reference to FIG. 11, a schematic block diagram illustrating a method for kinematic road gradient estimation is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1110. The method 1110 corresponds to the kinematic road gradient estimation portion of the overall RGE arbitration method 410 illustrated in FIG. 4. Block 1112 corresponds to block 412 of FIG. 4. At block 1112 the vehicle system receives inputs including acceleration (axs, ays, azs), yaw rate (r), and vehicle speed (Vx). In one embodiment, the vehicle system only receives the longitudinal acceleration (axs) and vehicle speed. In other embodiments, the vehicle system also receives GPS information (e.g., vehicle speed, acceleration, and elevation). The input may be received directly from sensors described with respect to FIG. 1, or indirectly over the CAN bus.
Blocks 1114 and 1116 correspond to the RGEkin block 424 of FIG. 4. At block 1114, the vehicle system determines extended Kalman Filter (EKF) parameters Yk and uk, which represent a system output and input respectively. The input parameters are based on vehicle speed (Vx), longitudinal acceleration (axs), and an acceleration offset axoff. The parameters may also be based on optional state information for embodiments using GPS information. At block 1116, the vehicle system estimates RGEkin using the EKF and determines the kinematic quality factor (QFkin).
The (axs, azs and r) inputs are provided to an RGEkin input compensation block 1118, which corresponds to block 430 of FIG. 4. At block 1118, the vehicle system determines the acceleration offset αxoff. The acceleration offset is provided as a feed-forward input compensation value for the kinetic road gradient estimation of blocks 114 and 116. The determination of the acceleration offset is described in detail with respect to FIGS. 14-17.
The vehicle system estimates a kinematic road gradient based on the kinematic vehicle frame equations 1-4 and uses an integral type filtering strategy. By ignoring high order angular terms, the vehicle frame kinematic equations may be simplified as shown in equations 7-9 below:




  
    
      
        
          
            a
            x
            s
          
          +
          
            
              l
              c
            
            ⁢
            
              r
              .
            
          
          -
          
            (
            
              
                
                  V
                  .
                
                x
              
              -
              
                r
                *
                
                  V
                  y
                
              
            
            )
          
        
        =
        
          
            
              -
              
                (
                
                  
                    
                      V
                      .
                    
                    x
                  
                  -
                  
                    r
                    *
                    
                      V
                      y
                    
                  
                
                )
              
            
            ⁢
            
              
                θ
                2
              
              2
            
          
          +
          
            g
            ⁡
            
              (
              
                
                  α
                  r
                
                +
                θ
              
              )
            
          
        
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          
            a
            y
            s
          
          +
          
            
              l
              d
            
            ⁢
            
              r
              .
            
          
          -
          
            (
            
              
                
                  V
                  .
                
                y
              
              +
              
                r
                *
                
                  V
                  x
                
              
            
            )
          
        
        =
        
          
            
              -
              
                (
                
                  
                    
                      V
                      .
                    
                    y
                  
                  +
                  
                    r
                    *
                    
                      V
                      x
                    
                  
                
                )
              
            
            ⁢
            
              
                φ
                2
              
              2
            
          
          -
          
            g
            ⁡
            
              (
              
                
                  φ
                  v
                
                +
                
                  φ
                  r
                
              
              )
            
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          
            a
            z
            s
          
          -
          g
        
        =
        
          
            
              
                V
                .
              
              x
            
            ⁢
            θ
          
          -
          
            
              g
              2
            
            ⁢
            
              
                (
                
                  
                    α
                    r
                  
                  +
                  θ
                
                )
              
              2
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




An equation for longitudinal acceleration ({dot over (V)}x) may be obtained by applying small angle approximation to equation 8, and is shown below:




  
    
      
        
          
            
              
                
                  V
                  .
                
                x
              
              =
              
              ⁢
              
                a
                x
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  a
                  x
                  s
                
                +
                
                  
                    l
                    c
                  
                  ⁢
                  
                    r
                    .
                  
                
                +
                
                  r
                  *
                  
                    V
                    y
                  
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  θ
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  
                    α
                    r
                  
                
              
            
          
        
        
          
            
              ≈
              
              ⁢
              
                
                  a
                  x
                  s
                
                +
                
                  
                    l
                    c
                  
                  ⁢
                  
                    r
                    .
                  
                
                +
                
                  r
                  *
                  
                    V
                    y
                  
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  θ
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    α
                    r
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  a
                  x
                  s
                
                +
                
                  a
                  x
                  off
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    α
                    r
                  
                
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




where axoff represents an acceleration offset.

The vehicle system estimates kinematic road gradient using a Kalman filter, according to one or more embodiments. The Kalman filter operates recursively on streams of noisy input data to provide a statistically optimal estimate of the underlying system state. Basically, the Kalman filter is an algorithm which uses a series of measurements that are observed over time and contain random variations (noise), and provides estimates of unknown variables that are more accurate than the measurements alone. A Kalman filter may be thought of as an approach to estimate the state of a “black box”, based on measurable inputs and outputs.
The Kalman filter algorithm estimates a state using a two-step process including a time update prediction step, and a measurement update step. In the prediction step, the Kalman filter provides estimates of the current state variables, along with their uncertainties. The measurement update step is a form of feedback control. Once the Kalman filter observes the next measurement, the estimates are updated using a weighted average, where more weight is given to estimates with higher certainty.
The weights are calculated from the covariance, which is a measure of the estimated uncertainty of the prediction of the system's state. The Kalman filter process is repeated every time step, where the new estimate and its covariance informs the prediction used in the following iteration. The Kalman filter is a recursive estimator. This means that only the estimated state from the previous time step and the current measurement are used to compute the estimate for the current state.
The Kalman Filter provides estimates for linear system models with additive independent white noise in both the transition and the measurement systems. However, most systems are nonlinear. An extended Kalman filter (EKF) is used to provide estimates of nonlinear systems. The EKF adapts techniques including multivariate Taylor Series expansions, to linearize calculations about a working point.
The vehicle system defines an extended state (axg). Since gravity is a constant value (9.8 m/s2), a first order dynamic of the extended state is equal to zero ({dot over (a)}xg=0). The extended state models the vehicle acceleration dynamic that is excited by variation in the road gradient, and the road gradient is relatively high at low frequencies.
The overall system model for estimation is shown below, and is referenced as equation 11:




  
    
      
        
          
        
        ⁢
        
          
            X
            k
          
          =
          
            
              
                f
                ⁡
                
                  (
                  
                    
                      X
                      
                        k
                        -
                        1
                      
                    
                    ,
                    
                      u
                      k
                    
                  
                  )
                
              
              +
              
                
                  ω
                  k
                
                ⁢
                
                  

                
                [
                
                  
                    
                      
                        
                          V
                          x
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          a
                          x
                          g
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                
                ]
              
            
            =
            
              
                
                  
                     
                     
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            
                              
                                V
                                x
                              
                              ⁡
                              
                                (
                                
                                  k
                                  -
                                  1
                                
                                )
                              
                            
                            -
                            
                              
                                
                                  a
                                  x
                                  g
                                
                                ⁡
                                
                                  (
                                  
                                    k
                                    -
                                    1
                                  
                                  )
                                
                              
                              ⁢
                              
                                τ
                                s
                              
                            
                          
                        
                      
                      
                        
                          
                            
                              a
                              x
                              g
                            
                            ⁡
                            
                              (
                              
                                k
                                -
                                1
                              
                              )
                            
                          
                        
                      
                    
                    ]
                  
                
                +
                
                  
                
                ⁢
                
                  
                    [
                    
                      
                        
                          1
                        
                      
                      
                        
                          0
                        
                      
                    
                    ]
                  
                  ⁢
                  
                    (
                    
                      
                        
                          a
                          x
                          s
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                      +
                      
                        
                          a
                          x
                          off
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                    )
                  
                  ⁢
                  
                    τ
                    s
                  
                
                +
                
                  
                    [
                    
                      
                        
                          
                            
                              ω
                              v
                            
                            ⁡
                            
                              (
                              k
                              )
                            
                          
                        
                      
                      
                        
                          
                            
                              ω
                              a
                            
                            ⁡
                            
                              (
                              k
                              )
                            
                          
                        
                      
                    
                    ]
                  
                  ⁢
                  
                    

                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    X
                    k
                  
                
              
              =
              
                
                  
                    
                      [
                      
                        
                          
                            V
                            x
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                        ⁢
                        
                          
                            a
                            x
                            g
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                      
                      ]
                    
                    T
                  
                  ⁢
                  
                    

                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    u
                    k
                  
                
                =
                
                  
                    
                      (
                      
                        
                          
                            a
                            x
                            s
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                        +
                        
                          
                            a
                            x
                            off
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                      
                      )
                    
                    ⁢
                    
                      τ
                      s
                    
                    ⁢
                    
                      

                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      Y
                      k
                    
                  
                  =
                  
                    
                      
                        [
                        
                          1
                          ⁢
                          
                            
                          
                          ⁢
                          0
                        
                        ]
                      
                      ⁢
                      
                        X
                        k
                      
                    
                    =
                    
                      
                        
                          V
                          x
                        
                        ⁢
                        
                          

                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          ω
                          k
                        
                      
                      =
                      
                        
                          [
                          
                            
                              
                                ω
                                v
                              
                              ⁡
                              
                                (
                                k
                                )
                              
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              
                                ω
                                a
                              
                              ⁡
                              
                                (
                                k
                                )
                              
                            
                          
                          ]
                        
                        T
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




where Xk is the state at time step (k), and Xk−1 is the state at a previous time step. uk is a system input and Yk is a system output. ωk is the process noise, ωv(k) is the process noise associated with the velocity, and ωa(k) is the process noise associated with acceleration. Vx is the longitudinal speed. axg (k) is the extended state acceleration due to gravity, axs(k) is the measured longitudinal acceleration, and axoff (k) is the acceleration offset. τs is the sampling time for the above discrete time implementation model, and T indicates a matrix transform.

Based on the system model shown in equations 11, an observer is designed to estimate the extended state (axg), and the kinematic road gradient is estimated according to equation 12 as shown below:




  
    
      
        
          
            α
            ^
          
          r
        
        =
        
          
            arcsin
            ⁡
            
              (
              
                
                  a
                  x
                  g
                
                g
              
              )
            
          
          =
          
            RGE
            kin
          
        
      
    
    
      
        (
        12
        )
      
    
  




The observer is designed using a standard Extended Kalman Filter (EKF) system model as represented by the function shown above in equation 11, and equation 13, both of which are shown below:

Xk=f(Xk−1,uk)+ωk 

Yk=CXk+ek  (13)

where C is a matrix, which is a relatively simple matrix in equation 11, and ek is an unknown noise.

An EKF is used for nonlinear systems. The nonlinear system model is linearized around the current state at every time step, and provides a transition matrix Fk. The obtained transition matrix Fk is then used to complete the steps of the standard EKF recursions. These recursions are described by two update steps: a prediction step and a measurement update step. In the prediction step, the system input is used to predict the future state of the system ({circumflex over (X)}k|k−1) as shown in equations 14 and 15 below:

{circumflex over (X)}k|k−1=f({circumflex over (X)}k−1|k−1,uk)  (14)

Pk|k−1=FkPk−1|k−1FkT+Qk  (15)

where Pk|k−1 is the estimated error covariance and Qk is the process noise covariance which is equal to E[ωk2]. Fk is the Jacobian and is provided by equation 16 as shown below:





  
    
      
        
          F
          k
        
        =
        
          
            
              ∂
              f
            
            
              ∂
              x
            
          
          ❘
          
            
              x
              ^
            
            
              
                k
                -
                1
              
              ❘
              
                k
                -
                1
              
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




A measurement update step is provided by equations 17-19 as shown below:

Kk=Pk|k−1CT(CPk|k−1CT+Rk)−1  (17)

{circumflex over (X)}k|k={circumflex over (X)}k|k−1+Kk(Yk−C{circumflex over (X)}k|k−1)  (18)

Pk|k=(I−KkC)Pk|k−1  (19)

where Kk is the EKF gain, and Rk is the measurement noise covariance which is equal to E[ek2].

Generally with an EKF system, the process noise covariance (Q) and the measurement noise covariance (R) are not known a priori, or independent of experience. However, well conditioned signals, such as Kk and Pk|k will converge to values that are dependent on Q and R. Therefore, the converged value of Kk may be obtained from predetermined data and it will be used directly for real-time estimation.
In one or more embodiments, the vehicle system includes vertical displacement in its estimation of kinematic road grade. Vertical displacement (Z) may be provided by the navigation system of FIG. 1 based on GPS information. The governing vehicle vertical kinetic equation in the global inertia frame is calculated according to equation 20 as shown below:

Ż=Vx sin αr  (20)

The EKF system model of equation (11) is modified to include vertical displacement Z(k), as shown in equation 21 below:




  
    
      
        
          [
          
            
              
                
                  
                    V
                    x
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                
                  Z
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                
                  
                    a
                    x
                    g
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                
                  
                    a
                    x
                    θ
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        
                          V
                          x
                        
                        ⁡
                        
                          (
                          
                            k
                            -
                            1
                          
                          )
                        
                      
                      -
                      
                        
                          
                            a
                            x
                            g
                          
                          ⁡
                          
                            (
                            
                              k
                              -
                              1
                            
                            )
                          
                        
                        ⁢
                        
                          T
                          s
                        
                      
                      -
                      
                        
                          
                            a
                            x
                            θ
                          
                          ⁡
                          
                            (
                            
                              k
                              -
                              1
                            
                            )
                          
                        
                        ⁢
                        
                          τ
                          s
                        
                      
                    
                  
                
                
                  
                    
                      
                        Z
                        ⁡
                        
                          (
                          
                            k
                            -
                            1
                          
                          )
                        
                      
                      +
                      
                        
                          
                            V
                            x
                          
                          ⁡
                          
                            (
                            
                              K
                              -
                              1
                            
                            )
                          
                        
                        ⁢
                        
                          1
                          g
                        
                        ⁢
                        
                          
                            a
                            x
                            g
                          
                          ⁡
                          
                            (
                            
                              k
                              -
                              1
                            
                            )
                          
                        
                      
                    
                  
                
                
                  
                    
                      
                        a
                        x
                        g
                      
                      ⁡
                      
                        (
                        
                          k
                          -
                          1
                        
                        )
                      
                    
                  
                
                
                  
                    
                      
                        a
                        x
                        θ
                      
                      ⁡
                      
                        (
                        
                          k
                          -
                          1
                        
                        )
                      
                    
                  
                
              
              ]
            
            +
            
              
                [
                
                  
                
                ⁢
                
                  
                    
                      1
                    
                  
                  
                    
                      0
                    
                  
                  
                    
                      0
                    
                  
                  
                    
                      0
                    
                  
                
                ]
              
              ⁢
              
                (
                
                  
                    
                      a
                      x
                      s
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                  +
                  
                    
                      a
                      x
                      cmp
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                
                )
              
              ⁢
              
                τ
                s
              
            
            +
            
              
                [
                
                  
                    
                      
                        
                          ω
                          v
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          ω
                          z
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          ω
                          ∝
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          ω
                          θ
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                  
                
                ]
              
              ⁢
              
                

              
              ⁢
              
                
              
              ⁢
              
                X
                k
              
            
          
          =
          
            
              
                
                  [
                  
                    
                      
                        V
                        x
                      
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      Z
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        a
                        x
                        g
                      
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        a
                        x
                        θ
                      
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                  
                  ]
                
                T
              
              ⁢
              
                

              
              ⁢
              
                
              
              ⁢
              
                u
                k
              
            
            =
            
              
                
                  (
                  
                    
                      
                        a
                        x
                        s
                      
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                    +
                    
                      
                        a
                        x
                        cmp
                      
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  T
                  s
                
                ⁢
                
                  

                
                ⁢
                
                  
                
                ⁢
                
                  Y
                  k
                
              
              =
              
                
                  CX
                  k
                
                =
                
                  
                    
                      [
                      
                        
                          
                            1
                          
                          
                            0
                          
                          
                            0
                          
                          
                            0
                          
                        
                        
                          
                            0
                          
                          
                            1
                          
                          
                            0
                          
                          
                            0
                          
                        
                      
                      ]
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      X
                      k
                    
                  
                  =
                  
                    
                      [
                      
                        
                          
                            V
                            k
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          Z
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                      
                      ]
                    
                    T
                  
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




where axcmp represents an GPS based acceleration offset, and axcmp is equal to (lc{dot over (r)}+r*Vy). axcmp is provided as feedforward input compensation in embodiments using the GPS information. Equation 21 differs from equation 11, in that equation 21 includes two extended states: axg (k) and axθ(k), where axθ(k) is the acceleration due to the vehicle pitch angle.

With reference to FIG. 12, a method for estimating kinematic road gradient is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1210. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 1210 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers.
At operation 1212 the vehicle system initializes estimator parameter values. The estimator parameter values include Xk, uk, and Yk. The estimator parameter values are initialized by setting them equal to default values. At operation 1214 the vehicle system determines if the kinematic estimator status (STATUSkin) is on. The estimation status factor is defined as a Boolean variable (TRUE/FALSE), where TRUE indicates that the status is on, and FALSE indicates that the status is not on. A status value of TRUE indicates that the estimation method is actively monitoring the vehicle dynamic states (e.g., Vx, axs, ays, azs, and r) and the road gradient information and its estimation channel are updating. If the determination at operation 1214 is positive, the vehicle system proceeds to operation 1216.
At operation 1216 the vehicle system evaluates the present vehicle speed (Vx) to determine if Vx is greater than or equal to a threshold vehicle entry speed (Vspdentry). In one embodiment, the threshold vehicle entry speed is approximately 7 kph (4.3 mph). If the determination at operation 1216 is negative, the vehicle system proceeds to operation 1218. At operation 1218 the vehicle system evaluates Vx to determine if it is less than a threshold vehicle exit speed (Vspdexit). In one embodiment, the threshold vehicle exit speed is approximately 3 kph (1.9 mph). If the determination at operation 1218 is positive, then the vehicle system proceeds to operation 1220.
At operation 1220, the vehicle system sets the kinematic estimation status (STATUSkin) to pause. By setting STATUSkin to pause, the vehicle system stops the estimation and retains the current RGEkin value, rather than resetting the value to a default value. At operation 1222 the vehicle system evaluates the present speed (Vx) to determine if Vx is less than a threshold vehicle standstill speed (Vspdstand). The vehicle standstill speed corresponds to the very low vehicle speed referenced in column 6 of Table A. In one embodiment, the threshold vehicle standstill speed is approximately 1 kph (0.6 mph). If the determination at operation 1222 is positive, the vehicle system proceeds to operation 1224 and sets the kinematic estimation status (STATUSkin) to off (FALSE). If the determination at operation 1216 is positive, or the determination at operation 1218 is negative, then the vehicle system proceeds to operation 1226.
Operations 1226-1236 relate to the extended Kalman filter operations, which are described above with reference to equations 11-21. First, at operation 1226 the vehicle system computes a compensated estimator input (uk). At operation 1228 the vehicle system obtains a new system output (Yk). At operation 1230 the vehicle system updates the EKF state (Xk) based on output prediction error, as described above with reference to equations 14-16. At operation 1232 the vehicle system updates the EKF parameter values (Xk, uk, and Yk). At operation 1234, the vehicle system predicts the next system dynamic output Yk. Then at operation 1236 the vehicle system provides the kinematic road gradient estimation.
If the determination at operation 1214 is negative, the vehicle system proceeds to operation 1238. At operation 1238 the vehicle system evaluates the present vehicle speed (Vx) to determine if Vx is greater than or equal to the threshold vehicle entry speed (Vspdentry), which is similar to the determination at operation 1216. In one embodiment, the threshold vehicle entry speed is approximately 7 kph (4.3 mph). If the determination at operation 1238 is positive, the vehicle system proceeds to operation 1240 and sets the kinematic estimation status (STATUSkin) to on (TRUE). After operation 1240 the vehicle system proceeds to operation 1226 to begin the EKF operations.
After operations 1224 and 1244, and after negative determinations at operations 1222 and 1238, the vehicle system proceeds to operation 1242. After operation 1242 the vehicle system returns to operation 1214 for another iteration.
FIGS. 13A-13D are graphs illustrating the impact of the method 1210 for estimating kinematic road gradient when a vehicle stops on an incline. FIGS. 13A-13D include 4 graphs of measurement data taken over a common period of time and are generally referenced by numeral 1310. FIG. 13A illustrates brake pressure (Pbrk) which is measured within the master cylinder according to one embodiment. In other embodiments, Pbrk is measured within the brake lines or at the wheel assemblies. FIG. 13B illustrates vehicle speed (Vx). FIG. 13C illustrates the gear selection (PRNDL), where a value of 3 corresponds to the drive (D) gear selection. FIG. 13C also illustrates the state of the brake pedal (Sbp), where a positive value indicates that the brake is applied, and a value of zero indicates that the brake is released. FIG. 13D illustrates the measured longitudinal acceleration value (axs), along with the static road grade estimation (RGEst) and the kinematic road grade estimation (RGEkin).
At time t1, the vehicle is accelerating, as indicated by the increasing slope of Vx. At time t2 the vehicle starts decelerating as indicated by the decreasing slope of Vx. At time t3, the vehicle stops. At time t4 the vehicle begins accelerating again, and at time t5 the vehicle begins decelerating again. The kinematic road grade estimation (RGEkin) converges relatively quickly, as referenced by numerals 1312 and 1314.
As described with respect to operations 1222 and 1224, the kinematic estimation status (STATUSkin) is set to off, if the vehicle speed decreases below a vehicle standstill threshold speed (e.g., 1 kph). Such conditions are illustrated by RGEkin and referenced by numeral 1320. After t3, the estimation status is set to off as depicted by the lack of any perceptible noise on RGEkin at 1320. Additionally once the vehicle stops, the status of the static road gradient estimation (STATUSst) and the RGEst converges as referenced by numeral 1322. As shown at 1320, even though the kinematic status is off, RGEkin is still locked at an accurate road grade estimate and both RGEst and RGEkin are providing generally equal estimates of 17.56% grade.
Existing methods for estimating road gradient using a kinematic approach use a derivative filtering strategy. For example, some existing methods estimate road gradient based on the following vehicle longitudinal kinematic equation:

{dot over (u)}=ax+rv−g sin θ  (22)

where {dot over (u)} is the derivative of vehicle longitudinal speed, a, is the vehicle longitudinal acceleration, r is the yaw rate, v is the vehicle lateral speed. However, such an approach amplifies any noise on the signal due to the derivative calculation of {dot over (u)}.

As such the vehicle system provides advantages over the existing methods by estimating kinematic road gradient using an EKF, which provides an integral type filtering strategy. Such an integrating strategy does not amplify noise on the signals, and therefore provides improved estimation accuracy with relatively fast convergence.
The vehicle system estimates kinematic road gradient using a feedforward input compensation corresponding to an acceleration offset (axoff). The acceleration offset is calculated according to equation 23, as shown below:

axoff=lc{dot over (r)}+rVy−g sin θ  (23)

where lc is the lateral offset, r is the yaw rate, Vy is the lateral velocity and θ is the vehicle pitch angle. In one or more embodiments, the vehicle system does not receive these values (lc, Vy, and θ) as inputs. Rather the vehicle system estimates these values from other inputs using a number of compensation strategies. In one or more embodiments, the vehicle system estimates the lateral velocity (Vy) term of equation 23 using a Coriolis Effect compensation strategy, the lateral offset (lc) term using a sensor lateral position offset effect compensation strategy, and the vehicle pitch angle (θ) using a vehicle body relative pitch effect compensation strategy.

The vehicle system estimates the rVy term of equation 23 using a Coriolis Effect compensation strategy which relates to vehicle turning maneuvers. The inertial sensors 52 of FIG. 1, provide the yaw rate, however the lateral velocity (Vy) is estimated by the vehicle system. Two computation methods are provided for estimating the lateral velocity (Vy): a first method, which is an estimation based method, and a second method, which is a derivation based method. The two methods are provided based on the lateral dynamic equation 8 and a vehicle cornering stiffness parameter, respectively. Furthermore, the two methods are integrated to provide the estimated lateral velocity (Vy) under continuous varying vehicle driving conditions.
The first method (the estimation based method) for estimating lateral velocity (Vy) is obtained by applying small angle approximation to Equation 8, and is based on the measured lateral acceleration (ays). Lateral velocity (Vy) is calculated according to equation 24 as shown below:

{dot over (V)}y=ays−ld{dot over (r)}−r*Vx+gφv+gφr  (24)

The road bank angle (φr) is unknown. Therefore it is treated as a disturbance and removed from the calculation. The vehicle body roll angle φv may be approximated by the measured lateral acceleration as shown in equation 25:

φv≈Kφ(ays−ld{dot over (r)})  (25)

where Kφ represents the vehicle body roll stiffness.

Equation 24 may be simplified as shown in equation 26 below:

{dot over (V)}y=(1+KØg)(ays−ld{dot over (r)})−r*Vx=Kay(ays−ld{dot over (r)})−r*Vx  (26)

By integrating both sides of equation 26, an equation is obtained for calculating the estimated lateral velocity ({circumflex over (V)}y(k), as shown by equation 27:

{circumflex over (V)}y(k)={circumflex over (V)}y(k−1)+[Kay(ays−ld{dot over (r)})−r*Vx]*τs  (27)

The measured lateral acceleration (ays) is generally subjected to acceleration disturbances and bank angle impact. Therefore the Vy estimation method described in equations 24-27 is suited for short time intervals of vehicle dynamic maneuvers in which large magnitudes of lateral acceleration (ays) and yaw rate (r) are present, due to the low frequency biased acceleration components.
Unlike the first method, the second method (the derivation based method) for estimating the lateral velocity (Vy) does not depend on integration. A steady state gain (Gr2β) from the vehicle yaw rate (r) to a vehicle side slip angle (β) may be obtained from a linear vehicle model as shown in equation 28:




  
    
      
        
          G
          
            r2
            ⁢
            
              
            
            ⁢
            β
          
        
        =
        
          
            
              l
              r
            
            
              V
              x
            
          
          -
          
            
              
                M
                r
              
              
                C
                
                  α
                  ⁢
                  
                    
                  
                  ⁢
                  r
                
              
            
            ⁢
            
              V
              x
            
          
        
      
    
    
      
        (
        28
        )
      
    
  




where Mr is the vehicle mass distributed to the rear axle and Cαr is a lumped rear axle cornering stiffness.

As the lateral velocity (Vy) becomes approximately equal to the product of the longitudinal velocity (Vx) and the vehicle side slip angle (β), an equation for calculating the lateral velocity may be obtained as shown below:

Vy=VxGr2βr=(lr−kUVx2)*r=Kr(Vx)*r  (29)

where kU is the cornering stiffness and is equal to Mr divided by Cαr.

Therefore the lateral velocity may be estimated by equation 30 as shown below:

{circumflex over (V)}y=r*(lr−kUVx2)=Kr(Vx)*r  (30)

A Coriolis effect compensation equation is provided below:




  
    
      
        
          r
          *
          
            
              V
              ^
            
            y
          
        
        =
        
          
            
              r
              2
            
            *
            
              (
              
                
                  l
                  r
                
                -
                
                  
                    
                      M
                      r
                    
                    
                      C
                      
                        α
                        ⁢
                        
                          
                        
                        ⁢
                        r
                      
                    
                  
                  ⁢
                  
                    V
                    x
                    2
                  
                
              
              )
            
          
          =
          
            
              r
              2
            
            *
            
              (
              
                
                  l
                  r
                
                -
                
                  
                    k
                    U
                  
                  ⁢
                  
                    V
                    x
                    2
                  
                
              
              )
            
          
        
      
    
    
      
        (
        31
        )
      
    
  




The second method as described with reference to equations 28-31 is suited for steady state vehicle turning maneuvers where the side slip angle (β) is small.
The vehicle system implements a comprehensive RGEkin structure that incorporates both the first method (estimation based) and the second method (derivation based). Then, the vehicle system selects which of the two methods to apply, based on vehicle dynamic conditions reflected by the yaw rate (r). The general design concept is to use the derivation based Vy from equation 30, in primarily forward driving and steady state turning situations but to switch to the estimation based Vy from equation 27, in the presence of large vehicle lateral dynamics (speed and acceleration). Since the estimation method is not expected to be carried out for a long period of time, the estimated Vy must be reset to the derived state whenever a steady state condition may be observed during a continuous vehicle handling maneuver. The estimated Vy may be further bounded by predetermined upper and lower limits. The output Vy is provided by a low pass filter.
With reference to FIG. 14, a method for determining a vehicle driving situation is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1410. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 1410 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers.
The method 1410 includes evaluating an indexed lateral dynamic (ILD), which correspond to a present vehicle driving situation. The ILD is defined as:

ILD=|r|+a*|{dot over (r)}|  (32)

where the variable (a) is a scaling parameter. The vehicle system compares the ILD to two thresholds values: a steady state threshold (SST) corresponding to the derivation method, and a dynamic handling threshold (DHT) corresponding to the estimation method, where (SST<DHT). The determination is further made together with two timers, a steady state timer (Tss) and a dynamic driving timer (Tdd) that count the duration a vehicle has stayed in each situation. Function LPF(x) indicates a general low pass filter function. A moving average filter is an example of a low pass filter that may be used to smooth out discontinuous output transition effects.

In operation 1412 the vehicle system determines if the vehicle is in the steady state driving situation. The vehicle system receives a flag (SSDflag) that is associated with the driving situation. A SSDflag that is equal to one (1) indicates a present steady state driving situation, and SSDflag that is equal to two (2) indicates a present dynamic handling situation. Thus, at operation 1412, the system evaluates the SSDflag to determine if SSDflag=2. If the determination at operation 1412 is positive, the vehicle system proceeds to operation 1414.
At operation 1414, the vehicle system evaluates the indexed lateral dynamic (ILD) to determine if the ILD is less than a steady state threshold (SST). The SST corresponds to the derivation method for determining the lateral velocity (Vy) as described above with reference to equations 28-31. If the determination at operation 1414 is negative, the vehicle system proceeds to operation 1416 and resets a steady state timer (Tss) by setting T equal to zero (0). If the determination at operation 1414 is positive, the vehicle system proceeds to operation 1418.
At operation 1418 the vehicle system increments the steady state timer (Tss) by counting up by one execution cycle time (ts). At operation 1420 the vehicle system evaluates the steady state timer (Tss) to determine if Tss is greater than or equal to a settling time threshold (tsssettle). If the determination at operation 1420 is positive, the vehicle system proceeds to operation 1422 and resets the dynamic driving timer (Tdd). Then, at operation 1424 the vehicle system sets the vehicle driving situation to steady state by setting the SSDflag equal to one (1).
If the determination at operation 1412 is negative (the vehicle is not driving in steady state conditions), the vehicle system proceeds to operation 1426. At operation 1426 the vehicle system evaluates the indexed lateral dynamic (ILD) to determine if the ILD is greater than the dynamic handling threshold (DHT). The DHT corresponds to the estimation method for determining the lateral velocity (Vy) as described above with reference to equations 24-27. If the determination at operation 1426 is negative, the vehicle system proceeds to operation 1428 and resets a dynamic driving timer (Tdd) by setting Tdd equal to zero (0). If the determination at operation 1426 is positive, the vehicle system proceeds to operation 1430.
At operation 1430 the vehicle system increments the dynamic driving timer (Tdd) by counting up by one execution cycle time (ts). At operation 1432 the vehicle system evaluates the dynamic driving timer (Tdd) to determine if Tdd is greater than or equal to an excitation time threshold (tddexcite). If the determination at operation 1432 is positive, the vehicle system proceeds to operation 1434 and resets the steady state timer (Tss). Then, at operation 1436 the vehicle system sets the vehicle driving situation to dynamic driving by setting the SSDflag equal to two (2). After operations 1416, 1424, 1428 and 1436, and after a negative determination at operations 1420 or 1432, the vehicle system proceeds to operation 1438 and then returns to operation 1412 for another iteration.
With reference to FIG. 15, a method for estimating a vehicle lateral velocity (Vy) based on a present vehicle driving situation is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1510. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 1510 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The method 1510 arbitrates, by selecting a method for estimating lateral velocity (Vy), from one of the two methods described above with reference to equations 24-31; either the integration estimation based method or the derivation estimation based method
In operation 1512 the vehicle system determines if the vehicle is in the dynamic driving situation. The vehicle system evaluates the SSDflag to determine if SSDflag=2. If the determination at operation 1512 is positive, the vehicle system proceeds to operation 1514, and estimates the lateral velocity (i) using the integration estimation provided in equation 27. After operation 1514 the vehicle system proceeds to operation 1516. At operation 1516 the vehicle system filters the estimated lateral velocity (i) using a low pass filter (LPF), and updates a lateral velocity output (Vy). If the determination at operation 1512 is negative, the vehicle system proceeds to operation 1518. At operation 1518 the vehicle system estimates the lateral velocity (i) using the derivation estimation provided in equation 30. After operation 1518, the vehicle system proceeds to operation 1516, filters the estimated lateral velocity ({circumflex over (V)}y) using a low pass filter (LPF), and updates the lateral velocity output (Vy) based on the derivation estimation. After operation 1516, the vehicle system proceeds to operation 1520 and provides the updated lateral velocity (Vy) output, and then returns to operation 1512.
In one or more embodiments, the vehicle system estimates the lateral offset (lc) term using a sensor lateral position offset effect compensation strategy. The lateral offset (lc) is defined by the following equation:

lc=lcst+htφv  (33)

where lcst is a static sensor installation offset. The variable htφv is a dynamic vehicle body roll component, where ht is the distance from the sprung mass C.G. to the roll center. The vehicle roll angle φv may be obtained from a roll dynamic model or from suspension deflections (in embodiments having suspension deflection sensors). In one or more embodiments, φv at steady state is derived from the yaw rate (r) at steady state as shown below:

φv=Gr2φr  (34)

Combining equations 33 and 34 and multiplying each side by the derivative of yaw rate ({dot over (r)}) provides the following equation:

lc{dot over (r)}=lcst{dot over (r)}+htGr2φ*r*{dot over (r)}  (35)

In a dynamic driving situation, the body roll angle (φv) is approximated from the lateral acceleration (ays) according to equation 36 as shown below:

φv=Gay2φ(ays−ld{dot over (r)})  (36)

Combining equations 35 and 36 provides the following equation:

lc{dot over (r)}=(lcst=htGay2φld){dot over (r)}+htGay2φays  (37)

where Gay2φ is equal to Kphi.

With reference to FIG. 16, a method for estimating the body roll angle (φv) based on a present vehicle driving situation is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1610. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 1610 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The method 1610 arbitrates, by selecting a method for estimating the body roll angle (φv), from one of the two methods described above with reference to equations 33-37; either using lateral acceleration or yaw rate.
In operation 1612 the vehicle system determines if the vehicle is in the dynamic driving situation. The vehicle system evaluates the SSDflag to determine if SSDflag=2. If the determination at operation 1612 is positive, the vehicle system proceeds to operation 1614, and estimates the body roll angle ({circumflex over (φ)}v) based on lateral acceleration as shown by equation 36. After operation 1614 the vehicle system proceeds to operation 1616. At operation 1616 the vehicle system filters the estimated body roll angle ({circumflex over (φ)}v) using a low pass filter (LPF), and updates a body roll angle output (φv). If the determination at operation 1612 is negative, the vehicle system proceeds to operation 1618. At operation 1618 the vehicle system estimates the body roll angle ({circumflex over (φ)}v) based on the yaw rate as provided in equation 34. After operation 1618, the vehicle system proceeds to operation 1616, filters the estimated body roll angle ({circumflex over (φ)}v) using a low pass filter (LPF), and updates the body roll angle output ({circumflex over (φ)}v). After operation 1616, the vehicle system proceeds to operation 1620 and provides the updated body roll angle output ({circumflex over (φ)}v), and then returns to operation 1612.
In one or more embodiments, the vehicle system estimates the vehicle body pitch angle (θ) using one of three compensation strategies. The different strategies correspond to the availability of sensors and accuracy. For example, a first strategy estimates the pitch angle (θ) based on the vehicle longitudinal speed (Vx) and the vehicle longitudinal acceleration (axs). A second strategy, estimates the pitch angle (θ) based on Vx, axs, and a vertical acceleration azs. A third strategy estimates the pitch angle (θ) based on Vx, axs, azs and an estimated acceleration due to gravity component axg provided by the extended Kalman filter. Generally, the more input signals used and the more accurate the source information is, the better the compensation strategy is at cancelling the coupling effect of the axg estimation, and the better the overall road gradient estimation.
The first compensation strategy for estimating the vehicle body pitch angle (θ), is referenced as dynamic pitch compensation. For this strategy, the vehicle system estimates the pitch angle (θ) based on the vehicle longitudinal speed (Vx) and the vehicle longitudinal acceleration (axs). A low frequency component of θ may be approximated by a first or second order linear dynamic system. The thrust to pitch vehicle dynamic is generally caused by the longitudinal forces exerted on the vehicle body according to the following equation:




  
    
      
        
          θ
          ∝
          
            
              K
              p
            
            ⁢
            
              a
              x
              r
            
          
        
        =
        
          
            
              
                
                  M
                  s
                
                ⁡
                
                  (
                  
                    
                      K
                      f
                    
                    +
                    
                      K
                      r
                    
                  
                  )
                
              
              ⁢
              
                h
                cg
              
            
            
              
                K
                f
              
              ⁢
              
                K
                r
              
              ⁢
              
                L
                2
              
            
          
          ⁢
          
            a
            x
            r
          
        
      
    
    
      
        (
        38
        )
      
    
  




where θ is proportional to Kpaxr, and Ms is a vehicle sprung mass, which is the vehicle mass and payload above the suspension, and hcg is the height of the vehicle CG above ground in the vertical direction. Kf and Kr are predetermined suspension stiffness at the front axle and the rear axle, respectively. axr is an equivalent reactive acceleration and arx is equal to ax+g sin αr. An initial thrust to the pitch dynamic model is based on axr and shown below:

θraw=Kp({dot over (V)}x+axg)  (39)

A first order relative pitch dynamic is provided by equation 40 as shown below:




  
    
      
        
          θ
          +
          
            
              
                D
                pitch
              
              
                K
                pitch
              
            
            ⁢
            
              θ
              .
            
          
        
        =
        
          θ
          raw
        
      
    
    
      
        (
        40
        )
      
    
  




which has a transfer function of:





  
    
      
        
          T
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            K
            pitch
          
          
            
              K
              pitch
            
            +
            
              
                D
                pitch
              
              ⁢
              S
            
          
        
      
    
    
      
        (
        41
        )
      
    
  




where Dpitch represents the general damping effect of the suspension and Kpitch represents the suspension stiffness to the pitch thrust from the longitudinal force.

FIG. 17 depicts a diagram for selecting Dpitch and Kpitch using a phase plane gain scheduling rule, which is generally represented by numeral 1710. The diagram 1710 illustrates eight regions, which are centered about an origin of a Cartesian coordinate system. The x-axis corresponds to {dot over (V)}x+axg and the y-axis corresponds to the pitch angle (θ). The pitch angle θ is bounded by θmax and θmin as shown. Generally, in the normal region, there is a linear relationship, however in stiffer region there is a non linear relationship. The values are slightly different between pitching forward and rearward.
The vehicle system uses the second compensation strategy to estimate the pitch angle (θ) based on Vx, axs, and a vertical acceleration (azs). The inertial sensors 52 of FIG. 1, provide the vertical acceleration (azs) according to one or more embodiments.
For convenience variable (Ψ) may be defined as shown below:

Ψ=azs+lc{dot over (r)}−({dot over (V)}x−r*Vy)  (42)

The following equation can be obtained by combining equations 7, 9, and 42 and substituting (αr+0):




  
    
      
        
          
            a
            z
            s
          
          -
          g
        
        =
        
          
            
              
                V
                .
              
              x
            
            ⁢
            θ
          
          -
          
            
              1
              
                2
                ⁢
                g
              
            
            [
            
              
                Ψ
                2
              
              +
              
                2
                ⁢
                
                  Ψ
                  ⁡
                  
                    (
                    
                      
                        
                          V
                          .
                        
                        x
                      
                      -
                      
                        r
                        *
                        
                          V
                          y
                        
                      
                    
                    )
                  
                
                ⁢
                
                  
                    θ
                    2
                  
                  2
                
              
              +
              
                
                  
                    (
                    
                      
                        
                          V
                          .
                        
                        x
                      
                      -
                      
                        r
                        *
                        
                          V
                          y
                        
                      
                    
                    )
                  
                  2
                
                ⁢
                
                  
                    θ
                    4
                  
                  4
                
              
            
            ]
          
        
      
    
    
      
        (
        43
        )
      
    
  




Equation 43 may be simplified by ignoring the high order term with θ4, as shown below:

Ψ({dot over (V)}x−r*Vy)θ2−2g{dot over (V)}xθ+Ψ2+2g(azs−g)=0  (44)

An equation for calculating the relative pitch angle (θ) according to the second strategy is obtained by solving for the root of equation 44, as shown in equation 45:




  
    
      
        θ
        =
        
          
            g
            ⁢
            
              
            
            ⁢
            
              
                V
                .
              
              x
            
            ⁢
            
              
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      V
                      .
                    
                    x
                    2
                  
                
                -
                
                  
                    Ψ
                    ⁡
                    
                      (
                      
                        
                          
                            V
                            .
                          
                          x
                        
                        -
                        
                          r
                          *
                          
                            V
                            y
                          
                        
                      
                      )
                    
                  
                  ⁡
                  
                    [
                    
                      
                        Ψ
                        2
                      
                      +
                      
                        2
                        ⁢
                        
                          g
                          ⁡
                          
                            (
                            
                              
                                a
                                z
                                s
                              
                              -
                              g
                            
                            )
                          
                        
                      
                    
                    ]
                  
                
              
            
          
          
            Ψ
            ⁡
            
              (
              
                
                  V
                  x
                
                -
                
                  r
                  *
                  
                    V
                    y
                  
                
              
              )
            
          
        
      
    
    
      
        (
        45
        )
      
    
  




The vehicle system uses the third compensation strategy to estimate the pitch angle (θ) based on Vx, axs, and gravity component axg provided by the extended Kalman filter. An approximation equation is provided below:




  
    
      
        
          
            sin
            ⁢
            
              
            
            ⁢
            θ
          
          ≈
          
            
              
                
                  a
                  z
                  s
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                
              
              ⁢
              
                ∝
                r
              
            
            
              a
              x
              s
            
          
        
        =
        
          
            
              a
              z
              s
            
            -
            
              
                
                  g
                  2
                
                -
                
                  
                    (
                    
                      a
                      x
                      g
                    
                    )
                  
                  2
                
              
            
          
          
            a
            x
            s
          
        
      
    
    
      
        (
        46
        )
      
    
  




The pitch angle (θ) may be iteratively identified together with an estimated road gradient (αr) through regression based on the following equation:




  
    
      
        {
        
          
            
              
                
                  
                    
                      a
                      _
                    
                    x
                    off
                  
                  ⁡
                  
                    (
                    
                      k
                      +
                      1
                    
                    )
                  
                
                =
                
                  
                    
                      l
                      c
                    
                    ⁢
                    
                      r
                      .
                    
                  
                  +
                  
                    r
                    *
                    
                      V
                      y
                    
                  
                  -
                  
                    g
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        
                          a
                          z
                          s
                        
                        -
                        
                          
                            
                              g
                              2
                            
                            -
                            
                              
                                (
                                
                                  
                                    
                                      a
                                      ^
                                    
                                    x
                                    g
                                  
                                  ⁡
                                  
                                    (
                                    k
                                    )
                                  
                                
                                )
                              
                              2
                            
                          
                        
                      
                      
                        a
                        x
                        s
                      
                    
                  
                
              
            
          
          
            
              
                
                  
                    
                      a
                      ^
                    
                    x
                    g
                  
                  ⁡
                  
                    (
                    
                      k
                      +
                      1
                    
                    )
                  
                
                =
                
                  
                    
                      
                        a
                        _
                      
                      x
                      off
                    
                    ⁡
                    
                      (
                      
                        k
                        +
                        1
                      
                      )
                    
                  
                  -
                  
                    
                      
                        a
                        _
                      
                      x
                      off
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                  +
                  
                    
                      
                        a
                        ^
                      
                      x
                      g
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        47
        )
      
    
  




where âxg(k=0)=axg(0 from the Kalman filter output based on an initial compensation value of axoff(t). This value may be zero (0), any other admissible value or the regression result from the previous cycle. āxoff(k+=0)=axoff (t).

Equations for estimating the road gradient ({circumflex over (α)}r) and the relative pitch angle ({circumflex over (θ)}) according to the third compensation strategy are provided by equations 48 and 49, as shown below:




  
    
      
        
          
            α
            ^
          
          r
        
        =
        
          arcsin
          (
          
            
              
                
                  a
                  ^
                
                x
                g
              
              ⁡
              
                (
                ∞
                )
              
            
            g
          
          )
        
      
    
    
      
        (
        48
        )
      
    
  
  
    
      
        
          θ
          ^
        
        =
        
          arcsin
          (
          
            
              
                
                  l
                  c
                
                ⁢
                
                  r
                  .
                
              
              +
              
                r
                *
                
                  V
                  y
                
              
              -
              
                
                  
                    a
                    _
                  
                  x
                  off
                
                ⁡
                
                  (
                  ∞
                  )
                
              
            
            g
          
          )
        
      
    
    
      
        (
        49
        )
      
    
  




where âxg(∞) and āxoff(∞) indicate the converged value of âxg(k) and āxoff(k) through the regression process as k>∞. The convergence generally occurs within 5-10 iterations.

As such the vehicle system provides advantages over the existing methods by estimating kinematic road gradient using input compensated acceleration offset to improve accuracy to remove uncertainties and offsets caused by lateral vehicle dynamics and their coupling effects.
The vehicle system estimates a dynamic road gradient (RGEdyn) based on input corresponding to estimated vehicle mass, vehicle speed, wheel torques (brake torque and powertrain output torque), and other drag forces acting on the vehicle (e.g., aerodynamic drag, rolling resistance and the road gradient load). Generally the dynamic RGE has lower quality conditions as compared to static and kinematic RGE because it involves more signals in its estimation, and the noise and error in each signal will be propagated in the final RGE output. However, the dynamic RGE does not rely on the longitudinal acceleration signal (axs), and therefore it has different signal coverage from that of the static and the kinematic estimations. Although dynamic road gradient estimation is secondary, the vehicle mass estimation portion is a useful parameter to many vehicle systems.
The vehicle system estimates vehicle mass independent of the road gradient estimate using a recursive least square (RLS) algorithm under certain vehicle conditions. Generally, the vehicle mass only changes when an object (e.g., a passenger) is added to or removed from the vehicle. Such a change usually only occurs when the vehicle is stopped. The vehicle system estimates vehicle mass while the vehicle is moving using vehicle dynamic equations, and independent of the road gradient estimation. By estimating vehicle mass and road gradient independently, or in parallel, the vehicle system improves the accuracy of each estimate over existing methods. The vehicle system estimates vehicle mass using an event-seeking strategy, whereby the mass is estimated at vehicle conditions where the signal to noise ratio of certain inputs is high. The vehicle system then estimates vehicle mass using a recursive least square (RLS) algorithm when the sought event occurs. The vehicle system re-initializes certain parameters in the mass estimation algorithm at vehicle conditions corresponding to potential mass changes (e.g., when the vehicle stops) in order to increase the sensitivity of the algorithm.
The vehicle system estimates a vehicle mass using dynamic vehicle frame equations that are based on Newton's Second Law (ΣF=ma) as provided below:

Me{dot over (V)}x=Fpwt−Fbrk−Fdrag−Fg  (50)

where Me is the estimated vehicle mass without payload (passengers or cargo), {dot over (V)}x is the derivative of vehicle speed, Fpwt is the propulsion force at the output of the powertrain, Fbrk is the brake force, Fdrag is the drag force and Fg is the force due to gravity.

Equation 50 may be rewritten to solve for Me and further separated into equations based on whether or not the vehicle is in gear, as shown below:




  
    
      
        
          M
          e
        
        =
        
          {
          
            
              
                
                  
                    M
                    v
                  
                  +
                  
                    
                      J
                      w
                    
                    
                      R
                      w
                      2
                    
                  
                  +
                  
                    
                      
                        r
                        p
                        2
                      
                      ⁢
                      
                        η
                        p
                        2
                      
                      ⁢
                      
                        J
                        p
                      
                    
                    
                      R
                      w
                      2
                    
                  
                
              
              
                
                  In
                  ⁢
                  
                    
                  
                  ⁢
                  gear
                
              
            
            
              
                
                  
                    M
                    v
                  
                  +
                  
                    
                      J
                      w
                    
                    
                      R
                      w
                      2
                    
                  
                
              
              
                
                  Not
                  ⁢
                  
                    
                  
                  ⁢
                  in
                  ⁢
                  
                    
                  
                  ⁢
                  gear
                
              
            
          
        
      
    
    
      
        (
        51
        )
      
    
  




where Mv is estimated vehicle mass, which includes the vehicle mass without payload (Me) and the vehicle payloads. The term “in gear” generally references whether or not the powertrain is connected to the wheels. Jw and Jp are the moment of inertia at the wheels and the powertrain, respectively. Rw is the wheel radius. rp is the lumped powertrain torque transfer ratio, and ηp is the lumped torque transfer efficiency. Jp, rp, and ηp vary depending on transmission gear ratio. These values also vary for hybrid electric vehicles having multiple propulsion devices (e.g., a motor and an engine).

An inertia dependent vehicle mass (Mi) is calculated according to equation 52, as shown below:




  
    
      
        
          M
          i
        
        =
        
          {
          
            
              
                
                  
                    
                      J
                      w
                    
                    +
                    
                      
                        r
                        p
                        2
                      
                      ⁢
                      
                        η
                        p
                        2
                      
                      ⁢
                      
                        J
                        p
                      
                    
                  
                  
                    R
                    w
                    2
                  
                
              
              
                
                  In
                  ⁢
                  
                    
                  
                  ⁢
                  gear
                
              
            
            
              
                
                  
                    J
                    w
                  
                  
                    R
                    w
                    2
                  
                
              
              
                
                  Not
                  ⁢
                  
                    
                  
                  ⁢
                  in
                  ⁢
                  
                    
                  
                  ⁢
                  gear
                
              
            
          
        
      
    
    
      
        (
        52
        )
      
    
  




An equivalent longitudinal tire force (Fpwt) acting on the wheels may be calculated based on the total powertrain propulsion torque (Tpwt—whl) or the total regenerative braking torque (Tbrk—whl) according to equation 53 as shown below:




  
    
      
        
          F
          pwt
        
        =
        
          
            
              T
              pwt_whl
            
            
              R
              w
            
          
          =
          
            
              
                1
                
                  R
                  w
                
              
              ⁢
              
                (
                
                  
                    T
                    eng_whl
                  
                  +
                  
                    T
                    mot_whl
                  
                
                )
              
            
            =
            
              
                
                  
                    
                      r
                      pe
                    
                    ⁢
                    
                      η
                      pe
                    
                  
                  
                    R
                    w
                  
                
                ⁢
                
                  T
                  eng
                
              
              +
              
                
                  
                    
                      r
                      pm
                    
                    ⁢
                    
                      η
                      pm
                    
                  
                  
                    R
                    w
                  
                
                ⁢
                
                  T
                  mot
                
              
            
          
        
      
    
    
      
        (
        53
        )
      
    
  




where Teng and Tmot are the engine and electric motor output torque, and Teng—whl and Tmot—whl are their corresponding torques at the wheels. rpe, rpm, ηpe, and ηpm are lumped torque transfer ratio and efficiency parameters of the engine, and the motor respectively.

An equivalent longitudinal brake force (Fbrk) acting on the wheels may be calculated based on the total frictional braking torque according to equation 54 as shown below:




  
    
      
        
          F
          brk
        
        =
        
          
            1
            
              R
              w
            
          
          ⁢
          
            K
            brk
          
          ⁢
          
            P
            brk
          
        
      
    
    
      
        (
        54
        )
      
    
  




where Pbrk is the effective brake line pressure in the wheel brake chambers at the wheel assemblies. This is an equivalent term evaluated on all wheels, because individual wheel brake chambers may have different applied brake line pressure. Kbrk is the frictional brake effectiveness factor that determines how much frictional brake torque can be realized over the braking pads given a certain brake pressure. The value of Kbrk varies with respect to brake system conditions such as temperature, humidity, etc. In one embodiment, Kbrk corresponds to a predetermined constant.

A drag force (Fdrag) may be calculated according to equation 55 as shown below:




  
    
      
        
          F
          drag
        
        =
        
          
            
              1
              2
            
            ⁢
            K
            ⁢
            
              
            
            ⁢
            
              ρ
              a
            
            ⁢
            
              A
              d
            
            ⁢
            
              C
              d
            
            ⁢
            
              V
              x
              2
            
          
          +
          
            
              M
              v
            
            ⁢
            
              gf
              r
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              α
              r
            
          
        
      
    
    
      
        (
        55
        )
      
    
  




where ρa and Cd are constant aerodynamic parameters, Ad is a vehicle cross sectional area, and fr is the rolling resistance coefficient.

A load on the vehicle due to gravity (Fg) may be calculated according to equation 56 as shown below:

Fg=Mvg sin αr  (56)

where αr is the road gradient.

Equation 56 may be converted to the torque domain as shown in equation 57 below:

Tg=Mvg Rw sin αr  (57)

where Rw is the wheel radius.

A first method is provided for estimating vehicle mass concurrently with road gradient using kinematic equations using a least-square type of estimator based on the following parameterized model:




  
    
      
        
          
            
              y
              =
              
              ⁢
              
                
                  V
                  .
                
                x
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      F
                      pwt
                    
                    -
                    
                      F
                      brk
                    
                    -
                    
                      F
                      drag
                    
                  
                  
                    M
                    e
                  
                
                -
                
                  
                    
                      M
                      v
                    
                    
                      M
                      e
                    
                  
                  ⁢
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    α
                    r
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        F
                        pwt
                      
                      -
                      
                        F
                        brk
                      
                      -
                      
                        F
                        drag
                      
                    
                    )
                  
                  ⁢
                  
                    1
                    
                      
                        M
                        v
                      
                      +
                      
                        M
                        i
                      
                    
                  
                
                -
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      
                        M
                        v
                      
                      ⁢
                      
                        
                      
                      ⁢
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      
                        α
                        r
                      
                    
                    
                      
                        m
                        v
                      
                      +
                      
                        M
                        i
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  φ
                  T
                
                ⁢
                γ
              
            
          
        
      
    
    
      
        (
        58
        )
      
    
  




The estimation results for vehicle mass (Mv) and road gradient (αr) based on the model of equation 58 are obtained as follows:




  
    
      
        
          M
          v
        
        =
        
          
            1
            
              γ
              1
            
          
          -
          
            M
            i
          
        
      
    
    
      
        (
        59
        )
      
    
  
  
    
      
        
          α
          r
        
        =
        
          arc
          ⁢
          
            
          
          ⁢
          
            sin
            (
            
              
                γ
                2
              
              
                1
                -
                
                  
                    γ
                    1
                  
                  ⁢
                  
                    M
                    i
                  
                
              
            
            )
          
        
      
    
    
      
        (
        60
        )
      
    
  




However solving for two unknowns (αr and Mv) simultaneously adds additional error to the estimates.
A second approach for estimating the vehicle mass using kinematic equations is to separate the vehicle mass estimate from the road gradient estimate based on the following mathematical model:




  
    
      
        
          
            a
            x
          
          +
          
            g
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              α
              
                r
                ⁢
                
                  
                
              
            
          
        
        =
        
          
            
              1
              
                M
                e
              
            
            ⁢
            
              (
              
                
                  F
                  pwt
                
                -
                
                  F
                  brk
                
                -
                
                  F
                  drag
                
              
              )
            
          
          +
          
            
              
                
                  M
                  e
                
                -
                
                  M
                  v
                
              
              
                M
                e
              
            
            ⁢
            g
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              α
              r
            
          
        
      
    
    
      
        (
        61
        )
      
    
  




Equation may be simplified as follows, as described above with reference to equation 10, and reproduced below:

{dot over (V)}x=ax=axs+axoff−g sin αr  (10)

where axoff represents an acceleration offset and is equal to (lc{dot over (r)}+r*Vy−g sin θ). lc is the lateral offset, r is the yaw rate, Vy is the lateral velocity and θ is the vehicle pitch angle. In one or more embodiments, the vehicle system does not receive these values (lc, Vy, and θ) as inputs. Rather the vehicle system estimates these values from other inputs using a number of compensation strategies.

The following equation is obtained by combining equations 61 and 10:




  
    
      
        
          
            a
            x
            s
          
          +
          
            a
            x
            off
          
        
        =
        
          
            
              1
              
                M
                e
              
            
            ⁢
            
              (
              
                
                  F
                  pwt
                
                -
                
                  F
                  brk
                
                -
                
                  F
                  drag
                
              
              )
            
          
          +
          
            a
            
              δ
              ⁢
              
                
              
              ⁢
              g
            
          
        
      
    
    
      
        (
        62
        )
      
    
  




The variable (aδg) represents the acceleration disturbance excited from road gradient variation, and may be calculated according to equation 63 as shown below:




  
    
      
        
          a
          
            δ
            ⁢
            
              
            
            ⁢
            g
          
        
        =
        
          
            
              
                M
                e
              
              -
              
                M
                v
              
            
            
              M
              e
            
          
          ⁢
          
            
          
          ⁢
          g
          ⁢
          
            
          
          ⁢
          sin
          ⁢
          
            
          
          ⁢
          
            α
            r
          
        
      
    
    
      
        (
        63
        )
      
    
  




where Mv is the sum of the vehicle mass Me and vehicle payloads.

The acceleration disturbance (aδg) contributes to the estimation error. When the vehicle payload is small, the aδg may be ignored. aδg may also be regressively compensated after an estimated vehicle mass is in steady state and converged. The acceleration offset (axoff) may be determined and used as feedforward input compensation as described above with reference to FIGS. 14-17.
Equation 62 may be rearranged to solve for Me and Mv as shown below:




  
    
      
        
          M
          e
        
        =
        
          
            
              
                F
                pwt
              
              -
              
                F
                brk
              
              -
              
                F
                drag
              
            
            
              
                a
                x
                s
              
              +
              
                a
                x
                off
              
              -
              
                a
                
                  δ
                  ⁢
                  
                    
                  
                  ⁢
                  g
                
              
            
          
          =
          
            
              
                
                  
                    T
                    pwt_whl
                  
                  -
                  
                    
                      K
                      brk
                    
                    ⁢
                    
                      P
                      brk
                    
                  
                
                
                  R
                  w
                
              
              -
              
                F
                drag
              
            
            
              
                a
                x
                s
              
              +
              
                a
                x
                off
              
              -
              
                a
                
                  δ
                  ⁢
                  
                    
                  
                  ⁢
                  g
                
              
            
          
        
      
    
    
      
        (
        64
        )
      
    
  
  
    
      
        
          M
          v
        
        =
        
          
            
              
                
                  
                    T
                    pwt_whl
                  
                  -
                  
                    
                      K
                      brk
                    
                    ⁢
                    
                      P
                      brk
                    
                  
                
                
                  R
                  w
                
              
              ⁢
              
                F
                drag
              
            
            
              
                a
                x
                s
              
              +
              
                a
                x
                off
              
              -
              
                a
                
                  δ
                  ⁢
                  
                    
                  
                  ⁢
                  g
                
              
            
          
          -
          
            M
            i
          
        
      
    
    
      
        (
        65
        )
      
    
  




The drag force (Fdrag) and the road gradient (αr) are compensated in order to obtain vehicle mass estimates. However, it is difficult to accurately compensate these terms due to signal noise. These terms exhibit slow dynamic variations as compared to the variations exhibited by the vehicle longitudinal dynamic. Therefore these terms may be treated as unknown constants in a short period estimation event. A difference based estimation equation may be derived as follows:

axi(t)Me=Fwhl(t)−Fdrag(t)+Meaδg(t)  (66)

axi(t+Δt)Me=Fwhl(t+Δt)−Fdrag(t+Δt)+Me aδg(t+Δt)  (67)

where the longitudinal acceleration (axi) at time index (i) is equal to (axs+axoff) and the force at the wheels (Fwhl) is equal to (Tpwt—whl−KbrkPbrk)/Rw. The variable (t) is a time instance within a qualified estimation event. At is a predefined time interval parameter.

As Fdrag (t)≈Fdrag(t+Δt) and aδg (t)≈aδg (t+Δt), an incremental vehicle longitudinal dynamic model may be obtained by taking the difference between equations 66 and 67 as shown below:




  
    
      
        
          M
          v
        
        =
        
          
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  F
                  whl
                
              
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  a
                  x
                  i
                
              
            
            ⁢
            
              ❘
              
                Δ
                ⁢
                
                  
                
                ⁢
                t
              
            
            ⁢
            
              -
              
                M
                i
              
            
          
          =
          
            
              
                
                  
                    F
                    whl
                  
                  ⁡
                  
                    (
                    
                      t
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                    
                    )
                  
                
                -
                
                  
                    F
                    whl
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
              
                
                  
                    a
                    x
                    i
                  
                  ⁡
                  
                    (
                    
                      t
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                    
                    )
                  
                
                -
                
                  
                    a
                    x
                    i
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
            
            -
            
              M
              i
            
          
        
      
    
    
      
        (
        68
        )
      
    
  




where ΔFwhl and Δaxi are the difference between the values of Fwhl and axi evaluated in the time interval Δt within a qualified estimation event.

Estimating vehicle mass (Mv) using such an approach isolates the vehicle mass estimate from noise factors such as: aerodynamic noise, rolling resistance and road gradient. As Δt→0, the equation 68 may be simplified to a derivative format, according to equation 69 as shown below:




  
    
      
        
          M
          v
        
        =
        
          
            
              
                
                  F
                  .
                
                whl
              
              
                
                  a
                  .
                
                x
                i
              
            
            -
            
              M
              i
            
          
          =
          
            
              
                1
                
                  R
                  w
                
              
              ⁢
              
                
                  
                    (
                    
                      
                        T
                        pwt_whl
                      
                      -
                      
                        
                          K
                          brk
                        
                        ⁢
                        
                          P
                          brk
                        
                      
                    
                    )
                  
                  ⁢
                  ′
                
                
                  
                    (
                    
                      
                        a
                        x
                        s
                      
                      +
                      
                        a
                        x
                        off
                      
                    
                    )
                  
                  ⁢
                  ′
                
              
            
            -
            
              M
              i
            
          
        
      
    
    
      
        (
        69
        )
      
    
  




With reference to FIG. 18, a schematic block diagram illustrating a method for vehicle mass estimation is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1810. The method 1110 is included within the dynamic road gradient estimation (RGEdyn) block 426 of FIG. 4, according to one or more embodiments. At block 1812 the vehicle system receives inputs including acceleration (a) yaw rate (r), vehicle speed (Vx), brake pressure (Pbrk) powertrain torque (Tpwt) and braking torque (Tbrk). The input may be received directly from sensors described with respect to FIG. 1, or indirectly over the CAN bus.
At block 1814, the vehicle system determines dynamic mass input parameters Yk and φk which represent a system output, and a system input respectively. The input parameters are based on the calculated vehicle speed (Vx) the longitudinal acceleration (axs), a powertrain force (Fpwt), a braking force (Fbrk), and an acceleration offset (axoff). Fpwt and Fbrk are calculated based on Tpwt and Tbrk. At block 1816, the vehicle system estimates the dynamic mass (My). At block 1818, the vehicle system determines the acceleration offset axoff. The acceleration offset is provided as feed-forward RGE input compensation. At block 1820 the vehicle system estimates road gradient αr.
The vehicle system estimates the dynamic vehicle mass (My) at block 1816 using a recursive least square (RLS) estimation according to equation 65, 68 or 69, according to one or more embodiments. The RLS estimation is based on the following parameterized estimation model:

y=φ*θ  (70)

A recursive least square estimation (RLS) estimation algorithm is developed as shown below:




  
    
      
        
          L
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            
              P
              ⁡
              
                (
                
                  k
                  -
                  1
                
                )
              
            
            ⁢
            
              φ
              ⁡
              
                (
                k
                )
              
            
          
          
            1
            +
            
              
                
                  φ
                  2
                
                ⁡
                
                  (
                  k
                  )
                
              
              ⁢
              
                P
                ⁡
                
                  (
                  
                    k
                    -
                    1
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        71
        )
      
    
  
  
    
      
        
          
            θ
            ^
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            
              θ
              ^
            
            ⁡
            
              (
              
                k
                -
                1
              
              )
            
          
          +
          
            
              L
              ⁡
              
                (
                k
                )
              
            
            ⁡
            
              [
              
                
                  y
                  ⁡
                  
                    (
                    k
                    )
                  
                
                -
                
                  
                    φ
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                  *
                  
                    
                      θ
                      ^
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        72
        )
      
    
  
  
    
      
        
          P
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            [
            
              1
              -
              
                
                  L
                  ⁡
                  
                    (
                    k
                    )
                  
                
                ⁢
                
                  φ
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            ]
          
          ⁢
          
            P
            ⁡
            
              (
              
                k
                -
                1
              
              )
            
          
        
      
    
    
      
        (
        73
        )
      
    
  




where {circumflex over (θ)} (k) is the estimated version of parameter θ at k-th execution cycle. The RLS estimator is a filtering strategy that averages the estimated data updated based on the error between the measurement and the predicted model output. Such an estimation strategy may be used when the parameter to be estimated is constant in nature. To initiate the estimation, P(0) is set equal to 1.

The parameters of the RLS equations 71-73 vary depending on which vehicle mass estimation equation is used, (e.g., equation 65, 68, or 69).
The following parameters are used when the RLS equations (71-73) are calculated based on the vehicle mass estimation of equation 65:




  
    
      
        
          
            
              θ
              =
              
                M
                e
              
            
          
        
        
          
            
              φ
              =
              
                
                  a
                  x
                  s
                
                +
                
                  a
                  x
                  off
                
              
            
          
        
        
          
            
              y
              =
              
                
                  
                    T
                    pwt_whl
                  
                  -
                  
                    
                      K
                      brk
                    
                    ⁢
                    
                      P
                      brk
                    
                  
                
                
                  
                    R
                    w
                  
                  -
                  
                    F
                    drag
                  
                
              
            
          
        
      
    
    
      
        (
        74
        )
      
    
  




The following parameters are used when the RLS equations (71-73) are calculated based the difference based vehicle mass estimation of equation 68:

θ=Me 

φ=Δaxi|Δt 

y=ΔFwhl|Δt  (75)

Alternatively, the following parameters are used when the RLS equations (71-73) are calculated based the differential based algorithm of equation 69:




  
    
      
        
          φ
          =
          
            
              
                a
                .
              
              x
              i
            
            =
            
              
                (
                
                  
                    a
                    x
                    s
                  
                  +
                  
                    a
                    x
                    off
                  
                
                )
              
              ⁢
              ′
            
          
        
        ⁢
        
          

        
        ⁢
        
          y
          =
          
            
              
                F
                .
              
              whl
            
            =
            
              
                
                  (
                  
                    
                      T
                      
                        pwt
                        whl
                      
                    
                    -
                    
                      
                        K
                        brk
                      
                      ⁢
                      
                        P
                        brk
                      
                    
                  
                  )
                
                ⁢
                ′
              
              
                R
                w
              
            
          
        
      
    
    
      
        (
        76
        )
      
    
  




The vehicle system uses the RLS to estimate vehicle mass at certain qualified vehicle events. The mass estimator is developed based on the condition when vehicle dynamic is predominantly longitudinal. Thus the vehicle system searches for conditions under which the vehicle longitudinal inertia effect is significantly dominating the related noise factors, or the signal to noise ratio is high. These qualified events generally occur during vehicle braking/deceleration or during vehicle acceleration.
In one or more embodiments, the vehicle system detects a qualified vehicle braking/deceleration event when each of the following conditions are satisfied:

V≧C—DynMassBrkEntry;  1.

axi≦C—aRGEAxNeg;  2.

|r|≦C_small—yaw;  3.

and

(Fwhl<C—FBrk and {dot over (F)}whl<C—DFBrk and {umlaut over (F)}whl<C—D2FBrk)or(Fwhl<C—FBrk2 and {dot over (F)}whl<C—DFBrk2).  4.

In one or more embodiments, the vehicle system detects a qualified vehicle acceleration event when each of the following conditions are satisfied:

V≧C—DynMassPropEntry;  1.

4≧C—aRGEAxPos;  2.

|r|≦C_small—yaw;  3.

and

(Fwhl>C—FProp and {dot over (F)}hi<C—DFPropUp and {dot over (F)}whl>C—DFPropDn and {umlaut over (F)}hi<C—D2FPropUp and {umlaut over (F)}whl<C—D2FPropDn)or(Fwhl>C—FProp2 and Pwhl>C—DFProp2 and {umlaut over (F)}whl>C—D2FProp2)  4.

In the above conditions, a variable that begins with “C” indicates a calibration parameter. The calibration parameters are predetermined values designed to target a specific estimation event. The conditions provide a tunable construction of qualified predominantly longitudinal vehicle motion. The more strict the parameters are tuned towards dynamic significance, the higher the estimation accuracy that may be achieved. However, such strict parameters will reduce the available estimation time interval, since a qualified event will be less frequent. Thus a qualified estimation event depends on a certain level of driving aggressiveness. The more forceful or rapid the vehicle acceleration or braking is, the more dominant the inertia effect will be excited in the vehicle dynamic. The tuning of the qualified estimation parameters is a tradeoff between estimation accuracy and event occurrence, since aggressive driving behavior is less common during a normal driving process.
With reference to FIG. 19, a method for estimating the vehicle mass is illustrated in accordance with one or more embodiments and is generally referenced by numeral 1910. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 1910 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The method 1910 is applicable to both acceleration and braking/deceleration events and corresponds to the event conditions listed above.
In operation 1912 the vehicle system evaluates the vehicle speed (Vx) to determine if Vx is greater than an estimation entry speed threshold. Operation 1912 corresponds to the first event condition. The estimated entry speed threshold corresponds to a calibrated brake entry speed threshold (C_DynMassBrkEntry) for the braking event; and to a calibrated propulsion entry speed threshold (C_DynMassPropEntry) for the acceleration event. If the determination at operation 1912 is positive, the vehicle system proceeds to operation 1914.
At operation 1914 the vehicle system updates the estimation model input (φ) and output (Y) based on the vehicle mass estimation algorithm selected, as described with reference to equations 74-76. At operation 1916 the vehicle system computes the estimation condition signals ({dot over (F)}whl and {umlaut over (F)}whl), as described with reference to equations 75 and 76.
At operation 1918 the vehicle system evaluates the longitudinal acceleration (axi) and the yaw rate (r) to determine if the vehicle dynamic is primarily longitudinal with a sufficiently small yaw rate. Operation 1918 corresponds to the second and third event conditions. For a braking event, the vehicle dynamic is primarily longitudinal if the longitudinal acceleration (axi) is less than or equal to a calibrated negative acceleration threshold (C_aRGEAxNeg). For an acceleration event, the vehicle dynamic is primarily longitudinal if the longitudinal acceleration (axi) is greater than or equal to a calibrated positive acceleration threshold (C_aRGEAxPos). For both a braking event and an acceleration event, the yaw rate (r) is sufficiently small if the magnitude of the yaw rate (|r|) is less than or equal to a calibrated small yaw rate threshold (C_small_yaw). If the determination at operation 1918 is positive, the vehicle system proceeds to operation 1920.
At operation 1920 the vehicle system evaluates the force at the wheels (Fwhl) along with higher order derivatives of the force at the wheels ({dot over (F)}whl and {umlaut over (F)}whl) to determine if a qualified estimation event is present. Operation 1920 corresponds to the fourth event condition.
For a braking event, the fourth condition is satisfied if: the wheel force (Fwhl) is less than a calibrated brake force (C_Fbrk); and the derivative of the wheel force ({dot over (F)}whl) is less than a derivative of the calibrated brake force (C_DFBrk); and a second derivative of the wheel force ({umlaut over (F)}whl) is less than a second derivative of the calibrated brake force (C_D2FBrk). Alternatively, the fourth braking condition may also be satisfied if the wheel force (Fwhl) is less than a second calibrated brake force (C_Fbrk2); and the derivative of the wheel force ({dot over (F)}whl) is less than a derivative of the second calibrated brake force (C_DFBrk2).
For an acceleration event, the fourth condition is satisfied if: the wheel force (Fwhl) is greater than a calibrated propulsion force (C_FProp); and the derivative of the wheel force ({dot over (F)}whl) is less than a derivative of the calibrated uphill propulsion force (C_DFPropUp); and the derivative of the wheel force ({dot over (F)}whl) is greater than a derivative of the calibrated downhill propulsion force (C_DFPropDn); and a second derivative of the wheel force ({umlaut over (F)}whl) is less than a second derivative of the calibrated uphill propulsion force (C_D2FPropUp); and the second derivative of the wheel force ({umlaut over (F)}whl) is less than a second derivative of the calibrated downhill propulsion force (C_D2FPropDn). Alternatively, the fourth acceleration condition may also be satisfied if the wheel force (Fwhl) is greater than a second calibrated propulsion force (C_FProp2); and the derivative of the wheel force ({umlaut over (F)}whl) is greater than a derivative of the second calibrated propulsion force (C_DFProp2); and the second derivative of the wheel force ({umlaut over (F)}whl) is greater than a second derivative of the second calibrated propulsion force (C_D2FProp2). If the determination at operation 1920 is positive, the vehicle system proceeds to operation 1922.
At operation 1922 the vehicle system executes an RLS vehicle mass estimation algorithm (e.g., using equations 71-73). At operation 1924 the vehicle system estimates vehicle mass (e.g., using equation 65, 68 or 69). After operation 1924, or in response to negative determinations at operation 1912, 1918 or 1920, the vehicle system proceeds to operation 1926, and then returns to operation 1912.
Vehicle mass is unlikely to change when the vehicle is moving. The vehicle system continues to monitor the vehicle state and carries out estimation updates in the presence of a qualified estimation event. The estimated vehicle mass is averaged throughout all the estimation instances with the assumption that the vehicle mass is constant and the estimation noise is zero mean Gaussian. However, vehicle mass may change when the vehicle is moving slowly and when it is stopped for a certain period of time. Under such vehicle conditions, the vehicle system adjusts the sensitivity of the RLS estimator to a more sensitive state, such that a change in vehicle mass is readily detected.
With reference to FIG. 20, a method for resetting vehicle mass estimation parameters is illustrated in accordance with one or more embodiments and is generally referenced by numeral 2010. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 2010 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The method 2010 resets vehicle mass estimation parameters under certain vehicle conditions for adjusting the sensitivity of the RLS estimator.
At operation 2012 the vehicle system evaluates the vehicle speed (Vx) to determine if the vehicle speed is less than a vehicle stop speed threshold. If the determination at operation 2012 is positive, the vehicle system proceeds to operation 2014 and increments, or counts up the vehicle stop/standstill timer (Tmr_stand) by one cycle time and the stop/standstill timer will start recording the vehicle stopped duration. At operation 2016 the vehicle system resets a drive-off timer (Tmr_drive), by setting Tmr_drive equal to zero.
At operation 2018 the vehicle system evaluates a total count of qualified estimation events (Cntr_EstTot), to determine if Cntr_EstTot is greater than a calibrated estimated event threshold (C_QuaEstEvnt). If the determination at operation 2018 is positive, the vehicle system proceeds to operation 2020 and updates a RLS estimation memory (Vmass_memo) by setting Vmass_memo equal to a present vehicle mass estimate (Vmass_est). Vmass_est corresponds to the vehicle mass (Vmass_rls) estimated at operation 1924 of FIG. 19, and the estimated vehicle mass (Mv) of equations 65, 68 and 69, according to one or more embodiments. After operation 2020, or in response to a negative determination at operation 2018, the vehicle system proceeds to operation 2022.
At operation 2022 the vehicle system evaluates the stop/standstill timer (Tmr_stand) to determine if Tmr_stand is greater than a calibrated standstill timer (C_StandstilTm). If the determination at operation 2022 is positive, the vehicle system proceeds to operation 2024 and resets the RLS parameter and mass estimation boundaries to their original values. For example, the vehicle system may set the RLS estimation parameter P(k) (from equation 73) to its initial value P(0). At operation 2026 the vehicle system resets the total qualified estimation instance counter (Enter_EstTot) by setting Enter_EstTot equal to zero. If the determination at operation 2012 is negative, the vehicle system proceeds to operation 2028.
At operation 2028, the vehicle system evaluates the vehicle speed (Vx) to determine if Vx is greater than a drive-off speed threshold. If the determination at operation 2028 is positive, the vehicle system proceeds to operation 2030 and increments, or counts up the vehicle drive-off Timer (Tmr_drive) by one cycle time, and the drive-off timer will start recording the vehicle moving duration. At operation 2032 the vehicle system evaluates the drive-off timer (Tmr_drive) and the stop/standstill timer (Tmr_stand) to determine if both Tmr_drive is greater than a calibrated drive-off timer (C_StandstilTm); and if Tmr_stand is greater than zero. If the determination at operation 2032 is positive, the vehicle system proceeds to operation 2034. At operation 2034 the vehicle system resets a stop/standstill timer (Tmr_stand), by setting Tmr_stand equal to zero. After operations 2026, and 2034, or in response to a negative determination at operations 2022, 2028 or 2032, the vehicle system proceeds to operation 2336, and then returns to operation 2012.
A brief or transient qualified estimation event, is insufficient for vehicle mass estimation. Generally, a transient event does not satisfy the second condition related to longitudinal inertia. However, even if such a transient event does satisfy the second condition, the vehicle system still removes a transient event from the estimation results because it will take a certain period of estimation presence before the estimated vehicle mass starts to converge. The estimated vehicle mass is updated when a sufficient number of estimation instances have be observed.
With reference to FIG. 21, a method for updating vehicle mass estimation parameters is illustrated in accordance with one or more embodiments and is generally referenced by numeral 2110. The method is implemented using software code contained within the VSC, according to one or more embodiments. In other embodiments, the method 2110 is implemented in other vehicle controllers, or distributed amongst multiple vehicle controllers. The method 2110 updates vehicle mass estimation state when the qualified estimation event exceeds a predetermined time interval.
At operation 2112 the vehicle system determines if a qualified estimation event is present. Operation 2112 corresponds to operations 1912-1920 of FIG. 19. If the determination at operation 2112 is positive the vehicle system proceeds to operation 2114, and computes Vmass_rls as described in operation 1924. At operation 2116 the vehicle system increments an estimation instance counter (Cntr_est) by one cycle time. At operation 2118 the vehicle system evaluates the estimation instance counter (Cntr_est) to determine if Cntr_est is greater than or equal to a calibrated time interval threshold (C_EstEvnt4Update), which corresponds to a sufficient instance of estimation for an estimation event. If the determination at operation 2118 is positive, the vehicle system proceeds to operation 2120 and updates an estimation update state (Vmass_update) by setting Vmass_update equal to the present vehicle mass estimation state (Vmass_rls).
At operation 2122 the vehicle system evaluates the estimation instance counter (Cntr_est) to determine if Cntr_est is equal to the calibrated time interval threshold (C EstEvnt4Update). If the determination at operation 2122 is negative, the vehicle system proceeds to operation 2124 and increments the estimation instance counter (Cntr_est) by one cycle time. If the determination at operation 2122 is positive, the vehicle system increments the accumulative counter (Cntr_EstTot) by Cntr_est, which counts up at each instance in qualified estimation events.
At operation 2128 the vehicle system evaluates the accumulated count of qualified estimation events (Cntr_EstTot) to determine if Cntr_EstTot is greater than a calibrated qualified estimation event threshold (C_QuaEstEvnt). If the determination at operation 2128 is positive, the vehicle system proceeds to operation 2130 and sets the dynamic mass estimation output (Vmass_est) equal to the RLS estimator state (Vmass_rls). If the determination at operation 2128 is negative, the vehicle system proceeds to operation 2132 and sets the dynamic mass estimation output (Vmass_est) equal to the stored estimation memory value (Vmass_memo).
If the determination at operation 2112 is negative, a qualified estimation event is not present, then the vehicle system proceeds to operation 2134. At operation 2134 the vehicle system decreases the estimation instance counter (Cntr_est) by one cycle time. At operation 2136 the vehicle system evaluates the estimation instance counter (Cntr_est) to determine if Cntr_est is less than or equal to zero. If the determination at operation 2136 is positive, the vehicle system proceeds to operation 2138 and sets the present vehicle mass estimation state (Vmass_rls) equal to the estimation update state (Vmass_update). After operations 2130, 2132, 2138, and in response to a negative determination at operation 2136, the vehicle system proceeds to operation 2140, and then returns to operation 2112.
With reference to FIG. 22, the vehicle system limits the vehicle mass estimate using dynamic boundaries according to one or more embodiments. The presence of a single inaccurate vehicle mass estimation instance may distort the estimation state. Therefore the vehicle system establishes boundaries to avoid biased and unreasonable estimation results caused by certain noise factors in the estimation process. As shown in the illustrated embodiment, the vehicle system sets an upper boundary (Ubd) and a lower boundary (Lbd). According to one or more embodiments, the upper boundary (Ubd) corresponds to a maximum gross vehicle weight (GVW), which includes the vehicle, passengers and cargo; and the lower boundary (Lbd) corresponds to a value that is less than the vehicle “curb” weight, which only includes the vehicle. As shown in FIG. 22, the Ubd corresponds to a GVW of approximately 2,500 Kg, and the Lbd corresponds to a Curb weight of approximately 1,700 Kg at initial conditions, with a difference of approximately 800 Kg at (time=0s).
The boundaries are dynamic in that they converge towards the estimated vehicle mass. Since vehicle mass generally does not change while the vehicle is moving, the boundaries narrow or converge toward a moving average of the estimated vehicle mass such that the difference between the two boundaries decreases as more estimation information becomes available to filter out unreasonable vehicle mass estimates. For example, as shown in the illustrated embodiment, at time 150 s, the Ubd has decreased to approximately 2,200 Kg, and the Lbd has increased to approximately 1,750 Kg, with a difference of approximately 450 Kg.
The estimation boundaries will be reset to the original upper and lower limits when the reset condition is satisfied as described above with reference to FIG. 20. Additionally, the boundaries will converge while the vehicle is accumulating qualified estimation updates. As described above with reference to FIG. 21, the vehicle system includes a counter (Cntr_est) that increments every time a qualified acceleration or deceleration event occurs. A mass estimation event flag (Mevent) is illustrated on FIG. 22, and scaled with the other signals such that Mevent is equal to “NQ” when no qualified event is detected, Mevent is equal to “QA” when a qualified acceleration event is detected (e.g., at point 2212); and Mevent is equal to “QB” when a qualified braking event is detected (e.g., at point 2214). Once a qualified event is detected, the counter starts counting up, once per execution cycle. Once the total number of counts exceeds a predetermined value the vehicle system sets the initial estimated vehicle mass equal to the RLS estimated mass as referenced by numeral 2216.
With reference to FIG. 23, as the vehicle system accumulates more vehicle mass information, the estimated mass converges upon an actual vehicle mass. FIG. 23 depicts test results taken for a vehicle having a known mass, which is represented by the actual mass waveform (Mass_actual). FIG. 23 also includes the vehicle speed (Vx), and the qualified event flag (Mevent) waveform. Mevent is equal to “NQ” when no qualified event is detected, Mevent is equal to “QA” when a qualified acceleration event is detected; and Mevent is equal to “QB” when a qualified braking event is detected. The qualified braking events correspond to vehicle deceleration as referenced by numeral 2312. The initial estimated vehicle mass (Vmass_est) is referenced by numeral 2314. The estimated vehicle mass (Vmass_est) converges towards the actual mass (Mass_actual) as additional information is accumulated, as referenced by numeral 2316.
FIG. 24 is another graph illustrating vehicle mass estimation. Although the vehicle mass estimate (Vmass_est) converges toward the actual mass (Mass_actual) over time as additional information is detected, the vehicle mass estimate may adjust away from the actual mass when there is noise present on the input signals, as referenced by numeral 2412. For example, in one or more embodiments, vehicle mass is estimated based on a drag force (Fdrag) as described above with reference to equation 65. The accuracy of the drag force estimate decreases under certain vehicle conditions (e.g., high wind). The boundaries described with reference to FIG. 22 reduce the impact of such noisy conditions, by limiting the estimated vehicle mass.
The vehicle system may estimate a dynamic road gradient based on the estimated vehicle mass, for example by using equation 60. Additionally the estimated vehicle mass may be used for other vehicle systems, such as brake apply and release detection (BARD) for micro-HEVs having start/stop functionality.
As such, the vehicle system provides advantages by estimating the vehicle mass independent of the road gradient estimation using an event-seeking strategy. The mass is estimated at vehicle conditions, or qualified events, where the signal to noise ratio of certain inputs is high. The vehicle system then estimates vehicle mass using a recursive least square (RLS) algorithm when the sought event occurs. The vehicle system re-initializes certain parameters in the mass estimation algorithm at vehicle conditions corresponding to potential mass changes (e.g., when the vehicle stops) in order to increase the sensitivity of the algorithm. The vehicle system estimates a dynamic road gradient (RGEdyn) based on the estimated vehicle mass. Generally the dynamic RGE has lower quality conditions as compared to static and kinematic RGE because it involves more signals in its estimation, and the noise and error in each signal will be propagated in the final RGE output. However, the dynamic RGE does not rely on the longitudinal acceleration signal and therefore it has different signal coverage from that of the static and the kinematic estimations. Additionally, the vehicle system estimates the vehicle mass independent of the road gradient estimation, thereby isolating noise in RGEdyn from the vehicle mass estimate.
As described above, with reference to FIG. 4, the vehicle system determines a quality factor for each of the static, kinematic and dynamic RGE algorithms, which are referenced by the variables: QFst, QFkin, QFdyn respectively. The vehicle system also determines an overall RGE quality factor (QFRGE). Since each RGE algorithm relies on different input signals, they each have different ranges and availability.
The estimation quality factor for the RGE output as well as that for each RGE algorithm (static, kinematic, and dynamic) corresponds to a value between 0-3, and is represented by variables (QFRGE, QFst, QFkin, and QFdyn). The quality factor evaluates both the availability and the accuracy of the estimated information. The quality factor (QF) factor is selected from a range of 0-3, where a quality factor of three (3) indicates full quality where the estimation condition is satisfied. A quality factor of two (2) indicates degraded quality when the estimation condition is not fully satisfied but the estimation state does not deviate too far from its setpoint. A quality factor of one (1) indicates that no qualified estimation has been achieved near the instantaneous spot and that the estimation state has not achieved sufficient convergence, a quality factor of zero indicates estimation not available.
The vehicle system evaluates the RGE input signal quality and processing, which are referenced by numerals 414 and 412 in FIG. 4. The vehicle longitudinal acceleration is an input signal 414 and is measured from an inertia sensor 52 (shown in FIG. 1). The signal is regarded as a good signal if its quality factor (QF) is equal to three (3) and the signal is within a specified operating range, e.g. +/−2 g. The input signal is locked to its latest high quality value when its quality degrades and it will be reset to a default value when the degraded quality persists for longer than a time threshold. Other input signals like, brake torque, powertrain torque, vehicle speed, yaw rate, etc. are all processed in the same manner.
The static RGE algorithm estimates road gradient based on the longitudinal acceleration input (axs). Static RGE is suited for vehicle standstill conditions, such as engine start/stop functionality when a micro-hybrid is stopped on an incline. The vehicle system evaluates the static RGE quality, based on the acceleration sensor measurement signal's stability. The sensor stability is indicated by the signals oscillation and variation level as compared to thresholds on its filtered first order derivative. The highest quality is achieved when the vehicle is in standstill and the vehicle body is stabilized such that the accelerometer output is near a converged steady state. A higher oscillation level will degrade the estimation quality and become unavailable when the preceding kinematic estimation fails to set the static estimation's initial value in the current vehicle stop event.
The kinematic RGE algorithm estimates RGE based on input corresponding to the vehicle's kinematic properties including the relationship between vehicle speed, accelerations and yaw rate. The algorithm processes the input using a Kalman Filter. The kinematic RGE algorithm is suited for normal vehicle motion conditions (e.g., Vx>5 kph) because some of the underlying calculations have low accuracy at low vehicle speeds. Additionally, a Kalman filter is a dynamic filter, a good estimation is usually not available until the estimation state has converged. The transient time depends on the difference between the initial state and the true road grade state as well as the estimation conditions like: signal frequency component richness and signal to noise ratio. The vehicle system evaluates the quality of the kinematic RGE based on assessing the quality of its input signals, its estimation state, vehicle state as well as the estimation settle down time period.
The dynamic RGE algorithm 426 estimates RGE based on input corresponding to vehicle speed, wheel torques (brake torque and powertrain output torque), and other drag forces acting on the vehicle (e.g., aerodynamic drag, rolling resistance and the road gradient load). Generally the dynamic RGE has lower quality conditions as compared to static and kinematic RGE because it involves more signals in its estimation, and the noise and error in each signal will be propagated in the final RGE output. Similar to the kinematic algorithm, the quality the dynamic RGE estimation is also evaluated over its input signal, estimation state, vehicle states, estimation qualification conditions, and estimation convergence/settle-down time period.
The signals that are provided to the Kalman filter based estimator (e.g., RGEkin) are grouped into two sets: the primary signals and the supportive signals. The primary signals are used to compute the inputs and outputs of the estimator system equation. The supportive signals are used to adjust the estimator mode or parameters to be adaptive to different driving situations.
When any of the primary signals are not ok, (e.g., QF<3), the estimation result will start deviating from the true value. The vehicle system evaluates the impact of the degraded signal to the estimation state by counting the time period during which degraded primary signals are present, which is referenced by variable Tdgd. The vehicle system also interpolates a degradation impact factor Cdg, with respect to Tdgd. Such an interpolation may be based on predetermined data stored in a look up table. In general, the principle of determining cdgi is that the longer the presence of the degraded primary signals, the worse the impact on the estimation results. An estimation quality index variable (qestc) is calculated according to equation 77 as shown below:

qestest=Tdgd*Cdgi  (77)

Alternatively, a pre-calibrated value of qest may be used and the estimation process will be halted in the presence of degraded primary signals to the estimator.
After the primary signals all resume good quality, the estimation does not converge back to the true estimation target state immediately. The estimation results can be used confidently after a converging phase occurs. In order to represent the converging phase, a fading factor, fdgi, is designed that takes a value between 0 and 1 to decrease the value of qest that has been accumulated from a previous degraded signal event such that:

qest0=qest0*fdgi  (78)

This computation is repeated in each control implementation cycle to gradually decrease the value of qest down to zero.
At each moment, the occurrence of a present impact increasing can occur simultaneously with the decreasing of the impact from a previous event. Thus, the final estimation quality index is determined by the maximum value among qest0 and qestc, i.e. qest=max(qest0,qestc).
With reference to FIG. 25, a method for estimating the quality of a Kalman estimator, or an RLS dynamic mass estimator, is illustrated in accordance with one or more embodiments and is generally referenced by numeral 2510.
In operation 2512 the vehicle system evaluates a road gradient estimation (RGE) signal to determine if the estimation is locked. An RGE is locked when it is not updating. The RGE may be locked to its latest high quality value when its quality degrades. If the determination at operation 2512 is positive the vehicle system proceeds to operation 2514 and locks the quality factor of the evaluated estimation (QF_est) to its current value. The variable QF_est is a generic value corresponding to the quality factor for kinematic or dynamic road gradient estimations. When the estimation algorithm does not satisfy a qualified estimation condition and the estimation is locked, its associated quality evaluation is locked as well until a next estimation event will happen. If the determination at operation 2512 is negative, the vehicle system proceeds to operation 2516.
At operation 2516 the vehicle system evaluates the road gradient estimation (RGE) signal to determine if the estimation is active. RGE is active, if it is updating. If the determination at operation 2516 is negative, the vehicle system proceeds to operation 2518 and sets the quality factor to a value corresponding to “no data”, (e.g., QF=0). If the determination at operation 2516 is positive, the vehicle system proceeds to operation 2520.
At operation 2520, the vehicle system compares the evaluated quality index (qest), as described above with reference to equations 77 and 78, to an OK threshold parameter (Q_est_ok), to determine if q_est is less than Q_est_ok; and there is no supportive signal degradation. If the determination at operation 2520 is positive, the vehicle system proceeds to operation 2522 and sets the quality factor to a value corresponding to “OK”, (e.g., QF=3). If the determination at operation 2520 is negative the vehicle system proceeds to operation 2524.
At operation 2524, the vehicle system compares the evaluated quality index (qt) to a FAULT threshold parameter (Q_est_degrade), to determine if q_est is greater than or equal to Q_est_degrade. If the determination at operation 2524 is positive, the vehicle system proceeds to operation 2524 and sets the quality factor to a value corresponding to “FAULT”. If the determination at operation 2524 is negative the vehicle system proceeds to operation 2528 and sets the quality factor to a value corresponding to “DEGRADED”, (e.g., QF=2). After operations 2514, 2518, 2522, 2526 and 2528, the vehicle system proceeds to operation 2530, and then returns to operation 2512.
With reference to FIGS. 26A-26C, the vehicle system monitors variation in the kinematic road gradient estimation (RGEkin). FIGS. 26A-26C include 3 graphs of measurement data taken over a common period of time and are generally referenced by numeral 2610. FIG. 26A illustrates vehicle speed (Vx), FIG. 26B illustrates an actual road gradient (RGEref) along with an estimated road gradient (RGEEST). The RGEEST corresponds to an arbitrated output (RGEout) of RGEst, RGEkin, and RGEdyn according to one or more embodiments. FIG. 26C illustrates a calculated percent error between RGEEST and RGEref.
The Kalman filter based estimation, such as that used for RGEkin, assumes that the road gradient state is stable and is close to a constant state in a short period of time or distance. This assumption is generally true, because most roads do not have rapid changes in road gradient. However, as shown in FIGS. 26B and 26C, large estimation error may occur at transitions between different road gradients. For example, the road gradient changes from approximately 0% to approximately 7%, as referenced by numeral 2612. A generally large estimation error is present at this transition, as referenced by numeral 2614.
The vehicle system includes a method, (not shown) for monitoring an estimation variation to determine the estimation quality (accuracy) based on the magnitude of the road gradient variation. In this method, a derivative of the estimated road gradient is calculated (RĠEEST). A transition period is identified, when the derivative (RĠEEST) is greater than a threshold estimation and such condition is true for longer than a predetermined period of time, the transition period is identified. Subsequently, the estimation quality (QF) will be reduced to a degraded level (QF=2) if no other degradation factor is also present simultaneously. Such degradation will be cancelled when the derivative signal settles down to a low value for another settling time period. In general, this method informs other vehicle systems that are using the RGE output, that the estimation is not fully accurate in the transition period of time when road gradient changes significantly in a short period of time.
The road gradient estimations are not available at very low vehicle speeds, as shown in column 6, of Table A. Although none of the estimates are available at very low vehicle speeds, (e.g., less than 5 mph), the vehicle system still provides an accurate road gradient estimate during such conditions. At low vehicle speeds the vehicle speed provided by or derived from a sensor signal, (e.g., wheel speed sensors or GPS), is not accurate. The vehicle system stops the estimation at such speeds, to avoid propagating this inaccuracy in the road gradient estimation, (e.g., RGEkin and RGEdyn), and locks the estimator state to its latest qualified estimation result.
A vehicle may travel an extended distance at such very low vehicle speeds. Therefore the actual road gradient may change after the RGE is stopped. In such conditions, the accuracy of the locked RGE output will decrease. Therefore the vehicle system downgrades the associated estimation quality accordingly. The vehicle system defines two vehicle speed thresholds, (Vstp and Vexit), which represent a vehicle stop speed threshold, and a vehicle exit speed threshold, respectively. Where Vstp is greater than Vexit. When the vehicle speed (Vx) is less than Vstp, The vehicle system (using a Kalman filter (or a least square type of filter) based estimation algorithm) halts or stops the estimation by locking all of the estimation states. However, after recording the data at the estimation stop moment, the estimation algorithm continues, before the vehicle speed (Vx) further decreases down to Vexit, or before Vx increases above Vstp. When operating in the speed range between Vexit<Vx<Vstp, the new estimation output (RGEout) will continue to compare the locked estimation output. If the difference between current RGE state and the locked RGE is greater than a first difference threshold, the estimation quality will be downgraded to a lower level indicating that the estimation output may not be accurate enough in this new situation. However, if the difference between them is greater than a second difference threshold, the estimation quality will be downgraded to a “No Data” level (e.g., QF=0) indicating that the estimation result cannot be used. After the vehicle speed (Vx) decreases further below Vexit, the estimation will be fully stopped and no more estimation activity will be carried out. Once the vehicle speed (Vs) increases above the Vstp threshold, the locked and recorded estimation states will be reloaded back to the estimation algorithm to resume normal estimation function.
The vehicle system stops the road gradient estimator once the vehicle speed (Vx) is less than Vexit, in one or more embodiments, because an estimation at such conditions would not be accurate. Vexit may be referenced as a “crawling” speed. However, the vehicle may continue to travel at such a crawling speed.
With reference to FIGS. 27A-27G, the kinematic or dynamic road gradient estimator will be available after the vehicle speed increases above a higher speed level and the static RGE estimator can take the responsibility after the vehicle is fully stopped. FIGS. 27A-27G include 7 graphs of measurement data taken over a common period of time and are generally referenced by numeral 2710. FIG. 27B illustrates vehicle speed (Vx); FIG. 27E illustrates the static quality factor (QFst); FIG. 27F illustrates the kinematic quality factor (QFkin); and FIG. 27G illustrates the quality factor of the arbitrated road gradient estimation (QFRGE). The static RGE is available when the vehicle is stopped. A vehicle speed of zero is referenced by numeral 2712 and a corresponding high quality RGEst is referenced by numeral 2714. Additionally the kinematic RGE is available when the vehicle is moving in normal motion as referenced by numeral 2716 and a corresponding high quality RGEkin is referenced by numeral 2718. The arbitrated RGE combines RGEst and RGEkin to provide a high quality output over most vehicle conditions, as referenced by numeral 2720. However, at certain very low vehicle speeds, as referenced by numeral 2722, there may be no high quality RGE available, as generally referenced by numeral 2724.
However, there may be no RGE algorithm available in the low speed range between 0 and Vexit. Thus, the road gradient variation while the vehicle is moving at such low speed cannot be detected and estimated. In order to assure RGE's robustness in reporting road gradient information, the vehicle system includes a conservative countermeasure to avoid a road gradient estimation being used without providing information regarding the quality of the RGE.
The vehicle system includes a method, (not shown) for monitoring vehicle speed (Vx) in a crawling event, according to one or more embodiments. A crawling event occurs when Vx is less than v_exit (e.g., less than 5 kph) and the current estimator quality evaluation is OK (e.g., QF=3). The vehicle system calculates a crawling distance (Dc), by integrating the vehicle speed. The vehicle system downgrades the quality factor to degraded (QF=2) when Dc is greater than a first distance threshold. The vehicle system downgrades the quality factor to no data (QF=0) when Dc is greater than a second distance threshold. The vehicle system resets the crawling distance (Dc) to zero, when any RGE estimator starts providing a qualified output. Meanwhile, the quality of RGE will be arbitrated to the current active estimator's quality evaluation result to replace the quality evaluated in the crawling event. The vehicle system will restart the crawling distance computation, when the quality evaluation (QF) downgrades in the speed range (v_exit<Vx<v_stp). Additionally, the vehicle system restarts the Dc computation when the continued estimation result deviates from the locked estimation result more than a certain threshold, even if the vehicle speed is greater than Vexit.
As such the vehicle system provides advantages over existing system by evaluating the quality of input signals, road gradient estimations and mass estimations and providing quality factors associated with each input and estimation. The vehicle system arbitrates or selects an estimation based in part on the quality factors. The vehicle system then provides a road gradient estimation and a vehicle mass estimation with a corresponding quality factor based on the selected estimation. Other vehicle systems may use the road gradient and vehicle mass estimations differently, depending on the associated quality factors.
As demonstrated by the embodiments described above, the vehicle, vehicle system, and method provide advantages over the prior art by arbitrating or selecting a road gradient estimation from a number of different estimations (static, kinematic, and dynamic). By arbitrating between these different estimations, the vehicle system provides a road gradient estimation over a wide range of vehicle conditions. The vehicle system also provides advantages over existing systems by estimating static road gradient using dynamically filtering of the longitudinal acceleration, and by using another estimation for the initial value; which reduces the time delay before a quality static road gradient estimated value may be provided. The vehicle system also provides advantages over the existing methods by estimating kinematic road gradient using an EKF, which provides an integral type filtering strategy. Such an integrating strategy does not amplify noise on the signals, and therefore provides improved estimation accuracy with relatively fast convergence. The vehicle system also provides advantages over the existing methods by estimating kinematic road gradient using input compensated acceleration offset to improve accuracy to remove uncertainties and offsets caused by lateral vehicle dynamics and their coupling effects. The vehicle system provides advantages by estimating vehicle mass independent of the road gradient estimation using an event-seeking strategy. The vehicle system re-initializes certain parameters in the mass estimation algorithm at vehicle conditions corresponding to potential mass changes (e.g., when the vehicle stops) in order to increase the sensitivity of the algorithm. The vehicle system estimates a dynamic road gradient (RGEdyn) based on the estimated vehicle mass. The vehicle system also provides advantages over existing system by evaluating the quality of input signals, road gradient estimations and mass estimations and providing quality factors associated with each input and estimation. The vehicle system arbitrates or selects an estimation based in part on the quality factors. The vehicle system then provides a road gradient estimation and a vehicle mass estimation with a corresponding quality factor based on the selected estimation. Other vehicle systems may use the road gradient and vehicle mass estimations differently, depending on the associated quality factors.
While the best mode has been described in detail, those familiar with the art will recognize various alternative designs and embodiments within the scope of the following claims. Additionally, the features of various implementing embodiments may be combined to form further embodiments of the invention. While various embodiments may have been described as providing advantages or being preferred over other embodiments or prior art implementations with respect to one or more desired characteristics, those of ordinary skill in the art will recognize that one or more features or characteristics may be compromised to achieve desired system attributes, which depend on the specific application and implementation. These attributes may include, but are not limited to: cost, strength, durability, life cycle cost, marketability, appearance, packaging, size, serviceability, weight, manufacturability, ease of assembly, etc. The embodiments described herein that are described as less desirable than other embodiments or prior art implementations with respect to one or more characteristics are not outside the scope of the disclosure and may be desirable for particular applications. Additionally, the features of various implementing embodiments may be combined to form further embodiments of the invention.

## Claims

1. A vehicle comprising:
a controller configured to generate output indicative of a vehicle mass estimation using a recursive least square estimation based on a longitudinal acceleration, an acceleration offset and a wheel torque when at least one of the longitudinal acceleration, a vehicle speed and a yaw rate indicate an occurrence of a qualified event;
wherein the vehicle mass estimation is independent of a road gradient estimation; and
wherein the acceleration offset includes at least one of a lateral velocity, a lateral offset, and a vehicle pitch angle.

2. The vehicle of claim 1 wherein the qualified event comprises at least one of an acceleration event and a braking event.

3. The vehicle of claim 2 wherein the acceleration event further comprises:
input indicative of the vehicle speed being greater than a propulsion entry speed threshold;
input indicative of the longitudinal acceleration being greater than a positive acceleration threshold; and
input indicative of a magnitude of the yaw rate being less than a yaw rate threshold.

4. The vehicle of claim 2 wherein the braking event further comprises:
input indicative of the vehicle speed being greater than a brake entry speed threshold;
input indicative of the longitudinal acceleration being less than a negative acceleration threshold; and
input indicative of a magnitude of the yaw rate being less than a yaw rate threshold.

5. The vehicle of claim 1 further comprising:
an inertial sensor for providing a signal indicative of the longitudinal acceleration.

6. The vehicle of claim 1 wherein the controller is further configured to:
set the output indicative of the vehicle mass estimation to a default value until the occurrence of a predetermined number of qualified events.

7. The vehicle of claim 1 wherein the controller is further configured to:
limit the vehicle mass estimation using a dynamic upper boundary having an initial value based on a maximum gross vehicle weight.

8. The vehicle of claim 7 wherein the controller is further configured to:
limit the vehicle mass estimation using a dynamic lower boundary having an initial value based on a vehicle curb weight.

9. The vehicle of claim 8 wherein at least one of the dynamic upper boundary and the dynamic lower boundary is configured to converge in response to an accumulation of qualified events.

10. The vehicle of claim 1 wherein the controller is further configured to:
generate output indicative of a dynamic road gradient estimation based on the vehicle mass estimation, the vehicle speed and the wheel torque.

11. A vehicle system comprising:
a controller configured to:
generate output indicative of a vehicle mass estimation based on a longitudinal acceleration and a wheel torque in response to a qualified event, and
generate output indicative of a dynamic road gradient estimation based on a vehicle speed, a wheel torque and the vehicle mass estimation.

12. The vehicle system of claim 11 wherein the controller is further configured to:
set the output indicative of the vehicle mass estimation to a default value until the occurrence of a predetermined number of qualified events.

13. The vehicle system of claim 11 wherein the controller is further configured to:
increment an event count in response to detecting a qualified event; and
generate the output indicative of the vehicle mass estimation when the event count is greater than a qualified event threshold.

14. The vehicle system of claim 11 wherein the controller is further configured to:
limit the vehicle mass estimation using a dynamic upper boundary having an initial value based on a maximum gross vehicle weight, and a dynamic lower boundary having an initial value based on a vehicle curb weight.

15. The vehicle system of claim 14 wherein at least one of the dynamic upper boundary and the dynamic lower boundary is configured to converge in response to an accumulation of qualified events.

16. The vehicle system of claim 11 wherein the controller is further configured to:
provide a dynamic quality factor that is indicative of an availability of the dynamic road gradient estimation.

17. The vehicle system of claim 11 wherein the vehicle mass estimation includes a vehicle payload.

18. A method for estimating a dynamic road gradient comprising:
generating output indicative of a vehicle mass estimation based on a longitudinal acceleration and a wheel torque in response to a qualified event; and
generating output indicative of a dynamic road gradient estimation based on a vehicle speed, a wheel torque and the vehicle mass estimation.

19. The method of claim 18 further comprising:
setting the output indicative of the vehicle mass estimation to a default value until the occurrence of a predetermined number of qualified events.