JP2016099862A - プラント制御装置 - Google Patents
プラント制御装置 Download PDFInfo
- Publication number
- JP2016099862A JP2016099862A JP2014237547A JP2014237547A JP2016099862A JP 2016099862 A JP2016099862 A JP 2016099862A JP 2014237547 A JP2014237547 A JP 2014237547A JP 2014237547 A JP2014237547 A JP 2014237547A JP 2016099862 A JP2016099862 A JP 2016099862A
- Authority
- JP
- Japan
- Prior art keywords
- target value
- value
- model
- cand
- plant
- 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.)
- Pending
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/12—Improving ICE efficiencies
Landscapes
- Combined Controls Of Internal Combustion Engines (AREA)
- Feedback Control In General (AREA)
- Output Control And Ontrol Of Special Type Engine (AREA)
Abstract
【課題】リファレンスガバナによる修正目標値の決定手法において、予測値を計算するためのモデルを区分的アフィンARXモデル(PWARXモデル)で表現する場合においても、評価関数の勾配を安定的に算出するための手法を提供する。
【解決手段】サブモデルの切り替えに差異が生じた場合には、微小摂動Δを小さい値(例えば、Δ:=Δ/2)に変更して、制御出力予測値z^(wcand+Δ),z^(wcand−Δ)を再計算する。微小摂動Δを変更することで、サブモデルの切り替えタイミングが揃うようになるので、勾配∇の探索時に異常値が算出されることを回避して、修正目標値の演算の収束性や安定性を高めることができる。
【選択図】図9
【解決手段】サブモデルの切り替えに差異が生じた場合には、微小摂動Δを小さい値(例えば、Δ:=Δ/2)に変更して、制御出力予測値z^(wcand+Δ),z^(wcand−Δ)を再計算する。微小摂動Δを変更することで、サブモデルの切り替えタイミングが揃うようになるので、勾配∇の探索時に異常値が算出されることを回避して、修正目標値の演算の収束性や安定性を高めることができる。
【選択図】図9
Description
本発明は、プラント制御装置に関する。
従来、例えば特許文献1には、プラントの制御量の出力値を目標値に近づけるようにフィードバック制御によって当該プラントの制御入力を決定するフィードバックコントローラと、当該プラントと当該フィードバックコントローラとを含む閉ループシステムのモデルを用いて当該プラントの特定状態量の将来の予測値を計算し、当該予測値と当該特定状態量に課せられた制約とに基づいて当該フィードバックコントローラに与えられる目標値を修正するリファレンスガバナと、を備えるプラント制御装置が開示されている。
このプラント制御装置において、リファレンスガバナは、具体的に、最急降下法を用いた繰り返し計算によって、特定状態量zの将来の予測値z^および制約z−と、修正目標値候補wcandと、オリジナルの目標値rとを用いて表される評価関数J(wcand)の最適値を探索する。また、リファレンスガバナは、この評価関数J(wcand)を最適化する修正目標値候補wcandを、制御量の最終的な修正目標値wとして決定する。このようにしてオリジナルの目標値rが修正され、最終的な修正目標値wがフィードバックコントローラに入力される。
リファレンスガバナによる最終的な修正目標値wの決定手法は、次の通りである。先ず、オリジナルの目標値rに基づいて、修正目標値候補wcandを複数用意する。続いて、修正目標値候補wcandから所定の微小値Δだけ離れた2つの近傍値wcand+Δ,wcand−Δに対する予測値z^(wcand+Δ),z^(wcand−Δ)を、上述した予測モデルによって予測する。続いて、評価関数J(wcand)の勾配∇(={J(wcand+Δ)−J(wcand−Δ)}/2Δ)を計算する。最終的な修正目標値wは、勾配∇の計算を有限回繰り返したときに、勾配∇を最小化する修正目標値候補wcandとして決定される。
秋田敏和ら、「PWARXモデルによる前方車追従行動のモデル化と解析」、第50回自動制御連合講演会、2007年11月24日〜25日
ところで、例えば非特許文献1に開示されている区分的アフィンARXモデル(PWARXモデル)は、離散的な状態遷移に応じて複数のサブモデル(ARXモデル)を切り替えるシステム表現手法であり、非線形性が強くモデル化の難しいエンジンのようなプラントを扱うのに有用である。そこで、本発明者は、上述したリファレンスガバナによる修正目標値wの決定手法において、予測値z^を計算するための予測モデルをこのPWARXモデルで表現することを検討している。
しかし、PWARXモデルによって予測値z^(w+Δ),z^(w−Δ)を予測する場合、サブモデルの切り替えタイミングが揃わない可能性がある。例えば、ある区分領域iからその隣の区分領域jへの遷移(サブモデルの切り替え)が、予測値z^(w+Δ)の予測では3ステップ目で起こるのに対して、予測値z^(w−Δ)の予測では2ステップ目で起こるというケースが生じ得る。そうすると、2ステップ目や3ステップ目で求めた予測値z(w+Δ)と予測値z(w−Δ)の差が大きくなり、上述した勾配∇の計算式の分母が急変して異常値が算出されてしまう虞がある。
本発明は、上述のような課題に鑑みてなされたものである。即ち、予測値zを計算するためのモデルをPWARXモデルで表現する場合においても、評価関数J(w)の勾配∇を安定的に算出するための手法を提供することを目的とする。
本発明は、上記の目的を達成するため、プラント制御装置であって、
プラントの制御量の出力値をその目標値に近づけるようにフィードバック制御によって前記プラントの制御入力を決定するフィードバックコントローラと、
前記プラントと前記フィードバックコントローラとを含む閉ループシステムのモデルを用いて前記プラントの特定状態量の将来の予測値を計算すると共に、計算した予測値と、前記特定状態量に課せられた制約と、オリジナルの目標値と、修正目標値の候補とを用いて表される評価関数の最適値を最急降下法による繰り返し計算によって探索し、前記最適値を用いて前記フィードバックコントローラに与えられる前記目標値を修正するリファレンスガバナと、を備え、
前記モデルは、前記特定状態量の拡大状態ベクトルに基づいて特定される区分領域に応じて、複数のサブモデルを切り替え可能なPWARXモデルによって表現され、
前記リファレンスガバナは、前記オリジナルの目標値に基づいて前記修正目標値の候補を用意し、用意した修正目標値の候補から正方向および負方向に所定値だけ離れた2つの近傍値のそれぞれを前記PWARXモデルに適用して前記特定状態量の予測値を2つ求め、求めた2つの予測値に基づいて前記修正目標値の候補に対する前記評価関数の勾配を求め、求めた勾配を最小化する前記修正目標値の候補を前記最適値として決定するように構成され、
前記リファレンスガバナは、前記2つの近傍値を前記PWARXモデルに適用する際に、前記2つの近傍値間で前記サブモデルを切り替えるタイミングが揃わない場合は、前記所定値をより小さい値に変更するように構成されていることを特徴とする。
プラントの制御量の出力値をその目標値に近づけるようにフィードバック制御によって前記プラントの制御入力を決定するフィードバックコントローラと、
前記プラントと前記フィードバックコントローラとを含む閉ループシステムのモデルを用いて前記プラントの特定状態量の将来の予測値を計算すると共に、計算した予測値と、前記特定状態量に課せられた制約と、オリジナルの目標値と、修正目標値の候補とを用いて表される評価関数の最適値を最急降下法による繰り返し計算によって探索し、前記最適値を用いて前記フィードバックコントローラに与えられる前記目標値を修正するリファレンスガバナと、を備え、
前記モデルは、前記特定状態量の拡大状態ベクトルに基づいて特定される区分領域に応じて、複数のサブモデルを切り替え可能なPWARXモデルによって表現され、
前記リファレンスガバナは、前記オリジナルの目標値に基づいて前記修正目標値の候補を用意し、用意した修正目標値の候補から正方向および負方向に所定値だけ離れた2つの近傍値のそれぞれを前記PWARXモデルに適用して前記特定状態量の予測値を2つ求め、求めた2つの予測値に基づいて前記修正目標値の候補に対する前記評価関数の勾配を求め、求めた勾配を最小化する前記修正目標値の候補を前記最適値として決定するように構成され、
前記リファレンスガバナは、前記2つの近傍値を前記PWARXモデルに適用する際に、前記2つの近傍値間で前記サブモデルを切り替えるタイミングが揃わない場合は、前記所定値をより小さい値に変更するように構成されていることを特徴とする。
本発明によれば、PWARXモデルを用いて、修正目標値の候補から正方向および負方向に所定値だけ離れた2つの近傍値から特定状態量の将来の予測値を計算して評価関数の勾配を求める場合において、サブモデルを切り替えるタイミングが当該2つの近傍値間で揃わない場合は、当該所定値をより小さい値に変更できるので、勾配探索時に異常値が算出されることを回避して、修正目標値の演算の収束性や安定性を高めることができる。
本発明の実施の形態のプラント制御装置は、車両動力プラントであるディーゼルエンジンの制御装置である。図1は、本実施の形態に係る制御装置が適用されるシステムの構成を示す図である。図1に示すように、ディーゼルエンジンの本体2には4つの気筒が直列に備えられ、気筒ごとにインジェクタ8が設けられている。エンジン本体2には吸気マニホールド4と排気マニホールド6が取り付けられている。
吸気マニホールド4にはエアクリーナ20から取り込まれた新気が流れる吸気通路10が接続されている。吸気通路10にはターボ過給機14のコンプレッサが取り付けられている。このコンプレッサの下流にはインタークーラ22が備えられ、インタークーラ22の下流にはディーゼルスロットル24が設けられている。排気マニホールド6にはエンジン本体2からの排気を大気中に放出するための排気通路12が接続されている。排気通路12にはターボ過給機14のタービンが取り付けられている。ターボ過給機14は可変容量型であって、タービンには可変ノズル16が備えられている。
図1に示すシステムは、排気系から吸気系へ排気を再循環させるEGR装置を備えている。EGR装置は、吸気通路10におけるディーゼルスロットル24の下流と排気マニホールド6とをEGR通路30によって接続する高圧ループEGR装置である。EGR通路30にはEGR弁32が設けられている。但し、EGR装置は、吸気通路10におけるコンプレッサの上流と、排気通路12におけるタービンの下流とを、EGR通路30とは別のEGR通路によって接続する低圧ループEGR装置であってもよい。
図1に示すECU(Electronic Control Unit)40が本実施の形態に係る制御装置に相当する。ECU40は、RAM(ランダムアクセスメモリ)、ROM(リードオンリメモリ)、CPU(マイクロプロセッサ)等を備えている。ECU40は、各種センサの信号を取り込み処理する。各種センサには、エンジン回転速度を検出する回転速度センサ42や、アクセルペダルの開度に応じた信号を出力するアクセルペダル開度センサ44などが含まれている。ECU40は、取り込んだ各センサの信号を処理して所定の制御プログラムに従ってアクチュエータを操作する。ECU40によって操作されるアクチュエータには、可変ノズル16、ディーゼルスロットル24、EGR弁32などが含まれている。
本実施の形態において、ECU40は、ディーゼルエンジンの過給圧・EGR率制御を実行する。過給圧・EGR率制御における制御入力(操作量)は可変ノズル開度、EGR弁開度およびディーゼルスロットル開度であり、制御出力(特定状態量)は過給圧とEGR率である。ここで、過給圧とEGR率にはハード上或いは制御上の制約が課せられている。ECU40は、過給圧とEGR率がそれぞれの制約を満たし、尚且つ、それぞれの目標値に追従させるように制御入力を決定する。
図2は本実施の形態に係る制御装置が有する目標値追従制御構造を示す図である。なお、図2に示す目的値追従制御構造は、ECU40のROMに格納された制御プログラムに従いCPUが動作することで仮想的に実現される構成である。この目標値追従制御構造は、目標値マップ(MAP)、リファレンスガバナ(RG)およびフィードバックコントローラ(FBC)を備えている。
目標値マップは、ディーゼルエンジン(DE)の運転条件を示す外生入力d=[エンジン回転速度;燃料噴射量]が与えられると、ディーゼルエンジンの制御量の目標値r=[EGR率目標値;過給圧目標値]を出力する。
リファレンスガバナは、目標値rが与えられると、制御出力(特定状態量)z=[EGR率;過給圧]に関する制約が満たされるように目標値rを修正し、修正目標値w=[EGR率修正目標値;過給圧修正目標値]を出力する。リファレンスガバナの詳細については後述する。
フィードバックコントローラは、リファレンスガバナから修正目標値wが与えられると、ディーゼルエンジンの状態量x=[EGR率;過給圧]を修正目標値wに近づけるように、フィードバック制御によってディーゼルエンジンの制御入力u=[ディーゼルスロットル開度;EGR弁開度;可変ノズル開度]を決定する。フィードバックコントローラの仕様に限定はなく、公知のフィードバックコントローラを用いることができる。
図3は図2に示す目標値追従制御構造を等価変形して得られたフィードフォワード構造である。図2において破線で囲まれた閉ループシステムは既に設計済みであるとして、図3に示すフィードフォワード構造では1つのモデル(P)とされている。閉ループシステムのモデルは次のモデル式(1)で表される。式(1)において、f,hはモデル式の関数である。また、kは離散時間ステップを表している。
ところで、上述したPWARXモデルは、非線形システムの表現能力の高さと、データクラスタリング理論やシステム同定理論との親和性が高く、非常によく用いられるシステムモデル形式である。上記閉ループシステムに対するPWARXモデルは、次式(2)で記述される。式(2)において、sはサブモデルの切り替え総数であり、θ1,・・・θsは各サブモデルの係数であり、Bw,Bdは既知の係数である。また、x−(k)は拡大状態ベクトル(x−(k)=[z(k−1),z(k−2),…,z(k−Nz),1]T)である(Nzは整数、Tは転置記号)。
リファレンスガバナは、上記の式(1),(2)で表される予測モデルを用いてディーゼルエンジンの制御出力zの将来の予測値z^を計算する。上述したように、本実施の形態における制御出力zはEGR率および過給圧であり、制御出力zには制約が課せられている。制御出力zがその上限値z−以下であることが制御出力zに課せられた制約である。制御出力予測値z^の計算には、状態量xおよび外生入力dに加えて修正目標値wが用いられる。リファレンスガバナは、制御出力予測値z^と制御出力上限値z−とに基づき、次式(3)で表される評価関数J(w)を用いて修正目標値wを計算する。
式(3)に示す評価関数J(w)の右辺第1項は修正目標値候補wcandを変数とする目的関数である。この目標関数はオリジナルの目標値rと修正目標値候補wcandとの距離が小さいほど小さな値を取るように構成されている。評価関数J(w)の右辺第2項はペナルティ関数である。ペナルティ関数は制御出力予測値z^が制約z−に抵触する場合に目的関数にペナルティを加えるように構成されている。ペナルティ関数には、ペナルティに重みを付けるための重み定数ρが設定されている。なお、右辺第2項のz^(k+i|k)は時刻kの時点での情報に基づく時刻k+iの時点の制御出力予測値を表し、Nhは予測ホライズン(予測ステップ数)を表している。
また、本実施の形態においては、勾配法(具体的には最急降下法)を用いて式(3)の評価関数J(w)を最小化する修正目標値wを探索している。この勾配演算では、修正目標値候補wcandに対して微小摂動Δを加えたときのJ(w+Δ)とJ(w−Δ)の差を2Δで割ることで勾配∇を求めている。つまり、勾配∇={J(w+Δ)−J(w−Δ)}/2Δである。図4は、評価関数J(w)の勾配∇を説明するための図である。図4に示すように、修正目標値候補wcandの勾配∇(wcand)は、修正目標値候補wcandの近傍のwcand+Δ,wcand−Δの評価関数J(wcand+Δ),J(wcand−Δ)を通る直線の傾きとして表される。なお、図4においては、評価関数J(w)の構造が一次元の場合を示しているが、本実施の形態では目標値rがEGR率と過給圧の2つあるので、実際の評価関数J(w)の構造は二次元となる。
図5は、本実施の形態に係るリファレンスガバナアルゴリズムを示す図である。図5に示すように、本実施の形態では、修正目標値候補wcandに対して、現在のディーゼルエンジンの運転条件に基づく閉ループシステムの将来予測と、この予測結果を用いた評価関数J(w)の演算とが有限回反復される。これにより、式(3)に示す評価関数J(w)を最小にする修正目標値候補wcandを選択し、選択した修正目標値候補wcandを最終的な修正目標値wとして決定する。
ところで、式(3)に示す評価関数の右辺第2項の制御出力予測値z^(k+i|k)は、修正目標値候補wcandの近傍値wcand+Δ,wcand−Δのそれぞれに対して計算される。ここで、式(2)に示したように、PWARXシステムは、x-(k)が属する区分領域に応じて使用するサブモデルの切り替えが行われる。そのため、微小摂動Δが大きい場合には、制御出力予測値z^(wcand+Δ),z^(wcand−Δ)の計算中にサブモデルの切り替えタイミングが揃わなくなる可能性がある。
図6乃至図8は、制御出力予測値z^(wcand+Δ),z^(wcand−Δ)の計算中に発生する問題を説明するための図である。なお、図6では、4ステップ先までの将来予測をしているものとする。図6に示すように、wcand+Δに対する将来予測では、サブモデルの切り替えが2ステップ目で発生しているのに対し、wcand−Δに対する将来予測では、サブモデルの切り替えが3ステップ目で発生している。このように、サブモデルの切り替えタイミングが揃わなくなると、図7に示すように、制約抵触量がwcand+Δとwcand−Δで大きく異なることになる。制約抵触量は式(3)の右辺第2項中のz^(k+i|k)−z−に相当するものである。そのため、制約抵触量がwcand+Δとwcand−Δで大きく異なれば、評価関数J(wcand+Δ)と評価関数J(wcand−Δ)も大きく異なることになるので、修正目標値候補wcandの勾配∇(wcand)の勾配∇の計算式の分母が急変する。そうすると、図8に示すような異常値が算出されてしまう。
そこで、本実施の形態では、サブモデルの切り替えに差異が生じた場合には、微小摂動Δを小さい値(例えば、Δ:=Δ/2)に変更して、制御出力予測値z^(wcand+Δ),z^(wcand−Δ)を再計算する。図9は、微小摂動をΔからΔ/2に変更した場合の、変更前後の制御出力予測値z^(wcand+Δ),z^(wcand−Δ)の推移を示した図である。なお、図9では、図6同様、4ステップ先までの将来予測をしているものとする。図9に示すように、微小摂動Δを変更することで、サブモデルの切り替えタイミングが揃うようになるので、勾配∇の探索時に異常値が算出されることを回避して、修正目標値の演算の収束性や安定性を高めることができる。
なお、本発明は上述の実施の形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。例えば、図2に示した目標値追従制御構造は、ディーゼルエンジンが低圧ループEGR装置と高圧ループEGR装置とを備える場合には、図10の(a)〜(d)に示すような制御入力と制御出力との組み合わせにも適用することができる。図10の(a)および(b)では、可変ノズル開度(VN開度)やディーゼルスロットル開度(Dスロ開度)の他に、低圧ループEGR装置のEGR弁開度(LPL−EGR弁開度)と高圧ループEGR装置のEGR弁開度(HPL−EGR弁開度)とが制御入力に含まれている。図10の(c)および(d)では、EGR率の代わりに、低圧ループEGR装置のEGR量(LP−EGR量)と高圧ループEGR装置のEGR量(HP−EGR量)とが制御出力に含まれている。
2 エンジン本体
16 可変ノズル
24 ディーゼルスロットル
32 EGR弁
40 ECU
16 可変ノズル
24 ディーゼルスロットル
32 EGR弁
40 ECU
Claims (1)
- プラントの制御量の出力値をその目標値に近づけるようにフィードバック制御によって前記プラントの制御入力を決定するフィードバックコントローラと、
前記プラントと前記フィードバックコントローラとを含む閉ループシステムのモデルを用いて前記プラントの特定状態量の将来の予測値を計算すると共に、計算した予測値と、前記特定状態量に課せられた制約と、オリジナルの目標値と、修正目標値の候補とを用いて表される評価関数の最適値を最急降下法による繰り返し計算によって探索し、前記最適値を用いて前記フィードバックコントローラに与えられる前記目標値を修正するリファレンスガバナと、を備え、
前記モデルは、前記特定状態量の拡大状態ベクトルに基づいて特定される区分領域に応じて、複数のサブモデルを切り替え可能なPWARXモデルによって表現され、
前記リファレンスガバナは、前記オリジナルの目標値に基づいて前記修正目標値の候補を用意し、用意した修正目標値の候補から正方向および負方向に所定値だけ離れた2つの近傍値のそれぞれを前記PWARXモデルに適用して前記特定状態量の予測値を2つ求め、求めた2つの予測値に基づいて前記修正目標値の候補に対する前記評価関数の勾配を求め、求めた勾配を最小化する前記修正目標値の候補を前記最適値として決定するように構成され、
前記リファレンスガバナは、前記2つの近傍値を前記PWARXモデルに適用する際に、前記2つの近傍値間で前記サブモデルを切り替えるタイミングが揃わない場合は、前記所定値をより小さい値に変更するように構成されていることを特徴とするプラント制御装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014237547A JP2016099862A (ja) | 2014-11-25 | 2014-11-25 | プラント制御装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014237547A JP2016099862A (ja) | 2014-11-25 | 2014-11-25 | プラント制御装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2016099862A true JP2016099862A (ja) | 2016-05-30 |
Family
ID=56077187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014237547A Pending JP2016099862A (ja) | 2014-11-25 | 2014-11-25 | プラント制御装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2016099862A (ja) |
-
2014
- 2014-11-25 JP JP2014237547A patent/JP2016099862A/ja active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5737537B2 (ja) | スライディングモードコントローラ | |
EP3511555B1 (en) | Plant control system | |
WO2014033858A1 (ja) | プラント制御装置 | |
JP2010216305A (ja) | 内燃機関の過給圧制御装置 | |
JP2014127083A (ja) | 車両動力プラントの制御装置 | |
JP2018067046A (ja) | プラント制御装置 | |
Nakada et al. | An application study of online reference governor to boost pressure control for automotive diesel engines | |
JP2016115252A (ja) | プラント制御装置 | |
EP3147488B1 (en) | Internal-combustion-engine control device | |
JP2016169688A (ja) | 内燃機関の制御装置 | |
JP2016099862A (ja) | プラント制御装置 | |
JP2016130877A (ja) | プラント制御装置 | |
JP2016098771A (ja) | 内燃機関の制御装置 | |
JP2017106333A (ja) | プラント制御装置 | |
JP2018063586A (ja) | プラント制御装置 | |
JP2015197087A (ja) | 内燃機関の制御装置 | |
JP2016091041A (ja) | プラント制御装置 | |
JP2017020357A (ja) | 内燃機関の制御装置 | |
JP2011043154A (ja) | 制御装置 | |
JP2014074987A (ja) | プラント制御装置 | |
EP2749760B1 (en) | Vehicle power plant control apparatus | |
JP2016157310A (ja) | プラント制御装置 | |
JP2012067662A (ja) | エンジン制御プログラム及び装置 | |
JP2011001871A (ja) | 制御装置 | |
JP2018080628A (ja) | プラント制御装置 |