JP4460493B2 - 圧延機におけるワークロールのサーマルクラウン計算方法 - Google Patents
圧延機におけるワークロールのサーマルクラウン計算方法 Download PDFInfo
- Publication number
- JP4460493B2 JP4460493B2 JP2005172016A JP2005172016A JP4460493B2 JP 4460493 B2 JP4460493 B2 JP 4460493B2 JP 2005172016 A JP2005172016 A JP 2005172016A JP 2005172016 A JP2005172016 A JP 2005172016A JP 4460493 B2 JP4460493 B2 JP 4460493B2
- Authority
- JP
- Japan
- Prior art keywords
- work roll
- roll
- rolled material
- heat transfer
- amount
- 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.)
- Active
Links
Images
Landscapes
- Control Of Metal Rolling (AREA)
Description
熱間圧延において板クラウンおよび板形状設定計算に使用するサーマルクラウン予測を高精度に算出するための計算方法に関して、同予測モデルに使用するロールバイト内熱伝達係数を高精度に算出するための方法に関する。
熱間圧延において高精度な板クラウン形状制御を実現するためには、圧延材の幅方向板厚分布に直接的に影響を及ぼすワークロールサーマルクラウン予測の高精度化が重要である。これまで、サーマルクラウン予測に関しては、例えば、特許文献１や特許文献２に示すように、熱間圧延で実用的な計算時間でかつ高精度な予測を行うモデルの提案が行われている。
上記予測モデルでサーマルクラウンを計算する際に必須であるワークロールへの入熱および抜熱にかかわる熱伝達係数は、サーマルクラウン量やワークロール温度の変化からその量を表現するように、あらかじめ設定値として与えている。圧延条件(板厚スケジュール、圧延速度など)やワークロール温度によって、ワークロールへの入熱が異なることは知られており、これにかかわる熱伝達係数を設定値として扱った場合、サーマルクラウン予測に誤差を生じる。
これに対して、例えば特許文献３に示すように、ロールバイト内の熱伝達係数を接触弧長の関数として与え、サーマルクラウン予測モデルに反映した予測方法が提案されている。
特許文献３による場合、ロールバイト内の熱伝達係数を、圧延ワークロール円周方向の表面温度とワークロール温度の関係から算出し、それを圧延条件で統計的に処理し、ロールバイト内における熱伝達係数のモデル化を行っている。
ロールバイト内では、ワークロールと圧延材とは密接に接触しており、ワークロールへの入熱量を評価する場合、すなわち、ロールバイト内での熱伝達係数を正確に算出するためには、圧延材からワークロールへの熱移動量を正確に計算したうえで、ロールバイト内の熱伝達係数を算出する必要がある。
しかしながら、特許文献３による場合、ワークロール表面のみの温度のみに着目し、圧延材からワークロールへの熱移動量を正確に評価したとは言い難く、そのため、サーマルクラウン予測にも誤差が生じるという問題があった。
一方、ロールバイト内でのワークロールと圧延材との接触時間は、極短時間であることから、ワークロールおよび圧延材における熱移動量変化(温度変化)は、例えば、非特許文献１(鉄鋼協会出版の『板圧延の理論と実際』の第6章で示されているように、ワークロールでは円周方向および圧延材では長手方向の熱伝導を無視しても、計算しても十分な精度であることが知られており、このことから、ワークロール半径方向および圧延材板厚方向への熱移動量計算(温度計算)は、非定常１次元熱伝導計算で十分であることが記載されている。
特開平07-80517号公報
特開平07-32041号公報
特開平2002-346620号公報
『板圧延の理論と実際』日本鉄鋼協会出版、1984
本発明は、上述した従来技術の問題を解決し、板クラウンおよび板形状設定計算に使用するサーマルクラウン予測を高精度に算出するための計算方法を提供することを目的としている。
上記目的を達成するために本発明の要旨とするところは、以下のとおりである。
（１）圧延機におけるワークロールのサーマルクラウン計算方法において、
前記予測モデルに用いるワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数を、ワークロールと圧延材とが接触状態にあるとし、ロールバイト内での圧延材からワークロールへの熱移動量(ロール半径方向温度分布)を算出し、ついで、この算出した熱移動量に基づくロールバイト入側と出側でのワークロール半径方向の熱移動の変化量(ロール内温度分布の変化量)から算出することを特徴とする圧延機におけるワークロールのサーマルクラウン計算方法である。
前記予測モデルに用いるワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数を、ワークロールと圧延材とが接触状態にあるとし、ロールバイト内での圧延材からワークロールへの熱移動量(ロール半径方向温度分布)を算出し、ついで、この算出した熱移動量に基づくロールバイト入側と出側でのワークロール半径方向の熱移動の変化量(ロール内温度分布の変化量)から算出することを特徴とする圧延機におけるワークロールのサーマルクラウン計算方法である。
（２）前記ロールバイト内での圧延材からワークロールへの熱移動量(ロール半径方向温度分布)を、ワークロールと圧延材との境界における摩擦発熱および圧延材スケール層、圧延材加工発熱、およびワークロールおよび圧延材熱物性値、ワークロールと圧延材が接触している時間を考慮して算出することを特徴とする上記（１）の圧延機におけるワークロールのサーマルクラウン計算方法である。
以上説明してきたように、本発明による場合、実用的な計算時間で、かつ高精度に圧延中のロールバイト内熱伝達係数を算出でき、高精度なサーマルクラウン予測が可能となることから、圧延中の圧延材の板クラウンおよび形状を安定的に製造することができる。
本発明は、サーマルクラウン予測の高精度化のために、ワークロールへの入熱にかかわる熱伝達係数を精度よく算出する方法について提供するものである。
そのため、サーマルクラウン予測モデルに関しては、例えば、特許文献１に示すように、異なる物性値を有する２層、すなわち芯材と外層材とで構成されたワークロールのロール内温度分布を圧延操業中に計算する際に、ロール半径方向温度分布の数式表現として、芯材と外層材との境界で、双方の熱伝導率の違いにより、熱流束の連続性から必然的に生じるロール半径方向温度勾配の不連続を表現可能な多項式を用いてサーマルクラウン計算を行うモデルを使用すればよく、そのときに使用するワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数を、本発明に基づく計算値に置き換えればよい。
また、本発明で計算される上記ワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数は、上記で示したサーマルクラウン予測モデルに限定されることはなく、他の予測モデルにもこの計算方法が適用できる。
上記の[背景技術]でも述べたように、ロールバイト内でワークロールへの入熱を評価する場合、圧延材からワークロールへの熱移動量を正確に計算する必要がある。
また、ロールバイト内でのワークロールと圧延材との接触時間は、極短時間であることから、ワークロールおよび圧延材における熱移動量変化(温度変化)は、例えば、非特許文献１の第6章で示されているように、ワークロールでは円周方向および圧延材では長手方向の熱伝導を無視しても、計算しても十分な精度であることが知られており、このことから、ワークロール半径方向および圧延材板厚方向への熱移動量計算(温度計算)は、非定常１次元熱伝導計算で十分であることが記載されている。
本発明によるワークロールへの熱移動量の計算は、ワークロールと圧延材との非定常１次元熱伝導計算を採用して算出する。このとき、本発明ではワークロールと圧延材は高圧で接触していることから、熱移動量計算は、ワークロールと圧延材は完全に接触状態にあるとして計算した。
本発明で実施した計算手順を以下に示す。図２には圧延中のロールバイト内でのワークロールと圧延材との接触状態を示した模式図、および本発明の計算方法を説明するための図を示す。このとき、上記非定常１次元熱伝導計算を行う場合、図２(a)で示したロールバイト入側での符号１で示す１線が圧延時間経過とともにロールバイト出側の符号５で示す５線になるという計算を行えばよい。
このとき、図２(b)に示したワークロールと圧延材との接触時間Δｔは、圧延条件から算出することができ、例えば次式で表現される。
ここで、Ｒ’(m)はヒッチコックのロール偏平半径、ＶＲ(m/sec)はワークロールの周速度、Ｒ(m)はワークロール半径、Ｈ、ｈ(m)は入出側板厚、ｐ(kgf/m)は圧延線荷重、Ｅ(kgf/m2)、νはワークロールのヤング率、ポアソン比である。
図２(b)に基づいた非定常１次元熱伝導方程式を式(1)に、境界条件を式(2)に示す。
ここで、ρ(y) (kg/m3)は密度、ｃ(y)(J/(kg・℃))は比熱、λ(y)(W/(m・℃))は熱伝導率、Ｔ(℃)は温度、ｔ(sec)は時間、ｑf(W/m２)は単位時間あたり、単位面積あたりの摩擦発熱、ｑｐ(W/m3)は単位時間あたり、単位体積あたりの加工発熱となり、ｙはワークロール半径および圧延材板厚方向であり、ここではｙ≧０でワークロール、ｙ＜０で圧延材を意味することになる。
ここで、Ｒはワークロール半径、Ｈ(t)は圧延時間とともに変化する板厚、λRおよびλSは、ワークロールと圧延材との境界面(y=0)でのロール側熱伝導率および圧延材側熱伝導率、ＴR、ＴSは上記境界面での温度である。
式(1)、(2)におけるρ(y)ｃ(y)、λ(y)は、
ｙ≧０でワークロールの熱物性値、
ｙ＜０で圧延材の熱物性値を使用すればよい。
ｙ≧０でワークロールの熱物性値、
ｙ＜０で圧延材の熱物性値を使用すればよい。
また、ワークロールと圧延材との境界にスケール層がある場合、その熱物性値をｙ＝０~−hs(ｈsはスケール厚)の位置に入れて計算すればよい。
さらに、式(２)のｑfおよびｑｐは、例えば、非特許文献１の第６章に記載されている方法で計算でき、例えば、次式で表現できる。
ｑf ＝ ｋf・In(Ｈ/ｈ)/Δｔ (5)
ｑｐ ＝ μ・ｋf・Ｖr (6)
ここで、 ｋfは２次元変形抵抗、μは摩擦係数、Ｖrはワークロールと圧延材との相対速度の絶対値の平均値である。
ｑｐ ＝ μ・ｋf・Ｖr (6)
ここで、 ｋfは２次元変形抵抗、μは摩擦係数、Ｖrはワークロールと圧延材との相対速度の絶対値の平均値である。
また、図２(a)における符号１で示す１線のロールバイト入側温度の圧延材側は、例えば、非特許文献１の第６章で記載されている既設モデルからの圧延材板厚方向温度分布の計算値、ワークロール側は、例えば、特許文献１で計算されるアイドル終了時(ワークロールと圧延材とが接触する直前における特許文献１の計算値)のワークロール半径方向温度分布を用いればよい。上記熱伝導方程式、境界条件、ロールバイト入側の温度初期条件などを用いれば、ロールバイト出側でのワークロール半径方向および圧延材板厚方向の温度分布を計算することができる。ここで、本発明では温度計算に有限要素解析を用いたが、差分法など熱伝導計算が可能な手法を用いてもよい。図２(c)には、ロールバイト入側と出側とのワークロール半径および圧延材板厚方向温度分布の計算結果を示す。
このとき、説明を簡単にするために、図２(b)での時間分割は４としているが計算精度維持のためには、少なくとも２０分割以上が好ましい。
ついで、ワークロールへの熱移動量(温度分布)の計算結果から、ワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数αＲ(W/(m2・℃))は、図２(c)に示したロールバイト入側と出側でのワークロール半径方向の熱移動の変化量(ロール内温度分布の変化量)を用いれば、次式で算出することができる。このとき、ワークロール半径方向温度分布は、本発明では有限要素解析で行われることから、ロール半径方向にｎ分割されたかたちで計算される。前記熱伝達係数は、その計算結果を用いて計算すればよい。
ここで、ＴR ｅ(y)およびＴR ｄ(y)はワークロール入側および出側での半径方向温度分布、ρRおよびｃRはワークロールの密度および比熱、ＴS ｅおよびＴR ｅはロールバイト入側での圧延材およびワークロールの表面温度、ΔyiはI分割目の要素長さ、Δｔは式(1)で示したワークロールと圧延材との接触時間である。また、前記ロールバイト入側と出側でのワークロール半径方向の熱移動の変化量は、式(７)の(ＴR d(y)i-ＴR e(y)i)項にあたることがわかる。
本発明を用いてロールバイト内の熱伝達係数を計算すれば、図２(c)で示すようなロールバイト内での急激な熱移動量を正確に表現できることから、圧延条件の違い(板厚、圧下率、圧延材温度、材質(変形抵抗)など)から必然的に生じるロールバイト内での熱伝達係数を正確に表現することが可能となることがわかる。
本発明の具体的な実施形態を(実施例)で詳しく説明する。
図３には本発明を適用したサーマルクラウン計算の計算フローを示す。具体的には、図３に示すように、発明の実施するための最良の形態で説明したように、ロールバイト内での非定常１次元熱伝導計算によるワークロールへの熱移動量の計算結果から、ロールバイト内での熱伝達係数を算出する。ついで、本発明では、例えば特許文献１で示したサーマルクラウン予測モデルの境界条件であるロールバイト内の熱伝達係数を、上記本発明から計算される熱伝達係数に置き換えてサーマルクラウン計算を行うことになる。
図１(a)には、本発明を熱間仕上圧延機列のF4スタンド（F4−stand）を適用したときのロールバイト内熱伝達係数の推移を示す。図１(b)には、ワークロールセンタ位置の熱膨張量の推移を実測値(F4出側板厚計より計算した値)と本発明を適用して予測したそれとの比較結果および圧延終了後の４分後にワークロール胴長方向のサーマルクラウン量の実測値と本発明を適用したサーマルクラウン予測との比較結果を示す。圧延条件の違い(図１(c))から、図１(a)示すように、本発明で計算したロールバイト内熱伝達係数は大きく異なることがわかる。本発明を適用すれば、上述したように圧延条件が異なることで生じるロールバイト内熱伝達係数の違いをサーマルクラウン予測モデルの中に反映することができることから、図１(b)に示すように、圧延初期から終了まで高精度にワークロールセンタの熱膨張量を表現でき、かつワークロール胴長方向プロフィルとして比較した結果でも高精度に予測可能であることがわかる。すなわち、本発明は、圧延条件の違い(板厚、圧下率、圧延材温度、材質(変形抵抗)など)から必然的に生じるロールバイト内での熱伝達係数を正確に表現可能であることがわかる。
本発明の実施形態を熱間仕上圧延機列に適用した結果を示したが、本発明は厚板圧延や冷間圧延にも適用可能であることは言うまでもない。
Claims (2)
- 圧延機におけるワークロールのサーマルクラウン計算方法において、
前記予測モデルに用いるワークロールと圧延材との接触部(ロールバイト内)の熱伝達係数を、ワークロールと圧延材とが接触状態にあるとし、ロールバイト内での圧延材からワークロールへの熱移動量(ロール半径方向温度分布)を算出し、ついで、この算出した熱移動量に基づくロールバイト入側と出側でのワークロール半径方向の熱移動の変化量(ロール内温度分布の変化量)から算出することを特徴とする圧延機におけるワークロールのサーマルクラウン計算方法。 - 前記ロールバイト内での圧延材からワークロールへの熱移動量(ロール半径方向温度分布)を、ワークロールと圧延材との境界における摩擦発熱および圧延材スケール層、圧延材加工発熱、およびワークロールおよび圧延材熱物性値、ワークロールと圧延材が接触している時間を考慮して算出することを特徴とする請求項１記載の圧延機におけるワークロールのサーマルクラウン計算方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005172016A JP4460493B2 (ja) | 2005-06-13 | 2005-06-13 | 圧延機におけるワークロールのサーマルクラウン計算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005172016A JP4460493B2 (ja) | 2005-06-13 | 2005-06-13 | 圧延機におけるワークロールのサーマルクラウン計算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2006346682A JP2006346682A (ja) | 2006-12-28 |
| JP4460493B2 true JP4460493B2 (ja) | 2010-05-12 |
Family
ID=37643064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005172016A Active JP4460493B2 (ja) | 2005-06-13 | 2005-06-13 | 圧延機におけるワークロールのサーマルクラウン計算方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4460493B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112238145B (zh) * | 2019-07-19 | 2022-09-02 | 北京宇轩智能科技有限公司 | 智能轧钢轧机精确咬钢信号处理方法 |
| CN114798755B (zh) * | 2022-04-26 | 2023-02-17 | 东北大学 | 一种分析冷连轧过程中热凸度对板形影响的方法 |
-
2005
- 2005-06-13 JP JP2005172016A patent/JP4460493B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2006346682A (ja) | 2006-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5839582B2 (ja) | 鋳型設計方法および鋳型 | |
| Zhang et al. | A semi-analytical method for the springback prediction of thick-walled 3D tubes | |
| JP5391205B2 (ja) | 制御装置 | |
| Zhang et al. | Modeling of rolling force of ultra-heavy plate considering the influence of deformation penetration coefficient | |
| Arif et al. | Roll deformation and stress distribution under thermo-mechanical loading in cold rolling | |
| Schmidtchen et al. | Fast numerical simulation of symmetric flat rolling processes for inhomogeneous materials using a layer model− part I: basic theory | |
| JP4460493B2 (ja) | 圧延機におけるワークロールのサーマルクラウン計算方法 | |
| Zhao et al. | Investigation of the interfacial heat transfer coefficient of sheet aluminum alloy 5083 in warm stamping process | |
| Wen | On a new concept of rotary draw bend-die adaptable for bending tubes with multiple outer diameters under non-mandrel condition | |
| Wen et al. | Investigations on the interfacial heat transfer coefficient during hot stamping of ultra-high strength steel with Al-Si coating | |
| Rout et al. | Finite element modeling of hot rolling: Steady-and unsteady-state analyses | |
| Arif et al. | Investigations on strain hardening during cutting of heat-resistant austenitic stainless steel | |
| JP7021608B2 (ja) | 連続鋳造鋳片の中心固相率推定方法 | |
| Hadała | Implementation of the heat balance in the finite element solution to the temperature field of the plastically deformed material | |
| Ragab et al. | Evaluation of estimates of roll separating force in bar rolling | |
| Kowalik et al. | Numerical and experimental estimation of forces during longitudinal rolling process of shaft formation | |
| Minutolo et al. | Dimensional analysis of a new type of groove for steel rebar rolling | |
| Steinheimer et al. | Thermal influences during rotary draw bending of tubes from stainless steel | |
| JP4710237B2 (ja) | 厚板鋼板の変形量予測方法及びその製造方法 | |
| Khan et al. | Thermal analysis of a cold rolling process—a numerical approach | |
| JP5644334B2 (ja) | 熱鍛型摩耗量予測装置 | |
| Hsu et al. | A three-dimensional inverse problem of estimating the surface thermal behavior of the working roll in rolling process | |
| JPS59125203A (ja) | 粗圧延鋼板の温度制御方法 | |
| JP6416676B2 (ja) | ロールプロフィール修正量の決定方法 | |
| Muvunzi et al. | Heat transfer in a hot stamping process: a review |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070904 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100127 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20100209 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100212 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 4460493 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130219 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130219 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140219 Year of fee payment: 4 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |