棒材轧制过程数学模型(共19页)

1、1轧制过程数学模型1.1轧制工艺参数模型随着科学技术的发展，计算机已广泛应用于轧钢生产过程的控制，促使轧钢生产向自动化、高速和优质方向发展。电子计算机在线控制生产过程，不仅仅只是电子计算机本身的硬件和软件的作用，更重要的是控制系统和各种各样的数学模型，正因为有适合轧钢生产的各种数学模型，才有可能实现电子计算机对整个轧钢生产各个环节的控制，获得高精度的产品。线材连轧生产过程的主要内容基本上可归纳为尺寸变化和温度变化两大类性质极不相同但又相互紧密联系的物理过程，涉及的数学模型主要是轧制工艺参数的制定、各环节的温降变化、产品质量控制及实现线材连轧生产的可靠性等。在线材连轧生产过程中，准确地计算（预估

2、）各个环节的温度变化是实现计算机控制的重要前提，这是因为轧件各道次的变形阻力、轧制压力、轧制力矩的准确确定与温度是分不开的，而各机架轧制压力的预估精度将直接关系到设备的使用安全等。 下面分别讨论线材连轧生产过程中的温降模型、变形阻力模型、轧制力与轧制力矩模型。1.1.1延伸系数及孔型尺寸计算模型在制订棒线材轧制工艺时，当坯料和产品断面面积F0和Fn给定之后，总延伸系数就可唯一确定：其中：n总轧制道次；i某一道次的延伸系数；Fi某一道次的轧件断面面积。对于箱型孔，轧件断面面积可通过下式计算：箱形孔示意图对于椭圆孔，轧件断面面积可通过下式计算：椭圆孔示意图对于圆孔，轧件断面面积可通过下式计算：圆孔

3、示意图1.1.2前滑模型孔型轧制时的前滑率计算可采用筱篬或斋藤提出的实验模型。两者都认为前滑仅是轧件、孔型几何尺寸的函数。斋藤模型以平均工作辊径定义前滑，当道次变形量较小时会出现负前滑的计算结果；筱篬模型改用孔型槽底处的最小辊径定义前滑，即前滑Sf为： Sf=V1/VR-1 （1.1）其中：V1 ，VR轧件出口速度及孔型槽底处的轧辊线速度。在孔型中轧制时，前滑值取平均值，其计算式为式中 变形区中性角的平均值； 咬入角的平均值； 摩擦角，一般为2127度； 轧辊工作直径的平均值； 轧件轧前、轧后高度的平均值；1.1.3轧件温降模型轧件在轧制过程中的温度变化，是由辐射、传导、对流引起的温降和金属变

4、形所产生的温升合成的，可用下式表示： （1.5）以上四项起作用的是辐射损失和金属变形热所产生的温升。各项温度变化的计算按下式进行：1、由于辐射引起的温降计算 （1.6）式中：Tf辐射引起的温降，；F轧件的散热表面积，m2；t冷却时间，s；T轧件表面绝对温度，K。2、由于传导引起的温降计算 （1.7）式中：TZ传导引起的温降，； 钢材的导热系数，1.255KJ/(m·h·)；Fz轧件与导热系体的接触面积，m2，对于轧辊Fz2lcbc×106；lc轧件与轧辊的接触弧长，mm；bc轧件轧前与轧后的平均宽度，mm；c0钢材平均比热容，在若杂货温度西热轧温度下取c00.62

5、7KJ/Kg· tz传导时间，s； hc轧件轧前与轧后的平均高度，mm。3、由于对流引起的温降计算 （1.8）式中 Td对流引起的温降，； T轧件表面绝对温度，K； T0环境绝对温度，K； V0轧件的移动速度，m/s； t对流时间，s； r轧件表面的相对黑度，r0.8；Tf同时间内的辐射温降，；4、由于变形热产生的温升计算 （1.9）式中Tb变形热产生的温升，； A该道次所需变形功，根据公式A=pVln(H/h)； P平均单位压力，MPa，粗略估计可用p=(ty0t75)×b/1500计算； V轧件体积，mm3H、h轧件轧前、轧后高度，mm； a系数，表明被轧件吸收的变形能

