JP2011149312A - Control device of internal combustion engine - Google Patents

Control device of internal combustion engine Download PDF

Info

Publication number
JP2011149312A
JP2011149312A JP2010010225A JP2010010225A JP2011149312A JP 2011149312 A JP2011149312 A JP 2011149312A JP 2010010225 A JP2010010225 A JP 2010010225A JP 2010010225 A JP2010010225 A JP 2010010225A JP 2011149312 A JP2011149312 A JP 2011149312A
Authority
JP
Japan
Prior art keywords
surging
phase
standard deviation
calculated
cylinder
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2010010225A
Other languages
Japanese (ja)
Other versions
JP5351064B2 (en
Inventor
Hideji Takamiya
秀治 高宮
Futoshi Nishioka
太 西岡
Masayoshi Nishino
正嘉 西野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Honda Motor Co Ltd
Original Assignee
Honda Motor Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2010010225A priority Critical patent/JP5351064B2/en
Publication of JP2011149312A publication Critical patent/JP2011149312A/en
Application granted granted Critical
Publication of JP5351064B2 publication Critical patent/JP5351064B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/12Improving ICE efficiencies

Landscapes

  • Valve Device For Special Equipments (AREA)
  • Output Control And Ontrol Of Special Type Engine (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a technique for accurately preventing surging without increasing cost. <P>SOLUTION: A rotation speed parameter (OMG) in accordance with a rotation speed of the internal combustion engine is detected, a torque parameter of each cylinder (MFJUD) is calculated, and a standard deviation (σ<SB>-</SB>TRQ) of the torque parameter of each cylinder is calculated. An allowable acceleration map of a power transmission mechanism, and a standard deviation map defining the standard deviation corresponding to the allowable acceleration for each equivalent inertial weight of a vehicle are stored. The allowable acceleration of the corresponding surging is calculated by referring to the allowable acceleration map based on the detected gear ratio. The corresponding standard deviation is calculated as a threshold (σ<SB>-</SB>TRQ<SB>-</SB>SUR) by referring to the standard deviation map based on the calculated allowable acceleration and the equivalent inertial weight of the vehicle. It is determined that there is a possibility of occurrence of surging when the calculated standard deviation is larger than the threshold, and a variable valve train is controlled so that an overlap period of an air intake/exhaust valve is reduced. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

この発明は、内燃機関の制御装置に関し、より具体的には、内燃機関のサージングを防止しつつ、より最適な位相を実現するように可変動弁機構を制御する制御装置に関する。   The present invention relates to a control device for an internal combustion engine, and more particularly, to a control device that controls a variable valve mechanism so as to realize a more optimal phase while preventing surging of the internal combustion engine.

下記の特許文献1には、エンジンの燃焼不安定状態を検出するためのラフネスセンサを設け、該ラフネスの発生時において、少なくとも低回転低負荷時には、吸排気ポートのオーバーラップが小さくなる方向に、吸排気ポートの開閉タイミングの補正を行うことが開示されている。   In Patent Document 1 below, a roughness sensor for detecting an unstable combustion state of the engine is provided, and at the time of occurrence of the roughness, at least at a low rotation and low load, the overlap of the intake and exhaust ports becomes smaller. It is disclosed that the opening / closing timing of the intake / exhaust port is corrected.

また、下記の特許文献2には、所定のサンプリング母数分のトルク変動量から標準偏差または不偏分散を算出し、該算出した標準偏差または不偏分散に基づいて、燃焼状態を示すパラメータとしての燃焼ラフネス値を算出することが記載されている。   Further, in Patent Document 2 below, a standard deviation or unbiased variance is calculated from a torque fluctuation amount corresponding to a predetermined sampling parameter, and combustion as a parameter indicating a combustion state is calculated based on the calculated standard deviation or unbiased variance. It is described that the roughness value is calculated.

特開昭61−185601号公報JP-A 61-185601 特開2001−123879号公報JP 2001-123879 A

吸排気弁の開閉タイミングによっては、内燃機関においてサージングが発生するおそれがある。サージングは、エンジンの不整燃焼により、車体が前後に周期的に振動する現象をいう。サージングが発生すると、車両のドライバビリティが低下するため、サージングを防止することが望まれる。   Depending on the opening / closing timing of the intake / exhaust valve, surging may occur in the internal combustion engine. Surging is a phenomenon in which the vehicle body periodically oscillates back and forth due to irregular combustion of the engine. When surging occurs, the drivability of the vehicle decreases, so it is desirable to prevent surging.

上記の特許文献1では、燃焼の不安定状態を検出するセンサを別途設けており、このような付加的なセンサは、コスト高を招くおそれがある。また、センサ出力のばらつきによって、所望の精度を維持するのが困難となるおそれがある。   In the above-mentioned Patent Document 1, a sensor for detecting an unstable combustion state is provided separately, and such an additional sensor may increase the cost. In addition, it may be difficult to maintain desired accuracy due to variations in sensor output.

また、サージングが発生する条件は、車両の等価慣性重量や動力伝達機構の変速比に応じて変化する。上記の文献では、これらのパラメータについて何ら考慮されていない。   In addition, the conditions under which surging occurs change according to the equivalent inertia weight of the vehicle and the gear ratio of the power transmission mechanism. In the above-mentioned document, no consideration is given to these parameters.

また、最近の内燃機関においては、可変動弁機構を搭載し、吸気弁の開閉タイミング(位相)やそのリフト量を制御することが行われている。このような可変動弁機構においては、位相を進角することにより、ポンピングロスを低減し、よって燃費向上を図ることができる。他方、位相を進角しすぎると、吸気弁と排気弁が共に開いているオーバーラップ期間が大きくなり、サージングが発生するおそれがある。   In recent internal combustion engines, a variable valve mechanism is mounted to control the opening / closing timing (phase) of the intake valve and its lift amount. In such a variable valve mechanism, the pumping loss can be reduced by advancing the phase, thereby improving the fuel consumption. On the other hand, if the phase is advanced too much, the overlap period during which both the intake valve and the exhaust valve are open becomes large, and surging may occur.

したがって、この発明の一つの目的は、コスト高を招くことなく、より良好な精度でサージングを防止することができる手法を提供することである。また、この発明の他の目的は、サージングを防止しつつ、最適な進角量を持つ位相となるように、可変動弁機構を制御する手法を提供することである。   Accordingly, one object of the present invention is to provide a technique capable of preventing surging with better accuracy without incurring high costs. Another object of the present invention is to provide a method for controlling a variable valve mechanism so as to achieve a phase having an optimum advance amount while preventing surging.

この発明の一つの側面によると、車両に搭載され、吸気弁と排気弁のオーバーラップ期間を変更可能な可変動弁機構を有する内燃機関の制御装置が、内燃機関の回転速度に応じた回転速度パラメータ(OMG)を検出する回転速度パラメータ検出手段と、前記回転速度パラメータに基づいて、気筒別の出力トルクを表す気筒別トルクパラメータ(MFJUD)を算出する手段と、前記算出された気筒別トルクパラメータの標準偏差(σ_TRQ)を算出する手段と、動力伝達機構の変速比毎に、サージングによる振動の許容加速度を規定した許容加速度マップ(図10(a))、および車両の等価慣性重量毎に、該許容加速度に対応する気筒別トルクパラメータの標準偏差を規定した標準偏差マップ(図10(b))を記憶する記憶手段と、前記車両に搭載された動力伝達機構について検出された現在の変速比に基づいて前記許容加速度マップを参照し、対応するサージングの許容加速度を求める手段と、前記求めた許容加速度と前記車両の等価慣性重量に基づいて前記標準偏差マップを参照し、対応する標準偏差をしきい値(σ_TRQ_SUR)として求める手段と、前記算出された標準偏差が前記求めたしきい値より大きければ、サージングが発生する可能性がある状態(F_SURGE)と判定する手段と、サージングが発生する可能性がある状態と判定されたならば、前記オーバーラップ期間が小さくなるように、前記可変動弁機構を制御する制御手段と、を備える。   According to one aspect of the present invention, a control device for an internal combustion engine that is mounted on a vehicle and has a variable valve mechanism that can change an overlap period of an intake valve and an exhaust valve is provided with a rotational speed corresponding to the rotational speed of the internal combustion engine. Rotation speed parameter detection means for detecting a parameter (OMG), means for calculating a cylinder specific torque parameter (MFJUD) representing output torque for each cylinder based on the rotation speed parameter, and the calculated cylinder specific torque parameter A standard deviation (σ_TRQ) of the power transmission mechanism, a permissible acceleration map (FIG. 10A) defining permissible acceleration of vibration due to surging for each speed ratio of the power transmission mechanism, and for each equivalent inertia weight of the vehicle, Storage means for storing a standard deviation map (FIG. 10B) that defines the standard deviation of the cylinder-specific torque parameter corresponding to the allowable acceleration; The means for obtaining the permissible acceleration of the corresponding surging by referring to the permissible acceleration map based on the current speed ratio detected for the power transmission mechanism mounted on the vehicle, the permissible acceleration obtained and the equivalent inertia of the vehicle Means for obtaining the corresponding standard deviation as a threshold (σ_TRQ_SUR) by referring to the standard deviation map based on weight, and surging can occur if the calculated standard deviation is greater than the obtained threshold. And a control means for controlling the variable valve mechanism so that the overlap period is reduced if it is determined that surging is likely to occur. .

この発明によれば、内燃機関の回転速度に基づいて算出された気筒別トルクパラメータに基づいて、サージングが発生する可能性のある状態かどうかを判定するので、付加的なセンサを設けることなく、該サージングが発生する可能性のある状態かどうかを判定することができる。また、車両の等価慣性重量と動力伝達機構の変速比を考慮して該判定を行うので、良好な精度でサージングの発生を判定することができる。よって、サージングを防止しつつ、より最適な位相となるよう可変動弁機構を制御することができる。   According to the present invention, it is determined based on the cylinder-specific torque parameter calculated based on the rotational speed of the internal combustion engine whether surging is likely to occur, so without providing an additional sensor, It can be determined whether or not the surging is likely to occur. In addition, since the determination is performed in consideration of the equivalent inertia weight of the vehicle and the gear ratio of the power transmission mechanism, it is possible to determine the occurrence of surging with good accuracy. Therefore, it is possible to control the variable valve mechanism so as to achieve a more optimal phase while preventing surging.

この発明の一実施形態によると、内燃機関の運転状態に基づいて、吸気バルブの最適位相(CAINCMDM)を求める手段を備え、前記制御手段は、前記サージングが発生する可能性がある状態と判定されたならば、前記吸気バルブの位相を、前記最適位相から第1の所定量(CAIN_PRE)だけ遅角させる遅角制御を、前記可変動弁機構を介して行うことにより、前記オーバーラップ期間を小さくする。こうして、サージングの発生の可能性があると判定されたときには位相が遅角されるので、サージングを防止することができる。   According to one embodiment of the present invention, it is provided with means for obtaining the optimum phase (CAINCMDM) of the intake valve based on the operating state of the internal combustion engine, and the control means is determined to be in a state where the surging may occur. If this is the case, the overlap period can be reduced by performing, via the variable valve mechanism, retard control that retards the phase of the intake valve by a first predetermined amount (CAIN_PRE) from the optimum phase. To do. Thus, since it is determined that there is a possibility of occurrence of surging, the phase is retarded, so surging can be prevented.

この発明の一実施形態によると、前記制御手段は、さらに、前記サージングが発生する可能性がある状態と判定されたことに応じた前記遅角制御の後に、該サージングが発生する可能性のある状態が解消されたならば、前記吸気バルブの位相を、前記第1の所定量だけ遅角させた位相から、第2の所定量(ΔCAIN)ずつ前記最適位相に向けて進角させる進角制御を、前記可変動弁機構を介して行う。このように、一旦遅角された位相は、サージングの発生のおそれがないと判断される限り、最適位相に向けて進角されるので、サージングを防止しつつ、より最適な位相を実現して、ポンピングロスを低減し、燃費向上を図ることができる。   According to an embodiment of the present invention, the control unit may further generate the surging after the retardation control according to the determination that the surging is likely to occur. If the condition is resolved, the advance angle control advances the phase of the intake valve from the phase retarded by the first predetermined amount toward the optimum phase by a second predetermined amount (ΔCAIN). Is performed via the variable valve mechanism. In this way, the phase once retarded is advanced toward the optimum phase as long as it is determined that there is no risk of surging. Therefore, a more optimal phase can be realized while preventing surging. , Pumping loss can be reduced and fuel consumption can be improved.

この発明の一実施形態によると、前記制御手段は、さらに、前記進角制御によって進角された前記吸気バルブの位相に対し、前記最適位相で上限リミット処理する。このような上限リミット処理を行うことにより、進角された位相が、最適位相を超えることが防止される。   According to an embodiment of the present invention, the control means further performs upper limit processing at the optimum phase with respect to the phase of the intake valve advanced by the advance angle control. By performing such an upper limit process, the advanced phase is prevented from exceeding the optimum phase.

本発明のその他の特徴及び利点については、以下の詳細な説明から明らかである。   Other features and advantages of the present invention will be apparent from the detailed description that follows.

この発明の一実施例に従う、エンジンおよびその制御装置を概略的に示す図。The figure which shows schematically the engine and its control apparatus according to one Example of this invention. この発明の一実施例に従う、可変動弁機構による吸排気バルブの開閉特性を示す図。The figure which shows the opening / closing characteristic of the intake / exhaust valve by a variable valve mechanism according to one Example of this invention. この発明の一実施例に従う、吸気バルブの位相に対する、気筒別トルクパラメータの標準偏差および体感ポイントをプロットしたグラフを示す図。The figure which shows the graph which plotted the standard deviation and the sensory point of the torque parameter according to cylinder with respect to the phase of an intake valve according to one Example of this invention. この発明の一実施例に従う、吸気バルブの位相が遅角されている場合(a)と進角されている場合(b)における、気筒別トルクパラメータの時間的推移を示す図。The figure which shows the time transition of the cylinder-specific torque parameter in the case where the phase of the intake valve is retarded (a) and the case where it is advanced (b) according to one embodiment of the present invention. この発明の一実施例に従う、変速比毎に、体感ポイントに対応するサージングの振動による加速度をプロットしたグラフを示す図。The figure which shows the graph which plotted the acceleration by the vibration of the surging corresponding to a bodily sensation point for every gear ratio according to one Example of this invention. この発明の一実施例に従う、車両の等価慣性重量毎に、標準偏差に対応するサージングによる振動の加速度をプロットしたグラフを示す図。The figure which shows the graph which plotted the acceleration of the vibration by the surging corresponding to a standard deviation for every equivalent inertial weight of the vehicle according to one Example of this invention. この発明の一実施例に従う、制御装置のブロック図。The block diagram of the control apparatus according to one Example of this invention. この発明の一実施例に従う、目標位相の推移の一例を示す図。The figure which shows an example of transition of the target phase according to one Example of this invention. この発明の一実施例に従う、しきい値設定部のより詳細なブロック図。The more detailed block diagram of the threshold value setting part according to one Example of this invention. この発明の一実施例に従う、(a)変速比に対するサージングによる振動の許容加速度を規定したマップ、および(b)車両の等価慣性重量毎に、サージングによる振動の許容加速度に対応する標準偏差(しきい値)を規定したマップを示す図。According to one embodiment of the present invention, (a) a map defining allowable acceleration due to surging with respect to the gear ratio, and (b) standard deviation corresponding to the allowable acceleration due to surging for each equivalent inertia weight of the vehicle. The figure which shows the map which prescribed | regulated (threshold value). この発明の第1の形態に従う、トルクパラメータを算出する処理のフローチャート。The flowchart of the process which calculates a torque parameter according to the 1st form of this invention. この発明の一実施例に従う、トルクパラメータ算出のために用いられるパラメータの推移を示す図。The figure which shows transition of the parameter used for torque parameter calculation according to one Example of this invention. この発明の一実施例に従う、トルクパラメータ算出のために用いられるパラメータの推移を示す図。The figure which shows transition of the parameter used for torque parameter calculation according to one Example of this invention. この発明の一実施例に従う、トルクパラメータ算出のために用いられるパラメータの推移を示す図。The figure which shows transition of the parameter used for torque parameter calculation according to one Example of this invention. この発明の一実施例に従う、エンジンの往復運動部品の作動に起因する慣性力トルクの算出手法を説明するための図。The figure for demonstrating the calculation method of the inertial force torque resulting from the action | operation of the reciprocating motion part of an engine according to one Example of this invention. この発明の一実施例に従う、単一気筒慣性力トルク、合成慣性トルク、および慣性力回転速度を示す図。The figure which shows the single cylinder inertial force torque, the synthetic | combination inertial torque, and the inertial force rotational speed according to one Example of this invention. この発明の一実施例に従う、気筒間のトルクパラメータを示す図。The figure which shows the torque parameter between cylinders according to one Example of this invention. この発明の一実施例に従う、サージング防止のための位相補正プロセスのフローチャート。5 is a flowchart of a phase correction process for preventing surging according to an embodiment of the present invention. この発明の一実施例に従う、サージング防止のための位相補正プロセスのフローチャート。5 is a flowchart of a phase correction process for preventing surging according to an embodiment of the present invention. この発明の第1の形態の変形例に従うトルクパラメータを算出するプロセスのフローチャート。The flowchart of the process which calculates the torque parameter according to the modification of the 1st form of this invention. この発明の第2の形態に従うトルクパラメータを算出するプロセスのフローチャート。The flowchart of the process which calculates the torque parameter according to the 2nd form of this invention. この発明の第2の形態の変形例に従うトルクパラメータを算出するプロセスのフローチャート。The flowchart of the process which calculates the torque parameter according to the modification of the 2nd form of this invention. クランク角センサ出力に含まれる外乱の影響を説明するための図。The figure for demonstrating the influence of the disturbance contained in a crank angle sensor output. この発明の一実施例に従う、燃焼相関関数の例を示す図。The figure which shows the example of a combustion correlation function according to one Example of this invention. この発明の一実施例に従う、燃焼相関関数の他の例を示す図。The figure which shows the other example of a combustion correlation function according to one Example of this invention. この発明の一実施例に従う、トルクパラメータの実測値のばらつきを示す図。The figure which shows the dispersion | variation in the measured value of a torque parameter according to one Example of this invention. この発明の第3の形態に従うトルクパラメータを算出するプロセスのフローチャート。The flowchart of the process which calculates the torque parameter according to the 3rd form of this invention.

次に図面を参照してこの発明の実施の形態を説明する。図1は、この発明の実施形態に従う、内燃機関(以下、エンジンと呼ぶ)およびその制御装置の全体的な構成図である。   Next, an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is an overall configuration diagram of an internal combustion engine (hereinafter referred to as an engine) and its control device according to an embodiment of the present invention.

電子制御ユニット(以下、「ECU」)という)10は、入出力インターフェース、中央演算処理装置(CPU)、およびメモリを備えるコンピュータである。メモリには、車両の様々な制御を実現するためのコンピュータ・プログラムおよび該プログラムの実施に必要なデータを格納することができる。本発明に従う様々な制御のためのプログラム、および該プログラムの実行の際に用いるデータおよびマップは、メモリに格納されている。ECU10は、車両の各部から送られてくるデータを入出力インターフェースを介して受け取って演算を行い、制御信号を生成し、これを、該入出力インターフェースを介してエンジンの各部を制御するために送る。   An electronic control unit (hereinafter referred to as “ECU”) 10 is a computer including an input / output interface, a central processing unit (CPU), and a memory. The memory can store a computer program for realizing various controls of the vehicle and data necessary for executing the program. Programs for various controls according to the present invention, and data and maps used in executing the programs are stored in a memory. The ECU 10 receives data sent from each part of the vehicle via the input / output interface, performs an operation, generates a control signal, and sends this to control each part of the engine via the input / output interface. .

エンジン12は、この実施形態では6気筒を備えており、図には、そのうちの一つの気筒が概略的に示されている。エンジン12は、吸気バルブ14を介して吸気管16に連結され、排気バルブ18を介して排気管20に連結されている。ECU10からの制御信号に従って燃料を噴射する燃料噴射弁22が、吸気管16に設けられている。代替的に、燃料噴射弁22を、燃焼室24に設けてもよい。   The engine 12 includes six cylinders in this embodiment, and one of the cylinders is schematically shown in the drawing. The engine 12 is connected to an intake pipe 16 via an intake valve 14 and is connected to an exhaust pipe 20 via an exhaust valve 18. A fuel injection valve 22 that injects fuel in accordance with a control signal from the ECU 10 is provided in the intake pipe 16. Alternatively, the fuel injection valve 22 may be provided in the combustion chamber 24.

エンジン12は、吸気管16から吸入される空気と、燃料噴射弁22から噴射される燃料との混合気を、燃焼室24に吸入する。燃料室24には、ECU10からの点火時期信号に従って火花を飛ばす点火プラグ26が設けられている。点火プラグ26による火花により、混合気は燃焼する。この燃焼により混合気の体積は増大し、ピストン28を下方に押し下げる。ピストン28の往復運動は、クランク軸30の回転運動に変換される。4サイクルエンジンでは、エンジンのサイクルは、吸入、圧縮、燃焼、および排気行程からなる。ピストン28は、1サイクルにつき2往復する。   The engine 12 sucks an air-fuel mixture of air sucked from the intake pipe 16 and fuel injected from the fuel injection valve 22 into the combustion chamber 24. The fuel chamber 24 is provided with a spark plug 26 that discharges a spark in accordance with an ignition timing signal from the ECU 10. The air-fuel mixture is combusted by the spark from the spark plug 26. This combustion increases the volume of the air-fuel mixture and pushes the piston 28 downward. The reciprocating motion of the piston 28 is converted into the rotational motion of the crankshaft 30. In a four-cycle engine, the engine cycle consists of intake, compression, combustion, and exhaust strokes. The piston 28 reciprocates twice per cycle.

連続可変動弁機構31は、この実施形態では、可変リフト機構32および可変位相機構33から構成される。可変リフト機構32は、ECU10からの制御信号に従って、吸気バルブ14のリフト量を連続的に変更することができる機構である。可変リフト機構32は、任意の既知の手法により実現することができる。例えば、カム、リフト可変リンク、アッパーリンク、ロアリンクから構成され、ロアリンクの角度をアクチュエータなどで変更して、バルブの最大リフト量を制御する手法が提案されている(たとえば、特開2004−036560号を参照)。   In this embodiment, the continuously variable valve mechanism 31 includes a variable lift mechanism 32 and a variable phase mechanism 33. The variable lift mechanism 32 is a mechanism that can continuously change the lift amount of the intake valve 14 in accordance with a control signal from the ECU 10. The variable lift mechanism 32 can be realized by any known method. For example, a method has been proposed in which a cam, a lift variable link, an upper link, and a lower link are configured, and the maximum lift amount of the valve is controlled by changing the angle of the lower link with an actuator or the like (for example, Japanese Patent Laid-Open No. 2004-2004) 036560).

