JP6753754B2 - Engine system controller - Google Patents

Engine system controller Download PDF

Info

Publication number
JP6753754B2
JP6753754B2 JP2016197186A JP2016197186A JP6753754B2 JP 6753754 B2 JP6753754 B2 JP 6753754B2 JP 2016197186 A JP2016197186 A JP 2016197186A JP 2016197186 A JP2016197186 A JP 2016197186A JP 6753754 B2 JP6753754 B2 JP 6753754B2
Authority
JP
Japan
Prior art keywords
value
state
deviation
target
control
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.)
Active
Application number
JP2016197186A
Other languages
Japanese (ja)
Other versions
JP2018059444A (en
Inventor
英正 高山
英正 高山
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.)
Hino Motors Ltd
Original Assignee
Hino Motors 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 Hino Motors Ltd filed Critical Hino Motors Ltd
Priority to JP2016197186A priority Critical patent/JP6753754B2/en
Publication of JP2018059444A publication Critical patent/JP2018059444A/en
Application granted granted Critical
Publication of JP6753754B2 publication Critical patent/JP6753754B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

  • Supercharger (AREA)
  • Output Control And Ontrol Of Special Type Engine (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Description

本発明は、ディーゼルエンジン、ディーゼルスロットル、EGR装置、及び、可変容量型ターボチャージャーを備えたエンジンシステムに適用されるエンジンシステムの制御装置に関する。 The present invention relates to a diesel engine, a diesel throttle, an EGR device, and an engine system control device applied to an engine system including a variable displacement turbocharger.

従来から、ディーゼルエンジンを用いたエンジンシステムにおいては、NOxの低減や燃費の向上を図るために様々な技術が用いられている。例えば特許文献1に記載のエンジンシステムでは、ディーゼルエンジンの排気側から吸気側へ排気ガスの一部を還流させる排気再循環(EGR:Exhaust Gas Recirculation)装置やターボチャージャーに流入する排気ガスの流速を変更可能な可変容量型ターボチャージャーが搭載されている。 Conventionally, in an engine system using a diesel engine, various technologies have been used in order to reduce NOx and improve fuel efficiency. For example, in the engine system described in Patent Document 1, the flow velocity of the exhaust gas flowing into the exhaust gas recirculation (EGR: Exhaust Gas Recirculation) device or turbocharger that recirculates a part of the exhaust gas from the exhaust side to the intake side of the diesel engine is determined. It is equipped with a variable capacity turbocharger that can be changed.

特開2003−193875号公報Japanese Unexamined Patent Publication No. 2003-193875

近年では、さらなるNOxの低減や燃費の向上を図るために、より高い精度のもとでエンジンシステムを制御することが求められている。
本発明は、ディーゼルエンジンを備えたエンジンシステムを高い精度のもとで制御することのできるエンジンシステムの制御装置を提供することを目的とする。
In recent years, in order to further reduce NOx and improve fuel efficiency, it is required to control the engine system with higher accuracy.
An object of the present invention is to provide a control device for an engine system capable of controlling an engine system including a diesel engine with high accuracy.

上記課題を解決するエンジンシステムの制御装置は、エンジンシステムの制御対象を制御するエンジンシステムの制御装置であって、前記エンジンシステムは、ディーゼルエンジン、ディーゼルスロットル、EGR装置、及び、可変容量型ターボチャージャーを含むものであり、前記制御対象は、前記ディーゼルスロットル、前記EGR装置のEGR弁、及び、前記可変容量型ターボチャージャーの可変ノズルを含むものであり、前記ディーゼルエンジンの状態量に関する複数のパラメーターの値を取得する取得部と、前記複数のパラメーターの一部である複数の制御パラメーターの各々について、前記取得部の取得値に基づく目標値を演算する目標値演算部と、前記複数の制御パラメーターの各々について前記目標値と前記取得部の取得値との偏差に基づく目標偏差を演算する目標偏差演算部と、前記目標偏差を構成要素に含む偏差ベクトルとゲイン行列とを乗算することで制御対象に対する制御指示値を演算する指示値演算部とを備え、前記制御パラメーターは、前記ディーゼルエンジンが吸入する作動ガスにおけるEGRガスの割合である吸気EGR率と前記可変容量型ターボチャージャーによって過給された吸入空気が流れる通路内の圧力であるブースト圧とを含む。 The control device of the engine system that solves the above problems is a control device of the engine system that controls the controlled object of the engine system, and the engine system is a diesel engine, a diesel throttle, an EGR device, and a variable displacement turbocharger. The controlled object includes the diesel throttle, the EGR valve of the EGR device, and the variable nozzle of the variable displacement turbocharger, and includes a plurality of parameters relating to the state quantity of the diesel engine. A target value calculation unit that calculates a target value based on the acquisition value of the acquisition unit for each of the acquisition unit that acquires the value and the plurality of control parameters that are a part of the plurality of parameters, and the plurality of control parameters. For each, the target deviation calculation unit that calculates the target deviation based on the deviation between the target value and the acquisition value of the acquisition unit, and the deviation vector including the target deviation as a component and the gain matrix are multiplied to control the control target. It includes an instruction value calculation unit that calculates a control instruction value, and the control parameters are an intake EGR rate that is a ratio of EGR gas to a working gas that is sucked by the diesel engine and a suction supercharged by the variable displacement turbocharger. Includes boost pressure, which is the pressure in the passage through which air flows.

上記構成によれば、吸気EGR率及びブースト圧の各々の目標偏差を構成要素に含む偏差ベクトルとゲイン行列とを乗算することにより制御対象の制御指示値が演算される。そして、制御対象には、吸気EGR率及びブースト圧の双方に影響力を持つディーゼルスロットル、EGR装置のEGR弁、及び、可変容量型ターボチャージャーの可変ノズルが設定されている。すなわち、上記構成によれば、吸気EGR率及びブースト圧の各々の目標偏差、及び、吸気EGR率及びブースト圧の各々が制御対象に与える影響力を考慮したうえで制御対象への制御指示値が演算される。これにより、他の制御パラメーターの影響を考慮しつつ、制御パラメーターの各々が目標値となるように制御対象の制御指示値が演算される。その結果、エンジンシステムを高い精度のもとで制御することができる。 According to the above configuration, the control instruction value to be controlled is calculated by multiplying the deviation vector including the target deviations of the intake EGR rate and the boost pressure in the components and the gain matrix. A diesel throttle having an influence on both the intake EGR rate and the boost pressure, an EGR valve of the EGR device, and a variable nozzle of the variable displacement turbocharger are set as control targets. That is, according to the above configuration, the control instruction value for the control target is determined in consideration of the target deviations of the intake EGR rate and the boost pressure, and the influence of each of the intake EGR rate and the boost pressure on the control target. It is calculated. As a result, the control instruction value of the control target is calculated so that each of the control parameters becomes a target value while considering the influence of other control parameters. As a result, the engine system can be controlled with high accuracy.

上記エンジンシステムの制御装置は、前記ディーゼルエンジンの状態量に関するパラメーターであって前記制御対象の制御により値が変化する状態パラメーターについて、前記取得部の取得値のもとで前記ディーゼルエンジンが定常状態にあるときの値である理想値を演算する理想値演算部と、前記状態パラメーターについて前記理想値と前記取得部の取得値との偏差である状態偏差を演算する状態偏差演算部とをさらに備え、前記偏差ベクトルは、前記目標偏差と前記状態偏差とを構成要素に含んでいることが好ましい。 The control device of the engine system puts the diesel engine into a steady state under the acquired value of the acquisition unit for the state parameter which is a parameter related to the state quantity of the diesel engine and whose value changes depending on the control of the controlled object. It further includes an ideal value calculation unit that calculates an ideal value that is a value at a certain time, and a state deviation calculation unit that calculates a state deviation that is a deviation between the ideal value and the acquisition value of the acquisition unit for the state parameter. It is preferable that the deviation vector includes the target deviation and the state deviation as components.

上記構成によれば、制御対象の制御指示値は、制御パラメーターの目標偏差に加えて、状態パラメーターの偏差である状態偏差を考慮して演算される。これにより、制御指示値は、制御パラメーターが目標値となるような値であり、かつ、ディーゼルエンジンが定常状態に近づくような値に設定される。その結果、排気ガスの性能が最もよい定常状態により近い状態にエンジンシステムを制御することができる。 According to the above configuration, the control instruction value of the control target is calculated in consideration of the state deviation, which is the deviation of the state parameter, in addition to the target deviation of the control parameter. As a result, the control instruction value is set to a value such that the control parameter becomes a target value and the diesel engine approaches a steady state. As a result, the engine system can be controlled to a state closer to the steady state with the best exhaust gas performance.

上記エンジンシステムの制御装置において、前記目標偏差演算部は、前記制御パラメーターの各々について偏差の積算値を前記目標偏差として演算することが好ましい。
上記構成によれば、各制御パラメーターの偏差の積算値が目標偏差に設定されるから、各制御パラメーターを目標値に偏差なく到達させるサーボ系を構成することができる。
In the control device of the engine system, it is preferable that the target deviation calculation unit calculates the integrated value of deviations for each of the control parameters as the target deviation.
According to the above configuration, since the integrated value of the deviation of each control parameter is set to the target deviation, it is possible to configure a servo system that causes each control parameter to reach the target value without deviation.

上記エンジンシステムの制御装置は、前記複数のパラメーターの推定値をモデルを用いて演算する状態量推定部をさらに備え、前記取得部は、前記状態量推定部が演算した前記推定値を取得することが好ましい。 The control device of the engine system further includes a state quantity estimation unit that calculates estimated values of the plurality of parameters using a model, and the acquisition unit acquires the estimated values calculated by the state quantity estimation unit. Is preferable.

取得部は、状態量推定部がモデルを用いて演算した推定値をエンジンの状態量に関するパラメーターの値として取得する。そのため、例えば、取得部がエンジンの状態量に関するパラメーターの値としてセンサーの観測値を取得する場合に比べて、より高い精度のもとでエンジンシステムを制御することができる。 The acquisition unit acquires the estimated value calculated by the state quantity estimation unit using the model as the value of the parameter related to the engine state quantity. Therefore, for example, the engine system can be controlled with higher accuracy than when the acquisition unit acquires the observed value of the sensor as the value of the parameter related to the state quantity of the engine.

上記エンジンシステムの制御装置は、前記ディーゼルエンジンの状態量に関するパラメーターを観測する観測部を備え、前記状態量推定部は、前記制御対象への制御指示値を入力ベクトルの構成要素、前記観測部の観測値を観測ベクトルの構成要素として用いるカルマンフィルター理論により前記推定値を演算するものであることが好ましい。そして、前記観測部が観測するパラメーターは、前記ブースト圧、前記ディーゼルエンジンが吸入する空気の量である吸入空気量、インテークマニホールド内の温度であるインマニ温度、これらディーゼルエンジンに対する吸気側の状態量を含み、また、前記可変容量型ターボチャージャーにおけるタービンの回転数であるタービン回転数、すなわちディーゼルエンジンに対する吸気側の状態と排気側の状態とのバランスを示す状態量を含むとよい。 The control device of the engine system includes an observation unit for observing parameters related to the state quantity of the diesel engine, and the state quantity estimation unit inputs a control instruction value to the control target as a component of an input vector of the observation unit. It is preferable that the estimated value is calculated by the Kalman filter theory using the observed value as a component of the observation vector. The parameters observed by the observation unit are the boost pressure, the intake air amount which is the amount of air taken in by the diesel engine, the intake manifold temperature which is the temperature in the intake manifold, and the state amount on the intake side with respect to these diesel engines. It is also preferable to include the turbine rotation speed, which is the rotation speed of the turbine in the variable displacement turbocharger, that is, the state amount indicating the balance between the state on the intake side and the state on the exhaust side with respect to the diesel engine.

上記構成によれば、状態量推定部の演算する推定値を、エンジンの今現在の状態量を示す値にさらに近づけることができる。そして、該推定値を利用して各制御対象の制御指示値が演算されることにより、エンジンシステムをさらに高い精度のもとで制御することができる。 According to the above configuration, the estimated value calculated by the state quantity estimation unit can be made closer to the value indicating the current state quantity of the engine. Then, the engine system can be controlled with higher accuracy by calculating the control instruction value of each control target using the estimated value.

エンジンシステムの制御装置の一実施形態を搭載したエンジンシステムの概略構成を示す図。The figure which shows the schematic structure of the engine system which carries out one Embodiment of the control device of an engine system. カルマンフィルター理論の一例を示すブロック線図。A block diagram showing an example of Kalman filter theory. エンジンシステムの制御装置の一例を示す機能ブロック図。The functional block diagram which shows an example of the control device of an engine system. 制御演算部の一例を示す機能ブロック図。The functional block diagram which shows an example of the control calculation part. 状態量推定部の一例を示す機能ブロック図。The functional block diagram which shows an example of the state quantity estimation part.

図1〜図5を参照して、エンジンシステムの制御装置の一実施形態について説明する。まず、図1を参照してエンジンシステムの全体構成について説明する。
図1に示すように、エンジンシステムは、軽油を燃料とするディーゼルエンジン10(以下、単にエンジン10という。)を備えている。エンジン10のシリンダーブロック11には6つのシリンダー12が形成されている。各シリンダー12においては、吸入した作動ガスに対してインジェクター13から燃料が噴射され、作動ガスと燃料との混合気が燃焼する。こうした混合気の燃焼が所定の順番で各シリンダー12において行われることにより、エンジン10のクランクシャフト10aが駆動される。
An embodiment of the control device of the engine system will be described with reference to FIGS. 1 to 5. First, the overall configuration of the engine system will be described with reference to FIG.
As shown in FIG. 1, the engine system includes a diesel engine 10 (hereinafter, simply referred to as an engine 10) that uses light oil as fuel. Six cylinders 12 are formed in the cylinder block 11 of the engine 10. In each cylinder 12, fuel is injected from the injector 13 with respect to the sucked working gas, and the mixture of the working gas and the fuel is burned. The crankshaft 10a of the engine 10 is driven by the combustion of the air-fuel mixture in each cylinder 12 in a predetermined order.

シリンダーブロック11には、各シリンダー12に作動ガスを供給するためのインテークマニホールド14と、各シリンダー12からの排気ガスが流入するエキゾーストマニホールド15とが接続されている。 An intake manifold 14 for supplying working gas to each cylinder 12 and an exhaust manifold 15 into which exhaust gas from each cylinder 12 flows are connected to the cylinder block 11.

インテークマニホールド14に接続される吸気通路16には、上流側から順に、図示されないエアクリーナー、ターボチャージャー17のコンプレッサー18、インタークーラー19が取り付けられている。また、吸気通路16には、インタークーラー19の下流側であって、かつ、後述するEGR通路25との接続部分よりも上流側に、吸気通路16の流路断面積を変更可能なディーゼルスロットル20(以下、単にスロットル20という。)が取り付けられている。 An air cleaner (not shown), a compressor 18 of a turbocharger 17, and an intercooler 19 are attached to the intake passage 16 connected to the intake manifold 14 in this order from the upstream side. Further, the intake passage 16 has a diesel throttle 20 (a diesel throttle 20) whose flow path cross-sectional area of the intake passage 16 can be changed on the downstream side of the intercooler 19 and on the upstream side of the connection portion with the EGR passage 25 described later. Hereinafter, the throttle 20) is simply attached.

