SUSPENSION STROKE RELATED VALUE ESTIMATION DEVICE AND METHOD

20220118812 · 2022-04-21

Assignee

Inventors

Cpc classification

International classification

Abstract

A control unit that controls control currents supplied to shock absorbers that generate damping forces according to the control currents determines coefficients of two functions that functionally identify an equivalent damping coefficient and an equivalent spring constant of each suspension based on relationships between the control current supplied to each shock absorber and the coefficients of the two functions that change according to the control current and a frequency of a relative vibration between a sprung and an unsprung of a vehicle, and calculates a relative displacement or a relative velocity between the sprung and the unsprung based on a vertical acceleration detected by a detection device and the two functions in which the coefficients are determined.

Claims

1. A suspension stroke related value estimation device for a vehicle which comprises a damping force generator configured to generate a damping force that damps a relative vibration between a sprung and an unsprung at each wheel position according to a supplied control current, a control unit configured to control a control current supplied to the damping force generator, and a detection device configured to detect a vertical acceleration of the sprung, and is configured to calculate at least one of a relative displacement and a relative velocity between the sprung and the unsprung as a stroke related value for a corresponding suspension, wherein the control unit includes a storage device that stores relationships between a control current and coefficients of two functions for functionally identifying an equivalent damping coefficient and an equivalent spring constant of each suspension that change according to the control current and a frequency of the relative vibration, the control unit is configured to determine the coefficients of the two functions based on the control current supplied to the damping force generator and the relationships, and calculate at least one of the relative displacement and relative velocity between the sprung and the unsprung based on a vertical acceleration of the sprung detected by the detection device and the two functions in which the coefficients are determined.

2. The suspension stroke related value estimation device for a vehicle according to claim 1, wherein the control unit is configured to specify a function for functionally identifying the equivalent spring constant and a function for functionally identifying the equivalent damping coefficient based on the determined coefficients, and calculate a relative displacement between the sprung and the unsprung as a value obtained by dividing a product of a mass of the sprung and a vertical acceleration of the sprung by a sum of a product of the specified function for functionally identifying the equivalent damping coefficient and the Laplace operator and the specified function for functionally identifying the equivalent spring constant.

3. The suspension stroke related value estimation device for a vehicle according to claim 2, wherein the control unit is configured to regard a transfer function from a vertical acceleration of the sprung to a relative displacement between the sprung and the unsprung as a product of a first transfer function from a vertical acceleration of the sprung to an intermediate variate and a second transfer function from the intermediate variate to a relative displacement between the sprung and the unsprung to calculate a relative displacement between the sprung and the unsprung.

4. The suspension stroke related value estimation device for a vehicle according to claim 3, wherein the function for functionally identifying the equivalent damping coefficient is a function of at least a first-order lag of a basic equivalent damping coefficient, and the function for functionally identifying the equivalent spring constant is a function of at least a first-order lag and at least a first-order lead of a basic equivalent spring constant.

5. The suspension stroke related value estimation device for a vehicle according to claim 4, wherein the control unit is configured to calculate a relative displacement between the sprung and the unsprung as a sum of the intermediate variate, a product of a sum of a time constant of the function of the first-order lag and a delay time constant of the function of the first-order lag and first-order lead and a derivative value of the intermediate variate, and a product of the time constant of the function of the first-order lag, the delay time constant of the function of the first-order lag and first-order lead and a second derivative value of the intermediate variate.

6. The suspension stroke related value estimation device for a vehicle according to claim 1, wherein the control unit is configured to specify a function for functionally identifying the equivalent spring constant and a function for functionally identifying the equivalent damping coefficient based on the determined coefficients, and calculate a relative velocity between the sprung and the unsprung as a value obtained by dividing a product of a mass of the sprung, a vertical acceleration of the sprung and the Laplace operator by a sum of a product of the specified function for functionally identifying the equivalent damping coefficient and the Laplace operator and the specified function for functionally identifying the equivalent spring constant.

7. The suspension stroke related value estimation device for a vehicle according to claim 6, wherein the control unit is configured to regard a transfer function from a vertical acceleration of the sprung to a relative velocity between the sprung and the unsprung as a product of a third transfer function from a vertical acceleration of the sprung to an intermediate variate and a fourth transfer function from the intermediate vanate to a relative velocity between the sprung and the unsprung to calculate a relative velocity between the sprung and the unsprung.

8. The suspension stroke related value estimation device for a vehicle according to claim 7, wherein the function for functionally identifying the equivalent damping coefficient is a function of at least a first-order lag of a basic equivalent damping coefficient, and the function for functionally identifying the equivalent spring constant is a function of at least a first-order lag and at least a first-order lead of a basic equivalent spring constant.

9. The suspension stroke related value estimation device for a vehicle according to claim 8, wherein the control unit is configured to calculate a relative velocity between the sprung and the unsprung as a sum of a product of a sum of a time constant of the function of the first-order lag and a delay time constant of the function of the first-order lag and first-order lead and a derivative value of the intermediate variate, and a product of the time constant of the function of the first-order lag, the delay time constant of the function of the first-order lag and first-order lead and a second derivative value of the intermediate variate.

