(57)【要約】【課題】鋳片全体の表面温度を予め定めた目標温度に制御する際の精度を高めることが可能な、連続鋳造機の２次冷却制御方法を提供する。【解決手段】鋳片表面温度測定工程と、鋳造速度把握工程と、トラッキング面設定工程と、鋳片目標温度設定工程と、温度固相率推定工程と、熱伝達係数推定工程と、伝熱凝固モデルパラメータ修正工程と、将来予測面設定工程と、将来予測工程と、将来温度影響係数予測工程と、鋳片表面参照温度算出工程と、最適化問題係数行列算出工程と、最適化問題求解工程と、冷却水量変更工程とを有し、該冷却水量変更工程で冷却ゾーン毎の冷却水量の変更を繰り返すことにより、鋳造中の任意の時刻において各トラッキング面が２次冷却制御対象の冷却ゾーン出口まで移動する間に、将来予測面の、将来予測面位置における鋳片の表面温度を、鋳片目標温度設定工程で定めた鋳片の表面温度の目標値に制御する、連続鋳造機の２次冷却制御方法とする。【選択図】図３
2014140862.tif 000002
【特許請求の範囲】
【請求項１】
連続鋳造機の鋳型から引き抜かれた鋳片を冷却する２次冷却帯を、前記鋳片の鋳造方向に複数の冷却ゾーンへと分割し、前記鋳片へ向けて噴射される冷却水量を各冷却ゾーンで制御することにより、前記鋳片の表面温度を制御する方法において、
予め定めたストランド内の温度測定点における前記鋳片の表面温度を、前記鋳片の鋳造中に測定する、鋳片表面温度測定工程と、
前記連続鋳造機の鋳造速度を把握する、鋳造速度把握工程と、
前記鋳片の断面内温度、前記鋳片の表面温度、および、前記鋳片の固相率分布を計算する対象であるトラッキング面を、鋳型内湯面位置から少なくとも２次冷却制御対象の冷却ゾーン出口までの領域で、予め定めた間隔で設定する、トラッキング面設定工程と、
前記トラッキング面における、前記鋳片の表面温度の目標値を定める、鋳片目標温度設定工程と、
鋳造が進むことにより、前記トラッキング面が前記鋳片の鋳造方向へ予め定めた間隔だけ進む毎に、伝熱方程式に基づく伝熱凝固モデルにより、前記鋳造方向に垂直な前記鋳片の断面内温度、前記鋳片の表面温度、および、前記鋳片の固相率分布を算出して更新する、温度固相率推定工程と、
前記伝熱凝固モデルで用いる前記鋳片の表面の熱伝達係数を、前記冷却水量を含む鋳造条件を用いて算出する、熱伝達係数推定工程と、
前記鋳片表面温度測定工程で測定された前記鋳片の表面温度と、前記温度固相率推定工程で推定された前記鋳片の表面温度との差を用いて、前記伝熱凝固モデルにおける鋳造条件に対するパラメータを修正する、伝熱凝固モデルパラメータ修正工程と、
前記トラッキング面設定工程で設定された前記トラッキング面の集合の中から、予め定めた鋳造方向に一定の間隔で、将来時刻における前記鋳片の表面温度、前記鋳造方向に垂直な前記鋳片の断面内温度、および、前記鋳片の固相率分布を予測する将来予測面を設定する、将来予測面設定工程と、
鋳造が進むことによって、任意の前記将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む間に、鋳造速度が現在時刻から変化しないと仮定して、それぞれの前記将来予測面が前記将来予測面位置に到達したときの前記鋳片の表面温度、前記鋳造方向に垂直な前記鋳片の断面内温度、および、前記鋳片の固相率分布を、前記将来予測面設定工程で用いた前記間隔毎に、前記伝熱凝固モデルを用いて繰り返し予測して更新する、将来予測工程と、
鋳造が進むことによって、任意の前記将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む毎に、鋳造速度が現在時刻から変化しないと仮定して、前記各冷却ゾーンの冷却水量がステップ関数状に変化した場合の、それぞれの前記将来予測面が前記将来予測面位置に到達するまでに通過する、各トラッキング面位置における前記鋳片の表面温度を予測し、該予測した前記鋳片の表面温度と、前記将来予測工程で予測した前記鋳片の表面温度との偏差を求め、該偏差を用いて、ステップ関数状に変化する前記冷却水量に対する変化影響係数を求める、将来温度影響係数予測工程と、
前記鋳片目標温度設定工程で設定した前記鋳片の表面温度の目標値と、前記将来温度影響係数予測工程で予測した、前記将来予測面が前記将来予測面位置に到達した時点における前記鋳片の表面温度の予測値との間の値である、時間に応じて決定される参照目標温度を算出する、鋳片表面参照温度算出工程と、
現在時刻における前記各冷却ゾーンの冷却水量を決定変数とし、前記将来予測工程および前記将来温度影響係数予測工程の各々においてそれぞれの前記将来予測面が通過した各将来予測面位置における将来温度影響係数、および、前記鋳片表面参照温度算出工程で算出した前記参照目標温度と前記将来予測工程で予測した前記鋳片の表面温度との偏差を算出し、それぞれの前記将来予測面で算出した該偏差の和を最小化する最適化問題の２次計画問題とし、該２次計画問題における決定変数に対する係数行列を算出する、最適化問題係数行列算出工程と、
前記２次計画問題を数値的に解くことにより、ステップ関数状に変化する前記冷却水量の変更量の、現在時刻における最適値を求める、最適化問題求解工程と、
前記最適値を、現在の冷却ゾーンの冷却水量へと加えることにより冷却水量を変更する、冷却水量変更工程と、を有し、
前記冷却水量変更工程で、前記冷却水量の変更を繰り返すことにより、鋳造中の任意の時刻において各トラッキング面が前記２次冷却制御対象の冷却ゾーン出口まで移動する間に、前記将来予測面の、前記将来予測面位置における前記鋳片の表面温度を、前記鋳片目標温度設定工程で定めた前記鋳片の表面温度の目標値に制御することを特徴とする、連続鋳造機の２次冷却制御方法。
【請求項２】
連続鋳造機の鋳型から引き抜かれた鋳片を冷却する２次冷却帯を、前記鋳片の鋳造方向に複数の冷却ゾーンへと分割し、前記鋳片へ向けて噴射される冷却水量を各冷却ゾーンで制御することにより、前記鋳片の表面温度を制御する装置であって、
予め定めたストランド内の温度測定点における前記鋳片の表面温度を、前記鋳片の鋳造中に測定する、鋳片表面温度測定部と、
前記連続鋳造機の鋳造速度を把握する、鋳造速度把握部と、
前記鋳片の断面内温度、前記鋳片の表面温度、および、前記鋳片の固相率分布を計算する対象であるトラッキング面を、鋳型内湯面位置から少なくとも２次冷却制御対象の冷却ゾーン出口までの領域で、予め定めた間隔で設定する、トラッキング面設定部と、
前記トラッキング面における、前記鋳片の表面温度の目標値を定める、鋳片目標温度設定部と、
鋳造が進むことにより、前記トラッキング面が前記鋳片の鋳造方向へ予め定めた間隔だけ進む毎に、伝熱方程式に基づく伝熱凝固モデルにより、前記鋳造方向に垂直な前記鋳片の断面内温度、前記鋳片の表面温度、および、前記鋳片の固相率分布を算出して更新する、温度固相率推定部と、
前記伝熱凝固モデルで用いる前記鋳片の表面の熱伝達係数を、前記冷却水量を含む鋳造条件を用いて算出する、熱伝達係数推定部と、
前記鋳片表面温度測定部で測定された前記鋳片の表面温度と、前記温度固相率推定部で推定された前記鋳片の表面温度との差を用いて、前記伝熱凝固モデルにおける鋳造条件に対するパラメータを修正する、伝熱凝固モデルパラメータ修正部と、
前記トラッキング面設定部で設定された前記トラッキング面の集合の中から、予め定めた鋳造方向に一定の間隔で、将来時刻における前記鋳片の表面温度、前記鋳造方向に垂直な前記鋳片の断面内温度、および、前記鋳片の固相率分布を予測する将来予測面を設定する、将来予測面設定部と、
鋳造が進むことによって、任意の前記将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む間に、鋳造速度が現在時刻から変化しないと仮定して、それぞれの前記将来予測面が前記将来予測面位置に到達したときの前記鋳片の表面温度、前記鋳造方向に垂直な前記鋳片の断面内温度、および、前記鋳片の固相率分布を、前記将来予測面設定部で用いた前記間隔毎に、前記伝熱凝固モデルを用いて繰り返し予測して更新する、将来予測部と、
鋳造が進むことによって、任意の前記将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む毎に、鋳造速度が現在時刻から変化しないと仮定して、前記各冷却ゾーンの冷却水量がステップ関数状に変化した場合の、それぞれの前記将来予測面が前記将来予測面位置に到達するまでに通過する、各トラッキング面位置における前記鋳片の表面温度を予測し、該予測した前記鋳片の表面温度と、前記将来予測部で予測した前記鋳片の表面温度との偏差を求め、該偏差を用いて、ステップ関数状に変化する前記冷却水量に対する変化影響係数を求める、将来温度影響係数予測部と、
前記鋳片目標温度設定部で設定した前記鋳片の表面温度の目標値と、前記将来温度影響係数予測部で予測した、前記将来予測面が前記将来予測面位置に到達した時点における前記鋳片の表面温度の予測値との間の値である、時間に応じて決定される参照目標温度を算出する、鋳片表面参照温度算出部と、
現在時刻における前記各冷却ゾーンの冷却水量を決定変数とし、前記将来予測部および前記将来温度影響係数予測部の各々においてそれぞれの前記将来予測面が通過した各将来予測面位置における将来温度影響係数、および、前記鋳片表面参照温度算出部で算出した前記参照目標温度と前記将来予測部で予測した前記鋳片の表面温度との偏差を算出し、それぞれの前記将来予測面で算出した該偏差の和を最小化する最適化問題の２次計画問題とし、該２次計画問題における決定変数に対する係数行列を算出する、最適化問題係数行列算出部と、
前記２次計画問題を数値的に解くことにより、ステップ関数状に変化する前記冷却水量の変更量の、現在時刻における最適値を求める、最適化問題求解部と、
前記最適値を、現在の冷却ゾーンの冷却水量へと加えることにより冷却水量を変更する、冷却水量変更部と、を有し、
前記冷却水量変更部で、前記冷却水量の変更を繰り返すことにより、鋳造中の任意の時刻において各トラッキング面が前記２次冷却制御対象の冷却ゾーン出口まで移動する間に、前記将来予測面の、前記将来予測面位置における前記鋳片の表面温度を、前記鋳片目標温度設定部で定めた前記鋳片の表面温度の目標値に制御することを特徴とする、連続鋳造機の２次冷却制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、連続鋳造機の２次冷却帯において、一部または全部の鋳片鋳造方向若しくは幅方向の表面温度分布を制御する、連続鋳造機の２次冷却制御方法および２次冷却制御装置に関する。
【背景技術】
【０００２】
鋼の連続鋳造においては、例えば垂直曲げ型連続鋳造機では、垂直な鋳型から引き出した鋳片を一旦湾曲させた後、一定湾曲半径で引抜き、その後矯正部において曲がりをなくした状態の鋳片として抽出し、切断する。ところが、ストランド（「鋳型＋２次冷却帯群＋ローラー群を有する引抜き装置」のセットを意味する。以下において同じ。）の曲げ部においては鋳片の下側表面に、矯正部においては鋳片上側表面に引っ張り応力がかかるため、鋳片表面の温度が脆化域とよばれる範囲にある場合、横ひび割れと呼ばれる表面割れ疵が発生することがある。このため、ストランドの曲げ部および矯正部において、鋳片表面部温度が上記脆化域を回避するように、冷却水量分布を適切に設定することが必要である。冷却水量分布の適切な設定は、例えば、一定鋳造速度の場合には冷却ゾーン水量分布を、事前にシミュレーションなどで適正な値に定めておくことで達成できる。
【０００３】
ところが、連々鋳における次の取鍋到着が遅れる場合には、連々鋳が中断しないように鋳造速度を所定値より低下させて到着を待つため、鋳造速度を操業中に変化させる必要がある。このとき、変更中の鋳造速度について、事前に鋳造速度に対して設定した各ゾーン水量を補間して設定する従来のカスケード水量制御では、鋳片の鋳型湯面から切断までの時間に対する冷却履歴が乱れ、表面の横ひび割れなどの鋳片品質不良が発生する。
【０００４】
また、鋳片表面のスケール付着などの影響により、冷却水量と表面の熱伝達係数との関係が事前のシミュレーションで仮定したものから変化する場合がある。このような場合にも鋳片表面温度が脆化域に入り、横ひび割れが発生することがある。
【０００５】
このような問題に対し、これまでに、いわゆるモデル予測制御による制御方法が開示されている。例えば、特許文献１には、引抜き鋳片を定間隔毎にトラッキングし、各トラッキング面の温度分布を伝熱モデルに基づき逐次的に計算し、鋳片引抜き軌跡をいくつかのゾーンに分割した各ゾーンの出側における計算温度と実測温度との関係から学習された熱伝達係数によって上記モデルを修正し、上記軌跡に沿って設けられた測温点における各トラッキング面の温度分布を上記修正モデルに基づいて一定時刻毎に予測するとともに、当該位置における目標温度と予測温度の差から求めたフィードフォワード水量と、実測温度と目標温度の差から求めたフィードバック水量と、を合計した水量を鋳片に散布する表面温度制御方法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開昭５７−１５４３６４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
特許文献１に開示されているフィードフォワード水量の算出方法では、冷却ゾーンに存在するトラッキング点毎に、各々が冷却ゾーン出口の測温点に到達した時点での温度を予測し、各トラッキング点が測温点に到達した時の温度予測値が目標値に一致する予測水量密度を求め、さらに当該冷却ゾーンの全トラッキング面について予測水量密度の重みつき平均値をフィードフォワード水量とする。この技術では、鋳型側の冷却ゾーンから順に、フィードフォワード水量を求める手続きと、この手続きで求めたフィードフォワード水量を用いた当該冷却ゾーンにおける温度分布の再計算を行って再計算温度を求める手続きとを行い、再計算温度を、下流側に隣接している冷却ゾーンの入口における初期温度とする手続きを繰返し、全冷却ゾーンの冷却水量を決定する。ところがこの技術では、再計算温度を下流側に隣接している冷却ゾーン入口における初期温度としても、下流側に隣接している冷却ゾーンの入口以外のトラッキング点の温度計算（再計算温度を求めた冷却ゾーンの下流側に隣接している冷却ゾーンよりもさらに下流側に存在している冷却ゾーンにおけるトラッキング点の温度計算）には、フィードフォワード水量の影響が表れない。したがって、特許文献１に開示されている技術では、温度予測計算において、上流側の水量変化が正しく反映されるまでの所要時間が長くなり、場合によっては水量がハンチングするなどの問題が発生していた。その結果、鋳片全体の表面温度を予め定めた目標温度に制御する際の精度が低下しやすかった。
【０００８】
そこで、本発明は、鋳片全体の表面温度を予め定めた目標温度に制御する際の精度を高めることが可能な、連続鋳造機の２次冷却制御方法および２次冷却制御装置を提供することを課題とする。
【課題を解決するための手段】
【０００９】
本発明の第１の態様は、連続鋳造機の鋳型から引き抜かれた鋳片を冷却する２次冷却帯を、鋳片の鋳造方向に複数の冷却ゾーンへと分割し、鋳片へ向けて噴射される冷却水量を各冷却ゾーンで制御することにより、鋳片の表面温度を制御する方法において、予め定めたストランド内の温度測定点における鋳片の表面温度を、鋳片の鋳造中に測定する鋳片表面温度測定工程と、連続鋳造機の鋳造速度を把握する鋳造速度把握工程と、鋳片の断面内温度、鋳片の表面温度、および、鋳片の固相率分布を計算する対象であるトラッキング面を、鋳型内湯面位置から少なくとも２次冷却制御対象の冷却ゾーン出口までの領域で、予め定めた間隔で設定するトラッキング面設定工程と、トラッキング面における、鋳片の表面温度の目標値を定める鋳片目標温度設定工程と、鋳造が進むことにより、トラッキング面が鋳片の鋳造方向へ予め定めた間隔だけ進む毎に、伝熱方程式に基づく伝熱凝固モデルにより、鋳造方向に垂直な鋳片の断面内温度、鋳片の表面温度、および、鋳片の固相率分布を算出して更新する温度固相率推定工程と、伝熱凝固モデルで用いる鋳片の表面の熱伝達係数を、上記冷却水量を含む鋳造条件を用いて算出する熱伝達係数推定工程と、鋳片表面温度測定工程で測定された鋳片の表面温度と、温度固相率推定工程で推定された鋳片の表面温度との差を用いて、伝熱凝固モデルにおける鋳造条件に対するパラメータを修正する伝熱凝固モデルパラメータ修正工程と、トラッキング面設定工程で設定されたトラッキング面の集合の中から、予め定めた鋳造方向に一定の間隔で、将来時刻における鋳片の表面温度、鋳造方向に垂直な鋳片の断面内温度、および、鋳片の固相率分布を予測する将来予測面を設定する将来予測面設定工程と、鋳造が進むことによって、任意の将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む間に、鋳造速度が現在時刻から変化しないと仮定して、それぞれの将来予測面が将来予測面位置に到達したときの鋳片の表面温度、鋳造方向に垂直な鋳片の断面内温度、および、鋳片の固相率分布を、将来予測面設定工程で用いた間隔毎に、伝熱凝固モデルを用いて繰り返し予測して更新する将来予測工程と、鋳造が進むことによって、任意の将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む毎に、鋳造速度が現在時刻から変化しないと仮定して、各冷却ゾーンの冷却水量がステップ関数状に変化した場合の、それぞれの将来予測面が将来予測面位置に到達するまでに通過する、各トラッキング面位置における鋳片の表面温度を予測し、該予測した鋳片の表面温度と、将来予測工程で予測した鋳片の表面温度との偏差を求め、該偏差を用いて、ステップ関数状に変化する冷却水量に対する変化影響係数を求める将来温度影響係数予測工程と、鋳片目標温度設定工程で設定した鋳片の表面温度の目標値と、将来温度影響係数予測工程で予測した、将来予測面が将来予測面位置に到達した時点における鋳片の表面温度の予測値との間の値である、時間に応じて決定される参照目標温度を算出する鋳片表面参照温度算出工程と、現在時刻における各冷却ゾーンの冷却水量を決定変数とし、将来予測工程および将来温度影響係数予測工程の各々においてそれぞれの将来予測面が通過した各将来予測面位置における将来温度影響係数、および、鋳片表面参照温度算出工程で算出した参照目標温度と将来予測工程で予測した鋳片の表面温度との偏差を算出し、それぞれの将来予測面で算出した該偏差の和を最小化する最適化問題の２次計画問題とし、該２次計画問題における決定変数に対する係数行列を算出する最適化問題係数行列算出工程と、上記２次計画問題を数値的に解くことにより、ステップ関数状に変化する冷却水量の変更量の、現在時刻における最適値を求める最適化問題求解工程と、該最適値を、現在の冷却ゾーンの冷却水量へと加えることにより冷却水量を変更する冷却水量変更工程と、を有し、該冷却水量変更工程で冷却水量の変更を繰り返すことにより、鋳造中の任意の時刻において各トラッキング面が２次冷却制御対象の冷却ゾーン出口まで移動する間に、将来予測面の、将来予測面位置における鋳片の表面温度を、鋳片目標温度設定工程で定めた鋳片の表面温度の目標値に制御することを特徴とする、連続鋳造機の２次冷却制御方法である。
【００１０】
本発明の第２の態様は、連続鋳造機の鋳型から引き抜かれた鋳片を冷却する２次冷却帯を、鋳片の鋳造方向に複数の冷却ゾーンへと分割し、鋳片へ向けて噴射される冷却水量を各冷却ゾーンで制御することにより、鋳片の表面温度を制御する装置であって、予め定めたストランド内の温度測定点における鋳片の表面温度を、鋳片の鋳造中に測定する鋳片表面温度測定部と、連続鋳造機の鋳造速度を把握する鋳造速度把握部と、鋳片の断面内温度、鋳片の表面温度、および、鋳片の固相率分布を計算する対象であるトラッキング面を、鋳型内湯面位置から少なくとも２次冷却制御対象の冷却ゾーン出口までの領域で、予め定めた間隔で設定するトラッキング面設定部と、トラッキング面における、鋳片の表面温度の目標値を定める鋳片目標温度設定部と、鋳造が進むことにより、トラッキング面が鋳片の鋳造方向へ予め定めた間隔だけ進む毎に、伝熱方程式に基づく伝熱凝固モデルにより、鋳造方向に垂直な鋳片の断面内温度、鋳片の表面温度、および、鋳片の固相率分布を算出して更新する温度固相率推定部と、伝熱凝固モデルで用いる鋳片の表面の熱伝達係数を、冷却水量を含む鋳造条件を用いて算出する熱伝達係数推定部と、鋳片表面温度測定部で測定された鋳片の表面温度と、温度固相率推定部で推定された鋳片の表面温度との差を用いて、伝熱凝固モデルにおける鋳造条件に対するパラメータを修正する伝熱凝固モデルパラメータ修正部と、トラッキング面設定部で設定されたトラッキング面の集合の中から、予め定めた鋳造方向に一定の間隔で、将来時刻における鋳片の表面温度、鋳造方向に垂直な鋳片の断面内温度、および、鋳片の固相率分布を予測する将来予測面を設定する将来予測面設定部と、鋳造が進むことによって、任意の将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む間に、鋳造速度が現在時刻から変化しないと仮定して、それぞれの将来予測面が将来予測面位置に到達したときの鋳片の表面温度、鋳造方向に垂直な鋳片の断面内温度、および、鋳片の固相率分布を、将来予測面設定部で用いた間隔毎に、伝熱凝固モデルを用いて繰り返し予測して更新する将来予測部と、鋳造が進むことによって、任意の将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む毎に、鋳造速度が現在時刻から変化しないと仮定して、各冷却ゾーンの冷却水量がステップ関数状に変化した場合の、それぞれの将来予測面が将来予測面位置に到達するまでに通過する、各トラッキング面位置における鋳片の表面温度を予測し、該予測した鋳片の表面温度と、将来予測部で予測した鋳片の表面温度との偏差を求め、該偏差を用いて、ステップ関数状に変化する冷却水量に対する変化影響係数を求める将来温度影響係数予測部と、鋳片目標温度設定部で設定した鋳片の表面温度の目標値と、将来温度影響係数予測部で予測した、将来予測面が将来予測面位置に到達した時点における鋳片の表面温度の予測値との間の値である、時間に応じて決定される参照目標温度を算出する鋳片表面参照温度算出部と、現在時刻における各冷却ゾーンの冷却水量を決定変数とし、将来予測部および将来温度影響係数予測部の各々においてそれぞれの将来予測面が通過した各将来予測面位置における将来温度影響係数、および、鋳片表面参照温度算出部で算出した参照目標温度と将来予測部で予測した鋳片の表面温度との偏差を算出し、それぞれの将来予測面で算出した該偏差の和を最小化する最適化問題の２次計画問題とし、該２次計画問題における決定変数に対する係数行列を算出する最適化問題係数行列算出部と、上記２次計画問題を数値的に解くことにより、ステップ関数状に変化する冷却水量の変更量の、現在時刻における最適値を求める最適化問題求解部と、該最適値を、現在の冷却ゾーンの冷却水量へと加えることにより冷却水量を変更する冷却水量変更部と、を有し、該冷却水量変更部で、冷却水量の変更を繰り返すことにより、鋳造中の任意の時刻において各トラッキング面が２次冷却制御対象の冷却ゾーン出口まで移動する間に、将来予測面の、将来予測面位置における鋳片の表面温度を、鋳片目標温度設定部で定めた鋳片の表面温度の目標値に制御することを特徴とする、連続鋳造機の２次冷却制御装置である。
【発明の効果】
【００１１】
本発明によれば、鋳片全体の表面温度を、予め定めた目標温度に常に一致するように制御することが可能な、連続鋳造機の２次冷却制御方法および２次冷却制御装置を提供することができる。その結果、いかなる鋳造速度でも、また鋳造速度が鋳造中に変化した場合でも、連続鋳造機の曲げセグメントや矯正セグメントにおいて、表面温度を鋼の脆化域を回避するように制御することが可能になる。したがって、本発明によれば、表面疵による欠陥のない鋳片を製造することが可能になる。
【図面の簡単な説明】
【００１２】
【図１】連続鋳造機９および冷却制御装置１０を説明する図である。
【図２】鋳造方向に垂直な鋳片断面の分割と格子点の例を示す図である。
【図３】本発明の冷却制御方法を説明する図である。
【図４】各将来予測面がその下流側に隣接する将来予測面位置まで移動する間に、表面温度を評価するトラッキング面の位置と、温度を予測する相対時刻との関係を説明する図である。
【図５】冷却制御装置１０に備えられている各部の関係とやりとりする情報を説明するブロック線図である。
【図６】鋳造速度低下時に、本発明の冷却制御方法を適用した場合の結果を表すグラフである。
【図７】鋳造速度低下時に、従来のカスケード水量制御を適用した場合の結果を表すグラフである。
【図８】鋳造中に第３冷却ゾーンの出口目標温度を変更した場合に、本発明の冷却制御方法で冷却水量を調節して表面温度を制御した結果を表すグラフである。
【図９】第４冷却ゾーンのスプレー熱伝達係数が低下した場合に、本発明の冷却制御方法で冷却水量を調節して鋳片表面温度を制御した結果を表すグラフである。
【発明を実施するための形態】
【００１３】
以下、本発明の実施の形態について説明する。なお、以下に説明する形態は本発明の例示であり、本発明は以下に説明する形態に限定されない。
【００１４】
図１は、本発明を実施する連続鋳造機９、および、本発明にかかる連続鋳造機の２次冷却制御装置（以下において、「冷却制御装置」ということがある。）１０を説明する図である。図１では、連続鋳造機９および冷却制御装置１０を簡略化して示している。
本発明を実施する連続鋳造機９では、外側が凝固したストランドをロール対で挟んで支持しながら、駆動装置を備えたピンチロールによって、鋳型１からストランドが所定の引抜き速度（鋳造速度）で引抜かれる。符号４は溶鋼メニスカスである。鋳造方向に所定の間隔をあけて配置された隣接する支持ロールの間には、鋳片５へ向けて冷却水を散布するミストスプレー２（またはスプレー２）の噴出口が設置される。散布される冷却水の流量は冷却水配管に設置した流量調整弁３により制御される。流量調整弁３の開度は、冷却制御装置１０から与えられる水量指示値に基づいて調節される。冷却水配管は、鋳片５の鋳造方向長さを複数個に区分した冷却ゾーン（冷却ゾーン境界線６によって区分された冷却ゾーン）に対応して設置されるので、ストランド内の鋳造方向冷却水量分布は、冷却ゾーンごとに制御される。以下の説明において、鋳型直下の冷却ゾーンから順に、第１冷却ゾーン、第２冷却ゾーン、…ということがある。
【００１５】
ストランド内における鋳片５の温度および固相率の分布は、鋳型内湯面から最終ロール出側まで鋳造方向に一定間隔で設置した計算点で、鋳片５に垂直な断面を設定し、各断面内の温度および固相率分布を、各計算点における冷却条件を反映した熱伝達係数の境界条件のもとで離散化した熱伝導方程式を解くことで計算する。熱伝導方程式の初期条件には、計算対象位置に存在する断面の上流側に隣接する断面の温度および固相率の計算結果を設定し、当該上流側に隣接する計算点から対象計算位置へ、鋳片引抜きにより断面が移動するまでの計算を繰り返すことにより、鋳片全体の温度および固相率を計算することができる。
【００１６】
熱伝導方程式の離散化には、例えば図２に示した直交する格子の二次元モデルを用いる。各格子点（ｉ、ｊ）における温度Ｔｉｊ、単位質量あたりのエンタルピーＨｉｊ、および、単位質量あたりの固相率ｆｉｊを変数とし、各格子点（ｉ、ｊ）における物性定数を、温度依存性を考慮して密度ρｉｊ、比熱Ｃｉｊ、および、熱伝導率λｉｊとして表す。このとき、エンタルピーＨｉｊ、温度Ｔｉｊ、および、固相率ｆｉｊの関係は、式（１）で表される。
【００１７】
【数１】
2014140862.tif 000003
【００１８】
時間刻みΔｔの間に、鋳造方向位置ｚからｚ＋Δｚまで引き抜かれる断面のエンタルピーＨｉｊおよび固相率ｆｉｊの分布の時間変化は、離散化した熱伝導方程式（２）、（４）、（７）、初期条件式（３）、および、境界条件式（５）、（６）、（８）、（９）を用いて表される。以下の式において、上付き添え字ｚは鋳造方向位置を表し、鋳型内湯面位置をｚ＝０とする。熱伝導方程式における時間刻みΔｔは、鋳造方向の断面設置刻みΔｚと時刻ｔ−１における鋳造速度ｖ（ｔ−１）を用いて、Δｔ＝Δｚ／ｖ（ｔ−１）に変換する。鋳片表面からの抜熱は、鋳片５へ向けて散布された冷却水による冷却、ロールとの接触、および、放射等鋳造方向断面位置による冷却方法の違いを考慮した境界条件を反映して設定し、ここでは、式（５）および式（８）に示した、外部を代表する温度ＴＥと表面温度Ｔｉｊｚとの差の１次式で表した時の熱伝達係数ＫｘまたはＫｙで代表した。
【００１９】
【数２】
2014140862.tif 000004
【００２０】
【数３】
2014140862.tif 000005
【００２１】
上記式（２）において、ｑｉ＋１／２、ｊｚは鋳造方向位置ｚ−１における幅方向の格子点（ｉ、ｊ）から格子点（ｉ＋１、ｊ）への熱流束であり、幅方向内部をｉ＝２、…、Ｉとするとき、下記式（４）で表される。
【００２２】
【数４】
2014140862.tif 000006
【００２３】
なお、上記式（１）におけるＬｉｊは、格子点（ｉ、ｊ）における凝固潜熱λｉ＋１／２、ｊ＝（λｉ＋１、ｊ＋λｉｊ）／２である。上記式（２）におけるΔｘｉは格子点（ｉ−１／２、ｊ）から格子点（ｉ＋１／２、ｊ）までの距離であり、上記式（２）におけるΔｙｉは格子点（ｉ、ｊ−１／２）から格子点（ｉ、ｊ＋１／２）までの距離である。また、幅方向境界条件は、短辺表面をｉ＝１とするとき、鋳造方向位置ｚ−１における熱伝達係数Ｋｘおよび外部代表温度ＴＥを用いて下記式（５）で表される。
【００２４】
【数５】
2014140862.tif 000007
【００２５】
また、幅方向中央線上では、幅方向中央線上をｉ＝Ｉ＋１とするとき、下記式（６）で表される対称境界条件を仮定する。
【００２６】
【数６】
2014140862.tif 000008
【００２７】
また、上記式（２）において、ｑｚｉ、ｊ＋１／２は厚み方向の格子点（ｉ、ｊ）から格子点（ｉ、ｊ＋１）への熱流束であり、厚み方向内部をｊ＝２、…、Ｊとするとき、下記式（７）で表される。
【００２８】
【数７】
2014140862.tif 000009
【００２９】
なお、λｉ、ｊ＋１／２＝（λｉ、ｊ＋１＋λｉｊ）／２である。上記式（７）において、Δｙは格子点（ｉ、ｊ）から格子点（ｉ、ｊ＋１）までの距離である。また、厚み方向境界条件は、長辺表面をｊ＝１とするとき、鋳造方向位置ｚ−１における熱伝達係数Ｋｙおよび外部代表温度ＴＥを用いて下記式（８）で表される。
【００３０】
【数８】
2014140862.tif 000010
【００３１】
また、厚さ方向中央線上では、厚み中央線上をｊ＝Ｊ＋１とするとき、下記式（９）で表される対称境界条件を仮定する。
【００３２】
【数９】
2014140862.tif 000011
【００３３】
鋳造方向位置ｚ＋ΔｚにおけるエンタルピーＨｉｊｚ＋Δｚを算出した後、完全液相のｆｉｊｚ＋Δｚ＝０または完全固相のｆｉｊｚ＋Δｚ＝１の場合には、上記式（１）に各々の値を代入することにより、温度Ｔｉｊｚ＋Δｚを求める。一方、０＜ｆｉｊｚ＋Δｚ＜１の場合には、液相中の溶質濃度で定まる状態図で表される液相線温度ＴＬ（Ｃｋ）（Ｃｋは溶質成分ｋの濃度）に一致するが、Ｓｃｈｅｉｌの式などで知られるように液相中の溶質濃度は固相率に依存するので、下記式（１０）で表されるモデルを使用し、当該式（１０）と上記式（１）とを連立した方程式の解として、ｆｉｊｚ＋ΔｚおよびＴｉｊｚ＋Δｚを求める。
【００３４】
【数１０】
2014140862.tif 000012
【００３５】
ミストスプレー２から散布された冷却水が衝突する鋳片の表面から流出する熱流束が下記式（１１）で表されるとき、熱伝達係数ｋは下記式（１２）で求める。
【００３６】
【数１１】
2014140862.tif 000013
【００３７】
【数１２】
2014140862.tif 000014
ここで、ＴＳは表面温度［℃］、Ｄｗは表面水量密度［ｌ／ｍ２］、νａはミストスプレー空気流速［ｍ／ｓ］であり、α、β、γ、及び、ｃは各々定数である。
【００３８】
冷却制御装置１０は、鋳片５の引抜き速度と、タンディッシュ内における溶鋼温度と、冷却水温とを用いて、温度評価点における鋳片表面温度の予測値を求め、この予測値と各冷却ゾーン内において予め定めた温度評価点における鋳片表面温度の目標値との偏差と、冷却水量とにより定められる評価関数を最小化するように、各冷却ゾーンの冷却水量の最適値を算出する。本発明にかかる連続鋳造機の２次冷却制御方法（以下において、「本発明の冷却制御方法」ということがある。）では、一回の制御周期内で行う、以下に説明する計算を繰り返すことにより、各トラッキング面における鋳片表面温度を、予め定めた鋳片表面温度の目標値に制御する。本発明の冷却制御方法を説明する図３を参照しつつ、本発明の冷却制御方法について、以下に説明する。
【００３９】
図３に示したように、本発明の冷却制御方法は、鋳片表面温度測定工程（Ｓ１）と、鋳造速度把握工程（Ｓ２）と、トラッキング面設定工程（Ｓ３）と、鋳片目標温度設定工程（Ｓ４）と、温度固相率推定工程（Ｓ５）と、熱伝達係数推定工程（Ｓ６）と、伝熱凝固モデルパラメータ修正工程（Ｓ７）と、将来予測面設定工程（Ｓ８）と、将来予測工程（Ｓ９）と、将来温度影響係数予測工程（Ｓ１０）と、鋳片表面参照温度算出工程（Ｓ１１）と、最適化問題係数行列算出工程（Ｓ１２）と、最適化問題求解工程（Ｓ１３）と、冷却水量変更工程（Ｓ１４）と、を有している。
【００４０】
鋳片表面温度測定工程（以下において、「Ｓ１」ということがある。）は、予め定めたストランド内の鋳片表面上の温度測定点における鋳片表面温度を、鋳造中に、鋳片表面温度計７を用いて測定する工程である。
【００４１】
鋳造速度把握工程（以下において、「Ｓ２」ということがある。）は、鋳造速度測定ロール８を用いて、連続鋳造機９の鋳片引抜速度（鋳造速度）を逐次測定することにより、鋳造速度を把握する工程である。このほか、Ｓ２は、例えば、冷却制御装置１０の上位計算機（不図示）から、鋳造速度の設定値に関するデータを受信することにより、鋳造速度を把握する工程、とすることもできる。
【００４２】
トラッキング面設定工程（以下において、「Ｓ３」ということがある。）は、鋳片断面内温度、鋳片表面温度、および、固相率分布を計算する対象であるトラッキング面を、鋳型内湯面位置から少なくとも２次冷却制御対象の冷却ゾーン出口までの領域で、予め定めた間隔で設定する工程である。
【００４３】
鋳片目標温度設定工程（以下において、「Ｓ４」ということがある。）は、Ｓ３で設定したトラッキング面における鋳片表面温度の目標値を定める工程である。
【００４４】
温度固相率推定工程（以下において、「Ｓ５」ということがある。）は、鋳造が進むことにより、Ｓ３で定めたトラッキング面が鋳片の鋳造方向へ予め定めた間隔だけ進む毎に、伝熱方程式に基づく伝熱凝固モデルにより、鋳造方向に垂直な鋳片断面内温度、鋳片表面温度、および、固相率分布を算出して更新する工程である。
Ｓ５では、鋳片の鋳造方向に一定間隔で設定した垂直な断面における温度および固相率分布の、前回制御周期からの変更量を、鋼が凝固する際の変態発熱を考慮した熱伝導方程式を解くことにより算出する。
より具体的には、現在時刻をｔとし、上記式（２）乃至式（１０）を時刻ｔ−１と時刻ｔとの間の変数間の関係式とみなして、鋳型内湯面に隣接する計算点から２次冷却制御対象の冷却ゾーン出口までの各計算点における断面の温度および固相率分布を計算する。
【００４５】
熱伝達係数推定工程（以下において、「Ｓ６」ということがある。）は、伝熱凝固モデルで用いる鋳片表面の熱伝達係数（上記式（５）及び式（８）で表される熱伝達係数）を、現在時刻ｔにおける伝熱凝固モデルパラメータの推定値と、時刻ｔ−１における冷却水量等の鋳造条件を用いて算出する工程である。
【００４６】
伝熱凝固モデルパラメータ修正工程（以下において、「Ｓ７」ということがある。）は、Ｓ１で測定された鋳片の表面温度と、Ｓ５で推定された鋳片表面温度との差を用いて、伝熱凝固モデルにおける鋳造条件に対するパラメータを修正する工程である。
【００４７】
伝熱凝固モデルのパラメータの修正は、Ｓ１で測定された鋳片の表面温度とＳ５で推定された鋳片表面温度の推定値との誤差に補正係数をかけた値を、モデルパラメータ修正量として伝熱凝固モデルパラメータに加えることによって行う。鋳片の表面温度の測定点（以下において、「測温点」または「測温位置」ということがある。）が複数ある場合、補正係数は行列またはベクトルで表される。伝熱凝固モデルのパラメータの修正に用いる補正係数は、推定対象のパラメータ毎に以下の手順で求める。
【００４８】
１）補正対象のパラメータについて、現在の値から微小に変更した値を設定する。
２）予め定めた時間Ｔａを現在からさかのぼり、現在時刻ｔにおいて測温位置ｚｋにある断面が時刻ｔ−Ｔａにあった位置ｚｋ（ｔ−Ｔａ）における温度および固相率の断面内分布を初期値とする。そして、時刻ｔ−Ｔａにおける位置ｚｋ（ｔ−Ｔａ）からから現在時刻ｔにおける測温位置ｚｋまでの冷却条件の履歴を与えて、上記式（２）乃至（１０）の計算を繰返すことにより、現在時刻ｔにおいてパラメータを微小変更した場合の、測温点における温度推定値を算出する。上記遡及時間範囲Ｔａは、補正対象パラメータが測温位置ｚｋにある断面の状態に影響を及ぼす範囲に限定すればよい。
３）各パラメータ修正量に対する温度変化量の関係を表す線型関係式を、下記手順で求める。
パラメータθｌをΔθｌだけ変更したときに、Ｓ５で推定した表面温度Ｔｋ（ｔ）に対し、上記２）で算出した表面温度推定値がＴｋ＋ΔＴｋｌに変化したとすると、ΔＴｋｌは下記式（１３）で表すことができる。
【００４９】
【数１３】
2014140862.tif 000015
式（１３）におけるＡａｋｌの推定値は下記式（１４）で表される。
【００５０】
【数１４】
2014140862.tif 000016
なお、Ａａｋｌをｋ行ｌ列の成分とする行列をＡａと書くと、全修正対象パラメータによる測温点における表面温度への影響を合わせた温度変化推定値は、Δθｌを第ｌ成分とするベクトルΔθ＝［Δθ１  Δθ２  …  Δθｌ］Ｔを用いてＡａΔθと表される。
【００５１】
パラメータの最適修正量は、下記式（１５）で表される、各測温点の温度測定値Ｔａｋ（ｔ）とＴｋ（ｔ）との偏差ψａｋ（ｔ）を並べたベクトルψａ（ｔ）を、修正後パラメータによる温度変化ＡａΔθが、数値的計算誤差やデータのばらつきを考慮して最もよく近似するように決定する。
【００５２】
【数１５】
2014140862.tif 000017
【００５３】
すなわち、ΔＡａをゲイン行列Ａａの各成分の誤差を表す行列とするとき、
【００５４】
【数１６】
2014140862.tif 000018
を最小化する値を求める。ただし、<ｘ>は変数ｘの期待値を表す。
【００５５】
Ｊの最小値は解析的に解くことができ、Ｊを最小化するパラメータ修正量Δθ（ｔ）は下記式（１７）で表される。
【００５６】
【数１７】
2014140862.tif 000019
ただし、<ΔＡａ>＝０とする。ゲイン行列からなる<ΔＡａＴΔＡａ>は、ゲイン行列の各成分の相関が０であると仮定すれば、対角成分ΔＡａｉｉの分散を各々同じ位置の対角成分とする行列で表されるので、プロセスなどの知識により予め定めておく。
以上のようにして求めたパラメータ修正量Δθ（ｔ）を現在のパラメータに加えた
【００５７】
【数１８】
2014140862.tif 000020
を、次回時刻以降の制御操作量算出に用いる。
【００５８】
将来予測面設定工程（以下において、「Ｓ８」ということがある。）は、Ｓ３で設定したトラッキング面の集合の中から、予め定めた鋳造方向に一定の間隔で、将来時刻における鋳片表面温度、鋳片断面内温度、および、固相率分布を予測する将来予測面を設定する工程である。
【００５９】
将来予測工程（以下において、「Ｓ９」ということがある。）は、鋳造が進むことによって、Ｓ８で設定した任意の将来予測面が現在時刻から下流側に隣接する将来予測面位置まで進む間に、鋳造速度が現在時刻から変化しないと仮定して、Ｓ８で設定した各将来予測面が上記下流側に隣接する将来予測面位置に到達したときの鋳片表面温度、鋳片断面内温度、および、固相率分布を、Ｓ８で定めた間隔（伝熱計算間隔）毎に上記伝熱凝固モデルを用いて繰り返し予測して更新する工程である。Ｓ９では、現在時刻における鋳造速度、各冷却ゾーンの冷却水量、および、Ｓ７で修正した伝熱凝固モデルのパラメータの値を用いて、鋳片表面温度、鋳片断面内温度、および、固相率分布を予測する。予測計算の初期値には、Ｓ５で求めた現在時刻ｔにおける各将来温度予測面の鋳片表面温度、鋳片断面内温度、および、固相率分布の値を用いる。
【００６０】
図４は、Ｓ８で設定した各将来予測面がその下流側に隣接する将来予測面位置まで移動する間に、表面温度を評価するトラッキング面の位置と、温度を予測する相対時刻との関係を説明する図である。図４では、「●」で示した時刻に表面温度が予測されることを示している。図４に示した、複数の「●」を結んだ斜めの直線の傾きは、現在時刻ｔにおける鋳造速度ｖ（ｔ）に相当する。Ｓ９では、将来予測面ｉのトラッキング面位置ｚｉにおける鋳片表面温度予測値を、将来予測温度Ｔｐｒｅｄｉｊとする。
【００６１】
将来温度影響係数予測工程（以下において、「Ｓ１０」ということがある。）は、鋳造が進むことによりＳ８で設定した将来予測面が現在時刻からその下流側に隣接する将来予測面位置まで進む毎に、鋳造速度が現在時刻から変化しないと仮定して、各冷却ゾーンの冷却水量がステップ関数状に変化した場合の、各将来予測面がその下流側に隣接する将来予測面位置に到達するまでに通過する各トラッキング面位置における鋳片表面温度を予測し、この予測した鋳片表面温度とＳ９で予測した鋳片表面温度との偏差を求め、この偏差を用いて、ステップ関数状に変化する冷却水量に対する変化影響係数を求める工程である。
Ｓ１０では、各冷却ゾーンｋについて、現在時刻ｔで各冷却水量ｑｋ（ｔ）をステップ状にΔｑｋだけ変更した場合に、将来予測面ｉがその鋳造方向下流側に隣接する将来予測面の位置ｚｊに到達したときの鋳片表面温度Ｔｋｉｊを予測し、Ｓ９で求めたＴｐｒｅｄｉｊとの間の偏差ΔＴｋｉｊ（ｔ）＝Ｔｋｉｊ−ＴｐｒｅｄｉｊとΔｑｋとの関係を
【００６２】
【数１９】
2014140862.tif 000021
と表した時の係数Ｍｋｉｊを将来温度影響係数として求める。Ｓ１０では、各将来予測面について、ｊ行ｋ列成分に将来温度影響係数Ｍｋｉｊを並べた表面温度変化ゲイン行列Ｍｉを算出する。
【００６３】
鋳片表面参照温度算出工程（以下において、「Ｓ１１」ということがある。）は、Ｓ４で設定した鋳片表面温度の目標値と、Ｓ１０で予測した、将来予測面が将来予測面位置に到達した時点における鋳片表面温度の予測値との間の値である、時間に応じて決定される中間目標値（Ｓ１０の予測計算を繰り返すたびにＳ４で設定した鋳片表面温度の目標値に漸近する温度）である参照目標温度を算出する工程である。
Ｓ１１では、例えば、現在時刻において第ｉ冷却ゾーンの入り口にある断面の温度評価点ｚｊにおける参照目標温度Ｔｒｅｆｉｊは、下記式（２０）に示したように、将来予測温度Ｔｐｒｅｄｉｊと目標温度Ｔｔｇｔｊとの間を時間ｔｉｊの指数関数に従う比で内分する温度として定めることができ、Ｓ１１は、時間の関数で表される参照目標温度軌道Ｔｒｅｆｉｊ（ｔ）を求める工程、とすることができる。
【００６４】
【数２０】
2014140862.tif 000022
ここで、Ｔｒは予め定めた減衰パラメータに相当する時定数である。
【００６５】
最適化問題係数行列算出工程（以下において、「Ｓ１２」ということがある。）は、現在時刻ｔにおける各冷却ゾーンの冷却水量を決定変数とし、Ｓ９およびＳ１０の各々において各将来予測面が通過した各将来予測面位置における将来温度影響係数と、参照目標温度と鋳片表面将来予測温度との偏差を算出し、算出した当該偏差の各将来予測面に関する和を最小化する最適化問題の２次計画問題とし、この２次計画問題における決定変数に対する係数行列を算出する工程である。
Ｓ１２では、Ｓ１１の評価時刻ｔにおける各評価位置ｚｊの鋳片表面温度応答Ｔｐｒｅｄｉｊ（ｔ）＋ΔＴｉｊ（ｔ）と参照目標温度軌道Ｔｒｅｆｉｊ（ｔ）との偏差の重み付き二乗和と、各冷却ゾーンにおける冷却水量の変更ステップ幅Δｑｋの二乗和との合計を評価関数とし、この評価関数を最小化するΔｑ＝［Δｑ１  Δｑ２  …  ΔｑＫ］Ｔを求める。評価関数は下記式（２１）で表される。
【００６６】
【数２１】
2014140862.tif 000023
ここで、Ｔｐｒｅｄｉ、Ｔｒｅｆｉ、および、ΔＴｉは、それぞれ、式（２２）、式（２３）、および、式（２４）で表される。
【００６７】
【数２２】
2014140862.tif 000024
【００６８】
【数２３】
2014140862.tif 000025
【００６９】
【数２４】
2014140862.tif 000026
評価関数の温度偏差の項は、Ｓ１０で求めたゲイン行列を用いて下記式（２５）のように書き換えることができ、さらに、冷却水量の変更ステップ幅Δｑｋに無関係な項を除けば、上記評価関数の最小化は、下記式（２６）で表されるＪ’の最小化と等価である。
【００７０】
【数２５】
2014140862.tif 000027
【００７１】
【数２６】
2014140862.tif 000028
【００７２】
Ｊ’の最小化は、Δｑを決定変数とする２次計画問題である。ＱはＩ×Ｉ次元の非負定行列、ＲはＫ×Ｋ次元の正定行列である。例えば、Ｑには対角成分が負でない定数である対角行列などを用い、Ｒには対角成分が正の定数である対角行列などを用いる。さらに、冷却水量の変更ステップ幅の上限および下限や、冷却水量の上限および下限などに基づく制約条件を加えることにより、ミストスプレー２における物理的な制約を反映することができる。
【００７３】
最適化問題求解工程（以下において、「Ｓ１３」ということがある。）は、Ｓ１２における２次計画問題を数値的に解くことにより、現在時刻におけるΔｑの最適値Δｑ＊を求める工程である。上記２次計画問題は凸２次計画問題であるため、Δｑに制約がない場合、最適解Δｑ＊は、下記式（２７）で求められる。また、Δｑに制約がある場合には、有効制約法などを用いることにより、容易に最適解Δｑ＊を求めることができる。
【００７４】
【数２７】
2014140862.tif 000029
【００７５】
冷却水量変更工程（以下において、「Ｓ１４」ということがある。）は、Ｓ１３で求めた最適解Δｑ＊を、現在の冷却ゾーンの冷却水量ｑ（ｔ）へと加えることにより
【００７６】
【数２８】
2014140862.tif 000030
に変更する。このようにして変更された冷却水量ｑ（ｔ＋１）は、次回の制御周期で用いられる。
【００７７】
Ｓ１乃至Ｓ１４を有する本発明の冷却制御方法によれば、表面温度を評価するトラッキング面の鋳造方向下流側に隣接している冷却ゾーンの入口以外の位置にも、冷却水量の変更の影響をすぐに反映することができるので、鋳片全体の表面温度を、予め定めた目標温度に常に一致するように制御することが可能になる。したがって、本発明の冷却制御方法によれば、鋳片全体の表面温度を予め定めた目標温度に制御する際の精度を高めることが可能になる。鋳片全体の表面温度を精度良く目標温度に制御することにより、いかなる鋳造速度でも、また鋳造速度が鋳造中に変化した場合でも、連続鋳造機の曲げセグメントや矯正セグメントにおいて、表面温度を鋼の脆化域を回避するように制御することが可能になるので、表面疵による欠陥のない鋳片を製造することが可能になる。
【００７８】
以上説明した本発明の冷却制御方法は、例えば、図５に示した冷却制御装置１０を用いて実施することができる。図５に示したように、冷却制御装置１０は、鋳片表面温度測定部として機能する鋳片表面温度計７と、鋳造速度把握部として機能する鋳造速度測定ロール８と、トラッキング面設定部１０ａと、鋳片目標温度設定部１０ｂと、温度固相率推定部１０ｃと、熱伝達係数推定部１０ｄと、伝熱凝固モデルパラメータ修正部１０ｅと、将来予測面設定部１０ｆと、将来予測部１０ｇと、将来温度影響係数予測部１０ｈと、鋳片表面参照温度算出部１０ｉと、最適化問題係数行列算出部１０ｊと、最適化問題求解部１０ｋと、冷却水量変更部１０ｌと、を有している。上述のように、鋳片表面温度計７はＳ１で用いられ、鋳造速度測定ロール８はＳ２で用いられる。また、トラッキング面設定部１０ａではＳ３が、鋳片目標温度設定部１０ｂではＳ４が、温度固相率推定部１０ｃではＳ５が、熱伝達係数推定部１０ｄではＳ６が、伝熱凝固モデルパラメータ修正部１０ｅではＳ７が、それぞれ行われる。さらに、将来予測面設定部１０ｆではＳ８が、将来予測部１０ｇではＳ９が、将来温度影響係数予測部１０ｈではＳ１０が、鋳片表面参照温度算出部１０ｉではＳ１１が、それぞれ行われ、最適化問題係数行列算出部１０ｊではＳ１２が、最適化問題求解部１０ｋではＳ１３が、冷却水量変更部１０ｌではＳ１４が、それぞれ行われる。したがって、冷却制御装置１０を用いることにより、本発明の冷却制御方法を実施することができる。したがって、本発明によれば、鋳片全体の表面温度を、予め定めた目標温度に常に一致するように制御することが可能な、連続鋳造機の２次冷却制御装置を提供することができる。
【実施例】
【００７９】
スラブ用連続鋳造機において、鋳型出口直下の第１冷却ゾーンから最終の第１０冷却ゾーンまでを対象に、本発明を適用した実施例を以下に示す。
温度目標値は、鋳造速度一定と仮定して、各冷却ゾーン水量を最適化した場合のストランド伝熱凝固計算による、トラッキング面位置における鋳片表面温度計算値を用いた。本実施例で使用した連続鋳造機は、鋳片幅２３００ｍｍ、鋳片厚３００ｍｍ、鋳型内メニスカス位置から２次冷却帯出口までの距離２８．５ｍのスラブ用連続鋳造機である。本実施例における伝熱計算の更新間隔は２５ｍｍ、トラッキング面の間隔は１２５ｍｍ、将来温度予測面の間隔は１．２５ｍとした。トラッキング面では、鋳片の断面を長辺中心線および短辺中心線で分割した４分の１断面（図２参照）を、厚み方向に２０分割および幅方向に４０分割して、上記伝熱凝固モデルによる計算を行った。
なお、鋳片の鋳片表面温度の測定は、第４冷却ゾーン出側の、メニスカスから５．２５ｍ離れた位置で行い、鋳片長辺面中央で、放射温度計にて測定を行った。
【００８０】
［実施例１］
鋳造中に鋳込速度を２５％減少させた場合に本発明の冷却制御方法を適用した（実施例１）。実施例１における、各冷却ゾーン出口での鋳片幅方向中央部表面温度、各冷却ゾーンにおける冷却水量、および、鋳造速度と時間（時刻）との関係についての結果を、図６に示す。鋳造速度を０．８ｍ／ｍｉｎから０．６ｍ／ｍｉｎに急に減少させ、その５分後に０．８ｍ／ｍｉｎに戻した場合、実施例１における、各冷却ゾーン出口の鋳片表面温度と目標温度との二乗誤差平方根は、１２℃から１８℃の間であった。
一方、鋳造中に鋳込速度を２５％減少させた場合に従来の水量カスケード制御を適用したとき（比較例）の結果を、図７に示す。比較例では、実施例１と同じ条件で鋳造速度を変化させたにもかかわらず、各冷却ゾーン出口の鋳片表面温度と目標温度との二乗誤差平方根は、１７℃から２４℃であった。図６および図７に示したように、特に、鋳造速度を０．８ｍ／ｍｉｎから０．６ｍ／ｍｉｎへと低減した後、および、鋳造速度を０．６ｍ／ｍｉｎから０．８ｍ／ｍｉｎに戻した後における第１冷却ゾーンから第５冷却ゾーンまでの冷却水量の制御を比較すると、図６に示した実施例１では、図６に示した比較例よりも、第１冷却ゾーンから第５冷却ゾーンの冷却水量が冷却ゾーン出口の鋳片表面温度と目標温度との差が少なくなるように、好適な形でずれている様子が確認された。この結果から、本発明によれば、鋳造速度を変更しても、鋳片の表面温度を目標温度に高精度に制御可能であることが確認された。
【００８１】
［実施例２］
鋳造中に第３冷却ゾーンの温度目標値を２０℃低下させるように変更した場合に本発明の冷却制御方法を適用した（実施例２）。なお、この目標温度とは将来予測工程で予測される鋳片表面温度が近付くべき目標値である。実施例２における、鋳片表面温度の実績値、目標温度、冷却水量、および、鋳造速度と時間との関係についての結果を、図８に示す。
図８に示したように、温度目標値を低下させた後、第３冷却ゾーンの冷却水量は次第に増加させた結果、第３冷却ゾーンの出口における鋳片表面温度は、２０℃低下させた変更後の目標温度に漸近した。これに対し、温度目標値を低下させた後に、第４冷却ゾーンの冷却水量を若干低減することにより、第４冷却ゾーンの入口における鋳片温度の低下を補償した。その結果、第４冷却ゾーンの出口における鋳片表面温度の変化幅は３℃に抑制された。すなわち、本発明によれば、鋳片の表面温度を目標温度に高精度に制御可能であることが確認された。
なお、実施例２では、第３冷却ゾーンよりも鋳造方向の上流側に位置する第１冷却ゾーンや第２冷却ゾーンにおける冷却水量、および、温度には変化がなかった。そのため、第１冷却ゾーンおよび第２冷却ゾーンの結果の図示は省略し、第３冷却ゾーンおよび第４冷却ゾーンの結果のみを図示した。
【００８２】
［実施例３］
事前の冷却水量計算で設定された冷却水量で冷却すると、第４冷却ゾーンの出口における鋳片表面温度が目標温度よりも１６℃高くなると予想された際に、本発明の冷却制御方法によって実際の熱伝達係数を逐次推定しながら第４冷却ゾーンの冷却水量を調整した（実施例３）。実施例３における、鋳片表面温度の実績値、目標温度、冷却水量、および、鋳造速度と時間との関係についての結果を、図９に示す。
図９に示したように、第４冷却ゾーンでは冷却水量を当初の設定値よりも増大させるように制御され、その結果、第４冷却ゾーンの出口における鋳片表面温度を目標値に一致させることができた。この結果から、本発明によれば、鋳片の表面温度を目標温度に高精度に制御可能であることが確認された。
なお、実施例３では、第３冷却ゾーンよりも鋳造方向の上流側に位置する第１冷却ゾーンや第２冷却ゾーンにおける冷却水量、および、温度には変化がなかった。そのため、第１冷却ゾーンおよび第２冷却ゾーンの結果の図示は省略し、第３冷却ゾーンおよび第４冷却ゾーンの結果のみを図示した。
【符号の説明】
【００８３】
１…鋳型
２…ミストスプレー
３…流量調整弁
４…溶鋼メニスカス
５…鋳片
６…冷却ゾーン境界線（入口または出口位置）
７…鋳片表面温度計
８…鋳造速度測定ロール
９…連続鋳造機
１０…冷却制御装置
１０ａ…トラッキング面設定部
１０ｂ…鋳片目標温度設定部
１０ｃ…温度固相率推定部
１０ｄ…熱伝達係数推定部
１０ｅ…伝熱凝固モデルパラメータ修正部
１０ｆ…将来予測面設定部
１０ｇ…将来予測部
１０ｈ…将来温度影響係数予測部
１０ｉ…鋳片表面参照温度算出部
１０ｊ…最適化問題係数行列算出部
１０ｋ…最適化問題求解部
１０ｌ…冷却水量変更部
【図１】
2014140862.tif 000031
【図２】
2014140862.tif 000032
【図３】
2014140862.tif 000033
【図４】
2014140862.tif 000034
【図５】
2014140862.tif 000035
【図６】
2014140862.tif 000036
【図７】
2014140862.tif 000037
【図８】
2014140862.tif 000038
【図９】
2014140862.tif 000039