Method for reducing drive shaft vibration of eco-friendly vehicle

Abstract

A method for reducing vibration of a drive shaft of an eco-friendly vehicle includes calculating a model velocity of the drive shaft, obtaining a vibration component based on a deviation between an actual velocity of the drive shaft and the calculated model velocity, and generating a vibration reduction compensation torque for reduction in vibration of the drive shaft from the vibration component.

Claims

1. A method for reducing vibration of a drive shaft of an eco-friendly vehicle, the method comprising: acquiring an actual velocity of the drive shaft from which a torque of a driving motor driving the vehicle is output; calculating a model velocity of the drive shaft; obtaining a vibration component based on a deviation between the acquired actual velocity of the drive shaft and the calculated model velocity; and generating a vibration reduction compensation torque for reducing vibration of the drive shaft from the vibration component, wherein the step of calculating includes: calculating a drive shaft output requirement torque based on a torque delivered to the drive shaft; estimating a drive shaft input torque, which is inputted to the drive shaft, using the actual velocity of the drive shaft; estimating a disturbance torque using the drive shaft output requirement torque and the drive shaft input torque; calculating a drive shaft model input torque, which is obtained by adding the disturbance torque to the drive shaft output requirement torque, using the estimated disturbance torque; and calculating the model velocity using a drive shaft model which receives the drive shaft model input torque as an input.

2. The method according to claim 1, wherein the drive shaft output requirement torque is calculated by subtracting a sum of a vehicle load torque and a brake torque required for a brake of a driving wheel from a sum of torques required for the driving motor, an engine, and a starter-generator directly connected to the engine, respectively, in a hybrid electric vehicle (HEV) of the eco-friendly vehicle using the driving motor and the engine as driving sources of the vehicle.

3. The method according to claim 1, wherein the drive shaft output requirement torque is calculated by subtracting a sum of a vehicle load torque and a brake torque required for a brake of a driving wheel from a torque required for the driving motor in an electric (EV) of the eco-friendly vehicle vehicle using the driving motor as a driving source of the vehicle.

4. The method according to claim 1, wherein the step of estimating the drive shaft input torque comprises estimating the drive shaft input torque using a transfer function of the drive shaft model and the actual velocity (co) of the drive shaft, wherein the drive shaft input torque (T.sub.acc) is estimated using the following Equation:
T.sub.acc=Q(s)/G.sub.m(s), in which Q(s) is a transfer function of a filter, and G.sub.m(s) is the transfer function of the drive shaft model.

5. The method according to claim 4, wherein the transfer function Q(s) of the filter is set such that an order of a numerator is less than or equal to an order of a denominator in Q(s)/G.sub.m(s).

6. The method according to claim 4, wherein a primary low pass filter (LPF) is used as the filter, and the transfer function of the filter is Q(s)=1/(s+1), in which T denotes a time constant of the filter.

7. The method according to claim 4, wherein the transfer function of the drive shaft model is G.sub.m(s)=1/J.sub.ms, in which J.sub.m denotes a moment of inertia of the drive shaft.

8. The method according to claim 4, wherein the disturbance torque (d) is estimated by using the filter (Q) and using the following Equation:
d=Q(s)T.sub.netT.sub.acc, in which Q(s) denotes the transfer function of the filter, T.sub.net denotes the drive shaft output requirement torque, and T.sub.acc denotes the drive shaft input torque, and a filter having the same transfer function as a transfer function of a filter used to estimate the drive shaft input torque is applied to the filter (Q) for estimating the disturbance torque (d).

9. The method according to claim 1, wherein the step of estimating the disturbance torque comprises calculating the disturbance torque to be a difference value between the drive shaft output requirement torque and the drive shaft input torque.

10. The method according to claim 9, wherein a filter for eliminating a noise component is applied to the drive shaft output requirement torque to estimate the disturbance torque (d) using the following Equation:
d=Q(s)T.sub.netT.sub.acc, in which Q(s) denotes a transfer function of the filter, T.sub.net denotes the drive shaft output requirement torque, and T.sub.acc denotes the drive shaft input torque.

11. The method according to claim 10, wherein a primary low pass filter (LPF) is used as the filter, and the transfer function of the filter is Q(s)=1/(s+1), in which denotes a time constant of the filter.

12. The method according to claim 1, wherein the calculating of the model velocity comprises calculating the model velocity (.sub.m) using the following Equation:
.sub.m=G.sub.m(s)T.sub.m, in which G.sub.m(s) denotes a transfer function of the drive shaft model, and T.sub.m denotes the drive shaft model input torque.

13. The method according to claim 12, wherein the transfer function of the drive shaft model is G.sub.m(s)=1/J.sub.ms, in which J.sub.r, denotes a moment of inertia of the drive shaft.