10. A suspension stroke related value estimation method for a vehicle which calculates at least one of a relative displacement and a relative velocity between a sprung and an unsprung of a vehicle as a suspension stroke related value using a damping force generator configured to generate a damping force that damps a relative vibration between the sprung and the unsprung according to a supplied control current, a storage device that stores relationships between a control current and coefficients of two functions that functionally identify an equivalent damping coefficient and an equivalent spring constant of each suspension that change according to the control current and a frequency of the relative vibration, and a detection device that detects a vertical acceleration of the sprung, comprising: a step of acquiring a vertical acceleration of the sprung detected by the detection device, a step of acquiring a control current supplied to the damping force generator, a step of determining coefficients of the two functions based on an acquired control current and the relationships stored in the storage device, and a step of calculating at least one of a relative displacement and a relative velocity between the sprung and the unsprung based on an acquired vertical acceleration of the sprung and the two functions for which the coefficients determined.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0037] FIG. 1 is a schematic configuration diagram showing a first embodiment of a damping control apparatus to which the suspension stroke related value estimation device according to the present disclosure is applied.

[0038] FIG. 2 is a flowchart showing a damping control routine according to the first embodiment.

[0039] FIG. 3 is a block diagram of calculation of a relative displacement z.sub.si performed in step 240 in the flowchart shown in FIG. 2.

[0040] FIG. 4 is a flowchart showing a damping control routine according to the second embodiment.

[0041] FIG. 5 is a block diagram of calculation of a relative velocity dz.sub.si performed in step 440 of the flowchart shown in FIG. 4.

[0042] FIG. 6 is a diagram showing an example of a frequency response of an equivalent damping coefficient c.sub.e of a suspension for various control currents 1.

[0043] FIG. 7 is a diagram showing an example of a frequency response of an equivalent spring constant k.sub.e of a suspension for various control currents 1.

[0044] FIG. 8 is a diagram showing a single-wheel model of a vehicle for explaining a principle of estimation in the suspension stroke related value estimation device and method of the present disclosure.

[0045] FIG. 9 is a map for determining a coefficient c.sub.o based on a control current I.

[0046] FIG. 10 is a map for determining a coefficient a.sub.c based on a control current I.

[0047] FIG. 11 is a map for determining a coefficient k.sub.o based on a control current I.

[0048] FIG. 12 is a map for determining a coefficient a.sub.k based on a control current I.

[0049] FIG. 13 is a map for determining a coefficient b.sub.k based on a control current I.

[0050] FIG. 14 is a diagram showing an example of a gain characteristic of transfer from a true value to an estimated value of a stroke velocity of a suspension in the embodiment (thick solid line) and in a prior art (thin solid line).

[0051] FIG. 15 is a diagram showing an example of a phase characteristic of transfer from a true value to an estimated value of a stroke velocity of a suspension in the embodiment (thick solid line) and in a prior art (thin solid line).

DETAILED DESCRIPTION

[Principle of the Present Disclosure Adopted in the Embodiment]

[0052] In order to facilitate the understanding of the present disclosure, a principle of the suspension stroke related value estimation device and the method of the present disclosure will be described prior to the description of the embodiment.

<Equivalent Damping Coefficient and Equivalent Spring Constant>

[0053] Although not shown in the drawing, a random vibration test was performed in which a control current I supplied to a variate damping force shock absorber was set to various values and each wheel was vibrated by a four-wheel shaker. A vertical displacement z.sub.w and a vertical velocity dz.sub.w of each wheel were measured. The broken lines in FIGS. 6 and 7 show examples of an equivalent damping coefficient c.sub.e(I) and an equivalent spring constant ke(I) calculated according to the following equations (1) and (2) based on the test results, respectively. The horizontal axis of FIGS. 6 and 7 is a frequency (logarithmic scale) of a relative vibration between a sprung and an unsprung. Therefore, it can be seen that the equivalent damping coefficient c.sub.e(I) and the equivalent spring constant k.sub.e(I) change according to the control current I and the frequency of the relative vibration.

[00001] c e ( I ) = F w dz w ( 1 ) k e ( I ) = F w z w ( 2 )

[0054] The solid lines in FIGS. 6 and 7 show an equivalent damping coefficient c.sub.e(I) and an equivalent spring constant k.sub.e(I), respectively in which the test results shown by the broken lines are approximated to a gentle curve. The following equations (3) and (4) show functions of the an equivalent damping coefficients c.sub.2e(I) and an equivalent spring constant k.sub.2e(I) that function-identify the equivalent damping coefficient c.sub.e(I) and the equivalent spring constant k.sub.e(I) shown by the solid lines in FIGS. 6 and 7, respectively. In the equations (3) and (4) and equations described below, s is the Laplace operator. Therefore, the equations (3) and (4) are functions for functionally identifying the equivalent damping coefficient and the equivalent spring constant, respectively of the suspension, which change according to the control current I and the frequency of the relative vibration.

[00002] c 2 e ( I ) = 1 1 + a c ( I ) s c 0 ( I ) ( 3 ) k 2 e ( I ) = 1 + b k ( I ) s 1 + a k ( I ) s k 0 ( I ) ( 4 )

