(57)【要約】      （修正有） 【課題】コンピュータにより、熟練作業者による注湯作業に可及的に近づけることが可能な、取鍋の傾動による自動注湯の制御方法を提供する。【解決手段】注湯プロセスを遂行するプログラムで制御されるサーボモータにより取鍋を傾動させて鋳型に注湯するに当たり、所望の注湯流量パターンによって鋳型に注湯すべくサーボモータを制御する方法であって、サーボモータに印加する入力電圧から取鍋による注湯流量までの数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量とサーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定し、この注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理し、これにより、サーボモータに印加する入力電圧を獲得し、この獲得した入力電圧に基づきサーボモータを制御することを特徴とする。【選択図】図２
2008290148.tif 000002
【特許請求の範囲】
【請求項１】
注湯プロセスを遂行するプログラムを予め設定したコンピュータによって制御されるサーボモータにより取鍋を傾動させて鋳型に注湯するに当たり、前記取鍋による所望の注湯流量パターンによって鋳型に注湯すべく前記サーボモータを制御する方法であって、
前記サーボモータに印加する入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、
この作成した数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量と前記サーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定し、この注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理し、
これにより、前記サーボモータに印加する入力電圧を獲得し、
この獲得した入力電圧に基づき前記サーボモータを制御することを特徴とする自動注湯制御方法。
【請求項２】
請求項１に記載の自動注湯制御方法において、
前記注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理したのち、前記サーボモータに印加する入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、
この作成した数理モデルの逆問題を解き、かつ、目標注湯流量をＦＦ制御器で処理し、
ゲインスケジュールドＰＩ制御器で処理した結果とＦＦ制御器で処理した結果を足し合わせることにより、前記サーボモータに印加する入力電圧を獲得し、
この獲得した入力電圧に基づき前記サーボモータを制御することを特徴とする自動注湯制御方法。
【請求項３】
請求項１または２に記載の自動注湯制御方法において、
前記取鍋の傾動中に溶湯金属を含めた取鍋の重心が変動し、その重力方向の加速度が前記重量計測器の計測データに重畳される問題に対して、重心変動による加速度を打ち消し、重量計測器の計測データから重心変動の影響を取り除いて溶湯金属の鋳込み重量を計測することを特徴とする自動注湯制御方法。
【請求項４】
請求項１または２に記載の自動注湯制御方法において、
前記サーボモータへの入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、
この作成した数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量と前記サーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて、前記取鍋から流出する単位時間あたりの溶湯金属の重量を実時間で推定することを特徴とする自動注湯制御方法。
【請求項５】
請求項１〜４のうちいずれか１項に記載の自動注湯制御方法において、
前記取鍋は、円筒形状のもの、または扇形のものであることを特徴とする自動注湯制御方法。
【請求項６】
注湯プロセスを遂行するプログラムを予め設定したコンピュータによって制御されるサーボモータにより取鍋を傾動させて鋳型に注湯するに当たり、前記取鍋による所望の注湯流量パターンによって鋳型に注湯すべく前記サーボモータを制御するシステムであって、
前記サーボモータへの入力電圧から前記取鍋による注湯流量までの数理モデルを作成する数理モデル作成手段と、
作成した数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量と前記サーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定する注湯流量推定手段と、
この注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理するゲインスケジュールドＰＩ制御器処理手段と、
を含むことを特徴とする自動注湯装置のサーボモータの制御システム。
【請求項７】
注湯プロセスを遂行するプログラムを予め設定したコンピュータによって制御されるサーボモータにより取鍋を傾動させて鋳型に注湯するに当たり、前記取鍋による所望の注湯流量パターンによって鋳型に注湯すべく前記サーボモータを制御するための制御プログラムを記憶した記憶媒体であって、
前記サーボモータに印加する入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、
この作成した数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量と前記サーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定し、この注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理し、
これにより、前記サーボモータに印加する入力電圧を獲得し、
この獲得した入力電圧に基づき前記サーボモータを制御することを特徴とする取鍋用傾動制御プログラムを記憶した記憶媒体。
【請求項８】
請求項７に記載の取鍋用傾動制御プログラムを記憶した記憶媒体において、
注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理したのち、前記サーボモータに印加する入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、
この作成した数理モデルの逆問題を解き、かつ、目標注湯流量をＦＦ制御器で処理し、
ゲインスケジュールドＰＩ制御器で処理した結果とＦＦ制御器で処理した結果を足し合わせることにより、前記サーボモータに印加する入力電圧を獲得し、
この獲得した入力電圧に基づき前記サーボモータを制御することを特徴とする取鍋用傾動制御プログラムを記憶した記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自動注湯制御方法、自動注湯装置のサーボモータの制御システムおよび取鍋用傾動制御プログラムを記憶した記憶媒体に係り、より詳しくは、注湯プロセスを遂行するために予めプログラムを設定されたコンピュータによって制御されるサーボモータにより取鍋を傾動させて鋳型に溶湯を注入するに当たり、所望の注湯流量パターンによって溶湯を鋳型に注入すべく前記サーボモータを制御する方法、自動注湯装置のサーボモータの制御システムおよび取鍋用傾動制御プログラムを記憶した記憶媒体に関する。
【背景技術】
【０００２】
鋳造工場における注湯のように極めて危険でかつ最悪の作業から労働者を解放すべく、注湯プロセスの機械化・自動化が、近年行われるようになってきている。そして、従来、このための装置としては、取鍋と、取鍋を駆動する駆動手段と、取鍋の重量を検出する検出手段と、予め取鍋が傾動されたときの取鍋内の重量の変動割合を記憶しておき、前記検出手段からの信号に対応して取鍋の傾動速度を補正し、前記駆動手段に補正後の傾動速度信号を送信する記憶演算装置とを具備したものがある（例えば、特許文献１参照）。
【特許文献１】特開平６-７９１９号
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、このように構成された従来の自動注湯装置においては、駆動手段等に係る情報の記憶演算装置への入力が、現実的にはティーチング＆プレイバック方式により行われているため、不適切な取鍋傾動速度や注湯状況の変化に対応できず、この結果、鋳型に注入される溶湯が流量不足になったり、注湯時にほこり・のろなどの不純物が鋳型内に飲み込まれたりして、鋳物の品質低下を招くなどの問題があった。
【０００４】
本発明は上記の事情に鑑みてなされたもので、その目的は、予めプログラムを設定されたコンピュータにより、熟練作業者による注湯作業に可及的に近づけることが可能な、取鍋の傾動による自動注湯の制御方法、自動注湯装置のサーボモータの制御システムおよび取鍋用傾動制御プログラムを記憶した記憶媒体を提供することにある。
【課題を解決するための手段】
【０００５】
上記の目的を達成するため、本発明における自動注湯制御方法は、注湯プロセスを遂行するプログラムを予め設定したコンピュータによって制御されるサーボモータにより取鍋を傾動させて鋳型に注湯するに当たり、前記取鍋による所望の注湯流量パターンによって鋳型に注湯すべく前記サーボモータを制御する方法であって、前記サーボモータに印加する入力電圧から前記取鍋による注湯流量までの数理モデルを作成し、この作成した数理モデルの逆問題を解き、かつ、重量計測器によって計測される溶湯金属の鋳込み重量から重心変動の影響を取り除いた重量と前記サーボモータに印加する入力電圧から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定し、この注湯流量と目標注湯流量をゲインスケジュールドＰＩ制御器で処理し、これにより、前記サーボモータに印加する入力電圧を獲得し、この獲得した入力電圧に基づき前記サーボモータを制御することを特徴とする。
【０００６】
なお、本発明に利用する数理モデル法とは、プロセスの熱収支・物質収支・化学反応・制限条件などの式を解いて、利益・コストなどコンピュータ制御の目的とする関数をだし、その最大・最小を求めてそれが達成できるように制御を行う方法である。
またなお、取鍋は重心付近で支持してある。
【発明の効果】
【０００７】
上記の説明から明らかなように本発明によれば、プログラムを予め設定されたコンピュータにより、熟練作業者による注湯作業に可及的に近づけた状態で取鍋によって自動注湯を行うことが可能になるうえに、注湯流量を推定してフィードバックするため，目標注湯流量が変動する場合や，注湯工程で外乱要素が存在する場合においても所望の注湯流量を精度良く実現できるなどの優れた実用的効果を奏する。
【発明の実施の形態】
【０００８】
以下、本発明を適用した自動注湯装置の実施例について図１〜図１４に基づき詳細に説明する。図１に示すように、本自動注湯装置は、円筒形状の取鍋１と、この取鍋１を傾動させるサーボモータ２と、サーボモータの出力軸の回転運動を直線運動に変換する２組のボールねじ機構３・４により前記取鍋１および前記サーボモータ２を垂直方向および水平方向へそれぞれ移動させる移動手段５と、前記取鍋１内の溶湯の重量を検出するロードセル（図示せず）と、コンピュータを利用して前記サーボモータ２および前記２組のボールねじ機構３・４の動作を演算しかつ制御するコントロールシステム６と、で構成してある。
【０００９】
そして、前記取鍋１は、これの重心位置に前記サーボモータ２の出力軸を連結させてその重心位置で傾動可能に支持してあって、重心位置を中心にして鋳型の湯口に対して傾動・反傾動するようになっている。
なお、重心位置を中心にして傾動するようにすることにより、前記サーボモータ２にかかる負荷が大きくなることを防ぐことができる。
【００１０】
また、前記移動手段５は、鋳型の湯口に正確に注湯すべく前記取鍋１を傾動に連動させて前後移動および昇降させ、その出湯口先端を仮想回転軸として固定出湯点を得ることができるよう作動する。
【００１１】
このように構成したものは、サーボモータ２に印加する入力電圧による取鍋1の傾動と、取鍋1の傾動によって取鍋1から流出する溶湯の流量に関する数理モデルを作成し、この作成した数理モデルの逆問題を解くことにより前記サーボモータ２に印加する入力電圧を獲得し、かつ、重量計測器としてのロードセルによって計測される鋳込み重量から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて注湯流量を推定し、この推定した注湯流量をゲインスケジュールドＰＩ制御器で処理して前記サーボモータ２に印加する入力電圧を獲得し、この獲得した入力電圧に基づきコントロールシステム６を介して取鍋1の傾動を制御する。
【００１２】
すなわち、取鍋1の注湯時の縦断面図である図２において、取鍋1の傾動角度をθ［deg］、取鍋1の傾動中心である出湯口より下部の溶湯体積（濃い網掛け部）をＶs（θ）［ｍ3］、出湯口に対する水平面の面積（濃い網掛け部と薄い網掛け部の境界上の面積）をＡ（θ）［ｍ２］、出湯口より上部の溶湯体積（薄い網掛け部）をＶｒ［ｍ3］、上部溶湯の高さをｈ［ｍ］、取鍋1から流出する溶湯の流量をｑ［ｍ3/ｓ］とすると、注湯時における時刻ｔ［ｓ］からΔｔ［ｓ］後の取鍋内溶湯の収支式は下記の式（1）のようになる。
Ｖｒ（ｔ）＋Ｖs（θ（ｔ））
＝Ｖｒ（ｔ＋Δｔ）＋Ｖs（θ（ｔ＋Δｔ））＋ｑ（ｔ）Δｔ  （1）
【００１３】
式（1）を溶湯体積Ｖｒ［ｍ3］についてまとめ、Δｔ→０とすると下記の式（２）となる。
【００１４】
【数１】
2008290148.tif 000003
【００１５】
また、取鍋1の傾動角速度ω［deg/ｓ］を下記の式（３）とする。
ω（ｔ）＝ｄθ（ｔ）/ｄｔ      （３）    
よって、式（３）を式（２）に代入すると、下記の式（４）が得られる。  
【００１６】
【数２】
2008290148.tif 000004
【００１７】
また、出湯口より上部の溶湯体積Ｖｒ［ｍ3］は下記の式（５）で表すことができる。
【００１８】
【数３】
2008290148.tif 000005
ここで、面積Ａｓ［ｍ２］は、図３に示す出湯口水平面からの高さｈｓ［ｍ］における溶湯水平面積を示す。
【００１９】
また、面積Ａｓ［ｍ２］を出湯口水平面の面積Ａ［ｍ２］と面積Ａ［ｍ２］に対する面積変化量ΔＡｓ［ｍ２］に分割すると、溶湯体積Ｖｒ［ｍ3］は下記の式（６）となる。
【００２０】
【数４】
2008290148.tif 000006
【００２１】
また、取鍋1を含む一般的な取鍋においては、面積変化量ΔＡｓ［ｍ２］は出湯口水平面の面積Ａ［ｍ２］に対して微小であるから、下記の式（７）が得られる。
【００２２】
【数５】
2008290148.tif 000007
したがって、式（６）は下記の式（８）と示すことができる。
Ｖｒ（ｔ）≒Ａ（θ（ｔ））ｈ（ｔ）  （８）
よって、式（８）より下記の式（９）が得られる。    
ｈ（ｔ）≒Ｖｒ（ｔ）/Ａ（θ（ｔ））  （９）  
【００２３】
また、ベルヌーイの定理を用いて、出湯口より上部の溶湯高さｈ［ｍ］から溶湯流量ｑ［ｍ3/ｓ］までを下記の式（10）で示す。
【００２４】
【数６】
2008290148.tif 000008
ここで、ｈｂ［ｍ］は図４に示すように取鍋1内の溶湯の上面からの溶湯深さ、Ｌｆ［ｍ］は溶湯深さｈｂ［ｍ］における出湯口の幅、ｃは流量係数、ｇは重力加速度をそれぞれ示す。
【００２５】
また、式（４）、式（９）および式（１０）より注湯流量モデルの基礎式は下記の式（１１）および式（１２）となる。
【００２６】
【数７】
2008290148.tif 000009
【００２７】
【数８】
2008290148.tif 000010
【００２８】
また、出湯口に対する水平面の面積Ａ（θ）［ｍ２］は取鍋1の傾動角度θ［deg］に対して変動する。したがって、式（１１）および式（１２）の注湯流量モデルは、システム行列、入力行列および出力行列が取鍋1の傾動角度に依存して変動する非線形パラメータ変動モデルとなる。
【００２９】
図５は本自動注湯装置の第１実施例における注湯プロセスのブロック線図を示し、図５において、モータモデルは下記の式（１６）の１次遅れ系で示される。
ｄω（ｔ）/ｄｔ＝−ω（ｔ）/Ｔｍ＋Ｋｍｕ（ｔ）/Ｔｍ      （１６）  
ここで、Ｔｍ［ｓ］は時定数、Ｋｍ［deg/ｓＶ］はゲイン定数をそれぞれ示す。本自動注湯装置ではＴｍ＝０．００６［ｓ］、Ｋｍ＝２４．５８［deg/ｓＶ］である。
【００３０】
また、ロードセルの動特性を考慮してロードセルＰＬを下記の式（１７）で示す。
ｄｗＬ/ｄｔ＝−ｗＬ（ｔ）/ＴＬ＋ｗ（ｔ）/ＴＬ      （１７）
ここで、ｗ［Ｋｇ］は取鍋1から流出した液体の流出重量、ｗＬ［Ｋｇ］はロードセルで計測される計測重量、ＴＬ［ｓ］はロードセルの応答遅れを示す時定数である。本自動注湯装置ではステップ応答法により時定数を同定した結果、ＴＬ＝０．１０［ｓ］となった。
【００３１】
式（１１）および式（１２）に示す注湯流量モデルにおいて、図６は取鍋1の各傾動角度θ［deg］に対する出湯口面積水平Ａ（θ）［ｍ２］と出湯口下部の溶湯(液体)体積Ｖｓ（θ）［ｍ3］を示す。図６において、（ａ）は取鍋1の傾動角度θ［deg］に対する出湯口面積水平Ａ（θ）［ｍ２］、（b）は取鍋1の傾動角度θ［deg］に対する出湯口下部の溶湯(液体)体積Ｖｓ（θ）［ｍ3］を示す。
【００３２】
次に、上述のようにして求めた注湯流量モデルを用いて、逆モデルによる注湯流量フィードフォワード制御を構築する。
なお、フィードフォワード制御とは、制御対象に加える操作量を予め決められた値に調節することにより、出力が目標値になるようにする制御法であって、制御対象の入出力関係や外乱の影響などが明確な場合には性能の良い制御を行うことができる。
【００３３】
図７は、所望の注湯流量パターンｑｒｅｆ［ｍ3/ｓ］を実現するためサーボモータ２へ印加する制御入力ｕ［Ｖ］を導出するシステムにおける制御系のブロック線図を示す。ここで、サーボモータ２の逆モデルＰｍ-1は下記の式（１８）により示される。
【００３４】
【数９】
2008290148.tif 000011
【００３５】
式（１１）および式（１２）に示す注湯流量モデルの基礎式に対する逆モデルを導出する。ベルヌーイの定理である式（１０）より出湯口上部の溶湯高さｈ［ｍ］に対する注湯流量ｑ［ｍ3/ｓ］を求めることができる。取鍋1の形状から考えられる出湯口上部の最大溶湯高さｈｍａｘ［ｍ］をｎ分割したときの分割幅をΔｈ［ｍ］とし、各々の溶湯高さをｈｉ＝ｉΔｈ（ｉ＝0、…ｎ）で示す。したがって、溶湯高さｈ＝［ｈ0ｈ1…ｈｎ］Ｔに対する注湯流量ｑ＝［ｑ0ｑ1…ｑｎ］Ｔを下記の式（１９）に示す。
ｑ＝ｆ（ｈ）  （１９）
ここで、関数ｆ（ｈ）は式（１０）に示すベルヌーイの定理である。したがって、式（１９）の逆関数は下記の式（２０）となる。
ｈ＝ｆ−１（ｑ）      （２０）
【００３６】
この式（２０）は式（１９）をＬｏｏｋｕｐ  Ｔａｂｌｅで表現し、入出力関係を逆にすることで表すことができる。
ここで、分割間隔ｑｉ→ｑｉ+1、ｈｉ→ｈｉ+1は線形補間により近似する。分割幅が小さいほど、高精度に注湯流量ｑ［ｍ3/ｓ］と出湯口上部の溶湯高さｈ［ｍ］の関係を表現できる。実装可能な範囲で分割幅を小さくすることが望まれる。
【００３７】
所望の注湯流量パターンｑｒｅｆ［ｍ3/ｓ］を実現する出湯口上部の溶湯高さｈｒｅｆ［ｍ］は式（２０）より下記の式（２１）となる。
ｈｒｅｆ（ｔ）＝ｆ−１（ｑｒｅｆ（ｔ））      （２１）
【００３８】
また、出湯口上部の溶湯高さｈｒｅｆ［ｍ］における出湯口上部の溶湯体積Ｖｒｅｆ［ｍ3］は、式（９）を用い下記の式（２２）で示す。
Ｖｒｅｆ（ｔ）＝Ａ（（θ（ｔ））ｈｒｅｆ（ｔ）      （２２）
【００３９】
次に、式（２２）で得られた出湯口上部の溶湯体積Ｖｒｅｆ［ｍ3］と所望の注湯流量パターンｑｒｅｆ［ｍ3/ｓ］を、式（１１）の注湯流量モデルの基礎式に代入して、下記の式（２３）に示す所望の注湯流量パターンを実現する取鍋1の傾動角速度ωｒｅｆ［deg/ｓ］を導出する。
【００４０】
【数１０】
2008290148.tif 000012
【００４１】
まず、式（１９）から式（２３）を順に解き、得られた取鍋1の傾動角速度ωｒｅｆ［deg/ｓ］を式（１８）に代入することにより、所望の注湯流量パターンｑｒｅｆ［ｍ3/ｓ］を実現すべくサーボモータ２へ印加する制御入力ｕ［Ｖ］を得ることができる。
【００４２】
式（２２）より得られた出湯口上部の溶湯体積Ｖｒｅｆ［ｍ3］と所望の注湯流量パターンｑｒｅｆ［ｍ3/ｓ］を式（２３）に代入すると、所望の注湯流量パターンを実現する取鍋1の傾動角速度ωｒｅｆ［deg/ｓ］が得られる。そして、得られた取鍋1の傾動角速度ωｒｅｆ［deg/ｓ］を、式（１８）のサーボモータ２の逆モデルに代入すると、サーボモータ２へ印加する制御入力ｕ［Ｖ］を得ることができる。
【００４３】
式（１１）、式（１２）および式（１７）に示す取鍋傾動式注湯装置の注湯プロセスに対して、図５は、注湯流量数理モデルの逆モデルによる注湯流量フィードフォワード制御と、ゲインスケジュールドＰＩ制御による注湯流量フィードバック制御を併合した２自由度制御システムを示す。
そして，フィードバック部は，ロードセルによって計測される鋳込み重量から拡張カルマンフィルタに基づく指数減衰型オブザーバを用いて，注湯流量を推定し，この推定した注湯流量をゲインスケジュールドＰＩ制御器で処理することで、外乱要素が存在する環境下でも所望の注湯流量パターンに追従させる注湯流量制御システムを構築する。 
【００４４】
フィードフォワード部は目標値追従を良好にし、フィードバック部は定常誤差および外乱の影響を除去する役割をもつ。また、（１１）、（１２）式の注湯流量モデルは、注湯流量に対して非線形特性を持つため，非線形特性にも対応できるよう，フィードバックコントローラに，注湯流量に応じて、比例ゲインと積分ゲインが変動するゲインスケジュールドＰＩ制御系を用いている。
【００４５】
上述の注湯流量２自由度制御システムを、対象液体を水とした取鍋傾動式自動注水装置に適用した際の実験結果を図８に示す。ここで，本実験における外乱要素は，事前に取鍋内液体容量から求めた出湯開始角度より，＋2[deg]傾動した際に取鍋内の液体が流出を開始する角度誤差とした。
【００４６】
図８において，破線は目標注水流量パターンであり，実線は，本発明による注湯流量２自由度制御システムにて実験を行った注水流量の結果であり，一点鎖線は，注湯流量フィードフォワード制御の注水流量の結果である。この結果より，注湯流量２自由度制御システムでは，目標注水流量パターンが変動した際にも実際の注水流量が追従でき，また，外乱要素が存在した際にも，精度良く追従できていることが確認された。
【００４７】
次に、本自動注湯装置の第２実施例として、前記取鍋１の傾動時の重心変動に対するロードセル補償法を適用した傾動式自動注湯装置について述べる。上述の第１実施例の傾動式自動注湯装置は，溶湯落下位置を安定化させるために，取鍋１を出湯口中心に回転するように取鍋１を傾動に合わせて上下，前後に動作している。この取鍋１の上下動による取鍋１の重心変動により，上下方向に力が生じ，ロードセルによって計測される鋳込み重量へ影響を与える。これにより，真の鋳込み重量を得ることができなくなる。そこで、本第２実施例の傾動式自動注湯装置では、ロードセルによる計測注湯鋳込み重量を用いて注湯流量推定を行なうが，取鍋１の重心変動によって推定精度が低下する。そこで，高精度推定を可能とするために，取鍋１の重心変動によるロードセルへの影響を除去するロードセル補償法を構築する。図９にロードセル補償法のブロック線図を示す。ここで，GMvは取鍋重力方向動作のモータモデル，GLvは重力方向加速度からロードセルへの影響までを表したロードセルモデルである。
【００４８】
ロードセル補償法で用いるロードセルモデルは式(２７)に示すような二次遅れ系，垂直方向のモータモデルは式（２６）に示す一次遅れ系のモデルを用いる。Ｋmz［ｍＶ／ｓ］はモータゲイン，Ｔmzs［ｓ］はモータ時定数，Ｋl［ｋｇｓ２／ｍ］はロードセルゲイン，ωnl［ｒａｄ／ｓ］はロードセル固有振動数，ζlはロードセル減衰係数である。ここで，同定実験の結果Ｋmz＝０．０８２８［ｍ／ｓＶ］、Ｔmz＝０．００７［ｓ］となった。
GMv（ｓ）＝Kmz／（1＋Tmzs）                      （２６）
GLv（ｓ）＝Klωnls／（s2＋2ζlωnl＋ω2nl）        （２７）
【００４９】
また，ロードセルモデルのパラメータは，Ｋl=０．１８４，ωnl=０．７５０，ζｌ＝７．４４となった。ロードセル補償法を用いてロードセルにより計測された計測鋳込み重量から，取鍋１の上下動加速度により励起された力の除去したときの結果を図１０に示す。実験結果より，シミュレーションの鋳込み重量と補償結果が一致している。したがって，このロードセル補償法を流量推定システムに適用することにより，より高性能な注湯流量推定が可能となる。
【００５０】
次に、注湯流量推定方法について述べる。
拡張カルマンフィルタに基づいた指数減衰型オブザーバについて説明すると、離散時間系の拡張カルマンフィルタ［文献：Ｋ．Ｒｅｉｆ，Ｒ．Ｕｎｂｅｈａｕｅｎ，Ｔｈｅ  Ｅｘｔｅｎｄｅｄ  Ｋａｌｍａｎ  Ｆｉｌｔｅｒａｓ  ａｎ  Ｅｘｐｏｎｅｎｔｉａｌ  Ｏｂｓｅｒｂｅｒ  ｆｏｒ  Ｎｏｎｌｉｎｅａｒ  Ｓｙｓｔｅｍｓ，ＩＥＥＥ，Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，Ｖｏｌ．４７，Ｎｏ．８，（１９９９），ｐｐ．２３２４−２３２８．］を用いて，指数減衰型オブザーバの構築を行なう。以下にアルゴリズムを示す。対象とするシステムを式(２８)，(２９)に示す。
zn+1＝f（ｚｎ，ｘｎ）        （２８）
yn＝h（ｚｎ）              （２９）
ここで，n∈N0は離散時間，ｚｎ∈Rq，ｘｎ∈Rq，ｙｎ∈Rmはそれぞれ，状態変数，入力，出力である。また，関数ｆ，ｈはC1関数であると仮定する。式(２８),(２９)により，オブザーバは式(３０)，(３１)となる。ここで，オブザーバゲインKnは時変のｑ×m行列である。
【００５１】
【数１１】
2008290148.tif 000013
【００５２】
また，推定された状態量
【００５３】
【数１２】
2008290148.tif 000014
はそれぞれ，ｐｒｉｏｒｉ  ｅｓｔｉｍａｔｅ，ｐｏｓｔｅｒｉｏｒｉ  ｅｓｔｉｍａｔｅ  と呼ばれる。
【００５４】
このとき，オブザーバゲインKnを拡張カルマンフィルタのカルマンゲイン更新アルゴリズムを用いて更新する。拡張カルマンフィルタのカルマンゲインの更新アルゴリズムを式(３２)〜(３８)に示す。ここで，Qはｑ×ｑ正定対称行列，Ｒはｍ×ｍ正定対称行列，α≧１の実数である。
【００５５】
【数１３】
2008290148.tif 000015
【００５６】
Q，Rはそれぞれ拡張カルマンフィルタでは，システムノイズ，観測ノイズの共分散行列を意味している。また，αは収束の度合いを調節するパラメータである。α＝1のとき，拡張カルマンフィルタとなる。
以上の拡張カルマンフィルタの更新アルゴリズムを用いたオブザーバは，指数減衰型オブザーバとなることが文献［Ｋ．Ｒｅｉｆ，Ｒ．Ｕｎｂｅｈａｕｅｎ，Ｔｈｅ  Ｅｘｔｅｎｄｅｄ  Ｋａｌｍａｎ  Ｆｉｌｔｅｒａｓ  ａｎ  Ｅｘｐｏｎｅｎｔｉａｌ  Ｏｂｓｅｒｂｅｒ  ｆｏｒ  Ｎｏｎｌｉｎｅａｒ  Ｓｙｓｔｅｍｓ，ＩＥＥＥ，Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，Ｖｏｌ．４７，Ｎｏ．８，（１９９９），ｐｐ．２３２４−２３２８．］により証明されている。
【００５７】
次に、この離散時間型拡張カルマンフィルタを用い，注湯流量推定システムの構築を行う。
前記取鍋１の傾動角速度からロードセルによる計測注湯鋳込み重量までのシステムを考える。まず，連続時間系の注湯システムである式(１１)、(１２)、(１６)、（１７）の微分方程式を差分方程式へと変換する。求めた差分方程式を式（３９），（４０）に示す。
ここで，ｔ＝ｎｋｓ，ｔｓ［ｓ］はサンプリング時間，ｎはサンプリング番号ｎ＝１，２，３，・  ・  ・とする。
【００５８】
【数１４】
2008290148.tif 000016
【００５９】
ｙ（ｎ）＝ｗｌ（ｎ）        （４０）
ａｆ、ｂｆは式（４１）、式（４２）で示される。
【００６０】
【数１５】
2008290148.tif 000017
【００６１】
ｂｆ（ｎ）＝−∂Ｖｓ（θ（ｎ））／∂θ（ｎ）    （４２）
である。
【００６２】
式（３９）、式（４０）のシステムモデルに対して、指数減衰型オブザーバーを構築する。
式（３９）、式（４０）を式（３０）、式（３１）で表現すると式（４３）〜式（４６）となる。
ｚｎ＝［Ｖｒ（ｎ）  ｗ（ｎ）  ｗｌ（ｎ）］Ｔ      （４３）
ｘ（ｎ）＝ω（ｎ）                              （４４）
【００６３】
【数１６】
2008290148.tif 000018
【００６４】
ｈ（ｚｎ）＝ｗｌ（ｎ）        （４６）
【００６５】
【数１７】
2008290148.tif 000019
【００６６】
式（４７）を用いてシミュレーションより求めたwlと実際の実験結果を比較し，vの分散を求める。傾動開始角度と流出開始角度が3[deg]ずれていたときまでの推定が可能であることを目標とする。そこで，3[deg]の初期角度のずれをシステムノイズとすることで，初期角度のずれを考慮した注湯流量推定システムを設計する。図１０に3［deg]ずれた時の実験結果とシステムノイズを考慮した式（４７）より得られるwlを示す。各システムノイズの分散は，3[deg]ずれた時の流出重量の結果に近づくように，Σvq＝１．０×１０−１０［ｍ６／ｓ２］，Σvw＝１．０×１０−１２［m6］，Σｖｗｌ＝１．０×１０-12［ｍ６］とした。図１１に示すように，システムノイズを付加することにより，シミュレーションでの鋳込み重量が，３［ｄｅｇ］の初期角度ずれが生じた時の実験結果に近づいていることが確認できる。
【００６７】
これより，システムノイズの共分散行列Q を式(４８)とする。
【００６８】
【数１８】
2008290148.tif 000020
【００６９】
構築した離散時間系拡張カルマンフィルタを用いた指数減衰型拡張オブザーバにより注湯流量を推定する。推定シミュレーション，および推定実験結果を図１２に示す。このときのオブザーバゲインを図１３に示す。ここで，オブザーバゲインをKn＝［Ｋｑ  Ｋｗ  Ｋｗｌ］Tとする。推定シミュレーションおよび実験結果より，精度の高い推定が可能であることが確認できる。
【００７０】
また，鋳造工場において，取鍋１内に溶湯を供給する場合，手作業により行なわれている。そのため，所望の容量だけ正確に注ぐことは非常に困難である。よって，取鍋１から溶湯が出始めるときの取鍋１の角度にばらつきが生じてしまう。ここで，取鍋１の内重量と取鍋１の形状が既知であれば，この溶湯が出始めるときの角度は計算により求めることができるが，取鍋１の内部は手作業で成形されており，正確な形状を得ることができず，溶湯流出開始角度は得ることは困難である。そこで，所望の注湯流出開始角度と実際の注湯流出開始角度がずれてしまう。したがって，取鍋傾動開始角度と注湯流出開始角度にずれが生じた場合の注湯流量推定実験を行う。
【００７１】
それぞれ初期角度26［ｄｅｇ］から1，3，5［ｄｅｇ］のずれが生じた場合の注湯流量推定結果を図１４に示す。図１４に示すように，3［ｄｅｇ］以上のずれが生じた場合，初期の推定誤差は大きくなるが，それ以降の注湯流量推定は可能であることが確認され，高精度に注湯流量を推定可能であることが確認される。実際の鋳造工場では，計算された取鍋の出湯開始角度と実際の出湯開始角度との誤差は2［ｄｅｇ］程度であることから，実際の運用においても注湯流量の推定精度は高いことが確認された。拡張カルマンフィルタを用いたオブザーバはシステムノイズ，観測ノイズを設定するだけで，オブザーバゲインをシステマチックに得ることが可能である。また，システムノイズの共分散行列を操作することにより，ある程度外乱が生じた場合でも所望の状態量を推定することが可能である。
【００７２】
なお、本自動注湯装置の第２実施例においては、取鍋１の重心変動によるロードセルへの影響を除去するロードセル補償法を構築するとしてあるが、このロードセルは、溶湯および取鍋の静止重量と、これらの上下動による加速度によって励起される力と、を同時に測定可能ならどの場所に設けても良く、例えば、取鍋を搭載して上下方向および水平方向へ移動させる移動部材に設けてもよい。
【図面の簡単な説明】
【００７３】
【図１】本発明を適用した自動注湯装置の一実施例を示す模式図である。
【図２】図１の自動注湯装置における取鍋の縦断面図である。
【図３】図２における要部拡大詳細図である。
【図４】取鍋の注湯口先端の斜視図である。
【図５】本自動注湯装置の第１実施例における注湯プロセスのブロック線図である。
【図６】取鍋1の各傾動角度θ［deg］に対する出湯口面積水平Ａ（θ）［ｍ２］と出湯口下部の溶湯(液体)体積Ｖｓ（θ）［ｍ3］を示すグラフである。
【図７】注湯流量フィードフォワード制御系のブロック図である。
【図８】注湯流量２自由度制御システムを対象液体を水とした取鍋傾動式自動注水装置へ適用した際の実験結果を示すグラフである。
【図９】ロードセル補償法のブロック線図である。
【図１０】ロードセル補償法を用いた計量鋳込み重量を示すグラフである。
【図１１】システムノイズを伴う注湯シミュレーションを示すグラフである。
【図１２】離散時間系拡張カルマンフィルタを用いた指数減衰型拡張オブザーバにより流量を推定した推定シミュレーション，および推定実験結果を示すグラフである。
【図１３】図１２におけるオブザーバゲインを示すグラフである。
【図１４】注湯出湯開始初期角度の誤差が生じた場合の注湯流量推定結果を示すグラフである。
【図１】
2008290148.tif 000021
【図２】
2008290148.tif 000022
【図３】
2008290148.tif 000023
【図４】
2008290148.tif 000024
【図５】
2008290148.tif 000025
【図６】
2008290148.tif 000026
【図７】
2008290148.tif 000027
【図８】
2008290148.tif 000028
【図９】
2008290148.tif 000029
【図１０】
2008290148.tif 000030
【図１１】
2008290148.tif 000031
【図１２】
2008290148.tif 000032
【図１３】
2008290148.tif 000033
【図１４】
2008290148.tif 000034