可変位相機構33は、ECU10からの制御信号に従って、吸気バルブ14の位相を連続的に変更する。可変位相機構は、任意の既知の手法により実現することができる。たとえば、電磁的に吸気バルブの位相を進角または遅角に制御する手法が提案されている(たとえば、特開2000―227033号を参照)。   The variable phase mechanism 33 continuously changes the phase of the intake valve 14 in accordance with a control signal from the ECU 10. The variable phase mechanism can be realized by any known method. For example, a method has been proposed in which the phase of the intake valve is electromagnetically controlled to advance or retard (see, for example, JP 2000-227033 A).

なお、代替的に、可変リフト機構32および可変位相機構33を一体的に構成してもよい。また、本願発明は、リフト量および位相を連続的に変更可能なこれら機構に限定されるわけではなく、リフト量および位相を段階状に変更可能な機構にも適用可能である。   Alternatively, the variable lift mechanism 32 and the variable phase mechanism 33 may be configured integrally. Further, the present invention is not limited to these mechanisms that can continuously change the lift amount and phase, but can also be applied to mechanisms that can change the lift amount and phase in stages.

ECU10には、エンジン12のクランク軸30の回転角度を検出するクランク角センサ33およびエンジン12の吸気バルブ14を駆動するカムが連結されたカム軸の回転角度を検出するカム角センサ34が接続されており、これらのセンサの検出値はECU10に供給される。クランク角センサ33は、所定のクランク角度毎に1パルス(CRKパルス)を発生し、該パルスにより、クランク軸30の回転角度位置を特定することができる。また、カム角センサ34は、エンジン12の特定の気筒の所定クランク角度位置でパルス(CYLパルス)と、各気筒の吸入行程開始時の上死点(TDC)でパルス(TDCパルス)を発生する。これらのパルスは、燃料噴射時期、点火時期等の各種の制御タイミングおよびエンジン回転数の検出に使用される。   The ECU 10 is connected to a crank angle sensor 33 that detects the rotation angle of the crankshaft 30 of the engine 12 and a cam angle sensor 34 that detects the rotation angle of the camshaft to which the cam that drives the intake valve 14 of the engine 12 is coupled. The detection values of these sensors are supplied to the ECU 10. The crank angle sensor 33 generates one pulse (CRK pulse) for each predetermined crank angle, and the rotational angle position of the crankshaft 30 can be specified by the pulse. The cam angle sensor 34 generates a pulse (CYL pulse) at a predetermined crank angle position of a specific cylinder of the engine 12 and a pulse (TDC pulse) at the top dead center (TDC) at the start of the intake stroke of each cylinder. . These pulses are used for various control timings such as fuel injection timing and ignition timing, and detection of engine speed.

なお、カム角センサ34により出力されるTDCパルスと、クランク角センサ33により出力されるCRKパルスとの相対関係から、吸気バルブ14のカム軸の実際の位相CAINが検出される。この実施例では、位相CAINは、最遅角をゼロとし、進角になるほど大きい値を持つ。   The actual phase CAIN of the cam shaft of the intake valve 14 is detected from the relative relationship between the TDC pulse output from the cam angle sensor 34 and the CRK pulse output from the crank angle sensor 33. In this embodiment, the phase CAIN has a value that increases as the angle of advance becomes zero with the most retarded angle being zero.

また、連続可変動弁機構31には、吸気バルブ14のリフト量LFTを検出するためのリフト量センサ35が設けられ、該センサは、ECU10に接続されている。リフト量センサ35は、任意の適切な手法により構成されることができる。たとえば、ポテンショメータにより、これらのセンサを構成することができる。該センサの検出値は、ECU10に送られる。   The continuously variable valve mechanism 31 is provided with a lift amount sensor 35 for detecting the lift amount LFT of the intake valve 14, and this sensor is connected to the ECU 10. The lift amount sensor 35 can be configured by any appropriate method. For example, these sensors can be configured by a potentiometer. The detection value of the sensor is sent to the ECU 10.

吸気管16内にはスロットル弁46が配置されている。スロットル弁46は、ECU10からの制御信号に応じてアクチュエータ(図示せず)によって駆動されるドライブバイワイヤ(drive by wire:DBW)式のスロットル弁である。   A throttle valve 46 is disposed in the intake pipe 16. The throttle valve 46 is a drive-by-wire (DBW) type throttle valve that is driven by an actuator (not shown) in accordance with a control signal from the ECU 10.

スロットル弁開度センサ48がスロットル弁46に設けられており、スロットル開度THに応じた信号をECU10に出力する。   A throttle valve opening sensor 48 is provided in the throttle valve 46 and outputs a signal corresponding to the throttle opening TH to the ECU 10.

吸気管16のスロットル弁46の上流側に、エアフローメータ50が設置されている。エアフローメータ50は、吸入空気量を示す電気信号をECU10に出力する。   An air flow meter 50 is installed upstream of the throttle valve 46 in the intake pipe 16. The air flow meter 50 outputs an electrical signal indicating the amount of intake air to the ECU 10.

吸気管16のスロットル弁46の下流には吸気管内圧力センサ52および吸気温センサ54が備えられ、それぞれ、吸気管内絶対圧PBおよび吸気温度TAを示す電気信号をECU10に出力する。また、大気圧センサ55がエンジン外部の任意の位置に設置されており、大気圧PAを示す電気信号をECU10に出力する。さらに、エンジン12の水温TWを検出するエンジン水温センサ56が備えられ、エンジンの水温を示す電気信号をECU10に出力する。   An intake pipe pressure sensor 52 and an intake air temperature sensor 54 are provided downstream of the throttle valve 46 of the intake pipe 16 and output electric signals indicating the intake pipe absolute pressure PB and the intake air temperature TA to the ECU 10, respectively. An atmospheric pressure sensor 55 is installed at an arbitrary position outside the engine, and outputs an electric signal indicating the atmospheric pressure PA to the ECU 10. Further, an engine water temperature sensor 56 for detecting the water temperature TW of the engine 12 is provided, and an electric signal indicating the engine water temperature is output to the ECU 10.

排気管20の触媒58の上流側にはLAF(linear air-fuel)センサ60が設置されている。LAFセンサ60は、リーンからリッチにわたる広範囲において排ガス中の酸素濃度に比例する信号をECU10に出力する。   A LAF (linear air-fuel) sensor 60 is installed upstream of the catalyst 58 in the exhaust pipe 20. The LAF sensor 60 outputs a signal proportional to the oxygen concentration in the exhaust gas to the ECU 10 over a wide range from lean to rich.

