[go: up one dir, main page]

JP2009208115A - Method and device for calculating parameter of rolling control, and rolling simulation device - Google Patents

Method and device for calculating parameter of rolling control, and rolling simulation device Download PDF

Info

Publication number
JP2009208115A
JP2009208115A JP2008053420A JP2008053420A JP2009208115A JP 2009208115 A JP2009208115 A JP 2009208115A JP 2008053420 A JP2008053420 A JP 2008053420A JP 2008053420 A JP2008053420 A JP 2008053420A JP 2009208115 A JP2009208115 A JP 2009208115A
Authority
JP
Japan
Prior art keywords
rolling
parameter
stand
control
parameters
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.)
Pending
Application number
JP2008053420A
Other languages
Japanese (ja)
Inventor
Akira Murakami
晃 村上
Makiyuki Nakayama
万希志 中山
Kenichi Sano
研一 佐野
Mitsuo Okamoto
充央 岡本
Kaoru Shigeno
馨 重野
Kazuaki Tanebe
和明 種部
Shigeto Koizumi
重人 小泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kobe Steel Ltd
Original Assignee
Kobe Steel Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kobe Steel Ltd filed Critical Kobe Steel Ltd
Priority to JP2008053420A priority Critical patent/JP2009208115A/en
Publication of JP2009208115A publication Critical patent/JP2009208115A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Control Of Metal Rolling (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To calculate optimally adjusted parameters of rolling control by identifying friction coefficients of the whole of a rolling machine and by performing, based on the obtained friction coefficients, a simulation to reproduce a rolling state. <P>SOLUTION: A dynamic model of rolling is previously prepared. Concerning a temporal change of rolling loads or roll gaps at a plurality of rolling stands when accelerating or decelerating a rolling material, the difference between a calculated value of the temporal change, which is obtained from the dynamic model of rolling, and an actual value of the temporal change is made minimum by optimizing a friction coefficient parameter or a deformation resistance parameter at each rolling stand. Using the dynamic model of rolling into which these parameters are incorporated, parameters of rolling control are calculated. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

本発明は、複数の圧延スタンドを備えた圧延機に対して、最適化された圧延制御パラメータを算出する方法及び算出する装置、ならびに、最適化された圧延制御パラメータを算出可能な圧延シミュレーション装置に関する。   The present invention relates to a method and a device for calculating an optimized rolling control parameter for a rolling mill having a plurality of rolling stands, and a rolling simulation device capable of calculating an optimized rolling control parameter. .

圧延機を制御する制御パラメータ(圧延制御パラメータ)に関し、最適に調整されたパラメータを算出する方法としては、例えば、特許文献1に開示されたものがある。
特許文献1に開示されたシミュレーション装置は、圧延機などの各種プラントをモデル化したモデル手段と、各種プラントを制御する制御装置をモデル化したブロック要素手段とを備えており、このシミュレーション装置を用いて、各種プラントのパラメータを変更するものとなっている。前述のモデル手段におけるパラメータの変更については、誤差信号が許容誤差以下になるように、プラントに精通した者がパラメータを変更するか、推論手段を設けてパラメータを変更するか、同定機構によりパラメータを変更することとしている。
As a method for calculating an optimally adjusted parameter regarding a control parameter (rolling control parameter) for controlling the rolling mill, for example, there is one disclosed in Patent Document 1.
The simulation apparatus disclosed in Patent Document 1 includes a model unit that models various plants such as a rolling mill, and a block element unit that models a control device that controls the various plants. The parameters of various plants are changed. Regarding the change of the parameter in the above-mentioned model means, a person familiar with the plant changes the parameter so that the error signal is less than the allowable error, or the inference means is provided to change the parameter, or the parameter is changed by the identification mechanism. I am going to change it.

一方、特許文献2には、単圧延スタンドからなる圧延機に対するロールギャップ補正装置が開示されている。この装置では、摩擦係数の変化による圧延スタンド出側の板厚変化を修正するにあたり、ニューラルネットワークを用いて摩擦係数のモデル式を推定し、この推定式により得られた値を用いて圧延スタンドのロールギャップ設定値(制御パラメータ)を補正している。
特開昭63−239508号公報 特開平7−246411号公報
On the other hand, Patent Document 2 discloses a roll gap correction device for a rolling mill composed of a single rolling stand. In this equipment, when correcting the change in the thickness of the rolling stand exit side due to the change in the friction coefficient, a model formula of the friction coefficient is estimated using a neural network, and the value obtained from this estimation formula is used to estimate the rolling stand. The roll gap setting value (control parameter) is corrected.
JP 63-239508 A JP-A-7-246411

しかしながら、複数の圧延スタンドを備えた圧延機で圧延材を加速又は減速する際の制御に、特許文献1,2の技術を適用するには、以下に述べるような課題が存在する。
例えば、特許文献1の技術を圧延材の加減速時に適用した場合、プラントに精通した者が有している「定性的な表現によるエキスパートな知識」に基づくパラメータ変更や推論手段によるパラメータ変更では、特許文献1にも記載のように、曖昧な情報であり、正確なプラント制御ができないことは明らかである。また、同定機構によりパラメータを変更する方法についても、本技術は線形系に限られ非線形系の同定は不可能である。
However, in order to apply the techniques of Patent Documents 1 and 2 to control when accelerating or decelerating a rolled material with a rolling mill having a plurality of rolling stands, there are problems as described below.
For example, when the technique of Patent Document 1 is applied at the time of acceleration / deceleration of the rolling material, parameter change based on "expert knowledge by qualitative expression" possessed by a person familiar with the plant or parameter change by inference means, As described in Patent Document 1, it is ambiguous information and it is clear that accurate plant control cannot be performed. Also, with regard to a method of changing parameters by an identification mechanism, the present technology is limited to a linear system and cannot identify a nonlinear system.

また、特許文献2の技術で、圧延荷重式を用いて圧延荷重を計算するためには、入側板厚、出側板厚、前方張力、後方張力のデータが必要であるが、センサがないためにデータが得られない場合が往々にしてある。この場合には、特許文献2の手法を適用することができない。また、データが得られた場合であっても、センサの取付け位置に起因するむだ時間や信号処理による遅れ時間の発生による影響があり、入側板厚、出側板厚、前方張力、後方張力のすべてのデータについて同時性、すなわち同一時刻のデータであることが保てない。つまり、特許文献2の圧延荷重式は静的であるため、同時性を保てないデータの誤差は考慮されておらず、圧延の加減速時のような過渡状態を精確には表現できない。   Moreover, in order to calculate a rolling load using the rolling load formula with the technique of patent document 2, the data of an entrance side plate thickness, an exit side plate thickness, front tension, and back tension are required, but since there is no sensor. Often there is no data available. In this case, the method of Patent Document 2 cannot be applied. Even when data is obtained, there are effects due to dead time caused by the sensor mounting position and delay time due to signal processing, and all of the inlet side thickness, outlet side thickness, front tension, and rear tension Cannot be maintained at the same time, that is, at the same time. That is, since the rolling load formula of Patent Document 2 is static, an error in data that cannot maintain simultaneity is not considered, and a transient state such as during rolling acceleration / deceleration cannot be accurately expressed.

また、特許文献2の技術は単圧延スタンドに対するものであり、複数の圧延スタンドを備えた圧延機に適用するには、他の圧延スタンドとの干渉が考慮されておらず技術的困難を伴う。
つまり、複数の圧延スタンドを備えた圧延機では、各圧延スタンドが圧延材や制御系を介して互いに干渉していて、例えば、上流側の圧延スタンドの出側板厚偏差は下流側の圧延スタンドの出側板厚偏差に影響する。このため、上流側の圧延スタンドの動特性を考慮する必要がある。また、下流側の板厚や張力の偏差が上流側に及び、この上流側の影響が再び下流側の圧延スタンドの出側板厚偏差となる場合もある。これは、多スタンドの圧延機全体が、圧延材によって一体に連結された系であり、速度サクセシブや板厚制御等の制御系によって一体に連結された系であるからである。このため、複数の圧延スタンドを備えた圧延機の場合には全体のシミュレーションが必要である。特に、仕上板厚の高精度シミュレーションのためには、上流側の圧延スタンドのシミュレーションの精度向上が重要である。
Moreover, the technique of patent document 2 is for a single rolling stand, and it is technically difficult to apply to a rolling mill having a plurality of rolling stands because interference with other rolling stands is not considered.
In other words, in a rolling mill equipped with a plurality of rolling stands, each rolling stand interferes with each other via a rolling material and a control system.For example, the exit side plate thickness deviation of the upstream rolling stand is that of the downstream rolling stand. Affects the thickness deviation on the delivery side. For this reason, it is necessary to consider the dynamic characteristics of the upstream rolling stand. Further, the deviation of the downstream thickness and tension may be on the upstream side, and this upstream influence may be the downstream thickness deviation of the downstream rolling stand again. This is because the entire multi-stand rolling mill is a system that is integrally connected by a rolled material, and is a system that is integrally connected by a control system such as speed-successive or plate thickness control. For this reason, in the case of a rolling mill provided with a plurality of rolling stands, an entire simulation is required. In particular, for the high-precision simulation of the finishing plate thickness, it is important to improve the accuracy of the simulation of the upstream rolling stand.

また、特許文献2の技術では、圧延に大きな影響を及ぼす摩擦係数に着目し、その誤差が小さくなるように摩擦係数モデルを用いて摩擦係数を同定している。しかしながら、圧延材の変形抵抗モデルの誤差や、制御系における張力の不感帯等の非線形性や、制御系の過渡応答は考慮されていない。そのため、圧延全体のシミュレーションとして、誤差が大になる可能性が大きい。
以上述べた課題ゆえに、複数の圧延スタンドを備えた圧延機で圧延材を加速及び/又は減速する際の制御に、特許文献1,2の技術を採用したとしても、最適に調整された圧延制御パラメータを算出することが困難な状況となっていた。
In the technique of Patent Document 2, attention is paid to a friction coefficient that greatly affects rolling, and the friction coefficient is identified using a friction coefficient model so that the error is reduced. However, errors in the deformation resistance model of the rolled material, non-linearity such as a dead zone of tension in the control system, and transient response of the control system are not considered. Therefore, there is a high possibility that the error becomes large as a simulation of the whole rolling.
Because of the problems described above, even when the techniques of Patent Documents 1 and 2 are adopted for the control when accelerating and / or decelerating the rolling material with a rolling mill equipped with a plurality of rolling stands, the optimally adjusted rolling control is adopted. It was difficult to calculate the parameters.

そこで、本発明は、上記問題点を鑑み、圧延材とこの圧延材を圧延する複数の圧延スタンドとこの圧延スタンドの板厚制御を行う制御系とを動的に表現する圧延動的モデルを用いて、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗を基に圧延状態を再現するシミュレーションを行って、最適に調整された圧延制御パラメータを算出する圧延制御パラメータの算出方法及び算出装置を提供することを目的とする。
加えて、前記圧延動的モデルを用いて、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗により高精度に圧延状態を再現する圧延シミュレーション装置を提供することを目的とする。
In view of the above problems, the present invention uses a rolling dynamic model that dynamically expresses a rolled material, a plurality of rolling stands that roll the rolled material, and a control system that controls the thickness of the rolling stand. The friction coefficient and / or deformation resistance of the entire rolling mill is identified, and the simulation is performed to reproduce the rolling state based on the obtained friction coefficient and / or deformation resistance, and the optimally adjusted rolling control parameters are calculated. An object of the present invention is to provide a rolling control parameter calculation method and calculation device.
In addition, a rolling simulation apparatus that identifies the friction coefficient and / or deformation resistance of the entire rolling mill using the rolling dynamic model and reproduces the rolling state with high accuracy by the obtained friction coefficient and / or deformation resistance. The purpose is to provide.

前記目的を達成するため、本発明においては以下の技術的手段を講じた。
すなわち、本発明に係る圧延制御パラメータの算出方法は、圧延材と該圧延材を圧延する複数の圧延スタンドと該圧延スタンドの板厚制御を行う制御系とを動的に表現する圧延動的モデルを構築し、前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデルでの計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化し、最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデルを用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出することを特徴とする。
In order to achieve the above object, the present invention takes the following technical means.
That is, the rolling control parameter calculation method according to the present invention is a rolling dynamic model that dynamically expresses a rolled material, a plurality of rolling stands that roll the rolled material, and a control system that controls the thickness of the rolling stand. And the calculated value and the actual value of the rolling dynamic model with respect to the time change of the rolling load and / or the time change of the roll gap at the time of acceleration and / or deceleration of the rolled material. The parameters related to the friction coefficient and / or the deformation resistance at each rolling stand are optimized so that the evaluation function including the difference is minimized, and the optimized parameters related to the friction coefficient and / or the deformation resistance are incorporated. Further, using the rolling dynamic model, the rolling control parameters including the setting parameters for the plate thickness control and / or the setting parameters for the rolling stand are used. And wherein the optimizing to calculate a meter.

なお、ここで言う圧延制御パラメータとは、板厚制御に関する設定パラメータ及び/又は圧延機に関する設定パラメータであり、加減速補償量、制御ゲイン、不感帯幅、張力目標値、チューニング率、ドループ率、これらの速度依存特性等、圧延状態に影響を与え、且つ設定可能なパラメータのことをいう。設定可能なパラメータとしては、数値だけではなく、テーブルや関数の形でもよい。また、評価関数が最小とは、符号を変更して最大化することを含む。以下同様である。
この圧延制御パラメータの算出方法によれば、圧延の過渡特性を表現可能な圧延動的モデルを用いて、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗を基に圧延状態を精確に再現するシミュレーションを行って、最適に調整された圧延制御パラメータを算出することができる。
Note that the rolling control parameter referred to here is a setting parameter related to sheet thickness control and / or a setting parameter related to a rolling mill, such as acceleration / deceleration compensation amount, control gain, dead band width, target tension value, tuning rate, droop rate, etc. This is a parameter that affects the rolling state and can be set, such as speed-dependent characteristics. The parameters that can be set are not only numerical values but also tables and functions. The minimum evaluation function includes changing the sign to maximize it. The same applies hereinafter.
According to this rolling control parameter calculation method, a rolling dynamic model capable of expressing rolling transient characteristics is used to identify the friction coefficient and / or deformation resistance of the entire rolling mill, and the obtained friction coefficient and / or A simulation that accurately reproduces the rolling state based on the deformation resistance can be performed to calculate the optimally adjusted rolling control parameter.

かかる圧延動的モデルを用いた圧延システム全体のシミュレーションは、実際には板厚や張力が測定不能な場合でも、それらの値を算出可能であり、摩擦係数及び/又は変形抵抗に関するパラメータや圧延制御パラメータ等の各種パラメータを調整できる。なお、ここで言う測定不能な場合とは、センサ自体が存在しない場合や、存在しても内部で計測値を利用するだけで外部に取り出せず利用できない場合や、精度が不十分で利用できない場合や、測定時刻に誤差が大きく利用できない場合等、精確なデータを利用できない場合をいうものとする。板厚や張力に測定の遅れ、例えば、板の移送距離による遅れやフィルタによる遅れがある場合、すなわち、実績データのすべての項目に同時性がない場合であっても、各種パラメータを調整できる。   The simulation of the entire rolling system using such a rolling dynamic model can calculate those values even when the plate thickness and tension cannot be measured in practice, and parameters related to the friction coefficient and / or deformation resistance and rolling control. Various parameters such as parameters can be adjusted. Note that the case where measurement is impossible here means that the sensor itself does not exist, or even if it exists, it cannot be used because it cannot be extracted outside simply by using the measured value inside, or it cannot be used due to insufficient accuracy. In addition, it means a case where accurate data cannot be used, such as when the measurement time has a large error and cannot be used. Various parameters can be adjusted even when there is a measurement delay in the plate thickness or tension, for example, a delay due to a plate transfer distance or a delay due to a filter, that is, there is no simultaneity of all items of the performance data.

また、本技術的手段を用いることで、実際の圧延機や圧延機を制御する制御部に対し特殊な操作を行うことなく、各種パラメータを同定することができ操業上安全である。
好ましくは、前記摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータの最適化の計算を、上流側の圧延スタンドから下流側の圧延スタンドへ順に行うとよい。
また、前記摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータの最適化の計算を、圧延スタンド毎に1回又は複数回実施するとよい。
前記摩擦係数に関するパラメータは、前記摩擦係数が変化しなくなる上限の圧延速度である上限圧延速度と、圧延速度0での摩擦係数を規定する変数とを備えていることが好ましい。
Further, by using this technical means, various parameters can be identified without performing any special operation on the actual rolling mill or a control unit for controlling the rolling mill, which is safe in operation.
Preferably, calculation of optimization of the parameter relating to the friction coefficient and / or the parameter relating to deformation resistance may be performed in order from the upstream rolling stand to the downstream rolling stand.
In addition, the calculation of the parameter relating to the friction coefficient and / or the parameter relating to the deformation resistance may be performed once or a plurality of times for each rolling stand.
The parameter relating to the friction coefficient preferably includes an upper limit rolling speed that is an upper limit rolling speed at which the friction coefficient does not change, and a variable that defines a friction coefficient at a rolling speed of zero.

また、本発明に係る圧延制御パラメータの算出装置は、圧延材と該圧延材を圧延する複数の圧延スタンドとをモデル化した圧延機モデル部、及び前記圧延スタンドの板厚制御を行う板厚制御モデル部を備えた圧延動的モデル部と、前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデル部での計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化する圧延状態パラメータ最適化部と、前記圧延状態パラメータ最適化部により最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデル部を用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出する圧延制御パラメータ最適化部と、を有することを特徴とする。   The rolling control parameter calculation apparatus according to the present invention includes a rolling mill model unit that models a rolled material and a plurality of rolling stands that roll the rolled material, and a plate thickness control that performs plate thickness control of the rolling stand. The rolling dynamic model part provided with a model part, and the rolling dynamic model part regarding the time change of the rolling load and / or the time change of the roll gap at the time of acceleration and / or deceleration of the rolled material A rolling state parameter optimizing unit that optimizes a parameter relating to a friction coefficient and / or a parameter relating to deformation resistance in each rolling stand so that an evaluation function including a difference between a calculated value and an actual value is minimized. The rolling dynamic mode in which the parameter relating to the friction coefficient and / or the parameter relating to deformation resistance optimized by the rolling state parameter optimization unit is incorporated. It is using pole tip, and having a a rolling control parameter optimization unit for calculating and optimizing the rolling control parameters including configuration parameters configuration parameters and / or rolling stands of the thickness control.

なお、ここで言う圧延状態パラメータとは、圧延機(圧延材と圧延スタンド)に関するパラメータであり、摩擦係数、変形抵抗、これらの速度依存特性等、圧延状態に影響を与え、且つ設定できないパラメータをいう。圧延状態パラメータは、数値だけではなく、テーブルや関数の形でもよい。
この圧延制御パラメータの算出装置を用いることで、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗を基に圧延状態を精確に再現するシミュレーションを行って、最適に調整された圧延制御パラメータを算出することができる。
In addition, the rolling state parameter said here is a parameter regarding a rolling mill (rolled material and rolling stand), and influences a rolling state, such as a friction coefficient, deformation resistance, and these speed-dependent characteristics, and parameters that cannot be set. Say. The rolling state parameter may be not only a numerical value but also a table or a function.
Using this rolling control parameter calculation device, the friction coefficient and / or deformation resistance of the entire rolling mill is identified, and a simulation is performed to accurately reproduce the rolling state based on the obtained friction coefficient and / or deformation resistance. Thus, the optimally adjusted rolling control parameter can be calculated.

また、本発明に係る圧延シミュレーション装置は、圧延材と該圧延材を圧延する複数の圧延スタンドとをモデル化した圧延機モデル部、及び前記圧延スタンドの板厚制御を行う板厚制御モデル部を備えた圧延動的モデル部と、前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデル部での計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化する圧延状態パラメータ最適化部と、前記圧延状態パラメータ最適化部により最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデル部を用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出する圧延制御パラメータ最適化部と、を有していて、前記圧延動的モデル部は、圧延状態パラメータ最適化部で最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータと、圧延制御パラメータ最適化部で最適化された圧延制御パラメータとを用いて、圧延の状況を再現できるように構成されていることを特徴とする。   The rolling simulation apparatus according to the present invention includes a rolling mill model unit that models a rolled material and a plurality of rolling stands that roll the rolled material, and a plate thickness control model unit that performs plate thickness control of the rolling stand. The rolling dynamic model section provided, and the rolling dynamic model section calculates the time change of the rolling load and / or the time change of the roll gap at the time of acceleration and / or deceleration of the rolled material. A rolling state parameter optimizing unit that optimizes a parameter relating to a friction coefficient and / or a parameter relating to deformation resistance in each rolling stand so that an evaluation function including a difference between a value and an actual value is minimized, and the rolling state parameter The rolling dynamic model incorporating a parameter relating to a friction coefficient and / or a parameter relating to deformation resistance optimized by an optimization unit A rolling control parameter optimization unit that optimizes and calculates a rolling control parameter including a setting parameter of the plate thickness control and / or a setting parameter of a rolling stand, and the rolling dynamic model unit is The rolling condition can be reproduced by using the friction coefficient parameter and / or deformation resistance parameter optimized by the rolling condition parameter optimization unit and the rolling control parameter optimized by the rolling control parameter optimization unit. It is comprised as follows.

この圧延シミュレーション装置を用いることで、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗を基に圧延状態を精確に再現することができる。   By using this rolling simulation apparatus, the friction coefficient and / or deformation resistance of the entire rolling mill can be identified, and the rolling state can be accurately reproduced based on the obtained friction coefficient and / or deformation resistance.

本発明にかかる圧延制御パラメータの算出方法及び算出装置を用いることで、圧延材の加減速時に圧延材の形状を悪化させることなく板厚精度を向上させることが可能な圧延制御パラメータを算出できる。
本発明にかかる圧延シミュレーション装置を用いることで、圧延機全体の摩擦係数及び/又は変形抵抗を同定し、得られた摩擦係数及び/又は変形抵抗により高精度に圧延状態を再現することができる。
By using the rolling control parameter calculation method and calculation apparatus according to the present invention, it is possible to calculate a rolling control parameter capable of improving the plate thickness accuracy without deteriorating the shape of the rolled material during acceleration / deceleration of the rolled material.
By using the rolling simulation apparatus according to the present invention, the friction coefficient and / or deformation resistance of the entire rolling mill can be identified, and the rolling state can be reproduced with high accuracy by the obtained friction coefficient and / or deformation resistance.

以下、本発明にかかる圧延制御パラメータの算出方法及び算出装置、圧延シミュレーション装置の実施の形態を、図を基に説明する。
本発明にかかる圧延制御パラメータの算出装置1(圧延シミュレーション装置)はコンピュータ上で実現されるシミュレーションモデルである。使用ツールまたは使用言語としては、例えば、The Mathworks社のMatlab/Simulink(Matlab、Simulinkは、The Mathworks社の登録商標)、C言語、Fortran言語が挙げられるが、これらに限定されるものではない。
Embodiments of a rolling control parameter calculation method, calculation device, and rolling simulation device according to the present invention will be described below with reference to the drawings.
The rolling control parameter calculation device 1 (rolling simulation device) according to the present invention is a simulation model realized on a computer. Examples of the tool or language used include, but are not limited to, The Mathworks Matlab / Simulink (Matlab, Simulink is a registered trademark of The Mathworks), C language, and Fortran language.

図1には、圧延制御パラメータの算出装置1のブロック図が示されている。
この図に示す如く、本算出装置1は、圧延スタンド2及びその圧延スタンド2で冷間圧延される圧延材3をモデル化した圧延機モデル部4と、圧延スタンド2のロール速度やロールギャップを操作する板厚制御モデル部5とを有している。圧延機モデル部4と板厚制御モデル部5とで圧延動的モデル部6が構成されている。
さらに、算出装置1は、圧延制御パラメータを最適化する圧延制御パラメータ最適化部7と、最適化された摩擦係数及び/又は変形抵抗を含む圧延状態パラメータを計算する圧延状態パラメータ最適化部8と、を備えている。
FIG. 1 shows a block diagram of a rolling control parameter calculation device 1.
As shown in this figure, the calculation device 1 includes a rolling stand 2 and a rolling mill model unit 4 that models a rolling material 3 that is cold-rolled by the rolling stand 2, and the roll speed and roll gap of the rolling stand 2. And a plate thickness control model unit 5 to be operated. The rolling mill model unit 4 and the plate thickness control model unit 5 constitute a rolling dynamic model unit 6.
Further, the calculation device 1 includes a rolling control parameter optimization unit 7 that optimizes rolling control parameters, and a rolling state parameter optimization unit 8 that calculates rolling state parameters including an optimized friction coefficient and / or deformation resistance. It is equipped with.

なお、ここで言う圧延制御パラメータとは、板厚制御に関する設定パラメータや圧延機に関する設定パラメータであり、加減速補償量、制御ゲイン、不感帯幅、張力目標値、チューニング率、ドループ率、これらの速度依存特性等、圧延状態に影響を与え、且つ設定可能なパラメータのことをいう。
圧延状態パラメータとは、圧延機(圧延材と圧延スタンド)に関するパラメータであり、摩擦係数、変形抵抗、これらの速度依存特性等、圧延状態に影響を与え、且つ設定できないパラメータをいう。
The rolling control parameters referred to here are setting parameters related to sheet thickness control and setting parameters related to the rolling mill, and include acceleration / deceleration compensation amount, control gain, dead band width, target tension value, tuning rate, droop rate, and these speeds. It refers to parameters that influence the rolling state and can be set, such as dependency characteristics.
The rolling state parameter is a parameter relating to a rolling mill (rolled material and rolling stand), and refers to a parameter that affects the rolling state and cannot be set, such as a friction coefficient, deformation resistance, and speed-dependent characteristics thereof.

図2には、圧延機モデル部4によりモデル化されている「圧延スタンド2及び圧延材3」の模式図が示されている。
圧延スタンド2は複数(#1〜#5)配備されていてタンデム型の圧延機9となっている。#1圧延スタンド2に元板が通された後、#2〜#4圧延スタンド2を通過する毎に圧下され、#5圧延スタンド2を出たところで所定の仕上げ板厚となるように、各圧延スタンド2のロール速度やロールギャップが制御されるものとなっている。
圧延機モデル部4は圧延スタンド2のモデルを有している。圧延スタンド2のモデルとしては、「板圧延の理論と実際,日本鉄鋼協会編,1984」のp111〜124に開示されたものが採用可能であり、式(1)〜式(9)に示すようなものとなっている。
FIG. 2 shows a schematic diagram of the “rolling stand 2 and rolled material 3” modeled by the rolling mill model unit 4.
A plurality of rolling stands 2 (# 1 to # 5) are provided, and a tandem rolling mill 9 is provided. After the base plate is passed through the # 1 rolling stand 2, each time it passes through the # 2 to # 4 rolling stand 2, it is squeezed down, and when it exits the # 5 rolling stand 2, The roll speed and roll gap of the rolling stand 2 are controlled.
The rolling mill model unit 4 has a model of the rolling stand 2. As the model of the rolling stand 2, those disclosed in p111 to 124 of "Theory and practice of plate rolling, edited by the Japan Iron and Steel Institute, 1984" can be adopted, as shown in equations (1) to (9). It has become a thing.

Figure 2009208115
Figure 2009208115

圧延材3の圧延スタンド2,2間の移送については、上流側に位置する圧延スタンド2の出側板厚が当該圧延スタンド2へ出側板速度で移動するとしている。
従って、圧延スタンド2,2間の圧延材3を複数の部分(サンプル板片)に分割すると共に、最適化計算のサンプリング周期をΔtsとすると、サンプリング周期ごとにサンプル板片がVout,i・Δtsだけ移動することになって、#i+1圧延スタンド2の位置に達したサンプル板片の板厚を#i+1スタンドの入側板厚Hi+1とすることができる。
まず、摩擦係数μiを最適化する場合には、圧延モデル式中の摩擦係数μiの初期値としては、実際の圧延機9から得られた実績値から推定または学習した値を設定する。例えば、特開平10−249421号公報「圧延機の制御装置」の図5(a)の如く、圧延速度が低速の場合は、摩擦係数が大きく、圧延速度がある所定速度以上の高速の場合は、摩擦係数が小さく且つ一定値となるように設定するとよい。なお、摩擦係数μiを最適化しない場合には、摩擦係数μiは初期値のままであり後述の最適化は行わない。
Regarding the transfer of the rolling material 3 between the rolling stands 2 and 2, the exit side plate thickness of the rolling stand 2 located on the upstream side is moved to the rolling stand 2 at the exit side plate speed.
Therefore, it becomes possible to divide the rolled material 3 between the mill stands 2, 2 into a plurality of portions (sample plate piece), the sampling period of optimization calculation When Delta] t s, the sample plate piece V out for each sampling period, i -By moving by Δt s, the plate thickness of the sample plate piece that has reached the position of the # i + 1 rolling stand 2 can be made the entry side plate thickness H i + 1 of the # i + 1 stand.
First, when optimizing the friction coefficient mu i as the initial value of the friction coefficient mu i during rolling model equation, set the estimated or learned value from the actual result value obtained from the rolling mill 9. For example, as shown in FIG. 5A of Japanese Patent Laid-Open No. 10-249421 “Control device for rolling mill”, when the rolling speed is low, the friction coefficient is large, and when the rolling speed is higher than a predetermined speed, The friction coefficient is preferably set to be small and constant. If the friction coefficient μ i is not optimized, the friction coefficient μ i remains the initial value, and optimization described later is not performed.

変形抵抗kiを最適化する場合には、変形抵抗kiの初期値としては、静的変形抵抗や歪み速度に依存した式等の公知の式、すなわち「板圧延の理論と実際、日本鉄鋼協会編、1984」のp188〜189に開示された式(7.51)〜式(7.53)のいずれかより求めることができる。なお、変形抵抗kiを最適化しない場合には、変形抵抗kiは初期値のままであり後述の最適化は行わない。
また、前回のサンプリング時刻のスタンド間張力が求められているとし、圧延スタンド2の入側・出側張力が実績値や実績パラメータから得られているため、式(3)から先進率を求めることができ、式(1),式(2)から、各圧延スタンド2の入側・出側板速を求めることができる。ここで、式(1),式(2)のロール速度は式(4)を用いて求めることができる。
When optimizing the deformation resistance k i as the initial value of the deformation resistance k i, known of the formula or the like depending on the static deformation resistance and the strain rate, i.e. Theory and Practice "plate rolling, the Japan Iron and Steel It can be obtained from any one of formulas (7.51) to (7.53) disclosed in p. When the deformation resistance k i is not optimized, the deformation resistance k i remains at the initial value, and optimization described later is not performed.
Also, assuming that the tension between the stands at the previous sampling time is obtained, and the entry and exit tensions of the rolling stand 2 are obtained from the actual values and actual parameters, the advanced rate is obtained from the equation (3). From the equations (1) and (2), the entry / exit plate speeds of the respective rolling stands 2 can be obtained. Here, the roll speed of Formula (1) and Formula (2) can be calculated | required using Formula (4).

さらに、式(7)を用いて電流を計算でき、式(8),式(9)を用いて現サンプリング時刻のスタンド間張力を求めることができる。スタンド出側板厚と圧延荷重は、式(5),式(6)を連立し、求めることができる。
なお、先進率(式(3))の具体的な式としては、「板圧延の理論と実際、日本鉄鋼協会編、1984」のp33〜35に開示されたBland&Fordの方法が挙げられる。
後進率については、(1+bi)・Hi=(1+fi)・hiのマスフロー一定則から求められる。
Furthermore, the current can be calculated using the equation (7), and the tension between the stands at the current sampling time can be obtained using the equations (8) and (9). The stand exit side plate thickness and the rolling load can be obtained by simultaneous equations (5) and (6).
In addition, as a specific formula of the advanced rate (formula (3)), the method of Brand & Ford disclosed in p33 to 35 of “Theory and practice of plate rolling, edited by Japan Iron and Steel Institute, 1984” can be mentioned.
The reverse rate is obtained from a constant mass flow law of (1 + b i ) · H i = (1 + f i ) · h i .

圧延荷重(式(6))の具体的な式は、「板圧延の理論と実際、日本鉄鋼協会編、1984」のp19〜36に開示された式(2.97),式(2.185),式(2.194)のHillの式から求めることができる。
電流(式(7))の具体的な式は、「板圧延の理論と実際、日本鉄鋼協会編、1984」のp20〜35の式(2.104),式(2.186),式(2.192)のHillの式のトルクを電流に変換して用いればよい。例えば、トルク定数で除することにより電流を求める。
The specific formula of the rolling load (formula (6)) is the formula (2.97) and formula (2.185) disclosed in p19 to 36 of "Theory and practice of plate rolling, edited by Japan Iron and Steel Institute, 1984". ), And Hill's equation (2.194).
The specific formula of the current (formula (7)) is the formula (2.104), formula (2.186), formula (p.20 to p35 of “Theory and practice of plate rolling, edited by the Japan Steel Association, 1984”) The torque of the Hill equation of 2.192) may be converted into a current and used. For example, the current is obtained by dividing by the torque constant.