[0055] The equation (3) is a function of a first-order lag of a basic equivalent damping coefficient c.sub.o(I), and a.sub.c(I) is a time constant of the first-order lag. The equation (4) is a function of a first-order lag and first-order lead of a basic equivalent spring constant k.sub.o(I), a.sub.k(I) is a time constant of the first-order lag, and b.sub.k(I) is a time constant of the first-order lead. Notably, a vertical displacement z.sub.w and a vertical velocity dz.sub.w of each wheel differ depending on specifications of the vehicle. However, the equivalent damping coefficient c.sub.2e(I) can be identified by a function of at least a first-order lag of the basic equivalent damping coefficient c.sub.o(I) regardless of the specifications of the vehicle. Similarly, the equivalent spring constant k.sub.2e(I) can be identified by a function of at least first-order lag and at least first-order lead of the basic equivalent spring constant k.sub.o(I) regardless of the specifications of the vehicle.

[0056] In the present application, the basic equivalent damping coefficient c.sub.o(I) and the time constant a.sub.c(I) are referred to as coefficients of the function of the equivalent damping coefficient c.sub.e(I), and the basic equivalent spring constant k.sub.o(I) and the time constants a.sub.k(I) and b.sub.k(I) are referred to as coefficients of the function of the equivalent spring constant k.sub.e(I). That is, the values other than the Laplace operator s in the functions of the equivalent damping coefficient c.sub.2e(I) and the equivalent spring constant k.sub.2e(I) represented by the equations (3) and (4) are referred to as coefficients of the two functions. The coefficients c.sub.o(I), a.sub.c(I), k.sub.o(I), a.sub.k(I) and b.sub.k(I) are functions of the control current I. and are obtained, for example, by referring to the maps shown in FIGS. 9 to 13. It is to be noted that in FIGS. 9 to 13, (I) of each coefficient is omitted, and (I) of each coefficient is also omitted in the following description as necessary.

<Vehicle Model>

[0057] FIG. 8 shows a single-wheel model of the vehicle 100 for explaining the principle of estimation in the suspension stroke related value estimation device and method of the present disclosure. In FIG. 8, the reference numerals 102 and 104 show the unsprung and the sprung of the vehicle 100, respectively. A suspension spring 106 and a damper (damping force generator or the like) 108 are provided between the unsprung 102 and the sprung 104. A spring 110 indicates a spring of the unsprung 102.

<Calculation of Relative Displacement>

[0058] As shown in FIG. 8, the masses of the unsprung 102 and the sprung 104 are represented by m.sub.1 and m.sub.2, respectively, and their vertical displacements are represented by z.sub.1 and z.sub.2, respectively. A spring constant of the spring 106 and a damping coefficient of the damper 108 are represented by k.sub.2 and c.sub.2, respectively. A spring constant of the spring 110 is represented by k.sub.1, and a vertical displacement of a ground contact portion 112 of the unsprung 102, that is, a road surface 114, is represented by z.sub.0.

[0059] Further, derivative values of the vertical displacements z.sub.1 and z.sub.2, that is, vertical velocities are represented by dz.sub.1 and dz.sub.2, respectively, and second-order derivative values of the vertical displacements z.sub.1 and z.sub.2, that is, vertical accelerations are represented by ddz.sub.1 and ddz.sub.2, respectively. It is assumed that an upward displacement is positive for z.sub.1 and z.sub.2, and an upward force is positive for forces generated by the spring 106, the damper 108 and the like.

[0060] Equations of motion for vertical movements of the sprung 104 and the unsprung 102 of the vehicle 100 are expressed by the following equations (5) and (6), respectively.


m.sub.2ddz.sub.2=k.sub.2e(z.sub.1−z.sub.2)+c.sub.2e(dz.sub.1−dz.sub.2)  (5)


m.sub.1ddz.sub.1=−k.sub.2e(z.sub.1−z.sub.2)−c.sub.2e(dz.sub.1−dz.sub.2)+k.sub.1(z.sub.0−z.sub.1)  (6)

[0061] From the above equation (5), a relative displacement z.sub.s between the sprung 104 and the unsprung 102 is expressed by the following equation (7). Therefore, by detecting a vertical acceleration ddz.sub.2 of the sprung 104 and determining the coefficients of the above equations (3) and (4) based on a control current I, the equivalent damping coefficient c.sub.2e and the equivalent spring constant k.sub.2e are specified. Thus, the relative displacement z.sub.s can be calculated according to the following equation (7) Notably, the mass m.sub.2 of the sprung 104 is known. Further, there is a certain relationship between the coefficients of the above equations (3) and (4) and the control current I, and the relationships can be experimentally obtained as described later.

[00003] z s = z 1 - z 2 = 1 k 2 e + c 2 e s m 2 ddz 2 ( 7 )

[0062] By substituting the above equations (3) and (4) into the above equation (7), the relative displacement z.sub.s between the sprung 104 and the unsprung 102 is expressed by the following equation (8).

[00004] z s = 1 1 + b k ( I ) s 1 + a k ( I ) s k 0 ( I ) + 1 1 + a c ( I ) s c 0 ( I ) s m 2 ddz 2 = ( 1 + a k s ) ( 1 + a c s ) ( 1 + b k s ) ( 1 + a c s ) k 0 + ( 1 + a k s ) c 0 s m 2 ddz 2 = 1 + ( a k + a c ) s + a k a c s 2 k 0 + [ ( b k + a c ) k 0 + c 0 ] s + ( a c b k k 0 + a k c 0 ) s 2 m 2 ddz 2 ( 8 )