エンジン12には、エンジンの回転を車輪に伝達するための動力伝達装置(図示せず)が接続されており、該装置は、変速動作を行う変速機と、該回転を伝達する各種のシャフト、および左右輪の回転速度を調整する作動歯車装置など(いずれも図示せず)を備えている。変速機は、任意の種類のものでよく、たとえば、手動変速機(MT)、自動MT、自動変速機(AT)、および無断連続変速機(CVT)のいずれかであることができる。   The engine 12 is connected to a power transmission device (not shown) for transmitting the rotation of the engine to the wheels. The device includes a transmission that performs a shifting operation, various shafts that transmit the rotation, And an operating gear device (not shown) for adjusting the rotational speeds of the left and right wheels. The transmission can be of any type and can be, for example, any of a manual transmission (MT), an automatic MT, an automatic transmission (AT), and a continuously variable transmission (CVT).

変速機は、シフトレバーのポジションが、P(パーキング)およびN(ニュートラル)にシフトされたときには、クラッチを解放することにより、エンジン12の回転の車輪への伝達を遮断する。シフトレバーのポジションがR(後退)およびD(ドライブ)にシフトされたときには、所望の変速比を実現するようクラッチを締結して、エンジンの回転を車輪に伝達する。たとえば、手動変速機や自動MTの場合には、変速機とエンジンの間に該クラッチが設けられており、該クラッチの締結および解放を介して、エンジンの回転の車輪への伝達を制御することができる。自動変速機で典型的に採用されるトルクコンバータ方式の場合には、副変速機内のクラッチの締結および解放を介して、エンジンの回転の車輪への伝達を制御することができる。これらのクラッチの締結および解放を検出するためのセンサすなわちクラッチスイッチ61が設けられており、該クラッチスイッチ61の出力信号は、ECU10に供給される。また、変速機には、変速機によって現在選択されている変速比を検出するための変速比センサ63が設けられており、該変速比センサ63の出力信号は、ECU10に供給される。   When the position of the shift lever is shifted to P (parking) and N (neutral), the transmission disengages the transmission of the rotation of the engine 12 to the wheels by releasing the clutch. When the position of the shift lever is shifted to R (reverse) and D (drive), the clutch is engaged to achieve a desired gear ratio, and the rotation of the engine is transmitted to the wheels. For example, in the case of a manual transmission or an automatic MT, the clutch is provided between the transmission and the engine, and the transmission of the engine rotation to the wheels is controlled through engagement and release of the clutch. Can do. In the case of a torque converter system typically employed in an automatic transmission, transmission of engine rotation to wheels can be controlled through engagement and disengagement of a clutch in the auxiliary transmission. A sensor for detecting engagement and disengagement of these clutches, that is, a clutch switch 61 is provided, and an output signal of the clutch switch 61 is supplied to the ECU 10. The transmission is provided with a transmission ratio sensor 63 for detecting a transmission ratio currently selected by the transmission, and an output signal of the transmission ratio sensor 63 is supplied to the ECU 10.

図2を参照すると、この発明の一実施形態に従う、吸排気バルブの開閉特性が示されている。縦軸はリフト量を示し、横軸はクランク角度を表す。連続可変動弁機構31により、吸気バルブ14の位相を連続的に変化させることができると共に、リフト量(吸気バルブが開かれる最大リフト量)を連続的に変化させることができる。   Referring to FIG. 2, the opening / closing characteristics of an intake / exhaust valve according to an embodiment of the present invention are shown. The vertical axis represents the lift amount, and the horizontal axis represents the crank angle. The continuously variable valve mechanism 31 can continuously change the phase of the intake valve 14 and can continuously change the lift amount (the maximum lift amount by which the intake valve is opened).

排気バルブ18のクランク角度に対するリフト量(以下、作動特性と呼ぶ)が、符合EX1によって示されている。実線で示されるIN11〜IN15は、吸気バルブ14の基準位相における作動特性を示している。この実施例では、該基準位相は、TDC(圧縮上死点)に対して最遅角の所に設定されており、該基準位相をゼロ度とする。IN15からIN11に向けて、リフト量は増大しており、エンジンの運転状態に応じて、IN11〜IN15の間で作動特性が切り換えられる。点線で示されるIN21〜IN25は、作動特性IN11〜IN15がそれぞれTDCに対して最進角された状態の作動特性を示している。   The lift amount with respect to the crank angle of the exhaust valve 18 (hereinafter referred to as the operating characteristic) is indicated by the symbol EX1. IN11 to IN15 indicated by solid lines indicate operating characteristics of the intake valve 14 at the reference phase. In this embodiment, the reference phase is set at the most retarded angle with respect to TDC (compression top dead center), and the reference phase is set to zero degrees. The lift amount increases from IN15 to IN11, and the operating characteristics are switched between IN11 to IN15 according to the operating state of the engine. IN21 to IN25 indicated by dotted lines indicate the operating characteristics in a state where the operating characteristics IN11 to IN15 are respectively advanced to the TDC.

図に示されるように、この可変動弁機構においては、位相を進角にするほど、吸気バルブの閉タイミングが早くなるので、ポンピングロスを低減することができ、よって燃費向上を図ることができる。他方、位相を進角にするほど、吸気バルブと排気バルブが共に開いているオーバーラップ期間が長くなる。そのため、排気の吹き返し等が生じやすくなり、サージングを誘起するおそれがある。   As shown in the figure, in this variable valve mechanism, the more the phase is advanced, the earlier the closing timing of the intake valve, so that the pumping loss can be reduced, and the fuel efficiency can be improved. . On the other hand, as the phase is advanced, the overlap period in which both the intake valve and the exhaust valve are open becomes longer. For this reason, exhaust blow-back is likely to occur, and surging may be induced.

したがって、本願発明では、サージングを防止することが可能な最適な進角量に、位相を制御する手法を提案する。まず、本願発明の当該手法に至った技術的根拠を、図3〜図6を参照して説明する。   Therefore, the present invention proposes a method of controlling the phase to the optimum advance amount that can prevent surging. First, the technical basis that led to the method of the present invention will be described with reference to FIGS.

図3は、車両の商品性を考察した実験結果の一例を示す。図の横軸は、可変動弁機構31の前述した吸気バルブ14の位相を示し、ゼロ度(最遅角)を基準として、右に行くほど進角量が大きいことを表している。符号71は、該位相に対する、後述する気筒毎の出力トルクを表す気筒別トルクパラメータの標準偏差の推移を示す。標準偏差の値が小さいほど、該気筒の出力トルクの変動量(ばらつき)が小さいことを示す。符号73は、該位相に対する体感ポイントの推移を示す。体感ポイントは、車両の商品性を表すポイント(点数)であり、専門の評価者が車両を運転した際に、体感によって、サージングをどの程度感じるかを表したものである。体感ポイントの値が高いほど、サージングが抑制されており、よって商品性が高いことを示す。   FIG. 3 shows an example of an experimental result considering the merchantability of the vehicle. The horizontal axis of the figure indicates the phase of the above-described intake valve 14 of the variable valve mechanism 31 and indicates that the advance amount is larger toward the right with reference to zero degree (most retarded angle). Reference numeral 71 indicates the transition of the standard deviation of the cylinder-specific torque parameter representing the output torque of each cylinder, which will be described later, with respect to the phase. The smaller the standard deviation value, the smaller the variation amount (variation) of the output torque of the cylinder. Reference numeral 73 indicates the transition of the sensation point with respect to the phase. The sensation point is a point (score) representing the merchantability of the vehicle, and represents how much surging is felt by the sensation when a specialized evaluator drives the vehicle. The higher the value of the sensation point, the more the surging is suppressed, and the higher the commercial value.

図に示すように、位相が遅角側にあるときには、体感ポイントは、ほぼ一定の値で推移しており、標準偏差も、安定した挙動を示している。しかしながら、位相の進角量が大きくなると、領域75で示されるように、体感ポイントが減少し始める。それに伴い、標準偏差の値が急激に上昇し始める。   As shown in the figure, when the phase is on the retard side, the sensation point changes at a substantially constant value, and the standard deviation also shows a stable behavior. However, as the phase advance amount increases, as shown by the region 75, the sensation points start to decrease. Along with this, the value of the standard deviation starts to rise rapidly.

すなわち、位相の進角量が大きくなりすぎると、体感によって認識可能なサージングが発生することがわかる。また、体感ポイントの減少に伴い、気筒別トルクパラメータの標準偏差の値が急激に上昇し始めることがわかる。   That is, it can be understood that if the phase advance amount becomes too large, surging that can be recognized by the bodily sensation occurs. It can also be seen that the standard deviation value of the cylinder-specific torque parameter starts to increase rapidly as the sensation points decrease.

また、図4を参照すると、所定の気筒について、(a)は、位相が基準位相(前述したように、ゼロ度)にあるときの気筒別トルクパラメータの値の時間的推移を示し、(b)は、位相が、基準位相に対して所定量だけ進角したときの、気筒別トルクパラメータの値の時間的推移を示す。位相が進角するほど、気筒別トルクパラメータの変動量が大きくなっているのがわかる。   Referring to FIG. 4, for a given cylinder, (a) shows the temporal transition of the value of the cylinder-specific torque parameter when the phase is at the reference phase (as described above, zero degree), and (b ) Shows the temporal transition of the value of the torque parameter for each cylinder when the phase is advanced by a predetermined amount with respect to the reference phase. It can be seen that the amount of variation in the torque parameter for each cylinder increases as the phase advances.

このように、サージングの観点からの商品性と、気筒別トルクパラメータの標準偏差との間には相関があり、位相の進角量が大きすぎると、該標準偏差が増大され、それに伴い、サージングの体感が増すことが判明した。したがって、気筒別トルクパラメータの標準偏差に基づいて、所定の商品性を実現するのに必要な体感ポイント値を維持することが困難なサージングの発生のおそれがある状態かどうかを判断することができる。このようなサージングの発生のおそれがある状態と判断された場合には、吸気バルブの位相を遅角させることにより、すなわちオーバーラップ期間が小さくなるように可変動弁機構を制御することにより、サージングを防止することができる。   Thus, there is a correlation between the merchantability from the viewpoint of surging and the standard deviation of the torque parameter for each cylinder. If the phase advance is too large, the standard deviation is increased. It became clear that the bodily sensation increased. Therefore, based on the standard deviation of the cylinder-specific torque parameter, it is possible to determine whether or not there is a possibility of surging that makes it difficult to maintain the sensation point value necessary to achieve a predetermined commercial value. . When it is determined that such a surging may occur, the surging is delayed by retarding the phase of the intake valve, that is, by controlling the variable valve mechanism so that the overlap period is reduced. Can be prevented.

この具体的な方法を説明すると、サージングの大きさと体感ポイントとの間の相関は予め決まっており、これが、図5に示されている。該相関は、変速比に応じて変化し、図には、一例として、変速比が3速、4速および5速のそれぞれに対応する場合についての相関が示されている。   Explaining this specific method, the correlation between the magnitude of surging and the sensation points is predetermined, and this is shown in FIG. The correlation changes in accordance with the gear ratio, and in the drawing, as an example, the correlation is shown for the case where the gear ratio corresponds to the third speed, the fourth speed, and the fifth speed.

この実施形態では、サージングによる振動の大きさを、該振動の加速度で表している。当然ながら、サージングの振動加速度が小さいほど、体感ポイント値は大きくなり、よって商品性が高くなる。所定の商品性を実現するのに必要な体感ポイントの下限値は予め決まっている。したがって、図に示すように、現在の変速比に基づき、該所望の体感ポイントの下限値P1に対して許容されるサージングの振動加速度の上限値A1が決定される。図の例では、変速比が3速に対応する場合の該上限値A1が決定されている。   In this embodiment, the magnitude of vibration due to surging is represented by the acceleration of the vibration. Naturally, the smaller the vibration acceleration of surging, the greater the sensation point value, and the higher the commercial value. The lower limit value of the sensation points necessary for realizing predetermined merchantability is determined in advance. Therefore, as shown in the figure, based on the current gear ratio, the upper limit value A1 of the surging vibration acceleration allowed for the lower limit value P1 of the desired sensation point is determined. In the example of the figure, the upper limit value A1 when the gear ratio corresponds to the third speed is determined.

該許容される振動加速度A1は、変速比が小さくなるほど(すなわち、高速用の変速比になるほど)、高くなる。変速比が小さくなるほど、車速が高くなるので、体感するサージングの振動の許容幅が大きくなるためである。   The allowable vibration acceleration A1 increases as the gear ratio decreases (that is, as the gear ratio for high speed increases). This is because the vehicle speed increases as the gear ratio decreases, and the allowable width of surging vibrations to be experienced increases.

また、図6は、車両の等価慣性重量ごとに、気筒別トルクパラメータの標準偏差と、サージングの振動加速度との相関を、実験やシミュレーション等によって調べた結果を示す。ここで、等価慣性重量(kg)は、エンジンや駆動系などの回転部分の回転慣性力を、車両が加速するときの抵抗として見かけの車両重量増加分として加える重量をいう。等価慣性重量は、車両毎に予め決まっている。   FIG. 6 shows the results of examining the correlation between the standard deviation of the cylinder-specific torque parameter and the surging vibration acceleration for each equivalent inertia weight of the vehicle through experiments, simulations, and the like. Here, the equivalent inertia weight (kg) refers to a weight that is added as an apparent vehicle weight increase as a resistance when the vehicle accelerates the rotational inertia force of a rotating part such as an engine or a drive system. The equivalent inertia weight is determined in advance for each vehicle.

標準偏差値が大きいほど、すなわち気筒の出力トルクのばらつきが大きいほど、サージングの振動加速度は大きくなる。また、等価慣性重量が大きいほど、標準偏差に対応する振動加速度は小さくなるが、これは、慣性力が高いほど、振動が起こりにくくなるためである。   The greater the standard deviation value, that is, the greater the variation in the output torque of the cylinder, the greater the surging vibration acceleration. Further, the greater the equivalent inertia weight, the smaller the vibration acceleration corresponding to the standard deviation. This is because the higher the inertial force, the less likely vibration will occur.

図6の相関関係を参照することにより、図5により決定された許容される振動加速度の上限値A1に対応する標準偏差を、車両の等価慣性重量に基づいて求め、こうして求めた標準偏差をしきい値に設定する。図には、符号77で示される等価慣性重量に基づいて求めたしきい値が示されている。該しきい値を超えるような値の標準偏差が実際に算出されたならば、商品性の観点から許容されない大きさのサージングが発生するおそれのある状態と判定することができる。   By referring to the correlation in FIG. 6, the standard deviation corresponding to the upper limit value A1 of the allowable vibration acceleration determined in FIG. 5 is obtained based on the equivalent inertia weight of the vehicle, and the standard deviation thus obtained is calculated. Set to threshold. In the figure, threshold values obtained based on the equivalent inertia weight indicated by reference numeral 77 are shown. If the standard deviation of the value that exceeds the threshold value is actually calculated, it can be determined that there is a possibility of surging having a size that is unacceptable from the viewpoint of merchantability.

