JP7529545B2 - 制御装置 - Google Patents

制御装置 Download PDF

Info

Publication number
JP7529545B2
JP7529545B2 JP2020193678A JP2020193678A JP7529545B2 JP 7529545 B2 JP7529545 B2 JP 7529545B2 JP 2020193678 A JP2020193678 A JP 2020193678A JP 2020193678 A JP2020193678 A JP 2020193678A JP 7529545 B2 JP7529545 B2 JP 7529545B2
Authority
JP
Japan
Prior art keywords
input
control
output
dynamic
static
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
JP2020193678A
Other languages
English (en)
Other versions
JP2022082233A (ja
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.)
Kobelco Construction Machinery Co Ltd
Hiroshima University NUC
Original Assignee
Kobelco Construction Machinery Co Ltd
Hiroshima University NUC
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 Kobelco Construction Machinery Co Ltd, Hiroshima University NUC filed Critical Kobelco Construction Machinery Co Ltd
Priority to JP2020193678A priority Critical patent/JP7529545B2/ja
Priority to CN202180076777.4A priority patent/CN116457731A/zh
Priority to PCT/JP2021/035961 priority patent/WO2022107456A1/ja
Priority to EP21894331.4A priority patent/EP4235316A4/en
Priority to US18/252,983 priority patent/US20230418241A1/en
Publication of JP2022082233A publication Critical patent/JP2022082233A/ja
Application granted granted Critical
Publication of JP7529545B2 publication Critical patent/JP7529545B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/021Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a variable is automatically adjusted to optimise the performance
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Program-controlled manipulators
    • B25J9/16Program controls
    • B25J9/1628Program controls characterised by the control loop
    • B25J9/163Program controls characterised by the control loop learning, adaptive, model based, rule based expert control
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16HGEARING
    • F16H61/00Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing
    • F16H61/38Control of exclusively fluid gearing
    • F16H61/40Control of exclusively fluid gearing hydrostatic
    • F16H61/46Automatic regulation in accordance with output requirements
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60YINDEXING SCHEME RELATING TO ASPECTS CROSS-CUTTING VEHICLE TECHNOLOGY
    • B60Y2200/00Type of vehicle
    • B60Y2200/40Special vehicles
    • B60Y2200/41Construction vehicles, e.g. graders, excavators
    • B60Y2200/412Excavators
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16HGEARING
    • F16H61/00Control functions within control units of change-speed- or reversing-gearings for conveying rotary motion ; Control of exclusively fluid gearing, friction gearing, gearings with endless flexible members or other particular types of gearing
    • F16H61/38Control of exclusively fluid gearing
    • F16H61/40Control of exclusively fluid gearing hydrostatic
    • F16H61/4148Open loop circuits
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/42Servomotor, servo controller kind till VSS
    • G05B2219/42018Pid learning controller, gains adapted as function of previous error
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/42Servomotor, servo controller kind till VSS
    • G05B2219/42036Adaptive control, adaptive nonlinear control
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/42Servomotor, servo controller kind till VSS
    • G05B2219/42038Real time adaptive control
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/42Servomotor, servo controller kind till VSS
    • G05B2219/42043Adapt regulator as function of its output

Landscapes

  • Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Feedback Control In General (AREA)

Description

本発明は、制御対象を制御する制御装置に関するものである。
油圧ショベルをはじめとする作業機械を制御する制御装置では自動運転等の複雑な制御が要求される。このような複雑な制御が要求される制御装置では制御構造が階層化され、上流のコントローラからの指示に従って、制御対象を直接的に制御する下流の制御ループが作動される。しかしながら、制御対象の入出力特性が大きく変動して下流の制御ループの入出力特性が大きく変動すると、上流のコントローラが適切に作用せず、制御対象を適切に制御できなくなる。そのため、下流の制御ループを理想的な入出力特性に維持できれば、上流のコントローラを下流の制御ループの理想的な入出力特性に基づいて設計できる。
そこで、下流の制御ループの入出力特性の変動を抑制する技術が提案されている。例えば、非特許文献1には、上流のコントローラから入力される制御入力(u´)及び制御対象から出力される制御出力(y)の理想的な入出力特性を示す入出力モデル(Pm)と、制御入力(u´)に対する入出力モデル(Pm)の理想出力(ym)及び制御出力(y)のモデル誤差を軽減するための補償入力(uc)を生成する補償器(D)と、制御入力(u´)から補償入力(uc)を減じて制御対象への実入力(u)を算出するする減算器とを備える制御装置が開示されている。
梅井、岡島、松永、浅井、「モデル誤差抑制補償器の多入出力システムに対する設計」、システム制御情報学会論文誌、Vol27 No.2、pp.67-72、2014
しかしながら、非特許文献1の技術では、制御対象の入出力特性が大きく変動した場合、補償器ではその変動を十分に補償できなくなり、上流のコントローラは制御対象を適切に制御することができなくなる。
本発明は、上記の課題を解決するためになされたものであり、制御対象の入出力特性が大きく変動しても、設計時のコントローラを用いて制御対象を適切に制御することができる制御装置を提供することを目的とする。
本発明の一態様に係る制御装置は、制御対象を制御する制御装置であって、前記制御対象から出力される制御出力と目標出力との偏差を零にするための制御入力を算出するコントローラと、静的パラメータと前記制御入力とに基づいて、前記制御対象の静的特性の変動を補償する静的補償入力を算出する静的補償器と、動的パラメータと前記制御出力とに基づいて前記制御対象の動的特性の変動を補償する動的補償入力を算出する動的補償器と、前記静的補償入力と前記動的補償入力とを合成して実入力を算出し、前記制御対象に入力する合成器と、前記制御入力と前記制御出力との理想的な入出力関係を規定する入出力モデルを用いて、前記制御入力に対応する理想出力を算出する理想出力算出器と、前記制御出力と前記理想出力との差が最小化されるように前記静的パラメータ及び前記動的パラメータをそれぞれ調整するパラメータ調整器とを備える。
本構成によれば、静的パラメータ及び動的パラメータが調整されており、固定されていないため、制御対象の入出力特性が大きく変動したとしても、制御装置の設計時のコントローラを用いて制御対象を適切に制御できる。すなわち、制御入力と制御出力との理想的な入出力特性を示す入出力モデルを用いて制御入力に対応する理想出力が算出され、理想出力と制御出力との差が最小化されるように、静的補償器が有する静的パラメータと動的補償器が有する動的パラメータとが調整されている。これにより、制御対象の入出力特性の特性が大きく変動しても、制御入力と制御出力との入出力特性は理想的な入出力特性に維持される。そのため、制御対象の入出力特性が大きく変動したとしても、設計時のコントローラを用いて制御対象を適切に制御できる。これにより、コントローラの設計が容易化され、制御装置の開発が円滑化される。
さらに、本構成によれば、制御対象の稼働中に算出される理想出力と制御出力とを用いて静的パラメータ及び動的パラメータが調整されている。そのため、制御対象を含む機器の動作を止めることなく稼働中に静的パラメータ及び動的パラメータを調整するオンライン調整が可能となる。なお、特許文献1では、補償器のパラメータの調整は予定されていないため、パラメータを調整するためには機器の動作を止める必要がある。
さらに、動的パラメータ及び制御出力に基づいて算出された動的補償入力によって制御入力が補正されているため、立ち上がり特性及び減衰特性等の制御対象の動的特性の変動を補償することができる。さらに、制御入力と静的パラメータとに基づいて算出された静的補償入力によって制御入力が補正されているため、動的補償入力の合成に伴う制御入力のスケールの変動等の制御特性の静的特性の変動を抑制できる。
上記制御装置において、前記動的パラメータは、比例ゲインと微分ゲインとを含み、前記動的補償器は、前記制御出力に前記比例ゲインを乗じた値と前記制御出力の微分値に前記微分ゲインを乗じた値との和を前記動的補償入力として算出し、前記パラメータ調整器は、前記制御出力と前記理想出力との差が最小化されるように前記静的パラメータ、前記前記比例ゲイン、及び前記微分ゲインをそれぞれ調整してもよい。
本構成によれば、制御出力に比例ゲインを乗じた値と制御出力の微分値に微分ゲインを乗じた値との和が動的補償入力として算出されているため、制御対象の動的特性の変動を正確に補償することができる。そして、このようにして調整された、静的パラメータ、比例ゲイン、及び微分ゲインを用いて制御入力が補正されている。そのため、制御対象の入出力特性が大きく変動したとしても、設計時のコントローラを用いて制御対象を適切に制御できる。
上記制御装置において、前記パラメータ調整器は、逐次最小二乗法を用いて前記静的パラメータ及び前記動的パラメータを調整してもよい。
本構成によれば、逐次最小二乗法を用いて静的パラメータ及び動的パラメータが調整されているため、静的パラメータ及び動的パラメータがリアルタイムに調整される。その結果、制御対象の入出力特性が変動した場合、速やかに静的パラメータ及び動的パラメータを調整することができ、制御対象を適切に制御できる。
上記制御装置において、前記合成器は、前記静的補償入力から前記動的補償入力を減じた値を前記実入力として算出してもよい。
本構成によれば、静的補償入力から動的補償入力を減じることで実入力が算出されているため、実入力を容易に算出できる。
上記制御装置において、前記制御対象は、油圧システムであり、油圧システムは、油圧ポンプと、前記油圧ポンプを駆動する電動機と、前記油圧ポンプから供給される作動油によって作動するアクチュエータとを備え、前記油圧システムは、前記アクチュエータの状態値を前記制御出力として出力し、前記合成器は、前記電動機の指令値を前記実入力として算出してもよい。
油圧システムは、例えば作動油の温度が変化して入出力特性が大きく変動することがある。例えば、油圧システムとして油圧ショベルを用いた場合、バケットが掬い上げた土砂により質量及び慣性が変化し、それに伴って油圧システムの入出力特性が大きく変動することがある。本構成では、このような変動に対して静的パラメータ及び動的パラメータが調整されるため、油圧システムを適切に制御することができる。
本発明によれば、制御対象の入出力特性が大きく変動しても、設計時のコントローラを用いて制御対象を適切に制御することができる。
本発明の実施の形態に係る制御装置の構成の一例を示すブロック図である。 制御装置の処理の一例を示すフローチャートである。 制御ループを構成するフィードバック系を示すブロック図である。 制御対象として用いられた油圧モータ制御システムの油圧回路図である。 油圧モータ制御システムに対する制御装置の適用結果を示すグラフである。 静的ゲイン、比例ゲイン、及び微分ゲインの調整結果を示すグラフである。 図6を拡大して示したグラフである。
以下、本発明の実施の形態について図面を参照して説明する。図1は、本発明の実施の形態に係る制御装置1の構成の一例を示すブロック図である。制御装置1は、実入力up(k)に応答して制御出力y(k)を出力する制御対象100を制御する。制御対象100は、例えば、油圧ショベル及び油圧クレーンを含む作業機械であってもよいし、自動車であってもよい。詳細には、制御対象100は、作業機械の油圧システムであってもよいし、自動車の自動ブレーキシステムであってもよい。参照符号の括弧内に記載されたkは時刻を示す。
制御装置1は、目標設定器2、減算器3、コントローラ4、静的補償器5、動的補償器6、減算器7(合成器の一例)、パラメータ調整器9、減算器8、理想出力算出器10、メモリ11、及び検出器12を含む。目標設定器2、減算器3、コントローラ4、静的補償器5、動的補償器6、減算器7、減算器8、パラメータ調整器9、及び理想出力算出器10は、例えばCPU又はASIC等のプロセッサで構成される。目標設定器2、メモリ11、及び検出器12は制御対象100の構成要素から省かれてもよい。
目標設定器2は、制御出力y(k)の目標である目標出力r(k)を設定する。目標出力r(k)は制御対象100に応じて種々の値が採用される。例えば、制御対象100が作業機械であれば、上部旋回体を旋回する電動機に対する目標回転数が採用される。例えば、制御出力y(k)が予め定められた目標パターンに追従するように制御対象100が制御される場合、目標出力r(k)は目標パターンを構成する目標値の時系列データが採用される。例えば、目標出力r(k)は、作業装置又は上部旋回体を動作させる操作レバーに入力された操作量に応じた目標値であってもよい。例えば、目標出力r(k)は予め定められた一定の値が採用されてもよい。
減算器3は、目標出力r(k)から制御出力y(k)を減じることで、偏差e(k)を算出する。
コントローラ4は、制御出力y(k)に基づいて、偏差e(k)を零にするための制御入力uc(k)を算出する。コントローラ4は背景技術で説明した上流のコントローラに相当する。コントローラ4は、例えばPID制御により偏差e(k)を零にするための制御入力uc(k)を算出すればよい。PID制御に用いられる数式としては、例えば後述の式(17)が採用される。なお、コントローラ4は、PID制御以外のP制御、PD制御、及びPI制御等の種々のフィードバック制御又はフィードフォワード制御を用いて制御入力uc(k)を算出してもよい。
静的補償器5は、静的ゲインf0(静的パラメータの一例)を制御入力uc(k)に乗じ、制御対象100の静的特性の変動を補償する静的補償入力を算出する。静的特性とは、時間に依存しない制御対象100の特性である。静的特性としては、例えば、制御出力y(k)がとり得るスケールが該当する。静的ゲインf0はこの静的特性の変動を補償するためのゲインである。例えば、動的補償器6により算出される動的補償入力が過大になると実入力up(k)が過小になって、制御出力y(k)の値が想定されるスケールよりも大幅に小さくなる。このような事態を回避するために、静的補償器5は、静的ゲインf0を制御入力uc(k)に乗じる。
動的補償器6は、動的ゲイン(動的パラメータの一例)と制御出力y(k)とに基づいて制御対象100の動的特性の変動を補償する動的補償入力を算出する。動的特性とは、例えば制御対象100の立ち上がり特性及び減衰特性等の時間に依存する特性である。動的ゲインは、このような動的特性の変動を補償するためのゲインである。動的ゲインには、例えば比例ゲインKp及び微分ゲインKが含まれる。動的補償器6は、例えばKp・y(k)+K・Δy(k)の演算式により、動的補償入力を算出する。ここで、Δy(k)はy(k)の微分を表す。
減算器7は、静的補償入力から動的補償入力を減じることにより、実入力up(k)を算出し、実入力up(k)を制御対象100に入力する。これにより、制御対象100の動的特性及び静的特性が補償されるように制御入力uc(k)が調整される。実入力up(k)は例えば、下記の式で表される。
up(k)=f0・uc(k)-Kp・y(k)-K・Δy(k)
上述の静的補償器5、動的補償器6、減算器7、及び制御対象100は制御ループ50を構成する。制御ループ50は、制御対象100を直接的に制御する下流の制御ループである。制御ループ50は、制御入力uc(k)に応答して制御出力y(k)を出力する。
理想出力算出器10は、制御入力uc(k)と制御出力y(k)との理想的な入出力関係を示す伝達関数である入出力モデルGm(z-1)を用いて、制御入力uc(k)に対応する理想出力yr(k)を算出する。理想的な入出力関係とは、コントローラ4の設計時における、制御入力uc(k)と制御出力y(k)との関係が該当する。以下、制御入力uc(k)と制御出力y(k)との関係を制御ループ50の入出力特性と呼ぶ。例えば、コントローラ4が、初期の制御対象100を含む初期の制御ループ50の入出力特性に基づいて設計されている場合、入出力モデルは初期の制御ループ50の入出力特性を有している。したがって、制御対象100の入出力特性が初期の特性から変化して、制御ループ50の入出力特性が初期の入出力特性から変化したとしても、理想出力算出器10は、初期の制御ループ50の入出力特性にしたがった理想出力yr(k)を算出できる。入出力モデルGm(z-1)は例えば後述の式(19)、(20)、(21)で表される。
減算器8は、制御出力y(k)から理想出力yr(k)を減じて差分Aを算出し、差分Aをパラメータ調整器9に入力する。
パラメータ調整器9は、減算器8から入力された差分Aが最小化されるように静的ゲインf0及び動的ゲイン(Kp、K)をそれぞれ調整する。パラメータ調整器9は例えば逐次最小二乗法により静的ゲインf0及び動的ゲイン(Kp、K)を算出してもよい。この場合、制御装置1のサンプリング時間に同期して、静的ゲインf0及び動的ゲイン(Kp、K)が調整される。すなわち、静的ゲインf0及び動的ゲイン(Kp、K)のオンライン調整が可能となる。逐次最小二乗法としては、後述の式(9)で示す評価関数Jを後述の式(10)~(16)を用いて最小化する手法が採用できる。
メモリ11は、例えばRAM又はフラッシュメモリで構成される。メモリ11は、制御出力y(k)及び理想出力yr(k)を記憶する。なお、メモリ11は、時刻kから数サンプル前までに算出された制御出力y(k)及び理想出力yr(k)を記憶してもよい。
検出器12は、例えばセンサで構成され、制御対象100から出力される制御出力y(k)を検出する。例えば、制御対象100がアクチュエータで構成される場合、検出器12は、アクチュエータの状態を検出するセンサで構成される。例えば、制御対象100が上部旋回体を旋回させる電動機で構成される場合、検出器12は電動機の回転数を検出するセンサで構成される。
次に、制御装置1の処理について説明する。図2は、制御装置1の処理の一例を示すフローチャートである。ステップS1において、減算器3は目標出力r(k)から制御出力y(k)を減じて偏差e(k)を算出する。
ステップS2において、コントローラ4は、偏差e(k)と制御出力y(k)とを式(17)に入力して制御入力uc(k)を算出する。
ステップS3において、理想出力算出器10は、制御入力uc(k)と式(19)で示す入出力モデルGm(z-1)とを乗じ、理想出力yr(k)を算出する。
ステップS4において、検出器12は、制御入力uc(k)に対する応答として制御ループ50から出力された制御出力y(k)を検出する。
ステップS5において、減算器8は、検出器12で検出された制御出力y(k)から理想出力yr(k)を減じて差分Aを算出する。
ステップS6において、パラメータ調整器9は差分Aが最小化されるように逐次最小二乗法を用いて静的ゲインf0及び動的ゲイン(Kp、K)を算出する。ステップS6が終了すると、処理はステップS1に戻る。これにより、逐次的に静的ゲインf0及び動的ゲイン(Kp、K)が調整される。
このように、制御装置1によれば、制御入力uc(k)と制御出力y(k)との理想的な入出力特性を示す入出力モデルGm(z-1)を用いて制御入力uc(k)に対応する理想出力yr(k)が算出され、理想出力yr(k)と制御出力y(k)との差分Aが最小化されるように、静的補償器5が有する静的ゲインf0と動的補償器6が有する動的ゲイン(Kp、K)とが調整される。これにより、制御対象100の入出力特性が大きく変動しても、制御入力uc(k)と制御出力y(k)との入出力特性は、コントローラ4の設計時における理想的な入出力特性に維持される。そのため、制御対象100の入出力特性が大きく変動しても、設計時のコントローラ4を用いて制御対象100を適切に制御できる。これにより、コントローラ4の設計が容易化され、制御装置1の開発が円滑化される。
なお、本発明は以下の変形例が採用できる。
(1)パラメータ調整器9は、データベース駆動型制御手法を用いて静的ゲインf0及び動的ゲイン(Kp、K)を調整してもよい。データベース駆動型制御手法は、データベースに記憶された過去に算出されたパラメータに基づいて、制御対象の現在の状態に適合するパラメータを算出する手法である。
この手法を採用する場合、制御装置1は、過去に算出された静的ゲインf0及び動的ゲイン(Kp、K)を記憶するデータベースをさらに備える。パラメータ調整器9は、制御対象100の現在の状態を示す要求点をメモリ11から取得する。要求点は例えば1サンプルから数サンプル前までの制御出力y(k)及び理想出力yr(k)を含む。パラメータ調整器9は、要求点とデータベースに記憶されたパラメータセットとのそれぞれの距離を算出し、距離が短い順にk個のパラメータセットを抽出する。パラメータセットは例えば1組の静的ゲインf0、比例ゲインKp、及び微分ゲインKを含む。パラメータ調整器9は、抽出したk個のパラメータセットのそれぞれについて距離が短いほど値が大きくなるように重み係数を求める。パラメータ調整器9は、算出した重み係数を用いてk個のパラメータセットを平均し、最終的なパラメータセットを算出し、最終的なパラメータセットを静的ゲインf0及び動的ゲイン(Kp、K)として算出する。
(2)動的補償器6が動的補償入力を算出する際に用いる演算式には、制御出力y(k)の2次の微分項と2次の微分ゲインとの積が含まれていてもよい。さらに、この演算式には、制御出力y(k)のi次の微分項とi次の微分ゲインとの積を、i=1からi=n(nは正の整数)まで加算した値が含まれていてもよい。
(実施例)
次に、本発明の実施例について説明する。まず、制御ループ50の設計について説明する。図3は、制御ループ50を構成するフィードバック系を示すブロック図である。このフィードバック系は、次式で表される。
ここで、up(k)、y(k)、uc(k)、Pはそれぞれ実入力、制御出力、制御入力、制御対象を表す。また、Δは差分演算子を表しており、後退演算子z-1を用いてΔ=1-z-1と表される。f0(k)、Kp(k)、K(k)はそれぞれパラメータを表している。パラメータ調整器9は、f0(k)、Kp(k)、K(k)のパラメータを逐次最小二乗法によりオンラインでチューニングする。逐次最小二乗法のメリットは計算コストが低いことにある。パラメータ調整器9は、操業データ(実入力up(k)、制御出力y(k))から静的補償器5及び動的補償器6のパラメータを算出する。
次に、操業データに基づくパラメータの調整法について説明する。
f0(k)=0でないと仮定すると、式(1)は以下のように変形される。
但し、式(3)において、θ1(k)、θ2(k)、θ3(k)は式(4)で表される。
また、理想的な制御ループ50の伝達関数を示す入出力モデルGm(z-1)に制御入力uc(k)を入力した際に得られる応答を理想出力yr(k,θ(k))とする。この場合、理想出力yr(k,θ(k))は式(5)で表される。
式(3)と式(5)の関係より次式が得られる。
評価関数Jは以下のように定義される。
但し、Nはデータの総数であり、評価関数Jを最小化することにより、制御出力y(k)が理想出力yr(k)に追従するようにパラメータθ(k)が調整される。よって、最適化されたパラメータを用いることにより、静的補償器5と動的補償器6と制御対象100とを含む制御ループ50の入出力特性と入出力モデルGm(z-1)の入出力特性とを一致させることが可能となる。
次に、式(9)の二乗和を最小化するため以下に示す逐次最小二乗法が適用される。
ωは忘却係数である。θ(k)及びψ(k)は次式で表される。
誤差共分散行列Γ(k)の初期値Γ(0)と推定値θ(k)の初期値θ(0)とは次式で定められる。
αはα>0を満たす任意の実数である。Iは3×3の単位行列である。θi(0)は任意の実数である。f0は0でないとの条件よりθi(0)は0でないと定める。
次に、本実施例に係る制御装置1の構成について説明する。制御装置1は図1で表される。
制御ループ50はフィードバック比例積分(F-PD)制御系で構成された下流の制御ループである。コントローラ4は、上流の制御ループである。コントローラ4は、固定の制御パラメータを有するPID(比例積分微分)制御系で構成される。
図1の構成では、制御ループ50の入出力特性と入出力モデルGm(z-1)の入出力特性とが一致するように静的補償器5及び動的補償器6のパラメータが調整されている。これにより、下流の制御ループ50が入出力モデルGm(z-1)と同等の入出力特性を有することになる。その結果、上流のコントローラ4を理想的な入出力モデルGm(z-1)に基づいて設計することが可能となる。
本実施例ではコントローラ4は、式(17)で示すPID制御系で構成される。
kcは比例ゲインを示し、TIは積分時間[s]を示し、TDは微分時間[s]を示す。
次に、実施例に係る制御装置1を油圧モータ制御システムに適用したシミュレーションについて説明する。
図4は、制御対象100として用いられた油圧モータ制御システムの油圧回路図である。油圧モータ制御システムは、電動機401(サーボモータ)、油圧ポンプ402、配管403、油圧モータ404、慣性体405、油圧センサ406、リリーフバルブ407、オイルタンク408、及び制御装置1を含む。
制御装置1は、アナログ電圧で表される指示回転数Npを電動機401に入力する。この指示回転数Npに従って電動機401は駆動する。油圧ポンプ402は、電動機401に接続されており、電動機401の動力によって回転する。油圧ポンプ402の回転により、オイルタンク408から配管403に作動油が流入する。油圧ポンプ402により吐出流量Qpで作動油は押し出され、配管403内に圧力Pmが発生する。この圧力Pmにより油圧モータ404は油圧動力Tで回転する。
油圧モータ404の回転数はNmで表され、慣性体405の慣性モーメント及び粘性減衰係数はそれぞれI、Dで表される。油圧モータ404の動力として使用された作動油は再びオイルタンク408に戻り、油圧ポンプ402によって再利用される。
制御装置1は、検出器12で回転数Nmを検出し、検出した回転数Nmが目標回転数になるように、指示回転数Npを算出する。指示回転数Npは実入力up(k)に相当し、回転数Nmは制御出力y(k)に相当する。
次に、油圧モータ制御システムをモデル化したシミュレーションモデルによる本実施例の検証について説明する。
本検証では、制御ループ50の理想的な入出力モデルGm(z-1)を以下のように設計した。
分母のP(z-1)は次式で表される。係数p1、p2は次式で表される。
Tsはサンプリング時間、σ,δはそれぞれ制御対象100の立ち上がり特性及び減衰特性等の動的パラメータを表している。これらの動的パラメータは、制御対象100の入出力特性に基づいて設計者が任意に設定する。
次に、シミュレーションの結果について説明する。図5は、油圧モータ制御システムに対する制御装置1の適用結果を示すグラフである。図6は、静的ゲインf0、比例ゲインKp、微分ゲインKの調整結果を示すグラフである。図7は、図6を拡大して示したグラフである。
図5に示すように、本シミュレーションでは、開始から40[s]が経過した時点において慣性体405の慣性モーメントが1/10に変更されている。ここでは、実機の運用上の都合により、理想出力yrが100から0に変動した場合、制御出力yはランプ状に減少されている。図5において、理想出力yr、制御出力y、制御入力uc、実入力up、圧力Pm、及び吐出流量Qpは正規化されている。
図5に示す過渡応答において、油圧モータ404の回転数Nmである制御出力yが慣性モーメントの変動後も理想出力yrに追従していることが分かる。これは、図6及び図7に示すように、慣性モーメントの変動に応じて、静的ゲインf0、比例ゲインKp、微分ゲインKがチューニングされているからである。
一方、図5に示すように、慣性モーメントの変動前後において、コントローラ4からの制御入力ucは大きく変動していないことが分かる。したがって、上流のコントローラ4は下流の制御ループ50の入出力特性の変動の影響を受けることなく動作していることが分かる。
このように、本実施例では、静的ゲインf0、比例ゲインKp、及び微分ゲインKがオンラインで調整されているため、制御対象100の特性が変動しても制御ループ50の入出力特性は、入出力モデルGm(z-1)の理想的な入出力特性と一致する。その結果、入出力モデルGm(z-1)に基づいて、コントローラ4を設計することができる。よって、コントローラ4の設計が容易になり、制御装置1の開発が円滑化される。さらに、下流の制御ループ50の入出力特性の変動の影響が制御装置1の全体に与える影響も抑制することが可能となる。
1 :制御装置
2 :目標設定器
3 :減算器
4 :コントローラ
5 :静的補償器
6 :動的補償器
7 :減算器
8 :減算器
9 :パラメータ調整器
10 :理想出力算出器
11 :メモリ
12 :検出器
50 :制御ループ
100 :制御対象
Gm :入出力モデル
:微分ゲイン
Kp :比例ゲイン
f0 :静的ゲイン

Claims (5)

  1. 制御対象を制御する制御装置であって、
    前記制御対象から出力される制御出力と目標出力との偏差を零にするための制御入力を算出するコントローラと、
    静的パラメータと前記制御入力とに基づいて、前記制御対象の静的特性の変動を補償する静的補償入力を算出する静的補償器と、
    動的パラメータと前記制御出力とに基づいて前記制御対象の動的特性の変動を補償する動的補償入力を算出する動的補償器と、
    前記静的補償入力と前記動的補償入力とを合成して実入力を算出し、前記制御対象に入力する合成器と、
    前記制御入力と前記制御出力との理想的な入出力関係を規定する入出力モデルを用いて、前記制御入力に対応する理想出力を算出する理想出力算出器と、
    前記制御出力と前記理想出力との差が最小化されるように前記静的パラメータ及び前記動的パラメータをそれぞれ調整するパラメータ調整器とを備える、
    制御装置。
  2. 前記動的パラメータは、比例ゲインと微分ゲインとを含み、
    前記動的補償器は、前記制御出力に前記比例ゲインを乗じた値と前記制御出力の微分値に前記微分ゲインを乗じた値との和を前記動的補償入力として算出し、
    前記パラメータ調整器は、前記制御出力と前記理想出力との差が最小化されるように前記静的パラメータ、前記前記比例ゲイン、及び前記微分ゲインをそれぞれ調整する、
    請求項1記載の制御装置。
  3. 前記パラメータ調整器は、逐次最小二乗法を用いて前記静的パラメータ及び前記動的パラメータを調整する、
    請求項1又は2記載の制御装置。
  4. 前記合成器は、前記静的補償入力から前記動的補償入力を減じた値を前記実入力として算出する、
    請求項1~3のいずれかに記載の制御装置。
  5. 前記制御対象は、油圧システムであり、
    前記油圧システムは、
    圧ポンプと、
    前記油圧ポンプを駆動する電動機と、
    前記油圧ポンプから供給される作動油によって作動するアクチュエータとを備え、
    前記油圧システムは、前記アクチュエータの状態値を前記制御出力として出力し、
    前記合成器は、前記電動機の指令値を前記実入力として算出する、
    請求項1~4のいずれかに記載の制御装置。