以上より、#1圧延スタンド2〜#5圧延スタンド2のシミュレーションが可能となる。これらの式に対して、圧延速度パターン、元板の板厚履歴や変形抵抗履歴を入力することで、仕上板厚、圧延荷重、ロールギャップの変動を再現可能となる。
実際の圧延スタンド2では、通常、各スタンドの圧延荷重を測定する荷重計は各圧延スタンド2に設置されている。しかし、板厚計や張力計は設置されていない場合があり、このような場合には、圧延スタンド2出側板厚や張力の値を知ることはできないが、上記したシミュレーションモデルを用いれば、測定できないスタンド出側板厚等を算出することが可能である。
From the above, simulation of # 1 rolling stand 2 to # 5 rolling stand 2 is possible. By inputting the rolling speed pattern, the sheet thickness history of the base plate, and the deformation resistance history with respect to these equations, it is possible to reproduce variations in the finished plate thickness, rolling load, and roll gap.
In an actual rolling stand 2, a load meter that measures the rolling load of each stand is usually installed in each rolling stand 2. However, there are cases where a thickness gauge and tension meter are not installed. In such a case, the thickness and tension value of the rolling stand 2 cannot be known. It is possible to calculate the stand outlet side plate thickness that cannot be calculated.

なお、圧延スタンド2のモデルには、最上流側圧延スタンド2の入側の動特性、及び最下流側圧延スタンド2の出側の動特性、例えば、巻き取り装置(図示せず)による圧延スタンド2出側の張力変化を現す関係式などを含めてもよい。また、圧延スタンド2のモータ制御や圧延スタンド2の剛性の制御など、ローカルな制御が行われている場合には、圧延スタンド2、圧延材3部にこれらの制御を含めるものとする。また、圧延スタンド2,2間の張力の動特性を含むものとする。
板厚制御モデル部5は、目標板厚設定部10、目標張力設定部11、板厚張力制御部12、加減速補償部13、圧延速度パターン決定部14から構成されている。
The rolling stand 2 model includes dynamic characteristics on the inlet side of the most upstream rolling stand 2 and dynamic characteristics on the outlet side of the most downstream rolling stand 2, for example, a rolling stand by a winding device (not shown). (2) A relational expression representing a change in tension on the outlet side may be included. When local control such as motor control of the rolling stand 2 or rigidity control of the rolling stand 2 is performed, these controls are included in the rolling stand 2 and the rolling material 3 part. The dynamic characteristics of the tension between the rolling stands 2 and 2 are included.
The plate thickness control model unit 5 includes a target plate thickness setting unit 10, a target tension setting unit 11, a plate thickness tension control unit 12, an acceleration / deceleration compensation unit 13, and a rolling speed pattern determination unit 14.

なお、張力の制御を行わない場合には、目標張力設定部11は不要である。また、板厚張力制御部12は、少なくとも#5圧延スタンド2の出側板厚(仕上板厚)を一定に制御するものであるが、それ以外の圧延スタンド2での出側板厚の全て又は一部を制御するもの、あるいは、スタンド間の張力の全て又は一部を制御するものを含む。
図1,図3に示すように、板厚張力制御部12には、目標板厚設定部10において予め設定された#1〜#5圧延スタンド2の出側板厚の目標値が入力される。同様に、板厚張力制御部12には、目標張力設定部11において予め設定された圧延スタンド2間張力の目標値が入力される。なお、図3において、#12スタンド間張力とは、#1圧延スタンド2と#2圧延スタンド2との間における圧延材3の張力を意味する。また、図3は、全スタンドの出側板厚と全スタンド間張力を制御する場合における板厚張力制御部12の一例を示したものである。
Note that the target tension setting unit 11 is not required when tension control is not performed. Further, the plate thickness tension control unit 12 controls at least the outlet side plate thickness (finished plate thickness) of the # 5 rolling stand 2 to be constant, but all or one of the outgoing side plate thicknesses in the other rolling stands 2 are also controlled. Those that control the part, or those that control all or part of the tension between the stands.
As shown in FIGS. 1 and 3, the target value of the exit side plate thickness of the # 1 to # 5 rolling stands 2 preset in the target plate thickness setting unit 10 is input to the plate thickness tension control unit 12. Similarly, a target value of the tension between the rolling stands 2 preset in the target tension setting unit 11 is input to the plate thickness tension control unit 12. In FIG. 3, the tension between # 12 stands means the tension of the rolling material 3 between the # 1 rolling stand 2 and the # 2 rolling stand 2. FIG. 3 shows an example of the plate thickness tension control unit 12 in the case of controlling the exit side plate thickness of all stands and the tension between all stands.

板厚張力制御部12の出力は、#1〜#5圧延スタンド2におけるロールギャップと、#1〜#4圧延スタンド2におけるロール速度である。
板厚張力制御部12において行われる制御としては、公知のものが採用可能であり、PI制御や最適制御を用いるとよい。例えば「冷延タンデムミル高精度板厚制御技術の開発」(電気学会研究会資料,金属産業研究会,MID−01〜9,2001年3月16日)の手法を用いてもよい。
さらに、各圧延スタンド2においては、出側板厚や圧延スタンド2,2間の張力が測定又は計算され、その結果は板厚張力制御部12にフィードバック入力される。
The outputs of the plate thickness tension control unit 12 are the roll gap in the # 1 to # 5 rolling stands 2 and the roll speed in the # 1 to # 4 rolling stands 2.
As the control performed in the plate thickness tension control unit 12, a publicly known one can be adopted, and PI control or optimum control may be used. For example, the method of “development of cold rolling tandem mill high-accuracy sheet thickness control technology” (The Institute of Electrical Engineers of Japan, Metal Industry Research Society, MID-01 to 9, March 16, 2001) may be used.
Further, in each rolling stand 2, the exit side plate thickness and the tension between the rolling stands 2 and 2 are measured or calculated, and the result is fed back to the plate thickness tension control unit 12.

また、本実施形態の場合、図1に示すように、基準となる圧延スタンド2(#5圧延スタンド)での圧延速度パターンを決定する圧延速度パターン決定部14が備えられている。
この圧延速度パターン決定部14では、1コイル(1本の圧延材3)を圧延するに際し、#5圧延スタンド2における圧延速度パターン(圧延速度の時間による変化)が決定される。例えば、図4に示すように、加速度を指定して台形パターンとしてもよい。すなわち、圧延材3の先端部においては、圧延速度を徐々に加速し、圧延材3の終端部においては、圧延速度を徐々に減速し、溶接点(当該圧延材と次圧延材との接続点)において最低速度とするようにしている。圧延材3の中途部では圧延速度を最高且つほぼ一定としている。なお、圧延速度を加減速する場合、板厚偏差が大きくなることは従来より明らかとなっている。これは、加減速時に圧延ロールと圧延材3との間の摩擦係数等が変化することが大きく関与している。
In the case of the present embodiment, as shown in FIG. 1, a rolling speed pattern determining unit 14 that determines a rolling speed pattern in the reference rolling stand 2 (# 5 rolling stand) is provided.
The rolling speed pattern determining unit 14 determines a rolling speed pattern (change in rolling speed with time) in the # 5 rolling stand 2 when rolling one coil (one rolled material 3). For example, as shown in FIG. 4, a trapezoid pattern may be specified by designating acceleration. That is, the rolling speed is gradually increased at the front end portion of the rolled material 3, and the rolling speed is gradually decreased at the end portion of the rolled material 3, and the welding point (the connection point between the rolled material and the next rolled material). ) Is set to the minimum speed. In the middle of the rolled material 3, the rolling speed is the highest and substantially constant. In addition, when accelerating / decelerating a rolling speed, it has become clear from the past that a plate | board thickness deviation becomes large. This is greatly related to a change in the coefficient of friction between the rolling roll and the rolled material 3 during acceleration / deceleration.

さらに、図1に示す如く、圧延材3が加速又は減速するに際し、板厚張力制御部12の出力に対する出力補正量を算出する加減速補償部13と、圧延制御パラメータを最適化する圧延制御パラメータ最適化部7と、が設けられている。
詳しくは、加減速補償部13の出力補正量としては、各圧延スタンド2のロールギャップとロール速度とを採用しており、加減速補償部13は、図5,図6に示すパターンに従って、ロールギャップの補正量、ロール速度の補正量を導き出す。
圧延制御パラメータ最適化部7は、圧延材3の加速時及び/又は減速時における圧延材3の板厚偏差を含む評価関数が小さくなるように、加減速補償量等の所定の圧延制御パラメータを最適化する機能を備えている。すなわち、圧延の過渡特性を表現可能な圧延動的モデルを用いて最適化し、最適化された圧延制御パラメータを出力するようになっている。最適化時には、評価関数を最小化する非線形最適化問題を解くことになる。
Further, as shown in FIG. 1, when the rolled material 3 is accelerated or decelerated, an acceleration / deceleration compensator 13 for calculating an output correction amount with respect to the output of the sheet thickness tension controller 12, and a rolling control parameter for optimizing the rolling control parameter. And an optimization unit 7.
Specifically, as the output correction amount of the acceleration / deceleration compensation unit 13, the roll gap and the roll speed of each rolling stand 2 are adopted, and the acceleration / deceleration compensation unit 13 rolls according to the patterns shown in FIGS. Gap correction amount and roll speed correction amount are derived.
The rolling control parameter optimization unit 7 sets predetermined rolling control parameters such as an acceleration / deceleration compensation amount so that an evaluation function including a thickness deviation of the rolled material 3 at the time of acceleration and / or deceleration of the rolled material 3 becomes small. Has the ability to optimize. That is, optimization is performed using a rolling dynamic model capable of expressing the transient characteristics of rolling, and optimized rolling control parameters are output. At the time of optimization, a nonlinear optimization problem that minimizes the evaluation function is solved.

一方、本発明にかかる圧延制御パラメータの算出装置1は、圧延状態パラメータ最適化部8を有している。この圧延状態パラメータ最適化部8において、圧延材3を加速及び/又は減速する際のシミュレーションを精確なものとするために、各圧延スタンド2での圧延荷重及び/又はロールギャップを評価した上で、摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適に調整し算出する。
具体的には、圧延材3を加速及び/又は減速するに際し、複数の圧延スタンド2での圧延荷重の時間変化と実績圧延荷重の時間変化との差、及び/又は、複数の圧延スタンド2でのロールギャップの時間変化と実績ロールギャップの時間変化との差を含む評価関数が最小となるようにする。評価関数として圧延荷重、ロールギャップを選定した理由としては、板厚は非線形性(例えば、張力のデッドバンド等)やシミュレーションの初期値によって影響を受けやすく、シミュレーションによる実圧延の再現が難しい場合があるものの、圧延荷重やロールギャップはその変化量が大きくその再現が容易で、摩擦係数や変形抵抗のパラメータを調整しやすいからである。
On the other hand, the rolling control parameter calculation device 1 according to the present invention includes a rolling state parameter optimization unit 8. In this rolling state parameter optimizing unit 8, in order to make accurate simulation when accelerating and / or decelerating the rolled material 3, the rolling load and / or roll gap in each rolling stand 2 is evaluated. The parameter relating to the friction coefficient and / or the parameter relating to the deformation resistance is optimally adjusted and calculated.
Specifically, when accelerating and / or decelerating the rolled material 3, the difference between the time change of the rolling load at the plurality of rolling stands 2 and the time change of the actual rolling load, and / or the plurality of rolling stands 2 The evaluation function including the difference between the time change of the actual roll gap and the time change of the actual roll gap is minimized. The reason for selecting the rolling load and roll gap as the evaluation function is that the plate thickness is easily affected by non-linearity (for example, dead band of tension) and the initial value of simulation, and it is difficult to reproduce actual rolling by simulation. However, the change in the rolling load and the roll gap is large and easy to reproduce, and it is easy to adjust the friction coefficient and deformation resistance parameters.

ここで一例として、圧延荷重の差に基づき、摩擦係数を調整する場合を説明する。調整の容易さを考慮し、減速時のデータにより調整するものとする。
圧延状態パラメータ最適化部8における摩擦係数モデルとしては、例えば、式(10),(11)を考える。
Here, as an example, a case where the friction coefficient is adjusted based on the difference in rolling load will be described. Considering the ease of adjustment, adjustment is made based on the data during deceleration.
As the friction coefficient model in the rolling state parameter optimization unit 8, for example, equations (10) and (11) are considered.

Figure 2009208115
Figure 2009208115

ここで、αiは、#i圧延スタンド2において圧延速度0での摩擦係数を規定する変数であり、Vout,i0は、摩擦係数の変化する圧延速度の上限値である。つまり、Vout,i0より圧延速度が速ければ摩擦係数はμi=μi0で変化しないとする。ここで、μi0は一定圧延速度でのデータに基づく同定により精確であるとして、以後の調整パラメータには含めない。
圧延スタンド2ごとに、αi,Vout,i0を調整する。なお、単スタンドの場合は1スタンドでの調整となる。また、μi0も未知パラメータとして、αi,Vout,i0と同時に調整することも可能である。
Here, α i is a variable that defines the friction coefficient at a rolling speed of 0 in the #i rolling stand 2, and V out, i0 is an upper limit value of the rolling speed at which the friction coefficient changes. That is, if the rolling speed is faster than V out, i0 , the friction coefficient is assumed not to change with μ i = μ i0 . Here, μ i0 is assumed to be accurate by identification based on data at a constant rolling speed, and is not included in the subsequent adjustment parameters.
For each rolling stand 2, α i , V out, i0 are adjusted. In the case of a single stand, the adjustment is made with one stand. Also, μ i0 can be adjusted as an unknown parameter simultaneously with α i , V out, i0 .

まず、圧延速度パターン決定部14に保持されている圧延速度パターンの時間履歴データを圧延動的モデル部6に入力する。そして、例えば、#1〜#4圧延スタンド2のαi,Vout,i0をαi=0.5,0.7,0.9,1.1の4通り、Vout,i0=600,800,1000(mpm)の3通りを取り得るとして、以後、最適化を行う。
また、#5圧延スタンド2の摩擦係数μ5は、式(11)で表されるとし、当該摩擦係数μ5を算出する。
First, the time history data of the rolling speed pattern held in the rolling speed pattern determining unit 14 is input to the rolling dynamic model unit 6. And, for example, α i , V out, i0 of # 1 to # 4 rolling stand 2 are four ways of α i = 0.5 , 0.7, 0.9, 1.1 , V out, i0 = 600, Assuming that three patterns of 800 and 1000 (mpm) can be taken, optimization is performed thereafter.
Further, the friction coefficient μ 5 of the # 5 rolling stand 2 is expressed by the equation (11), and the friction coefficient μ 5 is calculated.

Figure 2009208115
Figure 2009208115

#5圧延スタンド2のα5,Vout,50,cμ5に関しては、例えば、α5=0.5,0.7,0.9,1.1の4通り、Vout,50=600,800,1000(mpm)の3通り、cμ5=2.0,2.5,3.0,の3通りを取り得るとして、以後、最適化を行う。したがって、#1スタンド〜#4スタンドは、4×3=12通り、#5スタンドは、4×3×3=36通りの場合についてシミュレーションを行い、最適化することとなる。
式(10),式(11)のように設定することで、減速時、増速時の摩擦係数の変化を明示的に設定でき、仕上板厚変動がより精確に再現できるようになる。また、精確に再現された摩擦係数を用いることで、高速時と低速時のロールギャップ差を確実に設定でき、低速時の加減速補償量の調整が容易となる。
As for α 5 , V out, 50 , cμ 5 of # 5 rolling stand 2, for example, α 5 = 0.5 , 0.7, 0.9, 1.1, V out, 50 = 600, Assuming that three patterns of 800 and 1000 (mpm) and cμ 5 = 2.0, 2.5, and 3.0 can be taken, optimization is performed thereafter. Therefore, the simulation is performed and optimized for the case of # 1 stand to # 4 stand and 4 × 3 = 12 ways and the # 5 stand is 4 × 3 × 3 = 36 ways.
By setting as in Expression (10) and Expression (11), it is possible to explicitly set the change in the friction coefficient at the time of deceleration and at the time of acceleration, and the finished plate thickness fluctuation can be reproduced more accurately. In addition, by using the accurately reproduced friction coefficient, the difference in roll gap at high speed and low speed can be set reliably, and the acceleration / deceleration compensation amount at low speed can be easily adjusted.

次に、摩擦係数を調整する場合の評価関数について説明する。評価関数Jerrは、次式のようにしている。 Next, an evaluation function for adjusting the friction coefficient will be described. The evaluation function J err is expressed as follows.

Figure 2009208115
Figure 2009208115

ここで、Jerriは、#i圧延スタンド2における評価関数の成分である。
なお、板厚偏差が利用できる場合には、Jerriを以下のようにするとよい。
Here, J erri is a component of the evaluation function in the #i rolling stand 2.
In addition, when plate | board thickness deviation can be utilized, it is good to set Jerri as follows.

Figure 2009208115
Figure 2009208115

この場合には、板厚偏差を直接評価することができる。
式(13),式(14)において、wpi,whiは重み係数であり、ΔPi_actは実際の圧延機から得られた、減速前の一定速度時からの圧延荷重変化の実績値、ΔPi_simは圧延動的モデル部6の計算で得られる、減速前の一定速度時からの圧延荷重変化である。Δhi_actは#i圧延スタンド2の、減速前の一定速度時からの実績出側板厚変化であり、Δhi_simは圧延動的モデル部6の計算で得られる、減速前の一定速度時からの出側板厚変化である。
In this case, the thickness deviation can be directly evaluated.
In Expressions (13) and (14), w pi and w hi are weighting factors, and ΔP i _act is the actual value of the change in rolling load from a constant speed before deceleration obtained from an actual rolling mill. , ΔP i — sim is a change in rolling load from a constant speed before deceleration, obtained by calculation of the rolling dynamic model unit 6. Δh i _act is the actual thickness change of the #i rolling stand 2 from the constant speed before deceleration, and Δh i _sim is the constant speed before deceleration obtained by the calculation of the rolling dynamic model unit 6. This is the change in the thickness of the outlet side from the time.

sは評価の開始時間、teは評価の終了時間である。例えば、図4に示す如く、減速時であれば、tsを減速開始時とし、teを減速終了後において圧延材3間の溶接点が圧延スタンド2を通過し始める時とすることができる。
なお、加速時であれば、tsを加速開始前における圧延材3間の溶接点が圧延スタンド2を通過し終えた時とし、teを加速終了後一定速度となり圧延状態が整定した時とすることができる。
上述した手法は例示であり、これに限定はされない。例えば、加速時及び減速時を同時に評価する場合には、加速時に対する評価関数と減速時に対する評価関数とを加えたものを評価関数とすることができる。この場合には、評価する時間が増えるため、さらなる圧延状態パラメータの精度向上が期待できる。また、実績板厚としては、板厚計から直接測定されたものだけではなく、マスフロー板厚のように測定値から算出されたものを含む。さらに、評価関数の形は式(12)〜式(14)に限定されるものではなく、圧延荷重の計算値と実績値との差、すなわち、Pi_act−Pi_simを含んでいればよい。板厚についても、板厚の計算値と実績値との差、すなわち、hi_act−hi_simを含んでいればよい。例えば、一定速度時にPi_act=Pi_simの場合、式(13)の括弧内をPi_act−Pi_simとしてもよい。一定速度時にhi_act=hi_simの場合、式(14)の第2項括弧内をhi_act−hi_simとしてもよい。
t s is the evaluation start time, and t e is the evaluation end time. For example, as shown in FIG. 4, at the time of deceleration, t s can be set as the start of deceleration, and t e can be set as the time when the welding point between the rolled materials 3 starts to pass through the rolling stand 2 after the end of deceleration. .
Incidentally, if during acceleration, when the welding points between the strip 3 in the acceleration before starting the t s is the time it has finished passing through the rolling stand 2, the rolling state becomes acceleration end after a fixed rate t e has settled and can do.
The method described above is an example, and the present invention is not limited to this. For example, when evaluating at the time of acceleration and deceleration at the same time, an evaluation function can be obtained by adding an evaluation function for acceleration and an evaluation function for deceleration. In this case, since the time for evaluation increases, further improvement in accuracy of the rolling state parameter can be expected. In addition, the actual plate thickness includes not only those measured directly from the plate thickness meter but also those calculated from measured values such as mass flow plate thickness. Furthermore, the shape of the evaluation function is not limited Equation (12) to Formula (14), the difference between the calculated and actual values of the rolling load, i.e., contain P i _ act -P i _ sim It only has to be. The thickness is also a difference between the calculated and actual values of the thickness, i.e., it needs to include a h i _ act -h i _ sim . For example, when P i _act = P i _sim at a constant speed, the parentheses in the equation (13) may be P i _act− P i _sim . For h i _ act = h i _ sim at a constant speed, the second term in the parentheses of formula (14) may be h i _ act -h i _ sim .

以上述べた手順で調整された摩擦係数を基に圧延動的モデルを用いて、圧延制御パラメータを最適化し算出するようにしている。
言い換えるならば、圧延状態パラメータ最適化部8が算出した摩擦係数μiを圧延モデル式中の摩擦係数μiに適用し、圧延機モデル部4と、板厚制御モデル部5と、圧延制御パラメータ最適化部7での計算を行って、最適化された圧延制御パラメータを算出するようにしている。
図7には、本発明における摩擦係数等の圧延状態パラメータの調整方法、調整された圧延状態パラメータを用いた圧延制御パラメータの最適化方法のフローチャートが示されている。以降、圧延状態パラメータが摩擦係数である場合を例示し説明を行うが、本発明では、圧延状態パラメータとして変形抵抗が採用可能であるし、摩擦係数と変形抵抗との両者を採用可能である。
The rolling control parameters are optimized and calculated using the rolling dynamic model based on the friction coefficient adjusted by the procedure described above.
In other words, the friction coefficient mu i of the rolling condition parameter optimization unit 8 is calculated by applying the friction coefficient mu i during rolling model equation, a rolling mill model section 4, the sheet thickness control model section 5, the rolling control parameter The optimized rolling control parameter is calculated by performing the calculation in the optimization unit 7.
FIG. 7 shows a flowchart of a method for adjusting a rolling state parameter such as a friction coefficient and a method for optimizing a rolling control parameter using the adjusted rolling state parameter in the present invention. Hereinafter, the case where the rolling state parameter is a friction coefficient will be exemplified and described. However, in the present invention, deformation resistance can be adopted as the rolling state parameter, and both the friction coefficient and deformation resistance can be adopted.

まず、摩擦係数等の圧延状態パラメータの初期値として、実際の圧延機のデータから推定または学習した値をセットする。[S71]
次に、#1圧延スタンド2の1ケース目の計算として、圧延機モデル部4、板厚制御モデル部5を用いて、S71でセットした摩擦係数を基に、時刻tsからteまでを含む時間において圧延のシミュレーションを行う。[S72〜S74]
それらで得られた結果(圧延荷重、板厚、ロールギャップなど)を基に、評価関数(式(12))を計算する。[S75]
次に、2ケース目の摩擦係数を圧延状態パラメータ最適化部8が決定し、2ケース目の圧延シミュレーションを行い、評価関数を計算する。[S77,S73〜S75]
全ケースについて評価関数を計算した後、圧延状態パラメータ最適化部8が最小の評価関数を与える摩擦係数パラメータ(摩擦係数に関するパラメータ)α1,Vout,10を決定する。そして、以降の最適化のためのシミュレーションに用いる最適な摩擦係数のパラメータを得る。[S76,S78]
同様に、順に#2〜#5圧延スタンド2についても摩擦係数パラメータを決定し、新たな摩擦係数を算出する。#5圧延スタンド2については、α5,Vout,50,cμ5の摩擦係数パラメータを決定する。[S79,S710,S73〜S78]
このように、圧延スタンド毎に上流側から最適化を行うことで、圧延材の移送による上流側圧延スタンド2から下流側圧延スタンド2への圧延荷重や板厚等の影響を精確に再現でき、得られた圧延状態パラメータの精度が向上する。また、全スタンドの圧延状態パラメータを一括で最適化する場合よりも、離散的な最適化の場合はシミュレーションのケース数が少なく、連続的な最適化の場合は収束が速くなる傾向があり、最適化の計算時間が短くなるという効果もある。
First, a value estimated or learned from actual rolling mill data is set as an initial value of rolling state parameters such as a friction coefficient. [S71]
Then, as a case-th calculation of # 1 rolling stand 2, the rolling mill model section 4, with a thickness control model section 5, based on the coefficient of friction is set at S71, from the time t s until t e Simulation of rolling is performed in the time including. [S72 to S74]
An evaluation function (formula (12)) is calculated based on the results obtained (rolling load, sheet thickness, roll gap, etc.). [S75]
Next, the rolling condition parameter optimization unit 8 determines the friction coefficient of the second case, performs the rolling simulation of the second case, and calculates the evaluation function. [S77, S73 to S75]
After calculating the evaluation function for all cases, the rolling state parameter optimization unit 8 determines the friction coefficient parameters (parameters related to the friction coefficient) α 1 , V out, 10 that give the minimum evaluation function. And the parameter of the optimal friction coefficient used for the simulation for subsequent optimization is obtained. [S76, S78]
Similarly, a friction coefficient parameter is determined for the # 2 to # 5 rolling stands 2 in order, and a new friction coefficient is calculated. For the # 5 rolling stand 2, the friction coefficient parameters of α 5 , V out, 50 , cμ 5 are determined. [S79, S710, S73 to S78]
Thus, by performing optimization from the upstream side for each rolling stand, it is possible to accurately reproduce the influence of rolling load, sheet thickness, etc. from the upstream rolling stand 2 to the downstream rolling stand 2 due to the transfer of the rolled material, The accuracy of the obtained rolling state parameters is improved. In addition, the number of simulation cases is smaller in the case of discrete optimization than in the case of optimizing the rolling state parameters for all stands, and the convergence tends to be faster in the case of continuous optimization. There is also an effect that the calculation time of the conversion is shortened.

次に、終了条件の判断[S711]であるが、#1〜#5圧延スタンド2の摩擦係数パラメータの調整を少なくとも2回行ったか否かで判定している。これは、張力の影響や制御系の構成の影響による下流側スタンドの摩擦係数の変更により、上流側に影響がある場合に、2回目の調整で下流側からの影響が考慮された最適化が行われることとなり、有効である。終了条件を、パラメータの変更が無くなったとき、または、許容値以下になったときとしてもよい。評価関数の変化が許容値以下になったときや所定の回数(複数回)、例えば、5回目のパラメータ調整で終了としてもよい。   Next, the end condition determination [S711] is made based on whether or not the friction coefficient parameters of the # 1 to # 5 rolling stands 2 have been adjusted at least twice. This is because when the upstream side is affected by the change in the friction coefficient of the downstream stand due to the influence of the tension or the control system configuration, the second adjustment will be optimized in consideration of the influence from the downstream side. It will be done and is effective. The end condition may be set when the parameter is no longer changed or when the parameter is less than the allowable value. When the change in the evaluation function becomes equal to or less than the allowable value or a predetermined number of times (a plurality of times), for example, the fifth parameter adjustment may be ended.

終了条件を満たした場合には、摩擦係数パラメータが確定したものとして、それを基に、圧延機モデル部4、板厚制御モデル部5、圧延制御パラメータ最適化部7を用いて、圧延のシミュレーションを用いた最適化を行う。[S712]
かかる圧延のシミュレーションにより得られた圧延制御パラメータは、各圧延スタンド2における最適化された摩擦係数を考慮したものであるため、この圧延制御パラメータを実際の圧延機に適用することで、圧延材3の加減速時に圧延材3の形状を悪化させることなく板厚精度を向上させることができる。[S713]
また、以上述べた圧延シミュレーションを行うことで、実際の圧延機を用いた実験を行うことなく、圧延材3の加減速時における最適な圧延制御パラメータを得ることができるようになる。
When the end condition is satisfied, it is assumed that the friction coefficient parameter has been determined. Based on this, the rolling simulation is performed using the rolling mill model unit 4, the plate thickness control model unit 5, and the rolling control parameter optimization unit 7. Perform optimization using. [S712]
Since the rolling control parameter obtained by the simulation of rolling takes into account the optimized friction coefficient in each rolling stand 2, the rolling material 3 can be obtained by applying this rolling control parameter to an actual rolling mill. The plate thickness accuracy can be improved without deteriorating the shape of the rolled material 3 during acceleration / deceleration. [S713]
In addition, by performing the rolling simulation described above, it is possible to obtain optimum rolling control parameters at the time of acceleration / deceleration of the rolled material 3 without performing an experiment using an actual rolling mill.

#1〜#5圧延スタンド2の摩擦係数の最適化は、上述のような離散量に対してだけでなく連続量に対してでもよいし、上述のような総当り法以外に他の最適化手法でもよい。例えば、「最適設計ハンドブック,山川宏編,朝倉書店,2003」のp8〜p63に記載されている非線形計画法、分枝限定法、シミュレーテェッドアニーリング法、遺伝的アルゴリズム等の発見的手法でもよい。
図8には、全ての摩擦パラメータαi,Vout,i0,α5,Vout,50,cμ5を連続量とし、逐次2次計画法等の非線形計画法を用いて最適化する場合のフローチャートを示す。
The optimization of the friction coefficient of the # 1 to # 5 rolling stands 2 may be performed not only on the discrete amount as described above but also on the continuous amount, and other optimizations besides the brute force method as described above A technique may be used. For example, heuristic techniques such as nonlinear programming, branch and bound, simulated annealing, genetic algorithm, etc. described in p8-p63 of “Optimum Design Handbook, Hiroshi Yamakawa, Asakura Shoten, 2003” may be used. .
FIG. 8 shows a case in which all friction parameters α i , V out, i 0 , α 5 , V out, 50 , cμ 5 are set as continuous quantities and optimized using nonlinear programming such as sequential quadratic programming. A flowchart is shown.

評価関数としては、全圧延スタンド2,2,・・・での圧延荷重の評価を行う式(12),式(13)を用いている。
まず、摩擦係数の初期値として、実際の圧延機のデータから推定または学習した値をセットする。[S81]
次に、全ての圧延スタンド2,2,・・・について、セットした摩擦係数を用いて、時刻tsからteまで、圧延機モデル部4、板厚制御モデル部5、圧延制御パラメータ最適化部7を用いて、圧延シミュレーションを行う。[S82〜S84]
それらで得られた結果(圧延荷重、板厚、ロールギャップなど)を基に、評価関数(式(12))を計算する。[S85]
得られた評価関数が最小値に収束しているか否かを判定し、収束していない場合は、次の摩擦係数を圧延状態パラメータ最適化部8が決定し、再び圧延のシミュレーションを行い、評価関数を計算する。[S87,S82〜S85]
なお、逐次2次計画法(SQP)のような公知の非線形最適化手法を用いる際には、S35〜S87のステップが、公知である非線形最適化手法のアルゴリズム中で実行される。また、S87中で、次の圧延状態パラメータの探索のためにシミュレーションが行われる場合もある。
As the evaluation function, Expressions (12) and (13) for evaluating the rolling load at all rolling stands 2, 2,... Are used.
First, a value estimated or learned from actual rolling mill data is set as the initial value of the friction coefficient. [S81]
Then, all the rolling stands 2, 2, about ..., using the set coefficient of friction, from time t s to t e, mill model section 4, the plate thickness control model section 5, the rolling control parameter optimization Using part 7, a rolling simulation is performed. [S82 to S84]
An evaluation function (formula (12)) is calculated based on the results obtained (rolling load, sheet thickness, roll gap, etc.). [S85]
It is determined whether or not the obtained evaluation function has converged to the minimum value. If it has not converged, the next friction coefficient is determined by the rolling state parameter optimizing unit 8, and the rolling simulation is performed again for evaluation. Calculate the function. [S87, S82 to S85]
When a known nonlinear optimization method such as sequential quadratic programming (SQP) is used, steps S35 to S87 are executed in an algorithm of a known nonlinear optimization method. In S87, a simulation may be performed to search for the next rolling state parameter.

得られた評価関数が最小値に収束している場合は、得られた摩擦係数パラメータαi,Vout,i0を基に、圧延機モデル部4、板厚制御モデル部5、圧延制御パラメータ最適化部7を用いて、圧延のシミュレーションを行う。[S89]
かかる圧延のシミュレーションにより得られた圧延制御パラメータは、各圧延スタンド2における摩擦係数の調整を考慮したものであるため、この圧延制御パラメータを実際の圧延機に適用することで、圧延材3の加減速時に圧延材3の形状を悪化させることなく板厚精度を向上させることができる。[S810]
一方、圧延材3を加速及び/又は減速するに際し、複数の圧延スタンド2でのロールギャップの時間変化と実績ロールギャップの時間変化との差が最小となるように、各圧延スタンド2での摩擦係数に関する摩擦係数パラメータを最適化し、最適化された摩擦係数パラメータが組み込まれた圧延動的モデルを用いて、圧延制御パラメータを算出することもできる。
When the obtained evaluation function has converged to the minimum value, the rolling mill model unit 4, the plate thickness control model unit 5, and the optimum rolling control parameters are based on the obtained friction coefficient parameters α i , V out, i0. A rolling simulation is performed using the conversion unit 7. [S89]
Since the rolling control parameters obtained by the rolling simulation take into account the adjustment of the friction coefficient in each rolling stand 2, the rolling control parameters can be applied to the rolling mill 3 by applying the rolling control parameters to an actual rolling mill. The plate thickness accuracy can be improved without deteriorating the shape of the rolled material 3 during deceleration. [S810]
On the other hand, when accelerating and / or decelerating the rolled material 3, the friction at each rolling stand 2 is minimized so that the difference between the time change of the roll gaps at the plurality of rolling stands 2 and the time change of the actual roll gaps is minimized. It is also possible to calculate the rolling control parameter using a rolling dynamic model in which the friction coefficient parameter related to the coefficient is optimized and the optimized friction coefficient parameter is incorporated.

その場合、#i圧延スタンド2における評価関数の成分Jerriは、次式のようにするとよい。 In that case, the component J erri of the evaluation function in the #i rolling stand 2 is preferably expressed by the following equation.

Figure 2009208115
Figure 2009208115

また、板厚が計測又は計算できる圧延スタンド2については、次式のようにしてもよい。   Moreover, about the rolling stand 2 which can measure or calculate plate | board thickness, you may make it like following Formula.

Figure 2009208115
Figure 2009208115

式(15),式(16)で、wsi,whiは重み係数である。Δsi_actは、実際の圧延機から得られた、減速前の一定速度時からのロールギャップ変化の実績値であり、Δsi_simはシミュレーションでの、減速前の一定速度時からのロールギャップ変化である。Δhi_actは#i圧延スタンド2の、減速前の一定速度時からの実績出側板厚変化であり、Δhi_simはシミュレーションでの、減速前の一定速度時からの出側板厚変化である。tsは評価の開始時間、teは評価の終了時間である。なお、評価関数の形は、式(12),式(15),式(16)に限定されるものではなく、ロールギャップの計算値と実績値との差、すなわち、si_act−si_simを含んでいればよい。例えば、一定速度時にsi_act=si_simの場合、式(15),式(16)の括弧内をsi_act−si_simとしてもよい。 In equations (15) and (16), w si and w hi are weighting factors. Δs i _ act was obtained from the actual rolling mill, is the actual value of the roll gap change from the time of constant speed before the deceleration, Δs i _ sim is in the simulation, roll from the time of constant speed before the deceleration It is a gap change. Δh i _act is the actual change in the delivery side thickness of the #i rolling stand 2 from the constant speed before deceleration, and Δh i _sim is the change in the delivery side thickness from the constant speed before deceleration in the simulation. is there. t s is the evaluation start time, and t e is the evaluation end time. The form of the evaluation function is not limited to Expression (12), Expression (15), and Expression (16), and the difference between the calculated value of the roll gap and the actual value, that is, s i _act −s. i _ sim may if it contains. For example, in the case of s i _act = s i _sim at a constant speed, the parentheses in the equations (15) and (16) may be set as s i _act −s i _sim .

評価関数の式(13),式(14)と評価関数の式(15),式(16)とを同時に用いることもできる。例えば、圧延スタンド2ごとの評価関数の成分として、全スタンドで「圧延荷重とロールギャップ」又は「圧延荷重とロールギャップと板厚」を採用することが可能であるし、ある圧延スタンド2では「圧延荷重とロールギャップ」を採用し、別の圧延スタンド2では、「圧延荷重とロールギャップと板厚」を採用することが可能である。圧延荷重とロールギャップの両方を採用することにより、利用する情報量が増えることになり、最適化される摩擦係数や変形抵抗の精度がさらに向上する。   The evaluation function equations (13) and (14) and the evaluation function equations (15) and (16) can be used simultaneously. For example, “rolling load and roll gap” or “rolling load and roll gap and plate thickness” can be adopted for all the stands as components of the evaluation function for each rolling stand 2. "Rolling load and roll gap" are adopted, and in another rolling stand 2, "rolling load, roll gap and plate thickness" can be adopted. By adopting both the rolling load and the roll gap, the amount of information to be used is increased, and the accuracy of the optimized friction coefficient and deformation resistance is further improved.

以下、本発明の技術を用いて圧延シミュレーションを行った結果を示す。
図9は、減速時の圧延実績データであって、本発明の技術を採用していない従来例である。減速時に、圧延荷重偏差が発生すると共に仕上板厚偏差が生じていることがわかる。
図9の実績データと式(12),式(14)の評価関数とを用いて、#1〜#4圧延スタンド2については12通り、#5圧延スタンド2については、36通りの場合について離散的な最適化計算である総当たり法を実施した。
結果として、次の最適化された摩擦係数パラメータが得られた。
Hereafter, the result of having performed the rolling simulation using the technique of this invention is shown.
FIG. 9 shows rolling performance data at the time of deceleration, which is a conventional example that does not employ the technique of the present invention. It can be seen that a rolling load deviation and a finished sheet thickness deviation occur during deceleration.
Using the performance data in FIG. 9 and the evaluation functions of the equations (12) and (14), 12 types are used for the # 1 to # 4 rolling stands 2 and 36 types are used for the # 5 rolling stand 2 in a discrete manner. The brute force method, which is a typical optimization calculation.
As a result, the following optimized coefficient of friction parameter was obtained.

・α1=0.5,Vout,10=600(mpm)
・α2=0.5,Vout,20=600(mpm)
・α3=0.9,Vout,30=600(mpm)
・α4=1.1,Vout,40=800(mpm)
・α5=0.7,Vout,50=1000(mpm),cμ5=2.5

図10は、最適化された摩擦係数パラメータを用いたシミュレーション結果、すなわち、本発明の技術を適用した圧延シミュレーションの結果である。
・ Α 1 = 0.5, V out, 10 = 600 (mpm)
・ Α 2 = 0.5, V out, 20 = 600 (mpm)
・ Α 3 = 0.9, V out, 30 = 600 (mpm)
・ Α 4 = 1.1, V out, 40 = 800 (mpm)
Α 5 = 0.7, V out, 50 = 1000 (mpm), cμ 5 = 2.5

FIG. 10 shows a simulation result using the optimized coefficient of friction parameter, that is, a rolling simulation result to which the technique of the present invention is applied.

図9と図10とを比較すると、圧延荷重偏差の実績(図9(b))とシミュレーション(図10(b))がほぼ一致するように摩擦係数パラメータが調整されていることがわかる。仕上板厚偏差に関しては、仕上板厚偏差の実績(図9(c))とシミュレーション(図10(c))とにおける板厚変動のピーク発生時刻がほぼ一致しており、加えて、実績データのノイズ等の外乱を考慮すると、ピーク値もほぼ一致していることが確認できる。
次に、最適化された摩擦係数パラメータを用いて、圧延制御パラメータとして#5圧延スタンド2のロールギャップ加減速補償量を最適化する。本実施形態では、#5圧延スタンド2のロール速度加減速補償量を変更するやり方について精説を行ってはいないが、ロール速度加減速補償量を最適化変数に加えることで、ロールギャップ加減速補償量と同様に最適化可能である。
Comparing FIG. 9 and FIG. 10, it can be seen that the friction coefficient parameter is adjusted so that the actual rolling load deviation (FIG. 9 (b)) and the simulation (FIG. 10 (b)) substantially coincide. With respect to the finish plate thickness deviation, the peak occurrence time of the plate thickness variation in the finish plate thickness deviation record (FIG. 9 (c)) and the simulation (FIG. 10 (c)) is almost the same. In consideration of disturbances such as noise, it can be confirmed that the peak values are almost the same.
Next, using the optimized friction coefficient parameter, the roll gap acceleration / deceleration compensation amount of the # 5 rolling stand 2 is optimized as a rolling control parameter. In the present embodiment, the method of changing the roll speed acceleration / deceleration compensation amount of the # 5 rolling stand 2 is not elaborated, but the roll gap acceleration / deceleration is added by adding the roll speed acceleration / deceleration compensation amount to the optimization variable. It can be optimized in the same way as the compensation amount.

評価関数Jerr’は、式(17)のように定義する。 The evaluation function J err ′ is defined as in Expression (17).

Figure 2009208115
Figure 2009208115

式(17)において、Δh5_simは、既述の通り、圧延動的モデル部6で計算された、減速前の一定速度からの#5圧延スタンド2の出側板厚変化、すなわち仕上板厚変化である。tsは評価の開始時間、teは評価の終了時間である。例えば、圧延状態パラメータの最適化と同様に、図4に示す如く、減速時であれば、tsを減速開始時とし、teを減速終了後において圧延材3間の溶接点が圧延スタンド2を通過し始める時とすることができる。加速時であれば、tsを加速開始前における圧延材3間の溶接点が圧延スタンド2を通過し終えた時とし、teを加速終了後一定速度となり圧延状態が整定した(一定となった)時とすることができる。 In Expression (17), Δh 5 _sim is the change in the exit side plate thickness of the # 5 rolling stand 2 from the constant speed before deceleration calculated by the rolling dynamic model unit 6 as described above, that is, the finish plate thickness It is a change. t s is the evaluation start time, and t e is the evaluation end time. For example, as in the optimization of the rolling state parameter, as shown in FIG. 4, at the time of deceleration, t s is set at the start of deceleration, and t e is the welding point between the rolled materials 3 after the completion of deceleration at the rolling stand 2. It can be when it begins to pass. At the time of acceleration, t s is the time when the welding point between the rolled materials 3 before the start of acceleration has finished passing through the rolling stand 2, and t e is a constant speed after the end of acceleration, and the rolling state is stabilized (becomes constant). It can be time).

一例として、#5圧延スタンド2のロールギャップの加減速補償量を次式で決定する。

・Vout,5≦100(mpm): Δs5=Δs50
・100≦Vout,5≦Vout,50
Δs5=Δs50・(Vout,5−Vout,502 /(Vout,50−100)2
・Vout,50≦Vout,5: Δs5=0

Δs50とVout,50の2つのパラメータが最適化すべき圧延制御パラメータである。
As an example, the acceleration / deceleration compensation amount of the roll gap of the # 5 rolling stand 2 is determined by the following equation.

Vout, 5 ≦ 100 (mpm): Δs 5 = Δs 50
・ 100 ≦ V out, 5 ≦ V out, 50 :
Δs 5 = Δs 50 · (V out, 5 −V out, 50 ) 2 / (V out, 50 −100) 2
・ V out, 50 ≦ V out, 5 : Δs 5 = 0

Two parameters of Δs 50 and V out, 50 are rolling control parameters to be optimized.

簡単な一例として、Δs50=−60,−120(μm)の2通り、Vout,50=600,800(μm)の2通りの計4通りで離散的な最適化を行う。
総当たり法で評価関数を計算することで、Δs50=−120(μm)、Vout,50=800(μm)の解を得た。
図11は、上記の如く最適化された圧延制御パラメータΔs50,Vout,50を用いて、圧延シミュレーションを行った結果である。図10と図11との比較から、減速時の板厚変動が半減していることは明らかである。
As a simple example, discrete optimization is performed in two ways: two types of Δs 50 = −60, −120 (μm) and two types of V out, 50 = 600,800 (μm).
By calculating the evaluation function using the brute force method, a solution of Δs 50 = −120 (μm) and V out, 50 = 800 (μm) was obtained.
FIG. 11 shows the result of rolling simulation using the rolling control parameters Δs 50 and V out, 50 optimized as described above. From the comparison between FIG. 10 and FIG. 11, it is clear that the plate thickness variation during deceleration is halved.

以上、本発明にかかる圧延制御パラメータの算出方法及び算出装置、圧延シミュレーション装置は、上述した実施の形態に限定されるものではない。
すなわち、圧延機9は#1〜#5圧延スタンド2を有するものを例示して説明したが、圧延スタンド2の数はこれに限定されず、リバースミルのようなシングルスタンドミルにも適用可能である。また、圧延形態も冷間圧延に限られるものではなく、熱間圧延にも適用可能である。
また、本実施形態の場合、各圧延スタンド2における摩擦係数に関するパラメータを調整したが、変形抵抗に関するパラメータを調整するようにしてもよい。例えば、既述のように、「板圧延の理論と実際、日本鉄鋼協会編、1984」のp188〜189に開示された式(7.51)〜式(7.53)のいずれかを用いることとし、その変形抵抗をki_nominalとする。そして、変形抵抗を調整する場合の変形抵抗モデルとしては、例えば、

i=ki_nominal・(β1i+β2・Vout,i

を用いることができる(ただし、β1i,β2は調整すべき圧延状態パラメータ)。
As described above, the rolling control parameter calculation method, calculation apparatus, and rolling simulation apparatus according to the present invention are not limited to the above-described embodiments.
That is, although the rolling mill 9 has been described by exemplifying the one having the # 1 to # 5 rolling stands 2, the number of the rolling stands 2 is not limited to this and can be applied to a single stand mill such as a reverse mill. is there. Further, the rolling form is not limited to cold rolling, and is applicable to hot rolling.
Moreover, in the case of this embodiment, although the parameter regarding the friction coefficient in each rolling stand 2 was adjusted, you may make it adjust the parameter regarding deformation resistance. For example, as described above, any one of formulas (7.51) to (7.53) disclosed in p188 to 189 of "Theory and practice of plate rolling, edited by Japan Iron and Steel Institute, 1984" is used. And let the deformation resistance be k i _ nominal . And as a deformation resistance model when adjusting the deformation resistance, for example,

k i = k i _ nominal · (β1 i + β2 · V out, i )

Can be used (where β1 i and β2 are rolling state parameters to be adjusted).

なお、変形抵抗に関するパラメータと摩擦係数パラメータとを同時に最適化することもできる。この場合には、調整すべきパラメータが増加し、さらに細かな調整も可能となるため、圧延動的モデルの精度を高めることができる。
また、評価関数の形も、式(12)〜式(16)に限られるものではなく、圧延スタンド毎に個別の評価関数を採用してもよい。
It is also possible to simultaneously optimize the parameter relating to the deformation resistance and the friction coefficient parameter. In this case, the parameter to be adjusted increases and finer adjustment is possible, so that the accuracy of the rolling dynamic model can be increased.
Further, the form of the evaluation function is not limited to Expression (12) to Expression (16), and an individual evaluation function may be adopted for each rolling stand.

圧延制御パラメータの算出装置及び圧延シミュレーション装置のブロック図である。It is a block diagram of the calculation apparatus of a rolling control parameter, and a rolling simulation apparatus. 圧延機モデル部がモデル化している圧延機の模式図である。It is a schematic diagram of the rolling mill which the rolling mill model part has modeled. 板厚張力制御部の入出力を示した図である。It is the figure which showed the input / output of the plate | board thickness tension control part. #5圧延スタンドにおける圧延速度パターンを示した図である。It is the figure which showed the rolling speed pattern in a # 5 rolling stand. ロールギャップの補正量(加減速補償)を示した図である。It is the figure which showed the correction amount (acceleration / deceleration compensation) of the roll gap. ロール速度の補正量(加減速補償)を示した図である。It is the figure which showed the correction amount (acceleration / deceleration compensation) of roll speed. 圧延制御パラメータの算出装置及び圧延シミュレーション装置における計算フローチャートである(離散量の場合)。It is a calculation flowchart in the calculation apparatus of a rolling control parameter, and a rolling simulation apparatus (in the case of a discrete amount). 圧延制御パラメータの算出装置及び圧延シミュレーション装置における計算フローチャートである(連続量の場合)It is a calculation flowchart in the calculation apparatus of a rolling control parameter, and a rolling simulation apparatus (in the case of continuous quantity). 圧延実績を示した図である(従来例)。It is the figure which showed the rolling performance (conventional example). 最適化した摩擦係数を用いて圧延状況をシミュレーションした結果を示す図である。It is a figure which shows the result of having simulated the rolling condition using the optimized friction coefficient. 最適化した加減速補償量を用いて圧延状況をシミュレーションした結果を示す図である。It is a figure which shows the result of having simulated the rolling condition using the optimized acceleration / deceleration compensation amount.

符号の説明Explanation of symbols

1 圧延制御パラメータの算出装置
2 圧延スタンド
3 圧延材
4 圧延機モデル部
5 板厚制御モデル部
6 圧延動的モデル部
7 圧延制御パラメータ最適化部
8 圧延状態パラメータ最適化部
9 圧延機
10 目標板厚設定部
11 目標張力設定部
12 板厚張力制御部
13 加減速補償部
14 圧延速度パターン決定部
DESCRIPTION OF SYMBOLS 1 Rolling control parameter calculation apparatus 2 Rolling stand 3 Rolled material 4 Rolling mill model part 5 Sheet thickness control model part 6 Rolling dynamic model part 7 Rolling control parameter optimization part 8 Rolling state parameter optimization part 9 Rolling machine 10 Target plate Thickness setting unit 11 Target tension setting unit 12 Plate thickness tension control unit 13 Acceleration / deceleration compensation unit 14 Rolling speed pattern determination unit

Claims (6)

圧延材と該圧延材を圧延する複数の圧延スタンドと該圧延スタンドの板厚制御を行う制御系とを動的に表現する圧延動的モデルを構築し、
前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデルでの計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化し、
最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデルを用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出することを特徴とする圧延制御パラメータの算出方法。
Build a rolling dynamic model that dynamically expresses a rolled material, a plurality of rolling stands for rolling the rolled material, and a control system for controlling the thickness of the rolling stand,
Evaluation including the difference between the calculated value and the actual value in the rolling dynamic model with respect to the time change of the rolling load and / or the time change of the roll gap at the rolling stand at the time of acceleration and / or deceleration of the rolled material. Optimize parameters related to coefficient of friction and / or deformation resistance at each rolling stand so that the function is minimized,
Optimizing rolling control parameters including the setting parameters of the plate thickness control and / or the setting parameters of the rolling stand using the rolling dynamic model incorporating the parameters relating to the optimized coefficient of friction and / or parameters relating to deformation resistance A method for calculating a rolling control parameter, characterized by:
前記摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータの最適化の計算を、上流側の圧延スタンドから下流側の圧延スタンドへ順に行うことを特徴とする請求項1に記載の圧延制御パラメータの算出方法。   The calculation method of the rolling control parameter according to claim 1, wherein calculation of optimization of the parameter relating to the friction coefficient and / or the parameter relating to deformation resistance is sequentially performed from the upstream rolling stand to the downstream rolling stand. . 前記摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータの最適化の計算を、圧延スタンド毎に1回又は複数回実施することを特徴とする請求項1又は2に記載の圧延制御パラメータの算出方法。   The calculation method of the rolling control parameter according to claim 1 or 2, wherein the calculation of the parameter relating to the friction coefficient and / or the parameter relating to the deformation resistance is performed once or plural times for each rolling stand. 前記摩擦係数に関するパラメータは、前記摩擦係数が変化しなくなる上限の圧延速度である上限圧延速度と、圧延速度0での摩擦係数を規定する変数とを備えていることを特徴とする請求項1〜3のいずれかに記載の圧延制御パラメータの算出方法。   The parameter relating to the friction coefficient includes an upper limit rolling speed that is an upper limit rolling speed at which the friction coefficient does not change, and a variable that defines a friction coefficient at a rolling speed of 0. 4. A method for calculating a rolling control parameter according to any one of 3 above. 圧延材と該圧延材を圧延する複数の圧延スタンドとをモデル化した圧延機モデル部、及び前記圧延スタンドの板厚制御を行う板厚制御モデル部を備えた圧延動的モデル部と、
前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデル部での計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化する圧延状態パラメータ最適化部と、
前記圧延状態パラメータ最適化部により最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデル部を用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出する圧延制御パラメータ最適化部と、
を有することを特徴とする圧延制御パラメータの算出装置。
A rolling mill model unit that models a rolled material and a plurality of rolling stands that roll the rolled material, and a rolling dynamic model unit that includes a plate thickness control model unit that controls the plate thickness of the rolling stand,
Regarding the time change of the rolling load and / or the time change of the roll gap in the rolling stand at the time of acceleration and / or deceleration of the rolled material, the difference between the calculated value and the actual value in the rolling dynamic model part is included. A rolling state parameter optimizing unit that optimizes a parameter relating to a coefficient of friction and / or a parameter relating to deformation resistance in each rolling stand so as to minimize the evaluation function;
Using the rolling dynamic model unit in which the parameters relating to the friction coefficient and / or the parameters relating to deformation resistance optimized by the rolling state parameter optimizing unit are incorporated, the setting parameters for the plate thickness control and / or the rolling stand A rolling control parameter optimization unit that optimizes and calculates rolling control parameters including setting parameters;
An apparatus for calculating a rolling control parameter, comprising:
圧延材と該圧延材を圧延する複数の圧延スタンドとをモデル化した圧延機モデル部、及び前記圧延スタンドの板厚制御を行う板厚制御モデル部を備えた圧延動的モデル部と、
前記圧延材の加速時及び/又は減速時における前記圧延スタンドでの圧延荷重の時間変化及び/又はロールギャップの時間変化に関し、前記圧延動的モデル部での計算値と実績値との差を含む評価関数が最小となるように、各圧延スタンドでの摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータを最適化する圧延状態パラメータ最適化部と、
前記圧延状態パラメータ最適化部により最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータが組み込まれた前記圧延動的モデル部を用いて、前記板厚制御の設定パラメータ及び/又は圧延スタンドの設定パラメータを含む圧延制御パラメータを最適化し算出する圧延制御パラメータ最適化部と、を有していて、
前記圧延動的モデル部は、圧延状態パラメータ最適化部で最適化された摩擦係数に関するパラメータ及び/又は変形抵抗に関するパラメータと、圧延制御パラメータ最適化部で最適化された圧延制御パラメータとを用いて、圧延の状況を再現できるように構成されていることを特徴とする圧延シミュレーション装置。
A rolling mill model unit that models a rolled material and a plurality of rolling stands that roll the rolled material, and a rolling dynamic model unit that includes a thickness control model unit that controls the thickness of the rolling stand,
Regarding the time change of the rolling load and / or the time change of the roll gap in the rolling stand at the time of acceleration and / or deceleration of the rolled material, the difference between the calculated value and the actual value in the rolling dynamic model part is included. A rolling state parameter optimizing unit that optimizes a parameter related to a coefficient of friction and / or a parameter related to deformation resistance in each rolling stand so that the evaluation function is minimized;
Using the rolling dynamic model unit in which the parameters relating to the friction coefficient and / or the parameters relating to deformation resistance optimized by the rolling state parameter optimizing unit are incorporated, the setting parameters for the plate thickness control and / or the rolling stand A rolling control parameter optimization unit that optimizes and calculates rolling control parameters including setting parameters, and
The rolling dynamic model unit uses the parameters related to the friction coefficient and / or the deformation resistance optimized by the rolling state parameter optimization unit and the rolling control parameters optimized by the rolling control parameter optimization unit. A rolling simulation apparatus characterized by being configured to reproduce a rolling situation.
JP2008053420A 2008-03-04 2008-03-04 Method and device for calculating parameter of rolling control, and rolling simulation device Pending JP2009208115A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008053420A JP2009208115A (en) 2008-03-04 2008-03-04 Method and device for calculating parameter of rolling control, and rolling simulation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008053420A JP2009208115A (en) 2008-03-04 2008-03-04 Method and device for calculating parameter of rolling control, and rolling simulation device

Publications (1)

Publication Number Publication Date
JP2009208115A true JP2009208115A (en) 2009-09-17

Family

ID=41181778

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008053420A Pending JP2009208115A (en) 2008-03-04 2008-03-04 Method and device for calculating parameter of rolling control, and rolling simulation device

Country Status (1)

Country Link
JP (1) JP2009208115A (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103384572A (en) * 2011-03-18 2013-11-06 株式会社日立制作所 Rolling control device, rolling control method and rolling control program
CN103722022A (en) * 2013-12-29 2014-04-16 北京首钢自动化信息技术有限公司 Friction coefficient model optimizing system and method in rolling process
CN104307886A (en) * 2014-11-12 2015-01-28 中冶南方(武汉)自动化有限公司 Composite rolling mill clad-rolling large-lag AGC (automatic gain control) method
JP2015147249A (en) * 2014-01-09 2015-08-20 Jfeスチール株式会社 Rolling mill control method, rolling mill control device, and rolled material manufacturing method
CN104889175A (en) * 2015-05-27 2015-09-09 首钢京唐钢铁联合有限责任公司 Tension setting method for improving leveling stability and product surface quality
JP2016074025A (en) * 2014-10-09 2016-05-12 Jfeスチール株式会社 Control method and control device of rolling machine
TWI554340B (en) * 2014-09-10 2016-10-21 東芝三菱電機產業系統股份有限公司 Rolling simulation device
CN111159952A (en) * 2019-12-31 2020-05-15 唐山市德龙钢铁有限公司 A method for alleviating rough rolling wrinkles of rod and wire rod
CN113649420A (en) * 2021-07-28 2021-11-16 首钢京唐钢铁联合有限责任公司 Temper mill rolling force obtaining method and device
CN116140374A (en) * 2023-04-14 2023-05-23 太原科技大学 A comprehensive quality prediction and process control method for strip rolling process
CN116727459A (en) * 2023-02-17 2023-09-12 山西太钢不锈钢股份有限公司 Method for setting rolling low-thickness high-grade non-oriented silicon steel model
CN117161111A (en) * 2023-08-02 2023-12-05 武汉科技大学 Prediction method and system for cross stiffness of rolling mill
CN118268386A (en) * 2024-06-03 2024-07-02 太原理工大学 Rolling force determination method and device for dynamic rolling of rolled piece from thin to thick
CN119422115A (en) * 2023-04-07 2025-02-11 发那科株式会社 Processing simulation device and processing simulation method
CN119972826A (en) * 2025-01-16 2025-05-13 北京科技大学 A method for optimizing the process parameters of seamless steel pipe outer diameter

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103384572A (en) * 2011-03-18 2013-11-06 株式会社日立制作所 Rolling control device, rolling control method and rolling control program
CN103722022A (en) * 2013-12-29 2014-04-16 北京首钢自动化信息技术有限公司 Friction coefficient model optimizing system and method in rolling process
CN103722022B (en) * 2013-12-29 2015-07-01 北京首钢自动化信息技术有限公司 Friction coefficient model optimizing system and method in rolling process
JP2015147249A (en) * 2014-01-09 2015-08-20 Jfeスチール株式会社 Rolling mill control method, rolling mill control device, and rolled material manufacturing method
TWI554340B (en) * 2014-09-10 2016-10-21 東芝三菱電機產業系統股份有限公司 Rolling simulation device
JP2016074025A (en) * 2014-10-09 2016-05-12 Jfeスチール株式会社 Control method and control device of rolling machine
CN104307886A (en) * 2014-11-12 2015-01-28 中冶南方(武汉)自动化有限公司 Composite rolling mill clad-rolling large-lag AGC (automatic gain control) method
CN104889175B (en) * 2015-05-27 2017-04-26 首钢京唐钢铁联合有限责任公司 Tension setting method for improving leveling stability and product surface quality
CN104889175A (en) * 2015-05-27 2015-09-09 首钢京唐钢铁联合有限责任公司 Tension setting method for improving leveling stability and product surface quality
CN111159952A (en) * 2019-12-31 2020-05-15 唐山市德龙钢铁有限公司 A method for alleviating rough rolling wrinkles of rod and wire rod
CN113649420A (en) * 2021-07-28 2021-11-16 首钢京唐钢铁联合有限责任公司 Temper mill rolling force obtaining method and device
CN113649420B (en) * 2021-07-28 2023-10-24 首钢京唐钢铁联合有限责任公司 Temper mill rolling force obtaining method and device
CN116727459A (en) * 2023-02-17 2023-09-12 山西太钢不锈钢股份有限公司 Method for setting rolling low-thickness high-grade non-oriented silicon steel model
CN119422115A (en) * 2023-04-07 2025-02-11 发那科株式会社 Processing simulation device and processing simulation method
CN116140374A (en) * 2023-04-14 2023-05-23 太原科技大学 A comprehensive quality prediction and process control method for strip rolling process
CN117161111A (en) * 2023-08-02 2023-12-05 武汉科技大学 Prediction method and system for cross stiffness of rolling mill
CN118268386A (en) * 2024-06-03 2024-07-02 太原理工大学 Rolling force determination method and device for dynamic rolling of rolled piece from thin to thick
CN119972826A (en) * 2025-01-16 2025-05-13 北京科技大学 A method for optimizing the process parameters of seamless steel pipe outer diameter

Similar Documents

Publication Publication Date Title
JP2009208115A (en) Method and device for calculating parameter of rolling control, and rolling simulation device
CN106166566B (en) Hot finisher goes out side temperature control equipment and its control method
KR960007487B1 (en) Control device of continuous hot rolling mill
JP2018005544A (en) Plant controller, rolling controller, plant control method, and plant control program
JPH11707A (en) Rolling mill control method and rolling mill control device
JP5790636B2 (en) Rolled material meander control method, rolled material meander control device, rolled material meander control program, and rolled material manufacturing method
Prinz et al. Online parameter estimation for adaptive feedforward control of the strip thickness in a hot strip rolling mill
JP3902585B2 (en) Sheet shape control method in cold rolling
JP5380199B2 (en) Method of controlling tension in rolling device and rolling device
CN100515592C (en) Rolling mill control device and rolling mill control method
JP2007245204A (en) Method and apparatus for learning rolling load model
JP5272677B2 (en) Rolling mill tension control method and rolling mill tension control apparatus
JP2012121063A (en) Method and device for controlling tandem rolling mill
JP3501982B2 (en) Sheet width control device
CN115502217B (en) snake control device
KR100531145B1 (en) Sheet width control method in hot rolling
JPH0232041B2 (en)
JP2005111546A (en) Thickness / crown control method
Ozaki et al. Nonlinear receding horizon control of thickness and tension in a tandem cold mill with a variable rolling speed
Pittner et al. A simple rolling mill model with linear quadratic optimal controller
KR100293222B1 (en) Rolling load prediction method using fuzzy theory
JP3389903B2 (en) Metal strip rolling control method
KR101105900B1 (en) Heavy plate rolling load prediction method
JP3205175B2 (en) Strip width control method in hot rolling
JP3908702B2 (en) Sheet width control method for continuous rolling mill