このように、本願発明によれば、変速比および等価慣性重量を考慮して、商品性の観点から許容されないサージングが発生するおそれのある状態かどうかを判定するので、該判定の精度を向上させることができる。   As described above, according to the present invention, it is determined whether there is a possibility that surging that is unacceptable from the viewpoint of merchantability occurs in consideration of the gear ratio and the equivalent inertia weight, so that the accuracy of the determination is improved. be able to.

サージングが発生するおそれのある状態と判定されたならば、位相を遅角化し、オーバーラップ期間を小さくする。図3および図4を参照して述べたように、位相の進角量を減らすことにより、オーバーラップ期間が短くなり、標準偏差の上昇を抑制することができる。よって、サージングの振動加速度を、許容される振動加速度の上限値以下に制限することができる。結果として、所望の体感ポイント値を維持することができ、よって所望の商品性を確保することができる。   If it is determined that surging is likely to occur, the phase is retarded and the overlap period is reduced. As described with reference to FIGS. 3 and 4, by reducing the amount of phase advance, the overlap period is shortened, and an increase in standard deviation can be suppressed. Therefore, the vibration acceleration of surging can be limited to the upper limit value of the allowable vibration acceleration. As a result, a desired sensation point value can be maintained, and thus desired merchantability can be ensured.

図7は、上述した技術的根拠に基づく、本願発明の一実施例に従うエンジンの制御装置100の機能ブロック図である。これらの機能ブロックは、ECU10において実現される。なお、該制御装置100は、気筒毎に、サージング防止のための位相制御を行う。   FIG. 7 is a functional block diagram of engine control apparatus 100 according to one embodiment of the present invention based on the above-described technical basis. These functional blocks are realized in the ECU 10. The control device 100 performs phase control for preventing surging for each cylinder.

気筒別トルクパラメータ算出部111は、気筒毎に、クランク角センサ33の出力から得られるエンジンの回転速度に基づいて、該気筒の出力トルクを表すトルクパラメータMFJUDを算出する。この詳細なプロセスは、後述される。   The cylinder-specific torque parameter calculation unit 111 calculates, for each cylinder, a torque parameter MFJUD representing the output torque of the cylinder based on the engine speed obtained from the output of the crank angle sensor 33. This detailed process is described below.

標準偏差算出部113は、気筒毎に、気筒別トルクパラメータMFJUDの標準偏差σ_TRQを算出する。所定数(たとえば、100個)の制御周期にわたって、気筒別トルクパラメータ値MFJUDi(i=1〜100)を算出し、これらの気筒別トルクパラメータ値MFJUDiに基づいて、以下の式(1)により、標準偏差σ_TRQを算出することができる。ここで、nは、上記の所定数に対応する。また、MFJUDバーは、該所定数の制御周期におけるトルクパラメータMFJUDの平均値を表す。

Figure 2011149312
The standard deviation calculator 113 calculates the standard deviation σ_TRQ of the cylinder specific torque parameter MFJUD for each cylinder. The cylinder-specific torque parameter value MFJUDi (i = 1 to 100) is calculated over a predetermined number (for example, 100) of control cycles, and based on these cylinder-specific torque parameter values MFJUDi, A standard deviation σ_TRQ can be calculated. Here, n corresponds to the predetermined number. The MFJUD bar represents the average value of the torque parameter MFJUD in the predetermined number of control cycles.
Figure 2011149312

しきい値設定部115は、サージングが発生する可能性がある状態かどうかを判定するためのしきい値σ_TUR_SURを設定する。しきい値σ_TUR_SURは、図5および図6を参照して説明したように、変速比センサ63によって検出された現在の変速比と、車両の等価慣性重量とに基づいて設定される。しきい値設定部115の詳細なブロック図は、さらに後述される。   The threshold value setting unit 115 sets a threshold value σ_TUR_SUR for determining whether or not surging is likely to occur. The threshold value σ_TUR_SUR is set based on the current speed ratio detected by the speed ratio sensor 63 and the equivalent inertia weight of the vehicle, as described with reference to FIGS. A detailed block diagram of the threshold setting unit 115 will be further described later.

サージ判定部117は、気筒別トルクパラメータの標準偏差σと、設定されたしきい値σ_TUR_SURとを比較する。標準偏差σが、しきい値σ_TUR_SURより大きければ、許容されない大きさのサージング(前述したように、より正確には、商品性についての所定のレベルを超えたサージング)が発生する可能性のある状態と判定し、サージ判定フラグF_SURGEに値1を設定する。標準偏差σが、しきい値σ_TUR_SUR以下であれば、該許容されないサージングが発生する可能性はない状態と判定し、サージ判定フラグF_SURGEにゼロを設定する。   The surge determination unit 117 compares the standard deviation σ of the cylinder-specific torque parameter with the set threshold value σ_TUR_SUR. If the standard deviation σ is larger than the threshold value σ_TUR_SUR, an unacceptable magnitude surging (more precisely, surging beyond a predetermined level of merchantability as described above) may occur. And the value 1 is set to the surge determination flag F_SURGE. If the standard deviation σ is equal to or smaller than the threshold σ_TUR_SUR, it is determined that there is no possibility that the unacceptable surging will occur, and the surge determination flag F_SURGE is set to zero.

最適位相決定部121は、現在のエンジンの運転状態に基づいて、所定のマップを参照することにより、現在のエンジンの運転状態に最適な位相CAINCMDMを決定する。エンジンの運転状態として、任意の適切なパラメータを用いることができる。たとえば、現在のエンジン回転数および吸気バルブのリフト量に基づいて所定のマップを参照することにより、対応する位相を求めることで、最適位相CAINCMDMを決定することができる。なお、該マップ上で、最適位相CAINCMDMは、対応するエンジンの運転状態下でポンピングロスをなるべく低減させることができる進角量を持つように設定されている。   The optimum phase determination unit 121 determines the optimum phase CAINCMDM for the current engine operating state by referring to a predetermined map based on the current engine operating state. Any appropriate parameter can be used as the operating state of the engine. For example, the optimum phase CAINCMDM can be determined by obtaining a corresponding phase by referring to a predetermined map based on the current engine speed and the lift amount of the intake valve. On the map, the optimum phase CAINCMDM is set to have an advance amount that can reduce the pumping loss as much as possible under the operation state of the corresponding engine.

遅角量切り換え部123は、サージ判定フラグF_SURGEの値に従って、出力を切り換える。サージ判定フラグF_SURGEに値1が設定されているときには、遅角量ΔCAIN_SURに所定値CAIN_PREを設定して出力する。サージ判定フラグF_SURGEにゼロが設定されているときには、遅角量ΔCAIN_SURに、該所定値CAIN_PREから所定量ΔCAINを減算した値を設定する。ここで、ΔCAINの大きさは、CAIN_PREの大きさより小さい。   The retard amount switching unit 123 switches the output according to the value of the surge determination flag F_SURGE. When the value 1 is set in the surge determination flag F_SURGE, a predetermined value CAIN_PRE is set in the retard amount ΔCAIN_SUR and output. When the surge determination flag F_SURGE is set to zero, a value obtained by subtracting the predetermined amount ΔCAIN from the predetermined value CAIN_PRE is set to the retardation amount ΔCAIN_SUR. Here, the magnitude of ΔCAIN is smaller than the magnitude of CAIN_PRE.

位相補正部125は、最適位相CAINCMDMから、遅角量ΔCAIN_SURを減算することにより、補正済み位相CAINCMDXを算出する(CAINCMDX=CAINCMDM−ΔCAIN_SUR)。   The phase correction unit 125 calculates the corrected phase CAINCMDX by subtracting the retardation amount ΔCAIN_SUR from the optimum phase CAINCMDM (CAINCMDX = CAINCMDM−ΔCAIN_SUR).

目標位相算出部127は、補正済み位相CAINCMDXを、最適位相CAINCMDMで上限リミット処理し、目標位相CAINCMDを算出する。すなわち、補正済み位相CAINCMDXが最適位相CAINCMDMを超えたならば、該最適位相CAINCMDMを目標位相CAINCMDに設定し、補正済み位相CAINCMDXが最適位相CAINCMDM以下ならば、該補正済み位相CAINCMDXを目標位相CAINCMDに設定する。こうして、目標位相CAINCMDは、最適位相CAINCMDMを超えないよう制御される。吸気バルブ14の位相は、該目標位相CAINCMDに達するように、可変位相機構33を介して制御される。   The target phase calculation unit 127 performs upper limit processing on the corrected phase CAINCMDX with the optimum phase CAINCMDM, and calculates the target phase CAINCMD. That is, if the corrected phase CAINCMDX exceeds the optimal phase CAINCMDM, the optimal phase CAINCMDM is set to the target phase CAINCMD. If the corrected phase CAINCMDX is less than or equal to the optimal phase CAINCMDM, the corrected phase CAINCMDX is set to the target phase CAINCMD. Set. Thus, the target phase CAINCMD is controlled so as not to exceed the optimum phase CAINCMDM. The phase of the intake valve 14 is controlled via the variable phase mechanism 33 so as to reach the target phase CAINCMD.

ここで図8を参照すると、目標位相CAICMDの推移の一例が概略的に示されている。横軸は、基準位相(ゼロ度)に対する位相値を示す。この例では、最適位相CAINCMDMは一定であると仮定する。   Referring now to FIG. 8, an example of the transition of the target phase CAICMD is schematically shown. The horizontal axis indicates the phase value with respect to the reference phase (zero degree). In this example, it is assumed that the optimal phase CAINCMDM is constant.

制御周期nにおいて、サージングが発生するおそれのある状態が発生したために、サージ判定フラグF_SURGEに値1が設定され、それに応じて、最適位相CAINCMDMから、所定値CAIN_PREだけ遅角された位相が、目標位相となる。この遅角動作により、次の制御周期(n+1)では、サージングの発生のおそれがなくなり、サージ判定フラグF_SURGEにゼロが設定される。その結果、遅角量切り換え部123により、ΔCAINだけ、目標位相は進角される。   Since a state where surging may occur occurs in the control cycle n, the surge determination flag F_SURGE is set to a value of 1, and the phase delayed by the predetermined value CAIN_PRE from the optimal phase CAINCMDM accordingly. It becomes a phase. By this retarding operation, in the next control cycle (n + 1), there is no possibility of surging, and the surge determination flag F_SURGE is set to zero. As a result, the target phase is advanced by ΔCAIN by the retard amount switching unit 123.

ΔCAINずつの進角動作は、サージング発生のおそれがないと判断される限り行われる。これは、サージング発生のおそれがない限り、現在の位相をなるべく最適位相に近づけるためである。したがって、制御周期(n+2)においてもサージング発生のおそれがないと判断されれば、さらにΔCAINだけ、目標位相は進角されることとなる。この進角動作の結果、最適位相CAINCMDMに達したならば、サージ判定フラグF_SURGEに値1が設定されない限り、目標位相算出部127による上限リミット処理により、目標位相は該最適位相CAINCMDMに維持される。   The advance operation by ΔCAIN is performed as long as it is determined that there is no risk of surging. This is to make the current phase as close as possible to the optimum phase as long as there is no risk of surging. Accordingly, if it is determined that there is no possibility of occurrence of surging even in the control cycle (n + 2), the target phase is advanced by ΔCAIN. As a result of this advance operation, if the optimum phase CAINCMDM is reached, the target phase is maintained at the optimum phase CAINCMDM by the upper limit processing by the target phase calculation unit 127 unless the value 1 is set in the surge determination flag F_SURGE. .

このようにすることにより、サージングの発生を防止しつつ、できる限り最適位相に近づくように、位相を制御することができる。したがって、サージングの発生を防止しつつ、ポンピングロスを低減して、燃費向上を図ることができる。   By doing so, the phase can be controlled so as to be as close to the optimum phase as possible while preventing the occurrence of surging. Accordingly, it is possible to improve the fuel consumption by reducing the pumping loss while preventing the occurrence of surging.

図9は、図7のしきい値設定部115のより詳細なブロック図を示す。変速比取得部131は、変速比センサ63によって検出された現在の変速比を取得する。許容加速度算出部133は、図10の(a)に示すようなマップを、該取得した変速比に基づいて参照し、サージングの振動について許容される加速度(より正確には、許容加速度の上限値)を求める。当該マップは、図5に基づいて作成されており、予めメモリ等の記憶装置に記憶されている。なお、図の横軸の変速比は、右側へ向かうほど小さくなり、これは、右側に向かうほど高速用の変速比であることを表している。   FIG. 9 shows a more detailed block diagram of the threshold value setting unit 115 of FIG. The gear ratio acquisition unit 131 acquires the current gear ratio detected by the gear ratio sensor 63. The allowable acceleration calculation unit 133 refers to a map as shown in FIG. 10A based on the acquired gear ratio, and is allowed for the surging vibration (more precisely, the upper limit value of the allowable acceleration). ) The map is created based on FIG. 5 and stored in advance in a storage device such as a memory. Note that the gear ratio on the horizontal axis in the figure decreases toward the right side, and this indicates that the gear ratio for high speed increases toward the right side.

切り換え部135は、クラッチスイッチ61によって出力される、クラッチが締結されているか否か、すなわちエンジンの回転が車輪に伝達される状態か遮断されている状態かを示す信号を取得する。該信号が、クラッチが締結されていることを示すならば、許容加速度算出部133によって求めた、現在の変速比に対応する許容加速度を出力する。該信号が、クラッチが締結されていないことを示すならば、所定値を出力する。該所定値は、エンジンの回転が車輪に伝達されていない場合における、サージングの振動についての許容加速度であり、これは、予め決められることができる。   The switching unit 135 acquires a signal output by the clutch switch 61 indicating whether or not the clutch is engaged, that is, whether or not the rotation of the engine is transmitted to the wheels or is disconnected. If the signal indicates that the clutch is engaged, the allowable acceleration corresponding to the current gear ratio obtained by the allowable acceleration calculation unit 133 is output. If the signal indicates that the clutch is not engaged, a predetermined value is output. The predetermined value is an allowable acceleration with respect to surging vibration when the rotation of the engine is not transmitted to the wheel, and this can be determined in advance.