6、的相对部分，在T/Ty00.4时，当静力变形时（102s1）为0.92.6；当动力变形时（102s1）为1921； b强度极限，MPa； ty0钢材的熔点温度，K。取钢材的密度7.8，则得： （1.10）由于传导和对流引起的温降很小，甚至可以忽略不计。此时可以采用.采利柯夫方法计算在孔型中轧制和移送到下一孔型时间内，轧件温度得变化：（1.11）式中t0进入该孔型前得轧件温度，；l轧后轧件横截面周边长，mm；F轧后轧件横截面面积，mm2；t轧件冷却时间，s；t1在该孔型中金属温度得升高，；t1值按下式确定： t1=0.183Kmln（1.13）式中Km金属塑性变形抗力，MPa；延伸系数。1.1.

7、4变形抗力模型钢铁材料在热状态下的物理特性，与其温度、化学成分、应力、应变状态等诸多因素有关。目前在这方面的研究还不够充分，对于大多数钢种，只能给出离散数据的描述；但对于碳钢，平均变形抗力(MPa/mm2)可按以下模型计算： （1.14）其中：简单应力状态下的材料热变形抗力。 0.28exp（） （TTd） 0.28g（C,t）exp() (TTd)（1.15）式中：参数T，Td按下式计算： T= Td=0.95 (1.16)其中：t0轧件温度；C材料的碳含量百分数。式（2）中的函数g(C,t)为： （1.17）式（1）中的fm为考虑材料应变量等因素的影响系数： （1.18）在孔型设计时，式（

8、1）、式（5）中的平均应变和平均应变速率按下式计算： ln （1.19）其中：F0，FH轧件入口断面面积和轧件被孔型压掉部分的断面面积；F轧辊平均工作半径；N轧辊转数，r/min ；Lm接触弧平均长度。式（1）中的指数m由下式给出： (0.019C0.126)T0.075C0.05 （TTd） （0.081C0.154）T0.019C0.2070.027/（C0.32）（TTd） （1.20）该模型的适用范围为:材料的碳含量小于1.2，温度在(700-1200)范围内，平均应变小于0.7，平均应变速率在(0.1100)s1范围之内，该模型的最大优点是其数学上的完整性，利于实现计算机编程计算，但

9、也因此使其精确性有所降低，可用来进行预报。另外，为了适应计算机在线控制轧钢生产对变形阻力数学模型的要求和进一步提高计算精度，周纪华等采用碳钢和合金在高温、高速下测定得到的变形温度、变形速度和变形程度对变形阻力影响的大量实测数据而建立了非线性回归模型。它是以各种钢种为单位，得到各回归系数值，结构如下式。 （1.21）式中 ；基准变形阻力，即1000、0.4和10s1时的变形阻力（MPa）；变形温度（）；变形速度（s-1）；变形程度（对数应变）；， （1.23） （1.24）式中 v1轧件出口速度； l变形区长度； R轧辊半径； h0轧前高度。、回归系数，其值取决于钢种。各回归系数值按钢种的分类列

10、于表1-1至表1-14。表1-1普通碳钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a6Q215150. 0-2.7933.5560.2784-0.24600.42321.468Q235150.6-2.8783.6650.1861-0.12160.37951.402Q235-F140.3-2.9233.7210.3102-0.26590.45541.520表1-2优质碳素结构钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a608F136.13.3874.3120.51300.53200.58871.87908AL136.82.9993

11、.8180.35520.31860.49961.74210151.42.7713.5280.11470.03530.45371.59320152.72.6093.3210.20980.13320.38981.45445158.82.7803.5390.22620.15690.34171.379表1-3低合金钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a616Mn156.72.7233.4460.25450.21970.46581.56616MnCu160.12.4273.0900.06370.03870.40051.49914MnMoV177.12.6943.4

12、290.26160.24450.41571.49920Mn136.53.0573.8920.37430.31940.43371.51520MnSi163.02.4943.1740.06530.02380.42471.49210Ti161.22.5273.2170.15200.08390.40901.46015Ti171.12.0712.6400.14570.08400.36981.92610CrNi5MoV161.22.9223.7200.24510.20860.37521.36210CrNi2MoV153.12.9193.7160.26520.23790.40421.41928Cr2Ni2M

13、o154.83.0573.8920.22200.16970.37921.38430CrSiMo159.92.8333.6700.16270.09450.34541.33712Mn160.92.7443.4930.22700.18650.44331.54312MnNb164.52.6823.4140.12160.05080.40791.46312MnNb164.92.5323.2240.12090.04900.38461.42312MnNb164.12.6453.3670.18060.12870.40211.46712Mn*164.72.5413.2340.21860.18250.48011.5

14、2912MnNb*152.72.2702.8900.09440.03690.49271.57812MnNb*154.02.3633.0080.22640.19820.51931.63212MnNb*157.42.3422.9810.13780.09330.56701.844Y12CaS168.62.4543.1240.41220.44370.66382.364Y45CrCaS157.62.9333.7340.21180.13780.57282.034Si15AQ154.32.6213.3360.23940.20720.66312.403注：1.表示Nb含量不同。 2.*在进行变形抗力试验时，将

15、试件加热到1250，再冷却到8501150进行压缩得到变形抗力的试验数据。表1-4合金结构钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a609Mn2162.12.7103.4490.22480.17310.49421.67827SiMn176.13.1140.9650.38080.36440.57881.99420MnSiTi168.62.9593.7660.26310.22140.38111.45820MnV156.32.7103.4500.14990.07400.38411.34640MnSiV173.82.9973.8150.29260.23380.364

16、01.46315CrA165.62.8463.6230.21890.17720.55071.79920Cr148.42.4943.1750.18780.13090.43271.46940Cr153.42.8393.6140.17310.10500.35701.38338CrSiA186.63.2624.1530.19280.12580.18581.12230CrMnSiA168.73.0003.8180.29530.25080.42301.52535CrMnSiA189.73.1313.9850.18600.11260.18961.13020CrMnTi164.42.7293.4720.150

17、70.07080.39801.42315CrMoA170.92.7313.4770.14920.09060.60102.03320CrMoA163.32.7293.4730.12150.05830.62832.13835CrMoA184.23.2224.1010.15650.07650.19001.09342CrMo153.02.9823.7960.25080.18100.41601.48320CrMnMoA173.72.9253.7230.23860.19480.44191.49440CrMnMo173.03.1323.9870.19600.12310.41831.48412CrMoVA19

18、4.02.8243.5960.27290.25870.22821.07925Cr2MoVA172.73.0313.8580.12710.06460.44101.51045CrNiMoVA178.93.2714.1640.21470.15450.41671.52838CrMoAlA180.23.0913.9340.25400.21710.42581.49820MnVB163.72.7253.4700.16530.09280.54291.79720Mn2TiB162.02.7793.5370.09550.01280.31441.29740MnB155.63.0553.8900.16110.0646

19、0.23011.18045MnMoB163.22.9033.6950.28510.24560.49201.64340CrNiA160.83.1804.0480.23330.16830.38631.41518Cr2Ni4WA183.32.7953.5580.27530.26620.51351.66612CrNi4A165.92.8723.6560.25260.22000.52691.70323MnNiCrMo177.72.7923.5540.28800.26930.49311.650表1-5碳素工具钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a6T18A135.4

20、2.9473.7510.03400.10890.27811.312T10A141.83.0613.8960.00200.16980.61012.352T12A129.03.2164.0940.16910.04150.35911.564T13A133.13.5104.4680.28510.17710.43081.618表1-6合金工具钢变形抗力数学模型回归系数钢 种回归系数0/Mpaa1a2a3a4a5a69SiCr171.43.5864.5650.30330.23770.46861.739CrMn156.13.8684.9240.39570.36210.55102.123CrWMn140.83

21、.1864.0560.09790.01700.51731.916Cr6WV231.63.4314.3680.21420.15950.39381.6169Mn2157.62.8153.5830.32240.28120.57181.164W6Mo5cRv2290.12.6113.3240.29650.29060.24901.370表1-7弹簧钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a665Mn147.52.8403.6160.11530.01220.36661.49360Si2MnA172.83.2364.1180.25880.19450.59562.17550

22、CrVA159.53.2094.0850.29650.24050.44351.554表1-8轴承钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a6GCr9149.6-3.4864.4370.27170.18170.39891.574GCr15153.43.6854.6910.44800.41210.38881.589Cr2Ni4A202.82.9733.7850.31550.32140.48021.65420CrNiMoA170.23.1083.9560.25030.22150.47611.572表1-9不锈耐酸钢变形抗力数学模型回归系数钢 种回 归 系 数0/M