JP2020193678A 2020-11-20 2020-11-20 制御装置 Active JP7529545B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2020193678A JP7529545B2 (ja) 2020-11-20 2020-11-20 制御装置
CN202180076777.4A CN116457731A (zh) 2020-11-20 2021-09-29 控制装置
PCT/JP2021/035961 WO2022107456A1 (ja) 2020-11-20 2021-09-29 制御装置
EP21894331.4A EP4235316A4 (en) 2020-11-20 2021-09-29 Control device
US18/252,983 US20230418241A1 (en) 2020-11-20 2021-09-29 Control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020193678A JP7529545B2 (ja) 2020-11-20 2020-11-20 制御装置

Publications (2)

Publication Number Publication Date
JP2022082233A JP2022082233A (ja) 2022-06-01
JP7529545B2 true JP7529545B2 (ja) 2024-08-06

Family

ID=81708889

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020193678A Active JP7529545B2 (ja) 2020-11-20 2020-11-20 制御装置

Country Status (5)

Country Link
US (1) US20230418241A1 (ja)
EP (1) EP4235316A4 (ja)
JP (1) JP7529545B2 (ja)
CN (1) CN116457731A (ja)
WO (1) WO2022107456A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115157258A (zh) * 2022-07-25 2022-10-11 广东美的智能科技有限公司 多轴伺服系统及其控制方法、装置、主控制器、存储介质
JP7724557B2 (ja) * 2023-08-22 2025-08-18 国立大学法人広島大学 モデル誤差抑制補償器
CN119024680B (zh) * 2024-10-29 2025-01-24 中国科学院近代物理研究所 在线全局优化射频低电平系统反馈控制参数的方法及装置

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3545006B2 (ja) * 1992-04-06 2004-07-21 三菱電機株式会社 2自由度制御装置及び電動機のサーボ制御装置
JP2000307604A (ja) * 1999-04-23 2000-11-02 Nec Commun Syst Ltd Atmリンク切り換えシステム
JP2002312004A (ja) * 2001-04-18 2002-10-25 Nagoya Industrial Science Research Inst サーボ制御システム
JP2005330856A (ja) * 2004-05-19 2005-12-02 Denso Corp 自動車の制御装置
JP2006090206A (ja) * 2004-09-24 2006-04-06 Denso Corp 自動車の制御装置
CN100494698C (zh) * 2007-10-31 2009-06-03 东北大学 一种动态补偿液压伺服阀零漂的方法
JP5731933B2 (ja) * 2011-08-30 2015-06-10 川崎重工業株式会社 適応制御装置および適応制御方法ならびに射出成形機の制御装置および制御方法
JP2015148928A (ja) * 2014-02-06 2015-08-20 日立建機株式会社 建設機械
JP6279356B2 (ja) * 2014-03-10 2018-02-14 株式会社神戸製鋼所 作業機械の油圧駆動装置
US9809247B2 (en) * 2015-01-30 2017-11-07 Steering Solutions Ip Holding Corporation Motor control current sensor loss of assist mitigation for electric power steering
CN104898425B (zh) * 2015-05-19 2017-08-04 北京化工大学 基于静态前馈补偿的抗饱和内模控制系统设计方法
CN106773665A (zh) * 2015-11-23 2017-05-31 广东工业大学 一种基于李雅普诺夫模型的超声波电机自适应控制系统
GB2551701A (en) * 2016-06-21 2018-01-03 Univ Court Univ Of Edinburgh Control or processing system and method
RU2636603C1 (ru) * 2017-02-16 2017-11-24 Общество с ограниченной ответственностью "Ракурс-инжиниринг" Способ управления мощностью реактивных гидравлических турбин
US10181791B2 (en) * 2017-04-14 2019-01-15 Allegro Microsystems, Llc Converter digital control circuit with adaptive compensation
US11254416B2 (en) * 2019-05-16 2022-02-22 The Boeing Company Method and apparatus for remote angle calibration and measurement
CN110297425B (zh) * 2019-07-16 2022-09-16 华侨大学 一种参数带宽化与能量化的自适应抗扰控制方法
CN110543097B (zh) * 2019-08-06 2023-01-20 广东工业大学 一种基于模型参考自适应的压电陶瓷驱动器控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
尾島 康仁,"微分要素を含む非線形システムに対するデータ駆動型制御系の一設計",電気学会研究会資料 The Papers of Technical Meeting on "Control",IEE Japan,日本,一般社団法人電気学会,2016年04月04日,第1-6頁