しきい値算出部137は、図10の(b)のようなマップから、該制御装置100が搭載される車両の等価慣性重量に対応するマップを選択し(図の例では、3個の等価慣性重量値に対するマップが示されている)、切り換え部135から受け取った許容加速度の値に基づいて該選択したマップを参照し、対応する標準偏差値をしきい値σ_TRQ_SURとして求める。ここで、該車両についての等価慣性重量は、所定値として予め算出されて、メモリ等の記憶装置に記憶されている。また、当該マップは、図6に基づいて作成されており、予めメモリ等の記憶装置に記憶されている。   The threshold value calculation unit 137 selects a map corresponding to the equivalent inertia weight of the vehicle on which the control device 100 is mounted from the map as shown in FIG. 10B (in the example shown in the figure, three equivalents). A map for the inertia weight value is shown), the selected map is referred to based on the allowable acceleration value received from the switching unit 135, and the corresponding standard deviation value is obtained as the threshold σ_TRQ_SUR. Here, the equivalent inertia weight for the vehicle is calculated in advance as a predetermined value and stored in a storage device such as a memory. Further, the map is created based on FIG. 6 and is stored in advance in a storage device such as a memory.

こうして、現在の変速比および車両の等価慣性重量を考慮したしきい値が設定されるので、サージ判定の精度を向上させることができる。すなわち、所定の変速比用のサージング振動の許容加速度を、すべての変速比について採用すると、商品性を満たすことが困難であったり、サージングの判定が過剰に行われることになりうる。たとえば、変速比が4速のものを一律に用いると、実際の変速比が3速の場合には、3速の許容加速度の方が4速の許容加速度より低いため、サージング振動が体感されて商品性が低下するおそれがある。また、実際の変速比が5速の場合には、より厳しいしきい値となってしまい、サージングが発生する可能性があるとの判定が過剰に行われ、よって、位相を進角することができないおそれがある。現在の変速比に応じた許容加速度を算出することで、所望の商品性を維持しつつ、過剰なサージ判定を抑制して、位相を進角させることができる。   In this way, since the threshold value considering the current gear ratio and the vehicle's equivalent inertia weight is set, the accuracy of surge determination can be improved. In other words, if the allowable acceleration of surging vibration for a predetermined gear ratio is adopted for all gear ratios, it may be difficult to satisfy the merchantability or excessive surging may be determined. For example, if gears with a gear ratio of 4th are used uniformly, if the actual gear ratio is 3rd, the 3rd speed allowable acceleration is lower than the 4th speed allowable acceleration, so surging vibration is experienced. There is a risk that the merchantability is reduced. Further, when the actual gear ratio is 5th, the threshold becomes more severe, and it is excessively determined that surging may occur, so that the phase may be advanced. It may not be possible. By calculating the allowable acceleration according to the current gear ratio, it is possible to advance the phase while suppressing excessive surge determination while maintaining desired merchantability.

また、サージングの大きさは、前述したように、車両の等価慣性重量に依存して決まる。したがって、等価慣性重量に適したしきい値を設定することにより、サージ判定の精度を向上させることができる。また、等価慣性重量は、車両によって予め決まっているため、許容加速度に対するしきい値を、等価慣性重量毎に予め算出してマップ上に規定することができる。これにより、各種の車両に対して同じマップを使用することが可能となり、車両毎にマップを異ならせる必要がない。したがって、しきい値を設定するのに、工数を低減することができる。   Further, as described above, the magnitude of surging is determined depending on the equivalent inertia weight of the vehicle. Therefore, the accuracy of surge determination can be improved by setting a threshold value suitable for the equivalent inertia weight. Further, since the equivalent inertia weight is determined in advance by the vehicle, a threshold value for the allowable acceleration can be calculated in advance for each equivalent inertia weight and defined on the map. This makes it possible to use the same map for various types of vehicles, and there is no need for different maps for each vehicle. Therefore, man-hours can be reduced for setting the threshold value.

次に、図11〜図19を参照して、本願発明の一実施形態に従う、ECU1のCPUにより実行されるプロセスを説明する。   Next, with reference to FIGS. 11-19, the process performed by CPU of ECU1 according to one Embodiment of this invention is demonstrated.

図11は、気筒別トルクパラメータの算出プロセスを示し、具体的には、気筒別トルクパラメータ算出部111(図7)により実行される。該プロセスは、TDCパルスの発生に同期して実行される。   FIG. 11 shows a calculation process of the cylinder-specific torque parameter, which is specifically executed by the cylinder-specific torque parameter calculation unit 111 (FIG. 7). The process is performed in synchronization with the generation of the TDC pulse.

この実施例では、クランク角センサ33からのCRKパルスは、クランク角度6度毎に取得される。CRKパルスの発生時間間隔を示す時間パラメータCRME(i)については、クランク角度720度分のデータ(i=0〜ND−1,データ数NDは120個)がECU1のメモリに格納される。また、点火順の気筒識別番号をk(=1〜6)とし、1TDC期間(ここでは、クランク角度120度の期間)内のデータ数をNTDC(この例では、20)とすると、このプロセスの1回の実行で、パラメータiが(k−1)NTDCから(kNTDC−1)までの演算が行われる。たとえば今回の処理が1番目の気筒(k=1)に対応する演算を行うときは、パラメータiは0から(NTDC−1)までの値をとり、5番目の気筒(k=5)に対応する演算を行うときは、パラメータiは4NTDCから(5NTDC−1)までの値をとる。   In this embodiment, the CRK pulse from the crank angle sensor 33 is acquired every 6 degrees of the crank angle. Regarding the time parameter CRME (i) indicating the CRK pulse generation time interval, the data for the crank angle of 720 degrees (i = 0 to ND-1, the number of data ND is 120) is stored in the memory of the ECU 1. Further, if the cylinder identification number in the ignition order is k (= 1 to 6) and the number of data in one TDC period (here, the period of the crank angle is 120 degrees) is NTDC (20 in this example), this process In one execution, calculation is performed from parameter (k-1) NTDC to (kNTDC-1). For example, when the current process performs an operation corresponding to the first cylinder (k = 1), the parameter i takes a value from 0 to (NTDC-1) and corresponds to the fifth cylinder (k = 5). Parameter i takes a value from 4NTDC to (5NTDC-1).

ステップS11において、回転速度OMG(i)(rad/s)を、式(2)に従って算出する。回転速度OMGは、CRKパルスの発生時間間隔を角速度に変換したものであり、クランク軸30の回転する速度を表している。ここで、Dθは、時間パラメータCRMEを計測する角度間隔4π/NDであり、この実施例では、π/30(rad)である。
OMG(i)=Dθ/CRME(i) (2)
In step S11, the rotational speed OMG (i) (rad / s) is calculated according to equation (2). The rotational speed OMG is obtained by converting the CRK pulse generation time interval into an angular speed, and represents the speed at which the crankshaft 30 rotates. Here, Dθ is an angular interval of 4π / ND for measuring the time parameter CRME, and in this embodiment is π / 30 (rad).
OMG (i) = Dθ / CRME (i) (2)

図12(a)には、エンジン回転数NEが上昇していく時の時間パラメータCRMEの推移の一例が示されており、図12(b)には、それに対応する回転速度OMGが示されている。   FIG. 12 (a) shows an example of the transition of the time parameter CRME when the engine speed NE increases, and FIG. 12 (b) shows the corresponding rotational speed OMG. Yes.

ステップS12において、式(3)に従って720度フィルタ処理を実行し、フィルタ処理後回転速度OMGR(i)を算出する。
OMGR(i)=OMG(i)−(OMG(ND)−OMG(0))×Dθ×i/4π (3)
In step S12, a 720 degree filter process is executed according to equation (3), and a post-filter process rotational speed OMGR (i) is calculated.
OMGR (i) = OMG (i) − (OMG (ND) −OMG (0)) × Dθ × i / 4π (3)

720度フィルタ処理は、1サイクル期間(クランク角度720度)における線形変化分をキャンセルし、比較的周期の短い変動を抽出する。この処理は、エンジンの負荷側からエンジンに加わるトルク(エンジンにより駆動される車両のタイヤや補機から加わるトルク、エンジンの摺動部品の摩擦によるトルク等)に起因する回転変動成分を除くために行われる。図12(c)には、図12(b)に基づいて算出されたOMGRの推移が示されている。   The 720 degree filter process cancels the linear change in one cycle period (crank angle 720 degrees), and extracts a fluctuation with a relatively short cycle. This processing is performed to remove rotational fluctuation components caused by torque applied to the engine from the engine load side (torque applied to the tires and accessories driven by the engine, torque due to friction of sliding parts of the engine, etc.). Done. FIG. 12C shows the transition of the OMGR calculated based on FIG.

ステップS13において、式(4)に従って、相対回転速度OMGREFを算出する。
OMGREF(i)=OMGR(i)−OMGR((k-1)NTDC) (4)
In step S13, the relative rotational speed OMGREF is calculated according to the equation (4).
OMGREF (i) = OMGR (i) −OMGR ((k-1) NTDC) (4)

ここで、OMGR((k−1)NTDC)は基準回転速度であり、対象となる気筒の圧縮上死点(TDC)におけるフィルタ処理後回転速度に相当する。   Here, OMGR ((k−1) NTDC) is a reference rotational speed, and corresponds to the post-filtering rotational speed at the compression top dead center (TDC) of the target cylinder.

図13を参照すると、クランク角度が0〜720度の間の相対回転速度OMGREFの挙動の一例が示されている。#1〜#6は、点火順に6つの気筒を識別するのに付された気筒識別番号である。相対回転速度OMGREFは、圧縮上死点後、燃料と空気の混合気の燃焼により一旦上昇した後、減少する。なお、失火が発生すると、相対回転速度は、圧縮上死点後、上昇することなく下降することとなる。   Referring to FIG. 13, an example of the behavior of the relative rotational speed OMGREF when the crank angle is between 0 and 720 degrees is shown. # 1 to # 6 are cylinder identification numbers assigned to identify the six cylinders in the order of ignition. The relative rotational speed OMGREF once rises due to combustion of the fuel / air mixture after compression top dead center and then decreases. If misfire occurs, the relative rotational speed will decrease without increasing after compression top dead center.

相対回転速度OMGREFで囲まれた面積S(網掛けされた部分)は、燃焼によって発生するトルクを表している。したがって、気筒ごとに、相対回転速度OMGREFを1TDC期間(前述したように、この実施例では、クランク角度120度の期間)にわたって積算することにより得られる積算値を、その気筒で発生するトルクを表すトルクパラメータとして用いることができる。   An area S (shaded portion) surrounded by the relative rotational speed OMGREF represents torque generated by combustion. Therefore, for each cylinder, an integrated value obtained by integrating the relative rotational speed OMGREF over a period of 1 TDC (in this embodiment, a period of a crank angle of 120 degrees) represents the torque generated in that cylinder. It can be used as a torque parameter.

好ましくは、上記の積算値を算出する前に、ステップS14およびS15を実施して、慣性トルクの影響を相対回転速度から除去する。慣性トルクは、エンジン12の往復運動部品(ピストンおよびコンロッド)の質量、コンロッドの長さ、クランク半径、クランクプーリ等のエンジンの負荷側の回転部品の慣性力に基づくトルクである。相対回転速度には、慣性力に基づく成分が含まれているが、慣性トルクは、燃焼により生成される出力トルクに寄与するものではないので、これを除去するのが好ましい。   Preferably, before calculating the integrated value, steps S14 and S15 are performed to remove the influence of inertia torque from the relative rotational speed. The inertia torque is a torque based on the inertia force of the rotating parts on the load side of the engine such as the mass of the reciprocating motion parts (piston and connecting rod) of the engine 12, the length of the connecting rod, the crank radius, and the crank pulley. The relative rotational speed includes a component based on an inertial force, but the inertial torque does not contribute to the output torque generated by combustion, and therefore it is preferable to remove this.

ステップS14において、式(5)に従い、各気筒の圧縮上死点における慣性力回転速度OMGI(k)を算出する。
OMGI(k)=K・OMG((k−1)NTDC)/3I (5)
In step S14, the inertial force rotational speed OMGI (k) at the compression top dead center of each cylinder is calculated according to the equation (5).
OMGI (k) = K.OMG ((k-1) NTDC) / 3I (5)

ここで、Kは比例定数であり、Iは、クランクプーリ等の回転部品の慣性モーメントを示す。式(5)の根拠については、後述される。なお、変速機として自動変速機を搭載している場合には、該自動変速機のロックアップクラッチが係合しているか否かに応じて、慣性モーメントIの値を変更することが好ましい。これにより、ロックアップクラッチの係合/非係合にかかわらず、より正確なトルクパラメータの算出を行うことができる。   Here, K is a proportional constant, and I indicates the moment of inertia of a rotating component such as a crank pulley. The basis of equation (5) will be described later. When an automatic transmission is mounted as the transmission, it is preferable to change the value of the moment of inertia I depending on whether or not the lockup clutch of the automatic transmission is engaged. Thus, more accurate torque parameters can be calculated regardless of whether the lockup clutch is engaged or not.

図14(a)は、各気筒の圧縮上死点近傍において、基準回転速度の算出と同じタイミングで算出される慣性力回転速度OMGIの推移の一例を示す。   FIG. 14A shows an example of transition of the inertial force rotational speed OMGI calculated at the same timing as the calculation of the reference rotational speed in the vicinity of the compression top dead center of each cylinder.

ステップS15において、式(6)に従い、修正相対回転速度OMGREFM(i)を算出する。修正相対回転速度は、慣性トルクの影響が除去された相対回転速度である。
OMGREFM(i)=OMGREF(i)+OMGI(k) (6)
In step S15, the corrected relative rotational speed OMGREFM (i) is calculated according to the equation (6). The corrected relative rotational speed is a relative rotational speed from which the influence of inertia torque is removed.
OMGREFM (i) = OMGREF (i) + OMGI (k) (6)

図14(b)は、各気筒の修正相対回転速度の推移の一例を示す。   FIG. 14B shows an example of the transition of the corrected relative rotational speed of each cylinder.

ステップS16において、式(7)に従い、修正相対回転速度OMGREFMを積算して、トルクパラメータMFJUD(k)を算出する。

Figure 2011149312
In step S16, the corrected relative rotational speed OMGREFM is integrated according to the equation (7) to calculate the torque parameter MFJUD (k).
Figure 2011149312

図14(c)は、1TDC期間にわたって修正相対回転速度を積算することにより算出される各気筒のトルクパラメータMFJUDの推移の一例を示す。   FIG. 14C shows an example of the transition of the torque parameter MFJUD of each cylinder calculated by integrating the corrected relative rotational speed over 1 TDC period.

ステップS17において、気筒識別番号kが気筒数Nと等しいか否かを判断する。その答えがNoであるときは、気筒識別番号kを1だけインクリメントし(S18)、答えがYesであるときは、気筒識別番号kを1に戻す(S19)。   In step S17, it is determined whether the cylinder identification number k is equal to the number N of cylinders. When the answer is No, the cylinder identification number k is incremented by 1 (S18), and when the answer is Yes, the cylinder identification number k is returned to 1 (S19).

