調整された横方向および縦方向制御のためのモデル予測制御装置
１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスが開示される。本プロセスは、シート測定データを収集するために規則的な間隔で複数のシート特性を測定することを含む。シート測定データは、複数のシート特性測定値配列を確立するために処理され、その配列は、次いで、共通解像度にマップされる。共通解像度シート特性測定値配列は、１つのより大きい１次元共通解像度測定値配列へと結合される。共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列は、抄紙機プロセス・モデル状態オブザーバへの入力として使用されて、シート製造プロセスの推定される現在の内部状態が生成される。複数の未来のシート特性目標配列は結合されて、１つの目標配列となる。ウェブ製造プロセスの推定される現在の内部状態配列、および抄紙機プロセス・モデルは、シート特性の未来の予測配列を生成するために用いられる。シート特性の未来の予測配列と、目標配列と、オブジェクト関数の重みと、最後のアクチュエータ設定値と、ハード制約とが、オブジェクト関数に入力され、その関数を解くことによりシート製作プロセスの調整されたＭＤおよびＣＤ制御のためのアクチュエータ設定値の最適な変更を得ることができる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シート製作プロセスの制御に関連し、より詳細には、シート製作機における縦方向および横方向アクチュエータのオペレーションを調整するのための方法に関連する。
【背景技術】
【０００２】
シート製作機におけるシート特性の制御は、シートの特性を可能な限り目標値に近い値に保つことに関する。シート特性の制御には、２組の異なるアクチュエータが使用される。第１に、シート特性の横方向（ＣＤ：ｃｒｏｓｓ ｄｉｒｅｃｔｉｏｎ）平均値にだけ影響する縦方向（ＭＤ：ｍａｃｈｉｎｅ ｄｉｒｅｃｔｉｏｎ）アクチュエータがある。各ＭＤアクチュエータは、シート特性において異なる動的応答を有することができる。第２に、ＣＤ方向にシートを横断して配列されたＣＤアクチュエータがある。ＣＤアクチュエータの各アレイは、シート特性の平均値及びＣＤ形状の双方に影響することができる。ＣＤアクチュエータは、シート特性における異なる動的応答および異なる空間応答を有することができる。シート特性の全体的な制御問題は、非常に多変量になる。即ち、ＣＤアレイ中の１つのＣＤアクチュエータは、幾つかのシート特性において隣接するＣＤゾーンに影響し、また特定のシート特性を制御するよう意図したＣＤアクチュエータの平均的な効果は、幾つかのＭＤアクチュエータによっても影響される幾つかのシート特性の平均値に影響する可能性がある。その問題はまた、非常に大規模なものである。典型的な制御プロセスは、数千の出力（シート特性の測定値）、および数百の入力（アクチュエータの設定値）を有する可能性がある。プロセスはまた、特定の空間方向およびアクチュエータ内設定方向に制御することは困難または不可能である。
【０００３】
現在、従来の大部分のシート製作装置では、シート特性の制御を２つの制御問題に分離する。第１に、ＣＤ平均値は、シート特性のＣＤ平均値に対するＣＤアクチュエータの効果を利用せずに、ＭＤアクチュエータを利用するだけで制御される。第２に、シートを横断して配列されたＣＤアクチュエータは、シート特性の平均値まわりのＣＤ変動を制御するためだけに利用される。現在、ＭＤアクチュエータを調整するための明示的なハード制約処理とともに、モデル予測制御を用いるＭＤ制御方式が利用できる。
【０００４】
モデル予測制御を用いて１つまたは複数のシート特性を制御する、ＣＤアクチュエータ・アレイの最適に調整された制御は、例えば、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｈｅｎｄｅｒｓｏｎ Ｂ、およびＳｔｅｗａｒｔ Ｃ、「Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ａｎｄ ｍｕｌｔｉｖａｒｉａｂｌｅ ｃｏｎｔｒｏｌ ｏｆ ｓｕｐｅｒｃａｌｅｎｄｅｒｓ」、Ｃｏｎｔｒｏｌ Ｓｙｓｔｅｍｓ ２００２、２００２年６月、ストックホルム、スウェーデン、および、Ｂａｃｋｓｔｒｏｍ Ｊ． Ｕ、Ｇｈｅｏｒｇｈｅ Ｃ、Ｓｔｅｗａｒｔ Ｇ． Ｅ、Ｖｙｓｅ Ｒ．Ｎ、「Ｃｏｎｓｔｒａｉｎｅｄ ｍｏｄｅｌ ｐｒｅｄｉｃｔｉｖｅ ｃｏｎｔｒｏｌ ｆｏｒ ｃｒｏｓｓ ｄｉｒｅｃｔｉｏｎａｌ ｍｕｌｔｉ−ａｒｒａｙ ｐｒｏｃｅｓｓｅｓ」、Ｐｕｌｐ ＆ Ｐａｐｅｒ Ｃａｎａｄａ、Ｔ１２８ １０２：５ （２００１） などの記事で論じられている。
【０００５】
ＭＤアクチュエータおよびＣＤアクチュエータを調整する必要性は、本願の権利者が所有する２０００年７月２５日に発行された米国特許第６，０９４，６０４号で明らかにされている。この問題に対して提案される解決策はまた、互いに通信するアクチュエータにおける分散型局所化インテリジェント制御装置のシステムを含む‘６０４特許中に開示されている。
【発明の開示】
【課題を解決するための手段】
【０００６】
（発明の概要）
上記で概要を説明した問題に対処するために、本発明は、調整されたＭＤおよびＣＤ制御のための柔軟な大規模多変量モデル予測制御装置を提供し、その制御装置は、入力としてシート特性測定値の複数の配列を取り込み、複数の出力配列（アクチュエータ設定値）を生成する。その配列は、任意の次元のものとすることができる。１つのＭＤ配列は、１×１の配列と見なすことができる。任意の数の入力および出力配列とすることができる。本発明は、新しく最適に調整される設定値を均等な制御間隔で計算する。各シート特性に対して、ＣＤコンポーネントだけ、ＭＤコンポーネントだけ、またはＭＤ及びＣＤコンポーネントの双方に、制御することができる。本発明は、予測区間（ｐｒｅｄｉｃｔｉｏｎ ｈｏｒｉｚｏｎ）Ｈｐから未来の制御区間Ｈｃにおけるアクチュエータ設定値に対する、動的及び２次元応答を予測する。次いで、本発明は、未来の予測シート特性を可能な限り目標に近付ける未来の最適な設定値を計算する。制御装置はまた、アクチュエータに対する物理的制限を明示的に考慮に入れる。制御装置は、アクチュエータ設定値を困難な空間方向およびアクチュエータ内設定プロセス方向で発行することを回避することにより、２つのタイプの方向問題を処理する。これにより、閉ループの２次元ロバストな安定性を得ることができる。
【０００７】
従って、本発明は、１枚の材料シートを製作するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスを提供し、そのプロセスは、
シート測定データを収集するために規則的な間隔で複数のシート特性を測定するステップと、
そのシート測定データを、複数のシート特性測定値配列を確立するために処理するステップと、
１次元共通解像度測定値配列を確立するためにシートの特性測定値配列を処理するステップと、
シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
未来のシート特性の目標配列を確立するステップと、
シート製造プロセスの推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
シート特性の未来の予測配列、未来のシート特性目標配列、および前のアクチュエータ設定値配列を、可解なオブジェクト関数に入力して、シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含む。
【０００８】
本発明はまた、１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスを提供し、そのプロセスは、
シート測定データを収集するために規則的な間隔で複数のシート特性を測定するステップと、
複数のシート特性測定値配列を確立するために、シート測定データを処理するステップと、
シート特性測定値配列を共通解像度にマップするステップと、
共通解像度シート特性測定値配列を、１つのより大きい１次元共通解像度測定値配列に結合するステップと、
共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力することにより、シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
複数の未来のシート特性目標配列を１つの目標配列に結合するステップと、
シート製造プロセスの推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
シート特性の未来の予測配列、目標配列、オブジェクト関数の重み、最後のアクチュエータ設定値配列、およびハード制約を、オブジェクト関数に入力するステップと、
オブジェクト関数を解いて、シート製作プロセスの調整されたＭＤおよびＣＤ制御のために現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含む。
【０００９】
本発明は、ＣＤアクチュエータ・アレイおよびＭＤアクチュエータを最適に操作しかつ調整するように動作し、シート特性におけるＭＤおよびＣＤ変動を最小化することができる。
【００１０】
本発明は、ＭＤアクチュエータとＣＤアクチュエータ・アレイの間の対話を最適に調整する。本発明は更に、空間利得の低い方向に移動するコストを表現するためにオブジェクティト関数中に一般（ｇｅｎｅｒａｌ）重み付け関数を有する。本発明は更に、小さなアクチュエータ内設定方向に移動するコストを表現するための明示的な重み付け関数を有する。本発明は更に、ＣＤアクチュエータ・アレイ設定値平均の許容可能範囲を指定するハード制約を含む。本発明は、シート特性のＣＤコンポーネントだけ、ＭＤコンポーネントだけ、あるいはＣＤ及びＭＤコンポーネントの双方に制御するようにセットアップすることができる。
【００１１】
本発明のプロセスは、複数の分散された制御装置ではなく、１つの集中化された制御装置を使用することが好ましい。
【００１２】
本発明は、ハード・アクチュエータ制約を明示的に考慮に入れる。
【００１３】
なお、本発明の諸態様を、単なる例に過ぎないが、添付の図面中に示す。
【発明を実施するための最良の形態】
【００１４】
図１は、本発明のプロセスに従って制御可能なシート製作機の一実施例として典型的な抄紙機１２を示す。縦方向ＭＤは、シートが製造されているとき、シート製作機を介して搬送される方向２０として定義される。横方向ＣＤは、ＭＤに直角な方向２２である。図示の抄紙機による紙シートの製造プロセス全体は、まず、木材パルプをマシンのウェット・エンド１４にあるヘッド・ボックス１に供給することを含む。ヘッド・ボックス１は、移動ワイヤ１６上に、そのパルプを抄紙機の幅にわたって薄く配するように動作する。抄紙機１２の残りでは、製造中の紙シートがそれを加熱し加圧するローラ列を介して搬送されるとき水分が除去されることにより、紙が形成される。仕上がった紙シートは、最終的に、マシンのドライ・エンド１８にある格納リール１１に巻き取られる。
【００１５】
抄紙プロセスを制御するために、シート特性をコンスタントに測定する必要があり、シートの品質を保証するために抄紙機を調整しなければならない。その制御は、一般に、製造プロセスの様々な段階でシート特性を測定することによって、またその測定された情報を用いて抄紙機内のアクチュエータを調整し、所望の目標からのシート特性の任意の変動を補償することによって達成される。図１の抄紙機１２では、シート特性のＣＤプロファイルを表す測定値配列を提供するために、２つの走査測定装置６および１０が使用される。新しいＣＤプロファイルは、通常１０秒から３０秒の範囲の均等な走査またはサンプリング間隔で取得される。典型的な測定プロファイルの実施例は、（乾燥繊維の）重量、水分、紙厚（キャリパ）、光沢、および平滑性である。測定値配列（ＣＤプロファイル）は、様々なサイズを有しており、通常６００から２０００エレメントに及ぶ。
【００１６】
本発明のプロセスは、品質管理システム（ＱＣＳ）コンピュータ２５中のソフトウェア・アプリケーションによって実施されることが好ましい。ＱＣＳは、本発明のプロセスが利用するある範囲のシステム・サービスを提供する。例えば、ＱＣＳによって提供される主要なシステム・サービスは、測定装置とアクチュエータへの通信インターフェースである。図１で、通信インターフェースは、アクチュエータとセンサを相互接続するためのＬＡＮに類似したネットワーク２６を含む。他の主要なシステム・サービスは、本発明に対するヒューマン・マシン・インターフェース（ＨＭＩ）である。走査装置またはセンサの固定アレイなどの測定装置は、マシン幅にわたりシート特性の測定を行う。測定装置は、通常、搭載コンピュータを有し、そのコンピュータは、信号処理を行い、ＱＣＳコンピュータへの通信インターフェースを提供する。２つのタイプのアクチュエータがある。第１に、シートの幅全体に影響するだけの、即ち、シート特性の平均値を変える縦方向（ＭＤ）アクチュエータである。第２に、マシンの幅全体にわたるアクチュエータ・アレイである横方向ＣＤアクチュエータ・ビームがある。ＣＤアクチュエータ・ビームは、シート特性の平均値とシート特性のＣＤ形状との双方に影響する。アクチュエータは、通常、規制的な制御に加えて隣接するアクチュエータおよびＱＣＳゲートウェイと通信する搭載コンピュータを有しており、インテリジェントである。このような構成は、Ｓｐｉｎｎｅｒ他による、「Ｄｉｓｔｒｉｂｕｔｅｄ Ｉｎｔｅｌｌｉｇｅｎｃｅ ａｃｔｕａｔｏｒ ｃｏｎｔｒｏｌｌｅｒ ｗｉｔｈ ｐｅｅｒ−ｔｏ−ｐｅｅｒ ａｃｔｕａｔｏｒ ｃｏｍｍｕｎｉｃａｔｉｏｎ（ピアツーピアのアクチュエータ通信を用いる分散型知的アクチュエータ制御装置」と題する米国特許第５，７７１，１７５号に全体的に開示されており、その開示を参照により本明細書に組み込む。
【００１７】
図１の抄紙機におけるＭＤアクチュエータの一実施例は、ヘッド・ボックス１のシック（ｔｈｉｃｋ）・ストック・バルブ２であり、それは、到着するパルプの粘度（ｃｏｎｓｉｓｔｅｎｃｙ）を制御し、後に、製造中の紙シートのＭＤ重量およびＭＤ水分に影響する。図示されている他のＭＤアクチュエータは、ドライヤ・セクションのスチーム流量バルブ８であり、それは、ドライヤ・セクションのロールにより与えられる熱を調整し、後に、シートのＭＤ水分およびＭＤ紙厚に影響する。
【００１８】
図１で、抄紙機１２の横方向に延長されるＣＤアクチュエータ・アレイの一実施例は、ヘッド・ボックス１上にマウントされたスライス・アクチュエータ３のアレイであり、それは、ヘッド・ボックス開口部領域を調整するように動作し、後に、シートのＣＤ重量、水分、および紙厚に影響する。スライス・アクチュエータ３はまた、パルプ・フローの速度が一定に維持される場合、ＭＤ重量および水分に影響する。ＣＤスチーム・アクチュエータ４は、シートに蒸気を与え、ＭＤ／ＣＤ水分、およびＣＤ紙厚に影響する。そのＣＤスチーム・アクチュエータ・ビームがカレンダ・スタック中に位置する場合、それはまた、ＭＤおよびＣＤ光沢および平滑性に影響することになる。ＣＤリウェット（ｒｅｗｅｔ）・アクチュエータ５および７は、シートに細かい水のスプレーを加え、ＭＤ／ＣＤ水分およびＣＤ紙厚に影響する。ＣＤリウェット・アクチュエータが、カレンダ・スタックのすぐ前に位置する場合、リウェット・アクチュエータはまた、ＭＤ／ＣＤ光沢および平滑性に影響することも可能である。図１に示す最後のＣＤアクチュエータ・アレイは、カレンダ・スタック中の誘導加熱ビーム９である。誘導加熱ビームは、ＣＤ紙厚、ＭＣ／ＣＤ光沢および平滑性に影響する。
【００１９】
様々なアレイ中のＣＤアクチュエータは、不均一な間隔で配置することができ、通常、３０から２００エレメントの範囲である。
【００２０】
本発明のプロセスは、シート特性のすべての測定値配列を取り込み、各アクチュエータが各シート特性に対して有する影響を考慮に入れて、すべてのＭＤアクチュエータおよびＣＤアクチュエータのためにアクチュエータ設定値を最適に計算することを含む。
【００２１】
図２は、抄紙機プロセスに組み込まれた、本発明のプロセスの閉ループのブロック図を示す。本発明のプロセスは、波線３０でマーク付けされた境界で定義される。抄紙機制御プロセスは、３１で概略を示す。
【００２２】
最初に、少なくとも１つの特性測定値配列が、本発明のプロセスへの入力配列として提供される。図２のブロック図で、入力配列として、例えば、重量と、水分と、紙厚とを表す３つのシート特性配列ｙ１（ｋ）と、ｙ２（ｋ）と、ｙ３（ｋ）とが、ステップ３５で提供される。ｋは、現在のサンプリング・インスタントを示す。当業者であれば、本発明のプロセスに入力されるシート特性配列の数は、測定され制御されるシート特性にのみ依存することが明らかであろう。
【００２３】
各シート特性測定値配列は、異なる次元を有することができる。ステップ３６で、シート特性測定値配列は、通常、まず時間フィルタＦｉを用いてフィルタされ、周知のフィルタリング技法を用いて、シート特性中のノイズおよび制御できないＭＤ変動が除去される。その時間フィルタリングは、本発明の部分ではないので、それを、他のＱＣＳシステム・サービスと見なすことができる。フィルタされたシート特性測定値配列は、以下で述べる本発明の共通解像度マッピング・コンポーネント３９への入力となる。
【００２４】
ステップ３８で、入力配列は、まず、共通解像度Ｎｙｃにマップされる。共通解像度は、好ましくは、正確な２次元プロセス・モデルを取得するために、最も高いアクチュエータ解像度の３倍を超えるものとすべきである。共通解像度マッピング・コンポーネント３９により、得られた共通解像度配列ｙｆ１（ｋ）、ｙｆ２（ｋ）、およびｙｆ３（ｋ）中に偽測定情報が現れないようにすることができる。
【００２５】
次いで、ステップ４０で、共通解像度測定値配列は結合されて、次元が１×３Ｎｙｃの１次元配列ｙｆ（ｋ）となる。結合された共通解像度測定値配列ｙｆ（ｋ）およびアクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）は、次いで、状態オブザーバ・コンポーネント４２に送られる。状態オブザーバ・コンポーネント４２は、結合された測定値配列ｙｆ（ｋ）およびアクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）に基づき、抄紙機プロセスの推定される現在の内部状態を表す配列ｘ（ｋ）を生成する。
【００２６】
各シート特性測定値配列は、未来のシート特性目標配列ｙ１ｒｅｆ（ｋ＋ｊ）と、ｙ２ｒｅｆ（ｋ＋ｊ）と、ｙ３ｒｅｆ（ｋ＋ｊ）とにそれぞれ関連づけられる。未来の目標配列は、抄紙機のオペレータによって提供される情報に基づくＱＣＳシステム・サービスとして提供される。ｊ＞０は、未来のサンプリング・インスタンスを示す。ステップ４４で、シート特性の共通解像度測定値配列と同様に、未来のシート特性目標配列は結合されて、１つのより大きな目標配列ｙｒｅｆ（ｋ＋ｊ）となる。
【００２７】
シート特性コンポーネント・セレクタ・モジュール４６により、ユーザは、本発明の制御装置がシート特性のＣＤ／ＭＤコンポーネントを共に制御すべき場合、ＣＤコンポーネントだけ、またはＭＤコンポーネントだけを指定することができる。シート特性コンポーネント・セレクタ・モジュール４６によると、目標配列ｙｒｅｆ（ｋ＋ｊ）および共通解像度測定値配列ｙｆ（ｋ）を、所望モードを達成するように変更することができる。
【００２８】
推定される現在の状態配列ｘ（ｋ）と、結合された未来のシート特性目標配列ｙｒｅｆ（ｋ）と、アクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）とは、ＣＤＭＤ−ＭＰＣコア・モジュール４８への入力として使用される。抄紙機プロセスのモデル５０、ならびにオブジェクト関数の重みおよびハード制約５２はまた、ＣＤＭＤ−ＭＰＣコア・モジュール４８への入力として働く。その情報に基づき、ＣＤＭＤ−ＭＰＣコア・モジュール４８は、アクチュエータの物理的制限（ハード制約）の下に、可能な限りその目標にすべてのシート特性を近付けるように最適に調整された設定値を生成する。
【００２９】
最適に調整された設定値の計算は、以下のサブ関数で行われる。
【００３０】
推定される現在の状態の現在の内部状態ｘ（ｋ）およびプロセス・モデルに基づき、ＣＤＭＤＭＰＣ予測モジュールは、シート特性の未来の予測ｙｐ（ｋ＋ｊ）を生成する。ただし、ｊ＞０は未来のサンプリング・インスタンスを表す。抄紙機プロセス・モデルは、以下の状態空間形式（Ａ、Ｂ、Ｃ、Ｎｄ）で表されることが好ましい。
【００３１】
ｘ（ｋ＋１）＝Ａｘ（ｋ）＋ＢΔｕ（ｋ−Ｎｄ） （１）
ｙ（ｋ）＝Ｃｘ（ｋ）
ただし、ｋはサンプリング・インスタンスであり、Ａは、プロセスの動的な時間情報を含む状態遷移行列であり、Ｂは、プロセスの静的な空間情報を含む状態入力行列であり、Ｃは、状態出力行列であり、Ｎｄは、サンプル中のプロセス・トランスポート遅延である。代替的に、抄紙機モデルは、インパルス応答モデル、ステップ応答モデル、または伝達（ｔｒａｎｓｆｅｒ）関数モデルなど、他の形式で表すこともできる。抄紙機プロセス・モデルは、２次元プロセス・モデルを識別するための自動化ツールを用いて取得することが好ましい。このような自動化ツールは、Ｇｏｒｉｎｅｖｓｋｙ Ｄ．、Ｈｅａｖｅｎ Ｅ．Ｍ．、Ｇｈｅｏｒｇｈｅ Ｃ、「Ｈｉｇｈ ｐｅｒｆｏｒｍａｎｃｅ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ ｃｒｏｓｓ−ｄｉｒｅｃｔｉｏｎａｌ ｐｒｏｃｅｓｓｅｓ」、Ｃｏｎｔｒｏｌ ｓｙｓｔｅｍｓ １９９８、Ｐｏｖｏｒｏ、フィンランド、１９９８年９月、による参照中で論じられており、その開示を参照により本明細書に組み込む。
【００３２】
シート特性の未来の予測ｙｐ（ｋ＋ｊ）は、次に、未来の目標配列ｙｒｅｆ（ｋ＋ｊ）、オブジェクト関数の重みＱｉ、最後のアクチュエータ設定値ｕ（ｋ−１）、ハード制約、およびオブジェクト関数Ｊ（ｔ）と共にＱＰ定式化モジュールに渡される。オブジェクト関数Ｊ（ｔ）は、次の形式であることが好ましい。
【００３３】
【数１】
【００３４】
ｅ（ｋ＋ｊ）＝ｙｒｅｆ（ｋ＋ｊ）−ｙｐ（ｋ＋ｊ）は、シート特性における未来の予測されるエラーである。Ｑ１は、異なるシート特性間およびシートの異なるＣＤ位置間の相対的な重要性を指定する重み付け行列である。Ｑ１を用いると、例えば、水分を重量よりも重要なものとして指定でき、またシートの中央が、シートの端部よりも重要であると指定可能である。Ｑ２は、２つの連続するサンプル・インスタンス間のアクチュエータ設定値の大きな変更によるコストを指定する重み付け行列である。Ｍは、重み付け行列Ｑ３と併せて、ユーザに、アクチュエータ設定値プロファイルにおける異なる空間方向に対するコストを指定可能にする行列である。Ａおよびｂは、ハード制約を指定する制約行列である。空間的に低利得の方向には、閉ループ・システムの空間的にロバストな安定性を可能とするために、高いコストに割り当てる必要がある。低利得方向は、Ｓｔｅｗａｒｔ ＧＥ、Ｂａｃｋｓｔｒｏｍ Ｊ．Ｕ、Ｂａｋｅｒ Ｐ、Ｇｈｅｏｒｇｈｅ Ｃ、およびＶｙｓｅ Ｒ．Ｎ．、「Ｃｏｎｔｒｏｌｌａｂｉｌｉｔｙ ｉｎ ｃｒｏｓｓ−ｄｉｒｅｃｔｉｏｎａｌ ｐｒｏｃｅｓｓｅｓ：Ｐｒａｃｔｉｃａｌ ｒｕｌｅｓ ｆｏｒ ａｎａｌｙｓｉｓ ａｎｄ ｄｅｓｉｇｎ」、８７ｔｈ Ａｎｎｕａｌ Ｍｅｅｔｉｎｇ、ＰＡＰＴＡＣ、モントリオール、ケベック州、２００１年２月、による参照中で論じられている短い空間波長に相当し、その開示を参照により本明細書に組み込む。Ｑ４は、参照または目標設定値からずらしたアクチュエータ設定値のコストを指定する重み付け行列である。ＣＤアクチュエータ・アレイに対して、アクチュエータのエネルギ消費の観点から、またはシート製作機の動作可能性の観点から、関連するアクチュエータの目標値設定するのが普通である。Ｓは、重み付け行列Ｑ５と併せて、ユーザに、特定のアクチュエータ内設定方向にＣＤアクチュエータ・アレイおよびＭＤアクチュエータを移動するコストを指定可能にする行列である。ロバストな安定性を保証するために、低いアクチュエータ内設定利得方向へ移動するためには高いコストを割り当てる必要がある。特定のシート製作プロセスに対するアクチュエータ内設定方向性の影響は、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｈｅｎｄｅｒｓｏｎ Ｂ、およびＳｔｅｗａｒｔ Ｃ、「Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ａｎｄ ｍｕｌｔｉｖａｒｉａｂｌｅ ｃｏｎｔｒｏｌ ｏｆ ｓｕｐｅｒｃａｌｅｎｄｅｒｓ」、Ｃｏｎｔｒｏｌ Ｓｙｓｔｅｍｓ ２００２、２００２年６月、ストックホルム、 スウェーデン、による参照中で論じられており、その開示を参照により本明細書に組み込む。
【００３５】
本発明のプロセスで考慮されているハード制約は、次のものである。
【００３６】
１．例えば、オペレータ制御下のまたは故障している、本発明の制御下にないアクチュエータを制御装置の方へ移動してはならない。
【００３７】
２．アクチュエータの設定値は、その物理的な上限および下限内でなければならない。
【００３８】
３．１次および２次の曲げ限界（ＣＤアクチュエータ・ビームにのみ適用可能）。
【００３９】
４．アクチュエータ設定値平均を特定の制限で、または指定された範囲内で維持すること（ＣＤアクチュエータ・アレイにのみ適用可能）。
【００４０】
５．アクチュエータ設定値変更の最大値。
【００４１】
ＱＰ定式化モジュールはこれらの入力を取り込み、標準形式で２次計画法を作成する。
【００４２】
【数２】
【００４３】
ここで、Φはヘッセ行列であり、φはヤコビ行列であり、Ａおよびｂは、制約行列である。
【００４４】
式（３）の２次計画法は、Ｂａｒｔｌｅｔｔ Ｒ． Ａ、Ｂｉｅｇｌｅｒ Ｌ． Ｔ．、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｇｏｐａｌ Ｖ、「Ｑｕａｄｒａｔｉｃ ｐｒｏｇｒａｍｍｉｎｇ ａｌｇｏｒｉｔｈｍｓ ｆｏｒ ｌａｒｇｅ−ｓｃａｌｅ ｍｏｄｅｌ ｐｒｅｄｉｃｔｉｖｅ ｃｏｎｔｒｏｌ」、Ｊｏｕｒｎａｌ ｏｆ Ｐｒｏｃｅｓｓ Ｃｏｎｔｒｏｌ、１２（２００２） ７７５−７９５、による参照中で論じられている高度にカスタム化されたＱＰソルバ（ｓｏｌｖｅｒ）を用いて解かれる。２次計画法を解くことにより、シート製作プロセスの調整されたＭＤおよびＣＤ制御のためのアクチュエータ設定値Δｕ（ｔ）の最適な変更配列を得ることができる。
【００４５】
ステップ５４で、アクチュエータ設定値Δｕ（ｔ）の最適な変更配列は、次いで、ｕ（ｔ）を形成するためにアクチュエータ設定値の最後の配列ｕ（ｔ―１）に加えられ、ステップ５６で、それは、抄紙機プロセス中の異なるＭＤアクチュエータおよびＣＤアクチュエータ・アレイに送達するために、次いで、設定値ｕｉ（ｔ）に分割される。
【００４６】
本発明を、明確化し理解しやすくするために、例として詳細に説明してきたが、添付の請求項の範囲に含まれる幾つかの変更および修正を実施できることは明らかであろう。
【図面の簡単な説明】
【００４７】
【図１】本発明のプロセスに従って動作可能な典型的なシート製作機の概略図である。
【図２】本発明の諸プロセス・ステップを示すブロック図である。
【技術分野】
【０００１】
本発明は、シート製作プロセスの制御に関連し、より詳細には、シート製作機における縦方向および横方向アクチュエータのオペレーションを調整するのための方法に関連する。
【背景技術】
【０００２】
シート製作機におけるシート特性の制御は、シートの特性を可能な限り目標値に近い値に保つことに関する。シート特性の制御には、２組の異なるアクチュエータが使用される。第１に、シート特性の横方向（ＣＤ：ｃｒｏｓｓ ｄｉｒｅｃｔｉｏｎ）平均値にだけ影響する縦方向（ＭＤ：ｍａｃｈｉｎｅ ｄｉｒｅｃｔｉｏｎ）アクチュエータがある。各ＭＤアクチュエータは、シート特性において異なる動的応答を有することができる。第２に、ＣＤ方向にシートを横断して配列されたＣＤアクチュエータがある。ＣＤアクチュエータの各アレイは、シート特性の平均値及びＣＤ形状の双方に影響することができる。ＣＤアクチュエータは、シート特性における異なる動的応答および異なる空間応答を有することができる。シート特性の全体的な制御問題は、非常に多変量になる。即ち、ＣＤアレイ中の１つのＣＤアクチュエータは、幾つかのシート特性において隣接するＣＤゾーンに影響し、また特定のシート特性を制御するよう意図したＣＤアクチュエータの平均的な効果は、幾つかのＭＤアクチュエータによっても影響される幾つかのシート特性の平均値に影響する可能性がある。その問題はまた、非常に大規模なものである。典型的な制御プロセスは、数千の出力（シート特性の測定値）、および数百の入力（アクチュエータの設定値）を有する可能性がある。プロセスはまた、特定の空間方向およびアクチュエータ内設定方向に制御することは困難または不可能である。
【０００３】
現在、従来の大部分のシート製作装置では、シート特性の制御を２つの制御問題に分離する。第１に、ＣＤ平均値は、シート特性のＣＤ平均値に対するＣＤアクチュエータの効果を利用せずに、ＭＤアクチュエータを利用するだけで制御される。第２に、シートを横断して配列されたＣＤアクチュエータは、シート特性の平均値まわりのＣＤ変動を制御するためだけに利用される。現在、ＭＤアクチュエータを調整するための明示的なハード制約処理とともに、モデル予測制御を用いるＭＤ制御方式が利用できる。
【０００４】
モデル予測制御を用いて１つまたは複数のシート特性を制御する、ＣＤアクチュエータ・アレイの最適に調整された制御は、例えば、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｈｅｎｄｅｒｓｏｎ Ｂ、およびＳｔｅｗａｒｔ Ｃ、「Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ａｎｄ ｍｕｌｔｉｖａｒｉａｂｌｅ ｃｏｎｔｒｏｌ ｏｆ ｓｕｐｅｒｃａｌｅｎｄｅｒｓ」、Ｃｏｎｔｒｏｌ Ｓｙｓｔｅｍｓ ２００２、２００２年６月、ストックホルム、スウェーデン、および、Ｂａｃｋｓｔｒｏｍ Ｊ． Ｕ、Ｇｈｅｏｒｇｈｅ Ｃ、Ｓｔｅｗａｒｔ Ｇ． Ｅ、Ｖｙｓｅ Ｒ．Ｎ、「Ｃｏｎｓｔｒａｉｎｅｄ ｍｏｄｅｌ ｐｒｅｄｉｃｔｉｖｅ ｃｏｎｔｒｏｌ ｆｏｒ ｃｒｏｓｓ ｄｉｒｅｃｔｉｏｎａｌ ｍｕｌｔｉ−ａｒｒａｙ ｐｒｏｃｅｓｓｅｓ」、Ｐｕｌｐ ＆ Ｐａｐｅｒ Ｃａｎａｄａ、Ｔ１２８ １０２：５ （２００１） などの記事で論じられている。
【０００５】
ＭＤアクチュエータおよびＣＤアクチュエータを調整する必要性は、本願の権利者が所有する２０００年７月２５日に発行された米国特許第６，０９４，６０４号で明らかにされている。この問題に対して提案される解決策はまた、互いに通信するアクチュエータにおける分散型局所化インテリジェント制御装置のシステムを含む‘６０４特許中に開示されている。
【発明の開示】
【課題を解決するための手段】
【０００６】
（発明の概要）
上記で概要を説明した問題に対処するために、本発明は、調整されたＭＤおよびＣＤ制御のための柔軟な大規模多変量モデル予測制御装置を提供し、その制御装置は、入力としてシート特性測定値の複数の配列を取り込み、複数の出力配列（アクチュエータ設定値）を生成する。その配列は、任意の次元のものとすることができる。１つのＭＤ配列は、１×１の配列と見なすことができる。任意の数の入力および出力配列とすることができる。本発明は、新しく最適に調整される設定値を均等な制御間隔で計算する。各シート特性に対して、ＣＤコンポーネントだけ、ＭＤコンポーネントだけ、またはＭＤ及びＣＤコンポーネントの双方に、制御することができる。本発明は、予測区間（ｐｒｅｄｉｃｔｉｏｎ ｈｏｒｉｚｏｎ）Ｈｐから未来の制御区間Ｈｃにおけるアクチュエータ設定値に対する、動的及び２次元応答を予測する。次いで、本発明は、未来の予測シート特性を可能な限り目標に近付ける未来の最適な設定値を計算する。制御装置はまた、アクチュエータに対する物理的制限を明示的に考慮に入れる。制御装置は、アクチュエータ設定値を困難な空間方向およびアクチュエータ内設定プロセス方向で発行することを回避することにより、２つのタイプの方向問題を処理する。これにより、閉ループの２次元ロバストな安定性を得ることができる。
【０００７】
従って、本発明は、１枚の材料シートを製作するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスを提供し、そのプロセスは、
シート測定データを収集するために規則的な間隔で複数のシート特性を測定するステップと、
そのシート測定データを、複数のシート特性測定値配列を確立するために処理するステップと、
１次元共通解像度測定値配列を確立するためにシートの特性測定値配列を処理するステップと、
シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
未来のシート特性の目標配列を確立するステップと、
シート製造プロセスの推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
シート特性の未来の予測配列、未来のシート特性目標配列、および前のアクチュエータ設定値配列を、可解なオブジェクト関数に入力して、シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含む。
【０００８】
本発明はまた、１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスを提供し、そのプロセスは、
シート測定データを収集するために規則的な間隔で複数のシート特性を測定するステップと、
複数のシート特性測定値配列を確立するために、シート測定データを処理するステップと、
シート特性測定値配列を共通解像度にマップするステップと、
共通解像度シート特性測定値配列を、１つのより大きい１次元共通解像度測定値配列に結合するステップと、
共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力することにより、シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
複数の未来のシート特性目標配列を１つの目標配列に結合するステップと、
シート製造プロセスの推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
シート特性の未来の予測配列、目標配列、オブジェクト関数の重み、最後のアクチュエータ設定値配列、およびハード制約を、オブジェクト関数に入力するステップと、
オブジェクト関数を解いて、シート製作プロセスの調整されたＭＤおよびＣＤ制御のために現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含む。
【０００９】
本発明は、ＣＤアクチュエータ・アレイおよびＭＤアクチュエータを最適に操作しかつ調整するように動作し、シート特性におけるＭＤおよびＣＤ変動を最小化することができる。
【００１０】
本発明は、ＭＤアクチュエータとＣＤアクチュエータ・アレイの間の対話を最適に調整する。本発明は更に、空間利得の低い方向に移動するコストを表現するためにオブジェクティト関数中に一般（ｇｅｎｅｒａｌ）重み付け関数を有する。本発明は更に、小さなアクチュエータ内設定方向に移動するコストを表現するための明示的な重み付け関数を有する。本発明は更に、ＣＤアクチュエータ・アレイ設定値平均の許容可能範囲を指定するハード制約を含む。本発明は、シート特性のＣＤコンポーネントだけ、ＭＤコンポーネントだけ、あるいはＣＤ及びＭＤコンポーネントの双方に制御するようにセットアップすることができる。
【００１１】
本発明のプロセスは、複数の分散された制御装置ではなく、１つの集中化された制御装置を使用することが好ましい。
【００１２】
本発明は、ハード・アクチュエータ制約を明示的に考慮に入れる。
【００１３】
なお、本発明の諸態様を、単なる例に過ぎないが、添付の図面中に示す。
【発明を実施するための最良の形態】
【００１４】
図１は、本発明のプロセスに従って制御可能なシート製作機の一実施例として典型的な抄紙機１２を示す。縦方向ＭＤは、シートが製造されているとき、シート製作機を介して搬送される方向２０として定義される。横方向ＣＤは、ＭＤに直角な方向２２である。図示の抄紙機による紙シートの製造プロセス全体は、まず、木材パルプをマシンのウェット・エンド１４にあるヘッド・ボックス１に供給することを含む。ヘッド・ボックス１は、移動ワイヤ１６上に、そのパルプを抄紙機の幅にわたって薄く配するように動作する。抄紙機１２の残りでは、製造中の紙シートがそれを加熱し加圧するローラ列を介して搬送されるとき水分が除去されることにより、紙が形成される。仕上がった紙シートは、最終的に、マシンのドライ・エンド１８にある格納リール１１に巻き取られる。
【００１５】
抄紙プロセスを制御するために、シート特性をコンスタントに測定する必要があり、シートの品質を保証するために抄紙機を調整しなければならない。その制御は、一般に、製造プロセスの様々な段階でシート特性を測定することによって、またその測定された情報を用いて抄紙機内のアクチュエータを調整し、所望の目標からのシート特性の任意の変動を補償することによって達成される。図１の抄紙機１２では、シート特性のＣＤプロファイルを表す測定値配列を提供するために、２つの走査測定装置６および１０が使用される。新しいＣＤプロファイルは、通常１０秒から３０秒の範囲の均等な走査またはサンプリング間隔で取得される。典型的な測定プロファイルの実施例は、（乾燥繊維の）重量、水分、紙厚（キャリパ）、光沢、および平滑性である。測定値配列（ＣＤプロファイル）は、様々なサイズを有しており、通常６００から２０００エレメントに及ぶ。
【００１６】
本発明のプロセスは、品質管理システム（ＱＣＳ）コンピュータ２５中のソフトウェア・アプリケーションによって実施されることが好ましい。ＱＣＳは、本発明のプロセスが利用するある範囲のシステム・サービスを提供する。例えば、ＱＣＳによって提供される主要なシステム・サービスは、測定装置とアクチュエータへの通信インターフェースである。図１で、通信インターフェースは、アクチュエータとセンサを相互接続するためのＬＡＮに類似したネットワーク２６を含む。他の主要なシステム・サービスは、本発明に対するヒューマン・マシン・インターフェース（ＨＭＩ）である。走査装置またはセンサの固定アレイなどの測定装置は、マシン幅にわたりシート特性の測定を行う。測定装置は、通常、搭載コンピュータを有し、そのコンピュータは、信号処理を行い、ＱＣＳコンピュータへの通信インターフェースを提供する。２つのタイプのアクチュエータがある。第１に、シートの幅全体に影響するだけの、即ち、シート特性の平均値を変える縦方向（ＭＤ）アクチュエータである。第２に、マシンの幅全体にわたるアクチュエータ・アレイである横方向ＣＤアクチュエータ・ビームがある。ＣＤアクチュエータ・ビームは、シート特性の平均値とシート特性のＣＤ形状との双方に影響する。アクチュエータは、通常、規制的な制御に加えて隣接するアクチュエータおよびＱＣＳゲートウェイと通信する搭載コンピュータを有しており、インテリジェントである。このような構成は、Ｓｐｉｎｎｅｒ他による、「Ｄｉｓｔｒｉｂｕｔｅｄ Ｉｎｔｅｌｌｉｇｅｎｃｅ ａｃｔｕａｔｏｒ ｃｏｎｔｒｏｌｌｅｒ ｗｉｔｈ ｐｅｅｒ−ｔｏ−ｐｅｅｒ ａｃｔｕａｔｏｒ ｃｏｍｍｕｎｉｃａｔｉｏｎ（ピアツーピアのアクチュエータ通信を用いる分散型知的アクチュエータ制御装置」と題する米国特許第５，７７１，１７５号に全体的に開示されており、その開示を参照により本明細書に組み込む。
【００１７】
図１の抄紙機におけるＭＤアクチュエータの一実施例は、ヘッド・ボックス１のシック（ｔｈｉｃｋ）・ストック・バルブ２であり、それは、到着するパルプの粘度（ｃｏｎｓｉｓｔｅｎｃｙ）を制御し、後に、製造中の紙シートのＭＤ重量およびＭＤ水分に影響する。図示されている他のＭＤアクチュエータは、ドライヤ・セクションのスチーム流量バルブ８であり、それは、ドライヤ・セクションのロールにより与えられる熱を調整し、後に、シートのＭＤ水分およびＭＤ紙厚に影響する。
【００１８】
図１で、抄紙機１２の横方向に延長されるＣＤアクチュエータ・アレイの一実施例は、ヘッド・ボックス１上にマウントされたスライス・アクチュエータ３のアレイであり、それは、ヘッド・ボックス開口部領域を調整するように動作し、後に、シートのＣＤ重量、水分、および紙厚に影響する。スライス・アクチュエータ３はまた、パルプ・フローの速度が一定に維持される場合、ＭＤ重量および水分に影響する。ＣＤスチーム・アクチュエータ４は、シートに蒸気を与え、ＭＤ／ＣＤ水分、およびＣＤ紙厚に影響する。そのＣＤスチーム・アクチュエータ・ビームがカレンダ・スタック中に位置する場合、それはまた、ＭＤおよびＣＤ光沢および平滑性に影響することになる。ＣＤリウェット（ｒｅｗｅｔ）・アクチュエータ５および７は、シートに細かい水のスプレーを加え、ＭＤ／ＣＤ水分およびＣＤ紙厚に影響する。ＣＤリウェット・アクチュエータが、カレンダ・スタックのすぐ前に位置する場合、リウェット・アクチュエータはまた、ＭＤ／ＣＤ光沢および平滑性に影響することも可能である。図１に示す最後のＣＤアクチュエータ・アレイは、カレンダ・スタック中の誘導加熱ビーム９である。誘導加熱ビームは、ＣＤ紙厚、ＭＣ／ＣＤ光沢および平滑性に影響する。
【００１９】
様々なアレイ中のＣＤアクチュエータは、不均一な間隔で配置することができ、通常、３０から２００エレメントの範囲である。
【００２０】
本発明のプロセスは、シート特性のすべての測定値配列を取り込み、各アクチュエータが各シート特性に対して有する影響を考慮に入れて、すべてのＭＤアクチュエータおよびＣＤアクチュエータのためにアクチュエータ設定値を最適に計算することを含む。
【００２１】
図２は、抄紙機プロセスに組み込まれた、本発明のプロセスの閉ループのブロック図を示す。本発明のプロセスは、波線３０でマーク付けされた境界で定義される。抄紙機制御プロセスは、３１で概略を示す。
【００２２】
最初に、少なくとも１つの特性測定値配列が、本発明のプロセスへの入力配列として提供される。図２のブロック図で、入力配列として、例えば、重量と、水分と、紙厚とを表す３つのシート特性配列ｙ１（ｋ）と、ｙ２（ｋ）と、ｙ３（ｋ）とが、ステップ３５で提供される。ｋは、現在のサンプリング・インスタントを示す。当業者であれば、本発明のプロセスに入力されるシート特性配列の数は、測定され制御されるシート特性にのみ依存することが明らかであろう。
【００２３】
各シート特性測定値配列は、異なる次元を有することができる。ステップ３６で、シート特性測定値配列は、通常、まず時間フィルタＦｉを用いてフィルタされ、周知のフィルタリング技法を用いて、シート特性中のノイズおよび制御できないＭＤ変動が除去される。その時間フィルタリングは、本発明の部分ではないので、それを、他のＱＣＳシステム・サービスと見なすことができる。フィルタされたシート特性測定値配列は、以下で述べる本発明の共通解像度マッピング・コンポーネント３９への入力となる。
【００２４】
ステップ３８で、入力配列は、まず、共通解像度Ｎｙｃにマップされる。共通解像度は、好ましくは、正確な２次元プロセス・モデルを取得するために、最も高いアクチュエータ解像度の３倍を超えるものとすべきである。共通解像度マッピング・コンポーネント３９により、得られた共通解像度配列ｙｆ１（ｋ）、ｙｆ２（ｋ）、およびｙｆ３（ｋ）中に偽測定情報が現れないようにすることができる。
【００２５】
次いで、ステップ４０で、共通解像度測定値配列は結合されて、次元が１×３Ｎｙｃの１次元配列ｙｆ（ｋ）となる。結合された共通解像度測定値配列ｙｆ（ｋ）およびアクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）は、次いで、状態オブザーバ・コンポーネント４２に送られる。状態オブザーバ・コンポーネント４２は、結合された測定値配列ｙｆ（ｋ）およびアクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）に基づき、抄紙機プロセスの推定される現在の内部状態を表す配列ｘ（ｋ）を生成する。
【００２６】
各シート特性測定値配列は、未来のシート特性目標配列ｙ１ｒｅｆ（ｋ＋ｊ）と、ｙ２ｒｅｆ（ｋ＋ｊ）と、ｙ３ｒｅｆ（ｋ＋ｊ）とにそれぞれ関連づけられる。未来の目標配列は、抄紙機のオペレータによって提供される情報に基づくＱＣＳシステム・サービスとして提供される。ｊ＞０は、未来のサンプリング・インスタンスを示す。ステップ４４で、シート特性の共通解像度測定値配列と同様に、未来のシート特性目標配列は結合されて、１つのより大きな目標配列ｙｒｅｆ（ｋ＋ｊ）となる。
【００２７】
シート特性コンポーネント・セレクタ・モジュール４６により、ユーザは、本発明の制御装置がシート特性のＣＤ／ＭＤコンポーネントを共に制御すべき場合、ＣＤコンポーネントだけ、またはＭＤコンポーネントだけを指定することができる。シート特性コンポーネント・セレクタ・モジュール４６によると、目標配列ｙｒｅｆ（ｋ＋ｊ）および共通解像度測定値配列ｙｆ（ｋ）を、所望モードを達成するように変更することができる。
【００２８】
推定される現在の状態配列ｘ（ｋ）と、結合された未来のシート特性目標配列ｙｒｅｆ（ｋ）と、アクチュエータ設定値の過去の変更配列Δｕｄ（ｋ）とは、ＣＤＭＤ−ＭＰＣコア・モジュール４８への入力として使用される。抄紙機プロセスのモデル５０、ならびにオブジェクト関数の重みおよびハード制約５２はまた、ＣＤＭＤ−ＭＰＣコア・モジュール４８への入力として働く。その情報に基づき、ＣＤＭＤ−ＭＰＣコア・モジュール４８は、アクチュエータの物理的制限（ハード制約）の下に、可能な限りその目標にすべてのシート特性を近付けるように最適に調整された設定値を生成する。
【００２９】
最適に調整された設定値の計算は、以下のサブ関数で行われる。
【００３０】
推定される現在の状態の現在の内部状態ｘ（ｋ）およびプロセス・モデルに基づき、ＣＤＭＤＭＰＣ予測モジュールは、シート特性の未来の予測ｙｐ（ｋ＋ｊ）を生成する。ただし、ｊ＞０は未来のサンプリング・インスタンスを表す。抄紙機プロセス・モデルは、以下の状態空間形式（Ａ、Ｂ、Ｃ、Ｎｄ）で表されることが好ましい。
【００３１】
ｘ（ｋ＋１）＝Ａｘ（ｋ）＋ＢΔｕ（ｋ−Ｎｄ） （１）
ｙ（ｋ）＝Ｃｘ（ｋ）
ただし、ｋはサンプリング・インスタンスであり、Ａは、プロセスの動的な時間情報を含む状態遷移行列であり、Ｂは、プロセスの静的な空間情報を含む状態入力行列であり、Ｃは、状態出力行列であり、Ｎｄは、サンプル中のプロセス・トランスポート遅延である。代替的に、抄紙機モデルは、インパルス応答モデル、ステップ応答モデル、または伝達（ｔｒａｎｓｆｅｒ）関数モデルなど、他の形式で表すこともできる。抄紙機プロセス・モデルは、２次元プロセス・モデルを識別するための自動化ツールを用いて取得することが好ましい。このような自動化ツールは、Ｇｏｒｉｎｅｖｓｋｙ Ｄ．、Ｈｅａｖｅｎ Ｅ．Ｍ．、Ｇｈｅｏｒｇｈｅ Ｃ、「Ｈｉｇｈ ｐｅｒｆｏｒｍａｎｃｅ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ ｃｒｏｓｓ−ｄｉｒｅｃｔｉｏｎａｌ ｐｒｏｃｅｓｓｅｓ」、Ｃｏｎｔｒｏｌ ｓｙｓｔｅｍｓ １９９８、Ｐｏｖｏｒｏ、フィンランド、１９９８年９月、による参照中で論じられており、その開示を参照により本明細書に組み込む。
【００３２】
シート特性の未来の予測ｙｐ（ｋ＋ｊ）は、次に、未来の目標配列ｙｒｅｆ（ｋ＋ｊ）、オブジェクト関数の重みＱｉ、最後のアクチュエータ設定値ｕ（ｋ−１）、ハード制約、およびオブジェクト関数Ｊ（ｔ）と共にＱＰ定式化モジュールに渡される。オブジェクト関数Ｊ（ｔ）は、次の形式であることが好ましい。
【００３３】
【数１】
【００３４】
ｅ（ｋ＋ｊ）＝ｙｒｅｆ（ｋ＋ｊ）−ｙｐ（ｋ＋ｊ）は、シート特性における未来の予測されるエラーである。Ｑ１は、異なるシート特性間およびシートの異なるＣＤ位置間の相対的な重要性を指定する重み付け行列である。Ｑ１を用いると、例えば、水分を重量よりも重要なものとして指定でき、またシートの中央が、シートの端部よりも重要であると指定可能である。Ｑ２は、２つの連続するサンプル・インスタンス間のアクチュエータ設定値の大きな変更によるコストを指定する重み付け行列である。Ｍは、重み付け行列Ｑ３と併せて、ユーザに、アクチュエータ設定値プロファイルにおける異なる空間方向に対するコストを指定可能にする行列である。Ａおよびｂは、ハード制約を指定する制約行列である。空間的に低利得の方向には、閉ループ・システムの空間的にロバストな安定性を可能とするために、高いコストに割り当てる必要がある。低利得方向は、Ｓｔｅｗａｒｔ ＧＥ、Ｂａｃｋｓｔｒｏｍ Ｊ．Ｕ、Ｂａｋｅｒ Ｐ、Ｇｈｅｏｒｇｈｅ Ｃ、およびＶｙｓｅ Ｒ．Ｎ．、「Ｃｏｎｔｒｏｌｌａｂｉｌｉｔｙ ｉｎ ｃｒｏｓｓ−ｄｉｒｅｃｔｉｏｎａｌ ｐｒｏｃｅｓｓｅｓ：Ｐｒａｃｔｉｃａｌ ｒｕｌｅｓ ｆｏｒ ａｎａｌｙｓｉｓ ａｎｄ ｄｅｓｉｇｎ」、８７ｔｈ Ａｎｎｕａｌ Ｍｅｅｔｉｎｇ、ＰＡＰＴＡＣ、モントリオール、ケベック州、２００１年２月、による参照中で論じられている短い空間波長に相当し、その開示を参照により本明細書に組み込む。Ｑ４は、参照または目標設定値からずらしたアクチュエータ設定値のコストを指定する重み付け行列である。ＣＤアクチュエータ・アレイに対して、アクチュエータのエネルギ消費の観点から、またはシート製作機の動作可能性の観点から、関連するアクチュエータの目標値設定するのが普通である。Ｓは、重み付け行列Ｑ５と併せて、ユーザに、特定のアクチュエータ内設定方向にＣＤアクチュエータ・アレイおよびＭＤアクチュエータを移動するコストを指定可能にする行列である。ロバストな安定性を保証するために、低いアクチュエータ内設定利得方向へ移動するためには高いコストを割り当てる必要がある。特定のシート製作プロセスに対するアクチュエータ内設定方向性の影響は、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｈｅｎｄｅｒｓｏｎ Ｂ、およびＳｔｅｗａｒｔ Ｃ、「Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ａｎｄ ｍｕｌｔｉｖａｒｉａｂｌｅ ｃｏｎｔｒｏｌ ｏｆ ｓｕｐｅｒｃａｌｅｎｄｅｒｓ」、Ｃｏｎｔｒｏｌ Ｓｙｓｔｅｍｓ ２００２、２００２年６月、ストックホルム、 スウェーデン、による参照中で論じられており、その開示を参照により本明細書に組み込む。
【００３５】
本発明のプロセスで考慮されているハード制約は、次のものである。
【００３６】
１．例えば、オペレータ制御下のまたは故障している、本発明の制御下にないアクチュエータを制御装置の方へ移動してはならない。
【００３７】
２．アクチュエータの設定値は、その物理的な上限および下限内でなければならない。
【００３８】
３．１次および２次の曲げ限界（ＣＤアクチュエータ・ビームにのみ適用可能）。
【００３９】
４．アクチュエータ設定値平均を特定の制限で、または指定された範囲内で維持すること（ＣＤアクチュエータ・アレイにのみ適用可能）。
【００４０】
５．アクチュエータ設定値変更の最大値。
【００４１】
ＱＰ定式化モジュールはこれらの入力を取り込み、標準形式で２次計画法を作成する。
【００４２】
【数２】
【００４３】
ここで、Φはヘッセ行列であり、φはヤコビ行列であり、Ａおよびｂは、制約行列である。
【００４４】
式（３）の２次計画法は、Ｂａｒｔｌｅｔｔ Ｒ． Ａ、Ｂｉｅｇｌｅｒ Ｌ． Ｔ．、Ｂａｃｋｓｔｒｏｍ Ｊ、Ｇｏｐａｌ Ｖ、「Ｑｕａｄｒａｔｉｃ ｐｒｏｇｒａｍｍｉｎｇ ａｌｇｏｒｉｔｈｍｓ ｆｏｒ ｌａｒｇｅ−ｓｃａｌｅ ｍｏｄｅｌ ｐｒｅｄｉｃｔｉｖｅ ｃｏｎｔｒｏｌ」、Ｊｏｕｒｎａｌ ｏｆ Ｐｒｏｃｅｓｓ Ｃｏｎｔｒｏｌ、１２（２００２） ７７５−７９５、による参照中で論じられている高度にカスタム化されたＱＰソルバ（ｓｏｌｖｅｒ）を用いて解かれる。２次計画法を解くことにより、シート製作プロセスの調整されたＭＤおよびＣＤ制御のためのアクチュエータ設定値Δｕ（ｔ）の最適な変更配列を得ることができる。
【００４５】
ステップ５４で、アクチュエータ設定値Δｕ（ｔ）の最適な変更配列は、次いで、ｕ（ｔ）を形成するためにアクチュエータ設定値の最後の配列ｕ（ｔ―１）に加えられ、ステップ５６で、それは、抄紙機プロセス中の異なるＭＤアクチュエータおよびＣＤアクチュエータ・アレイに送達するために、次いで、設定値ｕｉ（ｔ）に分割される。
【００４６】
本発明を、明確化し理解しやすくするために、例として詳細に説明してきたが、添付の請求項の範囲に含まれる幾つかの変更および修正を実施できることは明らかであろう。
【図面の簡単な説明】
【００４７】
【図１】本発明のプロセスに従って動作可能な典型的なシート製作機の概略図である。
【図２】本発明の諸プロセス・ステップを示すブロック図である。
【特許請求の範囲】
【請求項１】
１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスであって、
規則的な間隔で複数のシート特性を測定し、シートの測定データを収集するステップと、
前記シート測定データを処理し、複数のシート特性測定値配列を確立するステップと、
前記シート特性測定値配列を処理し、１次元共通解像度測定値配列を確立するステップと、
前記シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
未来のシート特性目標配列を確立するステップと、
前記シート製造プロセスの前記推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
前記シート特性の未来の予測配列、前記未来のシート特性目標配列、および前のアクチュエータ設定値配列を、可解なオブジェクト関数に入力し、前記シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含むプロセス。
【請求項２】
請求項１に記載のプロセスであって、前記シート特性測定値配列を処理し、１次元共通解像度測定値配列を確立する前記ステップは、
前記シート特性測定値配列を共通解像度にマップするステップと、
前記共通解像度シート特性測定値配列を、より大きい１次元共通解像度測定値配列に結合するステップと、
を含む、プロセス。
【請求項３】
請求項２に記載のプロセスであって、前記シート特性測定値配列を共通解像度にマップする前記ステップは、最も高いアクチュエータ解像度の３倍を超える前記共通解像度を選択するステップを含む、プロセス。
【請求項４】
請求項１に記載のプロセスであって、前記シート製造プロセスの前記推定される現在の内部状態配列を生成する前記ステップは、前記共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力するステップを含む、プロセス。
【請求項５】
請求項１に記載のプロセスであって、未来のシート特性目標配列を確立する前記ステップは、複数の未来のシート特性目標配列を１つの目標配列に結合するステップを含む、プロセス。
【請求項６】
請求項１に記載のプロセスであって、前記シート特性の未来の予測配列、前記未来のシート特性目標配列、および前記前のアクチュエータ設定値配列を、シート特性のオブジェクト関数に入力する前記ステップは、オブジェクト関数の重みおよびハード制約を入力するステップを含む、プロセス。
【請求項７】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、以下の状態空間形式（Ａ、Ｂ、Ｃ、Ｎｄ）；
ｘ（ｋ＋１）＝Ａｘ（ｋ）＋ＢΔｕ（ｋ−Ｎｄ）
ｙ（ｋ）＝Ｃｘ（ｋ）
で表され、
式中、
ｋは、サンプリング・インスタンスであり、
ｘは、前記プロセスの前記推定される現在の内部状態配列であり、
Δｕは、アクチュエータ設定値の過去の変更配列であり、
Ａは、前記プロセスの動的な時間情報を含む状態遷移行列であり、
Ｂは、前記プロセスの静的な空間情報を含む状態入力行列であり、
Ｃは、状態出力行列であり、
Ｎｄは、サンプル中のプロセス・トランスポート遅延である、プロセス。
【請求項８】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、インパルス応答モデルの形式で表される、プロセス。
【請求項９】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、ステップ応答モデルの形式で表される、プロセス。
【請求項１０】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、伝達関数モデルの形式で表される、プロセス。
【請求項１１】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、２次元プロセス・モデルを識別するための自動化ツールを用いて生成される、プロセス。
【請求項１２】
請求項１に記載のプロセスであって、前記オブジェクト関数は、
【数１】
で表され、
式中、
ｅ（ｋ＋ｊ）＝ｙｒｅｆ（ｋ＋ｊ）−ｙｐ（ｋ＋ｊ）は、シート特性における未来の予測エラーであり、
Ｑ１は、前記シートの異なるシート特性間および異なるＣＤ位置間の相対的な重要性を指定する重み付け行列であり、
Ｑ２は、２つの連続するサンプル・インスタンス間における大きな前記アクチュエータ設定値変更によるコストを指定する重み付け行列であり、
Ｍは、重み付け行列Ｑ３と併せて、ユーザに、前記アクチュエータ設定値プロファイルにおける異なる空間方向に対するコストを指定可能にする行列であり、
Ｑ４は、基準または目標設定値からずれたアクチュエータ設定値のコストを指定する重み付け行列であり、
Ｓは、重み付け行列Ｑ５と併せて、ユーザに、特定のアクチュエータ内設定方向にＣＤアクチュエータ・アレイおよびＭＤアクチュエータを移動するコストを指定可能にする行列であり、
Ａおよびｂは、ハード制約を指定する制約行列である、プロセス。
【請求項１３】
請求項１に記載のプロセスであって、各ＭＤアクチュエータは、１×１配列と見なされる、プロセス。
【請求項１４】
請求項１に記載のプロセスであって、前記シート測定データを処理し、複数のシート特性測定値配列を確立する前記ステップは、
シート特性におけるノイズおよび制御できないＭＤ変動を除去するために、時間フィルタを用いて前記シート特性測定データのフィルタリングを実施するステップを含む、プロセス。
【請求項１５】
請求項１に記載のプロセスであって、
シート特性のＭＤコンポーネントとＣＤコンポーネントの何れを制御すべきか指定するステップを、
追加して含むプロセス。
【請求項１６】
１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスであって、
規則的な間隔で複数のシート特性を測定し、シートの測定データを収集するステップと、
前記シート測定データを処理し、複数のシート特性測定値配列を確立するステップと、
前記シート特性測定値配列を共通解像度にマップするステップと、
前記共通解像度シート特性測定値配列を、１つのより大きい１次元共通解像度測定値配列に結合するステップと、
前記共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力することにより、シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
複数の未来のシート特性目標配列を１つの目標配列に結合するステップと、
前記シート製造プロセスの前記推定される現在の内部状態配列および前記シート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
前記シート特性の未来の予測配列、前記目標配列、オブジェクト関数の重み、最後のアクチュエータ設定値配列、およびハード制約を、オブジェクト関数に入力するステップと、
前記オブジェクト関数を解いて、前記シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含むプロセス。
【請求項１】
１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスであって、
規則的な間隔で複数のシート特性を測定し、シートの測定データを収集するステップと、
前記シート測定データを処理し、複数のシート特性測定値配列を確立するステップと、
前記シート特性測定値配列を処理し、１次元共通解像度測定値配列を確立するステップと、
前記シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
未来のシート特性目標配列を確立するステップと、
前記シート製造プロセスの前記推定される現在の内部状態配列およびシート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
前記シート特性の未来の予測配列、前記未来のシート特性目標配列、および前のアクチュエータ設定値配列を、可解なオブジェクト関数に入力し、前記シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含むプロセス。
【請求項２】
請求項１に記載のプロセスであって、前記シート特性測定値配列を処理し、１次元共通解像度測定値配列を確立する前記ステップは、
前記シート特性測定値配列を共通解像度にマップするステップと、
前記共通解像度シート特性測定値配列を、より大きい１次元共通解像度測定値配列に結合するステップと、
を含む、プロセス。
【請求項３】
請求項２に記載のプロセスであって、前記シート特性測定値配列を共通解像度にマップする前記ステップは、最も高いアクチュエータ解像度の３倍を超える前記共通解像度を選択するステップを含む、プロセス。
【請求項４】
請求項１に記載のプロセスであって、前記シート製造プロセスの前記推定される現在の内部状態配列を生成する前記ステップは、前記共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力するステップを含む、プロセス。
【請求項５】
請求項１に記載のプロセスであって、未来のシート特性目標配列を確立する前記ステップは、複数の未来のシート特性目標配列を１つの目標配列に結合するステップを含む、プロセス。
【請求項６】
請求項１に記載のプロセスであって、前記シート特性の未来の予測配列、前記未来のシート特性目標配列、および前記前のアクチュエータ設定値配列を、シート特性のオブジェクト関数に入力する前記ステップは、オブジェクト関数の重みおよびハード制約を入力するステップを含む、プロセス。
【請求項７】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、以下の状態空間形式（Ａ、Ｂ、Ｃ、Ｎｄ）；
ｘ（ｋ＋１）＝Ａｘ（ｋ）＋ＢΔｕ（ｋ−Ｎｄ）
ｙ（ｋ）＝Ｃｘ（ｋ）
で表され、
式中、
ｋは、サンプリング・インスタンスであり、
ｘは、前記プロセスの前記推定される現在の内部状態配列であり、
Δｕは、アクチュエータ設定値の過去の変更配列であり、
Ａは、前記プロセスの動的な時間情報を含む状態遷移行列であり、
Ｂは、前記プロセスの静的な空間情報を含む状態入力行列であり、
Ｃは、状態出力行列であり、
Ｎｄは、サンプル中のプロセス・トランスポート遅延である、プロセス。
【請求項８】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、インパルス応答モデルの形式で表される、プロセス。
【請求項９】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、ステップ応答モデルの形式で表される、プロセス。
【請求項１０】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、伝達関数モデルの形式で表される、プロセス。
【請求項１１】
請求項１に記載のプロセスであって、前記シート機プロセス・モデルは、２次元プロセス・モデルを識別するための自動化ツールを用いて生成される、プロセス。
【請求項１２】
請求項１に記載のプロセスであって、前記オブジェクト関数は、
【数１】
で表され、
式中、
ｅ（ｋ＋ｊ）＝ｙｒｅｆ（ｋ＋ｊ）−ｙｐ（ｋ＋ｊ）は、シート特性における未来の予測エラーであり、
Ｑ１は、前記シートの異なるシート特性間および異なるＣＤ位置間の相対的な重要性を指定する重み付け行列であり、
Ｑ２は、２つの連続するサンプル・インスタンス間における大きな前記アクチュエータ設定値変更によるコストを指定する重み付け行列であり、
Ｍは、重み付け行列Ｑ３と併せて、ユーザに、前記アクチュエータ設定値プロファイルにおける異なる空間方向に対するコストを指定可能にする行列であり、
Ｑ４は、基準または目標設定値からずれたアクチュエータ設定値のコストを指定する重み付け行列であり、
Ｓは、重み付け行列Ｑ５と併せて、ユーザに、特定のアクチュエータ内設定方向にＣＤアクチュエータ・アレイおよびＭＤアクチュエータを移動するコストを指定可能にする行列であり、
Ａおよびｂは、ハード制約を指定する制約行列である、プロセス。
【請求項１３】
請求項１に記載のプロセスであって、各ＭＤアクチュエータは、１×１配列と見なされる、プロセス。
【請求項１４】
請求項１に記載のプロセスであって、前記シート測定データを処理し、複数のシート特性測定値配列を確立する前記ステップは、
シート特性におけるノイズおよび制御できないＭＤ変動を除去するために、時間フィルタを用いて前記シート特性測定データのフィルタリングを実施するステップを含む、プロセス。
【請求項１５】
請求項１に記載のプロセスであって、
シート特性のＭＤコンポーネントとＣＤコンポーネントの何れを制御すべきか指定するステップを、
追加して含むプロセス。
【請求項１６】
１枚の材料シートを製造するためのシート製作機における縦方向ＭＤおよび横方向ＣＤアクチュエータの調整された制御のためのプロセスであって、
規則的な間隔で複数のシート特性を測定し、シートの測定データを収集するステップと、
前記シート測定データを処理し、複数のシート特性測定値配列を確立するステップと、
前記シート特性測定値配列を共通解像度にマップするステップと、
前記共通解像度シート特性測定値配列を、１つのより大きい１次元共通解像度測定値配列に結合するステップと、
前記共通解像度測定値配列およびアクチュエータ設定値の過去の変更配列を、シート機プロセス・モデル状態オブザーバに入力することにより、シート製造プロセスの推定される現在の内部状態配列を生成するステップと、
複数の未来のシート特性目標配列を１つの目標配列に結合するステップと、
前記シート製造プロセスの前記推定される現在の内部状態配列および前記シート機プロセス・モデルを用いて、シート特性の未来の予測配列を生成するステップと、
前記シート特性の未来の予測配列、前記目標配列、オブジェクト関数の重み、最後のアクチュエータ設定値配列、およびハード制約を、オブジェクト関数に入力するステップと、
前記オブジェクト関数を解いて、前記シート製作プロセスの調整されたＭＤおよびＣＤ制御のために、現在のアクチュエータ設定値の最適な変更配列を得るステップと、
を含むプロセス。
【図１】
【図２】
【図２】
【公表番号】特表２００６−５２５４４７（Ｐ２００６−５２５４４７Ａ）
【公表日】平成１８年１１月９日（２００６．１１．９）
【国際特許分類】
【出願番号】特願２００６−５１４２７６（Ｐ２００６−５１４２７６）
【出願日】平成１６年５月５日（２００４．５．５）
【国際出願番号】ＰＣＴ／ＵＳ２００４／０１３８４２
【国際公開番号】ＷＯ２００４／０９９８８６
【国際公開日】平成１６年１１月１８日（２００４．１１．１８）
【出願人】（５００５７５８２４）ハネウェル・インターナショナル・インコーポレーテッド (1,504)
【Ｆターム（参考）】
【公表日】平成１８年１１月９日（２００６．１１．９）
【国際特許分類】
【出願日】平成１６年５月５日（２００４．５．５）
【国際出願番号】ＰＣＴ／ＵＳ２００４／０１３８４２
【国際公開番号】ＷＯ２００４／０９９８８６
【国際公開日】平成１６年１１月１８日（２００４．１１．１８）
【出願人】（５００５７５８２４）ハネウェル・インターナショナル・インコーポレーテッド (1,504)
【Ｆターム（参考）】
[ Back to top ]