エキゾーストマニホールド15には、排気通路21が接続されている。排気通路21には、コンプレッサー18に連結軸22を介して連結されたタービン23が取り付けられている。また、エキゾーストマニホールド15には、吸気通路16に接続されて排気ガスの一部をEGR(Exhaust Gas Recirculation)ガスとして吸気通路16に導入するEGR装置24のEGR通路25が接続されている。EGR通路25には、EGRクーラー26が取り付けられている。EGR通路25のうちでEGRクーラー26とエキゾーストマニホールド15とを接続する部分は、第1EGR通路25Aと第2EGR通路25Bとで構成される。EGR通路25には、EGRクーラー26の下流側に、EGR通路25の流路断面積を変更可能なEGR弁27が取り付けられている。シリンダー12には、EGR弁27が開状態にあるときに排気ガスと吸入空気との混合気体が作動ガスとして供給され、EGR弁27が閉状態にあるときに吸入空気が作動ガスとして供給される。 An exhaust passage 21 is connected to the exhaust manifold 15. A turbine 23 connected to the compressor 18 via a connecting shaft 22 is attached to the exhaust passage 21. Further, the exhaust manifold 15 is connected to the EGR passage 25 of the EGR device 24 which is connected to the intake passage 16 and introduces a part of the exhaust gas into the intake passage 16 as EGR (Exhaust Gas Recirculation) gas. An EGR cooler 26 is attached to the EGR passage 25. The portion of the EGR passage 25 that connects the EGR cooler 26 and the exhaust manifold 15 is composed of a first EGR passage 25A and a second EGR passage 25B. An EGR valve 27 capable of changing the flow path cross-sectional area of the EGR passage 25 is attached to the EGR passage 25 on the downstream side of the EGR cooler 26. A mixed gas of exhaust gas and intake air is supplied to the cylinder 12 as a working gas when the EGR valve 27 is in the open state, and intake air is supplied as a working gas when the EGR valve 27 is in the closed state. ..

ターボチャージャー17は、タービン23に可変ノズル28が配設された可変容量型ターボチャージャー(VNT:Variable Nozzle Turbo)である。可変ノズル28は、ステッピングモーターを備えたアクチュエーター29の駆動により開度が変更され、タービン23に流入する排気ガスの流路断面積を変更する。 The turbocharger 17 is a variable capacity turbocharger (VNT: Variable Nozzle Turbo) in which a variable nozzle 28 is arranged on a turbine 23. The opening degree of the variable nozzle 28 is changed by driving an actuator 29 provided with a stepping motor, and the cross-sectional area of the flow path of the exhaust gas flowing into the turbine 23 is changed.

エンジンシステムは、吸入空気量センサー31、ブースト圧センサー32、インマニ温度センサー33、エンジン回転数センサー34、タービン回転数センサー35、アクセル開度センサー36を備える。 The engine system includes an intake air amount sensor 31, a boost pressure sensor 32, an in-mani temperature sensor 33, an engine speed sensor 34, a turbine speed sensor 35, and an accelerator opening sensor 36.