14. A method for reducing vibration of a drive shaft of an eco-friendly vehicle comprising: acquiring current vehicle state information; acquiring an actual velocity of the drive shaft from which a torque of a driving motor driving the vehicle is output; determining a model velocity calculation scheme corresponding to one of a torque-based model velocity calculation scheme and a wheel velocity-based model velocity calculation scheme of calculating a model velocity based on wheel velocity information of the vehicle according to the acquired vehicle state information; calculating a model velocity of the drive shaft by changing a current model velocity calculation scheme to the determined model velocity calculation scheme when the determined model velocity calculation scheme is different from the current model velocity calculation scheme; obtaining a vibration component based on a deviation between the acquired actual velocity of the drive shaft and the calculated model velocity; and generating a vibration reduction compensation torque for reduction in vibration of the drive shaft from the vibration component.

15. The method according to claim 14, wherein the vehicle state information is a gear stage of a transmission, wherein the torque-based model velocity calculation scheme is determined when a target gear stage is lower than or equal to a set stage when the gear stage shifts, and the wheel velocity-based model velocity calculation scheme is determined when a current gear stage is higher than the set stage.

16. The method according to claim 15, wherein the model velocity calculation scheme is determined according to the target gear stage in response to gear shifting, and then the current model velocity calculation scheme is changed to the determined model velocity calculation scheme after the termination of gear shifting.

17. The method according to claim 16, wherein, when the model velocity calculation scheme is changed from the wheel velocity-based model velocity calculation scheme to the torque-based model velocity calculation scheme, an actual velocity of the drive shaft corresponding to a point in time at which the model velocity calculation scheme is changed is used as an initial model velocity in the changed torque-based model velocity calculation scheme.

18. The method according to claim 14, wherein the vehicle state information is a vehicle speed, a value of an accelerator pedal position sensor (APS), or a value of a brake pedal sensor (BPS).

19. The method according to claim 18, wherein the torque-based model velocity calculation scheme is determined when the vehicle speed is less than or equal to a reference vehicle speed, and the wheel velocity-based model velocity calculation scheme is determined when the vehicle speed exceeds the reference vehicle speed.

20. The method according to claim 18, wherein the torque-based model velocity calculation scheme is determined when a rate of change or the amount of change of the value of the APS is greater than or equal to a reference value, and the wheel velocity-based model velocity calculation scheme is determined when the rate of change or the amount of change of the value of the APS is less than the reference value.

21. The method according to claim 18, wherein the torque-based model velocity calculation scheme is determined when a rate of change or the amount of change of the value of the BPS is greater than or equal to a reference value, and the wheel velocity-based model velocity calculation scheme is determined when the rate of change or the amount of change of the value of the BPS is less than the reference value.

22. The method according to claim 18, wherein, when the model velocity calculation scheme is changed from the wheel velocity-based model velocity calculation scheme to the torque-based model velocity calculation scheme, the model velocity calculated using the wheel velocity-based model velocity calculation scheme corresponding to a point in time at which the model velocity calculation scheme is changed is used as an initial model velocity in the changed torque-based model velocity calculation scheme.

23. The method according to claim 18, wherein, when the model velocity calculation scheme is changed from the torque-based model velocity calculation scheme to the wheel velocity-based model velocity calculation scheme, an initial model velocity (.sub.m) in the changed wheel velocity-based model velocity calculation scheme is obtained using the following Equation:
.sub.m=.sub.m.sub.m.sub._.sub.init+.sub.m.sub._.sub.offset, in which .sub.m denotes a model velocity calculated using the wheel velocity-based model velocity calculation scheme corresponding to a point in time at which the model velocity calculation scheme is changed, .sub.m.sub._.sub.init denotes a difference value between model velocities calculated using the torque-based model velocity calculation scheme and the wheel velocity-based model velocity calculation scheme, respectively, at a point in time at which the model velocity calculation scheme is changed, and .sub.m.sub._.sub.offset denotes a value set to vary from 0 to .sub.m.sub._.sub.init with a constant slope or pattern and is a velocity determined to be a value according to .sub.m.sub._.sub.init.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The above and other features of the present disclosure will now be described in detail with reference to certain exemplary embodiments thereof illustrated in the accompanying drawings which are given hereinbelow by way of illustration only, and thus are not limitative of the present invention, and wherein:

(2) FIG. 1 is a block diagram illustrating a system configuration of an electric vehicle (EV) according to the related art;

(3) FIG. 2 is a block diagram illustrating a system configuration of a hybrid electric vehicle (HEV) according to the related art;

(4) FIG. 3 is a diagram illustrating a control system for reduction in drive shaft vibration of an eco-friendly vehicle according to an embodiment in the present disclosure;

(5) FIG. 4 is a flowchart illustrating a model velocity calculation process in a method for reducing drive shaft vibration of the eco-friendly vehicle according to an embodiment in the present disclosure;

(6) FIG. 5 is a diagram illustrating a configuration of a control system for reduction in drive shaft vibration of the eco-friendly vehicle capable of selecting a model velocity calculation scheme according to another embodiment in the present disclosure;

(7) FIG. 6 is a flowchart illustrating a process of selecting a model velocity calculation scheme and a model velocity according to a gear stage according to another embodiment in the present disclosure;

(8) FIG. 7 is a flowchart illustrating a process of selecting a model velocity calculation scheme and a model velocity according to a vehicle speed, a value of an accelerator pedal position sensor (APS), and a value of a brake pedal sensor (BPS) according to another embodiment in the present disclosure; and

(9) FIGS. 8 and 9 are diagrams illustrating a method of calculating a model velocity initial value at the time of change to a wheel velocity-based model velocity calculation scheme according to vehicle state information such as a vehicle speed, a value of an APS, and a value of a BPS according to another embodiment in the present disclosure.

(10) It should be understood that the appended drawings are not necessarily to scale, presenting a somewhat simplified representation of various preferred features illustrative of the basic principles of the invention. The specific design features of the present invention as disclosed herein, including, for example, specific dimensions, orientations, locations, and shapes will be determined in part by the particular intended application and use environment.

(11) In the figures, reference numbers refer to the same or equivalent parts of the present invention throughout the several figures of the drawing.

DETAILED DESCRIPTION

(12) Hereinafter reference will now be made in detail to various embodiments, examples of which are illustrated in the accompanying drawings and described below. While the invention will be described in conjunction with exemplary embodiments, it will be understood that present description is not intended to limit the invention to those exemplary embodiments. On the contrary, the invention is intended to cover not only the exemplary embodiments, but also various alternatives, modifications, equivalents, and other embodiments, which may be included within the spirit and scope of the invention as defined by the appended claims.

(13) The present disclosure relates to a method for reducing drive shaft vibration of an eco-friendly vehicle in which vibration generated from a drive shaft is reduced using a motor mechanically connected to the drive shaft. In particular, the present disclosure provides an improved method of calculating a model velocity such that an accurate vibration component can be extracted.

(14) As described in the foregoing, a vibration component of the drive shaft needs to be extracted to reduce vibration occurring from the drive shaft. Vibration reduction performance varies according to accuracy of extraction of the vibration component, and thus, it is important to accurately extract the vibration component.

(15) The vibration component of the drive shaft is extracted as below. An ideal model for the drive shaft, that is, a calculation model capable of calculating an ideal velocity (model velocity) of the drive shaft which does not consider vibration is designed. A model velocity corresponding to the drive shaft velocity excluding the vibration component is calculated using the calculation model, thereby extracting the vibration component using a difference between the calculated model velocity and an actual velocity corresponding to an actual drive shaft velocity.

(16) When the model velocity is accurately calculated as the ideal velocity of the drive shaft excluding the vibration component, an accurate vibration component may be extracted by calculating the difference between the two velocities. However, the calculated model velocity contains an error component when compared to the ideal velocity of the drive shaft excluding the vibration component.

(17) In order to eliminate the error, the difference between the model velocity and the actual velocity (actual drive shaft velocity) is obtained, and then error elimination control using a high pass filter (HPF), etc. is applied. An appropriate order of an error elimination controller needs to be determined based on a form (order) of the error component.

(18) In general, the order of the error elimination controller increases in proportion to the order of the error component, and occurrence of phase delay increases as the order of the error elimination controller increases. Thus, a different vibration component from actual vibration may be extracted.

(19) That is, when the order of the error component is minimized, an accurate vibration component may be extracted. The model velocity needs to be calculated to be a value as close to the ideal drive shaft velocity excluding the vibration component as possible in order to reduce the order of the error elimination controller including the filter, etc.

(20) In this regard, a disturbance torque applied to a vehicle may be observed using a disturbance observer when the model velocity is calculated, and the disturbance torque may be compensated when the model velocity is calculated, thereby minimizing the order of an error component contained in the model velocity. In this way, accuracy of extraction of the vibration component may be enhanced.

(21) The present disclosure may be applied to both the pure EV having the system configuration illustrated in FIG. 1 and the HEV having the system configuration illustrated in FIG. 2. FIGS. 1 and 2 will be referred to in description below with regard to the system configurations of the EV and the HEV.

(22) The drive shaft is a shaft from which a torque of the driving motor 13 is output in the systems of FIGS. 1 and 2. The drive shaft corresponds to an output shaft of a driving motor and an input shaft of a transmission. Here, a velocity of the drive shaft is the same as a velocity of the motor.

(23) Therefore, in description below, a drive shaft velocity may be replaced by a motor velocity, and a drive shaft model may be replaced by a motor model.

(24) The present disclosure will be described below with reference to FIGS. 3 and 4.

(25) FIG. 3 is a diagram illustrating a control system for reduction in drive shaft vibration of the eco-friendly vehicle according to an embodiment in the present disclosure, and FIG. 4 a flowchart illustrating a model velocity calculation process in a method for reducing drive shaft vibration of the eco-friendly vehicle according to an embodiment in the present disclosure.

(26) The control system for reduction in drive shaft vibration (that is, the anti-jerk control system) illustrated in FIG. 3 may be included in a MCU 17 illustrated in FIGS. 1 and 2.

(27) First, referring to FIG. 3, T.sub.1 is a torque required for a driving motor (first motor that drives the vehicle) (MG1) 13, and is a torque value subjected to torque compensation by a vibration reduction compensation torque T.sub.vib to reduce drive shaft vibration. The torque value may be a motor torque instruction value for driving the driving motor 13 to output the required torque.

(28) The motor torque instruction value T.sub.1 may be an instruction value delivered from a vehicle control unit (VCU) or a hybrid control unit (HCU) (not illustrated) to the MCU 17.

(29) In addition, in the present invention, the MCU 17 controls a torque output of the driving motor 13 using a final torque instruction value T.sub.1 calculated to compensate for the motor torque instruction value T.sub.1 by the vibration reduction compensation torque T.sub.vib such that vibration generated from the drive shaft can be reduced.

(30) The torque of T.sub.1 is output from the driving motor 13 according to a torque instruction which is compensated by the vibration reduction compensation torque T.sub.vib) and thus, T.sub.1 may be an actual output torque from the driving motor.

(31) In addition, T.sub.2 and T.sub.3 are torque values required for the engine 11 and the starter-generator (second motor directly connected to the engine) (MG2) 15, and may be torque output values of the engine 11 and the starter-generator 15 controlled according to instruction values thereof.

(32) However, in the EV, the engine 11 and the starter-generator 15 are not present, and thus, each of T.sub.2 and T.sub.3 becomes 0.

(33) In addition, T.sub.brake is a torque value required for a brake, and T.sub.load is a vehicle load torque generated by an angle of inclination of a road on which the vehicle is driven. T.sub.load is a torque value that can be calculated from the angle of inclination and a weight of the vehicle.

(34) Referring to FIG. 3, T.sub.2, T.sub.3, T.sub.brake, and T.sub.load are torques converted into torques of a shaft of the driving motor (MG1).

(35) In addition, referring to FIG. 3, T.sub.net is a torque applied to the drive shaft. This torque is a torque considering torques of T.sub.1, T.sub.2, T.sub.3, T.sub.brake, and T.sub.load, that is, output torques and load torques of torque sources (engine, driving motor, starter-generator, and brake) that apply torques to the drive shaft, and may be expressed as below.
T.sub.net=T.sub.1+T.sub.2+T.sub.3(T.sub.brake+T.sub.load)(1)

(36) Referring to FIG. 3, T.sub.acc is a torque obtained by adding an unknown disturbance torque d applied from the inside or outside of a driven vehicle to T.sub.net. T.sub.acc is an input torque of the drive shaft which actually accelerates/decelerates the drive shaft. T.sub.acc is delivered to the drive shaft such that a velocity co of the drive shaft is output.

(37) Referring to FIG. 3, G(s) represents a transfer function in an actual drive shaft 100.

(38) Here, except the disturbance torque d, all the above-described torques delivered to the drive shaft are known values. Thus, when an estimated value d of the disturbance torque as close to the actual disturbance torque d as possible can be obtained, a model velocity .sub.m may be accurately calculated.

(39) The actual drive shaft velocity is measured and obtained through a sensor, etc., and is used as a variable in a process of calculating the model velocity .sub.m and a process of extracting a vibration component .sub.vib.

(40) The vibration reduction compensation torque T.sub.vib for reduction of vibration generated from the drive shaft may be obtained through a process of calculating the model velocity .sub.m, a process of extracting the vibration component .sub.vib using an HPF or a band pass filter (BPF) from a difference between the calculated model velocity .sub.m and the measured actual velocity , and a process of calculating the vibration reduction compensation torque T.sub.vib based on the extracted vibration component .sub.vib and driving state information such as a driving mode, a gear stage, etc. of the vehicle.

(41) Here, the model velocity .sub.m is calculated by a drive shaft model velocity calculator 200 based on a torque component input to the drive shaft of the vehicle. As shown in FIG. 4, the model velocity .sub.m is obtained by calculating a drive shaft output requirement torque T.sub.net (S11), estimating a drive shaft input torque T.sub.acc input to the drive shaft using a measured actual velocity of the drive shaft (S12), estimating a disturbance torque d using the calculated drive shaft output requirement torque T.sub.net and the estimated drive shaft input torque T.sub.acc (S13), calculating a drive shaft model input torque T.sub.m considering a disturbance torque using the estimated disturbance torque d (514), and calculating the model velocity .sub.m using a drive shaft model 241 which receives the drive shaft model input torque T.sub.m as an input (S15).

