本发明专利技术涉及一种轧制力和轧制温度相互迭代的计算方法，初步计算轧制力后，根据所得轧制力计算相应的轧辊压扁半径和轧制温度，根据计算所得的轧制温度计算变形抗力、摩擦系数以达到根据轧制温度计算轧制力的目的，然后比较前后所计算的轧制力是否基本相等，若二者相差很小，也即达到了收敛条件，则迭代结束，取最后计算得到的轧制力、轧制温度作为计算结果。若不满足收敛条件，则根据所计算的轧制力再次计算轧辊压扁半径和轧制温度，然后根据轧制温度再次计算轧制力，并再次判断是否达到收敛条件。本发明专利技术考虑到了在现实生产过程中轧制力、轧制温度二者相互影响，存在着热力耦合的事实，通过数值迭代的方法，将轧制力和轧制温度的计算过程当做一个整体，还原轧制力和轧制温度相互影响的客观情况，以同时提高轧制力和轧制温度模型的预测精度。
【技术实现步骤摘要】
本专利技术涉及一种乳制力和乳制温度相互迭代的计算方法。
技术介绍
在钢板热连乳乳制过程中，乳制规程是制定热乳生产工艺最基本的内容之一。而 乳制规程的好坏，直接影响后续控制的难度和决定最终产品的质量，所以制定合理的乳制 规程对提高钢板带生产质量具有重大的意义。而要制定合理的乳制规程，准确的数学模型 是关键。其中，乳制力预测模型，乳制温度预测模型为乳制规程制定过程中的关键模型。准 确的预测在钢板乳制过程中乳制力和乳制温度的大小为合理制定乳制规程的最基本的保 障。在热连乳过程中，乳制温度直接影响着乳件在乳制过程中的变形抗力大小进而影响着 乳制力的大小，而另一方面，乳件因摩擦、塑性变形引起的温度变化多少又与乳制力大小直 接相关，所以，乳制温度和乳制力之间存在着热力耦合关系。而在现有技术中，在计算乳制 力和乳制温度时，都是将两个模型分开计算，人为的分离了乳制力和乳制温度的相互影响， 从而导致所预测的乳制力和乳制温度与实测值存在一定的误差。
技术实现思路
本专利技术的目的是为了克服现有技术的不足，提供了一种乳制力和乳制温度相互迭 代的计算方法。将将乳制力和温度计算当成一整体，还原乳制力和乳制温度相互影响的实 际情况，以达到同时提高乳制力和乳制温度预测精度的效果。 本专利技术是通过以下技术方案实现： 本专利技术所涉及的乳制力和乳制温度相互迭代的计算方法，实现步骤如下： 1、一种乳制力和乳制温度相互迭代的计算方法，其特征在于，包括以下步骤： 1)取计数器i=1，其中i为正整数且i彡1，计算考虑压扁后的变形区接触弧长 的水平投影& =^/^，_ ;预设定压扁后乳辊半径民'、平面变形抗力I以及摩擦系数yi 初始值，即Ri' =1?，1^=1(。，1^=y。，!?为工作辑原始半径，mm，Ah=H_h为压下量，H、h分别为板带入口、出口厚度，mm，K。为平面变形抗力初始值，MPa，y。为摩擦系数初始值； g+b 2)确定乳制力P1=BHQp，kN;式中，氧i =二^为变形区内乳件的平均宽度， A B、b分别为乳前、乳后板宽，mm;QP为平面应力状态系数； 3)根据步骤2)确定的乳制力Pi，确定压扁后乳辊半径心广卻+^)，mm;根 m 据压扁后乳親半径Ri+1 '计算接触弧长1#：^ ===^^TT; 4)根据步骤2)确定的乳制力Pi，计算乳制温度°C;其中，因塑性变形而产生 的温升°C;其中，p为乳件比重，kg/m3;c为乳件比热，X/(kg?m3)， nP为塑性变形热效率系数；摩擦温升°C;其中，nf为摩 擦热效率系数，f为前滑值，b为后滑值； 5)根据步骤4)所得到的乳制温度，计算变形抗力MPa， 根据变形抗力〇sl计算平面变形抗力K1+1=l. 155〇S1，其中，e为真应变；:&为变形速率，为热力学温度，K;al~a4为待定系数； 6)根据步骤4)所得到的乳制温度h，计算摩擦系数y1+1 = mKjO. 877-0. 0039^)，其中，&为与乳辊材质及表面状态有关修正系数，K2为与乳制材 料种类有关修正系数；K3为与乳制速度有关修正系数，K4为张力影响系数； 7)确定乳制力p1+1，其中 'Pi+iKi+iQpBmlj(i+i)， 8)判断是否满足收敛条件| (p1+1-Pl)/Pl | < 10 5; 9)若不满足收敛条件，则进行如下步骤： ①计蕖 ②计算乳制温度ti+1，其中，在计算乳制温度ti+1时， ③更新乳制力、乳制温度和压扁后乳辑半径，即Pi=pi+1，ti=ti+1，Ri+1' =Ri+2' ； 计算器i自增1，即i=i+1 然后返回步骤5); 10)若满足收敛条件，则乳制力p1+1作为最终乳制力，并根据此乳制力p1+1计算得 到最终乳制温度t1+1。 与现有技术相比，本专利技术有益的技术效果是：基于实际生产过程中乳制力和乳制 温度相互影响的客观事实，采用数值循环的方法，将乳制温度和乳制力的计算过程当做一 整体，来达到同时提高乳制力和乳制温度的技术效果，以减少计算误差，提高模型预测精 度。【附图说明】 图1为本专利技术的结构示意图。【具体实施方式】 为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对 本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并 不用于限定本专利技术。 请参阅图1，图1为本专利技术的计算流程图。 本专利技术所涉及的乳制力和乳制温度相互迭代的计算方法，实现步骤如下： 1、一种乳制力和乳制温度相互迭代的计算方法，其特征在于，包括以下步骤： 1)取计数器i=1，其中i为正整数且i彡1，计算考虑压扁后的变形区接触弧长 的水平投影4 ，mm;预设定压扁后乳辊半径民'、平面变形抗力I以及摩擦系数yi 初始值，即Ri' =1?，1^=1(。，1^=y。，!?为工作辑原始半径，mm，Ah=H_h为压下量，H、h分别为板带入口、出口厚度，mm，K。为平面变形抗力初始值，MPa，y。为摩擦系数初始值； 2)确定乳制力PfBUQpkN;式中，B=￥为变形区内乳件的平均宽度，B、 b分别为乳前、乳后板宽，mm;QP为平面应力状态系数； 3)根据步骤2)确定的乳制力Pi，确定压扁后乳辊半径4+1'=即+A)，mm;根 BnAh 据压扁后乳辊半径R1+1 '计算接触弧长#錢; 4)根据步骤2)确定的乳制力Pi，计算乳制温度ti，°C;其中，因塑性变形而产生 的温汧f°C;其中，P为乳件比重，kg/m3;c为乳件比热，J/(kg?m3)， nP为塑性变形热效率系数；摩擦温升°C;其中，nf为摩 擦热效率系数，f为前滑值，b为后滑值； 5)根据步骤4)所得到的乳制温度h，计算变形抗力根据变形抗力〇sl计算平面变形抗力K1+1=l. 155〇S1，其中，e为真应变；r为变形速率，为热力学温度，K;al~a4为待定系数； 6)根据步骤4)所得到的乳制温度h，计算摩擦系数y1+1 = mKjO. 877-0. 0039^)，其中，&为与乳辊材质及表面状态有关修正系数，K2为与乳制材 料种类有关修正系数；K3为与乳制速度有关修正系数，K4为张力影响系数； 7)确定乳制力p1+1，其中 'Pi+iKi+iQpBmlj(i+i)， 8)判断是否满足收敛条件| (p1+1-Pl)/Pl | < 10 5; 9)若不满足收敛条件，则进行如下步骤：①计算 ②计算乳制温度ti+1，其中，在计算乳制温度ti+1时， ③更新乳制力、乳制温度和压扁后乳辑半径，即Pi=pi+1，ti=ti+1，Ri+1' =Ri+2'； 计算器i自增1，即i=i+1 然后返回步骤5); 10)若满足收敛条件，则乳制力p1+1作为最终乳制力，并根据此乳制力p1+1计算得 到最终乳制温度t1+1。 以上所述仅为本专利技术的较佳实施例而已，并不用以限制本专利技术，凡在本专利技术的精 神和原则之内。【主权项】1. 一种乳制力和乳制温度相互迭代的计算方法，其特征在于，包括以下步骤： 1) 取计数器i = 1，其中i为正整数且i多1，计算考虑压扁后的变形区接触弧长的水 平投景，_ ;预设定压扁后乳辊半径民'、平面变形抗力Ki以及摩擦系数y ^刀 始值，即民' y。，!?为工作辑原本文档来自技高网...
【技术保护点】
一种轧制力和轧制温度相互迭代的计算方法，其特征在于，包括以下步骤：1)取计数器i＝1，其中i为正整数且i≥1，计算考虑压扁后的变形区接触弧长的水平投影mm；预设定压扁后轧辊半径Ri'、平面变形抗力Ki以及摩擦系数μi初始值，即Ri'＝R，Ki＝K0，μi＝μ0，R为工作辊原始半径，mm，Δh＝H‑h为压下量，H、h分别为板带入口、出口厚度，mm，K0为平面变形抗力初始值，MPa，μ0为摩擦系数初始值；2)确定轧制力Pi＝BmldiKiQp，kN；式中，为变形区内轧件的平均宽度，B、b分别为轧前、轧后板宽，mm；Qp为平面应力状态系数；3)根据步骤2)确定的轧制力Pi，确定压扁后轧辊半径mm；根据压扁后轧辊半径Ri+1'计算接触弧长4)根据步骤2)确定的轧制力Pi，计算轧制温度ti，℃；其中，因塑性变形而产生的温升℃；其中，ρ为轧件比重，kg/m3；c为轧件比热，J/(kg·m3)，ηP为塑性变形热效率系数；摩擦温升℃；其中，ηf为摩擦热效率系数，f为前滑值，b为后滑值；5)根据步骤4)所得到的轧制温度ti，计算变形抗力MPa，根据变形抗力σsi计算平面变形抗力Ki+1＝1.155σsi，其中，ε为真应变；为变形速率，s‑1；为热力学温度，K；a1～a4为待定系数；6)根据步骤4)所得到的轧制温度ti，计算摩擦系数μi+1＝K1K2K3K4(0.877‑0.0039ti)，其中，K1为与轧辊材质及表面状态有关修正系数，K2为与轧制材料种类有关修正系数；K3为与轧制速度有关修正系数，K4为张力影响系数；7)确定轧制力pi+1，其中，pi+1＝Ki+1QpBmld(i+1)；8)判断是否满足收敛条件|(pi+1‑pi)/pi|＜10‑5；9)若不满足收敛条件，则进行如下步骤：①计算Ri+2′=R(1+cPi+1BmΔh)]]>和ld(i+2)=ΔhRi+2′;]]>②计算轧制温度ti+1，其中，在计算轧制温度ti+1时，ΔtHF=2μi+1pi+1ld(i+2)×106(h1+h2)cρ×f2+b2f+bηf;]]>③更新轧制力、轧制温度和压扁后轧辊半径，即pi＝pi+1，ti＝ti+1，Ri+1'＝Ri+2'；计算器i自增1，即i＝i+1然后返回步骤5)；10)若满足收敛条件，则轧制力pi+1作为最终轧制力，并根据此轧制力pi+1计算得到最终轧制温度ti+1。...
【技术特征摘要】
【专利技术属性】
技术研发人员：范彬，熊志宏，谢强，尹旭妮，周超伦，
申请(专利权)人：湖南城市学院，
类型：发明
国别省市：湖南;43
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。