JP7508920B2 - Control device and design method - Google Patents
Control device and design method Download PDFInfo
- Publication number
- JP7508920B2 JP7508920B2 JP2020125530A JP2020125530A JP7508920B2 JP 7508920 B2 JP7508920 B2 JP 7508920B2 JP 2020125530 A JP2020125530 A JP 2020125530A JP 2020125530 A JP2020125530 A JP 2020125530A JP 7508920 B2 JP7508920 B2 JP 7508920B2
- Authority
- JP
- Japan
- Prior art keywords
- state
- command value
- control device
- controlled object
- trolley
- 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
Links
- 238000000034 method Methods 0.000 title claims description 23
- 238000013461 design Methods 0.000 title description 5
- 238000011156 evaluation Methods 0.000 claims description 22
- 230000014509 gene expression Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000002123 temporal effect Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 18
- 230000010365 information processing Effects 0.000 description 14
- 238000012545 processing Methods 0.000 description 7
- 238000005457 optimization Methods 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Landscapes
- Control And Safety Of Cranes (AREA)
- Feedback Control In General (AREA)
Description
本発明は、制御装置、及び設計方法に関する。 The present invention relates to a control device and a design method.
制御対象を制御するシステムとしては、例えば、フィードバック制御器やフィードフォワード制御器を用いる制御システムがある(例えば、特許文献1参照)。 Examples of systems that control a controlled object include control systems that use feedback controllers and feedforward controllers (see, for example, Patent Document 1).
特許文献1には、例えば、クレーンの吊荷の旋回振れを低減するフィードバック制御器の設計指針の一例は開示されているが、フィードフォワード制御器に関する設計指針は開示されていない。 For example, Patent Document 1 discloses an example of a design guideline for a feedback controller that reduces the swing of a load suspended by a crane, but does not disclose any design guideline for a feedforward controller.
本発明は、上記のような従来の問題に鑑みてなされたものであって、制御対象の状態を変化させる際にフィードフォワード制御器に適した指令値を出力させることができる制御装置を提供することを目的とする。 The present invention was made in consideration of the above-mentioned problems in the conventional technology, and aims to provide a control device that can output a command value suitable for a feedforward controller when changing the state of the controlled object.
前述した課題を解決する本発明の第1の態様は、制御対象を制御するための第1指令値を出力するフィードフォワード制御器と、前記制御対象の状態を第1状態から第2状態へと変化させる際に、前記制御対象のモデルを含む制約条件を用いて計算された、前記第1指令値の時間変化に関する第1プロファイルを記憶する記憶装置と、を備え、前記フィードフォワード制御器は、前記制御対象の状態を前記第1状態から前記第2状態へと変化させる際、前記第1プロファイルに基づいて、時間的に変化する前記第1指令値を出力する、制御装置である。 A first aspect of the present invention that solves the above-mentioned problems includes a feedforward controller that outputs a first command value for controlling a control object, and a storage device that stores a first profile of the time change of the first command value calculated using a constraint condition including a model of the control object when changing the state of the control object from a first state to a second state, and the feedforward controller is a control device that outputs the first command value that changes over time based on the first profile when changing the state of the control object from the first state to the second state.
本発明の第2の態様は、制御対象を制御するフィードフォワード制御器の指令値を計算する計算方法であって、前記制御対象の状態を第1状態から第2状態へと変化させる際、前記制御対象を示す状態空間モデルの状態変数と、前記指令値と、を決定変数とする評価関数を定式化するステップと、所定の制約条件を有する前記評価関数を用いて、前記制御対象を前記第1状態から前記第2状態へと変化させる際の前記指令値を計算するステップと、を含む計算方法。 A second aspect of the present invention is a calculation method for calculating a command value of a feedforward controller that controls a controlled object, the calculation method including the steps of: formulating an evaluation function in which, when changing the state of the controlled object from a first state to a second state, the command value and state variables of a state space model representing the controlled object are decision variables; and calculating the command value when changing the controlled object from the first state to the second state using the evaluation function having a predetermined constraint condition.
本発明によれば、制御対象の状態を変化させる際にフィードフォワード制御器に適した指令値を出力させることができる技術を提供することができる。 The present invention provides a technology that can output a command value appropriate for a feedforward controller when changing the state of a controlled object.
本明細書及び添付図面の記載により、少なくとも以下の事項が明らかとなる。 The following points become clear at least from the description in this specification and the accompanying drawings.
<<<クレーン10の概要について>>>
図1は、走行(x軸)、横行(y軸)、巻上(z軸)の3軸の動きが可能な天井クレーン10(以下、単にクレーンと称する。)の概要を示す図である。クレーン10は、全体がx軸方向に移動可能であるとともに、ガーダ(不図示)に沿ってy軸方向に移動可能なトロリー20と、トロリー20からの巻上げ、または巻下げ可能なワイヤに吊り下げられた吊荷21とを含む。
<<<<Overview of the
1 is a diagram showing an overview of an overhead crane 10 (hereinafter simply referred to as a crane) capable of movement along three axes: traveling (x-axis), traversing (y-axis), and hoisting (z-axis). The
以下、本実施形態では、トロリー20の座標を(x,y,0)で表現し、吊荷21の座標を(x~,y~,z~)で表現する。なお、「x~」、「y~」、「z~」の夫々は、下記の記号(1)~(3)と同じである。また、吊荷21の振れ幅θは、吊荷21を、xy平面に投影した角度φ、xz平面に投影した角度α、yz平面に投影した角度βで評価される。
<<<トロリー20のモデリング>>>
Hereinafter, in this embodiment, the coordinates of the
<<<Modeling of Trolley 20>>>
図2は、トロリー20及び吊荷21の構成の一例を示す図である。なお、トロリー20には、例えば、フィードバック制御器や加算器等が含まれているが、図2では、便宜上、トロリー20、及び吊荷21をモデリングする際の最低限の構成のみを図示している。
Figure 2 shows an example of the configuration of the
トロリー20は、フィードフォワード制御器30、及びモータ40を含んで構成される。フィードフォワード制御器30は、x軸方向の指令値ux、y軸方向の指令値uy、z軸方向の指令値uzをモータ40に対して出力する。なお、以下、x軸~y軸の3軸の指令値をまとめて、指令値uと記載する。
The
モータ40は、指令値uxに基づいて、クレーン10自体をx軸方向に移動させる。また、モータ40は、指令値uyに基づいて、トロリー20の車輪を回転させることにより、トロリー20をガーダ(不図示)に沿ってy軸方向に移動させる。さらに、モータ40は、指令値uzに基づいて、吊荷21が接続されたワイヤを、z軸方向に巻上げ、または巻下げを行う。
The
そして、トロリー20が、フィードフォワード制御器30からの指令値uで動作する場合、以下の式(1)~(3)が成立する。なお、以下、本実施形態では、「トロリー20の移動」には、クレーン10自体が移動することと、トロリー20がガーダを移動することを含むこととする。
なお、ここでは、クレーン10及びトロリー20は、一次遅れの特性を有し、モータ40のゲインを「K」、モータ40の時定数を「T」としている。このため、例えば、x軸方向のモータ40のゲインは、Kxとなり、x軸方向のモータ40の時定数は、Txと表現できる。また、式(1)の左辺の項は、「x」を時間tで2階微分した表現であり、右辺のxに関する記号は、「x」を時間tで1階微分した表現である。
When the
Here, the
<<<吊荷21のモデリング>>>
吊荷21に関する運動方程式を変形し、吊荷21の質量を「m」、重力加速度を「g」、x方向のダンピング係数を「Dα」、x方向のダンピング係数を「Dβ」とすると、吊荷21には、以下の式(4)、(5)が成立する。
<<<状態空間モデル>>>
式(1)~(5)は2階微分を含む式であるため、クレーン10の解析を容易に行うべく、ここでは、一階微分のみが含まれる表現に変換する。具体的には、式(1)~(5)を、以下の式(6)のように表現する。
このように、フィードフォワード制御器30によって制御されるクレーン10は、式(7)に示したような状態空間モデルで表現することができる。
<<<Modeling of Suspended
By transforming the equation of motion for the suspended
<<<<State Space Model>>>
Since the formulas (1) to (5) include second-order differentials, they are converted here into expressions including only first-order differentials in order to facilitate the analysis of the
In this way, the
<<<トロリー20の移動について>>>
図3は、静止状態のトロリー20を、位置Aから、搬送先である位置Dへと移動させる際の加減速パターンの一例を示す図である。例えば、トロリー20は、位置Aから、位置Bへと加速運転し、位置Bから位置Cへと等速運転し、位置Cから位置Dへと減速運転する。なお、位置Aの静止状態は、「第1状態」に相当し、位置Bの等速状態は、「第2状態」に相当し、静止状態及び等速状態はともに「定常状態」に相当する。
<<<<About the movement of the
3 is a diagram showing an example of an acceleration/deceleration pattern when the
ところで、図3に示すような加減速パターンでトロリー20を移動させる際、トロリーを、順次加速、等速、減速させる指令値uが、フィードフォワード制御器30から出力される必要がある。ただし、一般に、トロリー20を、所定の速度まで加速させる指令値u、または所定の速度から減速させる指令値uには無数のパターンがあるため、一意に指令値uを定めることは難しい。
Incidentally, when moving the
そこで、詳細は後述するが、本実施形態では、例えば、トロリー20が位置Bに到達した際に、吊荷21の振れ角や、振れ角速度がゼロになる等の各種制約条件の下、最適化問題を解く手法を用いて指令値uを定める。さらに、本実施形態では、最適化問題を解く手法以外の例として、トロリー20の物理特性と終端条件(制約条件)から指令値uの時間多項式を導出する方法を用いる。以下、それぞれの手法を実行する情報処理装置について説明する。
Thus, although details will be described later, in this embodiment, the command value u is determined using a method for solving an optimization problem under various constraint conditions, such as the sway angle and sway angular velocity of the suspended
<<<情報処理装置について>>>
図4は、情報処理装置50のハードウェア構成の一例を示す図である。情報処理装置50は、CPU(Central Processing Unit)60、メモリ61、記憶装置62、入力装置63、表示装置64、及び通信装置65を含むコンピュータである。
<<<<About information processing devices>>>
4 is a diagram showing an example of a hardware configuration of the
CPU60は、メモリ61や記憶装置62に格納されたプログラムを実行することにより、情報処理装置50における様々機能を実現する。
The
メモリ61は、例えばRAM(Random-Access Memory)等であり、プログラムやデータ等の一時的な記憶領域として用いられる。
記憶装置62は、CPU60によって実行されるプログラムやデータセット等の各種情報を格納する不揮発性の記憶装置である。
The
入力装置63は、ユーザによるコマンドやデータの入力を受け付ける装置であり、キーボード、タッチパネルディスプレイ上でのタッチ位置を検出するタッチセンサなどの入力インタフェースを含む。
The
表示装置64は、例えばディスプレイなどの装置であり、通信装置65は、ネットワークを介して他のコンピュータと各種プログラムやデータの受け渡しを行う。
==情報処理装置50の処理について==
The
Regarding the Processing of
<<最適化問題について>>
図5は、情報処理装置50のCPU60が実行する処理S10の一例を示すフローチャートである。ここでは、CPU60は、図3で示した、トロリー20を位置Aから位置Bへ移動する際、最適化問題を解くことにより状態変数xp及び指令値uを計算する。
<<About the optimization problem>>
5 is a flowchart showing an example of processing S10 executed by the
まず、CPU60は、10個の状態変数xpと、3軸の指令値uとの二乗和を示す、式(8)の評価関数を定式化する(S20)。
ここで、式(9)及び式(10)は、式(8)の評価関数から決定変数を求める際の制約条件である。なお、式(9)は、クレーン10の状態空間モデルを示す式であり、式(10)は、状態変数xpの各種条件を示す式である。
Here, equations (9) and (10) are constraint conditions when determining the decision variables from the evaluation function of equation (8). Note that equation (9) is an equation showing the state space model of the
式(9)において、xp(t)は、時刻tにおける状態ベクトルであり、xp(t+1)は、時刻tからは、1サンプル後の時刻の状態ベクトルである。また、A,Bは、上述した行列であり、u(t)は、時刻tにおける指令値(指令値ベクトル)である。 In equation (9), xp(t) is the state vector at time t, and xp(t+1) is the state vector one sample after time t. A and B are the matrices described above, and u(t) is the command value (command value vector) at time t.
式(10)に記載された条件としては、トロリー20が位置Aから、位置Bへ移動する際に許容される吊荷21のx軸方向の振れ幅の上下限値、x軸方向の振れは幅速度の(dΔx/dt)の上下限値が含まれる。また、式(10)には、吊荷21のy軸方向の振れは幅(Δy)の上下限値、y軸方向の振れは幅速度の(dΔy/dt)の上下限値が含まれる。さらに、式(10)は、状態ベクトルxpの終端条件と、トロリー20の指令値uの上下限値を含む。ここで、状態ベクトルxpの終端条件とは、例えば、トロリー20が位置Bに到達した際の状態ベクトルxp(例えば、トロリー20の位置、速度等)を含む。
The conditions described in formula (10) include the upper and lower limits of the swing width in the x-axis direction of the
なお、本実施形態では、上述した各種上下限値以外に、トロリー20が位置Bに到達した際に、吊荷21の振れ角や、振れ角速度がゼロになること、トロリー20が位置Bの等速運転時の速度となること、を制約条件として含む。
In addition to the various upper and lower limit values described above, this embodiment includes constraints such that when the
そして、CPU60は、式(9)及び式(10)の制約条件を用いて、式(8)の評価関数が最小となる決定変数(指令値u、状態変数xp)を計算する(S21:第2ステップ)。この結果、例えば、図6に示すよう、指令値u、状態変数xpの夫々が計算されることになる。したがって、図2に示すフィードフォワード制御器30が、図6に示す指令値uを出力すると、トロリー20の状態変数xp(例えば、トロリー20の速度、吊荷21の振れ幅、振れ幅速度)は、図6に示すよう変化する。
Then, the
ところで、図6は、位置Aの静止状態のトロリー20を加速させ、位置Bにおいて等速状態とするための指令値uの波形であるが、本実施形態では、図3に示すよう、トロリー20を位置Aから、位置Dへと移動させる。そして、位置Bから、位置Cまでは、トロリー20は等速運転を行うため、図7に示すよう、位置Bにおける指令値uが出力され続ける。なお、指令値uは、図6に示すよう、直線状に変化する波形ではないが、図7は、便宜上指令値uを直線として描いている。
FIG. 6 shows the waveform of the command value u for accelerating the
また、位置Cの等速状態(第2状態)から、トロリー20を減速させ、位置Dにて静止状態(第1状態)とするため、位置Cから位置Dにおいては、位置Aから位置Bへの指令値uの波形を反転した指令値uを、フィードフォワード制御器30は出力する。フィードフォワード制御器30が、このような指令値uを出力することにより、吊荷21の振れ幅、振れ幅速度を抑制しつつ、トロリー20を位置Aから位置Dへと適切に移動させることができる。なお、詳細は後述するが、計算された状態変数xpは、フィードバック制御器の目標値として用いることができる。
In addition, in order to decelerate the
<<最適化問題の解の近似式>>
図8は、情報処理装置50のCPU60が実行する処理S11の一例を示すフローチャートである。処理S11において、CPU60は、図5の処理S20,S21に加え、計算された状態変数xp、及び指令値uの近似式を計算する(S22)。例えば、CPU60は、図9に示すように、指令値uxの一点鎖線で示す近似式(関数)を、例えば最小二乗法を用いて計算する。
<<Approximate solution to optimization problem>>
Fig. 8 is a flow chart showing an example of processing S11 executed by the
したがって、フィードフォワード制御器30は、処理S22で計算された指令値uの近似式に基づいて、指令値uを出力することにより、吊荷21の振れ幅、振れ幅速度を抑制しつつ、トロリー20を位置Aから位置Dへと適切に移動させることができる。このような近似式を用いることにより、例えば、位置Aから位置Dへの指令値uに関する情報のデータサイズを小さくすることができる。
Therefore, the
<<近似式を用いた評価関数>>
図8では、評価関数の最適解を近似することとしたが、例えば、時間多項式で表現された指令値uと、状態変数xpとを決定変数として評価関数を定式化しても良い。具体的には、まず指令値uとして、式(11)の時間多項式を用いる。
8, the optimal solution of the evaluation function is approximated, but the evaluation function may be formulated by using, for example, the command value u expressed by a time polynomial and the state variable xp as decision variables. Specifically, first, the time polynomial of Equation (11) is used as the command value u.
そして、図10の処理S12において、CPU60は、式(11)で表現された指令値uと、10個の状態変数xpとの二乗和を示す、式(8)の評価関数を定式化する(S30)。なお、ここでは、決定変数は、xp、u、ai、bi、ci、nとなる。そして、CPU60は、式(8)の評価関数において、式(9)~(11)と、多項式の次数nの上下限値とを制約条件とし、決定変数を計算する(S31)。なお、多項式の次数nの上下限値は、以下の様に表現され、多項式のai、bi、ci、nは、「時間多項式の複数の変数」に相当する。
Then, in process S12 of Fig. 10, the
この結果、フィードフォワード制御器30は、処理S31で計算された指令値uの時間多項式に基づいて、指令値uを出力することにより、吊荷21の振れ幅、振れ幅速度を抑制しつつ、トロリー20を位置Aから位置Dへと適切に移動させることができる。
As a result, the
<<物理特性及び終端条件に基づく指令値>>
図5、図8、及び図10の処理S10~12では、最適化問題を解くことにより、指令値uを計算しているがこれに限られない。例えば図11に示すよう、CPU60は、クレーン10の運動方程式(つまり、動特性モデル)と、制約条件とを用いて、運動方程式の複数の解のうち、最小次数で表される時間多項式を解として計算する(S40)。
<<Command value based on physical characteristics and terminal conditions>>
In steps S10 to S12 in Fig. 5, Fig. 8, and Fig. 10, the command value u is calculated by solving an optimization problem, but this is not limited to this. For example, as shown in Fig. 11, the
なお、ここで、制約条件とは、いわゆる境界条件であり、例えば、位置Bにおけるモータ40の3軸方向の速度、吊荷21の3軸方向の振れ角、振れ角速度を含む。また、処理S40で実行される具体的な計算方法は、例えば、「よくわかる機械の制振設計:三好孝典 日刊工業新聞社」の第4章や「境界条件を考慮した有限時間多項式表現による制御入力の解析解導出、計測自動制御学会論文集:三好 孝典、寺嶋 一彦、Vol.36,No.6. pp.489-496、2000年」に開示されている。
The constraint conditions are so-called boundary conditions, and include, for example, the speeds of the
=====クレーン10の制御装置=====
図12は、クレーン10の構成の一例を示す図である。クレーン10は、モータ100、センサ110、及び制御装置120を含んで構成される。
=====
12 is a diagram showing an example of the configuration of the
モータ100は、上述したモータ40と同じであり、入力される指令値x軸、y軸方向にトロリー20を移動させるとともに、吊荷21の巻上げ、巻下げを行う。
The
センサ110は、例えばトロリー20のx軸の位置、y軸の位置、x軸方向の速度、y軸方向の速度、吊荷21のz軸の位置、z軸方向の速度等、10個の状態変数の夫々に対応する値を出力する。なお、本実施形態においては、センサ110が出力する10個の情報(位置x,y、速度dx/dt,dy/dt、・・)を「出力xps」とする。ここで、センサ110は、10個の速度等を出力することとしたが、速度や振れ角は、センサ110から取得せずに、他のセンサ(不図示)のセンサ情報から演算により求めてもよい。
The
制御装置120は、モータ100の動作を制御する装置であり、記憶装置200、フィードフォワード制御器(FF制御器)201、加算器202,204、及びフィードバック制御器(FB制御器)203を含んで構成される。
The
記憶装置200は、処理S10~S12の何れかで計算された指令値uの時間的なプロファイル(または近似式)を示す情報I1と、情報I1に対応する状態変数xpの時間的なプロファイル(または近似式)を示す情報I2とを記憶する。なお、以降、処理S10~S12の何れかで計算された指令値uを、フィードフォワード制御器201の「指令値uff」とする。また、指令値uffの時間的なプロファイルを示す情報I1は、「第1プロファイル」に相当し、状態変数xpの時間的なプロファイルを示す情報I2は、「第2プロファイル」に対応する。さらに、情報I1は、指令値uffの「時間変化に関する情報」であり、情報I2は、状態変数xpの「時間変化に関する情報」である。
The
フィードフォワード制御器201は、情報I1に基づいて、指令値uffを出力する。加算器202(第1加算器)は、情報I2から、センサ110からの10個の出力xpsに対応する10個の状態変数xpを選択する。そして加算器202は、選択した状態変数xpを目標値とし、目標値と、出力xpsとの差eを計算する。
The
フィードバック制御器203は、誤差eに基づいて、モータ40を駆動するための指令値ufbを出力する。
The
加算器204(第2加算器)は、フィードフォワード制御器201からの指令値uff(第1指令値)と、フィードバック制御器203からの指令値ufb(第2指令値)とを加算し、加算結果を指令値utとしてモータ100に出力する。この結果、モータ100は、指令値utに基づいて駆動される。
The adder 204 (second adder) adds the command value uff (first command value) from the
ところで、例えば、トロリー20を位置Aから位置Bへ移動させる際に用いられる情報I1,I2が、処理S10により計算された時間的なプロファイルに基づく場合、情報I1,I2は、図6で示した波形になる。この場合、フィードフォワード制御器201は、図6に示した指令値uffを出力する。そして、加算器202への目標値として、図6で示した10個の出力(トロリーや吊荷の位置x,y,z、速度dx/dt,dy/dt,dz/dt、振れ幅Δx,Δy、振れ幅速度dΔx/t,dΔy/t)が入力される。
For example, when the information I1 and I2 used when moving the
ここで、センサ110の出力xpsが、図6で示した4つの出力と一致している場合には、差eはゼロとなる。したがって、この場合には、フィードバック制御器203からの指令値ufbもゼロとなり、トロリー20は、指令値uffのみに基づいて動作することになる。つまり、指令値utは、指令値uffと等しくなる。
Here, if the output xps of the
一方、センサ110の出力xpsと、図6で示した10個の出力とがずれており、差eがゼロでない場合、フィードバック制御器203は、差eに応じた指令値ufbを出力する。この結果、トロリー20は、図6の波形(位置等)からのずれが小さくなるよう、フィードバック制御器203により制御される。したがって、いずれの場合であっても、トロリー20は、図6に示した位置、速度を保ちつつ、位置Aから位置Bへと移動することになる。
On the other hand, if the output xps of the
なお、位置Bから位置Cへの移動、及び位置Cから位置Dへの移動についても、位置Aから位置Bの移動と同様であるため、本実施形態では、吊荷21の振れ角、振れ角速度を抑制しつつ、適切にトロリー20を移動させることができる。また、ここでは、例えば処理S10により計算された時間的なプロファイルが、情報I1,I2である場合を説明したが、処理S11,S12により計算された結果(近似式)が、情報I1,I2である場合も同様である。
Note that movement from position B to position C and movement from position C to position D are similar to movement from position A to position B, so in this embodiment, the
<<<制御装置150>>>
図13は、クレーン11の実施例を示す図である。図13のクレーン11には、制御装置120の代わりに制御装置150が設けられている。なお、図12と、図13とでは同じ符号が付されたブロックは同じである。
<<<
Fig. 13 is a diagram showing an embodiment of the
制御装置150は、モータ100の動作を制御する装置であり、記憶装置200、フィードフォワード制御器201、加算器202,204、フィードバック制御器203、及び選択装置205を含んで構成される。図12と、図13とで同じ符号が付されたブロックは同じであるため、ここでは、選択装置205について説明する。
The
選択装置205は、選択情報に基づいて、トロリー20の位置、速度等の目標を示す目標値r、または、記憶装置200に記憶された情報I1,I2の何れかを出力する。選択装置205が目標値rを出力する場合、フィードフォワード制御器201、及びフィードバック制御器203は、トロリー20の目標値rに応じた指令値uff、ufbを出力する。この結果、トロリー20は、指定された目標値rで定められた動作をすることになる。
Based on the selection information, the
一方、選択装置205が情報I1,I2を出力する場合、上述したよう、吊荷21の振れ角等が小さくなるよう、フィードフォワード制御器201、フィードバック制御器203が指令値uff、ufbを出力する。このような選択装置205が設けられることにより、利用者は、柔軟にトロリー20の移動を制御することができる。
On the other hand, when the
<<その他>>
本実施形態では、制御対象の一例としてクレーン10について説明したが、制御対象は、フィードフォワード制御器により制御される対象であればクレーンに限られない。例えば、制御対象は、フィードフォワード制御器により庫内の温度を一定に保つ冷凍庫のコンプレッサーであっても良い。このような場合、フィードフォワード制御器は、例えば、庫内の温度を、所定の温度Ta(第1状態)から温度Tb(第2状態)に変化させるよう、温度制御器を制御する。そして、庫内の温度の上下限値等を制約条件とすることで、本実施形態と同様の手法を用いてフィードフォワード制御器の指令値を計算することができる。
<<Others>>
In this embodiment, the
=====まとめ=====
以上、本発明の一実施形態であるクレーン10について説明した。制御装置120のフィードフォワード制御器201は、情報I1に基づいて、トロリー20を静止状態(位置A)から等速状態(位置B)に変化させる際に、例えば吊荷21の振れ角等を抑制できる指令値uを出力することができる。また、フィードフォワード制御器201は、予め記憶装置200に格納された情報I1を用いているため、演算処理量が少なく、制御装置120のCPU(不図示)における計算負荷を小さくできる。
== ...
The
また、フィードフォワード制御器201は、位置Bから位置Cにおいて、位置Bにおける指令値uffを出力し続けることにより、トロリー20(移動する装置)を等速運転させることができる。
In addition, the
また、フィードフォワード制御器201は、位置Cから位置Dにおいて、位置Aから位置Bまでの指令値uffの波形を反転し、出力する。これにより、本実施形態では、吊荷21の振れ角等を抑制することができる。
In addition, the
また、制御装置120には、センサ110からの出力xpsと、出力xpsに応じた情報I2の値(目標値)との差eに基づいて、指令値ufbを出力するフィードバック制御器203が設けられている。したがって、フィードフォワード制御器201からの指令値uffに基づいて動作するトロリー20の位置等が、目標値(例えば、図6の状態変数xpの値)とずれた場合、フィードバック制御器203がずれを補正することができる。この結果、トロリー20を移動させる際、吊荷21の振れ角等を精度良く抑制することができる。
The
また、情報I1,I2として、例えば、図5に示すように、指令値uff及び状態変数xpを決定変数とする、式(8)の評価関数の最適解を用いることができる。 In addition, as information I1 and I2, for example, as shown in FIG. 5, the optimal solution of the evaluation function of equation (8) in which the command value u ff and the state variable x p are decision variables can be used.
また、情報I1,I2として、例えば、図8に示すように、指令値uff及び状態変数xpを決定変数とする、式(8)の評価関数の最適解の近似式を用いることができる。 In addition, as information I1 and I2, for example, as shown in FIG. 8, an approximation equation of the optimal solution of the evaluation function of equation (8) can be used, in which the command value u ff and the state variable x p are decision variables.
また、情報I1,I2として、例えば、図10に示すように、時間多項式で表現された指令値uff及び状態変数xpを決定変数とする、式(8)の評価関数の最適解を用いることができる。 In addition, as information I1 and I2, for example, as shown in FIG. 10, the optimal solution of the evaluation function of equation (8) can be used, in which the command value u ff and the state variable x p expressed as a time polynomial are the decision variables.
また、例えば、図11に示すように、評価関数を定式化せず、例えばクレーン10の運動方程式の複数の解のうち、最小次数で表現される時間多項式を解としても良い。
Also, for example, as shown in FIG. 11, the evaluation function may not be formulated, and the solution may be, for example, a time polynomial expressed with the smallest degree among multiple solutions of the equation of motion of the
また、本実施形態では、クレーン10に用いられるフィードフォワード制御器について詳細を説明したが、制御対象は、クレーン10に限られない。具体的には、サーボモータや恒温槽のコンプレッサー等、フィードフォワード制御器が用いられる制御対象に適用することができる。
In addition, in this embodiment, the feedforward controller used in the
また、本実施形態では、例えば、処理S10~S12に示したように、状態変数xp及び指令値uffを決定変数とする評価関数を用いて、指令値uff(つまり、処理S10~S12の指令値u)を計算することができる。なお、指令値uffは、例えば、評価関数の最適解や、最適解を近似した式であっても良い。 In this embodiment, for example, as shown in steps S10 to S12, the command value uff (i.e., the command value u in steps S10 to S12) can be calculated using an evaluation function in which the state variable xp and the command value uff are decision variables. Note that the command value uff may be, for example, an optimal solution of the evaluation function or an equation that approximates the optimal solution.
上記の実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。また、本発明は、その趣旨を逸脱することなく、変更や改良され得るとともに、本発明にはその等価物が含まれるのはいうまでもない。 The above embodiments are intended to facilitate understanding of the present invention, and are not intended to limit the scope of the present invention. Furthermore, the present invention may be modified or improved without departing from the spirit of the present invention, and it goes without saying that the present invention includes equivalents.
10,11 クレーン
20 トロリー
21 吊荷
30,201 フィードフォワード制御器(FF制御器)
40,100 モータ
50 情報処理装置
60 CPU
61 メモリ
62,200 記憶装置
63 入力装置
64 表示装置
65 通信装置
110 センサ
120,150 制御装置
202,204 加算器
203 フィードバック制御器(FB制御器)
205 選択装置(SEL)
10, 11
40, 100
61
205 Selection device (SEL)
Claims (9)
前記制御対象の状態を第1状態から第2状態へと変化させる際に、前記制御対象のモデルを含む制約条件を用いて計算された、前記第1指令値の時間変化に関する第1プロファイルを記憶する記憶装置と、
前記制御対象からの出力と、目標値との差に応じた第2指令値を出力するフィードバック制御器と、
前記第1指令値と、前記第2指令値とを加算し、加算結果に基づいて前記制御対象を制御する加算器と、
を備え、
前記フィードフォワード制御器は、前記制御対象の状態を前記第1状態から前記第2状態へと変化させる際、前記第1プロファイルに基づいて、時間的に変化する前記第1指令値を出力し、
前記制御対象は、前記第1指令値に基づいて移動する装置であり、
前記第1状態は静止状態、第2状態は、等速運転状態であり、
前記フィードフォワード制御器は、前記制御対象の状態が前記等速運転状態になると、前記第1プロファイルに基づいて、前記等速運転状態の前記第1指令値を出力し、
前記制約条件は、前記制御対象の状態空間モデルと、前記状態空間モデルの状態変数の条件を含み、
前記記憶装置は、前記制御対象の状態を前記静止状態から前記等速運転状態へと変化させる際に、前記制約条件を用いて計算された、前記状態変数の時間的な変化を示す第2プロファイルを更に記憶し、
前記フィードバック制御器は、前記第2プロファイルに基づいて、前記状態変数のうち、前記制御対象からの出力に対応する値を、前記目標値として出力する、
制御装置。 a feedforward controller that outputs a first command value for controlling a controlled object;
a storage device that stores a first profile regarding a time change of the first command value calculated using a constraint condition including a model of the control object when changing a state of the control object from a first state to a second state;
a feedback controller that outputs a second command value according to a difference between an output from the controlled object and a target value;
an adder that adds the first command value and the second command value and controls the controlled object based on a result of the addition;
Equipped with
the feedforward controller outputs the first command value which changes over time based on the first profile when changing a state of the controlled object from the first state to the second state;
the controlled object is a device that moves based on the first command value,
The first state is a stationary state, and the second state is a constant speed operation state,
the feedforward controller outputs, when a state of the controlled object becomes the constant speed operation state, the first command value for the constant speed operation state based on the first profile;
the constraint condition includes a state space model of the controlled object and a condition on a state variable of the state space model;
the storage device further stores a second profile indicating a change over time of the state variable calculated using the constraint condition when the state of the controlled object is changed from the stationary state to the constant speed driving state;
the feedback controller outputs, as the target value, a value of the state variable corresponding to an output from the controlled object based on the second profile.
Control device.
前記フィードフォワード制御器は、前記制御対象の状態を前記等速運転状態から前記静止状態へ変化させる際、前記第1プロファイルに基づいて、前記制御対象の状態を前記静止状態から前記等速運転状態へと変化させる際の前記第1指令値の時間的な波形が反転するよう、前記第1指令値を出力する、
制御装置。 The control device according to claim 1 ,
the feedforward controller outputs the first command value based on the first profile when changing the state of the controlled object from the constant speed operation state to the stationary state such that a temporal waveform of the first command value when changing the state of the controlled object from the stationary state to the constant speed operation state is inverted.
Control device.
前記制御対象からの出力と、目標値との差に応じた第2指令値を出力するフィードバック制御器と、
前記第1指令値と、前記第2指令値とを加算し、加算結果に基づいて前記制御対象を制御する加算器と、
を更に備える
制御装置。 The control device according to claim 1 or 2 ,
a feedback controller that outputs a second command value according to a difference between an output from the controlled object and a target value;
an adder that adds the first command value and the second command value and controls the controlled object based on a result of the addition;
The control device further comprises:
前記第1及び第2プロファイルは、前記第1指令値及び前記状態変数を決定変数とした評価関数の最適解である、
制御装置。 The control device according to claim 1 ,
The first and second profiles are optimal solutions of an evaluation function in which the first command value and the state variables are decision variables.
Control device.
前記第1及び第2プロファイルは、前記第1指令値及び前記状態変数を決定変数とした評価関数の最適解の近似式である、
制御装置。 The control device according to claim 1 ,
the first and second profiles are approximation expressions of an optimal solution of an evaluation function in which the first command value and the state variables are decision variables;
Control device.
前記第1及び第2プロファイルは、複数の変数を含む時間多項式で表現された前記第1指令値と、前記状態変数と、を決定変数とした評価関数の最適解である、
制御装置。 The control device according to claim 1 ,
The first and second profiles are optimal solutions of an evaluation function in which the first command value expressed by a time polynomial including a plurality of variables and the state variables are decision variables.
Control device.
前記制約条件は、前記モデルに対応する前記制御対象の運動方程式を含み、
前記第1プロファイルは、前記運動方程式の複数の解のうち最小の次数を有する時間多項式で表現される解である、
制御装置。 The control device according to claim 1 ,
the constraint condition includes an equation of motion of the control object corresponding to the model;
The first profile is a solution expressed by a time polynomial having a minimum degree among a plurality of solutions of the equation of motion.
Control device.
前記制御対象は、トロリー及び吊荷を有するクレーンであり、
前記制約条件は、前記トロリーの速度に応じた前記第1指令値に関する条件と、前記吊荷の振れ幅に関する条件と、前記吊荷の振れ幅速度に関する条件と、を含む、
制御装置。 The control device according to any one of claims 1 to 7 ,
the control object is a crane having a trolley and a load,
The constraint conditions include a condition regarding the first command value according to the speed of the trolley, a condition regarding the swing width of the load, and a condition regarding the swing width speed of the load.
Control device.
前記制御対象の状態を静止状態から等速状態へと変化させる際、前記制御対象を示す状態空間モデルの状態変数と、第1指令値と、を決定変数とする評価関数を定式化するステップと、
所定の制約条件を有する前記評価関数を用いて、前記制御対象を前記静止状態から前記等速状態へと変化させる際の前記第1指令値を計算するステップと、
前記評価関数に基づいて、前記状態変数のうち、前記制御対象からの出力に対応する値を、目標値として計算するステップと、
前記制御対象からの出力と、前記目標値との差に応じた第2指令値を計算するステップと、
前記第1指令値と、前記第2指令値とを加算して前記指令値を計算するステップと、
を含み、
前記制約条件は、前記制御対象の状態空間モデルと、前記状態空間モデルの状態変数の条件を含む、
計算方法。 A method for calculating a command value for controlling a controlled object, comprising the steps of:
formulating an evaluation function having state variables of a state space model representing the control object and a first command value as decision variables when changing a state of the control object from a stationary state to a constant velocity state;
calculating the first command value when changing the controlled object from the stationary state to the constant velocity state by using the evaluation function having a predetermined constraint condition;
calculating, as a target value, a value of the state variables corresponding to an output from the controlled object based on the evaluation function;
calculating a second command value corresponding to a difference between an output from the controlled object and the target value;
calculating the command value by adding the first command value and the second command value;
Including,
the constraint condition includes a state space model of the controlled object and a condition of a state variable of the state space model;
Method of calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020125530A JP7508920B2 (en) | 2020-07-22 | 2020-07-22 | Control device and design method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020125530A JP7508920B2 (en) | 2020-07-22 | 2020-07-22 | Control device and design method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022021748A JP2022021748A (en) | 2022-02-03 |
JP7508920B2 true JP7508920B2 (en) | 2024-07-02 |
Family
ID=80220658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020125530A Active JP7508920B2 (en) | 2020-07-22 | 2020-07-22 | Control device and design method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7508920B2 (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002091570A (en) | 2000-09-20 | 2002-03-29 | Yaskawa Electric Corp | Servo control method |
JP2002354857A (en) | 2001-05-22 | 2002-12-06 | Yaskawa Electric Corp | Motor controller |
JP2004272883A (en) | 2003-02-20 | 2004-09-30 | Mitsubishi Electric Corp | Servo control device |
JP2006525928A (en) | 2003-05-30 | 2006-11-16 | リープヘル−ヴェルク ネンツィング ゲーエムベーハー | Crane or excavator with optimal movement guidance for handling rope-loaded luggage |
JP2009153250A (en) | 2007-12-19 | 2009-07-09 | Ricoh Co Ltd | Drive apparatus, image forming apparatus |
JP2013206460A (en) | 2012-03-29 | 2013-10-07 | Mitsubishi Electric Corp | Optimization of system control based on solving multipoint boundary value problem |
JP2020030797A (en) | 2018-08-20 | 2020-02-27 | 富士電機株式会社 | Controller designing device and designing method |
-
2020
- 2020-07-22 JP JP2020125530A patent/JP7508920B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002091570A (en) | 2000-09-20 | 2002-03-29 | Yaskawa Electric Corp | Servo control method |
JP2002354857A (en) | 2001-05-22 | 2002-12-06 | Yaskawa Electric Corp | Motor controller |
JP2004272883A (en) | 2003-02-20 | 2004-09-30 | Mitsubishi Electric Corp | Servo control device |
JP2006525928A (en) | 2003-05-30 | 2006-11-16 | リープヘル−ヴェルク ネンツィング ゲーエムベーハー | Crane or excavator with optimal movement guidance for handling rope-loaded luggage |
JP2009153250A (en) | 2007-12-19 | 2009-07-09 | Ricoh Co Ltd | Drive apparatus, image forming apparatus |
JP2013206460A (en) | 2012-03-29 | 2013-10-07 | Mitsubishi Electric Corp | Optimization of system control based on solving multipoint boundary value problem |
JP2020030797A (en) | 2018-08-20 | 2020-02-27 | 富士電機株式会社 | Controller designing device and designing method |
Also Published As
Publication number | Publication date |
---|---|
JP2022021748A (en) | 2022-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5762625B2 (en) | Trajectory control device | |
US7348748B2 (en) | Motorized system and method of control | |
US11640557B2 (en) | Machine learning device, numerical control system, and machine learning method | |
Boscariol et al. | Robust point-to-point trajectory planning for nonlinear underactuated systems: Theory and experimental assessment | |
US7818087B2 (en) | Method and control device for guiding the movement of a movable machine element of a machine | |
JP5366840B2 (en) | Trajectory control device | |
Zhang et al. | Time-optimal path tracking for robots under dynamics constraints based on convex optimization | |
WO2007046257A1 (en) | Motor control device, method, and program storage medium | |
Petko et al. | Trajectory tracking controller of the hybrid robot for milling | |
JP7508920B2 (en) | Control device and design method | |
JPH06274228A (en) | Numerical control device | |
CN110788859B (en) | Controller parameter universe self-adaptive adjustment system | |
JP2012137990A (en) | Numerical control unit, movement control method, movement control program, and storage medium | |
JP2023159641A (en) | Path generation method of generating path of robot using command from external apparatus and robot control system | |
Haas | Set point optimisation for machine tools | |
JPH0695294B2 (en) | Trajectory interpolation method for industrial robot | |
JP2003015715A (en) | Command value generating method and command value generation system and prgraming device | |
JP7502135B2 (en) | Robot program adjustment device | |
CN115128996B (en) | Error compensation method, device, equipment and medium for numerical control machining | |
Zhang et al. | An adaptive type-2 fuzzy sliding mode tracking controller for a robotic manipulator | |
JP7460750B2 (en) | numerical control device | |
JP6800384B1 (en) | Positioning control device and positioning method | |
JP2021144438A (en) | Controller system and control method thereof | |
KR101425057B1 (en) | Slight vibration continuous robot motion implementing method using real time motion profile generation algorithm for multi-articulated robot | |
JP2858696B2 (en) | Preview control method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230614 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240216 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240305 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240408 |
|
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: 20240521 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240603 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7508920 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |