一种中厚板轧件断面形状的计算方法技术

一种中厚板轧件断面形状的计算方法，属于轧制技术领域，包括以下步骤：（１）对上支撑辊辊身长度进行划分，计算得到轧件的等份宽度和左右两侧剩余部分，求解上辊系变形；（２）对所求解的上辊系变形，判断上工作辊轴线位移是否满足精度要求，如果满足要求转入步骤（３），如果不满足重新求解上辊系变形，（３）计算下工作辊挠曲位移，（４）计算轧件断面厚度分布，轧件的断面厚度分布可近似用二次曲线来表示，计算得到的数据利用最小二乘法进行回归即可确定表示轧件断面形状的二次曲线的系数。本发明专利技术能够自适应轧制过程中轧件宽度的变化，对轧件使用较少的划分份数仍可以得到满意的计算结果，提高了运算速度。

【技术实现步骤摘要】

本专利技术属于轧制
，涉及了一种中厚板轧件断面形状高精度的求解方法。
技术介绍
中厚板轧件的断面形状直接影响了轧件的跑偏、侧弯，是板型预测、分析的重要影响因素。断面形状与轧制过程中轧机的辊系变形密切相关，传统的辊系变形分析方法基本思想是将轧辊离散成若干个单元，同时将轧辊所承受的载荷及变形也按相同的单元离散化。 轧机轧制过程中辊系的变形主要包括四个部分①与轧件接触引起的工作辊压扁②工作辊轴线的挠度③工作辊与支撑辊之间的压扁④支撑辊轴线的挠度。采用弹性基础梁模型，将支撑辊看成弹性基础，工作辊与支撑辊之间看成具有弹性常数为K1的弹簧，轧件引起的工作辊压扁看作是具有弹性常数K2弹簧。 考虑到中厚板轧机上、下辊系的相似性，为简化起见，以上辊系为研究对象，上辊系(包括上工作辊、上支撑辊和轧机牌坊)模型的建立如附图说明图1所示，PL，PR-左右两侧轧制力；YBL，YBR-上支撑辊左右两侧轴承中心位置在垂直方向的位移；YWL，YWR-上工作辊左右两侧轴承中心位置在垂直方向的位移；JL，JR-左右两侧弯辊力；p(j)-轧件和上工作辊之间的压力分布；q(j)-上工作辊与上支撑辊之间的压力分布；L2-上支撑辊辊身长度；L4-上支撑辊两侧轴承中心距离；L′4-上工作辊两侧轴承中心距离；L5-上支撑辊辊身边缘与上支撑辊轴承中心的距离；L′5-上支撑辊辊身边缘与上工作辊轴承中心的距离。 传统求解轧件断面形状的方法是先将上支撑辊全长分割为m个等份，认为上工作辊与上支撑辊之间的压力分布q(j)和上工作辊与轧件之间的压力分布p(j)在每个分割区是均布的。上辊系的力平衡方程和变形的位移协调方程如下 1)作用于上工作辊、上支撑辊垂直方向的力平衡方程 上工作辊 上支撑辊 式中，J为上工作辊弯辊力，JL为上工作辊左侧弯辊力，JR为上工作辊右侧弯辊力，N是板边处的分割区域编号，NL是板边处的左侧分割区域编号，NR是板边处的右侧分割区域编号，YR是上支撑辊轴承中心垂直方向的位移，YRL是上支撑辊左侧轴承中心垂直方向的位移，YRL是上支撑辊右侧轴承中心垂直方向的位移，K3是刚度系数，考虑轴承、压下螺丝及牌坊的变形，下标L和R分别表示左侧和右侧的值。 2)上工作辊和上支撑辊接触部分的位移协调方程 其中，K1为工作辊和支撑辊之间的弹性常数，YWL和YWR分别是上工作辊左右两侧轴承中心位置在垂直方向的位移；YR(i)和YW(i)是第i分割区中心处上支撑辊和上工作辊轴线位移；αB(i，j)和αW(i，j)分别为上支撑辊和上工作辊轴线位移影响函数；RCB(i)和RCW(i)分别为上支撑辊和上工作辊第i单元的轧辊凸度。 3)力矩平衡方程 上工作辊 上支撑辊 使用以上公式(1)～(7)求解辊系变形时，由于对上支撑辊辊身全长分割为m个等份并未考虑到轧件的宽度是否能正好等于Δz的整数倍，只能四舍五入近似求解，随着轧制过程中轧件宽展的变化，计算结果会发生跳变，影响了计算精度。
技术实现思路
针对利用传统方法计算轧件断面形状所遇到的问题，本专利技术提出了一种对轧件宽度划分后左右两侧剩余部分进行处理的方法，本方法可以自适应宽度的变化，消除了使用传统方法为减少边部的影响提高划分数量而影响计算速度的问题，提高了计算速度和计算精度，可进行工程上的在线应用。 本专利技术采用的技术方案包括以下步骤 (1)对上支撑辊辊身长度进行划分，计算得到轧件的等份宽度和左右两侧剩余部分，按照公式(2)、(3)、(4)、(7)、(13)、(14)和(15)生成计算矩阵，求解上辊系变形，包括计算上工作辊与上支撑辊之间的压力q(1)，q(2)，...，q(m)、上支撑辊两侧轴头位移YRL和YBR及上工作辊两侧轴头位移YWL和YWR；上辊系包括上工作辊、上支撑辊和轧机牌坊。 先将上支撑辊辊身全长分割为m个等份，对轧件等份划分后的两侧剩余部分单独进行处理，叠加到原有的辊系平衡方程和变形的位移协调方程中，不改变原有模型结构，使计算过程统一化。板宽方向的划分情况见图2所示。N1，N2-轧件占据辊身划分单元的开始、结束编号；Δz-每个划分单元的长度；ΔL-划分后左侧剩余的长度；ΔR-划分后右侧剩余的长度。 上支撑辊身长为L2，划分为m个等份，每个等份的宽度为 假设轧件钢板宽度为W，咬钢后偏离轧辊中心线ΔS，按照图2的划分情况，轧件占据辊身划分单元的开始、结束编号N1、N2以及划分后两侧剩余部分ΔL、ΔR。 其中，λL、λR为整数， 将得到的轧件划分后两侧剩余部分加入到辊系平衡方程和变形的位移协调方程中，即对轧件划分后的部分分为两类，一类为Δz的整数倍，使用传统的辊系平衡方程和变形的位移协调方程求解，另一类为两侧划分后的剩余部分，单独处理后叠加到原来方程中。对以上相关公式(1)、(5)、(6)进行修改，修改后的公式见(13)、(14)和(15)。 式中K1为工作辊和支撑辊之间的弹性常数；PL，PR-左右两侧轧制力；YBL，YBR-上支撑辊左右两侧轴承中心位置在垂直方向的位移；YWL，YWR-上工作辊左右两侧轴承中心位置在垂直方向的位移；JL，JR-左右两侧弯辊力；p(j)-轧件和上工作辊之间的压力分布；q(j)-上工作辊与上支撑辊之间的压力分布；L2-上支撑辊辊身长度；L4-上支撑辊两侧轴承中心距离；L′4-上工作辊两侧轴承中心距离；L5-上支撑辊辊身边缘与上支撑辊轴承中心的距离；L′5-上支撑辊辊身边缘与上工作辊轴承中心的距离；YB(i)和YW(i)是第i分割区中心处上支撑辊和上工作辊轴线位移；αB(i，j)和αW(i，j)分别为上支撑辊和上工作辊轴线位移影响函数；RCB(i)和RCW(i)分别为上支撑辊和上工作辊第i单元的轧辊凸度； 以上公式中包括的未知数为q(1)，q(2)，...，q(m)，YBL，YBR，YWL，YWR，方程总数为m+4个。将上述m+4个方程式列成矩阵形式，利用高斯消去法进行求解，即可求出与上辊系变形相关的变量。上工作辊与上支撑辊之间的压力q(1)，q(2)，...，q(m)、上支撑辊两侧轴头位移YBL和YBR及上工作辊两侧轴头位移YWL和YWR； (2)对所求解的上辊系变形，判断上工作辊轴线位移是否满足精度要求，如果满足要求转入步骤(3)，如果不满足重新求解上辊系变形。 上辊系变形的求解过程的收敛条件为ε＜δ，即最近两次求解的上工作辊轴线位移差小于指定的值，其中 δ取1.0×10-5mm (3)计算下工作辊挠曲位移，在板坯轧制过程中，下辊系(包括下工作辊和下支撑辊)会发生对称于轧制中心线的挠曲变形，如图3所示，h(i)-初始辊缝分布；h(i)′-轧制过程中的辊缝分布，YW(i)-上工作辊轴线位移；Δz-单元划分长度。根据计算出的上工作辊轴线位移为YW(i)，去除上工作辊轴线倾斜位移后即可得到下工作辊的挠曲位移 (4)计算轧件断面厚度分布，轧件的断面厚度分布可近似用二次曲线来表示，将公式(18)计算得到的数据利用最小二乘法进行回归即可确定表示轧件断面形状的二次曲线的系数。 根据设定辊缝h(i)、本文档来自技高网...

【技术保护点】
一种中厚板轧件断面形状的计算方法，其特征是该方法包括以下步骤：　（１）对上支撑辊辊身长度进行划分，计算得到轧件的等份宽度和左右两侧剩余部分，求解上辊系变形，包括计算上工作辊与上支撑辊之间的压力ｑ（１），ｑ（２），…，ｑ（ｍ）、上支撑辊两侧轴头位移Ｙ↓［ＢＬ］和Ｙ↓［ＢＲ］及上工作辊两侧轴头位移Ｙ↓［ＷＬ］和Ｙ↓［ＷＲ］；　Ｋ↓［３Ｌ］Ｙ↓［ＢＬ］＋Ｋ↓［３Ｒ］Ｙ↓［ＢＲ］＝＊ｑ（ｊ）Δｚ　（２）　Ｙ↓［Ｂ］（ｉ）＋ｑ（ｉ）／Ｋ↓［１］－Ｙ↓［Ｗ］（ｉ）＝－Ｒ↓［ＣＷ］（ｉ）＋Ｒ↓［ＣＢ］（ｉ）／２　（３）　Ｙ↓［Ｂ］（ｉ）＝＊ｑ（ｊ）Δｚα↓［Ｂ］（ｉ，ｊ）＋（Ｙ↓［ＢＲ］－Ｙ↓［ＢＬ］）［Ｌ↓［５］＋（ｉ－０．５）Δｚ］／Ｌ↓［４］＋Ｙ↓［ＢＬ］　（４）　Ｋ↓［３Ｒ］Ｌ↓［４］Ｙ↓［ＢＲ］＝＊ｑ（ｊ）Δｚ［Ｌ↓［５］＋（ｊ－０．５）Δｚ］　（７）　＊ｑ（ｊ）Δｚ＝＊ｐ（ｊ）Δｚ＋ｐ（Ｎ↓［１］－１）ΔＬ＋ｐ（Ｎ↓［２］＋１）ΔＲ＋Ｊ↓［Ｌ］＋Ｊ↓［Ｒ］　（１３）　Ｙ↓［Ｗ］（ｉ）＝＊ｐ（ｊ）Δｚα↓［Ｗ］（ｉ，ｊ）＋ｐ（Ｎ↓［１］－１）ΔＬα（ｉ，Ｎ↓［１］－１）＋ｐ（Ｎ↓［２］＋１）ΔＲα（ｉ，Ｎ↓［２］＋１）－＊ｑ（ｊ）Δｚα↓［Ｗ］（ｉ，ｊ）＋（Ｙ↓［ＷＲ］－Ｙ↓［ＷＬ］）［Ｌ′↓［５］＋（ｉ－０．５）Δｚ］／Ｌ′↓［４］＋Ｙ↓［ＷＬ］　（１４）　＊ｐ（ｊ）Δｚ［Ｌ′↓［５］＋（ｊ－０．５）Δｚ］＋ｐ（Ｎ↓［１］－１）ΔＬ［（Ｎ↓［１］－１）Δｚ－ΔＬ／２］＋ｐ（Ｎ↓［２］＋１）ΔＲ（Ｎ↓［２］Δｚ＋ΔＲ／２）＋Ｌ′↓［４］Ｊ↓［Ｒ］＝＊ｑ（ｊ）Δｚ［Ｌ′↓［５］＋（ｊ－０．５）Δｚ］　（１５）　式中：Ｋ↓［１］为工作辊和支撑辊之间的弹性常数；Ｐ↓［Ｌ］，Ｐ↓［Ｒ］－左右两侧轧制力；Ｙ↓［ＢＬ］，Ｙ↓［ＢＲ］－上支撑辊左右两侧轴承中心位置在垂直方向的位移；Ｙ↓［ＷＬ］，Ｙ↓［ＷＲ］－上工作辊左右两侧轴承中心位置在垂直方向的位移；Ｊ↓［Ｌ］，Ｊ↓［Ｒ］－左右两侧弯辊力；ｐ（ｊ）－轧件和上工作辊之间的压力分布；ｑ（ｊ）－上工作辊与上支撑辊之间的压力分布；Ｌ↓［２］－上支撑辊辊身长度；Ｌ↓［４］－上支撑辊两侧轴承中心距离；Ｌ′↓［４］－上工作辊两侧轴承中心距离；Ｌ↓［５］－上支撑辊辊身边缘与上支撑辊轴承中心的距离；Ｌ′↓［５］－上支撑辊辊身边缘与上工作辊轴承中心的距离；Ｙ↓［Ｂ］（ｉ）和Ｙ↓［Ｗ］（ｉ）是第ｉ...

【技术特征摘要】

【专利技术属性】
技术研发人员：何纯玉，吴迪，赵宪明，刘相华，胡贤磊，
申请(专利权)人：东北大学，
类型：发明
国别省市：89[中国|沈阳]