23、paa1a2a3a4a5a61Cr13179.41.9662.5410.05710.02950.52191.7812 Cr13184.32.3452.9860.00060.10600.59202.0103 Cr13196.83.1203.9720.28470.25270.46031.6354 Cr13213.43.2034.0770.27070.24320.52971.9041Cr18Ni9Ti224.62.2582.8750.35270.37450.32251.278表1-10硅钢变形抗力数学模型回归系数钢 种回 归 系 数0/Mpaa1a2a3a4a5a6D31093.732.7090.1

24、8800.08280.16131.150将此模型的计算结果与经典的变形抗力曲线图对比后发现，当变形速度在130s-1，变形温度在8501200之间时，结果能够很好的吻合曲线。因此在编制程序进行计算时采用的便是此模型。1.1.5轧制压力模型工程计算中经常采用如下简化的专用于孔型轧制的轧制压力公式计算轧制压力： （1.25）式中：平均变形抗力；接触投影面积；确定轧件与轧辊的接触面积，经常采用如下公式：用矩形箱形孔，方六角，六角方，方平椭圆，平椭圆方以及矩形平辊系统轧制时 （1.26）按方椭轧制方案时 （1.27）按椭椭，椭圆，圆椭，椭立椭和立椭椭轧制时 （1.28） 椭圆圆 （1.29）圆椭圆 （

25、1.30）载荷系数，针对各种孔型轧制情况的值回归模型为：式中：W考虑不同轧制条件的无量纲参数； 其中：分别为轧件如出口断面面积。1.1.6轧制力矩及功率模型轧制力矩计算公式为： （1.31）式中：轧制压力平均接触弧长度力臂系数力臂系数也采用对各种孔型轧制情况的回归模型：（1.34）轧制功率是单位时间所做的功，即： （1.35）式中：A变形功，KJ； t轧制时间，s。又由轧制所消耗的功与轧制力矩之间的关系为： （1.36）式中：角度，rad； 角速度，rad/s； R轧辊半径，mm； V轧辊线速度，m/s。得：将上式用工程上常用的参数和质量单位表示为：（KW） （1.37）式中 M轧制力矩，t&

26、#183;m；n轧辊转速，r/m。1.2工艺参数的计算 以产品规格10mm，钢种45#，终轧速度为16m/s，为例，各参数的计算过程如下。1.2.1轧件断面面积的计算在已知延伸系数的前提下，只需根据成品的断面面积依次往上推便可，对于10钢，成品断面面积F2075.39mm2，则F19F20×2075.39×1.22492.28 mm2，其它依次类推。在不知延伸系数的前提下，可根据孔型的尺寸通过几何换算的方法计算轧件的断面面积。 第一孔（箱形孔）的尺寸如下表所示h(mm)B(mm)b(mm)R(mm)r(mm)S(mm)117.0 173.2 152.0 23.4 18.0

27、25则其轧件截面面积0.2261.34418691mm2其它箱形孔略。第七孔（椭圆孔）的尺寸如下表所示h(mm)B(mm)R(mm)S(mm)52.8 97.2 65.86 10 则其轧件截面面积1.669.973849mm2其它椭圆孔略。第八孔（圆孔）的尺寸如下表所示D1(mm)B(mm)R(mm)(o)S(mm)60.0 65.82 30.00 306 则其轧件截面面积2923 mm2其它圆孔略。1.2.2轧件断面周长的计算1、对于箱型孔，轧件周长可按下式计算C=2(hk+Bk)则第一孔的周长C1=2×(117+173.2)=580.4mm，其余箱型孔雷同。2、对于椭圆孔，轧件周长可按下式计算则第五孔的周长277.41mm，其余椭圆孔雷同。3、对于圆孔，轧件周长可按下式计算C=bk则第六孔的周长C679×3.1415926248.19mm，其余圆孔雷同。1.2.3轧制速度的计算1、工作辊径的计算要设定轧辊转速需知道工作辊径，工作辊径按以下公式计算： 式中 DK轧辊工作辊径，mm； D0轧辊名义直径，mm； S辊缝，mm； hP轧件出口断面平均高度