[0063] Assuming that a transfer function from a vertical acceleration ddz.sub.2 of the sprung 104 to a relative displacement z.sub.s is a product of a first transfer function y/ddz.sub.2 from the vertical acceleration ddz.sub.2 to an intermediate variate y and a second transfer function z.sub.s/y from the intermediate variate y to the relative displacement z.sub.s, a transfer function from the vertical acceleration ddz.sub.2 to the relative displacement z.sub.s is expressed by the following equation (9).

[00005] z s ddz 2 = y ddz 2 .Math. z s y = 1 + ( a k + a c ) s + a k a c s 2 k 0 + [ ( b k + a c ) k 0 + c 0 ] s + ( a c b k k 0 + a k c 0 ) s 2 m 2 ( 9 )

[0064] It is assumed that the second transfer function z.sub.s/y from the intermediate variate y to the relative displacement z.sub.s is expressed by the following equation (10), and the first transfer function y/ddz.sub.2 from the vertical acceleration ddz.sub.2 to the intermediate variate y is expressed by the equation (11).

[00006] z s y = 1 + ( a k + a c ) s + a k a c s 2 ( 10 ) y ddz 2 = m 2 k 0 + [ ( b k + a c ) k 0 + c 0 ] s + ( a c b k k 0 + a k c 0 ) s 2 ( 11 )

[0065] Representing a derivative value and a second derivative value of the intermediate variate y as dy and ddy, respectively, the relative displacement z.sub.s is expressed by the following equation (12) corresponding to the above equation (10). Further, since the following equation (13) is established from the above equation (11), the second derivative ddy of the intermediate variate y is expressed by the following equation (14).

[00007] z s = y + ( a k + a c ) dy + a k a c ddy ( 12 ) m 2 ddz 2 = k 0 y + [ ( b k + a c ) k 0 + c 0 ] dy + ( a c b k k 0 + a k c 0 ) ddy ( 13 ) ddy = 1 a c b k k 0 + a k c 0 { m 2 ddz 2 - k 0 y - [ ( b k + a c ) k 0 + c 0 ] dy } ( 14 )

[0066] In the first embodiment described later, as will be described in detail later, a vertical acceleration ddz.sub.2 of the sprung 104 is detected, the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k are calculated based on a control current I, and a relative displacement (suspension stroke) z.sub.s is calculated using the above equations (12) and (14). If necessary, a relative velocity (suspension stroke velocity) dz.sub.s may be calculated as a derivative value of the relative displacement z.sub.s.

[0067] Further, since the following equation (15) corresponding to the above equation (7) is established, a relative velocity dz.sub.s can be calculated according to the equation (15).

[00008] dz s = dz 1 - dz 2 = s k 2 e + c 2 e s m 2 ddz 2 ( 15 )

[0068] Assuming that a transfer function from the vertical acceleration ddz.sub.2 to the relative velocity dz.sub.s of the sprung 104 is a product of a third transfer function y/ddz.sub.2 from the vertical acceleration ddz.sub.2 to the intermediate variate y and a fourth transfer function from the intermediate variate y to the relative velocity dz.sub.s, a transfer function from the vertical acceleration ddz.sub.2 to a relative velocity dz.sub.s is expressed by the following equation (16).

[00009] dz s ddz 2 = y ddz 2 .Math. dz s y = 1 + ( a k + a c ) s + a k a c s 2 k 0 + [ ( b k + a c ) k 0 + c 0 ] s + ( a c b k k 0 + a k c 0 ) s 2 m s s ( 16 )

[0069] It is assumed that the fourth transfer function dz.sub.s/y from the intermediate variate y to the relative velocity dz.sub.s is expressed by the following equation (17), and the third transfer function y/ddz.sub.2 from the vertical acceleration ddz.sub.2 to the intermediate variate y is expressed by the equation (18).

[00010] dz s y = 1 + ( a k + a c ) s + a k a c s 2 ( 17 ) y ddz 2 = m 2 s k 0 + [ ( b k + a c ) k 0 + c 0 ] s + ( a c b k k 0 + a k c 0 ) s 2 ( 18 )

[0070] The relative velocity dz.sub.s is expressed by the following equation (19) corresponding to the above equation (17). Further, since the following equation (20) is established from the above equation (18), the second derivative ddy of the intermediate variate y is expressed by the following equation (21).

[00011] dz s = y + ( a k + a c ) dy + a k a c ddy ( 19 ) m 2 sddz 2 = k 0 y + [ ( b k + a c ] k 0 + c 0 ] dy + ( a c b k k 0 + a k c 0 ) ddy ( 20 ) ddy = 1 a c b k k 0 + a k c 0 { m 2 sddz 2 - k 0 y - [ ( b k + a c ) k 0 + c 0 ] dy } ( 21 )

[0071] In the second embodiment described later, as will be described in detail later, a vertical acceleration ddz.sub.2 of the sprung 104 is detected, the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k are calculated based on a control current I. and a relative velocity (suspension stroke velocity) dz.sub.s is calculated using the above equations (19) and (21). If necessary, a relative displacement (suspension stroke) z.sub.s may be calculated as an integral value of the relative velocity dz.sub.s.