Also Published As

Publication number Publication date
CN116457731A (zh) 2023-07-18
EP4235316A1 (en) 2023-08-30
WO2022107456A1 (ja) 2022-05-27
EP4235316A4 (en) 2024-05-08
US20230418241A1 (en) 2023-12-28
JP2022082233A (ja) 2022-06-01

Similar Documents

Publication Publication Date Title
JP7529545B2 (ja) 制御装置
CN100562823C (zh) 一种复杂轨迹的轮廓控制方法
JP6165585B2 (ja) プラント制御装置の制御ゲイン最適化システム
Yao et al. Nonlinear adaptive robust backstepping force control of hydraulic load simulator: Theory and experiments
CN118971207B (zh) 一种涉网火电机组中调速系统的参数整定方法及系统
CN110023857A (zh) 伺服控制装置
JP2022112300A (ja) パラメータ調整装置及びパラメータ調整方法
JP2019053594A (ja) 制御装置及び機械学習装置
CN109324503B (zh) 基于鲁棒积分的多层神经网络电机系统控制方法
US9683462B2 (en) Cooling fluid flow control system for steam turbine system and program product
JP7003341B1 (ja) ダイクッション制御装置、ダイクッション制御方法、およびダイクッション制御プログラム
CN120701640A (zh) 一种适用复杂情况的液压升降系统同步控制系统及方法
WO2019053826A1 (ja) 圧延ラインの数学モデル算出装置および制御装置
Aguilar-Avelar et al. A feedback linearization controller for trajectory tracking of the Furuta pendulum
Abbott et al. Application of an optimal control synthesis strategy to an electro-hydraulic positioning system
Wang et al. Closed-loop simulation integrating finite element modeling with feedback controls in powder bed fusion additive manufacturing
JP4189716B2 (ja) モータ速度制御装置
Wos et al. Self-tuning controllers based on polynomial methods for electro-hydraulic servo drive
CN114755914A (zh) 基于imc-pid的航空发动机伺服控制器设计方法
JP2881873B2 (ja) デジタル適応制御装置
JP2565600B2 (ja) 圧延機に於ける蛇行制御方法
KR20220122254A (ko) 플랜트 제어 장치 및 그 제어 방법
CN111542789B (zh) 用设定值加权来闭环控制控制器的方法
KR950013243B1 (ko) 압연기 제어방법 및 장치
JP4036520B2 (ja) 単純適応制御装置

Legal Events

Date Code Title Description
RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7426

Effective date: 20201204

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230720

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240326

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240405

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: 20240716

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240725

R150 Certificate of patent or registration of utility model

Ref document number: 7529545

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150