JPS6199512A - 圧延ロ−ルの偏心状態計測方法 - Google Patents
圧延ロ−ルの偏心状態計測方法Info
- Publication number
- JPS6199512A JPS6199512A JP59220127A JP22012784A JPS6199512A JP S6199512 A JPS6199512 A JP S6199512A JP 59220127 A JP59220127 A JP 59220127A JP 22012784 A JP22012784 A JP 22012784A JP S6199512 A JPS6199512 A JP S6199512A
- Authority
- JP
- Japan
- Prior art keywords
- eccentricity
- roll
- rolling
- load
- amplitude
- 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.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
- B21B37/58—Roll-force control; Roll-gap control
- B21B37/66—Roll eccentricity compensation systems
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Control Of Metal Rolling (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
[産業上の利用分野]
本発明は、圧延ロールの偏心状態を計測する方法に関釘
るらのである。
るらのである。
[従米扶術1
一般に、圧延での板厚制御では、圧延荷重と7−クロー
ル間の間隙(ワークロールギャップ)とが相関的に制御
される。
ル間の間隙(ワークロールギャップ)とが相関的に制御
される。
いわゆるスキッドマーク等に起因する板厚変動では、モ
ニタしている圧延荷重が大きくなった場、合、ワークロ
ールギャップを狭めるようにオートマチックデーノフン
トロール(A G C)が働く。
ニタしている圧延荷重が大きくなった場、合、ワークロ
ールギャップを狭めるようにオートマチックデーノフン
トロール(A G C)が働く。
ところで、このオートマチソクデーノコントロールによ
る板厚制御は、ユーザ側の厳しい精度要求に十分に適合
できないことから、M厚精度に寄与する他の要素が種々
検討され、重要な要素として圧延ロールの偏心が指摘さ
れている。
る板厚制御は、ユーザ側の厳しい精度要求に十分に適合
できないことから、M厚精度に寄与する他の要素が種々
検討され、重要な要素として圧延ロールの偏心が指摘さ
れている。
かかる圧墓σ−ルの偏心は、圧延中の熱変形。
偏摩耗やワークロールの再研磨時の真円からのずれ等に
よってもたらされる。0−ル偏心に起因する圧延荷重の
増大に対してはワークロールギャップを広げるように制
御する必要があり、前記入キッドマーク等による板厚変
動に起因する制御特性とは逆方向の制御特性が要求され
る。この場合、圧風ロール偏心に起因する圧延荷重の変
動は、スキッドマーク等による板厚変動に起因する圧延
荷重に重畳され、重畳された圧延荷重として検出される
ので、単純に圧延荷重の変動に伴なうオートマチックゲ
ージコントロールを実行すると、圧延ロール偏心に起因
する変動分は、ワークロールギャップの@11g1に関
して逆方向に加算されることになり。
よってもたらされる。0−ル偏心に起因する圧延荷重の
増大に対してはワークロールギャップを広げるように制
御する必要があり、前記入キッドマーク等による板厚変
動に起因する制御特性とは逆方向の制御特性が要求され
る。この場合、圧風ロール偏心に起因する圧延荷重の変
動は、スキッドマーク等による板厚変動に起因する圧延
荷重に重畳され、重畳された圧延荷重として検出される
ので、単純に圧延荷重の変動に伴なうオートマチックゲ
ージコントロールを実行すると、圧延ロール偏心に起因
する変動分は、ワークロールギャップの@11g1に関
して逆方向に加算されることになり。
11++1#の精度が大幅に低下する。
このため、圧延ロールの偏心バラノータ(偏心振幅、偏
心位相)を計測らしくは推定する手法が幾つか提案され
ている。
心位相)を計測らしくは推定する手法が幾つか提案され
ている。
その代表的な手法は、圧延ロールの偏心に起因する圧延
荷重変動が圧延ロールの回転に同調した比較的短い周期
(比較的高い周波WL)で現われるのに対し、スキッド
マーク等による板厚変動に伴なう圧延荷重変動は比較的
長い周期(比較的低い周波数)で現われることか呟これ
らが重畳して検出される圧延荷重変動を圧延ロールの1
回転かそれ以上の回転範囲にわたってサンプリングし、
それらのサンプリング値を周波数解析(7−リエ解析)
にかけて、高周波数成分として得られる個々丁
の圧延ロールの偏心に起因する圧
延面m変動を分離抽出する手法である。
荷重変動が圧延ロールの回転に同調した比較的短い周期
(比較的高い周波WL)で現われるのに対し、スキッド
マーク等による板厚変動に伴なう圧延荷重変動は比較的
長い周期(比較的低い周波数)で現われることか呟これ
らが重畳して検出される圧延荷重変動を圧延ロールの1
回転かそれ以上の回転範囲にわたってサンプリングし、
それらのサンプリング値を周波数解析(7−リエ解析)
にかけて、高周波数成分として得られる個々丁
の圧延ロールの偏心に起因する圧
延面m変動を分離抽出する手法である。
しかしながら、上記の7一リエ解析手法は、演n量が膨
大で計算時■0も長時間を要し、@厚のオンライン制御
中に上記の手法を実行することは実際上不可能である。
大で計算時■0も長時間を要し、@厚のオンライン制御
中に上記の手法を実行することは実際上不可能である。
かかる問題を改善するため、特公昭56−22281号
公報には、上、下のバック7ツブロールの偏心に着目し
、第1回目のフーリエ解析と第2回目の7−リエ解析と
を上下バック7ノプロールの位相をずらせた状態で行な
って、第1回目と第2回目の7−リエ解析により求めら
れた相対振幅から上下バック7γブロールの振幅を分離
抽出するようにしたらのが提案されている。
公報には、上、下のバック7ツブロールの偏心に着目し
、第1回目のフーリエ解析と第2回目の7−リエ解析と
を上下バック7ノプロールの位相をずらせた状態で行な
って、第1回目と第2回目の7−リエ解析により求めら
れた相対振幅から上下バック7γブロールの振幅を分離
抽出するようにしたらのが提案されている。
しかしながら、この方式は、7一リエ解析手法が有する
本質的な問題を何ら解決しておらず、上。
本質的な問題を何ら解決しておらず、上。
下バツクアツプロールの偏心しか推定できず、しかも計
測途中で上、下バツクアツプロールの位相を変更する必
要があるため、オフラインでしか実行できず、さらに非
能率でもある。
測途中で上、下バツクアツプロールの位相を変更する必
要があるため、オフラインでしか実行できず、さらに非
能率でもある。
一方、特公昭56−43803号公報には、圧延荷重又
はワークロールギャップを検出し、この検出出力を遅延
時間rだけ遅延させた信号とロール回転基本波信号、そ
の高調波信号および上下ロールの回転基本波信号の和か
ら求まるビート波信号との相関を上記遅延時間τを種々
変えて各信号毎に求め、この相関結果からロール偏心量
および偏心位相を求めるようにした偏心計測方式が提案
されている。
はワークロールギャップを検出し、この検出出力を遅延
時間rだけ遅延させた信号とロール回転基本波信号、そ
の高調波信号および上下ロールの回転基本波信号の和か
ら求まるビート波信号との相関を上記遅延時間τを種々
変えて各信号毎に求め、この相関結果からロール偏心量
および偏心位相を求めるようにした偏心計測方式が提案
されている。
この方式は、7−リエ解析法によらずに、圧延ロールの
偏心を検出することができる利点があるが、基準バック
アップミールの基本周波数とその高調波および上、下バ
ツクアツプロールのビートのみを考慮しているに過ぎな
いため、モデルが単純化され、求められた偏心量等が必
ずしも、実際の値を表わすものでIよないおそれがあり
、高い精度を期待し難い。
偏心を検出することができる利点があるが、基準バック
アップミールの基本周波数とその高調波および上、下バ
ツクアツプロールのビートのみを考慮しているに過ぎな
いため、モデルが単純化され、求められた偏心量等が必
ずしも、実際の値を表わすものでIよないおそれがあり
、高い精度を期待し難い。
さらに、上記2つの手法は、例えばバックアップミール
1回転中に得られろ多点データをもとに演算を行なうの
で、1回の制御動作と次の制御動作との間の時間間隔が
長くなり、連続したIIIIIlが行なえない基本的な
欠点がある。
1回転中に得られろ多点データをもとに演算を行なうの
で、1回の制御動作と次の制御動作との間の時間間隔が
長くなり、連続したIIIIIlが行なえない基本的な
欠点がある。
[発明の目的J
本発明の目的は、圧延荷重らしくはワークロールギャッ
プの1回の測定ごとに、圧延ロール個々に偏心振幅並び
に偏心位相を適確に算出することができる圧延ロールの
偏心状態計測方法を提供することである。
プの1回の測定ごとに、圧延ロール個々に偏心振幅並び
に偏心位相を適確に算出することができる圧延ロールの
偏心状態計測方法を提供することである。
[発明の構成J
本発明の基本的な特徴は、圧延機を構成する全ての圧延
ロールについて、各圧延ロールが固有の偏心振幅および
偏心位相を有するものとして、偏心ベクトルとでもいう
べきベクトル量を導入し、全ての圧延ロールの偏心の相
乗結果として現われる圧延荷重変動らしくはワークロー
ルギャップ変動が、上記ベクトルの各成分によって線形
に展開できることに着目し、このベクトル量を逐次形同
定刑に基づいて1回の計測データが得られるごとに逐次
推定を行なうようにしたことにある。ベクトルの各成分
が得られると、各圧延ロールの偏心振幅および偏心位相
は、簡単な演算によって求めることができる。
ロールについて、各圧延ロールが固有の偏心振幅および
偏心位相を有するものとして、偏心ベクトルとでもいう
べきベクトル量を導入し、全ての圧延ロールの偏心の相
乗結果として現われる圧延荷重変動らしくはワークロー
ルギャップ変動が、上記ベクトルの各成分によって線形
に展開できることに着目し、このベクトル量を逐次形同
定刑に基づいて1回の計測データが得られるごとに逐次
推定を行なうようにしたことにある。ベクトルの各成分
が得られると、各圧延ロールの偏心振幅および偏心位相
は、簡単な演算によって求めることができる。
L発明の効果J
本発明によれば、1回の測定ごとに圧延ロールの偏心パ
ラメータ(偏心振幅、偏心位相)が得られるので、オン
ラインWill 11時には、フィードフォワード制御
を連続的に応答性よく行なうことができる。
ラメータ(偏心振幅、偏心位相)が得られるので、オン
ラインWill 11時には、フィードフォワード制御
を連続的に応答性よく行なうことができる。
また、本発明によれば、逐大形同定則を用いているので
、得られたロール偏心パラメータの推定値の精度を同時
に計算できる。このため、圧延運転開始前のキスロール
運転時に行なうオフラインのロール偏心パラメータ計測
において、所望の精度を得るのに必要な測定値の個数を
判定できる。
、得られたロール偏心パラメータの推定値の精度を同時
に計算できる。このため、圧延運転開始前のキスロール
運転時に行なうオフラインのロール偏心パラメータ計測
において、所望の精度を得るのに必要な測定値の個数を
判定できる。
さらに、かかる精度の収れん性を判定することによって
、採用したロール偏心モデルの妥当性が評価でき、必要
に応じて偏心モデルを切り換えることができる。
、採用したロール偏心モデルの妥当性が評価でき、必要
に応じて偏心モデルを切り換えることができる。
[実施例]
以下、添付の図面を参照して本発明の実施例を詳細に説
明する。
明する。
l 第1図1ユ示オよう
、ユ、オア施9□1よ、ワー、。
、ユ、オア施9□1よ、ワー、。
−ル1,2を上、下バックアップロール3,4で支持し
た圧延は5について考察する。
た圧延は5について考察する。
上バツクアップロール3に対しては、圧延荷重Pを検出
するためのロードセル6が設けられ、下バツクアツプロ
ール4に対しては、ワークミール1.2のa−ルギャッ
プを調節するためのロールギャップ調整手段7が設けら
れている。*た、上。
するためのロードセル6が設けられ、下バツクアツプロ
ール4に対しては、ワークミール1.2のa−ルギャッ
プを調節するためのロールギャップ調整手段7が設けら
れている。*た、上。
下ワークロールおよび上、下バック77ブa−ル1.2
,3.4の回転角θ1.θ21 θ、lθ、は、パルス
ジェネレータ8,9,10.11によって夫々検出され
るようになっている。
,3.4の回転角θ1.θ21 θ、lθ、は、パルス
ジェネレータ8,9,10.11によって夫々検出され
るようになっている。
上記ロードセル6によって検出される圧延荷重Pおより
上、下ワークロールおよび上、下バックアップロール1
.2,3.4の回転角θ3.θ21θ1.θ6、さらに
、ロールギャップll4g1手段7に対する制御出力へ
Scは、コンピュータに入力され、コンピュータは、こ
れら入力データに基づいて以下の処理を実行する。
上、下ワークロールおよび上、下バックアップロール1
.2,3.4の回転角θ3.θ21θ1.θ6、さらに
、ロールギャップll4g1手段7に対する制御出力へ
Scは、コンピュータに入力され、コンピュータは、こ
れら入力データに基づいて以下の処理を実行する。
〈データ処理〉
N] 信号処理ルーチン
データの読取りは、所定のサンプリングタイミング(こ
のタイミングとしては、例えば上バツク7ノプロール3
のパルスジェネレータ10のパルスを用いることができ
る。)ごとに行ない、読取ったデータについて、例えば
圧延荷重信号Pはその直を成分を除去する等の前処理を
行なう5いま、一般に、圧延機による圧延において、ワ
ークロールギャップが△Sだけ変化したときの圧延荷重
の変化1八Pは、 但し、M・・・圧延機の等価ミル定数 Q・・・被圧延材の塑性係数 となる、ただし、圧延運転前におけるキスロール時には
Q=■なので、MQ/(M+Q)=Mである。
のタイミングとしては、例えば上バツク7ノプロール3
のパルスジェネレータ10のパルスを用いることができ
る。)ごとに行ない、読取ったデータについて、例えば
圧延荷重信号Pはその直を成分を除去する等の前処理を
行なう5いま、一般に、圧延機による圧延において、ワ
ークロールギャップが△Sだけ変化したときの圧延荷重
の変化1八Pは、 但し、M・・・圧延機の等価ミル定数 Q・・・被圧延材の塑性係数 となる、ただし、圧延運転前におけるキスロール時には
Q=■なので、MQ/(M+Q)=Mである。
ロールギャップ変化量ΔSは、制御によるロールギャッ
プ変化成分ΔScと、各圧延ロールの偏心成分が重畳さ
れたロール偏心成分ΔSeとの合成量である。
プ変化成分ΔScと、各圧延ロールの偏心成分が重畳さ
れたロール偏心成分ΔSeとの合成量である。
ΔS=ΔSc−ΔSe −(2)
圧延荷重の変化量ΔPは、ロードセル6の出力から直流
成分を取り除いた信号の値より求まり、制御によるロー
ルギャップ変化成分ΔScは、ロールギャップ調整手段
71÷対する制御出力から検出することができる。
成分を取り除いた信号の値より求まり、制御によるロー
ルギャップ変化成分ΔScは、ロールギャップ調整手段
71÷対する制御出力から検出することができる。
したがって、圧延荷重変化ΔPのうち、ロール偏心に起
因する成分ΔPeは、(1)、(2)式より以下の通り
に求めることができる。
因する成分ΔPeは、(1)、(2)式より以下の通り
に求めることができる。
以上の処理は、信号処理ルーチン[+1にお〜1て行な
われる。
われる。
[ul 偏心パラメータの推定演算
いま、簡単のため、各ロールの偏心について高調波成分
を無視し、各ロールの偏心をロールの基本周波数θ;の
正弦波で近似すると、ロール偏心によるロールギャップ
変化ΔSeは、各ロールの偏心の重ね合せとして、次式
で表わすことができる。
を無視し、各ロールの偏心をロールの基本周波数θ;の
正弦波で近似すると、ロール偏心によるロールギャップ
変化ΔSeは、各ロールの偏心の重ね合せとして、次式
で表わすことができる。
ΔSe =、X Ai・5in(θ1+φi) ・
・・44)+=1 但し、 1 ・・ ロール番号 Ai・・・各ロールの偏心振幅 φ1・・ 各ロールの偏心位相 θ1 ・ 各ロールの回転角度 上記の式(4)において、各ロールの回転角θ1は、各
パルスノエネレータ8〜11の出力から既知であり、そ
の他、各ロールの角速度ωIを連続的に測定し、これを
時間的に積分することによっても求めることができる。
・・44)+=1 但し、 1 ・・ ロール番号 Ai・・・各ロールの偏心振幅 φ1・・ 各ロールの偏心位相 θ1 ・ 各ロールの回転角度 上記の式(4)において、各ロールの回転角θ1は、各
パルスノエネレータ8〜11の出力から既知であり、そ
の他、各ロールの角速度ωIを連続的に測定し、これを
時間的に積分することによっても求めることができる。
いま、k番目のサンプリングパルスに同期して得られた
ロール偏心による荷重変化分をΔPe(k)。
ロール偏心による荷重変化分をΔPe(k)。
各ロールの回転角をθ1(k)とすると、(4)式より
、i + n(
k)+ A1・SInφi ・eosθ1(k)]++
n(k) ・・・・・(6)と展開さ
れる。ここで、++(k)はΔPe(k)を測定する際
の観測雑音(測定誤差)である。
、i + n(
k)+ A1・SInφi ・eosθ1(k)]++
n(k) ・・・・・(6)と展開さ
れる。ここで、++(k)はΔPe(k)を測定する際
の観測雑音(測定誤差)である。
いま、既知の量である各ロールの回転角θ1(k)の正
弦および余弦をベクトル成分とするベクトルa(k)と
、未知量である偏心パラメータ(AI、φ1)の関数で
あられされる成分を有するベタレレマ(以下では、偏心
ベクトルという)とを以下の如く定義する。
弦および余弦をベクトル成分とするベクトルa(k)と
、未知量である偏心パラメータ(AI、φ1)の関数で
あられされる成分を有するベタレレマ(以下では、偏心
ベクトルという)とを以下の如く定義する。
上記(7)、(8)式で夫々定義されるベクトルを用い
ると、(6)式は以下の如くに表わされる。
ると、(6)式は以下の如くに表わされる。
ΔPe(k) = 7”(k) ・マ+n(k)
・−(9!この式(9)は、既知のベクトルaと未知の
偏心ベクトルマとの内積でロール偏心による圧延荷重成
分ΔPe(k)が表わされること、換言すれば、偏心ベ
クトル×の成分xiの線形関数で表わされることを意味
する。
・−(9!この式(9)は、既知のベクトルaと未知の
偏心ベクトルマとの内積でロール偏心による圧延荷重成
分ΔPe(k)が表わされること、換言すれば、偏心ベ
クトル×の成分xiの線形関数で表わされることを意味
する。
したがって、上記(9)式における偏心ベクトルXの推
定には逐次形量定刻を適用することができる。
定には逐次形量定刻を適用することができる。
いま、k番目のサンプリングパルスと同期した△
測定結果に基づく偏心ベクトルマの推定値マ(k)は、
逐次形量定刻にしたがえば、以下の式によって求めるこ
とができる。
逐次形量定刻にしたがえば、以下の式によって求めるこ
とができる。
T △
+7(kL[ΔPe(k)−1(k)・ マ(k−1)
]・・・・・・(10) 但し、ベクトルマ(k)は、ゲインベクトルと呼ばれ、 で与えられる。
]・・・・・・(10) 但し、ベクトルマ(k)は、ゲインベクトルと呼ばれ、 で与えられる。
また、B (k)は推定精度を表わす共分散行列であり
、次式で与えられる。
、次式で与えられる。
B(k)=−[B(k−11
A1
・・・・・(12)
なお、7(0)は推定値の初期値として、また共分散行
列B(0)は推定値の共分散行列の初期値として、夫々
適当に与えることができる。
列B(0)は推定値の共分散行列の初期値として、夫々
適当に与えることができる。
また、パラメータλ。、λ、lλ2は推定機能を14節
するためのパラメータであり、[λ。=λ、=λ、=I
Jとすれば、逐次形最小2水法による推定となり、[^
、=11 λ。=1/λ、=σe2(σe2は観測雑音
n(k)のパワー)] とすれば、カルマン・フィルタ
による推定となる。
するためのパラメータであり、[λ。=λ、=λ、=I
Jとすれば、逐次形最小2水法による推定となり、[^
、=11 λ。=1/λ、=σe2(σe2は観測雑音
n(k)のパワー)] とすれば、カルマン・フィルタ
による推定となる。
また、[λ2=11 λ。=λ、=γ(0,99S <
7く1)1と設定すれば、重み付き最小2乗法による推
定となり、ロールの熱膨張等による偏心状態の変化に対
する追従速度を早めることができる。
7く1)1と設定すれば、重み付き最小2乗法による推
定となり、ロールの熱膨張等による偏心状態の変化に対
する追従速度を早めることができる。
上記の演算;こより、偏心ベクトルh各成分xIの推定
値が求まると、各ロールの偏心振幅AI。
値が求まると、各ロールの偏心振幅AI。
偏心位相φ1は夫々以下の如く求めることができる。
Ai =r了第−W ・・・・・(13)以上のよう
にして、偏心パラメータが各ロールごとに推定されるこ
とになる。かがる推定演算は、偏心パラメータり推定演
算ルーチン[11]において行なわれる。
にして、偏心パラメータが各ロールごとに推定されるこ
とになる。かがる推定演算は、偏心パラメータり推定演
算ルーチン[11]において行なわれる。
[+111推定精度・モデル判定用
この判別は、上記共分散行列B (k)の性質を利用し
て、推定の精度や推定に用いたモデルの妥当(性を特定
する。
て、推定の精度や推定に用いたモデルの妥当(性を特定
する。
上記共分散行列B (k)の対角成分の和trB(k)
は、次の関係を満足することが知られている。
は、次の関係を満足することが知られている。
LrB(k)oCE [If ” (k)−7(k)
II l −−(15)但し、E([)は、「の期待値
を表わす。
II l −−(15)但し、E([)は、「の期待値
を表わす。
つまり、trB(k)は、(i心ベクトル7の推定値へ
7(k)の精度(分散)に対応し、これによって現在の
推定精度を確認することができる。
推定精度を確認することができる。
また、一般に、圧延荷重のす/プル数(k)が増加する
と、trB(k)は小さくなるので、このLrB (i
t)に対して予め許容誤差を適当に設定しておけば、許
容誤差以下になるサンプル数が判定でき、例えばキスロ
ール状態で行なうオフライ/時の偏心パラメータ推定で
は、許容誤差以下となるサンプル数に達した時点でサン
プリングを終了するようにすれば、オフラインでの計測
が短時間で無駄なく行なえる。
と、trB(k)は小さくなるので、このLrB (i
t)に対して予め許容誤差を適当に設定しておけば、許
容誤差以下になるサンプル数が判定でき、例えばキスロ
ール状態で行なうオフライ/時の偏心パラメータ推定で
は、許容誤差以下となるサンプル数に達した時点でサン
プリングを終了するようにすれば、オフラインでの計測
が短時間で無駄なく行なえる。
一方、上記共分散行列B (k)の対角成分和しrBf
k)が許容誤差以下に収れんしない場合がありうる。こ
れは、例えば、ロール中に全く間−の半径を有するロー
ルが2個以上あるような場合に生じうる。つまり、この
場合には、これら同一径のロールから生ずる偏心は全く
同じ周波数を持っているので圧延荷重信号からこれらの
偏心パラメータを分離して計測することが不可能となる
。
k)が許容誤差以下に収れんしない場合がありうる。こ
れは、例えば、ロール中に全く間−の半径を有するロー
ルが2個以上あるような場合に生じうる。つまり、この
場合には、これら同一径のロールから生ずる偏心は全く
同じ周波数を持っているので圧延荷重信号からこれらの
偏心パラメータを分離して計測することが不可能となる
。
実際に、trB(k)が小さい値に収れんしない場合に
は、キスロール状態での推定の時、慨において、同一径
のロールの偏心を合成した1つのロール偏心で近似する
縮小モデルに変更する。
は、キスロール状態での推定の時、慨において、同一径
のロールの偏心を合成した1つのロール偏心で近似する
縮小モデルに変更する。
上記の判定処理は、実際には材料の圧延を開始する前に
以下の1序で行なう。
以下の1序で行なう。
まず、ステップ(1)−1では、サンプリング回数kが
予め設定した回fiKに達したときに、tr13(K)
を予め設定したモデル判定用しきい値ε、と比較し、L
rB(K)>ε1のときは、収れん性がないものとして
、合成モデル(以下に詳述する)に切り替える。
予め設定した回fiKに達したときに、tr13(K)
を予め設定したモデル判定用しきい値ε、と比較し、L
rB(K)>ε1のときは、収れん性がないものとして
、合成モデル(以下に詳述する)に切り替える。
trB(K)≦C1のときは、最初に設定した初期モデ
ルのままに推定を続行する。
ルのままに推定を続行する。
また、ステップ(III)−2では、推定演算ごとに、
LrB(k)を予め設定した精度判定用のしきい値ε2
(εハε、)と比較し、しきい値ε2に達していない場
合には[1]、[口lと同様の処理による推定演算を続
行し、順次に推定値を更新する。
LrB(k)を予め設定した精度判定用のしきい値ε2
(εハε、)と比較し、しきい値ε2に達していない場
合には[1]、[口lと同様の処理による推定演算を続
行し、順次に推定値を更新する。
しきい値ε2以下にcrB(k)が低下したときは、推
定値の精度が十分であると特定して、キスロール状態で
のオフライン推定演算を終了する。
定値の精度が十分であると特定して、キスロール状態で
のオフライン推定演算を終了する。
[IVl 合成モデルによる偏心パラメータ演算例え
ば、ワークロール1. 2のロール径が全く同一の場合
には、これら2つのワーククールの偏心を合成して1つ
のロールの偏心に置換する。
ば、ワークロール1. 2のロール径が全く同一の場合
には、これら2つのワーククールの偏心を合成して1つ
のロールの偏心に置換する。
つまり、2つのロールの偏心量の和を1つの偏心量Al
25:n (θ1+φi)で置換したうえで、上記と同
様にして合成偏心パラメータを演算する。
25:n (θ1+φi)で置換したうえで、上記と同
様にして合成偏心パラメータを演算する。
[V] 偏心パラメータの確定
圧延運転開始前に、上記N]〜[IV]の処理によって
求めら八たオフライン推定値は、これを各o −ルの偏
心バラ/−夕としてコンビコータの所定アドレスにメモ
リされる。
求めら八たオフライン推定値は、これを各o −ルの偏
心バラ/−夕としてコンビコータの所定アドレスにメモ
リされる。
〈オンライン偏心制御〉
次に、オンライン偏心制御方式の一例を説明する。
ロール偏心状態が圧延中において不変であhば、オフラ
インで求めておいた偏心のパラメータを用いたフィード
フォワードによって、ロール偏心を補償するようなロー
ルI!1度制御を行なえばよいが、現実には、ロールの
熱膨張等によってロール偏心は変化するので、オンライ
ンでの修正が必要となる。
インで求めておいた偏心のパラメータを用いたフィード
フォワードによって、ロール偏心を補償するようなロー
ルI!1度制御を行なえばよいが、現実には、ロールの
熱膨張等によってロール偏心は変化するので、オンライ
ンでの修正が必要となる。
本推定方法を用いたオンライン制御には、始めに述べた
ように、連応性と多周波成分への披張性という特徴があ
るが、ここでは従来一般に行われている制御方式に準じ
た制御方法について説明する。
ように、連応性と多周波成分への披張性という特徴があ
るが、ここでは従来一般に行われている制御方式に準じ
た制御方法について説明する。
各時点にでのロール偏心によるロール間FA変動は、5
in(θ1(k)+φi off) −(16)(と
推定できる。
in(θ1(k)+φi off) −(16)(と
推定できる。
一方、オンラインでの推定は、ロール偏心のモデルとし
ては、Hi−ば基準バラ2フ/ブ0−ル3の周波数の正
弦波より成るモデルとし、AsとφSのみを推定すれば
よいような形の簡略モデルを用いる。パラメータの推定
は、基準ロール1回転分のデータを用いて第1図中の(
1)、([1)の処理と同様になされる。得られた推定
値は、次の1回転onとすると、次の1@転の闇の時、
慨kにおける、ヘ ロール偏心によるロール間隙変動の推定値ΔSeとなる
。
ては、Hi−ば基準バラ2フ/ブ0−ル3の周波数の正
弦波より成るモデルとし、AsとφSのみを推定すれば
よいような形の簡略モデルを用いる。パラメータの推定
は、基準ロール1回転分のデータを用いて第1図中の(
1)、([1)の処理と同様になされる。得られた推定
値は、次の1回転onとすると、次の1@転の闇の時、
慨kにおける、ヘ ロール偏心によるロール間隙変動の推定値ΔSeとなる
。
前に述べたように、ロール偏心状態の経時変化に対処す
るためには、オンライン方式は必要不可欠なものである
。しかし、オンライン方式として上記で述べた簡略モデ
ルによるものを泪いようとすると、各ロールの周波数成
分の重ね合せにより生ずるビート現象を表現するために
オフラインで求めた厳密なモデルによる推定の結果を用
いて補なう必要がある。そのため、制御には一般に、両
但し、0≦a≦1 が用いられる。
るためには、オンライン方式は必要不可欠なものである
。しかし、オンライン方式として上記で述べた簡略モデ
ルによるものを泪いようとすると、各ロールの周波数成
分の重ね合せにより生ずるビート現象を表現するために
オフラインで求めた厳密なモデルによる推定の結果を用
いて補なう必要がある。そのため、制御には一般に、両
但し、0≦a≦1 が用いられる。
ハ
Δ5e(k)を打ち消すために、時点にでのo −小開
度Δ5c(k)が △ Δ5c(k) ” Δ5e(k) −(1
9)となるようにフィード7fワード制御を行なう。
度Δ5c(k)が △ Δ5c(k) ” Δ5e(k) −(1
9)となるようにフィード7fワード制御を行なう。
〈実施結果〉
次に、本発明方法にしたがった実施結果について説明す
る。
る。
この実験に用いた各ロールの偏心振幅AI、偏心位相φ
1およびロール半径Riは表に示す通りであり、圧延速
度は300 (z/分)であった、上バツクアップロー
ル1回転ごとにパルスノエネレータ10からは、90個
のパルスが発せられ、このパルスをサンプリングタイミ
ングとして約1011 Seeごとに圧延荷重の測定を
行なった。なお、観測雑音のパワーは2(Lon’)と
仮定した。
1およびロール半径Riは表に示す通りであり、圧延速
度は300 (z/分)であった、上バツクアップロー
ル1回転ごとにパルスノエネレータ10からは、90個
のパルスが発せられ、このパルスをサンプリングタイミ
ングとして約1011 Seeごとに圧延荷重の測定を
行なった。なお、観測雑音のパワーは2(Lon’)と
仮定した。
第2図には、偏心による圧延荷重変化ΔPeの状態を示
す。
す。
この圧延荷重変化ΔPeをらとに、^。=2゜λ1=1
.λ、=O,Sとしたカルマンフィルタによって偏心パ
ラメータの逐次推定を行なった。
.λ、=O,Sとしたカルマンフィルタによって偏心パ
ラメータの逐次推定を行なった。
第3図には推定演算によって求められた各ロールの偏心
振幅を、また@4図には各ロールの偏心位相を夫々示す
。
振幅を、また@4図には各ロールの偏心位相を夫々示す
。
上記表と比較すれば明らかなように、推定値は真の値に
ごく短期間で収れんすることが分る。
ごく短期間で収れんすることが分る。
また、tjtJ5図からは、推定精度に関係するtrB
(k)も早期に10−5以下の値に達することが分る。
(k)も早期に10−5以下の値に達することが分る。
さらに、第6図、第7図には、(18)式において、α
=0(つまり、簡略モデルのみによる場合)と、α=1
(厳密モデルによるオフライン推定のみによる場合)
の圧延荷重変動ΔPeを示す。
=0(つまり、簡略モデルのみによる場合)と、α=1
(厳密モデルによるオフライン推定のみによる場合)
の圧延荷重変動ΔPeを示す。
簡略モデルでは、第6図に示すように、ロール偏心の影
響が十分除去されていないのに対し、第7図の場合には
、ロール偏心の経時変化が無い状態での実施例なので、
ロール偏心の影響がほぼ除去されている。
響が十分除去されていないのに対し、第7図の場合には
、ロール偏心の経時変化が無い状態での実施例なので、
ロール偏心の影響がほぼ除去されている。
以上の実施例では、簡単のためにロール偏心の、
imp″″”11°5”t、 *a*+:”
”61”埋から明らかなように、この種の高調波成分を
容易に推定演算に組み入れることができる。
imp″″”11°5”t、 *a*+:”
”61”埋から明らかなように、この種の高調波成分を
容易に推定演算に組み入れることができる。
また、上記の実施例では、4段の圧延機につし・で説明
したが、段数に制限されることはなく、より多段の圧延
機における各ロールの偏心パラメータを求めることがで
きることは〜1うまでもない。
したが、段数に制限されることはなく、より多段の圧延
機における各ロールの偏心パラメータを求めることがで
きることは〜1うまでもない。
さらに、ロール偏心に対するオンライン制御に関しては
、前記簡略モデルに限定されるものではなく、オンライ
ンで厳密モデルを用いて本発明方法にかかる偏心パラメ
ータの推定演算を行なって、求められた偏心パラ/−夕
に基づいて偏心制御を実行しうろことはいうまでもない
。
、前記簡略モデルに限定されるものではなく、オンライ
ンで厳密モデルを用いて本発明方法にかかる偏心パラメ
ータの推定演算を行なって、求められた偏心パラ/−夕
に基づいて偏心制御を実行しうろことはいうまでもない
。
第1図は本発明にがかるa−小偏心状態計測システムを
示す構成図、第2図〜第70は各々本発明にかかる推定
方法にしたがった計算結果を示すグラフで、第2図はロ
ール偏心に起因する圧延荷重変動分ΔPeを示し、第3
図、第4図は夫々推定方法によって得られるロール偏心
振幅および偏心位相の収れん状態を示し、Pt55図は
推定精度を表わすLrBの変化を示し、第6図、第7図
は夫々ロール偏心制御を、簡略モデルのみを用l・て行
なつた場合と、ff1ffiモデルのオフライン推定値
のみを用いて行なった場合を示す。
示す構成図、第2図〜第70は各々本発明にかかる推定
方法にしたがった計算結果を示すグラフで、第2図はロ
ール偏心に起因する圧延荷重変動分ΔPeを示し、第3
図、第4図は夫々推定方法によって得られるロール偏心
振幅および偏心位相の収れん状態を示し、Pt55図は
推定精度を表わすLrBの変化を示し、第6図、第7図
は夫々ロール偏心制御を、簡略モデルのみを用l・て行
なつた場合と、ff1ffiモデルのオフライン推定値
のみを用いて行なった場合を示す。
Claims (1)
- (1)圧延ロールの偏心に起因する圧延荷重の変動分を
検出し、該変動分を、各圧延ロールの偏心振幅および偏
心位相で表わされるベクトル成分で線形に展開し、これ
らベクトル成分を所定のタイミングごとに逐次形同定則
を用いて逐次推定して各圧延ロールの偏心振幅および偏
心位相を演算するようにしたことを特徴とする圧延ロー
ルの偏心状態計測方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP59220127A JPS6199512A (ja) | 1984-10-18 | 1984-10-18 | 圧延ロ−ルの偏心状態計測方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP59220127A JPS6199512A (ja) | 1984-10-18 | 1984-10-18 | 圧延ロ−ルの偏心状態計測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS6199512A true JPS6199512A (ja) | 1986-05-17 |
JPS649087B2 JPS649087B2 (ja) | 1989-02-16 |
Family
ID=16746329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP59220127A Granted JPS6199512A (ja) | 1984-10-18 | 1984-10-18 | 圧延ロ−ルの偏心状態計測方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS6199512A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998024567A1 (de) * | 1996-12-04 | 1998-06-11 | Voest-Alpine Industrieanlagenbau Gmbh | Verfahren zur kompensation der exzentrizität der stütz- und/oder arbeitswalzen in einem duo- oder quarto-walzgerüst |
JP2015000419A (ja) * | 2013-06-17 | 2015-01-05 | 株式会社神戸製鋼所 | 板クラウン予測モデルの決定方法 |
EP2644288A4 (en) * | 2010-11-22 | 2015-07-22 | Toshiba Mitsubishi Elec Inc | ROLLER CONTROL DEVICE |
-
1984
- 1984-10-18 JP JP59220127A patent/JPS6199512A/ja active Granted
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998024567A1 (de) * | 1996-12-04 | 1998-06-11 | Voest-Alpine Industrieanlagenbau Gmbh | Verfahren zur kompensation der exzentrizität der stütz- und/oder arbeitswalzen in einem duo- oder quarto-walzgerüst |
EP2644288A4 (en) * | 2010-11-22 | 2015-07-22 | Toshiba Mitsubishi Elec Inc | ROLLER CONTROL DEVICE |
JP2015000419A (ja) * | 2013-06-17 | 2015-01-05 | 株式会社神戸製鋼所 | 板クラウン予測モデルの決定方法 |
Also Published As
Publication number | Publication date |
---|---|
JPS649087B2 (ja) | 1989-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4126027A (en) | Method and apparatus for eccentricity correction in a rolling mill | |
KR900003970B1 (ko) | 다중 압연기의 롤 편심 제거 제어 방법 및 장치 | |
JP2001198612A (ja) | 周期的ゲージ偏差補正システム | |
US4222254A (en) | Gauge control using estimate of roll eccentricity | |
EP0015866B1 (en) | Method of controlling roll eccentricity of rolling mill and apparatus for performing the same method | |
JP2000288614A (ja) | 圧延機の板厚制御装置 | |
JPS6199512A (ja) | 圧延ロ−ルの偏心状態計測方法 | |
KR100397630B1 (ko) | 압연플랜트의 롤 편심제어방법 | |
CN114074120B (zh) | 设备控制装置、轧制控制装置、设备控制方法以及计算机可读记录介质 | |
CN111036686B (zh) | 一种冷轧机支撑辊偏心的预估补偿方法 | |
US4763273A (en) | Apparatus for detecting eccentricity of roll in rolling mill | |
US3892112A (en) | Rolling mill gauge control | |
JP2710863B2 (ja) | 圧延機の板厚制御方法 | |
JPH07214133A (ja) | 熱延鋼帯の捲取温度制御方法 | |
EP0992295A2 (en) | Method and device for the active compensation of periodic disturbances during hot or cold rolling | |
KR100848651B1 (ko) | 압연기의 상부 및 하부 편심 분리 검출 방법 | |
JPS59209413A (ja) | 自動板厚制御におけるゲインの決定方法 | |
US20050155402A1 (en) | Method and apparatus for detecting roll eccentricity utilizing pulse generator in rolling mill | |
JPS6216811A (ja) | ロ−ル偏心制御の位相遅れ補正方法 | |
JPS5897417A (ja) | ロ−ル偏心制御装置 | |
Ginzburg | CHA? TER12 | |
JPH06292916A (ja) | 厚鋼板圧延における板厚制御法 | |
SU1346288A1 (ru) | Устройство дл компенсации эксцентриситета валков прокатной клети | |
Lee et al. | Fault diagnosis of roll shape under the speed change in hot rolling mill | |
SU770588A1 (ru) | Способ устранения влияния эксцентриситета опорных валков и устройство для его осуществления 1 |