Embodiment

[0072] The present disclosure will now be described in detail with reference to the accompanying drawings.

<Configuration Common to the First and Second Embodiments>

[0073] As shown in FIG. 1, the suspension stroke related value estimation device 10 of the first and second embodiments is applied to a damping control apparatus 14 of a vehicle 11. The vehicle 11 has left and right front wheels 12FL and 12FR which are steering wheels, and left and right rear wheels 12RL and 12RR which are non-steering wheels. Further, the vehicle 11 has front wheel suspensions 18FL and 18FR that suspend the front wheels 12FL and 12FR from a vehicle body 16, respectively and rear wheel suspensions 18RL and 18RR that suspend the rear wheels 12RL and 12RR from the vehicle body 16, respectively.

[0074] The front wheel suspensions 18FL and 18FR include suspension arms 20FL and 20FR, respectively, and the rear wheel suspensions 18RL and 18RR include suspension arms 20RL and 20RR, respectively. In FIG. 1, only one suspension arm 20FL to 20RR is shown, but a plurality of these arms may be provided.

[0075] The front wheels 12FL and 12FR are rotatably supported around rotation axes (not shown) by wheel support members 22FL and 22FR, respectively, and the wheel support members 22FL and 22FR are connected to the vehicle body 16 by the suspension arms 20FL and 20FR, respectively. Similarly, the rear wheels 12RL and 12RR are rotatably supported around rotation axes (not shown) by wheel support members 22RL and 22RR, respectively, and the wheel support members 22RL and 22RR are connected to the vehicle body 16 by the suspension arms 20RL and 20RR, respectively.

[0076] The front wheel suspensions 18FL and 18FR include shock absorbers 24FL and 24FR which function as damping force generators and suspension springs 26FL and 26FR, respectively. Similarly, the rear wheel suspensions 18RL and 18RR include shock absorbers 24RL and 24RR which function as damping force generators and suspension springs 26RL and 26RR, respectively.

[0077] In the illustrated embodiment, each shock absorber is connected to the vehicle body 16 at an upper end and to the corresponding wheel support member at a lower end. Each suspension spring is mounted between the vehicle body 16 and the wheel support member via the corresponding shock absorber. However, each shock absorber and suspension spring may be arranged between the vehicle body 16 and the wheel support member or the suspension arm.

[0078] The suspensions 18FL to 18RR may be any type of suspension as long as the wheels 12FL to 12RR and the vehicle body 16 are allowed to be displaced in the vertical direction with respect to each other. Further, each suspension spring 26FL to 26RR may be arbitrary springs such as a compression coil spring and an air spring.

[0079] As can be seen from the above description, at least the vehicle body 16 constitutes a sprung of the vehicle 11, and at least the wheels 12FL to 12RR and the wheel support members 22FL to 22RR constitute an unsprung of the vehicle 11.

[0080] The shock absorbers 24FL to 24RR are shock absorbers of variate damping force type having actuators 28FL to 28RR, respectively. Each actuator is configured to vary an effective passage cross-sectional area of a variate orifice provided on a piston (not shown) to vary a damping coefficient by means of a control current supplied to the actuator being controlled by an electronic control unit 30. The shock absorbers 24FL to 24RR may be shock absorbers having any configuration known in the art as long as a damping coefficient can be changed by controlling a control current.

[0081] As shown in FIG. 1, the vehicle body 16 is provided with vertical acceleration sensors 32FL, 32FR, 32RL and 32RR, respectively that detect vertical accelerations ddz.sub.2i of the sprung at positions corresponding to the left and right front wheels 12FL and 12FR and the left and right rear wheels 12RL and 12RR. Each signal indicating the vertical acceleration ddz.sub.2i is input to the electronic control unit 30. In addition, i is fl, fr, ri and rr that mean the left front wheel, the right front wheel, the left rear wheel and the right rear wheel, respectively.

[0082] Although not shown in detail in FIG. 1, the electronic control unit 30 includes a microprocessor and a drive circuit. The microcomputer has a CPU, a ROM, a RAM, and an input/output port device, and has a general configuration in which these are connected to each other by a bidirectional common bus. A control program for controlling the shock absorbers 24FL to 24RR is stored in the ROM functioning as a storage device, and the shock absorbers 24FL to 24RR are controlled by the CPU according to the control program.

[0083] Further, the ROM functions as a storage device for storing maps for calculating coefficients such as c.sub.o and a.sub.c, that is, maps shown in FIGS. 9 to 13. Further, the ROM stores a map for calculating control currents to be supplied to the actuators 28FL to 28RR, that is, target control currents Iti based on relative velocities dz.sub.si. Coefficients such as c.sub.o and a.sub.c may be values common to all suspensions. However, in general, specifications of the front wheel suspensions 18FL and 18FR are the same as each other and specifications of the rear wheel suspensions 18RL and 18RR are the same as each other, but the specifications of the front wheel suspensions and the specifications of the rear wheel suspensions are different from each other. Therefore, in the first and second embodiments, coefficients such as c.sub.o and a.sub.c are calculated for the front wheels and the rear wheels.

<Damping Control in the First Embodiment>

[0084] As will be described in detail later, in the first embodiment, the electronic control unit 30 calculates a relative displacement z.sub.si and a relative velocity dz.sub.si between the sprung and the unsprung for each wheel according to the flowchart shown in FIG. 2 and a control program corresponding to the block diagram shown in FIG. 3. Further, the electronic control unit 30 calculates target control currents Iti based on the relative velocities dz.sub.si, and controls control currents Ii supplied to the actuators so that the control currents Ii become the corresponding target control currents Iti.

[0085] Next, the damping control in the first embodiment will be described with reference to the flowchart shown in FIG. 2. The damping control according to the flowchart shown in FIG. 2 is repeatedly executed by the CPU of the electronic control unit 30 at predetermined time intervals in the order of, for example, the left front wheel, the right front wheel, the left rear wheel, and the right rear wheel when an ignition switch (not shown in the figure) is ON. This also applies to the control by the flowchart shown in FIG. 4 described later.

[0086] First, in step 210, the CPU reads a signal indicating a vertical acceleration ddz.sub.2i of the sprung detected by corresponding one of the vertical acceleration sensors 32FL to 32RR and stores it in the RAM.

[0087] In step 220, the CPU reads a previous value Ifi of the control current supplied to corresponding one of the actuators 28FL to 28RR and stores it in the RAM.

[0088] In step 230, the CPU determines the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k by referring to the maps shown in FIGS. 9 to 13 based on the previous value Ifi of the control current, and stores the coefficients in the RAM. Notably, when the previous value Ifi of the control current is Iffl or Iffr, the coefficients for the corresponding front wheel are calculated, and when the previous value Ifi of the control current is Ifl or Ifr, the coefficients for the corresponding rear wheel are calculated.

[0089] In step 240, as will be described in detail later, the CPU calculates a relative displacement z.sub.si between the sprung and the unsprung based on the vertical acceleration ddz.sub.2i of the sprung and the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k according to the block diagram shown in FIG. 3.

[0090] In step 250, the CPU calculates a relative velocity dz.sub.si between the sprung and the unsprung as a derivative value of the relative displacement z.sub.si.

[0091] In step 260, the CPU calculates a control current to be supplied to corresponding one of the actuators 28FL to 28RR, that is, a target control current Iti, from a map not shown in the drawing based on the relative velocity dz.sub.si.

[0092] In step 270, the CPU controls a control current Ii supplied to corresponding one of the actuators 28FL to 28RR so that the control current becomes the corresponding target control current Iti, thereby controlling a damping coefficient of corresponding one of the shock absorbers 24FL to 24RR.

[0093] As shown in FIG. 3, a vertical acceleration ddz.sub.2i of the sprung is input to a multiplier block 40, and a product m.sub.2ddz.sub.2i of the vertical acceleration ddz.sub.2i and the mass m.sub.2 of the sprung, which is the output of the multiplier block 40, is input to a positive input terminal of an adder block 42. A product k.sub.oy of the coefficient k.sub.o and the intermediate variate y is input to a negative input terminal of the adder block 42 from a multiplier block 44 described later. Further, a product [(bk+ac) ko+co] dy of a value (bk+ac) ko+co in a square bracket on the right side of the above equation (14) and a derivative value dy of the intermediate variate y is input to the other negative input terminal of the addition block 42.

[0094] An output of the adder block 42 is a value in a curly brace on the right side of the above equation (14) and is input to a multiplier block 48 so that it is multiplied by a coefficient 1/(a.sub.cb.sub.kk.sub.o+a.sub.kc.sub.o) for a curly brace on the right side of the above equation (14) by the multiplier block 48. Therefore, an output of the multiplier block 48 is a value on the right side of the above equation (14), and thus represents a second-order derivative ddy of the intermediate variate y.

[0095] The output of the multiplier block 48 is input to an adder block 54 as the intermediate variate y via integrator blocks 50 and 52. Further, the output of the multiplier block 48 is input to a multiplier block 56, and a product a.sub.ka.sub.cddy of the second-order derivative value ddy of the intermediate variate y and a product a.sub.ka.sub.c of the coefficients, which is the output of the multiplier block 56, is input to the adder block 54. Further, an output of the integrator block 50, that is, a derivative value dy of the intermediate variate y is input to a multiplier block 58, and a product (a.sub.k+a.sub.c) dy of the derivative value dy of the intermediate variate y and a sum a.sub.k+a.sub.c of the coefficients, which is the output of the multiplier block 58, is input to the adder block 54.

[0096] Therefore, an output of the adder block 54 is a value on the right side of the above equation (12), and thus represents a relative displacement (suspension stroke) z.sub.si between the sprung and the unsprung. The output of the integrator block 50, that is, the derivative value dy of the intermediate variate y, is input to a multiplier block 46 which multiplies a value (b.sub.k+a.sub.c) k.sub.o+c.sub.o in the square bracket on the right side of the equation (14). The output of the integrator block 52, i.e. the intermediate variate y, is input to the above-mentioned multiplier block 44, which multiplies the coefficient k.sub.o.

[0097] As can be seen from the above description, the electronic control unit 30 cooperates with the vertical acceleration sensors 32FL to 32RR to function as a suspension stroke related value estimation device that estimates relative displacements z.sub.si and a relative velocities dz.sub.si between the sprung and the unsprung by executing steps 210 to 250.