(42) The drive shaft output requirement torque T.sub.net may be calculated by a drive shaft output requirement torque calculator 210 which receives torques T.sub.1, T.sub.2, T.sub.3, and T.sub.brake required for the driving motor 13, the engine 11, the starter-generator 15, and the brake (not illustrated) and a vehicle load torque T.sub.load as inputs (S11), and may be obtained by subtracting the vehicle load torque T.sub.load from torques required for torque sources of the vehicle which apply torques to the drive shaft.

(43) The torque sources of the vehicle may correspond to the driving motor 13, the engine 11, the starter-generator 15, and the brake. Here, the torque T.sub.brake required for the brake is a negative torque similarly to the load torque T.sub.load. Thus, the drive shaft output requirement torque T.sub.net may be calculated as in the following Equation.
T.sub.net=T.sub.1+T.sub.2+T.sub.3(T.sub.brake+T.sub.load)(2)
Here, the torque T.sub.1 required for the driving motor 13 may correspond to a torque instruction value for the driving motor, and all the torques T.sub.2, T.sub.3, and T.sub.brake required for the engine 11, the starter-generator 15, and the brake correspond to torque values converted into torques of the shaft of the driving motor (MG1).

(44) In this instance, the torques T.sub.2 and T.sub.3 required for the engine 11 and the starter-generator 15 may correspond to values obtained by converting values of an engine torque instruction and a starter-generator torque instruction into torques of the shaft of the driving motor (MG1), respectively, and the torque T.sub.brake required for the brake may correspond to a converted value of a brake torque which needs to be generated by a brake of a driving wheel.

(45) Next, the drive shaft input torque T.sub.acc using the measured actual velocity of the drive shaft (S12) is estimated by a drive shaft input torque estimator 220 of FIG. 3 which receives the actual velocity of the drive shaft as an input.

(46) As described in the foregoing, when the actual drive shaft input torque T.sub.acc is applied to the drive shaft, the drive shaft is rotated at a velocity of . When a transfer function in the drive shaft 100 is referred to as G(s), the actual velocity of the drive shaft and the drive shaft input torque T.sub.acc may be expressed by the following Equations.
=G(s)T.sub.acc(3)
T.sub.acc=/G(s)(4)

(47) In the above Equation (4), G(s) is the transfer function in the actual drive shaft 100. Thus, when a transfer function of an ideal drive shaft model which models the actual drive shaft 100, that is, the drive shaft model 241 designed to calculate the ideal model velocity .sub.m ignoring vibration is referred to as G.sub.m(s), the drive shaft input torque T.sub.acc may be estimated using G.sub.m(s) instead of G(s) in Equation 4.

(48) When the drive shaft is presumed to be a rigid body, and an estimated value of the drive shaft input torque is referred to as T.sub.acc, the estimated value of the drive shaft input torque may be expressed by the following Equation.
T.sub.acc/G.sub.m(s)=J.sub.ms, wherein G.sub.m(s)=1/J.sub.ms(5)

(49) In Equation (5), J.sub.m denotes a moment of inertia of the drive shaft which is a rigid body.

(50) Referring to Equation (5), the transfer function G.sub.m(s) for calculating the estimated value of the drive shaft input torque T.sub.acc applied to the drive shaft from the actual velocity of the drive shaft is a system in which an order of the numerator is greater than an order of the denominator. Thus, from a mathematical point of view, the transfer function G.sub.m(s) is calculated by differentiating the actual drive shaft velocity and multiplying by the moment of inertia, and thus is vulnerable to a noise component of the actual velocity of the drive shaft.

(51) In this regard, the denominator and the numerator may have the same order by using a filter Q. When a primary low pass filter LPF is used in the present embodiment, the drive shaft input torque T.sub.acc may be estimated as in the following Equation.

(52) $\begin{array}{cc}\begin{array}{c}{T}_{\mathrm{acc}}^{}=Q\left(s\right)/G_m\left(s\right)\\ =J_{m}s/\left(s+1\right),\end{array}& \left(6\right)\\ \mathrm{wherein}Q\left(s\right)=1/\left(s+1\right)& \end{array}$

(53) Q(s) is a transfer function of the LPF, and a time constant of the LPF Q is greater than a frequency of a vibration component such that a torque estimated by the vibration component can be excluded.

(54) The transfer function Q(s) of the filter is set such that an order of the numerator is less than or equal to an order of the denominator in Q(s)/G.sub.m(s) at all times, and the filter is additionally applied such that an order of the numerator is constantly less than or equal to an order of the denominator in the transfer function Q(s)/G.sub.m(s) for calculating the estimated value of the drive shaft input torque T.sub.acc from the actual velocity of the drive shaft , thereby ensuring robustness against the noise component.

(55) Next, the estimated value d of the disturbance torque is calculated by a disturbance torque estimator 230 which receives the drive shaft output requirement torque T.sub.net and the drive shaft input torque T.sub.acc as inputs, and may be calculated to be a difference between the drive shaft output requirement torque T.sub.net calculated by the drive shaft output requirement torque calculator 210 and the drive shaft input torque T.sub.acc calculated by the drive shaft input torque estimator 220.

(56) In this instance, the same filter Q as the filter used when the drive shaft input torque is estimated is applied to the drive shaft output requirement torque T.sub.net output from the drive shaft output requirement torque calculator 210 such that phase lag and magnitude change, which occur from the filter, similarly occur from the drive shaft output requirement torque, and the estimated value d of the disturbance torque is calculated using the drive shaft output requirement torque which is subjected to low-pass filtering through the filter Q.

(57) When the primary LPF is used, the disturbance torque d may be estimated as in the following Equation.
d=Q(s)T.sub.netT.sub.acc=T.sub.net/(s+1)J.sub.ms/(s+1)(7)

(58) When the disturbance torque d is estimated as described above, the drive shaft model input torque T.sub.m for calculation of the model velocity is calculated using the estimated disturbance torque d and the drive shaft output requirement torque T.sub.net calculated from the drive shaft output requirement torque calculator 210. In this instance, an equation for calculation of the drive shaft model input torque T.sub.m for calculation of the model velocity may be expressed by the following Equation.
T.sub.m=T.sub.m+d(8)

(59) As described above, when the drive shaft model input torque T.sub.m is calculated, the model velocity .sub.m is calculated by a velocity calculator 240 which receives the drive shaft model input torque T.sub.m as an input. In this instance, the model velocity .sub.m may be calculated by an equation below from the drive shaft model input torque T.sub.m using the transfer function G.sub.m(s) of the drive shaft model 241.
.sub.m=G.sub.m(s)T=T.sub.m/J.sub.ms(9)

(60) When the model velocity .sub.m is calculated as described above, a vibration component calculator 300 obtains a vibration component based on a deviation between the model velocity .sub.m and the actual velocity . In this instance, the vibration component .sub.vib can be calculated by applying an error elimination controller such as the HPF, etc. to the difference between the model velocity .sub.m and the actual velocity .

(61) Subsequently, after the vibration component .sub.vib is obtained, a vibration reduction compensation torque generator 400 calculates the vibration reduction compensation torque T.sub.vib based on the extracted vibration component .sub.vib and driving state information such as a driving mode, a gear stage, etc. of the vehicle.

(62) The vibration component .sub.vib using the model velocity .sub.m and the actual velocity and the vibration reduction compensation torque T.sub.vib based on the calculated vibration component .sub.vib and the driving state information of the vehicle may be calculated using the conventional method.

(63) When the vibration reduction compensation torque T.sub.vib is obtained as described above, the MCU 17 compensates for the torque instruction value T.sub.1 for the driving motor 13 by the vibration reduction compensation torque T.sub.vib, and controls a torque output of the driving motor 13 according to the compensated final torque instruction value T.sub.1.

(64) According to another embodiment, a plurality of model velocity calculation schemes may be applied, and one of the plurality of model velocity calculation schemes may be selected based on vehicle state information.

(65) A process of selecting a model velocity calculation scheme is additionally provided such that an advantage of each model velocity calculation scheme can be used according to a vehicle state. FIG. 5 illustrates a configuration of a control system for reduction in drive shaft vibration of the eco-friendly vehicle capable of selecting a model velocity calculation scheme as another embodiment.

(66) The embodiment of FIG. 5 is an embodiment allowing selection of a model velocity calculation scheme suitable for a current vehicle state from two model velocity calculation schemes, and one of the two model velocity calculation schemes is a model velocity calculation scheme described with reference to FIGS. 3 and 4, that is, a scheme in which a model velocity is calculated by the drive shaft model velocity calculator 200 of FIG. 3.

(67) In addition, the other of the two model velocity calculation schemes is a conventional model velocity calculation scheme.

(68) The model velocity calculation scheme according to the embodiment of FIGS. 3 and 4 is a scheme of calculating a model velocity .sub.m using a torque component, in which velocity delay does not occur when compared to an actual velocity of the drive shaft.

(69) In addition, the conventional model velocity calculation scheme is a scheme of calculating a model velocity .sub.m using a wheel velocity. The wheel velocity is a component generated when vibration of the drive shaft is attenuated. Thus, when the wheel velocity is used, an accurate model velocity can be calculated. In addition, the model velocity is calculated using only signal processing for the wheel velocity, and thus, a load factor is low.

(70) Therefore, the embodiment of FIG. 5 allows selection of a model velocity calculation scheme according to a vehicle state such that the above-described advantages of the respective model velocity calculation schemes can be used. In description below, the model velocity calculation scheme according to the embodiment of FIGS. 3 and 4 will be referred to as a torque-based model velocity calculation scheme, and the conventional model velocity calculation scheme will be referred to as a wheel velocity-based model velocity calculation scheme.

(71) In addition, a model velocity .sub.m obtained by the torque-based model velocity calculation scheme will be referred to as a torque-based model velocity, and a model velocity .sub.m obtained by the wheel velocity-based model velocity calculation scheme will be referred to as a wheel velocity-based model velocity.

(72) Referring to FIG. 5, reference numeral 200 denotes a torque-based model velocity calculator that calculates a model velocity according to the torque-based model velocity calculation scheme, and reference numeral 201 denotes a wheel velocity-based model velocity calculator that calculates a model velocity according to the wheel velocity-based model velocity calculation scheme.

(73) Reference numeral 202 denotes a model velocity selector that selects and changes a model velocity calculation scheme and outputs a model velocity calculated according to the selected model velocity calculation scheme in description below.

(74) First, vehicle state information may correspond to a gear stage of a transmission in a vehicle including the transmission, and a model velocity calculation scheme is selected according to the gear stage.

(75) In a low stage less than or equal to a set stage, a gear ratio of a drive shaft to a wheel is great, and thus a delay component of a wheel velocity may be great from a view point of the drive shaft. Therefore, the torque-based model velocity .sub.m is advantageously used. Conversely, in a high stage greater than the set stage, the wheel velocity-based model velocity .sub.m is advantageously used.

(76) FIG. 6 illustrates a process of selecting a model velocity calculation scheme and a model velocity according to a gear stage, and illustrates a process of changing a model velocity when a gear is shifted. First, in S21, the MCU 17 determines whether the gear is shifted based on a signal delivered from a transmission control unit (TCU) (not illustrated).

(77) When the gear is shifted, that is, when gear shifting is determined to be started, a model velocity calculation scheme to be used after gear shifting is determined based on a target gear stage of gear shifting in S22.

(78) When the model velocity calculation scheme according to the target gear stage is determined, whether a currently used model velocity calculation scheme needs to be changed is determined by comparing the determined model velocity calculation scheme after gear shifting with the currently used model velocity calculation scheme in S23.

(79) When the model velocity calculation scheme after gear shifting is different from the current model velocity calculation scheme, whether gear shifting is actually terminated is determined in S24. When actual termination of gear shifting is verified from a signal of a transmission controller, the current model velocity calculation scheme is changed to the new model velocity calculation scheme determined according to the target gear stage in S25.

(80) Here, when a model velocity calculation scheme is changed from the wheel velocity-based model velocity calculation scheme to the torque-based model velocity calculation scheme, an actual velocity of the drive shaft co at the time of change is used as a model velocity calculation initial value (initial model velocity .sub.m).

(81) When a model velocity calculation scheme is changed from the torque-based model velocity calculation scheme to the wheel velocity-based model velocity calculation scheme, the calculation scheme may be directly changed using a wheel velocity at the time of change.

(82) In other words, a model velocity .sub.m calculated by the changed wheel velocity-based model velocity calculation scheme is directly used as a model velocity initial value at the time of termination of gear shifting. Since a vibration reduction torque compensation function is deactivated at the time of gear shifting, the function is initialized at the time of termination of gear shifting, and then the function is activated again, an initial value of a model velocity at the time of changing the calculation scheme is not problematic.

(83) Next, at least one of a vehicle speed, a value of an accelerator pedal position sensor (APS), and a value of a brake pedal sensor (BPS) may correspond to vehicle state information that determines a model velocity calculation scheme as an example applicable to a vehicle which does not include a transmission.

(84) FIG. 7 illustrates a process of selecting a model velocity calculation scheme and a model velocity according to a vehicle speed, a value of the APS, and a value of the BPS. The MCU 17 receives vehicle state information such as a vehicle speed, a value of the APS, and a value of the BPS in S21, and determines whether a model velocity calculation scheme needs to be changed based on the received vehicle state information in S22.

(85) When vehicle speed information is used, the torque-based model velocity calculation scheme is used (torque-based model velocity is selected) under the condition of a low speed less than or equal to a predetermined reference vehicle speed, and the wheel velocity-based model velocity calculation scheme is used (wheel velocity-based model velocity is selected) under the condition of a high speed exceeding the reference vehicle speed.

(86) In addition, the APS is a sensor for detecting an operation state of an accelerator pedal, and the BPS is a sensor for detecting an operation state of a brake pedal. The vehicle is determined to be in a rapidly accelerated state when a rate of change (or the amount of change) of a value of the APS is greater than or equal to a first reference value, and is determined to be in a rapidly decelerated state when a rate of change (or the amount of change) of a value of the BPS is greater than or equal to a second reference value.