こうして、TDCパルスの発生に同期して、気筒毎に、該気筒の出力トルクを表すトルクパラメータMFJUDが算出される。   Thus, the torque parameter MFJUD representing the output torque of the cylinder is calculated for each cylinder in synchronization with the generation of the TDC pulse.

ここで、慣性力回転速度OMGIを算出するより詳細な手法を説明する。1つの気筒で発生する慣性力によるトルク(単一気筒慣性トルクと呼ぶ)TI1は、図15に示すようにコンロッド長をL、クランク半径をR、オフセットをe、クランク軸30の回転角速度をω、ピストン28及びコンロッドの合計質量をmとし、角度θおよびφを図のように定義すると、式(8)で与えられる。以下の式中の角度の単位はラジアン(rad)が用いられる。

Figure 2011149312
Here, a more detailed method for calculating the inertial force rotation speed OMGI will be described. As shown in FIG. 15, the torque TI1 due to the inertial force generated in one cylinder is L, the connecting rod length is R, the crank radius is R, the offset is e, and the rotational angular velocity of the crankshaft 30 is ω. When the total mass of the piston 28 and the connecting rod is m, and the angles θ and φ are defined as shown in the figure, the equation (8) is given. The unit of angle in the following formula is radians (rad).
Figure 2011149312

図16(a)は、式(8)により算出される単一気筒慣性トルクTI1を、クランク角度θの関数として示す。単一気筒慣性トルクTI1の位相を120度ずらして6気筒分を加算した合成慣性トルクTIは、図14(b)のように推移し、式(9)で近似することができる。
TI=−Asin3θ (9)
ここで、Aは、回転角速度ω(rad/s)の2乗に比例する係数である。
FIG. 16A shows the single cylinder inertia torque TI1 calculated by the equation (8) as a function of the crank angle θ. The combined inertia torque TI obtained by shifting the phase of the single cylinder inertia torque TI1 by 120 degrees and adding six cylinders changes as shown in FIG. 14B and can be approximated by the equation (9).
TI = −Asin3θ (9)
Here, A is a coefficient proportional to the square of the rotational angular velocity ω (rad / s).

一方、前述したようにIを慣性モーメントとすると、合成慣性トルクTIは、式(10)で与えられる。
TI=I×(dω/dt) (10)
On the other hand, as described above, when I is an inertia moment, the combined inertia torque TI is given by Expression (10).
TI = I × (dω / dt) (10)

式(9)および(10)から式(11)が得られ、これを回転角速度ωについて解くと、合成トルクTIに対応する慣性力回転速度ωIは、式(12)で与えられる。
−Asin3θ=I×(dω/dt) (11)
ωI=(Acos3θ×dt/dθ)/3I (12)
When Expression (11) is obtained from Expressions (9) and (10) and is solved with respect to the rotational angular speed ω, the inertial force rotational speed ωI corresponding to the combined torque TI is given by Expression (12).
−Asin3θ = I × (dω / dt) (11)
ωI = (Acos3θ × dt / dθ) / 3I (12)

図14(c)は、慣性力回転速度ωIの推移を示す。   FIG. 14C shows the transition of the inertial force rotation speed ωI.

よって、圧縮上死点での慣性力回転速度OMGIは、式(12)のθをゼロとして、式(13)により算出することができる。
OMGI=(A/3I)(1/OMG) (13)
Therefore, the inertial force rotational speed OMGI at the compression top dead center can be calculated by the equation (13) with θ in the equation (12) being zero.
OMGI = (A / 3I) (1 / OMG) (13)

係数Aは、回転速度OMGの2乗に比例するので、比例定数をKとすると、式(13)は、式(14)のように変形することができる。
OMGI=K・OMG/3I (14)
Since the coefficient A is proportional to the square of the rotational speed OMG, if the proportionality constant is K, the equation (13) can be transformed into the equation (14).
OMGI = K · OMG / 3I (14)

したがって、各気筒の圧縮上死点における慣性力回転速度OMGIは、前述した式(5)のように表されることができる。   Therefore, the inertial force rotational speed OMGI at the compression top dead center of each cylinder can be expressed as the above-described equation (5).

図16(c)に示すように、圧縮上死点(θ=0、120、240、・・・)での慣性力回転速度OMGIは最大の値となるので、相対回転速度OMGREFに慣性力回転速度OMGIを加算すること(これは、基準回転速度から慣性力回転速度OMGIを減算することと等価である)により、慣性力回転速度ωIの影響を除去した修正相対回転速度OMGREFMを得ることができる。なお、図16(c)に示す慣性力回転速度ωIの周期変動成分は、修正相対回転速度OMGREFMを1TDC期間(この実施例では、120度)にわたって積算することによりキャンセルされる。   As shown in FIG. 16 (c), the inertial force rotation speed OMGI at the compression top dead center (θ = 0, 120, 240,...) Is the maximum value, and therefore the inertial force rotation is equal to the relative rotation speed OMGREF. By adding the speed OMGI (this is equivalent to subtracting the inertial force rotational speed OMGI from the reference rotational speed), it is possible to obtain a corrected relative rotational speed OMGREFM that eliminates the influence of the inertial force rotational speed ωI. . Note that the periodic fluctuation component of the inertial force rotational speed ωI shown in FIG. 16C is canceled by integrating the corrected relative rotational speed OMGREFM over 1 TDC period (120 degrees in this embodiment).

図17は、或る運転状態におけるトルクパラメータMFJUDの値を示す。(a)は、気筒間のトルクパラメータのバラツキがほとんど無い状態を示し、(b)は、気筒#1が、他の気筒に比べて、トルクパラメータの偏差が大きい状態を示す。   FIG. 17 shows the value of the torque parameter MFJUD in a certain operating state. (A) shows a state in which there is almost no variation in torque parameters between cylinders, and (b) shows a state in which cylinder # 1 has a larger deviation in torque parameters than other cylinders.

次に、図18は、サージ判定に基づく目標位相の算出プロセスを示し、該プロセスは、より具体的には、標準偏差算出部113、しきい値設定部115、サージ判定部117、最適位相決定部121、遅角量切り換え部123、位相補正部125、および目標位相算出部127(図7)によって実行される。このプロセスは、たとえば所定の時間間隔で実行されることができる。   Next, FIG. 18 shows a target phase calculation process based on surge determination, and more specifically, the standard deviation calculation unit 113, threshold setting unit 115, surge determination unit 117, optimum phase determination. This is executed by the unit 121, the retard amount switching unit 123, the phase correction unit 125, and the target phase calculation unit 127 (FIG. 7). This process can be performed, for example, at predetermined time intervals.

ステップS101において、図11で算出された気筒別トルクパラメータMFJUDについて、直近の所定数iの制御周期で算出された気筒別トルクパラメータMFJUDiを取得する。iが100であれば、100個の気筒別トルクパラメータMFJUDi(i=1〜100)が取得される。これらの気筒別トルクパラメータMFJUDiに基づいて、前述した式(1)に示すように、該気筒別トルクパラメータの標準偏差σ_TRQを算出する。   In step S101, for each cylinder torque parameter MFJUD calculated in FIG. 11, the cylinder specific torque parameter MFJUDi calculated in the most recent predetermined number i of control cycles is acquired. If i is 100, 100 cylinder-specific torque parameters MFJUDi (i = 1 to 100) are acquired. Based on these cylinder-specific torque parameters MFJUDI, the standard deviation σ_TRQ of the cylinder-specific torque parameters is calculated as shown in the above-described equation (1).

ステップS102において、クラッチスイッチ61の検出結果を介して、クラッチが締結されているかどうか(すなわち、エンジン回転が車輪に伝達されている状態か遮断されている状態か)を判断する。クラッチが締結されていれば、変速比センサ63を介して検出された現在の変速比を取得し、該変速比に基づいて、図10(a)に示すようなマップを参照し、対応するサージングによる振動の許容加速度を求める(S103)。クラッチが締結されていなければ、予め設定された所定値を、該許容加速度に設定する(S104)。   In step S102, it is determined, based on the detection result of the clutch switch 61, whether or not the clutch is engaged (that is, whether the engine rotation is being transmitted to the wheel or disconnected). If the clutch is engaged, the current gear ratio detected via the gear ratio sensor 63 is acquired, and based on the gear ratio, a map as shown in FIG. The allowable acceleration of vibration due to is obtained (S103). If the clutch is not engaged, a predetermined value set in advance is set as the allowable acceleration (S104).

ステップS105において、図10(b)に示すように、車両の等価慣性重量に対応するマップを選択し、該選択したマップを、ステップS103またはS104で求めた許容加速度に基づいて参照し、対応する標準偏差値を、しきい値σ_TRQ_SURとして求める。   In step S105, as shown in FIG. 10B, a map corresponding to the equivalent inertia weight of the vehicle is selected, and the selected map is referenced based on the allowable acceleration obtained in step S103 or S104. A standard deviation value is obtained as a threshold value σ_TRQ_SUR.

ステップS106において、ステップS101で算出された標準偏差σ_TRQと、しきい値σ_TRQ_SURとを比較する。標準偏差σ_TRQがしきい値σ_TRQ_SURより大きければ、許容されない大きさのサージングの発生の可能性がある状態と判定し、ステップS107において、サージ判定フラグF_SURGEに値1を設定する。標準偏差σ_TRQがしきい値σ_TRQ_SUR以下であれば、許容されない大きさのサージングの発生の可能性はない状態とステップS108において、サージ判定フラグF_SURGEにゼロを設定する。   In step S106, the standard deviation σ_TRQ calculated in step S101 is compared with the threshold value σ_TRQ_SUR. If the standard deviation σ_TRQ is larger than the threshold σ_TRQ_SUR, it is determined that there is a possibility of occurrence of unacceptable magnitude surging, and a value 1 is set to the surge determination flag F_SURGE in step S107. If standard deviation σ_TRQ is equal to or less than threshold value σ_TRQ_SUR, zero is set to surge determination flag F_SURGE in a state where there is no possibility of occurrence of unacceptable magnitude surging.

図19に進み、ステップS111において、現在の運転状態に基づいて、所定のマップを参照することにより、位相の最適値(最適位相)CAINCMDMを決定する。ステップS112において、サージ判定フラグF_SURGEに値1が設定されているか否かを判断する。値1が設定されていれば、遅角量ΔCAIN_SURに、所定値CAIN_PREを設定する(S113)。サージ判定フラグF_SURGEにゼロが設定されていれば、遅角量ΔCAIN_SURに、(所定値CAIN_PRE―ΔCAIN)を設定する(S114)。ここで、ΔCAINは所定量(<CAIN_PRE)であり、図8を参照して前述したように、サージングの発生の可能性が解消されたときに各制御周期で最適位相に向けて進角させる量である。   Proceeding to FIG. 19, in step S111, an optimum value of phase (optimum phase) CAINCMDM is determined by referring to a predetermined map on the basis of the current operating state. In step S112, it is determined whether or not a value 1 is set in the surge determination flag F_SURGE. If the value 1 is set, a predetermined value CAIN_PRE is set to the retardation amount ΔCAIN_SUR (S113). If the surge determination flag F_SURGE is set to zero, (predetermined value CAIN_PRE−ΔCAIN) is set to the retardation amount ΔCAIN_SUR (S114). Here, ΔCAIN is a predetermined amount (<CAIN_PRE), and as described above with reference to FIG. 8, when the possibility of occurrence of surging is eliminated, the amount to be advanced toward the optimum phase in each control cycle It is.

ステップS115において、最適位相CAINCMDMから、遅角量ΔCAIN_SURを減算することにより、補正済み位相CAINCMDXを算出する。ステップS116において、補正済み位相CAINCMDXを、最適値CAINCMDMで上限リミット処理を行うことにより、目標位相CAINCMDを算出する。こうして、吸気バルブ14の位相は、該目標位相CAINCMDに達するように、可変動弁機構31を介して制御される。   In step S115, the corrected phase CAINCMDX is calculated by subtracting the retardation amount ΔCAIN_SUR from the optimum phase CAINCMDM. In step S116, the target phase CAINCMD is calculated by performing upper limit processing on the corrected phase CAINCMDX with the optimum value CAINCMDM. Thus, the phase of the intake valve 14 is controlled via the variable valve mechanism 31 so as to reach the target phase CAINCMD.

以下は、気筒別トルクパラメータ算出部111により実行されるトルクパラメータMFJUDを算出する手法の代替形態の説明である。   The following is an explanation of an alternative form of the method for calculating the torque parameter MFJUD executed by the cylinder specific torque parameter calculation unit 111.

図20は、図11に示すプロセス(第1の形態)の変形例を示し、ステップS15およびS16を、ステップS16a、S16bおよびS16cに変更したものである。   FIG. 20 shows a modification of the process (first form) shown in FIG. 11, in which steps S15 and S16 are changed to steps S16a, S16b and S16c.

ステップS16aでは、式(15)に従い、相対回転速度OMGREF(i)の積算値として、トルクパラメータMFJUDa(k)を算出する。

Figure 2011149312
In step S16a, the torque parameter MFJUDa (k) is calculated as an integrated value of the relative rotational speed OMGREF (i) according to the equation (15).
Figure 2011149312

ステップS16bでは、式(16)に従い、慣性力回転速度OMGIの積算値MFTH(k)を算出する。
MFTH(k)=−NTDC×OMGI(k) (16)
In step S16b, an integrated value MFTH (k) of inertial force rotation speed OMGI is calculated according to equation (16).
MFTH (k) = − NTDC × OMGI (k) (16)

ステップS16cでは、式(17)に従い、トルクパラメータMFJUDを算出する。
MFJUD(k)=MFJUDa(k)−MFTH(k) (17)
In step S16c, torque parameter MFJUD is calculated according to equation (17).
MFJUD (k) = MFJUDa (k) −MFTH (k) (17)

このように、この変形例では、修正相対回転速度を積算することに代えて、相対回転速度を積算した値と慣性力回転速度を積算した値とを用いてトルクパラメータMFJUDを算出する。   Thus, in this modification, instead of integrating the corrected relative rotational speed, the torque parameter MFJUD is calculated using the value obtained by integrating the relative rotational speed and the value obtained by integrating the inertial force rotational speed.

図21は、トルクパラメータを算出するプロセスの第2の形態のフローを示す。図11および図20を参照して説明した第1の形態では、時間パラメータCRMEを回転速度OMGに変換したが、この形態では、時間パラメータCRMEを速度パラメータとして使用する。以下のプロセス中で算出される相対時間パラメータCRMEREFの積算値は、図13に一例として示された相対回転速度OMGREFの積算値と同様の変化を呈するので、トルクを表すパラメータとして用いることができる。なお、以下に説明する点以外は、第1の形態と同じである。   FIG. 21 shows a flow of the second form of the process of calculating the torque parameter. In the first mode described with reference to FIGS. 11 and 20, the time parameter CRME is converted into the rotational speed OMG. In this mode, the time parameter CRME is used as the speed parameter. Since the integrated value of the relative time parameter CRMEREF calculated in the following process exhibits the same change as the integrated value of the relative rotational speed OMGREF shown as an example in FIG. 13, it can be used as a parameter representing torque. The first embodiment is the same as the first embodiment except for the points described below.

ステップS32において、式(18)に従って720度フィルタ処理を実行し、フィルタ処理後時間パラメータCRMER(i)を算出する。
CRMER(i)=CRME(i)−(CRME(0)−CRME(ND))×Dθ×i/4π (18)
In step S32, a 720 degree filter process is executed according to the equation (18), and a post-filter process time parameter CRMER (i) is calculated.
CRMER (i) = CRME (i) − (CRME (0) −CRME (ND)) × Dθ × i / 4π (18)

ステップS33において、式(19)に従い、相対時間パラメータCRMEREF(i)を算出する。ここで、CRMER((k−1)NTDC)は基準時間パラメータであり、対象となる気筒の圧縮上死点におけるフィルタ処理後時間パラメータに相当する。
CRMEREF(i)=CRMER((k−1)NTDC)−CRMER(i) (19)
In step S33, the relative time parameter CRMEREF (i) is calculated according to the equation (19). Here, CRMER ((k−1) NTDC) is a reference time parameter, and corresponds to a time parameter after filtering at the compression top dead center of the target cylinder.
CRMEREF (i) = CRMER ((k-1) NTDC) -CRMER (i) (19)

ステップS34において、式(20)に従い、慣性力時間パラメータCRMEI(k)を算出する。
CRMEI(k)=3I・CRME((k−1)NTDC)/K (20)
In step S34, the inertia force time parameter CRMEI (k) is calculated according to the equation (20).
CRMEI (k) = 3I.CRME ((k-1) NTDC) / K (20)

ステップS35において、式(21)に従い、修正相対時間パラメータCRMEREFM(i)を算出する。
CRMEREFM(i)=CRMEREF(i)−CRMEI(k) (21)
In step S35, the corrected relative time parameter CRMEREFM (i) is calculated according to the equation (21).
CRMEFM (i) = CRMEREF (i)-CRMEI (k) (21)

ステップS36において、式(22)に従い、修正相対時間パラメータCRMEREFMの積算値を算出することにより、トルクを表すトルクパラメータMFJUD(k)を算出する。

Figure 2011149312
In step S36, the torque parameter MFJUD (k) representing the torque is calculated by calculating the integrated value of the corrected relative time parameter CRMEREFM according to the equation (22).
Figure 2011149312

ステップS37では、気筒識別番号kが気筒数Nと等しいか否かを判断する。答えがNoであるときは、気筒識別番号kを1だけインクリメントし(S38)、答えがYesであるときは、気筒識別番号kを1に戻す(S39)。   In step S37, it is determined whether or not the cylinder identification number k is equal to the number N of cylinders. When the answer is No, the cylinder identification number k is incremented by 1 (S38), and when the answer is Yes, the cylinder identification number k is returned to 1 (S39).

図22は、図21に示すプロセスの変形例を示し、ステップS36を、ステップS36a、S36bおよびS36cに変更したものである。   FIG. 22 shows a modification of the process shown in FIG. 21, in which step S36 is changed to steps S36a, S36b, and S36c.

ステップS36aでは、式(23)により、相対時間パラメータの積算値MFJUDc(k)を算出する。

Figure 2011149312
In step S36a, the integrated value MFJUDc (k) of the relative time parameter is calculated by equation (23).
Figure 2011149312

ステップS36bでは、式(24)に従い、CRMEIの積算値MFTHa(k)を算出する。
MFTHa(k)=NTDC×CRMEI(k) (24)
In step S36b, the CRMEI integrated value MTFHa (k) is calculated according to the equation (24).
MFTHa (k) = NTDC × CRMEI (k) (24)

ステップS36cにおいて、式(25)に従い、トルクパラメータMFJUDを算出する。
MFJUD(k)=MFJUDc−MFTHa(k) (25)
In step S36c, torque parameter MFJUD is calculated according to equation (25).
MFJUD (k) = MFJUDc-MFTHa (k) (25)

この変形例では、修正相対時間パラメータCRMEREFMを積算することに代えて、相対時間パラメータCRMEREFの積算値と慣性力時間パラメータCRMEIの積算値とを用いて、トルクパラメータMFJUDを算出する。   In this modified example, instead of integrating the corrected relative time parameter CRMEREFM, the torque parameter MFJUD is calculated using the integrated value of the relative time parameter CRMEREF and the integrated value of the inertial force time parameter CRMEI.

トルクパラメータの算出は、様々な変形が可能である。たとえば、時間パラメータCRME(i)を式(2)に適用して回転速度OMGを算出したが、高回転時に算出精度が低下しないようにするため、式(26)により算出される5個の時間パラメータCRMEの積算値CRME30(i)を用いて回転速度OMGを算出するようにしてもよい。

Figure 2011149312
The torque parameter can be calculated in various ways. For example, the rotational speed OMG is calculated by applying the time parameter CRME (i) to the equation (2), but the five times calculated by the equation (26) are used so that the calculation accuracy does not decrease at the time of high rotation. The rotational speed OMG may be calculated using the integrated value CRME30 (i) of the parameter CRME.
Figure 2011149312

この場合、回転速度OMG(i)は、式(27)により算出される。ただし、回転速度の算出位相がずれるので、その分の位相補正を行うのが良い。
OMG(i)=5Dθ/CRME30(i) (27)
In this case, the rotational speed OMG (i) is calculated by Expression (27). However, since the calculated phase of the rotational speed is shifted, it is preferable to perform the phase correction accordingly.
OMG (i) = 5Dθ / CRME30 (i) (27)

また、上記の基準回転速度および基準時間パラメータは、圧縮上死点における回転速度および時間パラメータを用いたが、圧縮上死点近傍(例えば、±7.5度の範囲内)でもよい。ここで、7.5度は、回転速度パラメータのサンプリング周期が15度の場合に対応しており、一般的にサンプリング周期をθSPLとすると、±θSPL/2の範囲内でサンプリングされた回転速度パラメータを用いるのが良い。   Further, although the rotation speed and time parameters at the compression top dead center are used as the reference rotation speed and the reference time parameter, they may be near the compression top dead center (for example, within a range of ± 7.5 degrees). Here, 7.5 degrees corresponds to the case where the sampling period of the rotation speed parameter is 15 degrees. Generally, assuming that the sampling period is θSPL, the rotation speed parameter sampled within a range of ± θSPL / 2. It is good to use.

また、720度フィルタ処理は、上記の式(3)に代えて、下記の式(28)により行ってもよい。この式は、クランク角720度の期間の回転速度OMGの移動平均値OMGAVE(m)を用いて線形変化分をキャンセルするものである。ここで、mは、クランク角度720度の周期に対応する離散化時刻である。
OMGR(i)=OMG(i)
―(OMGAVE(m)−OMGAVE(m-1))×Dθ×i/4π (28)
Further, the 720 degree filter processing may be performed by the following equation (28) instead of the above equation (3). This equation cancels the linear change using the moving average value OMGAVE (m) of the rotational speed OMG during the period of the crank angle of 720 degrees. Here, m is a discretization time corresponding to a cycle of a crank angle of 720 degrees.
OMGR (i) = OMG (i)
― (OMGAVE (m) -OMGAVE (m-1)) × Dθ × i / 4π (28)

次に、気筒別トルクパラメータ算出部111により実行される、トルクパラメータを算出する手法の第3の形態について説明する。この形態では、クランク軸の捩れやクランク角センサによる時間パラメータCRMEの検出誤差等に起因する外乱の影響を排除することができる。   Next, a third embodiment of a method for calculating torque parameters, which is executed by the cylinder-specific torque parameter calculation unit 111, will be described. In this embodiment, it is possible to eliminate the influence of disturbance due to the crankshaft twist, the detection error of the time parameter CRME by the crank angle sensor, and the like.

図23(a)は、修正相対回転速度OMGREFMの実測値の一例を示しており、点線で囲んだ部分が、このような外乱の影響を受けた部分である。このような外乱の影響があると、トルクパラメータの算出に誤差が生じるおそれがある。   FIG. 23A shows an example of an actual measurement value of the corrected relative rotational speed OMGREFM, and a portion surrounded by a dotted line is a portion affected by such a disturbance. If there is an influence of such disturbance, an error may occur in calculation of the torque parameter.

そこで、この実施形態では、正常燃焼が行われ、かつクランク角センサの検出に影響を与える外乱がない場合の回転速度変化を近似する燃焼相関関数FCRを修正相対回転速度OMGREFMに乗算することにより、外乱の影響を排除する。図23(b)は、(a)に示す修正相対回転速度OMGREFMに関数FCRを乗算することにより算出したOMGREFMbを示しており、外乱が抑制されていることがわかる。   Therefore, in this embodiment, by multiplying the corrected relative rotational speed OMGREFM by the combustion correlation function FCR that approximates the rotational speed change when normal combustion is performed and there is no disturbance affecting the detection of the crank angle sensor, Eliminate the effects of disturbances. FIG. 23B shows OMGREFMb calculated by multiplying the corrected relative rotational speed OMGREFM shown in FIG. 23A by the function FCR, and it can be seen that the disturbance is suppressed.

図24は、相関関数FCRの一例を示し、式(29)で定義される。ここで、Nは気筒数であり、θは、特定気筒のピストンが上死点に位置する角度を基準としたクランク角度である(図15を参照)。なお、図24は、6気筒エンジンに対応する相関関数FCRを示している。
FCR={1−2cos(N・θ/2)}/2 (29)
FIG. 24 shows an example of the correlation function FCR, which is defined by Expression (29). Here, N is the number of cylinders, and θ is a crank angle based on the angle at which the piston of the specific cylinder is located at the top dead center (see FIG. 15). FIG. 24 shows a correlation function FCR corresponding to a 6-cylinder engine.
FCR = {1-2 cos (N · θ / 2)} / 2 (29)

相関関数FCRは、例えばエンジンの暖機後の定常運転状態において、正常燃焼時の各気筒の筒内圧を計測し、計測した気筒毎の筒内圧を加算することにより合成の筒内圧変化を算出し、その合成筒内圧変化を回転速度の変化に換算することにより、求めてもよい。図25は、そのようにして求めた関数FCRを示す。この関数FCRは、正常燃焼状態における回転速度変化波形を、最小値が0で最大値が1となるように正規化されている。   The correlation function FCR calculates the combined in-cylinder pressure change by, for example, measuring the in-cylinder pressure of each cylinder during normal combustion in the steady operation state after the engine is warmed up, and adding the measured in-cylinder pressure for each cylinder. Alternatively, it may be obtained by converting the change in the combined cylinder pressure into the change in the rotational speed. FIG. 25 shows the function FCR thus obtained. This function FCR is normalized so that the rotation speed change waveform in the normal combustion state is 0 at the minimum value and 1 at the maximum value.

図26(a)は、燃焼相関関数FCRによる相対回転速度の補正を行わない場合のトルクパラメータMFJUDのばらつきの範囲の例を示し、図26(b)は、この第3の実施形態におけるトルクパラメータMFJUDのばらつき範囲の例を示す。これらの図から明らかなように、相関関数FCRを用いた補正を行うことにより、トルクパラメータの算出精度が高まり、ばらつき範囲を減少させることができる(図の例では、約40%減少)。   FIG. 26A shows an example of the range of variation of the torque parameter MFJUD when the relative rotational speed is not corrected by the combustion correlation function FCR, and FIG. 26B shows the torque parameter in the third embodiment. An example of the variation range of MFJUD is shown. As is clear from these figures, by performing correction using the correlation function FCR, the calculation accuracy of the torque parameter can be increased and the variation range can be reduced (in the example of the figure, it is reduced by about 40%).

図27は、この第3の形態に従う、トルクパラメータを算出するプロセスのフローを示す。   FIG. 27 shows a process flow for calculating a torque parameter according to the third embodiment.

ステップS51〜53は、第1の形態に従う図11のステップS11〜S13と同じであるので、説明を省略する。   Steps S51 to S53 are the same as steps S11 to S13 of FIG. 11 according to the first embodiment, and thus description thereof is omitted.

ステップS54において、式(5)により算出される慣性力回転速度OMGI(k)を下記の式(30)に適用し、慣性力回転速度OMGIa(I)を算出する。第1の形態では、圧縮上死点での慣性力回転速度OMGI(k)をそのまま式(5)に適用して修正相対回転速度OMGREFMを算出したが、この形態では、各サンプリングタイミングにおける慣性力回転速度OMGIa(i)を算出して、相対回転速度OMGREFの修正を行う。
OMGIa(i)=OMGI(k)×{cos(N・Dθ・i/2)−1} (30)
In step S54, the inertial force rotation speed OMGI (k) calculated by the equation (5) is applied to the following equation (30) to calculate the inertial force rotation speed OMGIa (I). In the first form, the inertial force rotational speed OMGI (k) at the compression top dead center is applied to the equation (5) as it is to calculate the corrected relative rotational speed OMGREFM. In this form, the inertial force at each sampling timing is calculated. The rotational speed OMGIa (i) is calculated, and the relative rotational speed OMGREF is corrected.
OMGIa (i) = OMGI (k) × {cos (N · Dθ · i / 2) −1} (30)

ステップS55では、式(31)に、ステップS54で算出した慣性力回転速度OMGIa(i)を適用し、第1修正相対回転速度OMGREFMa(i)を算出する。
OMGREFMa(i)=OMGREF(i)−OMGIa(i) (31)
In step S55, the inertial force rotational speed OMGIa (i) calculated in step S54 is applied to the equation (31) to calculate the first corrected relative rotational speed OMGREFMa (i).
OMGREFMa (i) = OMGREF (i) -OMGIa (i) (31)

ステップS56では、ステップS55で算出した第1修正相対回転速度OMGREFMa(i)および式(32)により算出される相関関数FCR(i)を式(33)に適用し、第2修正相対回転速度OMGREFMb(i)を算出する。式(32)は、式(29)のθを、(Dθ・i)に置換したものである。
FCR(i)={1−2cos(N・Dθ・i/2)}/2 (32)
OMGREFMb(i)=OMGREFMa(i)×FCR(i) (33)
In step S56, the first corrected relative rotational speed OMGREFMa (i) calculated in step S55 and the correlation function FCR (i) calculated by the expression (32) are applied to the expression (33) to obtain the second corrected relative rotational speed OMGREFFMb. (I) is calculated. Expression (32) is obtained by replacing θ in Expression (29) with (Dθ · i).
FCR (i) = {1-2 cos (N · Dθ · i / 2)} / 2 (32)
OMGREFMb (i) = OMGREFMa (i) × FCR (i) (33)