吸入空気量センサー31は、コンプレッサー18の上流にて吸入空気の質量流量である吸入空気量Gaを観測する。ブースト圧センサー32は、スロットル20の下流であって、かつ、吸気通路16とEGR通路25との接続部分よりも上流にて、ターボチャージャー17によって過給された吸入空気の圧力であるブースト圧Pbを観測する。インマニ温度センサー33は、インテークマニホールド14内の温度であって、エンジン10が吸入する作動ガスの温度であるインマニ温度Timを観測する。エンジン回転数センサー34は、クランクシャフト10aの回転数であるエンジン回転数Neを観測する。タービン回転数センサー35は、ターボチャージャー17の連結軸22の回転数であるタービン回転数Ntを観測する。アクセル開度センサー36は、運転者が操作するアクセルペダルの踏み込み量であるアクセル開度ACCを観測する。上記各種センサー31〜36は、エンジン10の状態量(運転状態)に関するパラメーターを観測する観測部として機能可能であり、センサー群40(図3参照)を構成する。各種センサー31〜36の出力した信号は、エンジンシステムを統括制御する制御装置であるECU50に入力される。 The intake air amount sensor 31 observes the intake air amount Ga, which is the mass flow rate of the intake air, upstream of the compressor 18. The boost pressure sensor 32 is a boost pressure Pb which is the pressure of the intake air supercharged by the turbocharger 17 downstream of the throttle 20 and upstream of the connection portion between the intake passage 16 and the EGR passage 25. To observe. The intake manifold temperature sensor 33 observes the intake manifold temperature Tim, which is the temperature inside the intake manifold 14, which is the temperature of the working gas sucked by the engine 10. The engine speed sensor 34 observes the engine speed Ne, which is the speed of the crankshaft 10a. The turbine rotation speed sensor 35 observes the turbine rotation speed Nt, which is the rotation speed of the connecting shaft 22 of the turbocharger 17. The accelerator opening sensor 36 observes the accelerator opening ACC, which is the amount of depression of the accelerator pedal operated by the driver. The various sensors 31 to 36 can function as an observation unit for observing parameters related to the state quantity (operating state) of the engine 10, and constitute a sensor group 40 (see FIG. 3). The signals output from the various sensors 31 to 36 are input to the ECU 50, which is a control device that controls the engine system in an integrated manner.

ECU50は、プロセッサ、メモリ、入出力インターフェース等を有するマイクロコンピューターを中心に構成される。ECU50は、各種センサー31〜36の観測値を取得し、その取得した観測値を用いたカルマンフィルター理論により求めた状態量の推定値に基づいてスロットル20、EGR弁27、可変ノズル28といった制御対象60(図3参照)を制御する。 The ECU 50 is mainly composed of a microcomputer having a processor, a memory, an input / output interface, and the like. The ECU 50 acquires the observed values of various sensors 31 to 36, and the control target such as the throttle 20, the EGR valve 27, and the variable nozzle 28 is based on the estimated value of the state quantity obtained by the Kalman filter theory using the acquired observed values. 60 (see FIG. 3) is controlled.

ECU50の構成について説明するまえに、図2を参照してカルマンフィルター理論について説明する。
図2に示すように、カルマンフィルター理論において、制御対象に関する状態量を構成要素とする状態ベクトルをx(k)、制御対象に対する制御指示値を構成要素とする入力ベクトルをu(k)とした場合、状態方程式は、式(1)のように設定される。この入力ベクトルu(k)は、例えば、状態ベクトルx(k)、各種の運転条件におけるエンジンの理想的な状態量の値を構成要素とする理想ベクトルxss(k)、及び、制御対象に関する状態量の目標値を構成要素とする目標ベクトルr(k)に基づくフィードバック制御により演算される。また、予め設定されたセンサーの観測値を構成要素とする観測ベクトルをy(k)とした場合、観測方程式は、式(2)のように設定される。
Before explaining the configuration of the ECU 50, the Kalman filter theory will be described with reference to FIG.
As shown in FIG. 2, in the Kalman filter theory, the state vector whose component is the state quantity related to the controlled object is x (k), and the input vector whose component is the control instruction value for the controlled object is u (k). In this case, the equation of state is set as in Eq. (1). The input vector u (k) is, for example, a state vector x (k), an ideal vector xss (k) having an ideal state quantity value of the engine under various operating conditions as a component, and a state related to a controlled object. It is calculated by feedback control based on the target vector r (k) having the target value of the quantity as a component. Further, when the observation vector whose component is the observation value of the preset sensor is y (k), the observation equation is set as in the equation (2).

なお、kは、0以上の整数であって離散時刻である。また、Aはシステム行列、Bは入力行列、Cは観測行列であり、これらの行列は、例えば運動方程式といった第1原理モデリングやシミュレーションの結果等に基づいて予め設定される。また、w(k)はシステム雑音であり、v(k)は観測雑音である。これらシステム雑音w(k)及び観測雑音v(k)は、正規分布にしたがう互いに独立な白色雑音として取り扱われる。そのため、システム雑音w(k)及び観測雑音v(k)の平均値は0に設定される。また、システム雑音w(k)の分散σ 及び観測雑音v(k)の分散σ は、予め行ったシミュレーション等の結果に基づいて設定される。 Note that k is an integer of 0 or more and is a discrete time. Further, A is a system matrix, B is an input matrix, and C is an observation matrix, and these matrices are preset based on the results of first-principles modeling or simulation such as an equation of motion. Further, w (k) is system noise and v (k) is observation noise. These system noises w (k) and observed noises v (k) are treated as white noises that are independent of each other according to a normal distribution. Therefore, the average values of the system noise w (k) and the observed noise v (k) are set to 0. The dispersibility sigma v 2 of variance sigma w 2 and the observation noise v (k) of the system noise w (k) is set based on the results of simulations that were performed previously.

カルマンフィルター理論では、時刻k−1までの情報に基づいて時刻kにおける状態ベクトルx(k)の推定値である予測推定値x^(k)を演算する予測ステップと、時刻kまでの情報に基づいて時刻kにおける状態ベクトルx(k)の推定値であるフィルタリング推定値x^(k)を演算するフィルタリングステップとが行われる。システム雑音w(k)及び観測雑音v(k)が互いに独立な正規分布にしたがう白色雑音であることから、予測推定値x^(k)は式(3)で示され、フィルタリング推定値x^(k)は式(4)で示される。状態ベクトルx(k)の初期値x(0)は、予め行ったシミュレーションの結果等に基づいて設定される。式(4)において、G(k)はカルマンゲインである。カルマンゲインG(k)は、予め行ったシミュレーション等の結果に基づいて設定された一定値であってもよいし、時刻kごとに更新されてもよい。 In the Kalman filter theory, a prediction step for calculating a predicted estimated value x ^ (k), which is an estimated value of the state vector x (k) at time k, based on information up to time k-1, and information up to time k. The filtering step of calculating the filtering estimated value x ^ (k), which is the estimated value of the state vector x (k) at the time k, is performed based on the above. Since the system noise w (k) and the observed noise v (k) are white noises that follow a normal distribution independent of each other, the predicted estimated value x ^ (k) is expressed by Eq. (3), and the filtered estimated value x. ^ (K) is represented by the equation (4). The initial value x (0) of the state vector x (k) is set based on the result of a simulation performed in advance. In equation (4), G (k) is Kalman gain. The Kalman gain G (k) may be a constant value set based on the result of a simulation or the like performed in advance, or may be updated every time k.

カルマンゲインG(k)が更新される場合、予測ステップにおいて、式(5)に示すように、時刻k−1の事後誤差共分散行列P(k−1)やシステム行列Aに基づいて、時刻kにおける事前誤差共分散行列P(k)が演算される。一方、フィルタリングステップにおいては、式(6)に示すようにカルマンゲインG(k)が時刻kごとに更新される。また、フィルタリングステップにおいては、式(7)に示すように、時刻kにおける事前誤差共分散行列P(k)や更新されたカルマンゲインG(k)等に基づいて事後誤差共分散行列P(k)が演算される。 When the Kalman gain G (k) is updated, in the prediction step, as shown in equation (5), the time is based on the posterior error covariance matrix P (k-1) at time k-1 and the system matrix A. The prior error covariance matrix P (k) at k is calculated. On the other hand, in the filtering step, the Kalman gain G (k) is updated every time k as shown in the equation (6). Further, in the filtering step, as shown in the equation (7), the posterior error covariance matrix P (k) based on the pre-error covariance matrix P (k) at time k, the updated Kalman gain G (k), and the like. k) is calculated.

なお、事前誤差共分散行列P(k)の初期値P(0)は、予め行ったシミュレーション等の結果に基づき設定される。式(5)〜(7)において、Iは単位行列、Aはシステム行列Aの転置行列、Cは観測行列Cの転置行列である。Qは、観測雑音v(k)についての共分散行列であり、予め行ったシミュレーションの結果に基づき設定される。 The initial value P (0) of the prior error covariance matrix P (k) is set based on the result of a simulation or the like performed in advance. In Formula (5) ~ (7), I is a unit matrix, A T is the transposed matrix of the system matrix A, a C T is a transposed matrix of the observation matrix C. Q is a covariance matrix for the observed noise v (k) and is set based on the result of a simulation performed in advance.

すなわち、カルマンフィルター理論では、まず、時刻k−1までの情報に基づいて時刻kにおける状態ベクトルx(k)の推定値である予測推定値x^(k)が演算される。そして、時刻kにおいて観測ベクトルy(k)が入力されると、観測ベクトルy(k)から観測ベクトルy(k)の予測値である観測予測値(=Cx^(k))が減算され、続いて観測予測誤差(=y(k)−Cx^(k))が演算される。 That is, in the Kalman filter theory, first, the predicted estimated value x ^ (k), which is the estimated value of the state vector x (k) at the time k, is calculated based on the information up to the time k-1. Then, when the observation vector y (k) is input at time k, the observation prediction value (= Cx ^ (k)), which is the prediction value of the observation vector y (k), is subtracted from the observation vector y (k). Then, the observation prediction error (= y (k) -Cx ^ - (k)) is calculated.

次に、観測予測誤差に対してカルマンゲインG(k)が乗算されて予測推定値x^(k)についての補正項(=G(k)×(y(k)−Cx^(k)))が演算される。そして、予測推定値x^(k)に補正項が加算され、時刻kまでの情報に基づく状態ベクトルx(k)の推定値であるフィルタリング推定値x^(k)が演算される。 Next, the Kalman gain G (k) is multiplied by the observed prediction error, and the correction term (= G (k) × (y (k) −Cx ^ (k) for the predicted estimated value x ^ (k). ))) Is calculated. Then, the correction term is added to the predicted estimated value x ^ (k), and the filtering estimated value x ^ (k), which is the estimated value of the state vector x (k) based on the information up to the time k, is calculated.

次に、フィルタリング推定値x^(k)に対してシステム行列Aが乗算されたもの(=Ax^(k))と入力ベクトルu(k)に対して入力行列Bが乗算されたもの(=Bu(k))とが加算されることにより、時刻k+1における予測推定値x^(k+1)が演算される。予測推定値x^(k+1)は、時刻を一つ遅らせることにより、次の時刻kにおける予測推定値x^(k)として、入力ベクトルu(k)の演算や観測予測値(=Cx^(k))の演算に用いられる。 Next, the filtering estimate x ^ (k) multiplied by the system matrix A (= Ax ^ (k)) and the input vector u (k) multiplied by the input matrix B (=). By adding Bu (k)), the predicted estimated value x ^ (k + 1) at time k + 1 is calculated. The predicted estimated value x ^ - (k + 1) is set as the predicted estimated value x ^ - (k) at the next time k by delaying the time by one, and the operation of the input vector u (k) or the observed predicted value (= Cx) It is used for the calculation of ^ - (k)).

図3〜図5を参照してECU50の構成について詳しく説明する。
図3に示すように、ECU50は、スロットル20、EGR弁27、可変ノズル28といった制御対象60に対する制御指示値を構成要素とする入力ベクトルu(k)を演算する制御演算部51を備える。またECU50は、各種センサーの観測値を用いたカルマンフィルター理論を適用して、エンジン10の状態量を示すパラメーターを構成要素とする状態ベクトルx(k)の予測推定値x^(k)を演算する状態量推定部52を備える。
The configuration of the ECU 50 will be described in detail with reference to FIGS. 3 to 5.
As shown in FIG. 3, the ECU 50 includes a control calculation unit 51 that calculates an input vector u (k) having a control instruction value for a control target 60 such as a throttle 20, an EGR valve 27, and a variable nozzle 28 as a component. Further, the ECU 50 applies the Kalman filter theory using the observed values of various sensors to obtain a predicted estimated value x ^ (k) of the state vector x (k) having a parameter indicating the state quantity of the engine 10 as a component. A state quantity estimation unit 52 for calculation is provided.

制御演算部51は、取得部として、予測推定値x^(k)、及び、センサー群40の中から予め設定された複数のセンサーの観測値を構成要素とする観測ベクトルy(k)を取得する。制御演算部51は、予測推定値x^(k)及び観測ベクトルy(k)に基づき、状態パラメーターの理想値を構成要素とする理想ベクトルxss(k)と制御パラメーターの目標値を構成要素とする目標ベクトルr(k)とを演算する。制御演算部51は、理想ベクトルxss(k)、目標ベクトルr(k)、及び、予測推定値x^(k)に基づく状態フィードバック制御により入力ベクトルu(k)を演算する。制御演算部51は、入力ベクトルu(k)を制御対象60及び状態量推定部52に出力する。 As an acquisition unit, the control calculation unit 51 receives a predicted estimated value x ^ (k) and an observation vector y (k) whose components are observation values of a plurality of sensors preset from the sensor group 40. get. Based on the predicted estimated value x ^ (k) and the observation vector y (k), the control calculation unit 51 comprises an ideal vector xss (k) having an ideal value of a state parameter as a component and a target value of a control parameter as a component. The target vector r (k) to be calculated is calculated. The control calculation unit 51 calculates the input vector u (k) by state feedback control based on the ideal vector xss (k), the target vector r (k), and the predicted estimated value x ^ (k). The control calculation unit 51 outputs the input vector u (k) to the control target 60 and the state quantity estimation unit 52.

状態量推定部52は、観測ベクトルy(k)、制御演算部51から入力される入力ベクトルu(k)、及び、カルマンフィルター理論を用い、エンジン10の状態量を示す複数のパラメーターを構成要素とする予測推定値x^(k)を演算する。状態量推定部52は、演算した予測推定値x^(k)を制御演算部51に出力する。 The state quantity estimation unit 52 uses the observation vector y (k), the input vector u (k) input from the control calculation unit 51, and the Kalman filter theory, and constitutes a plurality of parameters indicating the state quantity of the engine 10. The predicted estimated value x ^ (k) is calculated. The state quantity estimation unit 52 outputs the calculated predicted estimated value x ^ (k) to the control calculation unit 51.

ECU50において、観測ベクトルy(k)の構成要素である観測パラメーターには、吸入空気量Ga、ブースト圧Pb、インマニ温度Tim、タービン回転数Nt、エンジン回転数Ne、及び、アクセル開度ACCが含まれている。 In the ECU 50, the observation parameters that are components of the observation vector y (k) include the intake air amount Ga, the boost pressure Pb, the intake manifold temperature Tim, the turbine speed Nt, the engine speed Ne, and the accelerator opening ACC. It has been.

状態ベクトルx(k)の構成要素には、目標ベクトルr(k)の構成要素である制御パラメーターとしての吸気EGR率ηi及びブースト圧Pbが含まれている。また、状態ベクトルx(k)の構成要素には、吸入空気量Ga、インマニ温度Tim、タービン回転数Ntが含まれている。吸気EGR率ηiは、作動ガス量Gwg(重量)におけるEGRガス量Gr(重量)の割合(=Gr/Gwg)を示す。また、状態ベクトルx(k)の構成要素には、インテークマニホールド14における作動ガスの密度である吸気密度ρim、エキゾーストマニホールド15における排気ガスの圧力である排気圧力Pem、エキゾーストマニホールド15における排気ガスの密度である排気密度ρem、及び、エンジン10が排出した排気ガス量Gex(重量)におけるEGRガス量Gr(重量)の割合(=Gr/Gex)である排気EGR率ηe等が含まれている。 The components of the state vector x (k) include an intake EGR rate ηi and a boost pressure Pb as control parameters that are components of the target vector r (k). Further, the components of the state vector x (k) include the intake air amount Ga, the intake manifold temperature Tim, and the turbine rotation speed Nt. The intake EGR rate ηi indicates the ratio (= Gr / Gwg) of the EGR gas amount Gr (weight) to the working gas amount Gwg (weight). Further, the components of the state vector x (k) include an intake density ρim, which is the density of the working gas in the intake manifold 14, an exhaust pressure Pem, which is the pressure of the exhaust gas in the exhaust manifold 15, and the density of the exhaust gas in the exhaust manifold 15. The exhaust density ρem, which is, and the exhaust EGR ratio ηe, which is the ratio (= Gr / Gex) of the EGR gas amount Gr (weight) to the exhaust gas amount Gex (weight) discharged by the engine 10, and the like are included.

入力ベクトルu(k)の構成要素には、スロットル20の制御指示値であるスロットル開度Ath、EGR弁27の制御指示値であるEGR弁開度Ar、可変ノズル28の制御指示値であるノズル開度Anzが含まれている。 The components of the input vector u (k) include a throttle opening Ath which is a control instruction value of the throttle 20, an EGR valve opening Ar which is a control instruction value of the EGR valve 27, and a nozzle which is a control instruction value of the variable nozzle 28. The opening degree Anz is included.