[0098] In the first embodiment, the step of acquiring a vertical acceleration of the sprung and the step of acquiring a previous value of the control current supplied to the damping force generator in the estimation method of the present disclosure are achieved by steps 210 and 220, respectively. Further, the step of determining the coefficients is accomplished by steps 230, and the step of calculating at least one of the relative displacement and the relative velocity between the sprung and the unsprung is accomplished by steps 240 and 250.

<Damping Control in the Second Embodiment>

[0099] As will be described in detail later, in the second embodiment, the electronic control unit 30 calculates a relative velocity dz.sub.si between the sprung and the unsprung for each wheel according to the flowchart shown in FIG. 4 and a control program corresponding to the block diagram shown in FIG. 5. Further, also in this embodiment, the electronic control unit 30 calculates a target control current Iti based on the relative velocity dz.sub.si, and controls the control current Ii supplied to each actuator so that the control current Ii becomes the corresponding target control current Iti.

[0100] Next, the damping control in the second embodiment will be described with reference to the flowchart shown in FIG. 4.

[0101] As can be seen from the comparison between FIGS. 4 and 2, the CPU executes steps 410 to 430 in the same manner as in steps 210 to 230, respectively. Therefore, the description of steps 410 to 430 will be omitted.

[0102] In step 440, as will be described in detail later, the CPU calculates a relative velocity dz.sub.si between the sprung and the unsprung based on the vertical acceleration ddz.sub.2i of the sprung and the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k according to the block diagram shown in FIG. 5. In FIG. 5, the blocks corresponding to the blocks shown in FIG. 3 are assigned the same reference numerals as those assigned in FIG. 3.

[0103] The CPU does not execute a step corresponding to step 250. Steps 460 and 470 are performed in the same manner as steps 260 and 270, respectively. Therefore, the description of steps 460 and 470 will be omitted.

[0104] As can be seen from the comparison between FIGS. 5 and 3, a differentiator block 38 is provided. A vertical acceleration ddz.sub.2i of the sprung is input to the differentiator block 38, and a product sddz.sub.2i of the Laplace operator s and the vertical acceleration ddz.sub.2i of the sprung, which is an output of the differentiator block 38, is input to a multiplier block 40. An output of the multiplier block 40, that is, a product m.sub.2sddz.sub.2i of the mass m.sub.2 of the sprung and the sddz.sub.2i is input to a positive input terminal of an adder block 42. Therefore, an output of the adder block 42 is a value in the curly brace on the right side of the above equation (21).

[0105] An output of a multiplier block 48 is a value on the right side of the above equation (21), and thus represents a second-order derivative ddy of the intermediate variate y. An output of an integrator block 50 is a derivative dy of the intermediate variate y, and an output of an integrator block 52 is the intermediate variate y. Therefore, an output of the adder block 54 is a value on the right side of the above equation (19), and thus represents a relative velocity (suspension stroke velocity) dz.sub.si between the sprung and the unsprung. In the second embodiment, a relative displacement (suspension stroke) z.sub.s may be calculated by integrating the relative velocity dz.sub.s.

[0106] As can be seen from the above description, the electronic control unit 30 cooperates with the vertical acceleration sensors 32FL to 32RR to function as a suspension stroke related value estimation device that estimates relative velocities dz.sub.si between the sprung and the unsprung by executing steps 410 to 440.

[0107] In the second embodiment, the step of acquiring a vertical acceleration of the sprung and the step of acquiring previous values of control currents supplied to the damping force generators in the estimation method of the present disclosure are achieved by steps 410 and 420, respectively. Further, the step of determining the coefficients is accomplished by step 430, and the step of calculating at least one of a relative displacement and a relative velocity between the sprung and the unsprung is accomplished by step 440.

[0108] Further, as can be seen from the above description, according to the first and second embodiments, relationships between the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k of the two functions that functionally identify an equivalent damping coefficient c.sub.e(I) and an equivalent spring constant k.sub.e(I) of each suspension that change according to a control current I and a frequency of the relative vibration and the control current I are stored in the storage device (ROM). The coefficients of the two functions are determined based on the control currents supplied to the shock absorbers 24FL to 24RR serving as the damping force generators and the above relationship. Further, at least one of a relative displacement z.sub.s and a relative velocity dz.sub.s between the sprung and the unsprung is calculated based on a vertical acceleration ddz.sub.2 of the sprung detected by the vertical acceleration sensors 32FL to 32RR serving as the detection device and the two functions in which the coefficients are determined.

[0109] Therefore, it is possible to specify two functions so that damping due to friction between suspension parts and due to rubber bushes and the like, and therefore non-linearity of a damping force is taken into consideration, and non-linearity of a spring force of the suspension is also taken into consideration. Therefore, even in a region where a frequency of the relative vibration between the sprung and the unsprung is relatively high, a suspension stroke related value can be estimated with high accuracy based on the vertical acceleration of the sprung.

[0110] For example, FIGS. 14 and 15 show examples of gain characteristics and phase characteristics of transfer from a true value to an estimated value of a stroke velocity of a suspension in the embodiments (thick solid line) and in the prior art (thin solid line), respectively. In FIGS. 14 and 15, the broken lines indicate that the gain characteristic and the phase characteristic vary according to a frequency of a relative vibration.