(87) When the vehicle is determined to be in the rapidly accelerated/decelerated state based on the value of the APS and the value of the BPS, the torque-based model velocity calculation scheme is used (torque-based model velocity is selected). Otherwise, the wheel velocity-based model velocity calculation scheme is used (wheel velocity-based model velocity is selected).

(88) When a model velocity calculation scheme is determined based on current vehicle state information such as a vehicle speed, a value of the APS, or a value of the BPS as described above, the determined model velocity calculation scheme is compared with a currently used model velocity calculation scheme to determine whether the model velocity calculation scheme needs to be changed.

(89) When the model velocity calculation scheme needs to be changed, the model velocity calculation scheme is changed to a new model velocity calculation scheme determined based on the current vehicle state information in S23.

(90) When the model velocity calculation scheme is changed from the wheel velocity-based model velocity calculation scheme to the torque-based model velocity calculation scheme in a process of changing a model velocity calculation scheme (changing a model velocity) as described above, a wheel velocity-based model velocity value at the time of change is used as an initial value of a torque-based model velocity.

(91) When the model velocity calculation scheme is changed from the torque-based model velocity calculation scheme to the wheel velocity-based model velocity calculation scheme, a value is determined according to a wheel velocity, and thus, an initial point cannot be set. Therefore, a discontinuous point of a model velocity may be generated at the time of change.

(92) A discontinuous point of a model velocity which is generated when the vibration reduction torque compensation function is continuously activated is a factor that degrades vibration component extraction performance, and thus, an initial point of a wheel velocity-based model velocity needs to be virtually set when a calculation scheme is changed. An equation for calculating a model velocity by setting a virtual initial point is given below.
.sub.m=.sub.m.sub.m.sub._.sub.init+.sub.m.sub._.sub.offset(10)

(93) .sub.m is a model velocity to which a virtual initial point is applied when a calculation scheme is changed, and the model velocity .sub.m of the virtual initial point is a model velocity at the time of change to the wheel velocity-based model velocity calculation scheme.

(94) .sub.m.sub._.sub.init is a difference value between the torque-based model velocity .sub.m and the wheel velocity-based model velocity .sub.m at the time of using an initial value (at the time of changing a calculation scheme), and .sub.m.sub._.sub.offset is a value which is set to vary from 0 to .sub.m.sub._.sub.init with a constant slope or pattern. For example, .sub.m.sub._.sub.offset may be calculated using a scheme illustrated in FIG. 8.

(95) Referring to FIG. 8, .sub.m.sub._.sub.offset is set to a determined velocity value for .sub.m.sub._.sub.init with respect to time.

(96) When an initial value .sub.m to be set is designated as a model velocity value of the wheel velocity-based calculation scheme at the time of changing a calculation scheme, it is possible to calculate a model velocity excluding a discontinuous point as in FIG. 9.

(97) As described in the foregoing, the present disclosure proposes an improved model velocity calculation scheme capable of extracting an accurate vibration component in a process of reducing vibration generated from a drive shaft of an eco-friendly vehicle.

(98) In the related art, a model velocity has been calculated using a wheel velocity or calculated based on a drive shaft input torque using a drive shaft model. When the model velocity is calculated using the wheel velocity, accuracy of the model velocity is high. However, an erroneous vibration component may be extracted since a model velocity is delayed when compared to an actual velocity at the time of acceleration/deceleration of a vehicle. When the model velocity is calculated based on the drive shaft input torque by designing the drive shaft model, the model velocity is not delayed at the time of acceleration/deceleration of the vehicle. However, accuracy of the model velocity is degraded due to inaccuracy of the model and a load torque, etc.

(99) On the other hand, according to the present disclosure, a drive shaft model is designed and a model velocity is calculated based on a drive shaft input torque, and thus the model velocity is not delayed at the time of acceleration/deceleration of a vehicle. At the same time, a disturbance torque delivered to the vehicle is observed through a disturbance observer and the model velocity is calculated based on the disturbance torque, and thus accuracy of the model velocity may be enhanced. Therefore, it is possible to more accurately extract a vibration component.

(100) According to a control method for reduction in drive shaft vibration of an eco-friendly vehicle of the present disclosure, a model velocity calculation scheme is improved such that an accurate vibration component can be extracted, thereby effectively reducing vibration generated from a drive shaft.

(101) In addition, one of a torque-based model velocity calculation scheme and a wheel velocity-based model velocity calculation scheme is selected and used when a model velocity is calculated, thereby having an effect of using advantages of both calculation schemes according to vehicle state.

(102) The invention has been described in detail with reference to exemplary embodiments thereof. However, it will be appreciated by those skilled in the art that changes may be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.