上述した状態ベクトルx(k)の構成要素には、状態パラメーターが含まれている。状態パラメーターは、本実施形態では、ブースト圧Pb、吸気密度ρim、吸気EGR率ηi、排気圧力Pem、排気密度ρem、排気EGR率ηe等である。このように状態パラメーターは、スロットル20、EGR弁27、可変ノズル28といった制御対象60の制御を通じてエンジン10の運転状態が変化すると値が変化するパラメーターである。また、制御パラメーターは、制御対象60を制御することでエンジン10を目標の運転状態へと導くためのパラメーターである。 The components of the state vector x (k) described above include state parameters. In the present embodiment, the state parameters are boost pressure Pb, intake density ρim, intake EGR rate ηi, exhaust pressure Pem, exhaust density ρem, exhaust EGR rate ηe, and the like. As described above, the state parameter is a parameter whose value changes when the operating state of the engine 10 changes through the control of the controlled object 60 such as the throttle 20, the EGR valve 27, and the variable nozzle 28. Further, the control parameter is a parameter for guiding the engine 10 to a target operating state by controlling the control target 60.

図4に示すように、制御演算部51は、状態量設定部51a、状態偏差演算部51b、目標偏差演算部51c、及び、指示値演算部51fを備えている。
状態量設定部51aは、理想値演算部として各状態パラメーターの理想値を構成要素とする理想ベクトルxss(k)を設定する。理想値は、各運転条件で排気ガスの性能が最もよい定常状態にエンジン10がある場合に各状態パラメーターがとるはずの値である。
As shown in FIG. 4, the control calculation unit 51 includes a state quantity setting unit 51a, a state deviation calculation unit 51b, a target deviation calculation unit 51c, and an indicated value calculation unit 51f.
The state quantity setting unit 51a sets an ideal vector xss (k) having an ideal value of each state parameter as a component as an ideal value calculation unit. The ideal value is a value that each state parameter should take when the engine 10 is in a steady state with the best exhaust gas performance under each operating condition.

例えば、状態量設定部51aは、エンジン回転数Neとインジェクター13への制御指示値である燃料噴射量Gfとに基づいて理想ベクトルxss(k)を演算する。燃料噴射量Gfは、エンジン回転数Neとアクセル開度ACCとに基づく要求トルクTrに基づいて演算され、入力ベクトルu(k)とは別個にインジェクター13に入力されてもよいし、入力ベクトルu(k)の一つとしてインジェクター13に入力されてもよい。入力ベクトルu(k)に燃料噴射量Gfが含まれる場合、状態量設定部51aが演算時に用いる燃料噴射量Gfは、前回の制御指示値である。状態量設定部51aは、予め行った実験やシミュレーションの結果に基づき、エンジン回転数Neと燃料噴射量Gfとに応じた理想値が規定された理想値マップを状態パラメーターごとに保持している。そして状態量設定部51aは、エンジン回転数Neと燃料噴射量Gfとに応じた理想値を各状態パラメーターの理想値マップから選択することにより理想ベクトルxss(k)を演算する。 For example, the state quantity setting unit 51a calculates the ideal vector xss (k) based on the engine speed Ne and the fuel injection amount Gf which is a control instruction value to the injector 13. The fuel injection amount Gf is calculated based on the required torque Tr based on the engine speed Ne and the accelerator opening ACC, and may be input to the injector 13 separately from the input vector u (k), or the input vector u. It may be input to the injector 13 as one of (k). When the input vector u (k) includes the fuel injection amount Gf, the fuel injection amount Gf used by the state amount setting unit 51a at the time of calculation is the previous control instruction value. The state quantity setting unit 51a holds an ideal value map for each state parameter in which ideal values corresponding to the engine speed Ne and the fuel injection amount Gf are defined based on the results of experiments and simulations performed in advance. Then, the state amount setting unit 51a calculates the ideal vector xss (k) by selecting the ideal value corresponding to the engine speed Ne and the fuel injection amount Gf from the ideal value map of each state parameter.

状態量設定部51aは、目標値演算部として各制御パラメーターの目標値を構成要素とする目標ベクトルr(k)を設定する。目標値は、エンジン10を目標の運転状態へ移行させるうえで各制御パラメーターがとるべき値である。 The state quantity setting unit 51a sets a target vector r (k) having a target value of each control parameter as a component as a target value calculation unit. The target value is a value that each control parameter should take in order to shift the engine 10 to the target operating state.

例えば、状態量設定部51aは、エンジン回転数Neと燃料噴射量Gfとに基づいて目標ベクトルr(k)を演算する。状態量設定部51aは、制御パラメーターごとに目標値を演算することにより目標ベクトルr(k)を演算する。状態量設定部51aは、予め行った実験やシミュレーションの結果に基づき、エンジン回転数Neと燃料噴射量Gfとに応じた目標値が規定された目標値マップを制御パラメーターごとに保持している。そして状態量設定部51aは、エンジン回転数Neと燃料噴射量Gfとに応じた目標値を各制御パラメーターの目標値マップから選択することにより目標ベクトルr(k)を演算する。 For example, the state quantity setting unit 51a calculates the target vector r (k) based on the engine speed Ne and the fuel injection amount Gf. The state quantity setting unit 51a calculates the target vector r (k) by calculating the target value for each control parameter. The state quantity setting unit 51a holds a target value map for each control parameter in which target values corresponding to the engine speed Ne and the fuel injection amount Gf are defined based on the results of experiments and simulations performed in advance. Then, the state amount setting unit 51a calculates the target vector r (k) by selecting the target value corresponding to the engine speed Ne and the fuel injection amount Gf from the target value map of each control parameter.

状態偏差演算部51bは、減算器であり、理想ベクトルxss(k)から予測推定値x^(k)を減算することにより、状態パラメーターの偏差を構成要素とする状態偏差ベクトルxe(k)を演算する。 The state deviation calculation unit 51b is a subtractor, and the state deviation vector xe (k) having the deviation of the state parameter as a component by subtracting the predicted estimated value x ^ (k) from the ideal vector xss (k). Is calculated.

目標偏差演算部51cは、減算器51dと積分回路51eとを有している。減算器51dは、目標ベクトルr(k)から予測推定値x^(k)を減算することにより、各制御パラメーターの偏差を構成要素とする制御偏差ベクトルre(k)を演算する。積分回路51eは、制御偏差ベクトルre(k)を積算することにより、制御パラメーターの偏差の積算値を構成要素とする目標偏差ベクトルxi(k)を演算する。 The target deviation calculation unit 51c has a subtractor 51d and an integration circuit 51e. The subtractor 51d calculates the control deviation vector re (k) having the deviation of each control parameter as a component by subtracting the predicted estimated value x ^ (k) from the target vector r (k). By integrating the control deviation vector re (k), the integration circuit 51e calculates the target deviation vector xi (k) having the integrated value of the deviation of the control parameters as a component.

指示値演算部51fは、ゲイン行列Kと、状態偏差ベクトルxe(k)と目標偏差ベクトルxi(k)とで構成される偏差ベクトルe(k)(=[xe(k) xi(k)])とを乗算することにより、入力ベクトルu(k)を構成するスロットル開度Ath、EGR弁開度Ar、ノズル開度Anzを演算する。このゲイン行列Kは、予め行った実験の結果やシミュレーションの結果に基づいて設定されるものであり、例えば最小二乗推定法等に基づいて設定される。 The indicated value calculation unit 51f has a deviation vector e (k) (= [xe (k) xi (k)] composed of a gain matrix K, a state deviation vector xi (k), and a target deviation vector xi (k). By multiplying with T ), the throttle opening Ath, the EGR valve opening Ar, and the nozzle opening Anz constituting the input vector u (k) are calculated. This gain matrix K is set based on the results of experiments performed in advance and the results of simulations, and is set based on, for example, the least squares estimation method.

すなわち、ECU50において、制御演算部51は、状態パラメーターについての偏差と制御パラメーターについての偏差とに基づいて入力ベクトルu(k)を演算し、その演算した入力ベクトルu(k)を制御対象60及び状態量推定部52に出力する。 That is, in the ECU 50, the control calculation unit 51 calculates the input vector u (k) based on the deviation of the state parameter and the deviation of the control parameter, and the calculated input vector u (k) is the control target 60 and Output to the state quantity estimation unit 52.

図5に示すように、状態量推定部52は、エンジン10の1サイクルの平均値を扱うエンジンモデル52aと、観測予測誤差を演算する誤差演算器52bと、カルマンフィルター理論による補正項を演算する補正項演算部52cと、を備える。 As shown in FIG. 5, the state quantity estimation unit 52 calculates the engine model 52a that handles the average value of one cycle of the engine 10, the error calculator 52b that calculates the observation prediction error, and the correction term by the Kalman filter theory. A correction term calculation unit 52c is provided.

エンジンモデル52aは、入力ベクトルu(k)、システム行列A、入力行列B、予測推定値x^(k)、及び、補正項演算部52cが演算した補正項に基づき、予測推定値x^(k+1)を演算する。エンジンモデル52aは、予測推定値x^(k+1)の時刻を一つ遅らせることにより、当該予測推定値x^(k+1)を次の時刻kの予測推定値x^(k)として取り扱う。また、エンジンモデル52aは、予測推定値x^(k)と観測行列Cとに基づき、観測予測値(=Cx^(k))を演算する。 The engine model 52a is based on the input vector u (k), the system matrix A, the input matrix B, the predicted estimated value x ^ (k), and the correction term calculated by the correction term calculation unit 52c, and the predicted estimated value x ^. - (K + 1) is calculated. Engine model 52a is predicted estimate x ^ - by delaying one time (k + 1), the predicted estimates x ^ - handled as (k) - a (k + 1) predicted estimate x ^ at the next time k .. Further, the engine model 52a calculates the observed predicted value (= Cx ^ (k)) based on the predicted estimated value x ^ (k) and the observation matrix C.