[0111] In the prior art, that is, in the estimation methods of Patent Documents 1 and 2 described above, in a region where a frequency of a relative vibration between a sprung and an unsprung is relatively high, as the frequency increases, the gain decreases and a phase shift between the true and estimated values increases. The decrease of the gain and the phase shift also occur in a region lower than a resonance frequency of the unsprung, and increase as the frequency of the relative vibration increases. Therefore, it is not possible to estimate a stroke velocity of a suspension with high accuracy in a region where the frequency of the relative vibration is relatively high.

[0112] On the other hand, according to the first and second embodiments, even in a region where a frequency of the relative vibration is high around a resonance frequency of the unsprung, the gain does not decrease and the phase shift between a true value and an estimated value does not increase. Therefore, a stroke velocity of a suspension can be estimated with high accuracy even in a region where the frequency of the relative vibration is relatively high.

[0113] Although not shown in the drawing, a gain characteristic and a phase characteristic of a transfer from a true value to an estimated value of a suspension stroke are also the same as the characteristics shown in FIGS. 14 and 15, respectively. Therefore, according to the first and second embodiments, a suspension stroke can be estimated with high accuracy even in a region where the frequency of the relative vibration is relatively high.

[0114] Further, according to the first and second embodiments, suspension stroke related values can be estimated without using an observer, and it is not necessary to spend a lot of time and effort to construct the observer, so suspension stroke related values can be estimated more easily than with conventional methods.

[0115] Further, according to the first and second embodiments, the equivalent damping coefficient c.sub.2e(I) and the equivalent spring constant k.sub.2e(1) can be functionally identified by the corresponding functions (equations (3) and (4), respectively), and the relationships between the coefficients c.sub.o, a.sub.c, k.sub.o, a.sub.k and b.sub.k of those functions and the control current I can be stored in the ROM serving as a storage device. Further, the coefficients of the two functions are determined based on the control current I and the above relationships, and at least one of a relative displacement z.sub.s and a relative velocity dz.sub.s between a sprung and an unsprung can be calculated based on a vertical acceleration ddz.sub.2 of the sprung and the two functions in which the coefficients are determined.

[0116] Further, according to the first and second embodiments, suspension stroke related values can be calculated by detecting a vertical acceleration ddz.sub.2 of the sprung by the vertical acceleration sensors 32FL to 32RR and acquiring the control currents supplied to the shock absorbers 24FL to 24RR.

[0117] In particular, according to the first embodiment, a relative displacement z.sub.s between a sprung and an unsprung can be calculated by regarding a transfer function z.sub.s/ddz.sub.2 from a vertical acceleration of the sprung to a relative displacement between the sprung and the unsprung as a product of the first transfer function y/ddz.sub.2 and the second function z.sub.s/y.

[0118] On the other hand, according to the second embodiment, a relative velocity dz.sub.s between a sprung and an unsprung can be calculated by regarding a transfer function dz.sub.s/ddz.sub.2 from a vertical acceleration of the sprung to a relative velocity between the sprung and the unsprung as a product of the third transfer function y/ddz.sub.2 and the fourth transfer function dz.sub.s/y.

[0119] Further, according to the first and second embodiments, a suspension stroke related value can be easily calculated by the functions of the intermediate variate y (equations (12) and (19)).

[0120] Although the present disclosure has been described in detail with reference to specific embodiments, it will be apparent to those skilled in the art that the present disclosure is not limited to the above-described embodiments, and various other embodiments are possible within the scope of the present disclosure.

[0121] For example, in the first and second embodiments, the equivalent damping coefficient c.sub.e(I) is a function represented by the equation (3), that is, a function of the first-order lag of the basic equivalent damping coefficient c.sub.o(I). However, the equivalent damping coefficient may be a function of second-order or higher-order lag. Similarly, the equivalent spring constant k.sub.e(I) is a function represented by the equation (4), that is, a function of the first-order lag and first-order lead of the basic equivalent spring constant k.sub.o(I). However, the equivalent spring constant may be a function of second-order or higher-order lag and second-order or higher-order lead. Again in these modification examples, the coefficients of each function are calculated based on a previous value Ifi of the control current.

[0122] Further, in the first and second embodiments, a vertical acceleration ddz.sub.2i of the sprung at the position of each wheel is detected by the vertical acceleration sensors 32FL to 32RR provided corresponding to each wheel. However, any one of the vertical acceleration sensors 32FL, 32FR, 32RL and 32RR may be omitted and a vertical acceleration at the position where the vertical acceleration sensor is omitted may be estimated based on the detected values of the three vertical acceleration sensors. Further, a vertical acceleration ddz.sub.21 of the sprung at the position of each wheel may be calculated based on a pitch rate of the sprung detected by a pitch rate sensor and a roll rate of the sprung detected by a roll rate sensor.

[0123] Further, in the first and second embodiments, a damping coefficient of each shock absorber is controlled based on a calculated relative velocity dz.sub.si between the sprung and the unsprung. However, a relative displacement z.sub.si and/or a relative velocity dz.sub.si calculated by the suspension stroke related value estimation device of the present disclosure may be used for other vehicle controls. The other vehicle controls include, for example, damping control by actively controlling a control force acting between a sprung and an unsprung, posture control to control posture of a sprung, and suspension control to control ride comfort of a vehicle and so on.