ステップS57では、式(34)に従ってトルクパラメータMFJUDを算出する。

Figure 2011149312
In step S57, torque parameter MFJUD is calculated according to equation (34).
Figure 2011149312

こうして、燃焼関数を用いることにより、クランク角センサの検出値に影響を与えうる外乱の影響を排除しつつトルクパラメータを算出することができる。また、式(32)に示す燃焼相関関数を用いることにより、燃焼相関関数値算出用のテーブル設定のための実験が不要となり、比較的簡単な演算でトルクパラメータの算出補正を行うことができる。   Thus, by using the combustion function, it is possible to calculate the torque parameter while eliminating the influence of disturbance that may affect the detection value of the crank angle sensor. Further, by using the combustion correlation function shown in the equation (32), an experiment for setting a table for calculating the combustion correlation function value becomes unnecessary, and the calculation correction of the torque parameter can be performed with a relatively simple calculation.

代替的に、図25に示す実測データに基づく相関関数を用いる場合には、図25に示す1周期分の関数値FCR(i)をパラメータiに応じて検索するFCRテーブルを予めメモリに格納しておき、ステップS56において、式(32)による演算に代えて、FCRテーブル検索を行う。実測データに基づく燃焼相関関数を用いることにより、燃焼相関関数に内燃機関の特性を反映するようトルクパラメータの算出を補正することができる。   Alternatively, when the correlation function based on the actual measurement data shown in FIG. 25 is used, an FCR table for searching the function value FCR (i) for one period shown in FIG. 25 according to the parameter i is stored in the memory in advance. In step S56, an FCR table search is performed instead of the calculation according to the equation (32). By using the combustion correlation function based on the actually measured data, the calculation of the torque parameter can be corrected so that the characteristic of the internal combustion engine is reflected in the combustion correlation function.

また、式(32)の演算も、コサイン関数を予めテーブルとしてメモリに記憶しておき、そのテーブルを検索することにより、相関関数値FCR(i)を算出するようにしてもよい。FCRを用いた補正は、第2の形態にも適用可能である。   Further, in the calculation of Expression (32), the correlation function value FCR (i) may be calculated by storing the cosine function in a memory in advance as a table and searching the table. Correction using the FCR is also applicable to the second embodiment.

上述した実施形態では、6気筒エンジンを例に説明したが、本願発明は、任意の数の気筒を有するエンジンに適用可能である。また、本願発明は、直接噴射式のエンジン、ディーゼルエンジン等のエンジンにも適用可能である。   In the above-described embodiment, a six-cylinder engine has been described as an example, but the present invention can be applied to an engine having an arbitrary number of cylinders. The present invention is also applicable to engines such as direct injection engines and diesel engines.

さらに、本発明は、汎用の(例えば、船外機等の)内燃機関に適用可能である。   Furthermore, the present invention is applicable to general-purpose internal combustion engines (for example, outboard motors).

10 ECU
12 エンジン
14 吸気バルブ
33 クランク角センサ
100 制御装置
10 ECU
12 Engine 14 Intake valve 33 Crank angle sensor 100 Control device

Claims (4)

車両に搭載され、吸気弁と排気弁のオーバーラップ期間を変更可能な可変動弁機構を有する内燃機関の制御装置であって、
内燃機関の回転速度に応じた回転速度パラメータを検出する回転速度パラメータ検出手段と、
前記回転速度パラメータに基づいて、気筒別の出力トルクを表す気筒別トルクパラメータを算出する手段と、
前記算出された気筒別トルクパラメータの標準偏差を算出する手段と、
動力伝達機構の変速比毎に、サージングによる振動の許容加速度を規定した許容加速度マップ、および車両の等価慣性重量毎に、該許容加速度に対応する気筒別トルクパラメータの標準偏差を規定した標準偏差マップを記憶する記憶手段と、
前記車両に搭載された動力伝達機構について検出された現在の変速比に基づいて前記許容加速度マップを参照し、対応するサージングの許容加速度を求める手段と、
前記求めた許容加速度と前記車両の等価慣性重量に基づいて前記標準偏差マップを参照し、対応する標準偏差をしきい値として求める手段と、
前記算出された標準偏差が前記求めたしきい値より大きければ、サージングが発生する可能性がある状態と判定する手段と、
前記サージングが発生する可能性がある状態と判定されたならば、前記オーバーラップ期間が小さくなるように、前記可変動弁機構を制御する制御手段と、
を備える制御装置。
A control device for an internal combustion engine that is mounted on a vehicle and has a variable valve mechanism that can change an overlap period of an intake valve and an exhaust valve,
A rotational speed parameter detecting means for detecting a rotational speed parameter according to the rotational speed of the internal combustion engine;
Means for calculating a cylinder-specific torque parameter representing an output torque for each cylinder based on the rotational speed parameter;
Means for calculating a standard deviation of the calculated cylinder-specific torque parameter;
A permissible acceleration map that prescribes the permissible acceleration of vibration due to surging for each gear ratio of the power transmission mechanism, and a standard deviation map that prescribes the standard deviation of the cylinder-specific torque parameter corresponding to the permissible acceleration for each equivalent inertia weight of the vehicle Storage means for storing
Means for referring to the allowable acceleration map based on a current gear ratio detected for a power transmission mechanism mounted on the vehicle and determining a corresponding surging allowable acceleration;
Means for referring to the standard deviation map based on the obtained allowable acceleration and the equivalent inertia weight of the vehicle, and obtaining a corresponding standard deviation as a threshold value;
Means for determining that surging may occur if the calculated standard deviation is greater than the determined threshold;
If it is determined that the surging is likely to occur, control means for controlling the variable valve mechanism so that the overlap period is reduced;
A control device comprising:
前記内燃機関の運転状態に基づいて、吸気バルブの最適位相を求める手段を備え、
前記制御手段は、前記サージングが発生する可能性がある状態と判定されたならば、前記吸気バルブの位相を、前記最適位相から第1の所定量だけ遅角させる遅角制御を、前記可変動弁機構を介して行うことにより、前記オーバーラップ期間を小さくする、
請求項1に記載の制御装置。
Means for determining an optimum phase of the intake valve based on the operating state of the internal combustion engine;
If it is determined that the surging is likely to occur, the control means performs a delay angle control for delaying the phase of the intake valve by a first predetermined amount from the optimum phase. Reducing the overlap period by performing via a valve mechanism;
The control device according to claim 1.
前記制御手段は、さらに、前記サージングが発生する可能性がある状態と判定されたことに応じた前記遅角制御の後に、該サージングが発生する可能性のある状態が解消されたならば、前記吸気バルブの位相を、前記第1の所定量だけ遅角させた位相から、第2の所定量ずつ前記最適位相に向けて進角させる進角制御を、前記可変動弁機構を介して行う、
請求項2に記載の制御装置。
The control means further, after the retardation control in response to the determination that the surging is likely to occur, if the state where the surging may occur is resolved, An advance angle control for advancing the phase of the intake valve from the phase delayed by the first predetermined amount toward the optimum phase by a second predetermined amount is performed via the variable valve mechanism.
The control device according to claim 2.
前記制御手段は、さらに、前記進角制御によって進角された前記吸気バルブの位相に対し、前記最適位相で上限リミット処理する、
請求項3に記載の制御装置。
The control means further performs an upper limit process at the optimum phase with respect to the phase of the intake valve advanced by the advance angle control.
The control device according to claim 3.
JP2010010225A 2010-01-20 2010-01-20 Control device for internal combustion engine Expired - Fee Related JP5351064B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010010225A JP5351064B2 (en) 2010-01-20 2010-01-20 Control device for internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010010225A JP5351064B2 (en) 2010-01-20 2010-01-20 Control device for internal combustion engine

Publications (2)

Publication Number Publication Date
JP2011149312A true JP2011149312A (en) 2011-08-04
JP5351064B2 JP5351064B2 (en) 2013-11-27

Family

ID=44536516

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010010225A Expired - Fee Related JP5351064B2 (en) 2010-01-20 2010-01-20 Control device for internal combustion engine

Country Status (1)

Country Link
JP (1) JP5351064B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2495511A (en) * 2011-10-12 2013-04-17 Ford Global Tech Llc A method and system for controlling the maximum torque of an internal combustion engine
GB2580979A (en) * 2019-02-04 2020-08-05 Delphi Automotive Systems Lux Method of torque and acceleration monitoring in a vehicle

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61185601A (en) * 1985-02-12 1986-08-19 Mazda Motor Corp Engine control device
JPH0325034A (en) * 1989-06-21 1991-02-01 Toyota Motor Corp Control device for vehicle driving system
JPH0492734A (en) * 1990-08-07 1992-03-25 Mazda Motor Corp Device and method for controlling on-vehicle computer
JPH051614A (en) * 1991-06-24 1993-01-08 Toyota Motor Corp Control device for internal combustion engine
JPH08121304A (en) * 1994-10-26 1996-05-14 Honda Motor Co Ltd Controller for internal combustion engine
JPH10148141A (en) * 1996-11-20 1998-06-02 Toyota Motor Corp Valve timing control device for internal combustion engine
JP2000227033A (en) * 1999-02-05 2000-08-15 Honda Motor Co Ltd Valve system control device for internal combustion engine
JP2001123879A (en) * 1999-10-28 2001-05-08 Denso Corp Combustion state detecting device for internal combustion engine
JP2001129367A (en) * 1999-08-26 2001-05-15 Nitto Denko Corp Spiral type gas-liquid contact membrane module
JP2001164949A (en) * 1999-12-06 2001-06-19 Nissan Motor Co Ltd Control device for variable valve system engine
JP2004036560A (en) * 2002-07-05 2004-02-05 Honda Motor Co Ltd Valve system for internal combustion engine
JP2004239151A (en) * 2003-02-05 2004-08-26 Toyota Motor Corp Control means of internal combustion engine
JP2005248909A (en) * 2004-03-08 2005-09-15 Denso Corp Combustion state detection device of internal combustion engine
JP2008075480A (en) * 2006-09-19 2008-04-03 Denso Corp Control device for engine
JP2008138579A (en) * 2006-12-01 2008-06-19 Denso Corp Variable valve timing control device for internal combustion engine
JP2009185788A (en) * 2008-02-08 2009-08-20 Denso Corp Vehicle control device

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61185601A (en) * 1985-02-12 1986-08-19 Mazda Motor Corp Engine control device
JPH0325034A (en) * 1989-06-21 1991-02-01 Toyota Motor Corp Control device for vehicle driving system
JPH0492734A (en) * 1990-08-07 1992-03-25 Mazda Motor Corp Device and method for controlling on-vehicle computer
JPH051614A (en) * 1991-06-24 1993-01-08 Toyota Motor Corp Control device for internal combustion engine
JPH08121304A (en) * 1994-10-26 1996-05-14 Honda Motor Co Ltd Controller for internal combustion engine
JPH10148141A (en) * 1996-11-20 1998-06-02 Toyota Motor Corp Valve timing control device for internal combustion engine
JP2000227033A (en) * 1999-02-05 2000-08-15 Honda Motor Co Ltd Valve system control device for internal combustion engine
JP2001129367A (en) * 1999-08-26 2001-05-15 Nitto Denko Corp Spiral type gas-liquid contact membrane module
JP2001123879A (en) * 1999-10-28 2001-05-08 Denso Corp Combustion state detecting device for internal combustion engine
JP2001164949A (en) * 1999-12-06 2001-06-19 Nissan Motor Co Ltd Control device for variable valve system engine
JP2004036560A (en) * 2002-07-05 2004-02-05 Honda Motor Co Ltd Valve system for internal combustion engine
JP2004239151A (en) * 2003-02-05 2004-08-26 Toyota Motor Corp Control means of internal combustion engine
JP2005248909A (en) * 2004-03-08 2005-09-15 Denso Corp Combustion state detection device of internal combustion engine
JP2008075480A (en) * 2006-09-19 2008-04-03 Denso Corp Control device for engine
JP2008138579A (en) * 2006-12-01 2008-06-19 Denso Corp Variable valve timing control device for internal combustion engine
JP2009185788A (en) * 2008-02-08 2009-08-20 Denso Corp Vehicle control device

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2495511A (en) * 2011-10-12 2013-04-17 Ford Global Tech Llc A method and system for controlling the maximum torque of an internal combustion engine
CN103047026A (en) * 2011-10-12 2013-04-17 福特环球技术公司 Method for controlling engine of motor vehicle
EP2581591A1 (en) * 2011-10-12 2013-04-17 Ford Global Technologies, LLC A Method for Controlling an Engine of a Motor Vehicle
US9031766B2 (en) 2011-10-12 2015-05-12 Ford Global Technologies, Llc Method for controlling an engine for a motor vehicle
GB2580979A (en) * 2019-02-04 2020-08-05 Delphi Automotive Systems Lux Method of torque and acceleration monitoring in a vehicle

Also Published As

Publication number Publication date
JP5351064B2 (en) 2013-11-27

Similar Documents

Publication Publication Date Title
US7877192B2 (en) Control for an internal-combustion engine
JP3991992B2 (en) Control device for internal combustion engine
JP3960339B2 (en) Intake air quantity variation detector
US8676478B2 (en) Engine stop control device
JP4529190B2 (en) Engine stop control device
JP4066971B2 (en) Control device for internal combustion engine
JP5844162B2 (en) Misfire detection device for internal combustion engine
WO2011086765A1 (en) Misfire detection device for internal combustion engine
JP2005273537A (en) Engine-output control system
JP2004150424A (en) Controlling device of internal combustion engine
JP2012127331A (en) Engine control unit
US10527522B2 (en) Misfire determination device and misfire determination method of internal combustion engine
JP5351064B2 (en) Control device for internal combustion engine
JP5305043B2 (en) Engine combustion state detection device
JP6286044B2 (en) Control device for internal combustion engine
JP2009281239A (en) Control device for internal combustion engine
JP2011157852A (en) Control device of internal combustion engine
US20040254714A1 (en) Engine control device
JP4927697B2 (en) Fuel property estimation device for internal combustion engine
JP4691142B2 (en) Misfire detection device for internal combustion engine
JP2004092603A (en) Combustion state prediction device of internal combustion engine
JP2004346869A (en) Control device of internal combustion engine
JP6271957B2 (en) Engine speed calculation device and engine control device
JP2014234757A (en) Mbt ignition timing computing device and engine controller using the same
JP2003056378A (en) Rotation sensor for crank shaft

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20121128

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130815

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130820

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130822

R150 Certificate of patent or registration of utility model

Ref document number: 5351064

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees