轧制过程数学模型

1轧制过程数学模型

1.1轧制工艺参数模型

随着科学技术的发展，计算机已广泛应用于轧钢生产过程的控制，促使轧钢生产向自动化、高速和优质方向发展。电子计算机在线控制生产过程，不仅仅只是电子计算机本身的硬件和软件的作用，更重要的是控制系统和各种各样的数学模型，正因为有适合轧钢生产的各种数学模型，才有可能实现电子计算机对整个轧钢生产各个环节的控制，获得高精度的产品。

线材连轧生产过程的主要内容基本上可归纳为尺寸变化和温度变化两大类性质极不相同但又相互紧密联系的物理过程，涉及的数学模型主要是轧制工艺参数的制定、各环节的温降变化、产品质量控制及实现线材连轧生产的可靠性等。

在线材连轧生产过程中，准确地计算（预估）各个环节的温度变化是实现计算机控制的重要前提，这是因为轧件各道次的变形阻力、轧制压力、轧制力矩的准确确定与温度是分不开的，而各机架轧制压力的预估精度将直接关系到设备的使用安全等。

下面分别讨论线材连轧生产过程中的温降模型、变形阻力模型、轧制力与轧制力矩模型。

1.1.1延伸系数及孔型尺寸计算模型

在制订棒线材轧制工艺时，当坯料和产品断面面积F 0和F n 给定之后，总延伸系数∑μ就可唯一确定：

n

n n i i n i F F

F F F F F F F F 011211021==

=-+∑ μμμμμ 其中：n ——总轧制道次；

μi ——某一道次的延伸系数； F i ——某一道次的轧件断面面积。

椭圆孔示意图

mB R F +-=)sin (2θθ

R

B 2arcsin

2=θ ⎪⎭⎫ ⎝

⎛

--=2cos 12θR h m

对于圆孔，轧件断面面积可通过下式计算：

圆孔示意图

απθ2-=

αθtan 422R R F +=

1.1.2前滑模型

孔型轧制时的前滑率计算可采用筱篬或斋藤提出的实验模型。两者都认为前滑仅是轧件、孔型几何尺寸的函数。斋藤模型以平均工作辊径定义前滑，当道次变形量较小时会出现负前滑的计算结果；筱篬模型改用孔型槽底处的最小辊径定义前滑，即前滑S f 为：

S f =V 1/V R -1 （1.1） 其中：

V 1 ，V R ——轧件出口速度及孔型槽底处的轧辊线速度。 在孔型中轧制时，前滑值取平均值f S ，其计算式为

()[]

1cos 1cos ++-=

h

h

D S f γγ

⎪⎪⎭

⎫

⎝⎛+=

βααγ212 D

h H --

=1cos α 式中 γ——变形区中性角的平均值； α——咬入角的平均值； β——摩擦角，一般为21~27度； D ——轧辊工作直径的平均值；

h H ,——轧件轧前、轧后高度的平均值； 1.1.3轧件温降模型

轧件在轧制过程中的温度变化，是由辐射、传导、对流引起的温降和金属变形所产生的温升合成的，可用下式表示：

b d Z f T T T T T ∆-∆+∆+∆=∆（1.5）

以上四项起作用的是辐射损失和金属变形热所产生的温升。各项温度变化的计算按下式进行：

1、由于辐射引起的温降计算

4

1000072.0⎪⎭⎫ ⎝⎛=∆T G Ft T f （1.6）

式中：

ΔT f ——辐射引起的温降，℃；

F ——轧件的散热表面积，m 2； t ——冷却时间，s ； T ——轧件表面绝对温度，K 。 2、由于传导引起的温降计算 c

Z

Z Z Gh c t F T 08.1λ=

∆（1.7）

式中：

ΔT Z ——传导引起的温降，℃；

λ——钢材的导热系数，λ≈1.255KJ/(m ·h ·℃)；

Fz ——轧件与导热系体的接触面积，m 2，对于轧辊Fz ＝2l c b c ×10－6； l c ——轧件与轧辊的接触弧长，mm ； b c ——轧件轧前与轧后的平均宽度，mm ；

c 0——钢材平均比热容，在若杂货温度西热轧温度下取

c 0＝0.627KJ/Kg ·℃

tz ——传导时间，s ；

hc ——轧件轧前与轧后的平均高度，mm 。

3、由于对流引起的温降计算

f r d T T T T t V T T T ∆⎪⎭⎫ ⎝⎛-⎥⎥⎦

⎤⎢⎢⎣

⎡+-∆4

05.22

301003.0ε（1.8） 式中

ΔTd ——对流引起的温降，℃； T ——轧件表面绝对温度，K ； T 0——环境绝对温度，K ； V 0——轧件的移动速度，m/s ； t ——对流时间，s ；

εr ——轧件表面的相对黑度，εr ≈0.8；

ΔTf ——同时间内的辐射温降，℃； 4、由于变形热产生的温升计算 G

c a A T b 0427)

1(-=∆（1.9）

式中

ΔTb ——变形热产生的温升，℃；

A ——该道次所需变形功，根据公式A=pVln(H/h)；

P ——平均单位压力，MPa ，粗略估计可用p=(t y0－t －75)×σb /1500计算；

V ——轧件体积，mm 3

H 、h ——轧件轧前、轧后高度，mm ；

a ——系数，表明被轧件吸收的变形能的相对部分，在T/T y0＞0.4时，

当静力变形时（102s －1）为0.9％～2.6％； 当动力变形时（102s －1）为19％～21％；

σb ——强度极限，MPa ； t y0——钢材的熔点温度，K 。

取钢材的密度γ＝7.8，则得：

)/ln()1(184.0h H a p T b -=∆（1.10）

由于传导和对流引起的温降很小，甚至可以忽略不计。此时可以采用А.И.采利柯夫方法计算在孔型中轧制和移送到下一孔型时间内，轧件温度得变化：

2732731000

0255.01000

3

3

100+⎪⎪⎭

⎫ ⎝⎛+∆++-

=∆t t F lt t T （1.11）

式中

t 0——进入该孔型前得轧件温度，℃； l ——轧后轧件横截面周边长，mm ； F ——轧后轧件横截面面积，mm 2； t ——轧件冷却时间，s ；

Δt 1——在该孔型中金属温度得升高，℃； Δt 1值按下式确定：

Δt 1=0.183K m ln μ（1.13） 式中

K m ——金属塑性变形抗力，MPa ； μ——延伸系数。 1.1.4变形抗力模型

钢铁材料在热状态下的物理特性，与其温度、化学成分、应力、应变状态等诸多因素有关。目前在这方面的研究还不够充分，对于大多数钢种，只能给出离散数据的描述；但对于碳钢，平均变形抗力(MPa/mm