誤差演算器52bは、観測ベクトルy(k)から観測予測値を減算することにより観測予測誤差(=y(k)−Cx^(k))を演算する。補正項演算部52cは、観測予測誤差にカルマンゲインG(k)を乗算することにより、予測推定値x^(k)についての補正項を演算し、その演算した補正項をエンジンモデル52aに出力する。 The error calculator 52b calculates the observation prediction error (= y (k) -Cx ^ - (k)) by subtracting the observation prediction value from the observation vector y (k). The correction term calculation unit 52c calculates the correction term for the predicted estimated value x ^ (k) by multiplying the observation prediction error by the Kalman gain G (k), and the calculated correction term is applied to the engine model 52a. Output.

上述した構成のECU50の作用について説明する。
ECU50は、状態パラメーターについて、予測推定値x^(k)、理想ベクトルxss(k)、及び、理想ベクトルxss(k)と予測推定値x^(k)との偏差で構成される状態偏差ベクトルxe(k)、これらを演算する。またECU50は、制御パラメーターについて、予測推定値x^(k)、目標ベクトルr(k)、及び、目標ベクトルr(k)と予測推定値x^(k)との偏差の積算値で構成される目標偏差ベクトルxi(k)、これらを演算する。そしてECU50は、ゲイン行列Kと、状態偏差ベクトルxe(k)と目標偏差ベクトルxi(k)とを構成要素とする偏差ベクトルe(k)とを乗算することで入力ベクトルu(k)を演算する。
The operation of the ECU 50 having the above configuration will be described.
The ECU 50 is a state in which the state parameter is composed of a predicted estimated value x ^ (k), an ideal vector xss (k), and a deviation between the ideal vector xss (k) and the predicted estimated value x ^ (k). Deviation vector xe (k), these are calculated. Further, the ECU 50 uses the integrated value of the predicted estimated value x ^ (k), the target vector r (k), and the deviation between the target vector r (k) and the predicted estimated value x ^ (k) for the control parameters. The target deviation vector xi (k) to be constructed is calculated. Then, the ECU 50 calculates the input vector u (k) by multiplying the gain matrix K by the deviation vector e (k) whose constituent elements are the state deviation vector xe (k) and the target deviation vector xi (k). To do.

このようにECU50では、スロットル開度Ath、EGR弁開度Ar、ノズル開度Anzに密接に関わり合う制御パラメーターである吸気EGR率ηi及びブースト圧Pbによって目標ベクトルr(k)が構成されている。そして、ゲイン行列Kと各制御パラメーターの偏差の積算値である積算偏差を含む偏差ベクトルe(k)との乗算によって入力ベクトルu(k)が演算される。そのため、入力ベクトルu(k)を構成するスロットル開度Ath、EGR弁開度Ar、ノズル開度Anzを演算する演算式の各々は、吸気EGR率ηiの積算偏差を含む項とブースト圧Pbの積算偏差を含む項とを有する。すなわち、スロットル開度Ath、EGR弁開度Ar、ノズル開度Anzの各々は、吸気EGR率ηiの積算偏差とブースト圧Pbの積算偏差とを考慮したうえで、これら吸気EGR率ηi及びブースト圧Pbの双方が目標値に近づく値に設定される。 As described above, in the ECU 50, the target vector r (k) is configured by the intake EGR rate ηi and the boost pressure Pb, which are control parameters closely related to the throttle opening Ath, the EGR valve opening Ar, and the nozzle opening Anz. .. Then, the input vector u (k) is calculated by multiplying the gain matrix K with the deviation vector e (k) including the integrated deviation which is the integrated value of the deviations of each control parameter. Therefore, each of the arithmetic expressions for calculating the throttle opening Ath, the EGR valve opening Ar, and the nozzle opening Anz constituting the input vector u (k) includes the term including the integrated deviation of the intake EGR rate ηi and the boost pressure Pb. It has a term including an integrated deviation. That is, each of the throttle opening Ath, the EGR valve opening Ar, and the nozzle opening Anz takes into consideration the integrated deviation of the intake EGR rate ηi and the integrated deviation of the boost pressure Pb, and these intake EGR rate ηi and boost pressure. Both Pb are set to values close to the target value.

その結果、例えば、ブースト圧Pbに基づくスロットル開度Athのフィードバック制御と吸気EGR率ηiに基づくEGR弁開度Ar及びノズル開度Anzのフィードバック制御とが互いに独立して行われる場合よりも高い精度のもとでエンジンシステムを制御することができる。 As a result, for example, the feedback control of the throttle opening Ath based on the boost pressure Pb and the feedback control of the EGR valve opening Ar and the nozzle opening Anz based on the intake EGR rate ηi are performed with higher accuracy than when they are performed independently of each other. The engine system can be controlled under.

また、偏差ベクトルe(k)が各状態パラメーターの偏差を構成要素とする状態偏差ベクトルxe(k)を含んでいる。そのため、開度Ath、Ar,Anzの各々を求める演算式は、状態パラメーターの偏差を含む項を状態パラメーターごとに有している。すなわち、開度Ath、Ar,Anzの各々は、各状態パラメーターの偏差を考慮して、各状態パラメーターの値を理想値に近づける値に設定される。つまり、各開度Ath、Ar,Anzは、制御パラメーターである吸気EGR率ηi及びブースト圧Pbの双方が目標値に近づく値であり、かつ、各状態パラメーターが理想値に近づく値に設定される。 Further, the deviation vector e (k) includes the state deviation vector xe (k) having the deviation of each state parameter as a component. Therefore, the arithmetic expression for obtaining each of the opening degree Ath, Ar, and Anz has a term including the deviation of the state parameter for each state parameter. That is, each of the opening degree Ath, Ar, and Anz is set to a value that brings the value of each state parameter closer to the ideal value in consideration of the deviation of each state parameter. That is, each opening degree Ath, Ar, and Anz is set to a value in which both the intake EGR rate ηi and the boost pressure Pb, which are control parameters, approach the target value, and each state parameter approaches the ideal value. ..

上記実施形態のECU50によれば、以下に列挙する効果が得られる。
(1)開度Ath、Ar,Anzの各々が制御パラメーターである吸気EGR率ηi及びブースト圧Pbの双方を目標値に近づける値に設定される。これにより、高い精度のもとでエンジンシステムを制御することができる。
According to the ECU 50 of the above embodiment, the effects listed below can be obtained.
(1) Each of the opening degree Ath, Ar, and Anz is set to a value that brings both the intake EGR rate ηi and the boost pressure Pb, which are control parameters, close to the target value. This makes it possible to control the engine system with high accuracy.

(2)また、開度Ath、Ar,Anzの各々は、吸気EGR率ηi及びブースト圧Pbの双方を目標値に近づけつつ、かつ、状態パラメーターの各々を理想値に近づける値に設定される。これにより、定常状態により近い状態でエンジン10を運転させることができる。その結果、排気ガスの性能の向上を図ることができる。 (2) Further, each of the opening degree Ath, Ar, and Anz is set to a value that brings both the intake EGR rate ηi and the boost pressure Pb close to the target value and each of the state parameters close to the ideal value. As a result, the engine 10 can be operated in a state closer to the steady state. As a result, the performance of the exhaust gas can be improved.

(3)目標偏差演算部51cは、制御偏差ベクトルre(k)を積算することで目標偏差ベクトルxi(k)を演算する。すなわち、目標偏差演算部51cは、吸気EGR率ηi及びブースト圧Pbの各々について偏差の積算値を目標偏差に設定する。これにより、各制御パラメーターを目標値に偏差なく到達させるサーボ系を構成することができる。 (3) The target deviation calculation unit 51c calculates the target deviation vector xi (k) by integrating the control deviation vector re (k). That is, the target deviation calculation unit 51c sets the integrated value of the deviation for each of the intake EGR rate ηi and the boost pressure Pb as the target deviation. This makes it possible to configure a servo system that allows each control parameter to reach the target value without deviation.

(4)状態量推定部52は、エンジンモデル52aが演算した時刻kの予測推定値x^(k)に基づいてカルマンフィルター理論による補正項を演算し、その補正項を用いて時刻k+1における予測推定値x^(k+1)を演算する。このようにカルマンフィルター理論を用いて予測推定値x^(k)を演算することで予測推定値x^(k)の精度が高まることから、エンジンシステムをより高い精度のもとで制御することができる。 (4) The state quantity estimation unit 52 calculates a correction term based on the Kalman filter theory based on the predicted estimated value x ^ (k) of the time k calculated by the engine model 52a, and uses the correction term at the time k + 1. Calculate the predicted estimated value x ^ (k + 1). By calculating the predicted estimated value x ^ - (k) using the Kalman filter theory in this way, the accuracy of the predicted estimated value x ^ - (k) is improved, so that the engine system is controlled with higher accuracy. can do.

(5)タービン回転数Ntは、吸気側の状態と排気側の状態とのバランスを示す1つの指標である。そのため、エンジン10の吸気側の状態を示す吸入空気量Ga、ブースト圧Pb、インマニ温度Timに加えて、タービン回転数Ntを観測値の1つとしてカルマンフィルター理論に適用することにより、予測推定値x^(k)の精度がさらに高まる。 (5) The turbine speed Nt is one index showing the balance between the state on the intake side and the state on the exhaust side. Therefore, by applying the turbine rotation speed Nt as one of the observed values to the Kalman filter theory in addition to the intake air amount Ga, boost pressure Pb, and intake manifold temperature Tim, which indicate the state of the intake side of the engine 10, the predicted estimated value The accuracy of x ^ - (k) is further improved.

なお、上記実施形態は、以下のように適宜変更して実施することもできる。
・状態量推定部52が各パラメーターの推定値を演算する方法としては、カルマンフィルター理論を用いた方法に限られない。例えば、予め行ったシミュレーションや実験の結果に基づいてエンジンモデル52aを構築し、このエンジンモデル52aに観測ベクトルy(k)及び入力ベクトルu(k)を適用することにより各パラメーターが演算される方法であってもよい。また例えば、予め行ったシミュレーションや実験の結果に基づいて構築された統計モデルに対して、観測ベクトルy(k)や入力ベクトルu(k)を適用することにより各パラメーターが演算される方法であってもよい。
The above embodiment can be modified as appropriate and implemented as follows.
-The method by which the state quantity estimation unit 52 calculates the estimated value of each parameter is not limited to the method using the Kalman filter theory. For example, a method in which an engine model 52a is constructed based on the results of simulations and experiments performed in advance, and each parameter is calculated by applying an observation vector y (k) and an input vector u (k) to the engine model 52a. It may be. Further, for example, it is a method in which each parameter is calculated by applying the observation vector y (k) and the input vector u (k) to a statistical model constructed based on the results of simulations and experiments performed in advance. You may.

・ECU50は、状態量推定部52が割愛された構成であってもよい。この場合、状態量設定部51aは、観測ベクトルy(k)の構成要素に基づいて、制御パラメーター及び状態パラメーターの各々の現在値及び目標値を演算する。こうした構成であっても上記(1)〜(3)に準ずる効果を得ることができる。 The ECU 50 may have a configuration in which the state quantity estimation unit 52 is omitted. In this case, the state quantity setting unit 51a calculates the current value and the target value of each of the control parameter and the state parameter based on the components of the observation vector y (k). Even with such a configuration, an effect similar to the above (1) to (3) can be obtained.

例えば、制御パラメーターである吸気EGR率ηiについては次のように構成することが可能である。すなわち、状態量設定部51aは、インテークマニホールド14における作動ガスの圧力であるインマニ圧Pim、インマニ温度Tim、エンジン回転数Ne、エンジン10の排気量D等を用いた状態方程式に基づいて作動ガス量Gwgを演算し、この作動ガス量Gwgから吸入空気量Gaを減算することによりEGRガス量Grを演算する。そして状態量設定部51aは、作動ガス量GwgとEGRガス量Grとに基づいて吸気EGR率ηiの演算値を演算する。また、状態量設定部51aは、エンジン回転数Neとアクセル開度ACCとに基づいて吸気EGR率ηiの目標値を演算する。そして、減算器51dは、これら吸気EGR率ηiの目標値と現在値との偏差を演算する。なお、EGRガス量Grがセンサーにより観測され、観測値そのものに基づいて吸気EGR率ηi(=Gr/(Ga+Gr))の現在値が求められてもよい。 For example, the intake EGR rate ηi, which is a control parameter, can be configured as follows. That is, the state amount setting unit 51a is based on a state equation using the intake manifold pressure Pim, the intake manifold temperature Tim, the engine rotation speed Ne, the exhaust amount D of the engine 10, and the like, which are the pressures of the working gas in the intake manifold 14. Gwg is calculated, and the EGR gas amount Gr is calculated by subtracting the intake air amount Ga from the working gas amount Gwg. Then, the state amount setting unit 51a calculates the calculated value of the intake EGR rate ηi based on the working gas amount Gwg and the EGR gas amount Gr. Further, the state amount setting unit 51a calculates a target value of the intake EGR rate ηi based on the engine speed Ne and the accelerator opening degree ACC. Then, the subtractor 51d calculates the deviation between the target value and the current value of the intake EGR rate ηi. The EGR gas amount Gr may be observed by the sensor, and the current value of the intake EGR rate ηi (= Gr / (Ga + Gr)) may be obtained based on the observed value itself.

・ECU50において、目標偏差演算部51cは、積分回路51eが割愛された構成であってもよい。こうした構成であっても、制御偏差ベクトルre(k)を目標偏差ベクトルxi(k)として演算することができる。また、目標偏差演算部51cは、制御偏差ベクトルre(k)を微分する微分回路を有していてもよい。こうした場合、目標偏差ベクトルxi(k)は、偏差の微分値を構成要素として有する。 In the ECU 50, the target deviation calculation unit 51c may have a configuration in which the integrator circuit 51e is omitted. Even with such a configuration, the control deviation vector re (k) can be calculated as the target deviation vector xi (k). Further, the target deviation calculation unit 51c may have a differentiating circuit for differentiating the control deviation vector re (k). In such a case, the target deviation vector xi (k) has a differential value of the deviation as a component.

・ECU50において、状態パラメーターの理想値の求め方は、エンジン回転数Ne、燃料噴射量Gf、及び、理想値マップを用いる方法に限られない。例えば、理想値は、観測ベクトルy(k)や予測推定値x^(k)から選択される構成要素の値を所定の演算式に代入することにより演算されてもよい。 -In the ECU 50, the method of obtaining the ideal value of the state parameter is not limited to the method using the engine speed Ne, the fuel injection amount Gf, and the ideal value map. For example, the ideal value may be calculated by substituting the value of the component selected from the observation vector y (k) or the predicted estimated value x ^ (k) into a predetermined calculation formula.

・ECU50において、制御パラメーターの目標値の求め方は、エンジン回転数Ne、燃料噴射量Gf、及び、目標値マップを用いる方法に限られない。例えば、目標値は、観測ベクトルy(k)や予測推定値x^(k)から選択される構成要素の値を所定の演算式に代入することにより演算されてもよい。 -In the ECU 50, the method of obtaining the target value of the control parameter is not limited to the method using the engine speed Ne, the fuel injection amount Gf, and the target value map. For example, the target value may be calculated by substituting the value of the component selected from the observation vector y (k) or the predicted estimated value x ^ (k) into a predetermined calculation formula.

・ECU50において、制御演算部51は、ゲイン行列Kと目標偏差ベクトルxi(k)とに基づく状態フィードバック制御によって制御指示値を演算してもよい。すなわち、制御演算部51は、状態パラメーターを用いることなく、制御パラメーターのみで状態フィードバック制御を行ってもよい。 In the ECU 50, the control calculation unit 51 may calculate the control instruction value by the state feedback control based on the gain matrix K and the target deviation vector xi (k). That is, the control calculation unit 51 may perform state feedback control using only the control parameters without using the state parameters.

・ブースト圧Pbは、ターボチャージャー17によって過給された吸入空気が流れる通路内の圧力であればよく、上記インマニ圧Pimをブースト圧Pbとしてもよい。
・制御パラメーターは、吸気EGR率ηiとブースト圧Pbだけに限られるものではなく、これら吸気EGR率ηiとブースト圧Pbに加えて、例えば吸入空気量Gaといった上述した状態パラメーターとして記載したパラメーターが含まれていてもよい。
The boost pressure Pb may be any pressure in the passage through which the intake air supercharged by the turbocharger 17 flows, and the intake manifold pressure Pim may be used as the boost pressure Pb.
The control parameters are not limited to the intake EGR rate ηi and the boost pressure Pb, but include the parameters described as the above-mentioned state parameters such as the intake air amount Ga in addition to the intake EGR rate ηi and the boost pressure Pb. It may be.

・状態パラメーターは、制御対象60の制御によって値が変化するパラメーターであればよい。そのため、上記実施形態に状態パラメーターとして記載したパラメーターの一部であってもよいし、上記実施形態に記載した状態パラメーターに加えて他のパラメーターが含まれていてもよい。 -The state parameter may be any parameter whose value changes depending on the control of the control target 60. Therefore, it may be a part of the parameters described as the state parameters in the above embodiment, or may include other parameters in addition to the state parameters described in the above embodiment.

10…ディーゼルエンジン、10a…クランクシャフト、11…シリンダーブロック、12…シリンダー、13…インジェクター、14…インテークマニホールド、15…エキゾーストマニホールド、16…吸気通路、17…ターボチャージャー、18…コンプレッサー、19…インタークーラー、20…ディーゼルスロットル、21…排気通路、22…連結軸、23…タービン、24…EGR装置、25…EGR通路、26…EGRクーラー、27…EGR弁、28…可変ノズル、29…アクチュエーター、31…吸入空気量センサー、32…ブースト圧センサー、33…インマニ温度センサー、34…エンジン回転数センサー、35…タービン回転数センサー、36…アクセル開度センサー、40…センサー群、50…ECU、51…制御演算部、51a…状態量設定部、51b…状態偏差演算部、51c…目標偏差演算部、51d…減算器、51e…積分回路、51f…指示値演算部、52…状態量推定部、52a…エンジンモデル、52b…誤差演算器、52c…補正項演算部、60…制御対象。 10 ... Diesel engine, 10a ... Crank shaft, 11 ... Cylinder block, 12 ... Cylinder, 13 ... Injector, 14 ... Intake manifold, 15 ... Exhaust manifold, 16 ... Intake passage, 17 ... Turbocharger, 18 ... Compressor, 19 ... Intercooler , 20 ... diesel throttle, 21 ... exhaust passage, 22 ... connecting shaft, 23 ... turbine, 24 ... EGR device, 25 ... EGR passage, 26 ... EGR cooler, 27 ... EGR valve, 28 ... variable nozzle, 29 ... actuator, 31 ... intake air volume sensor, 32 ... boost pressure sensor, 33 ... intake manifold temperature sensor, 34 ... engine rotation speed sensor, 35 ... turbine rotation speed sensor, 36 ... accelerator opening sensor, 40 ... sensor group, 50 ... ECU, 51 ... Control calculation unit, 51a ... state quantity setting unit, 51b ... state deviation calculation unit, 51c ... target deviation calculation unit, 51d ... subtractor, 51e ... integration circuit, 51f ... instruction value calculation unit, 52 ... state quantity estimation unit, 52a ... Engine model, 52b ... Error calculator, 52c ... Correction term calculation unit, 60 ... Control target.

Claims (4)

エンジンシステムの制御対象を制御するエンジンシステムの制御装置であって、
前記エンジンシステムは、ディーゼルエンジン、ディーゼルスロットル、EGR装置、及び、可変容量型ターボチャージャーを含むものであり、
前記制御対象は、前記ディーゼルスロットル、前記EGR装置のEGR弁、及び、前記可変容量型ターボチャージャーの可変ノズルを含むものであり、
前記ディーゼルエンジンの状態量に関する複数のパラメーターの値を取得する取得部と、
前記複数のパラメーターの一部である複数の制御パラメーターの各々について、前記取得部の取得値に基づく目標値を演算する目標値演算部と、
前記複数の制御パラメーターの各々について前記目標値と前記取得部の取得値との偏差に基づく目標偏差を演算する目標偏差演算部と、
前記ディーゼルエンジンの状態量に関するパラメーターであって制御対象の制御により値が変化する状態パラメーターについて、前記取得部の取得値のもとで前記ディーゼルエンジンが定常状態にあるときの値である理想値を演算する理想値演算部と、
前記状態パラメーターについて前記理想値と前記取得部の取得値との偏差である状態偏差を演算する状態偏差演算部と、
前記目標偏差および前記状態偏差を構成要素に含む偏差ベクトルとゲイン行列とを乗算することで前記制御対象に対する制御指示値を演算する指示値演算部とを備え、
前記制御パラメーターは、前記ディーゼルエンジンが吸入する作動ガスにおけるEGRガスの割合である吸気EGR率と前記可変容量型ターボチャージャーによって過給された吸入空気が流れる通路内の圧力であるブースト圧とを含む
エンジンシステムの制御装置。
An engine system control device that controls the control target of the engine system.
The engine system includes a diesel engine, a diesel throttle, an EGR device, and a variable displacement turbocharger.
The controlled object includes the diesel throttle, the EGR valve of the EGR device, and the variable nozzle of the variable capacity turbocharger.
An acquisition unit that acquires the values of a plurality of parameters related to the state quantity of the diesel engine, and
For each of the plurality of control parameters that are a part of the plurality of parameters, a target value calculation unit that calculates a target value based on the acquisition value of the acquisition unit, and a target value calculation unit.
A target deviation calculation unit that calculates a target deviation based on a deviation between the target value and the acquisition value of the acquisition unit for each of the plurality of control parameters.
For the state parameter that is a parameter related to the state amount of the diesel engine and whose value changes depending on the control of the controlled object, an ideal value that is a value when the diesel engine is in a steady state under the acquired value of the acquisition unit is set. The ideal value calculation unit to calculate and
A state deviation calculation unit that calculates a state deviation that is a deviation between the ideal value and the acquisition value of the acquisition unit for the state parameter.
And a command value calculator for calculating a control command value for the controlled object by multiplying the deviation vector and the gain matrix including the components of the target deviation and the state deviation,
The control parameters include an intake EGR rate, which is the ratio of EGR gas to the working gas sucked by the diesel engine, and a boost pressure, which is the pressure in the passage through which the intake air supercharged by the variable displacement turbocharger flows. Engine system controller.
前記目標偏差演算部は、前記制御パラメーターの各々について偏差の積算値を前記目標偏差として演算する
請求項に記載のエンジンシステムの制御装置。
The control device for an engine system according to claim 1 , wherein the target deviation calculation unit calculates an integrated value of deviations for each of the control parameters as the target deviation.
前記複数のパラメーターの推定値をモデルを用いて演算する状態量推定部をさらに備え、前記取得部は、前記状態量推定部が演算した前記推定値を取得する
請求項1または2に記載のエンジンシステムの制御装置。
The engine according to claim 1 or 2 , further comprising a state quantity estimation unit that calculates estimated values of the plurality of parameters using a model, and the acquisition unit acquires the estimated values calculated by the state quantity estimation unit. System control device.
前記ディーゼルエンジンの状態量に関するパラメーターを観測する観測部を備え、
前記状態量推定部は、前記制御対象への制御指示値を入力ベクトルの構成要素、前記観測部の観測値を観測ベクトルの構成要素として用いるカルマンフィルター理論により前記推定値を演算するものであり、
前記観測部が観測するパラメーターは、前記ブースト圧、前記ディーゼルエンジンが吸入する空気の量である吸入空気量、インテークマニホールド内の温度であるインマニ温度、及び、前記可変容量型ターボチャージャーにおけるタービンの回転数であるタービン回転数を含む
請求項に記載のエンジンシステムの制御装置。
It is equipped with an observation unit that observes parameters related to the state quantity of the diesel engine.
The state quantity estimation unit calculates the estimated value by the Kalman filter theory using the control instruction value for the controlled object as a component of the input vector and the observed value of the observation unit as a component of the observation vector.
The parameters observed by the observation unit are the boost pressure, the intake air amount which is the amount of air taken in by the diesel engine, the intake manifold temperature which is the temperature in the intake manifold, and the rotation of the turbine in the variable displacement turbocharger. The control device for an engine system according to claim 3 , which includes a turbine rotation speed which is a number.
JP2016197186A 2016-10-05 2016-10-05 Engine system controller Active JP6753754B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016197186A JP6753754B2 (en) 2016-10-05 2016-10-05 Engine system controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016197186A JP6753754B2 (en) 2016-10-05 2016-10-05 Engine system controller

Publications (2)

Publication Number Publication Date
JP2018059444A JP2018059444A (en) 2018-04-12
JP6753754B2 true JP6753754B2 (en) 2020-09-09

Family

ID=61909772

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016197186A Active JP6753754B2 (en) 2016-10-05 2016-10-05 Engine system controller

Country Status (1)

Country Link
JP (1) JP6753754B2 (en)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007154809A (en) * 2005-12-07 2007-06-21 Toyota Motor Corp Control unit for internal combustion engine
JP2010229972A (en) * 2009-03-30 2010-10-14 Daihatsu Motor Co Ltd Control device
JP5854662B2 (en) * 2011-06-28 2016-02-09 三菱重工業株式会社 Control device and method for internal combustion engine

Also Published As

Publication number Publication date
JP2018059444A (en) 2018-04-12

Similar Documents

Publication Publication Date Title
CN105626275B (en) Feedforward turbocharger control method for engine with supercharger
US6868329B2 (en) Control method of EGR system of engine
EP2292913B1 (en) Controlling exhaust gas recirculation in a turbocharged engine system
EP1864012B1 (en) Coordinated multivariable control of fuel and air in engines
US10161337B2 (en) Control device for internal combustion engine
CN109072772B (en) Physics-based vehicle turbocharger control techniques
JP2018178870A (en) Control device of engine
JP2005301764A (en) Controller using controlled object model
EP3880950B1 (en) Spark ignited engine load extension with low pressure exhaust gas recirculation and delta pressure valve
JP6753754B2 (en) Engine system controller
JP6855328B2 (en) Internal combustion engine throttle valve controller
JP6930902B2 (en) Valve controller
JP2011043156A (en) Control device
JP2019019802A (en) Control device for engine
JP2015222066A (en) Internal combustion engine control unit
JP2013142376A (en) Control device for internal combustion engine
JP2004108216A (en) Fuel injector of engine
JP2019094854A (en) Control device of internal combustion engine
JP6453122B2 (en) Control device for variable capacity turbocharger
JP6513440B2 (en) Controller of variable displacement turbocharger
JP2020002905A (en) Control device of internal combustion engine
JP2019157756A (en) Valve control device
JP6264317B2 (en) Control device for internal combustion engine
JP2019210843A (en) EGR control device
KR20200013580A (en) Control system

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190905

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200421

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200512

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200708

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200728

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200820

R150 Certificate of patent or registration of utility model

Ref document number: 6753754

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250