(57)【要約】  電動射出成形機の圧力制御に関し、高価で、煩雑な取付け作業、複雑な機械構造、更にはノイズ対策やゼロ点調整等の人手作業を必要とする圧力検出センサを使わずに、高精度で時間的に遅れの非常に小さい射出圧力検知方法が望まれていた。  本発明は、電動射出成形機の射出・保圧機構の数式モデルを表現した状態方程式を利用して導出した離散演算式を内蔵し、入力信号としてスクリュ位置信号とモータ電流指令或いはモータ実電流を受取り、上記状態方程式表現の状態変数の１つである射出圧力の推定値を出力する高ゲイン観測器を採用する。高ゲイン観測器は、圧力検出センサを必要とせず、高精度で、時間的遅れの非常に小さい射出圧力推定値を出力する。したがって、高ゲイン観測器の出力する射出圧力推定値は、射出圧力制御での射出圧力検出値のフィードバック信号として採用できる。
2010134359.tif 000002
【特許請求の範囲】
【請求項１】
電動射出成形機の圧力制御装置であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルを利用して導出した離散演算式を実行する高ゲイン観測器と、射出圧力を指令するための射出圧力指令器と、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或いはモータ実電流信号が入力される前記高ゲイン観測器が、内蔵する前記離散演算式を使用して算出して出力する射出圧力推定値と前記射出圧力指令器が出力する射出圧力指令とが入力されて、前記射出圧力指令と前記射出圧力推定値との差を出力する減算器と、前記減算器の出力を入力して、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出する圧力制御器と、を含む射出制御器と、
前記モータ電流指令信号が入力されるモータ制御器と、
を具備することを特徴とする電動射出成形機の圧力制御装置
【請求項２】
前記高ゲイン観測器が、前記スクリュ位置信号と前記モータ電流指令信号或いはモータ実電流信号を入力信号として下記(数１０４)の状態方程式で表現される前記射出軸系の数式モデルを利用して導出した下記(数１０５)の離散式(１２５)および(数１０６)の離散式(１２８)からなる前記離散演算式を実行して、射出圧力推定値ｘ^３および射出速度推定値ｘ^２を出力し、前記射出圧力推定値を前記減算器に入力し、前記減算器は前記射出圧力指令との差を算出して前記圧力制御器に入力し、前記圧力制御器は、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出して前記モータ制御器に出力することにより、サーボモータに所定のモータトルクを発生させて射出圧力を制御する電動射出成形機の圧力制御方法
【数１０４】
2010134359.tif 000106
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数
ｘ２：射出速度を最大射出速度で無次元化した状態変数
ｘ３：射出圧力を最大射出圧力で無次元化した状態変数
ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数
ｙ：計測できる状態変数ｘ１を表す出力変数
ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ：射出・保圧機構のモデル定数
ｆ(ｘ３)：射出圧力で決まる金型への樹脂流入量を決める関数
χ(ｘ２)、ψ(ｘ)：式(１２４)で表される非線形関数
【数１０５】
2010134359.tif 000107
ここで、η^１、η^２：状態変数ｘ２、ｘ３を推定するために導入した新状態変数
Ｉ２：２行２列の単位行列
ｋ：時間ｔｋを表す離散変数（ｋ＝０、１、２、・・・）
Δｔ：高ゲイン観測器の演算周期
ε：高ゲイン観測器でのパラメータで、一般に1より十分小さい正数
K１、Ｋ２：高ゲイン観測器のゲインを決めるパラメータで、行列A０の固有値の実数部が負になるように決める。
【数１０６】
2010134359.tif 000108
ここで、ｘ^２、ｘ^３：状態変数ｘ２、ｘ３の推定値
【請求項３】
前記高ゲイン観測器が、前記スクリュ位置信号と射出速度信号（モータ回転速度）と前記モータ電流指令信号或いはモータ実電流信号を入力信号として下記(数１０７)の状態方程式で表現される射出軸系の数式モデルを利用して導出した下記(数１０８)の離散式(１３２)および(数１０９)の離散式(１３５)からなる前記離散演算式を実行して、射出圧力推定値ｘ^３、射出速度推定値ｘ^２およびスクリュ位置推定値ｘ^１を出力し、前記射出圧力推定値を前記減算器に入力し、前記減算器は前記射出圧力指令との差を算出して前記圧力制御器に入力し、前記圧力制御器は、前記射出圧力推定値が前記射出圧力指令に追従するようにモータ電流指令信号を算出して前記モータ制御器に出力することにより、サーボモータに所定のモータトルクを発生させて射出圧力を制御する電動射出成形機の圧力制御方法
【数１０７】
2010134359.tif 000109
ここで、ｙ１、ｙ２：計測できる状態変数ｘ１及びｘ２をそれぞれ表す出力変数
その他の変数の意味は、(数１０４)と同じである。
【数１０８】
2010134359.tif 000110
ここで、η^１、η^２、η^３：状態変数ｘ１、ｘ２、ｘ３を推定するために導入した新状態変数
Ｉ３：３行３列の単位行列
ｋ：時間ｔｋを表す離散変数(ｋ＝０、１、２、・・・)
Δｔ：高ゲイン観測器の演算周期
ε：高ゲイン観測器でのパラメータで、一般に1より十分小さい正数
K０：高ゲイン観測器のゲインを決める行列で、行列A０の固有値の実数部が負になるように行列要素を決める。
【数１０９】
2010134359.tif 000111
ここで、ｘ^1、ｘ^２、ｘ^３：状態変数ｘ１、ｘ２、ｘ３の推定値
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電動射出成形機の圧力制御装置および圧力制御方法に関するものである。
【背景技術】
【０００２】
高精度な制御と優れた応答性を特長とするＡＣサーボモータは、モータに用いる永久磁石の性能向上とコストダウンの実現により大容量化が図られ、その結果、従来は油圧駆動であった中型射出成形機（型締力３５０ｔ以上）でもＡＣサーボモータを用いたＡＣサーボ駆動が適用されるようなった。【０００３】
射出成形機には、スクリュ回転により樹脂ペレットを溶融する可塑化機構、スクリュの前進移動により溶融樹脂を金型に高速射出し、高圧で保持する射出・保圧機構、及び金型を開閉する型開閉機構がある。ＡＣサーボ駆動はこれらのどの機構にも適用できるが、ここでは、射出・保圧機構でのＡＣサーボ技術を説明するために、射出・保圧機構の模式図を図３に示す。【０００４】
射出成形機には、地上に固定された図示しない射出成形機基盤上に、リニアガイドにより移動可能な図示しない移動装置台が設けられる。この移動装置台に図３に示す金型１以外の全部品が搭載される。この移動装置台を移動させることにより、射出シリンダ２の先端部を金型１に押し付けることができ、また、逆に金型１から射出シリンダ２の先端部を離すこともできる。図３は、金型１内に樹脂を射出するために射出シリンダ２の先端部を金型１に押し付けている状態を示す。【０００５】
移動装置台の上には、サーボモータ３、減速機４、ボールネジ５及び軸受６が固設されている。また、上記ボールネジのナット７、可動部８、スクリュ９及び圧力検出センサ１０が一体構造となっている。この一体構造は、ボールネジ５のナット７の移動により前後移動可能なように、可動部８がリニアガイド１１上に取り付けられている。【０００６】
サーボモータ３の回転運動は減速機４を介して倍力機構としてのボールネジ５に伝えられ、ボールネジ５の回転運動はボールネジ軸上のナット７の直線運動に変換され、可動部８を介してスクリュ９の前後移動とスクリュ９による溶融樹脂への加圧が実現される。また、スクリュ９の位置はサーボモータ軸に設けられたモータエンコーダ１２で検出される。スクリュ９が射出シリンダ２の先端部の溶融樹脂に与える圧力は、ボールネジのナット７と可動部８との間に挟まれて設置される圧力検出センサ（例えば、ロードセル）１０により検出される。溶融樹脂は、スクリュ９により金型１内のキャビティ１３に充填される。【０００７】
成形品を作る射出工程は、充填工程と保圧工程とからなる。充填工程では金型内に充填された樹脂温度を一様にするために、溶融樹脂を出来るだけ高速、短時間でキャビティ１３内に充填することが必要である。しかし、過度の高速化により、過度の射出圧力を発生させると成形不良の要因となるので、射出圧力はある所定の値以下という制約が課される。充填後の保圧工程では、溶融樹脂の冷却による樹脂収縮を補充するために一定時間圧力を掛ける必要がある。したがって、射出工程中の射出速度及び射出圧力に対しては、次の２つの要求課題が与えられる。【０００８】
（１）       充填工程中は、所定の射出速度を実現すると同時に射出圧力が所定の圧力以下であること。
（２）       保圧工程中は、所定の保持圧力を実現すると同時に射出速度が安全面の点から所定の速度以下であること。【０００９】
つまり、図４(ａ)に示す充填工程中(時間０〜ｔ１)では、所定の射出速度を実現するために図４(ｂ)のように与えられた射出速度指令を忠実に実行させるような射出速度制御がなされる。しかし、射出圧力は図４(ｃ)のように与えられた射出圧力指令以下に抑える制御がなされなければならない。ここで、図４(ｂ)、図４(ｃ)の縦軸に示す１００％はそれぞれ射出速度、射出圧力の最大値を示す。【００１０】
次に保圧工程(時間ｔ１〜ｔ２)に入ったときは、所定の射出圧力を実現するために図４(ｃ)のように与えられた射出圧力指令を忠実に実行させるような射出圧力制御がなされる。しかし、射出速度は図４(ｂ)のように与えられた射出速度指令以下に抑える制御がなされなければならない。【００１１】
図５は、前述の２つの要求課題(１)(２)(段落(０００８))を実現するための制御装置の内部構成を説明するためのブロック図である(特許文献１)。この制御装置は、射出制御器２０とモータ制御器（サーボアンプ）４０とで構成される。【００１２】
射出制御器２０について説明する。射出制御器２０は一定時間間隔Δｔ毎に制御演算を行い、制御指令を更新する。射出制御器２０は、射出速度指令器２１、変換器２２、パルス発生器２３、アナログ／デジタル（Ａ／Ｄ）変換器２５、射出圧力指令器２６、減算器２７、圧力制御器２８、デジタル／アナログ（Ｄ／Ａ）変換器２９から構成される。射出制御器２０には、前記の圧力検出センサ１０が接続される。【００１３】
射出速度指令器２１は、時間シーケンスに従った射出速度指令Ｖ＊ｉを変換器２２に出力する。変換器２２は、入力された射出速度指令Ｖ＊ｉに対してΔｔ時間に進むべきスクリュ９の移動量を、スクリュ移動指令Δｘ＊ｖとして次の式(１)により算出する。【００１４】
【数１】
2010134359.tif 000003
【００１５】
算出したスクリュ移動量指令Δｘ＊ｖをパルス発生器２３に出力する。【００１６】
パルス発生器２３は、スクリュ移動量指令Δｘ＊ｖに相当する数のパルスを含むパルス列２４を発生させる。そして、この発生させたパルス列２４をモータ制御器４０内のカウンタ４１に出力する。【００１７】
また、圧力検出センサ１０で検出された圧力を射出圧力センサ信号ＰｉとしてＡ／Ｄ変換器２５を介して射出制御器２０に取り込む。Ａ／Ｄ変換器２５は取り込んだ射出圧力センサ信号Ｐｉを減算器２７へ出力する。【００１８】
射出圧力指令器２６は、時間シーケンスに従った射出圧力指令Ｐ＊ｉを減算器２７に入力する。減算器２７は射出圧力指令Ｐ＊ｉと射出圧力センサ信号Ｐｉとの制御偏差ΔＰｉを次の式(２)により算出する。【００１９】
【数２】
2010134359.tif 000004
【００２０】
この後、算出した制御偏差ΔＰｉを圧力制御器２８に入力する。【００２１】
圧力制御器２８は、制御偏差ΔＰｉに対して比例積分微分(ＰＩＤ)制御演算を実行してモータ電流指令ｉ＊ｐを算出する。この算出したモータ電流指令ｉ＊ｐをＤ／Ａ変換器２９を介してモータ制御器４０に入力する。モータ制御器４０では、内臓のＡ／Ｄ変換器４２によりモータ電流指令ｉ＊ｐを取り込んだ後、このモータ電流指令ｉ＊ｐを比較器４３に入力する。【００２２】
次にモータ制御器（サーボアンプ）４０について説明する。モータ制御器４０はカウンタ４１及び４４、A／D変換器４２、比較器４３、減算器４５及び４８、位置制御器４６、微分回路４７、速度制御器４９、パルス幅変調制御(PWM)回路５０から構成される。このモータ制御器４０は前記サーボモータ３に接続されている。前記サーボモータ３には前記モータエンコーダ１２が取り付けられている。【００２３】
射出制御器２０からの速度指令用のパルス列２４を受けるモータ制御器４０内のカウンタ４１は、パルス列２４のパルス数を積算計数して、スクリュ９の位置指令ｘ＊を検出する。このｘ＊を減算器４５に入力する。また、カウンタ４４はモータエンコーダ１２からのパルスを受ける。カウンタ４４はこのパルス数を積算計数して、スクリュ９の実際の位置ｘを検出する。このｘを減算器４５に入力する。【００２４】
減算器４５は入力されたｘ＊及びｘから位置制御偏差(ｘ＊−ｘ)を算出した後、この結果を位置制御器４６に入力する。位置制御器４６は次の式(３)より速度指令ｖ＊を算出した後、その結果を減算器４８に入力する。【００２５】
【数３】
2010134359.tif 000005
【００２６】
ここで、ＫＰは位置制御器４６の比例ゲインである。【００２７】
モータエンコーダ１２はパルスを、カウンタ４４の他に微分回路４７にも入力する。微分回路４７はこのパルスにより実速度ｖを検出する。このｖを減算器４８に入力する。【００２８】
減算器４８は入力されたｖ＊及びｖから速度制御偏差(ｖ＊−ｖ)を算出して速度制御器４９に入力する。速度制御器４９は次の式(４)よりモータ電流指令ｉ＊ｖを算出し、比較器４３に入力する。【００２９】
【数４】
2010134359.tif 000006
【００３０】
ここで、ＫＰｖは速度制御器４９の比例ゲインであり、ＴＩｖは速度制御器４９の積分時間である。モータ制御器４０内では、位置制御ループのマイナループとして速度制御ループが形成されている。【００３１】
速度制御器４９からの出力であるモータ電流指令ｉ＊ｖと圧力制御器２８からの出力であるモータ電流指令ｉ＊ｐとが入力される比較器４３は、ｉ＊ｖとｉ＊ｐの小さい方の値を選択し、その値を最終的なモータ電流指令ｉ＊としてＰＷＭ回路５０に入力する。ＰＷＭ回路５０は与えられたモータ電流指令ｉ＊に基づいて所定の３相電圧をサーボモータ３に印加する。これにより、サーボモータ３はモータ電流ｉ＊で駆動される。比較器４３は位置制御ループにより決定されるモータ電流指令ｉ＊ｖを制御ループ外からのモータ電流指令ｉ＊ｐによって制限するために設けられている。【００３２】
次に、前記比較器４３の動作により前述した２つの要求課題(１)(２)(段落(０００８))が実現されることを、図４を用いて説明する。
まず、図４の充填工程の終了時間であるｔ１が、実際にキャビティ１３に樹脂が十分充填される時間よりも前に設定されていた場合を考える。このときは、キャビティ１３内にはまだ空隙が存在するので射出圧力は一般的に低くなっている。この状態から次の保圧工程に入ると、一般に射出圧力指令Ｐ＊ｉの方が実際の射出圧力よりも高いので、実際の射出圧力を高くするために圧力制御器２８から出力されるモータ電流指令ｉ＊ｐが大きくなる。【００３３】
ここで、モータ電流指令ｉ＊ｐが最終的なモータ電流指令として採用されたとすると、キャビティ１３内が未充填なので射出速度が上昇する。これにより、図４(ｂ)の射出速度指令を上回ってしまうおそれがある。しかし、比較器４３では常にｉ＊ｖとｉ＊ｐのうちで小さいほうの値が選択されるので、充填中に保圧制御に移行してｉ＊ｐがｉ＊ｖを上回っても、比較器４３においてｉ＊ｖが最終的なモータ電流指令として採用される。
つまり、射出圧力制御から射出速度制限に移行させることができ、前述の要求課題(２)を常に満足させることができる。【００３４】
次に、図４の充填工程の終了時間ｔ１が、実際にキャビティ１３に樹脂が十分充填される時間よりも後に設定されていた場合を考える。このときはキャビティ１３内に樹脂が十分充填されていても、まだ充填工程中であるので射出速度制御を実行する。しかし、この時点ではキャビティ１３が十分充填されているので実際のスクリュ速度はかなり低下している。そのため、速度制御器４９から出力されるモータ電流指令ｉ＊ｖが大きくなる。【００３５】
ここで、モータ電流指令ｉ＊ｖが最終的なモータ電流指令として採用されたとすると、既にキャビティ１３内に樹脂が十分充填されているので急激に圧力が上昇する。これにより、図４(ｃ)の射出圧力指令を上回ってしまうおそれがある。しかし、比較器４３では常にｉ＊ｖとｉ＊ｐのうちで小さい方の値が選択されるので、充填が完了し、ｉ＊ｖがｉ＊ｐを上回っても、比較器４３においてｉ＊ｐが最終的なモータ電流指令として採用される。つまり、射出速度制御から射出圧力制限に移行させることができ、前述の要求課題(１)を常に満足させることができる。【００３６】
以上述べた電動射出成形機の制御装置では、樹脂の圧力検出センサの装着を不可欠としていたが、特許文献２〜特許文献９は圧力検出センサを使わずに射出成形機での圧力制御を実現する考案である。【００３７】
特許文献２は、油圧式射出成形機を対象とし、樹脂圧力、樹脂温度、樹脂比容積の間の関係式を与える樹脂特性式を用いて、計測樹脂温度及び製品重量の目標値から決まる樹脂比容積を入力して必要となる樹脂圧力を算出する。さらに保圧工程移行時の初期金型温度と初期樹脂温度の計測値及び算出した前記樹脂圧力を用いて近似式より必要な保圧力設定値を求める。前記保圧力設定値を電圧指令値に変換して油圧サーボ弁アンプに入力し、油圧シリンダの油圧力により保圧設定値を実現している。【００３８】
特許文献３は、金型内キャビティの圧力を検出するために、キャビティに出入自在に設けられる突出しピンにキャビティ内圧力が掛かる構造で、突出しピンを支持するボールネジ機構を介してボールナットを回転駆動する突出し用電動機からなる。射出・保圧時のキャビティ内圧力を受けた突出しピンの位置変化を突出し用電動機に付属する回転位置検出器で検出し、射出前の初期突出し位置に突出しピンを保持するために必要な突出し用電動機の負荷電流を電流検出器で検出し、前記負荷電流が変換されてキャビティ内圧力を検知するものである。【００３９】
特許文献４は、金型内圧力を検出するのに、突出ピンを移動させる突出用サーボモータシステムに対してモータ速度信号とトルク指令信号を入力信号とし、負荷トルクを推定する外乱オブザーバを構成する。外乱オブザーバの演算内容が示されている。負荷トルク推定値から等価なキャビティ内圧力を得る。突出ピンを位置保持した状態でのサーボモータの駆動電流或はトルク指令値から直接キャビティ内圧力を得る方法も示す。【００４０】
特許文献５は、予め射出スクリュの駆動圧力と射出速度からキャビティ内樹脂圧力を推定するキャビティ内樹脂圧力推定関数を決定しておき、制御動作時には、前記キャビティ内樹脂圧力推定関数に成形品形状情報、樹脂情報、射出スクリュ駆動圧力及び射出速度のリアルタイム検出値を入力してキャビティ内樹脂圧力推定値をリアルタイムで取得し、前記キャビティ内樹脂圧力推定値と設定された基準圧力との偏差により充填速度を自動的に制御する方法である。高精度の推定関数を求める手順が示されている。【００４１】
特許文献６は射出用モータを対象に、樹脂圧推定値を求めるオブザーバを構成して、オブザーバ出力である樹脂圧推定値を樹脂圧検出値として用いる射出圧力制御装置を実現している。オブザーバは射出工程中のモータ回転速度、モータトルク指令及び射出機構全体の摩擦抵抗を入力信号として、モータ回転速度推定値及び樹脂圧力推定値を出力する。オブザーバで使う制御対象は次の２つを考えている。
（１）   射出用モータがスクリュに直結している場合
（２）   射出用モータが１段プーリ減速機を介してスクリュを駆動する場合
前記オブザーバモデル(１)での摩擦抵抗は、動摩擦抵抗と静止摩擦抵抗からなる。前記オブザーバモデル(２)での摩擦抵抗は、速度依存成分と荷重依存成分からなる動摩擦抵抗からなる。【００４２】
上記オブザーバモデル(１)での樹脂圧力は一定としている。上記オブザーバモデル(２)では、モータ回転速度推定値、樹脂圧力推定値の他に、負荷側プーリ回転速度推定値、ベルト張力推定値、スクリュが樹脂を押す力推定値も、オブザーバが出力する。ベルトを弾性体とし、樹脂圧力の時間変化は、負荷側プーリ回転速度及び回転加速度に比例し、さらにスクリュが樹脂を押す力に比例するとしている。スクリュが樹脂を押す力は一定としている。【００４３】
特許文献７は射出用モータを対象に、溶融樹脂から受ける反抗トルク(圧力推定値)を求めるオブザーバを構成して射出速度・圧力制御装置を実現している。オブザーバはモータ回転速度及びモータ電流指令値を入力信号として、モータ回転速度推定値、モータ軸と負荷軸の回転位置の差推定値及び外力トルク推定値を出力する。オブザーバの制御対象は、モータが1段プーリ減速機を介してスクリュを駆動するモデルを考えている。外力トルク推定値から換算した圧力を圧力検出値に使っている。【００４４】
特許文献８は、特許文献６の制御対象モデル(２)を用いたオブザーバを構成し、状態変数であるモータ回転速度、負荷側プーリ回転速度推定値、ベルト張力推定値及び樹脂圧力推定値をフィードバック信号として使って、樹脂圧力を圧力設定値に追従させる制御方式を考案している。また、モータトルク指令値を前記４つの状態フィードバックで決める制御方式を考案している。【００４５】
特許文献９はモータ回転速度とモータトルクから、モータが受ける負荷トルクを伝達関数の逆モデルから得る方法を考案している。逆モデルは高次微分作用を必要とする。【００４６】
圧力検出センサを使わずに、樹脂圧力を検知する特許文献２〜特許文献９の目的は、次のような不利を改善しようとする点で共通している。
（１）   高圧環境下で信頼性の高い圧力センサは高価になる。
（２）   金型内やスクリュノズル部への圧力センサ取付けは、特別な加工を施す必要があり、作業コストが無視できない。
（３）   電動モータから射出スクリュに至る射出軸系に取り付けるロードセルは、組み込むための機械構造を複雑にし、更には射出軸系の機械剛性の低下を招く。
（４）   歪みゲージを検出部に使用するロードセルでは、微弱なアナログ信号に対するノイズ対策が必要になり、また信号アンプのゼロ点調整やスパン調整等にも人手による作業が必要になる（特許文献１０）。【００４７】
他に、特許文献１１、特許文献１２は電動射出成形機の圧力制御に関する考案であるが、いずれも圧力検出センサの使用が不可欠である。特許文献１１は、圧力の制御をスクリュの変位の制御で行うという考え方に基づく式（明細書(１)式）で導入された無効速度ω１なる概念を導入して、モータトルクと圧力設定値とのギャップを生じさせる原因となる非線形損失による圧力損失を補償するパラメータとして無効速度ω１を使って正確な圧力制御を行うという考案である。外乱オブザーバが、ロードセルによる検出圧力と自らが出力する圧力推定値との差が入力されて、その差が０になるように無効速度推定値を出力する。特許文献１２は、特許文献１１の先願で、外乱オブザーバの構成が異なる。
【先行技術文献】
【特許文献】
【００４８】
【特許文献１】特許第３７８７６２７号（図５）
【特許文献２】特開平５−７７２９８号公報
【特許文献３】特開平６−８５６号公報
【特許文献４】特開平７−２９９８４９号公報
【特許文献５】特開平９−２７７３２５号公報
【特許文献６】再公表特許ＷＯ２００５／０２８１８１号公報
【特許文献７】特開２００６−１４２６５９号公報
【特許文献８】特開２００６−２５６０６７号公報
【特許文献９】特開２００８−２６５０５２号公報
【特許文献１０】特開２００３−２１１５１４号公報
【特許文献１１】特開平１０−２４４５７１号公報
【特許文献１２】特開平１０−４４２０６号公報
【非特許文献】
【００４９】
【非特許文献１】H.K.Khalil, Nonlinear Systems, 14.5 High-Gain Observers, Prentice-Hall,(2002), pp.610-625
【非特許文献２】B.D.O. Anderson andJ.B. Moore, Optimal Control, Linear Quadratic Methods, 7.2 DeterministicEstimator Design, Prentice-Hall, (1990), pp.168-178
【非特許文献３】A.M. Dabroom andH.K. Khalil, Discrete-time implementation of high-gain observers for numericaldifferentiation, Int. J. Control, Vol.72, No.17, (1999), pp.1523-1537
【非特許文献４】A.M. Dabroom andH.K. Khalil, Output Feedback Sampled-Data Control of Nonlinear Systems UsingHigh-Gain Observers, IEEE Trans. Automat. Contr., Vol.46, No.11, (2001),pp.1712-1725
【発明の概要】
【発明が解決しようとする課題】
【００５０】
圧力検出センサを使うことによる「背景技術」で述べた前記４つの不利(段落(００４６))を回避するために、圧力検出センサを使わずに電動射出成形機の制御装置に要求される機能（「背景技術」で述べた２つの要求課題(1)(2)(段落(０００８))）を達成できる制御装置を実現すること。
【課題を解決するための手段】
【００５１】
成形品を作る射出工程は充填工程と保圧工程からなるが、「背景技術」で述べた要求課題(１)より充填工程中も所定の射出圧力以下に保持することが求められることから、充填工程中の実際の射出圧力を時間的に遅れることなく検知する必要がある。更に要求課題(２)より保圧工程中に所定の保持圧力を実現する良好な保持圧制御を行うには、実際の射出圧力を時間的遅れなく検出する必要がある。したがって、射出工程中の圧力検知手段には時間的遅れのないことが要求される。【００５２】
次に検知圧力の誤差は、成形不良や安全運転の欠如につながることから、高精度の圧力検知が要求される。したがって、前記２つの必要条件
(A)時間的遅れが非常に小さい
(B)高精度である
を満たす圧力検知を実現するためには、高ゲイン観測器(非特許文献１)の手法を採用する。計測できる変数を入力して、すべての状態変数を推定する高ゲイン観測器が、前記２つの必要条件(A)(B)を満たしていることを簡単な制御対象モデルを使って説明する。式(５)は制御対象モデルの状態方程式を示す。【００５３】
【数５】
2010134359.tif 000007
【００５４】
ここで、ｘ１、ｘ２は状態変数、ｕは入力変数、ｙは出力変数である。φ(ｘ、ｕ)は変数ｘ、ｕからなる非線形関数である。例えば、ｘ１は位置変数、ｘ２は速度変数、ｕはモータ電流である。出力ｙ及び入力ｕは計測できるとし、状態ｘを推定する高ゲイン観測器は、式(６)で与えられる。【００５５】
【数６】
2010134359.tif 000008
【００５６】
ここで、ｘ^１、ｘ^２は、状態変数ｘ１、ｘ２の推定値を表す。Ｈ１、Ｈ２は高ゲイン観測器のゲイン定数で一般に１より大きな定数が与えられる。関数φ０は、高ゲイン観測器の演算に採用された関数φの公称（基準）関数を表す。式(６)の高ゲイン観測器を用いたときの推定値誤差ｘ〜１、ｘ〜２は、式(５)、式(６)より式(７)で与えられる。【００５７】
【数７】
2010134359.tif 000009
【００５８】
ここで、δは、実際には得られない真の関数φと高ゲイン観測器で採用した公称の関数φ０ との差、すなわち制御対象のモデル誤差と考えられる。次に１より十分小さい正のパラメータεを導入して、Ｈ１、Ｈ２を式(９)で与える。【００５９】
【数８】
2010134359.tif 000010
【００６０】
定数Ｈ１、Ｈ２には、式(９)から判るように大きなゲイン定数を採用することから高ゲイン観測器と呼ばれる。式(９)を使うと、式(７)は式(１０)で表される。【００６１】
【数９】
2010134359.tif 000011
【００６２】
推定値誤差ｘ〜１、ｘ〜２を式(１１)で表される新変数η１、η２に置き換える。【００６３】
【数１０】
2010134359.tif 000012
【００６４】
式(１１)より、式(１０)は式(１２)で表される。【００６５】
【数１１】
2010134359.tif 000013
【００６６】
ここで、パラメータεを十分小さくとれば、式(１２)より推定誤差η１、η２は制御対象のモデル誤差δから受ける影響を十分小さくできる。すなわち高ゲイン観測器を用いれば、状態変数に射出圧力を含める制御対象モデルを採用することにより、圧力検知に要求される前記必要条件(B)「高精度である」ことが満たされることが判る。【００６７】
次にモデル誤差の影響を無視すると、式(１２)は式(１３)、式(１４)で表される。【００６８】
【数１２】
2010134359.tif 000014
【００６９】
行列Aの共役複素数の固有値λ１、λ２の実数部Ｒｅ(λ１)＝Ｒｅ(λ２)が負になるようにＫ１、Ｋ２を決めると、式(１３)より推定値誤差η１、η２はその初期値η１０、η２０に対して式(１５)で与えられる。【００７０】
【数１３】
2010134359.tif 000015
【００７１】
ｔは時間変数で、C１(ｔ)〜C４(ｔ)はＫ１、Ｋ２で決まる一定振幅、一定周波数成分を表す。Ｒｅ(λ１)＜０で、パラメータεを１より十分小さくとれば、式(１５)より推定値誤差η１(ｔ)、η２(ｔ)は急速に０になることが判る。すなわち高ゲイン観測器を用いれば、圧力検知に要求される前記必要条件(A)「時間的遅れが非常に小さい」ことが満たされることが判る。【００７２】
式(６)の高ゲイン観測器では、すべての状態変数ｘ１、ｘ２の推定値を得たが、状態変数ｘ１は出力ｙとして計測できるので、状態変数ｘ２だけを推定すればよい。このときの高ゲイン観測器は、式(１６)で与えられる。【００７３】
【数１４】
2010134359.tif 000016
【００７４】
ここで、Hは高ゲイン観測器のゲイン定数で、１より大きな定数が与えられる。式(１６)は右辺に出力ｙの時間微分を含むので、直接、演算式としては使えないが、式(１６)で与えられる高ゲイン観測器は、前記２つの必要条件(A)(B)(段落(００５２))を満たすことを示す。式(５)の３番目の式より、式(１７)が得られる。【００７５】
【数１５】
2010134359.tif 000017
【００７６】
式(１６)、式(１７)より式(１８)が得られる。【００７７】
【数１６】
2010134359.tif 000018
【００７８】
式(５)の２番目の式を使うと、式(１８)より式(１９)を得る。【００７９】
【数１７】
2010134359.tif 000019
【００８０】
次に１より十分小さい正のパラメータεを導入してHを式(２１)で与える。【００８１】
【数１８】
2010134359.tif 000020
【００８２】
式(２１)を使うと、式(１９)は式(２２)で表される。【００８３】
【数１９】
2010134359.tif 000021
【００８４】
ここで、パラメータεを十分小さくとれば、式(２２)より推定値誤差ｘ〜２は、制御対象のモデル誤差δから受ける影響を十分小さくできる。すなわち高ゲイン観測器を用いれば、状態変数に射出圧力を含める制御対象モデルを採用することにより圧力検知に要求される前記必要条件(B)「高精度である」ことが満たされることが判る。【００８５】
次にモデル誤差δの影響を無視すると、式(２２)は式(２３)で表される。【００８６】
【数２０】
2010134359.tif 000022
【００８７】
式(２３)より推定値誤差ｘ〜２は、式(２４)で表される。【００８８】
【数２１】
2010134359.tif 000023
【００８９】
ここで、ｘ〜２０は推定値誤差ｘ〜２の初期値である。式(２４)より、パラメータεを１より十分小さくとれば、推定値誤差ｘ〜２(ｔ)は急速に０になることが判る。すなわち高ゲイン観測器を用いれば、圧力検知に要求される前記必要条件(A)「時間的遅れが非常に小さい」ことが満たされることが判る。式(１６)の高ゲイン観測器は、計測できる状態変数は推定せず、必要最小限の状態変数を推定するので式(６)の観測器より次数が低くなるので、低次元高ゲイン観測器と呼ばれる。【００９０】
次に式(１６)の演算を出力ｙの時間微分を使わずに行う方法を示す。式(２５)で与えられる新変数ｗ^を導入する。【００９１】
【数２２】
2010134359.tif 000024
【００９２】
式(２５)を使うと、式(１６)は式(２６)で与えられる。【００９３】
【数２３】
2010134359.tif 000025
【００９４】
式(２６)よりｗ^を計算し、式(２７)より推定値ｘ^２を求めることができる。【００９５】
【数２４】
2010134359.tif 000026
【００９６】
状態変数に射出圧力を含めた電動射出成形機の制御対象モデルに対して、高ゲイン観測器を適用する手順は後述の「実施例１」、「実施例２」で詳述する。
【発明の効果】
【００９７】
状態変数に射出圧力を含めた電動射出成形機の制御対象モデルに高ゲイン観測器を適用することにより、圧力検出センサを使用しないで、時間的遅れが非常に小さく、高精度な圧力検知が可能となる。これにより、圧力検出センサを使用しないで、電動射出成形機の制御装置に要求される機能を実現でき、さらに「背景技術」で述べた４つの不利(段落(００４６))を回避することができる。
【図面の簡単な説明】
【００９８】
【図１】本発明に係る圧力制御装置及び圧力制御方法の全体構成を示す実施例１の説明図である。
【図２】本発明に係る圧力制御装置及び圧力制御方法の全体構成を示す実施例２の説明図である。
【図３】電動射出成形機の従来の射出・保圧機構の構成を示す模式図である。
【図４】電動射出成形機の射出工程の時間的スケジュールを示す説明図である。
【図５】従来の電動射出成形機の圧力制御装置及び圧力制御方法の全体構成を示す説明図である。
【図６】本発明に係る電動射出成形機の射出・保圧機構の構成を示す模式図である。
【図７】本発明に係る実施例１での高ゲイン観測器の射出圧力推定シミュレーション結果の説明図である。
【図８】本発明に係る実施例１での高ゲイン観測器の射出速度推定シミュレーション結果の説明図である。
【図９】本発明に係る実施例２での高ゲイン観測器の射出圧力推定シミュレーション結果の説明図である。
【発明を実施するための形態】
【００９９】
以下、本発明の実施形態に係る電動射出成形機の制御装置を図面に基づいて説明する。
【実施例１】
【０１００】
図６は、圧力検出センサを使用しない射出・保圧機構の模式図を示す。図６は、圧力検出センサを除いて図３と同じ符号を持つ部品から構成されるので、図６の説明は、「背景技術」で述べた図３の説明(段落(０００４)〜(０００６))に代える。【０１０１】
図１は、本発明の一実施形態に係る高ゲイン観測器による射出圧力検知を電動射出成形機の制御装置に適用した例で、制御装置の内部構成を説明するためのブロック図である。この制御装置は、高ゲイン観測器３１を内蔵する射出制御器２０とモータ制御器(サーボアンプ)４０とから構成される。【０１０２】
射出制御器２０について説明する。射出制御器２０は一定時間間隔Δｔ毎に制御演算を行い、制御指令をモータ制御器４０に出力する。射出制御器２０は、射出速度指令器２１、変換器２２、パルス発生器２３、射出圧力指令器２６、減算器２７、圧力制御器２８、デジタル／アナログ（Ｄ／Ａ）変換器２９、アナログ／デジタル（Ａ／Ｄ）変換器３０及び高ゲイン観測器３１から構成される。【０１０３】
射出速度指令器２１は、時間シーケンスに従った射出速度指令Ｖ＊ｉを変換器２２に出力する。変換器２２は、入力された射出速度指令Ｖ＊ｉに対してΔｔ時間に進むべきスクリュ９の移動量を、スクリュ移動指令Δｘ＊ｖとして式(２８)により算出する。【０１０４】
【数２５】
2010134359.tif 000027
【０１０５】
算出したスクリュ移動量指令Δｘ＊ｖをパルス発生器２３へ出力する。パルス発生器２３は、スクリュ移動量指令Δｘ＊ｖに相当する数のパルスを含むパルス列２４を発生させる。このパルス列２４は、モータ制御器４０内のカウンタ４１へ出力される。【０１０６】
射出圧力指令器２６は、時間シーケンスに従った射出圧力指令Ｐ＊ｉを減算器２７へ出力する。モータ制御器４０内のモータ電流指令ｉ＊は、射出制御器２０内のＡ／Ｄ変換器３０を介して高ゲイン観測器３１に入力される。更にモータ制御器４０内のカウンタ４４の出力であるスクリュ９の位置信号ｘが高ゲイン観測器３１に入力される。高ゲイン観測器３１は、入力信号ｘ及びｉ＊を使って、射出軸系の数式モデルを利用して導出した内蔵する離散演算式を実行して、射出圧力推定値Ｐ^ｉ及び図示しない射出速度推定値ｖ^を出力する。【０１０７】
射出圧力推定値Ｐ^ｉは、減算器２７に入力される。減算器２７は射出圧力指令Ｐ＊ｉと射出圧力推定値Ｐ^ｉとの制御偏差ΔＰｉを式(２９)より算出する。【０１０８】
【数２６】
2010134359.tif 000028
【０１０９】
減算器２７は、算出した制御偏差ΔＰｉを圧力制御器２８へ出力する。【０１１０】
圧力制御器２８は、制御偏差ΔＰｉに対して比例積分微分(ＰＩＤ)制御演算を実行してモータ電流指令ｉ＊ｐを算出する。モータ電流指令ｉ＊ｐはＤ／Ａ変換器２９を介してモータ制御器４０へ出力される。モータ制御器４０ではＡ／Ｄ変換器４２を介してモータ電流指令ｉ＊ｐを取込んだ後、モータ電流指令ｉ＊ｐは比較器４３に入力される。【０１１１】
次にモータ制御器（サーボアンプ）４０について説明する。モータ制御器４０はカウンタ４１及び４４、A／D変換器４２、比較器４３、減算器４５及び４８、位置制御器４６、微分回路４７、速度制御器４９、パルス幅変調制御（PWM）回路５０から構成される。モータ制御器４０には、前記サーボモータ３が接続され、サーボモータ３には前記モータエンコーダ１２が取り付けられている。【０１１２】
射出制御器２０からの速度指令用のパルス列２４を受けるカウンタ４１は、パルス列２４のパルス数を積算計数して、スクリュ９の位置指令ｘ＊を算出し、減算器４５へ出力する。また、カウンタ４４はモータエンコーダ１２からのパルスを受けて積算計数して、スクリュ９の実際の位置ｘを検出する。このｘを減算器４５へ出力する。【０１１３】
減算器４５は、入力されたｘ＊及びｘから位置制御偏差(ｘ＊−ｘ)を算出し、位置制御器４６へ出力する。位置制御器４６は式(３０)より速度指令ｖ＊を算出し、減算器４８へ出力する。【０１１４】
【数２７】
2010134359.tif 000029
【０１１５】
ここで、ＫＰは位置制御器４６の比例ゲインである。モータエンコーダ１２はパルスを、カウンタ４４の他に微分回路４７へも出力する。微分回路４７は、このパルスにより実速度ｖを検出し、減算器４８へ出力する。【０１１６】
減算器４８は入力されたｖ＊及びｖから速度制御偏差(ｖ＊−ｖ)を算出して速度制御器４９へ出力する。速度制御器４９は式(３１)よりモータ電流指令ｉ＊ｖを算出し、比較器４３へ出力する。【０１１７】
【数２８】
2010134359.tif 000030
【０１１８】
ここで、ＫＰｖは速度制御器４９の比例ゲインであり、ＴＩｖは速度制御器４９の積分時間である。【０１１９】
速度制御器４９からのモータ電流指令ｉ＊ｖと圧力制御器２８からのモータ電流指令ｉ＊ｐとが入力される比較器４３は、ｉ＊ｖとｉ＊ｐの小さい方の値を選択し、その値を最終的なモータ電流指令ｉ＊としてＰＷＭ回路５０に入力する。ＰＷＭ回路５０は与えられたモータ電流指令ｉ＊に基づいて所定の３相電圧をサーボモータ３に印加する。これにより、サーボモータ３はモータ電流指令ｉ＊で駆動される。【０１２０】
比較器４３の動作により、「背景技術」で述べた電動射出成形機の制御装置に要求される前記2つの要求課題(１)(２)(段落(０００８))が実現されることは、「背景技術」の段落(００３２)〜(００３５)で詳述したので、ここでは述べない。【０１２１】
次に、位置信号ｘ及びモータ電流指令ｉ＊を入力して射出圧力推定値Ｐ^ｉ及び射出速度推定値ｖ^を出力する高ゲイン観測器３１の設計に必要になる射出軸系の数式モデルを図６を使って説明する。図６のモータ３の運動方程式は式(３２)で与えられる。【０１２２】
【数２９】
2010134359.tif 000031
【０１２３】
ここで、ＪＭはモータ本体慣性モーメント、ＪＧ１はモータ側減速歯車慣性モーメント、ωｍはモータ角速度、ＴＭはモータトルク、ｒ１はモータ側減速歯車半径、Ｆは減速機伝達力及びｔは時間である。ボールネジ５の運動方程式は式(３３)で与えられる。【０１２４】
【数３０】
2010134359.tif 000032
【０１２５】
ここで、ＪＳはボールネジ軸慣性モーメント、ＪＧ２は負荷側減速歯車慣性モーメント、ωsはボールネジ軸角速度、ｒ２は負荷側減速歯車半径及びＴａはボールネジ駆動トルクである。可動部８の運動方程式は式(３４)、式(３５)で与えられる。【０１２６】
【数３１】
2010134359.tif 000033
【０１２７】
ここで、Ｗは可動部重量、ｇは重力加速度、ｖはスクリュ（可動部）速度、ｘはスクリュ位置(射出開始時ｘ＝０)、Ｆａはボールネジ軸力、FＬはスクリュが樹脂から受ける負荷力、μは可動部―リニアガイド摩擦係数である。ボールネジ駆動トルクＴａとボールネジ軸力Ｆａの関係は式(３６)で与えられる。【０１２８】
【数３２】
2010134359.tif 000034
【０１２９】
ここで、ｌはボールネジリード及びηはボールネジ効率である。スクリュ速度ｖ、ボールネジ角速度ωs及びモータ角速度ωｍの関係は式(３７)で与えられる。【０１３０】
【数３３】
2010134359.tif 000035
【０１３１】
スクリュの受ける負荷力ＦＬは式(３８)で与えられる。【０１３２】
【数３４】
2010134359.tif 000036
【０１３３】
ここで、Ａｓはスクリュ断面積、Ｐｉは射出シリンダ貯留部（スクリュ先端部）の樹脂圧力を意味する射出圧力、Ｃｍｔは射出シリンダ粘性係数、αは速度べき乗数である。射出圧力の方程式は式(３９)、式(４０)で与えられる。【０１３４】
【数３５】
2010134359.tif 000037
【０１３５】
ここで、Ｖｉはシリンダ貯留部容積、Ｖｉ０はシリンダ貯留部容積初期値（射出開始時）、Ｑｉｎは金型への樹脂流入量及びβは樹脂体積弾性係数である。モータ特性は式(４１)で与えられる。【０１３６】
【数３６】
2010134359.tif 000038
【０１３７】
ここで、ＫＴはモータトルク係数及びｉｍはモータ電流である。式(３２)、式(３３)、式(３７)を使ってωs、Ｆを消去すると、式(４２)を得る。【０１３８】
【数３７】
2010134359.tif 000039
【０１３９】
式(３４)、式(３６)、式(３７)、式(４２)を使って、Ｔａ、Ｆａを消去すると式(４３)を得る。【０１４０】
【数３８】
2010134359.tif 000040
【０１４１】
式(４３)は、モータ軸に換算した射出軸系の運動方程式を表し、式(４４)はモータ軸換算等価慣性モーメントを表す。式(３５)、式(３７)より式(４５)が得られる。【０１４２】
【数３９】
2010134359.tif 000041
【０１４３】
式(３８)、式(４１)、式(４３)より、射出軸系の運動方程式は式(４６)で与えられる。【０１４４】
【数４０】
2010134359.tif 000042
【０１４５】
次に式(４０)は、式(４７)で表される。【０１４６】
【数４１】
2010134359.tif 000043
【０１４７】
ここで、ｘｍａｘは最大スクリュ位置である。式(３７)と式(４７)を使うと、式(３９)は式(４８)で表される。【０１４８】
【数４２】
2010134359.tif 000044
【０１４９】
次に変数の無次元化を行う。式(４５)の変数を無次元化すると、式(４９)が得られる。【０１５０】
【数４３】
2010134359.tif 000045
【０１５１】
ここで、ωｍａｘはモータ定格回転数、ｖｍａｘは最大射出速度である。【０１５２】
次に式(４６)の変数を無次元化すると、式(５０)が得られる。【０１５３】
【数４４】
2010134359.tif 000046
【０１５４】
ここで、ｉｍａｘはモータ定格電流である。Ｐｍａｘは最大射出圧力を表す。式(５０)では、次の関係式(５１)を使っている。【０１５５】
【数４５】
2010134359.tif 000047
【０１５６】
式(５０)は、式(５２)で表される。【０１５７】
【数４６】
2010134359.tif 000048
【０１５８】
ここで、ＴＭｍａｘ＝ＫＴｉｍａｘでモータ定格トルクを表す。【０１５９】
次に式(４８)の変数を無次元化すると、式(５３)が得られる。【０１６０】
【数４７】
2010134359.tif 000049
【０１６１】
ここで、Qｍａｘ＝Ａｓｖｍａｘで最大射出量を表す。式(５３)は式(５４)に書き換えられる。【０１６２】
【数４８】
2010134359.tif 000050
【０１６３】
ここで、金型への樹脂流入量[Ｑｉｎ／Ｑｍａｘ]は一般的に射出圧力[Ｐｉ／Ｐｍａｘ]の関数で、その関数関係はノズル形状、金型入口形状、キャビティ形状及び樹脂特性で決まる。【０１６４】
【数４９】
2010134359.tif 000051
【０１６５】
以上高ゲイン観測器３１の設計に必要な射出軸系の数式モデルは、式(４９)、式(５２)、式(５４)、式(５５)より、式(５６)、式(５７)、式(５８)として与えられる。【０１６６】
【数５０】
2010134359.tif 000052
【０１６７】
なお、金型に樹脂が完全に充填されたときには、式(５９)が成立つとする。【０１６８】
【数５１】
2010134359.tif 000053
【０１６９】
次に式(６０)で定義される状態変数ｘ１、ｘ２、ｘ３を導入する。【０１７０】
【数５２】
2010134359.tif 000054
【０１７１】
式(６１)で定義される入力変数ｕを導入する。ｕは計測できるとする。高ゲイン観測器３１での計算では、モータ実電流ｉｍはモータ電流指令ｉ＊に等しいと考える。【０１７２】
【数５３】
2010134359.tif 000055
【０１７３】
計測できる状態変数としてｘ１を選び、出力変数ｙを式(６２)で定義する。【０１７４】
【数５４】
2010134359.tif 000056
【０１７５】
式(５６)、式(５７)、式(５８)、式(６２)を状態方程式で表現すると、式(６３)、式(６４)となる。【０１７６】
【数５５】
2010134359.tif 000057
【０１７７】
χ(ｘ2)、ψ(ｘ)は非線形関数を表す。式(６３）、式(６４)をベクトルｘを使って、式(６７)、式(６８)で表現する。【０１７８】
【数５６】
2010134359.tif 000058
【０１７９】
新しい状態変数Ｘ１、Ｘ２を式(６９)で定める。【０１８０】
【数５７】
2010134359.tif 000059
【０１８１】
式(６７)、式(６８)は、式(６９)より式(７０)、式(７１)で表せる。【０１８２】
【数５８】
2010134359.tif 000060
【０１８３】
式(７０)は、式(７３)で表せる。【０１８４】
【数５９】
2010134359.tif 000061
【０１８５】
状態変数Ｘ１＝ｘ１は計測できるので、状態変数Ｘ１は推定する必要がない。したがって、高ゲイン観測器３１は計測できるスクリュ位置ｘ１とモータ電流指令ｕを入力して、状態変数推定値Ｘ^２を出力する。その推定値Ｘ^２は式(７４)で与えられる(非特許文献２)。Kは高ゲイン観測器３１のゲイン定数で一般的に１より大きな定数が与えられる。【０１８６】
【数６０】
2010134359.tif 000062
【０１８７】
φ０(Ｘ^２，ｙ，u)は高ゲイン観測器３１で使われるφ(Ｘ，u)の公称（基準）関数である。式(７４)を式(７５)のように書く。【０１８８】
【数６１】
2010134359.tif 000063
【０１８９】
新変数ｗ^を式(７６)で導入する。【０１９０】
【数６２】
2010134359.tif 000064
【０１９１】
式(７５)に式(７６)を使うと、推定値Ｘ^２は式(７７)、式(７８)より得られる。【０１９２】
【数６３】
2010134359.tif 000065
【０１９３】
次に１より十分小さい正のパラメータεを導入して、ゲインKを式(７９)で与える。【０１９４】
【数６４】
2010134359.tif 000066
【０１９５】
式(７２)と式(７９)を利用すると、式(７７)は式(８０)で表される。【０１９６】
【数６５】
2010134359.tif 000067
【０１９７】
ここで、χ０(ｗ＾１、ｙ)、ψ０(ｗ＾、ｙ)は高ゲイン観測器で使われるχ(Ｘ２)、ψ(Ｘ)の公称（基準）関数である。【０１９８】
式（８１）で表される新変数η^１、η^２を導入する。【０１９９】
【数６６】
2010134359.tif 000068
【０２００】
式(８１)を利用して式(８０)を書き直すと、式(８２)が得られる。【０２０１】
【数６７】
2010134359.tif 000069
【０２０２】
式(８１)より、式(８３)が得られる。【０２０３】
【数６８】
2010134359.tif 000070
【０２０４】
式(７９)は式(８４)で表される。【０２０５】
【数６９】
2010134359.tif 000071
【０２０６】
式(８３)と式(８４)を使うと、式(７８)は式(８５)で表される。【０２０７】
【数７０】
2010134359.tif 000072
【０２０８】
以上より高ゲイン観測器３１での状態変数推定値ｘ^２、ｘ^３を得る。計算手順は式(８２)、式(８５)より、式(８６)、式(８７)で表される。
（１）   計算手順１【０２０９】
【数７１】
2010134359.tif 000073
【０２１０】
（２）   計算手順２【０２１１】
【数７２】
2010134359.tif 000074
【０２１２】
計算手順１でη^1、η^2を求め、計算手順２でｘ^2、ｘ^3を求める。【０２１３】
次に「課題を解決するための手段」で述べた高ゲイン観測器３１に要求される２つの必要条件
（A）時間的遅れが非常に小さい
（B）高精度である
が満たされることを明らかにする。式(８６)で公称関数χ０(η^１、ｙ)、ψ０(η^、ｙ)ではなく、実際には得られない真の関数χ(η１、ｙ)、ψ(η，ｙ)を使用したときに得られる変数η１、η２は式(８８)で決まる。【０２１４】
【数７３】
2010134359.tif 000075
【０２１５】
推定値誤差η〜１＝η１―η^１、η〜2 ＝η2 ―η^2は、式(８６)、式(８８)より式(８９)で得られる。【０２１６】
【数７４】
2010134359.tif 000076
【０２１７】
パラメータεは１より十分小さいことから、式(８９)より推定誤差η〜１、η〜2は制御対象のモデル誤差δ１、δ２から受ける影響を十分小さくできることが判る。すなわち高ゲイン観測器３１を用いれば、式(８６)、式(８７)より得られる射出圧力推定値ｘ＾３ 及びモータ回転速度（射出速度）推定値ｘ^２は前記必要条件(B)「高精度である」ことが満たされる。【０２１８】
次にモデル誤差δ１、δ２の影響を無視すると、式(８９)は式(９０)で表される。【０２１９】
【数７５】
2010134359.tif 000077
【０２２０】
前記行列A０の共役複素数の固有値λ１、λ２の実数部Ｒｅ(λ１)＝Ｒｅ(λ２)が負になるようにK０、すなわちK１、K２を決めると、式(９０)より推定誤差η〜１、η〜2は、その初期値η〜１０、η〜２０に対して式(９１)で与えられる。【０２２１】
【数７６】
2010134359.tif 000078
【０２２２】
ｔは時間変数で、C１(ｔ)〜C４(ｔ)は、行列A０の要素で決まる一定振幅、一定周波数の成分を表す。Ｒｅ(λ１)＜０で、パラメータεが１より十分小さいことから、式(９１)より推定誤差η〜１、η〜2は急速に０になることが判る。すなわち高ゲイン観測器３１を用いれば、式(８６)、式(８７)より得られる射出圧力推定値ｘ^３及びモータ回転速度（射出速度）推定値ｘ^２は前記必要条件(A)「時間的遅れが非常に小さい」ことが満たされる。【０２２３】
射出制御器２０は一定時間間隔Δｔ毎に制御演算を行うので、高ゲイン観測器３１での演算式(８６)、(８７)を離散演算式に変換する（非特許文献３、非特許文献４）。【０２２４】
新たなパラメータαを導入して、演算周期Δｔを式(９２)で表す。【０２２５】
【数７７】
2010134359.tif 000079
【０２２６】
式(８６)に前進矩形近似を適用すると、時間微分を表すラプラス演算子ｓとｚ変換演算子ｚの間には式(９３)が成立つ。【０２２７】
【数７８】
2010134359.tif 000080
【０２２８】
式(８６)に式(９３)を適用すると、式(９４)が成立つ。【０２２９】
【数７９】
2010134359.tif 000081
【０２３０】
式(９４)を離散演算式で表現すると、式(９５)が得られる。【０２３１】
【数８０】
2010134359.tif 000082
【０２３２】
ここで、離散時間ｔｋでの推定値η^１(ｔｋ)をη^１(ｋ)で表現している。ｙ(ｋ)、ｕ(ｋ)についても同じである。χ0(ｋ)は式(６５)より得られ、ψ0(ｋ)は式(６３)の3番目の式から得られ、時間ｔｋでの値を示す。式(９５)は式(９７)で表せる。【０２３３】
【数８１】
2010134359.tif 000083
【０２３４】
ここで、Ｉ２は、２行２列の単位行列を表す。【０２３５】
式(８７)の離散演算式は式(９８)で与えられる。【０２３６】
【数８２】
2010134359.tif 000084
【０２３７】
高ゲイン観測器３１は、演算式(９７)、(９８)を一定時間間隔Δｔ毎に演算することにより射出圧力推定値ｘ^３(ｋ)及びモータ回転速度（射出速度）推定値ｘ^２(ｋ)を得ることができる。式(９７)、(９８)で実現する高ゲイン観測器３１は、計測できる状態変数ｘ１(ｋ)は推定せず、必要となる状態変数である射出圧力ｘ３(ｋ)及びモータ回転速度(射出速度)ｘ２(ｋ)を推定する低次元高ゲイン観測器である。【０２３８】
高ゲイン観測器３１から射出速度推定値ｘ^２(ｋ)が得られることから図１には示さないが、射出制御器２０で射出速度の監視に使うことができる。またモータ制御器４０内の微分回路４７を廃止して減算器４８に速度信号ｖの代わりに射出速度推定値ｘ^２(ｋ)を入力することも可能である。【０２３９】
電動射出成形機を対象に高ゲイン観測器３１を使ったときのシミュレーション計算結果を示す。制御対象のモデル数値は次の通りである。
最大スクリュ位置  ｘｍａｘ＝３７．２ｃｍ
最大射出速度  ｖｍａｘ＝１３．２ｃｍ／ｓｅｃ
最大射出圧力  Ｐｍａｘ＝１７６５２Ｎ／ｃｍ２
モータ定格回転数  ωｍａｘ＝２０９．４ｒａｄ／ｓｅｃ（２０００ｒｐｍ）
対象のモデル定数を使って、式(６３)で使われる係数ａ，ｂ，ｃ，ｄは式(９９)の値を用いた。なお、本計算では抵抗成分χ(ｘ２)＝０とした。【０２４０】
【数８３】
2010134359.tif 000085
【０２４１】
関数ｆ(Ｐ^ｉ／Ｐｍａｘ)は適切な関数形を用いた。ここでは、任意の樹脂圧設定値Ｐｓｅｔ (ｔ)の波形を実現するために、Ｐ^ｉ(ｔ)／Ｐｓｅｔ(ｔ)の値によって金型への樹脂流入量を決める油圧機器の比例電磁式リリーフ弁の特性を用いた。式(７９)で与えられる高ゲイン観測器３１のゲインＫはＫ１＝０．３４０，Ｋ２＝−０．００３１６，ε＝０．０１５として式(１００)で与えられる。演算周期Δｔ＝５ｍｓｅｃを採用した。【０２４２】
【数８４】
2010134359.tif 000086
【０２４３】
図７(ａ)は、圧力検出センサ１０を用いたときの時間軸ｔｉｍｅ(ｓｅｃ)に対する実際の射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形、図７(ｂ)は、そのときの高ゲイン観測器３１が算出した射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]の波形、図７(ｃ)は、圧力検出センサ１０を用いずに図１に示すように射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]をフィードバック信号として減算器２７に入力したときの実際の射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形である。図７(ｂ)の高ゲイン観測器３１の算出した射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]の波形は、図7(ａ)の実際の射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形とほとんど一致していることが判る。更に図７(ａ)に示す圧力検出センサを用いたときの射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形と図７(ｃ)に示す圧力検出センサ１０を用いずに高ゲイン観測器３１を用いたときの射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形は、ほとんど同じであることが判る。スクリュ位置[ｘ／ｘｍａｘ]＝０．５５のときに充填工程から保圧工程への切換を行っており、図４での時間ｔ１＝２．６秒である。高ゲイン観測器３１は、射出圧力を時間的遅れなく、高精度に推定できることが判る。したがって、高ゲイン観測器３１が算出した射出圧力推定値は、圧力検出センサ信号の代わりに、充填工程での射出圧力監視及び保圧工程での射出圧力フィードバック信号として使えることが判る。【０２４４】
図８(ａ)は、射出速度[ｖ／ｖｍａｘ]の実際の制御波形で、図８(ｂ)は、高ゲイン観測器３１が算出した射出速度推定値[ｖ^／ｖｍａｘ]の波形である。高ゲイン観測器３１は、射出速度に対しても非常に精度良く、時間的遅れなく推定していることが判る。
【実施例２】
【０２４５】
図２は、本発明の一実施形態に係る高ゲイン観測器による射出圧力検知を電動射出成形機の制御装置に適用した例で、制御装置の内部構成を説明するためのブロック図である。この制御装置は、高ゲイン観測器３２を内蔵する射出制御器２０とモータ制御器４０（サーボアンプ）とから構成される。【０２４６】
射出制御器２０について説明する。射出制御器２０は一定時間間隔Δｔ毎に制御演算を行い、制御指令をモータ制御器４０に出力する。射出制御器２０は、図１に示す「実施例１」と同じ符号をもつ部品と高ゲイン観測器３２とから構成される。「実施例１」と同じ符号を持つ部品の機能説明は「実施例１」での説明に代える(段落(０１０２)〜(０１１０））。【０２４７】
高ゲイン観測器３２について説明する。モータ制御器４０内のモータ電流指令ｉ＊が、射出制御器２０内のＡ／Ｄ変換器３０を介して高ゲイン観測器３２に入力される。モータ制御器４０内のカウンタ４４の出力であるスクリュ９の位置信号ｘが高ゲイン観測器３２に入力される。更にモータ制御器４０内の微分回路４７の出力であるスクリュ９の速度信号ｖが高ゲイン観測器３２に入力される。高ゲイン観測器３２は、入力信号ｘ、ｖ及びｉ＊を使って、射出軸系の数式モデルを利用して導出した内蔵する離散演算式を実行して、スクリュ位置推定値ｘ^、射出速度推定値ｖ^及び射出圧力推定値Ｐ^ｉを出力する。【０２４８】
次にモータ制御器（サーボアンプ）４０について説明する。モータ制御器４０は、図１に示す「実施例１」と同じ符号をもつ部品から構成されているので、モータ制御器４０の部品の機能説明は「実施例１」での説明に代える(段落(０１１１)〜(０１１９））。【０２４９】
次に高ゲイン観測器３２の設計に必要となる射出軸系の数式モデルについて述べる。前記数式モデルを表す状態方程式表現は、「実施例１」と同じで、式(６７)と同じ式(１０１)で与えられる。【０２５０】
【数８５】
2010134359.tif 000087
【０２５１】
状態変数ｘ１、ｘ２、ｘ３及び入力変数ｕは式(６０)、式(６１)で与えられる「実施例１」と同じである。【０２５２】
計測できる状態変数としてスクリュ位置ｘ１及びモータ回転速度ｘ２を選ぶので、出力変数ｙは式(１０２)で定義される。【０２５３】
【数８６】
2010134359.tif 000088
【０２５４】
高ゲイン観測器３２は、計測できるｙ１＝ｘ１、ｙ２＝ｘ２及びモータ電流指令ｕを入力して、すべての状態変数の推定値ｘ^１、ｘ^２、ｘ^３を出力する。その推定値ｘ^１、ｘ^２、ｘ^３の演算式は式(１０３)で与えられる(非特許文献２)。【０２５５】
【数８７】
2010134359.tif 000089
【０２５６】
χ０(ｙ２)、ψ０(ｘ^)は高ゲイン観測器３２で使われるχ(ｙ２)、ψ(ｘ)の公称（基準）関数である。Ｋは高ゲイン観測器３２のゲイン定数で、１より十分小さい正のパラメータεを導入して式(１０４)で与えられる。【０２５７】
【数８８】
2010134359.tif 000090
【０２５８】
新しい推定値変数η^を式(１０５)で導入する。【０２５９】
【数８９】
2010134359.tif 000091
【０２６０】
式(１０３)を新変数η^で書き直すと、式(１０６)、式(１０７)を得る。【０２６１】
【数９０】
2010134359.tif 000092
【０２６２】
以上より、高ゲイン観測器３２での状態推定値ｘ^１、ｘ^２、ｘ^３を得る手順は、式(１０６)及び式(１０５)より式(１０８)、式(１０９)で与えられる。
（１）   計算手順１【０２６３】
【数９１】
2010134359.tif 000093
【０２６４】
（２）   計算手順２【０２６５】
【数９２】
2010134359.tif 000094
【０２６６】
すなわち、計算手順１でη^を求め、計算手順２でｘ^を求める。【０２６７】
次に「課題を解決するための手段」で述べた高ゲイン観測器３２に要求される２つの必要条件
（A）時間的遅れが非常に小さい
（B）高精度である
が満たされることを示す。式(１０６)で公称関数χ0(ｙ２)、ψ０(η^)ではなく、実際には得られない真の関数χ(ｙ２)、ψ(η)を使用したときに得られる変数η１、η２、η３は式(１１０)で決まる。【０２６８】
【数９３】
2010134359.tif 000095
【０２６９】
推定誤差η〜＝η−η^は、式(１０６)、式(１１０)より、式(１１１)で得られる。【０２７０】
【数９４】
2010134359.tif 000096
【０２７１】
パラメータεは１より十分小さいことから、式(１１１)より推定誤差η〜は制御対象のモデル誤差δ１ ,δ２から受ける影響を十分小さくできることが判る。すなわち高ゲイン観測器３２を用いれば、式(１０８)、式(１０９)より得られる射出圧力推定値ｘ^３、モータ回転速度（射出速度）推定値ｘ^２、及びスクリュ位置推定値ｘ^１は前記必要条件(B)「高精度である」ことが満たされる。【０２７２】
次にモデル誤差δ１ ,δ２の影響を無視すると、式(１１１)は式(１１２)で表される。【０２７３】
【数９５】
2010134359.tif 000097
【０２７４】
前記行列Ａ０の共役複素数の固有値λ１、λ２の実数部Ｒｅ(λ１)＝Ｒｅ(λ１)が負になるように、実数固有値λ３＜０になるように行列Ｋ０を決めると、式(１１２)より推定値誤差η〜１、η〜２、η〜３は、その初期値η〜１０、η〜２０、η〜３０に対して式(１１３)で与えられる。【０２７５】
【数９６】
2010134359.tif 000098
【０２７６】
ｔは時間関数で、ｇｉ(ｉ＝１〜６)関数は行列Ａ０の要素及び誤差初期値η〜１０、η〜２０、η〜３０で決まる有限の関数である。Ｒｅ(λ１)＜０、λ３＜０で、パラメータεが１より十分小さいことから、式(１１３)より推定誤差η〜１、η〜２、η〜３は急速に０になることが判る。すなわち高ゲイン観測器３２を用いれば、式(１０８)、式(１０９)より得られる射出圧力推定値ｘ^３、モータ回転速度（射出速度）推定値ｘ^２、及びスクリュ位置推定値ｘ^１は前記必要条件(A)「時間的遅れが非常に小さい」ことが満たされる。【０２７７】
射出制御器２０は一定時間間隔Δｔ毎に制御演算を行うので、高ゲイン観測器３２での演算式(１０８)、(１０９)を離散演算式に変換する(非特許文献３、非特許文献４)。【０２７８】
演算周期Δｔを式(１１４)で表す。【０２７９】
【数９７】
2010134359.tif 000099
【０２８０】
式(１０８)に前進矩形近似を適用すると、ラプラス演算子ｓとｚ変換演算子ｚの間には式(１１５)が成立つ。【０２８１】
【数９８】
2010134359.tif 000100
【０２８２】
式(１０８)に式(１１５)を適用すると、式(１１６)が成立つ。【０２８３】
【数９９】
2010134359.tif 000101
【０２８４】
式(１１６)を離散演算式で表現すると、式(１１７)が得られる。【０２８５】
【数１００】
2010134359.tif 000102
【０２８６】
ここで、離散時間ｔｋでの推定値η^(ｔｋ)をη^(ｋ)で表現している。χ０(ｋ)は式(６５)より得られ、ψ0(ｋ)は式(６３)の3番目の式から得られる。式(１１７)は式(１１９)で表せる。【０２８７】
【数１０１】
2010134359.tif 000103
【０２８８】
ここで、Ｉ３は３行３列の単位行列である。【０２８９】
式(１０９)の離散演算式は式(１２０)で与えられる。【０２９０】
【数１０２】
2010134359.tif 000104
【０２９１】
電動射出成形機を対象に高ゲイン観測器３２を使ったときのシミュレーション計算結果を示す。制御対象のモデル数値は「実施例１」のそれと同じである。式(１０１)で使われる係数ａ、ｂ、ｃ、ｄは式(９９)の「実施例１」のそれと同じである。なお、本計算では抵抗成分χ(ｘ２)＝０とした。【０２９２】
式(１０４)で与えられる高ゲイン観測器３２のゲインＫは、ε＝０．１として式(１２１)で与えられる。演算周期Δｔ＝５ｍｓｅｃとした。【０２９３】
【数１０３】
2010134359.tif 000105
【０２９４】
図９(ａ)は圧力検出センサ１０を用いたときの時間軸ｔｉｍｅ(ｓｅｃ)に対する実際の射出圧力[Ｐｉ ／Ｐｍａｘ]の制御波形、図９(ｂ)は、そのときの高ゲイン観測器３２が算出した射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]の波形、図９(ｃ)は、圧力検出センサ１０を用いずに図２に示すように射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]をフィードバック信号として減算器２７に入力したときの実際の射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形である。図９(ｂ)の高ゲイン観測器３２の算出した射出圧力推定値[Ｐ^ｉ／Ｐｍａｘ]の波形は、図９(ａ)の実際の射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形とほとんど一致していることが判る。更に図９(ａ)に示す圧力検出センサ１０を用いたときの射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形と図９(ｃ)に示す圧力検出センサ１０を用いずに高ゲイン観測器３２を用いたときの射出圧力[Ｐｉ／Ｐｍａｘ]の制御波形は、ほとんど同じであることが判る。高ゲイン観測器３２は、射出圧力を時間的遅れなく、高精度に推定できることが判る。したがって、高ゲイン観測器３２が算出した射出圧力推定値は、圧力検出センサ信号の代わりに、充填工程での射出圧力監視及び保圧工程での射出圧力フィードバック信号として使えることが判る。図示はしないが、スクリュ位置及び射出速度の推定値も、射出圧力推定値と同様に非常に良好な結果が得られている。【０２９５】
「実施例１」及び「実施例２」での高ゲイン観測器は、いずれもモータ多軸駆動の電動射出成形機でも、圧力検出センサ信号の代わりに高ゲイン観測器の出力である射出圧力推定値を採用する圧力制御装置及び圧力制御方法として適用できる。
【産業上の利用可能性】
【０２９６】
電動射出成形機の圧力制御装置及び圧力制御方法において、圧力検出センサを使わずに、高ゲイン観測器が出力する射出圧力推定値を射出圧力検知信号として使うことにより、次の4つの不利を回避できる。
（１）   高圧環境下で信頼性の高い圧力センサは高価になる。
（２）   金型内やスクリュノズル部への圧力センサ取付けは、特別な加工を施す必要があり、作業コストが無視できない。
（３）   電動モータから射出スクリュに至る射出軸系に取り付けるロードセルは、組み込むための機械構造を複雑にし、更には射出軸系の機械剛性の低下を招く。
（４）   歪みゲージを検出部に使用するロードセルでは、微弱なアナログ信号に対するノイズ対策が必要になり、また信号アンプのゼロ点調整やスパン調整等にも人手による作業が必要になる（特許文献１０）。【０２９７】
更に、高ゲイン観測器の出力する射出圧力推定値は、高精度であり、時間的遅れも非常に小さいので、射出圧力の監視信号及び射出圧力制御の射出圧力フィードバック信号として使える。したがって、本発明による高ゲイン観測器による電動射出成形機の圧力制御装置及び圧力制御方法は、十分利用される価値を有すると考えられる。
【符号の説明】
【０２９８】
１      金型
２      射出シリンダ
３      サーボモータ
４      減速機
５      ボールネジ
６      軸受
７      ナット
８      可動部
９      スクリュ
１０    圧力検出センサ
１１    リニアガイド
１２    モータエンコーダ
１３    キャビティ
２０    射出制御器
２１    射出速度指令器
２２    変換器
２３    パルス発生器
２４    パルス列
２５    アナログ／デジタル（Ａ／Ｄ）変換器
２６    射出圧力指令器
２７    減算器
２８    圧力制御器
２９    デジタル／アナログ（Ｄ／Ａ）変換器
３０    アナログ／デジタル（Ａ／Ｄ）変換器
３１    高ゲイン観測器
３２    高ゲイン観測器                                       
４０    モータ制御器（サーボアンプ）
４１    カウンタ
４２    アナログ／デジタル（Ａ／Ｄ）変換器
４３    比較器
４４    カウンタ
４５    減算器
４６    位置制御器
４７    微分回路
４８    減算器
４９    速度制御器
５０    パルス幅変調制御（ＰＷＭ）回路
【図１】
2010134359.tif 000112
【図２】
2010134359.tif 000113
【図３】
2010134359.tif 000114
【図４】
2010134359.tif 000115
【図５】
2010134359.tif 000116
【図６】
2010134359.tif 000117
【図７】
2010134359.tif 000118
【図８】
2010134359.tif 000119
【図９】
2010134359.tif 000120
【手続補正書】
【提出日】平成22年2月16日(2010.2.16)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項２
【補正方法】変更
【補正の内容】
【請求項２】
電動射出成形機の圧力制御装置であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルを利用して導出した離散演算式を実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或いはモータ実電流信号を入力信号として下記(数１０４)の状態方程式で表現される前記射出軸系の数式モデルを利用して導出した下記(数１０５)の離散式(１２５)および(数１０６)の離散式(１２８)からなる前記離散演算式を実行して、射出圧力推定値ｘ^３および射出速度推定値ｘ^２を出力し、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値とを減算器に入力し、前記減算器は前記射出圧力指令と前記射出圧力推定値との差を算出して圧力制御器に入力し、前記圧力制御器は、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出してモータ制御器に出力し、前記モータ制御器はサーボモータに所定のモータトルクを発生させて射出圧力を制御する電動射出成形機の圧力制御方法
【数１０４】
2010134359.tif 000121
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数
ｘ２：射出速度を最大射出速度で無次元化した状態変数
ｘ３：射出圧力を最大射出圧力で無次元化した状態変数
ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数
ｙ：計測できる状態変数ｘ１を表す出力変数
ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ：射出・保圧機構のモデル定数
ｆ(ｘ３)：射出圧力で決まる金型への樹脂流入量を決める関数
χ(ｘ２)、ψ(ｘ)：式(１２４)で表される非線形関数
【数１０５】
2010134359.tif 000122
ここで、η^１、η^２：状態変数ｘ２、ｘ３を推定するために導入した新状態変数
Ｉ２：２行２列の単位行列
ｋ：時間ｔｋを表す離散変数（ｋ＝０、１、２、・・・）
Δｔ：高ゲイン観測器の演算周期
ε：高ゲイン観測器でのパラメータで、一般に1より十分小さい正数
K１、Ｋ２：高ゲイン観測器のゲインを決めるパラメータで、行列A０の固有値の実数部が負になるように決める。
【数１０６】
2010134359.tif 000123
ここで、ｘ^２、ｘ^３：状態変数ｘ２、ｘ３の推定値
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項３
【補正方法】変更
【補正の内容】
【請求項３】
電動射出成形機の圧力制御装置であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルを利用して導出した離散演算式を実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と射出速度信号（モータ回転速度）と前記サーボモータへの制御信号として印加されるモータ電流指令信号或いはモータ実電流信号を入力信号として下記(数１０７)の状態方程式で表現される射出軸系の数式モデルを利用して導出した下記(数１０８)の離散式(１３２)および(数１０９)の離散式(１３５)からなる前記離散演算式を実行して、射出圧力推定値ｘ^３、射出速度推定値ｘ^２およびスクリュ位置推定値ｘ^１を出力し、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値とを減算器に入力し、前記減算器は前記射出圧力指令と前記射出圧力推定値との差を算出して圧力制御器に入力し、前記圧力制御器は、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出してモータ制御器に出力し、前記モータ制御器はサーボモータに所定のモータトルクを発生させて射出圧力を制御する電動射出成形機の圧力制御方法
【数１０７】
2010134359.tif 000124
ここで、ｙ１、ｙ２：計測できる状態変数ｘ１及びｘ２をそれぞれ表す出力変数
その他の変数の意味は、(数１０４)と同じである。
【数１０８】
2010134359.tif 000125
ここで、η^１、η^２、η^３：状態変数ｘ１、ｘ２、ｘ３を推定するために導入した新状態変数
Ｉ３：３行３列の単位行列
ｋ：時間ｔｋを表す離散変数(ｋ＝０、１、２、・・・)
Δｔ：高ゲイン観測器の演算周期
ε：高ゲイン観測器でのパラメータで、一般に1より十分小さい正数
K０：高ゲイン観測器のゲインを決める行列で、行列A０の固有値の実数部が負になるように行列要素を決める。
【数１０９】
2010134359.tif 000126
ここで、ｘ^1、ｘ^２、ｘ^３：状態変数ｘ１、ｘ２、ｘ３の推定値
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０１３８
【補正方法】変更
【補正の内容】
【０１３８】
【数３７】
2010134359.tif 000127
【手続補正書】
【提出日】平成22年5月27日(2010.5.27)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１
【補正方法】変更
【補正の内容】
【請求項１】
電動射出成形機の圧力制御装置であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構において、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号が入力されて、射出圧力推定値を出力する高ゲイン観測器と、射出圧力を指令するための射出圧力指令を出力する射出圧力指令器と、前記射出圧力指令と前記射出圧力推定値とが入力されて、前記射出圧力指令と前記射出圧力推定値との差を出力する減算器と、前記減算器の出力を入力して、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出する圧力制御器と、を含む射出制御器と、
前記モータ電流指令信号が入力されるモータ制御器と、
を具備することを特徴とする電動射出成形機の圧力制御装置。
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項２
【補正方法】変更
【補正の内容】
【請求項２】
電動射出成形機の圧力制御方法であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルを利用して導出した離散演算式を実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号を入力信号として下記(数１０４)の状態方程式で表現される前記射出軸系の数式モデルを利用して導出した下記（数１０５）の離散式（１２５）および（数１０６）の離散式（１２８）からなる前記離散演算式を実行して、射出圧力推定値ｘ＾３および射出速度推定値ｘ＾２を出力し、減算器が、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値を入力信号として前記射出圧力指令と前記射出圧力推定値との差を出力し、圧力制御器が、前記減算器の出力を入力信号として前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を出力し、モータ制御器が、前記モータ電流指令信号を入力信号として、前記サーボモータに前記モータ電流指令信号に相当するモータトルク（モータ実電流）を発生させて前記射出圧力指令に等しい射出圧力を実現する電動射出成形機の圧力制御方法。
【数１０４】
2010134359.tif 000128
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数、  ｘ2:射出速度を最大射出速度で無次元化した状態変数、  ｘ3:射出圧力を最大射出圧力で無次元化した状態変数、  ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数、  ｙ：計測できる状態変数ｘ１を表す出力変数、  ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ：射出軸系の数式モデルの定数、  ｆ（ｘ３）：射出圧力で決まる金型への樹脂流入量を決める関数、χ（ｘ２）、ψ（ｘ）：式（１２４）で表される非線形関数。
【数１０５】
2010134359.tif 000129
ここで、ｋ：離散時間ｔｋを表す離散変数（ｋ＝０，１，２、・・・）、  η＾１（ｋ）、η＾２（ｋ）：状態変数ｘ２、ｘ３を推定するために導入した新状態変数η１、η２の離散時間ｔｋでの推定値η＾１（ｔｋ）、η＾２（ｔｋ）、  ｙ（ｋ）、ｕ（ｋ）：離散時間ｔｋでの出力変数値ｙ（ｔｋ）及び入力変数値ｕ（ｔｋ）、  ｘ＾２（ｋ）、ｘ＾３（ｋ）：離散時間ｔｋでの状態推定値ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）、  χ（ｋ）、ψ（ｋ）：離散時間ｔｋでの非線形関数値χ（ｔｋ）、ψ（ｔｋ）、I２：２行２列の単位行列、Δｔ：高ゲイン観測器の演算周期、ε：高ゲイン観測器でのパラメータで、一般に１より十分小さい正数、Ｋ１、Ｋ２：高ゲイン観測器のゲインを決めるパラメータで、行列Ａ０の固有値の実数部が負になるように決める。
【数１０６】
2010134359.tif 000130
ここで、ｘ＾２（ｋ）、ｘ＾３（ｋ）：離間時間ｔｋでの状態変数ｘ２、ｘ３の状態推定値ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）。
【手続補正３】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項３
【補正方法】変更
【補正の内容】
【請求項３】
電動射出成形機の圧力制御方法であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルを利用して導出した離散演算式を実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と射出速度信号（モータ回転速度信号）と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号を入力信号として下記(数１０７)の状態方程式で表現される前記射出軸系の数式モデルを利用して導出した下記（数１０８）の離散式（１３２）および（数１０９）の離散式（１３５）からなる前記離散演算式を実行して、射出圧力推定値ｘ＾３、射出速度推定値ｘ＾２およびスクリュ位置推定値ｘ＾１を出力し、減算器が、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値を入力信号として前記射出圧力指令と前記射出圧力推定値との差を出力し、圧力制御器が、前記減算器の出力を入力信号として前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を出力し、モータ制御器が、前記モータ電流指令信号を入力信号として、前記サーボモータに前記モータ電流指令信号に相当するモータトルク（モータ実電流）を発生させて前記射出圧力指令に等しい射出圧力を実現する電動射出成形機の圧力制御方法。
【数１０７】
2010134359.tif 000131
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数、  ｘ2:射出速度を最大射出速度で無次元化した状態変数、  ｘ3:射出圧力を最大射出圧力で無次元化した状態変数、  ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数、  ｙ１、ｙ２：計測できる状態変数ｘ１及びｘ２をそれぞれ表す出力変数、  ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ：射出軸系の数式モデルの定数、  ｆ（ｘ３）：射出圧力で決まる金型への樹脂流入量を決める関数、χ（ｘ２）、ψ（ｘ）：式（１３１）で表される非線形関数。
【数１０８】
2010134359.tif 000132
ここで、ｋ：離散時間ｔｋを表す離散変数（ｋ＝０、１、２、・・・）、  η＾１（ｋ）、η＾２（ｋ）、η＾３（ｋ）：状態変数ｘ１、ｘ２、ｘ３を推定するために導入した新状態変数η１、η２、η３の離散時間ｔｋでの推定値η＾１（ｔｋ）、η＾２（ｔｋ）、η＾３（ｔｋ）、  ｙ１（ｋ）、ｙ２（ｋ）、ｕ（ｋ）：離散時間ｔｋでの出力変数値ｙ１（ｔｋ）、ｙ２（ｔｋ）及び入力変数値ｕ（ｔｋ）、  ｘ＾３（ｋ）：離散時間ｔｋでの状態推定値ｘ＾３（ｔｋ）、  χ（ｋ）、ψ（ｋ）：離散時間ｔｋでの非線形関数値χ（ｔｋ），ψ（ｔｋ）、I３：３行３列の単位行列、Δｔ：高ゲイン観測器の演算周期、ε：高ゲイン観測器でのパラメータで、一般に１より十分小さい正数、Ｋ０：高ゲイン観測器のゲインを決める行列で、行列Ａ０の固有値の実数部が負になるように行列要素を決める。
【数１０９】
2010134359.tif 000133
ここで、ｘ＾１（ｋ）、ｘ＾２（ｋ）、ｘ＾３（ｋ）：離間時間ｔｋでの状態変数ｘ１、ｘ２、ｘ３の状態推定値ｘ＾１（ｔｋ）、ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）。
【手続補正書】
【提出日】平成22年8月2日(2010.8.2)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００５２
【補正方法】変更
【補正の内容】
【００５２】
次に検知圧力の誤差は、成形不良や安全運転の欠如につながることから、高精度の圧力検知が要求される。したがって、前記２つの必要条件
（A）    
時間的遅れが非常に小さい
（B）    
高精度である
を満たす圧力検知を実現するためには、高ゲイン観測器（非特許文献１）の手法を採用する。計測できる変数を入力して、すべての状態変数を推定する高ゲイン観測器が、前記２つの必要条件（A）（B）を満たしていることを簡単な制御対象モデルを使って説明する。式（５）は制御対象モデルの状態方程式および出力方程式を示す。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０１３２
【補正方法】変更
【補正の内容】
【０１３２】
【数３４】
2010134359.tif 000134
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０１３３
【補正方法】変更
【補正の内容】
【０１３３】
ここで、Aｓはスクリュ断面積、Pｉは射出シリンダ貯留部（スクリュ先端部）の樹脂圧力を意味する射出圧力、Cｍｔは射出シリンダ粘性係数、γは速度べき乗数である。射出圧力の方程式は式（３９）、式（４０）で与えられる。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】０１４４
【補正方法】変更
【補正の内容】
【０１４４】
【数４０】
2010134359.tif 000135
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】０１５３
【補正方法】変更
【補正の内容】
【０１５３】
【数４４】
2010134359.tif 000136
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】０１５７
【補正方法】変更
【補正の内容】
【０１５７】
【数４６】
2010134359.tif 000137
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】０１６６
【補正方法】変更
【補正の内容】
【０１６６】
【数５０】
2010134359.tif 000138
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】０１７３
【補正方法】変更
【補正の内容】
【０１７３】
計測できる状態変数としてｘ1を選び、出力変数ｙを式（６２）の出力方程式で定義する。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】０１７５
【補正方法】変更
【補正の内容】
【０１７５】
式（５６）、式（５７）、式（５８）、式（６２）を状態方程式と出力方程式で表現すると、式（６３）、式（６４）となる。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】０１７６
【補正方法】変更
【補正の内容】
【０１７６】
【数５５】
2010134359.tif 000139
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】０２３１
【補正方法】変更
【補正の内容】
【０２３１】
【数８０】
2010134359.tif 000140
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】０２５２
【補正方法】変更
【補正の内容】
【０２５２】
計測できる状態変数としてスクリュ位置ｘ1及びモータ回転速度ｘ２を選ぶので、出力変数ｙは式（１０２）の出力方程式で定義される。
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】０２８５
【補正方法】変更
【補正の内容】
【０２８５】
【数１００】
2010134359.tif 000141
【手続補正１４】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１
【補正方法】変更
【補正の内容】
【請求項１】
電動射出成形機の圧力制御装置であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルとして、前記スクリュの位置変数、射出速度変数および射出圧力変数の３変数からなる状態変数と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号を入力変数とする状態方程式と計測できる状態変数を出力変数とする出力方程式からなる連続時間系の数式モデルを採用し、前記数式モデルに対して前進矩形近似を適用して導出した離散演算式を一定時間間隔毎に実行する高ゲイン観測器と、射出圧力を指令するための射出圧力指令を出力する射出圧力指令器と、前記高ゲイン観測器が前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記モータ電流指令信号或は前記モータ実電流信号が入力されて内蔵する前記離散演算式を使用して算出して出力する射出圧力推定値と前記射出圧力指令器が出力する前記射出圧力指令とが入力されて、前記射出圧力指令と前記射出圧力推定値との差を出力する減算器と、前記減算器の出力を入力して、前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を算出する圧力制御器と、を含む射出制御器と、
前記モータ電流指令信号が入力されるモータ制御器と、
を具備することを特徴とする電動射出成形機の圧力制御装置
【手続補正１５】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項２
【補正方法】変更
【補正の内容】
【請求項２】
電動射出成形機の圧力制御方法であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルとして採用した下記（数１０４）の状態方程式（１２２）および出力方程式（１２３）に前進矩形近似を適用して導出した下記（数１０５）の離散演算式（１２５）および（数１０６）の離散演算式（１２８）を一定時間間隔毎に実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号を入力信号として射出圧力推定値ｘ＾３及び射出速度推定値ｘ＾２を出力し、減算器が、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値を入力信号として前記射出圧力指令と前記射出圧力推定値との差を出力し、圧力制御器が、前記減算器の出力を入力信号として前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を出力し、モータ制御器が、前記モータ電流指令信号を入力信号として、前記サーボモータに前記モータ電流指令信号に相当するモータトルクを発生させて前記射出圧力指令に等しい射出圧力を実現する電動射出成形機の圧力制御方法
【数１０４】
2010134359.tif 000142
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数、  ｘ2:射出速度を最大射出速度で無次元化した状態変数、  ｘ3:射出圧力を最大射出圧力で無次元化した状態変数、  ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数、  ｙ：計測できる状態変数ｘ１を表す出力変数、  ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ、γ：射出軸系の数式モデルの定数、  ｆ（ｘ３）：射出圧力で決まる金型への樹脂流入量を決める関数、χ（ｘ２）、ψ（ｘ）：式（１２４）で表される非線形関数
【数１０５】
2010134359.tif 000143
ここで、ｋ：離散時間ｔｋを表す離散変数（ｋ＝０，１，２、・・・）、  η＾１（ｋ）、η＾２（ｋ）：状態変数ｘ２、ｘ３を推定するために導入した新状態変数η１、η２の離散時間ｔｋでの推定値η＾１（ｔｋ）、η＾２（ｔｋ）、  ｙ（ｋ）、ｕ（ｋ）：離散時間ｔｋでの出力変数値ｙ（ｔｋ）及び入力変数値ｕ（ｔｋ）、  ｘ＾２（ｋ）、ｘ＾３（ｋ）：離散時間ｔｋでの状態推定値ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）、  χ（ｋ）、ψ（ｋ）：離散時間ｔｋでの非線形関数値χ（ｔｋ）、ψ（ｔｋ）、  I２：２行２列の単位行列、  Δｔ：高ゲイン観測器の演算周期、  ε：高ゲイン観測器でのパラメータで、一般に１より十分小さい正数、  Ｋ１、Ｋ２：高ゲイン観測器のゲインを決めるパラメータで、行列Ａ０の固有値の実数部が負になるように決める。
【数１０６】
2010134359.tif 000144
ここで、ｘ＾２（ｋ）、ｘ＾３（ｋ）：離間時間ｔｋでの状態変数ｘ２、ｘ３の状態推定値ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）
【手続補正１６】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項３
【補正方法】変更
【補正の内容】
【請求項３】
電動射出成形機の圧力制御方法であって、サーボモータの回転は減速機を介してボールネジに伝えられ、前記ボールネジの回転はボールネジ軸上のナットの直線運動に変換され、前記ナットにより駆動される可動部を介してスクリュが前進移動し、前記スクリュの前進により射出シリンダ先端部に貯留された溶融樹脂への加圧と金型内キャビティへの樹脂充填を実現する射出・保圧機構の運動を表現する射出軸系の数式モデルとして採用した下記（数１０７）の状態方程式（１２９）および出力方程式（１３０）に前進矩形近似を適用して導出した下記（数１０８）の離散演算式（１３２）および（数１０９）の離散演算式（１３５）を一定時間間隔毎に実行する高ゲイン観測器が、前記サーボモータの軸に設けたモータエンコーダで検出されるスクリュ位置信号と射出速度信号と前記サーボモータへの制御信号として印加されるモータ電流指令信号或はモータ実電流信号を入力信号として射出圧力推定値ｘ＾3、射出速度推定値ｘ＾２およびスクリュ位置推定値ｘ＾１を出力し、減算器が、射出圧力指令器が出力する射出圧力指令と前記射出圧力推定値を入力信号として前記射出圧力指令と前記射出圧力推定値との差を出力し、圧力制御器が、前記減算器の出力を入力信号として前記射出圧力推定値が前記射出圧力指令に追従するように前記モータ電流指令信号を出力し、モータ制御器が、前記モータ電流指令信号を入力信号として、前記サーボモータに前記モータ電流指令信号に相当するモータトルクを発生させて前記射出圧力指令に等しい射出圧力を実現する電動射出成形機の圧力制御方法
【数１０７】
2010134359.tif 000145
ここで、ｘ１：スクリュ位置を最大スクリュ位置で無次元化した状態変数、  ｘ2:射出速度を最大射出速度で無次元化した状態変数、  ｘ3:射出圧力を最大射出圧力で無次元化した状態変数、  ｕ：モータ電流指令或はモータ実電流をモータ定格電流で無次元化した入力変数、  ｙ１、ｙ２：計測できる状態変数ｘ１及びｘ２をそれぞれ表す出力変数、  ａ、ｂ、ｃ、ｄ、ｅ、ｈ、ｐ、γ：射出軸系の数式モデルの定数、  ｆ（ｘ３）：射出圧力で決まる金型への樹脂流入量を決める関数、  χ（ｘ２）、ψ（ｘ）：式（１３１）で表される非線形関数
【数１０８】
2010134359.tif 000146
ここで、ｋ：離散時間ｔｋを表す離散変数（ｋ＝０、１、２、・・・）、  η＾１（ｋ）、η＾２（ｋ）、η＾３（ｋ）：状態変数ｘ１、ｘ２、ｘ３を推定するために導入した新状態変数η１、η２、η３の離散時間ｔｋでの推定値η＾１（ｔｋ）、η＾２（ｔｋ）、η＾３（ｔｋ）、  ｙ１（ｋ）、ｙ２（ｋ）、ｕ（ｋ）：離散時間ｔｋでの出力変数値ｙ１（ｔｋ）、ｙ２（ｔｋ）及び入力変数値ｕ（ｔｋ）、  ｘ＾３（ｋ）：離散時間ｔｋでの状態推定値ｘ＾３（ｔｋ）、  χ（ｋ）、ψ（ｋ）：離散時間ｔｋでの非線形関数値χ（ｔｋ）、ψ（ｔｋ）、  I３：３行３列の単位行列、  Δｔ：高ゲイン観測器の演算周期、  ε：高ゲイン観測器でのパラメータで、一般に１より十分小さい正数、  Ｋ０：高ゲイン観測器のゲインを決める行列で、行列Ａ０の固有値の実数部が負になるように行列要素を決める。
【数１０９】
2010134359.tif 000147
ここで、ｘ＾１（ｋ）、ｘ＾２（ｋ）、ｘ＾３（ｋ）：離間時間ｔｋでの状態変数ｘ１、ｘ２、ｘ３の状態推定値ｘ＾１（ｔｋ）、ｘ＾２（ｔｋ）、ｘ＾３（ｔｋ）
【国際調査報告】
2010134359.tif 000148
2010134359.tif 000149
2010134359.tif 000150
2010134359.tif 000151