CN116371942A - 一种基于带钢横向强度不均的横向厚度分布预测方法 - Google Patents
一种基于带钢横向强度不均的横向厚度分布预测方法 Download PDFInfo
- Publication number
- CN116371942A CN116371942A CN202310652406.7A CN202310652406A CN116371942A CN 116371942 A CN116371942 A CN 116371942A CN 202310652406 A CN202310652406 A CN 202310652406A CN 116371942 A CN116371942 A CN 116371942A
- Authority
- CN
- China
- Prior art keywords
- strip steel
- transverse
- thickness distribution
- rolling
- transverse thickness
- 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.)
- Granted
Links
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 223
- 239000010959 steel Substances 0.000 title claims abstract description 220
- 238000009826 distribution Methods 0.000 title claims abstract description 110
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000005096 rolling process Methods 0.000 claims abstract description 91
- 238000004088 simulation Methods 0.000 claims abstract description 46
- 230000007246 mechanism Effects 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000008878 coupling Effects 0.000 claims abstract description 4
- 238000010168 coupling process Methods 0.000 claims abstract description 4
- 238000005859 coupling reaction Methods 0.000 claims abstract description 4
- 238000009864 tensile test Methods 0.000 claims abstract description 4
- 239000004973 liquid crystal related substance Substances 0.000 claims description 64
- 238000013000 roll bending Methods 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 18
- 230000009471 action Effects 0.000 claims description 13
- 101100379079 Emericella variicolor andA gene Proteins 0.000 claims description 12
- 238000002474 experimental method Methods 0.000 claims description 12
- 239000000835 fiber Substances 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 11
- 238000005452 bending Methods 0.000 claims description 10
- 230000001105 regulatory effect Effects 0.000 claims description 10
- 230000005489 elastic deformation Effects 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 101100001673 Emericella variicolor andH gene Proteins 0.000 claims description 3
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 3
- 230000007547 defect Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 description 9
- 238000005097 cold rolling Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 230000008844 regulatory mechanism Effects 0.000 description 1
- 238000005482 strain hardening Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B38/00—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product
- B21B38/04—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product for measuring thickness, width, diameter or other transverse dimensions of the product
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B38/00—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product
- B21B38/02—Methods or devices for measuring, detecting or monitoring specially adapted for metal-rolling mills, e.g. position detection, inspection of the product for measuring flatness or profile of strips
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明属于轧制过程自动化控制技术领域,具体涉及一种基于带钢横向强度不均的横向厚度分布预测方法,包括:截取部分带钢进行拉伸试验获取真实应力‑应变曲线;获取轧辊参数、轧制工艺参数以及轧制前后的带钢参数;建立关于带钢‑轧辊变形耦合分析的横向厚度分布仿真模型;利用横向厚度分布仿真模型进行模拟实验;构建板形执行机构的厚度调控功效系数计算模型,提取模拟实验稳定轧制阶段的带钢横向厚度分布数据,计算各板形执行机构的厚度调控功效系数;提取模拟实验稳定轧制阶段的带钢宽度数据和出口带钢横向厚度分布曲线,建立带钢横向厚度分布预测计算方程,输入板形执行机构的调控数值以获取对应横向厚度分布曲线。
Description
技术领域
本发明属于轧制过程自动化控制技术领域,涉及一种基于带钢横向强度不均的横向厚度分布预测方法。
背景技术
近年来冷轧带钢市场供求情况发生改变,国内外市场竞争日益激烈,冷轧带钢各项性能指标也朝着高精尖方向发展。板形是评价冷轧带钢质量的一个关键指标,主要从纵向和横向两个方面来衡量一般意义上所述的板形;纵向主要关注的是冷轧带钢平直度方面,而横向关注的是带钢板凸度以及边部减薄等问题,其中最常考虑的评价指标就是横向厚度分布,实际生产中若能降低带钢横向厚度偏差可提高带钢成材率,减少切损量。因此,想要研究带钢板形控制,有效的冷轧出口带钢横向厚度分布预测至关重要。
针对冷轧带钢横向厚度的控制问题,国内研究人员做了一些相关研究。中国期刊论文“基于ANSYS的冷轧带钢横向厚度的研究“以四辊板带轧机为研究对象,利用ANSYS软件建立了四辊冷轧板带的有限元模型,完整地模拟了板带的冷轧过程,根据仿真及实验结果,对比分析了板带冷轧过程中横向厚度的分布规律,为实际生产和理论研究提供了指导依据;中国博士论文“薄带钢冷轧过程的弹塑性有限元模拟”应用了四辊轧制模型对带钢轧后横截面厚度分布及凸度进行了模拟,分析了不同的轧制条件对轧后带钢横截面形状及中心、边部凸度的影响,得到了弯辊力、带钢变形抗力、张应力等对带钢板凸度及边部减薄的影响及弯辊力对改善板形特别是较宽带钢的显著作用。
上述研究所存在不足主要有三个方面:(1)带钢模型均设置为理想几何模型,但带钢横向厚度控制在实际的冷轧生产具有很强的复杂性,无论是来料还是轧制过程中的带钢横截面必为不规则截面,若在仿真过程中仍使用理想几何模型势必与实际生产情况不相贴合;(2)经过实验证明,带钢横向上的材料强度在实际生产过程中并非均匀分布,边部与中心处强度存在一定差别,这种差别进一步也会影响轧制力和轧辊的弹性变形情况,因此采用传统均质带钢模型会在带钢横向厚度分布的仿真分析中存在一定局限性,仿真得到的厚度调控功效系数曲线也存在一定的精度误差;(3)以往的有限元仿真仅给出了不同执行机构在单位变化量内变化对于带钢横向厚度分布的影响情况,未结合调控机理给出带钢横向厚度分布的快速预测方法。
发明内容
为解决上述技术问题,本发明的目的是提供一种基于带钢横向强度不均的横向厚度分布预测方法。
本发明提供一种基于带钢横向强度不均的横向厚度分布预测方法,包括:
步骤1:轧制前截取部分带钢,并对该段带钢横向不同位置处进行拉伸试验获取相应位置的真实应力-应变曲线;
步骤2:获取轧辊参数、轧制工艺参数以及轧制前后的带钢参数;
步骤3:根据步骤1获取的真实应力-应变曲线和步骤2获取的参数建立关于带钢-轧辊变形耦合分析的横向厚度分布仿真模型;
步骤4:制定针对不同板形执行机构的实验计划,利用横向厚度分布仿真模型对带钢轧制过程进行模拟实验;
步骤5:构建板形执行机构的厚度调控功效系数计算模型,提取步骤4中各模拟实验稳定轧制阶段的带钢横向厚度分布数据,计算各板形执行机构的厚度调控功效系数;
步骤6:提取步骤4中各模拟实验稳定轧制阶段的带钢宽度数据和出口带钢横向厚度分布曲线,基于出口带钢横向厚度分布曲线的拟合系数和厚度调控功效系数曲线的拟合系数建立带钢横向厚度分布预测计算方程,输入板形执行机构的调控数值以获取对应横向厚度分布曲线。
进一步的,所述步骤1具体为:
轧制前截取部分带钢,沿带钢宽度方向将带钢等间距分成N份带钢试样,对每份带钢试样进行拉伸实验,获取所有带钢试样的N条真实应力-应变曲线。
进一步的,获取真实应力-应变曲线具体为:
步骤1.1:定义真实应力-应变曲线由弹性变形阶段、均匀塑性变形阶段和局部塑性变形阶段组成;
步骤1.2:求出弹性变形阶段和均匀塑性变形阶段的分界点,即屈服点对应的屈服应力σs:
其中,P s是材料开始屈服时的载荷,由试验机载荷刻度盘读出;S 0是试样原始横截面面积;
步骤1.3:找出均匀塑性变形阶段各瞬时的流动应力即真实应力Y和对数应变∈:
其中,P为加载瞬间的载荷,由试验机载荷刻度盘读出;S为加载瞬间试样的横截面面积,由体积不变条件求出:
其中,l 0是试样标距的原始长度;Δl是试样标距的伸长量;同样塑性失稳点的应力和应变也由上式求出;
步骤1.4:到达均匀塑性变形阶段和局部塑性变形阶段的分界点,即塑性失稳点后试样开始出现颈缩,找出断裂时的真实应力Y K1及其对应的对数应变∈ K1:
其中,P K1是试样断裂时的载荷;S K1是试样断裂时的横截面面积;l K1是试样断裂时的标距总长度;
步骤1.5:使用下式对颈缩后的局部塑性变形阶段的曲线进行修正:
其中,Y K2是去除形状硬化后的真实应力;d是试件颈缩处直径;ρ是试件颈缩处外形曲率半径。
进一步的,所述步骤2具体为:
步骤2.1:采集轧辊参数,包括:轧辊直径、轧辊长度、轧辊密度、轧辊弹性模量以及轧辊泊松比;
步骤2.2:采集轧制工艺参数,包括:摩擦系数、轧制速度、前张力、后张力、压下率、工作辊弯辊力、中间辊弯辊力以及中间辊横移量;
步骤2.3:采集轧制前的带钢参数包括:带钢宽度和入口带钢横向厚度分布数据;将带钢两侧距边部100mm处设为边部减薄区,剩余部分设为中心区,并在带钢横向上以1.5mm为单位间隔取测点,分别测量中心区和两个边部减薄区每个测点的厚度值并连线获取入口带钢横向厚度分布曲线;采集轧制后的带钢参数包括:出口带钢横向厚度分布数据,轧制前后的带钢宽度视为相等。
进一步的,所述步骤3具体为:
步骤3.1:根据获取的轧辊数据建立轧机的三维弹塑性有限元模型;
步骤3.2:根据轧前带钢宽度数据和入口带钢横向厚度分布曲线建立带钢的三维弹塑性有限元模型;
步骤3.3:沿着宽度方向将带钢模型分成N个单元,并将步骤1中获取的N条真实应力-应变曲线写入对应位置的单元中,在仿真求解过程中,每个单元的变形过程均遵循对应的真实应力-应变曲线。
进一步的,所述步骤4具体为:
步骤4.1:根据获取的轧制工艺参数,制定针对于不同执行机构的实验计划,板形执行机构包括:中间辊横移、工作辊弯辊和中间辊弯辊;
步骤4.2:依据对照实验计划,更改横向厚度分布仿真模型中的中间辊横移量、工作辊弯辊力、中间辊弯辊力的取值,进行模拟仿真。
进一步的,所述步骤5具体为:
步骤5.1:建立以带钢宽度、长度、厚度为坐标轴的三维坐标系,以带钢的中心位置点作为坐标原点;将带钢假设为连续离散化的纵向纤维条,α点设为距带钢中心宽向位置x α 的离散化的纵向纤维条的离散点编号,设板形执行机构编号为i,其调节量为Δu i ;
步骤5.2:通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量:
通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量:
式中,H(x α )为选取的纵向纤维条的轧前横向厚度,h(x α )为选取的纵向纤维条的轧后横向厚度;Δh α,i (x α )为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量,为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量;
步骤5.3:通过带钢轧制前后横向厚度的相对改变量对厚度调控功效系数进行表示为:
式中,T α,i (x α )为第i种板形执行机构在x α 处的厚度调控功效系数;
步骤5.4:提取步骤4中各模拟实验稳定轧制阶段带钢横向厚度分布数据,根据步骤5.2计算获得位置x α 处的带钢轧制前后横向厚度的相对改变量,并根据步骤5.3计算获得各执行机构的厚度调控功效系数及其分布曲线。
进一步的,所述步骤6具体为:
步骤6.1:提取步骤4中各模拟实验稳定轧制阶段带钢宽度数据和出口带钢横向厚度分布曲线,对宽度数据进行归一化处理;并对得到的出口带钢横向厚度分布曲线使用五次多项式进行拟合,得到多项式拟合系数B 1、B 2、B 3、B 4、B 5;五次多项式如下:
其中,y为出口带钢横向厚度分布曲线;THK_Center为截距;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1];
步骤6.2:对步骤5.4中获取的三种板形执行机构的厚度调控功效系数曲线使用六次Legendre正交多项式进行拟合,获取三种板形执行机构的厚度调控功效系数曲线的拟合系数A 0、A 1、A 2、A 4、A 6,采用的六次Legendre正交多项式如下式:
其中,E(x)为带钢横向厚度分布曲线;A 0为常数项、A 1、A 2、A 4、A 6分别为厚度调控功效系数曲线的一次、二次、四次和六次项拟合系数,其绝对值大小表明厚度调控功效系数曲线中对一次、二次、四次和六次板形缺陷的调控分量;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1];e(x)为拟合误差;
步骤6.3:根据步骤2.3中划分的边部减薄区和中心区,以归一化的带钢宽度作为x轴、厚度预测值为y轴建立二维坐标系,建立带钢中心区每坐标点厚度预测值的计算方程为:
其中,H P 为计算得到的每坐标点厚度预测值;
其中,V WRB 为工作辊弯辊力,V IRB 为中间辊弯辊力,V IRS 为中间辊横移量,A 0W 、A 0I 和A 0S 都是常数项,A 1W 、A 2W 、A 4W 和A 6W 分别为工作辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1I 、A 2I 、A 4I 和A 6I 分别为中间辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1S 、A 2S 、A 4S 和A 6S 分别为中间辊横移调控功效系数曲线一次、二次、四次、六次项拟合系数;
步骤6.4: 建立带钢两个边部减薄区每坐标点厚度预测值的计算方程为:
其中,为带钢每个边部减薄区基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;/>和/>分别是带钢操作侧边部和驱动侧边部基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;H OS 和H DS 分别带钢中心区靠近操作侧最边部点厚度值和驱动侧最边部点厚度值;V e是边部减薄值;f t是边部减薄程度权重;f b是位置系数;
步骤6.5:对步骤6.1中得到的归一化宽度数据进行反归一化处理后分别代入步骤6.3和6.4的计算方程式,计算得到的每坐标点厚度值并获得带钢整体横向厚度分布曲线。
本发明的一种基于带钢横向强度不均的横向厚度分布预测方法,至少具有以下有益效果:
(1)首先,所提供的建立带钢模型完全依照采集的实际生产现场的入口带钢参数,使得带钢模型更加贴近实际。
(2)其次,本发明的建立的带钢模型的横向材料性能以现场实测的入口带钢横向厚度分布曲线建立,实现了基于带钢横向材料强度非均匀分布条件下的冷轧过程分析,进一步提高了出口带钢横向厚度分布的仿真精度。
(3)再次,采用本发明提供的方法可以预测冷轧出口带钢的横向厚度分布情况,为实际生产中带钢横向厚度分布的快速调整提供依据。
(4)最后,本发明采用的是基于轧制理论和三维弹塑性有限元仿真的模拟实验方法,可减少实际实验导致的设备、时间以及成本损耗。
附图说明
图1是本发明的一种基于带钢横向强度不均的横向厚度分布预测方法的流程图;
图2是真实应力-应变曲线示意图;
图3是在不同的工作辊弯辊力作用下,出口带钢的横向厚度分布曲线变化情况图;
图4是在不同的中间辊弯辊力作用下,出口带钢的横向厚度分布曲线变化情况图;
图5是在不同的中间辊横移量作用下,出口带钢的横向厚度分布曲线变化情况图;
图6是工作辊弯辊的厚度调控功效系数曲线;
图7是中间辊弯辊的厚度调控功效系数曲线;
图8是中间辊横移的厚度调控功效系数曲线。
具体实施方式
本实施例中以某厂的六辊冷轧机组为例,对其稳定轧制阶段的带钢横向厚度分布进行预测,轧机轧辊均为平辊。
如图1所示,本发明的一种基于带钢横向强度不均的横向厚度分布预测方法,包括:
步骤1:轧制前截取部分带钢,并对该段带钢横向不同位置处进行拉伸试验获取相应位置的真实应力-应变曲线;
具体实施时,轧制前截取部分带钢,沿带钢宽度方向将带钢等间距分成32份带钢试样,对每份带钢试样进行拉伸实验,获取所有带钢试样的32条真实应力-应变曲线,获取真实应力-应变曲线具体为:
步骤1.1:本实施例中获取的真实应力-应变曲线示意如图2所示,首先,定义真实应力-应变曲线由弹性变形阶段、均匀塑性变形阶段和局部塑性变形阶段组成。图中OC为弹性变形阶段,为均匀塑性变形阶段;/>和/>为局部塑性变形阶段;C为屈服点,是弹性变形阶段和均匀塑性变形阶段的分界点,此时的应力称为屈服应力;/>点为塑性失稳点,是均匀塑性变形阶段和局部塑性变形阶段的分界点,此时的应力称为抗拉强度,到达/>点后,试样开始出现颈缩;应力-应变曲线的末端/>和/>是破坏点,试样发生断裂,是塑性变形的终止点。
步骤1.2:求出屈服点对应的屈服应力σs:
其中,P s是材料开始屈服时的载荷,由试验机载荷刻度盘读出;S 0是试样原始横截面面积,单位m2。
步骤1.3:随后找出均匀塑性变形阶段各瞬时的流动应力即真实应力Y和对数应变∈:
其中,P为加载瞬间的载荷,由试验机载荷刻度盘读出;S为加载瞬间试样的横截面面积,由体积不变条件求出:
其中,l 0是试样标距的原始长度,单位mm;Δl是试样标距的伸长量,单位mm;同样塑性失稳点的应力和应变也由上式求出。
步骤1.4:到达塑性失稳点后试样开始出现颈缩,但颈缩后的局部塑性变形阶段的各加载瞬间横截面面积不能由体积不变条件求出,因此需找出断裂时的真实应力Y K1及其对应的对数应变∈ K1:
其中,P K1是试样断裂时的载荷;S K1是试样断裂时的横截面面积;l K1是试样断裂时的标距总长度。至此,可在真实应力-应变坐标系中确定出未修正的真实应力-应变曲线。
步骤1.5:由于颈缩出现后截面发生局部收缩,颈缩部位应力状态由单向拉应力转变为三向应力状态,增大拉伸应力会出现形状硬化,此时颈缩后的应力-应变曲线反映的是材料冷作硬化和形状硬化的总效应,此时需要去除形状硬化效应对真实应力-应变关系曲线的影响,使用下式对颈缩后的局部塑性变形阶段的曲线进行修正:
步骤2:获取轧辊参数、轧制工艺参数以及轧制前后的带钢参数,所述步骤2具体为:
步骤2.1:采集轧辊参数,包括:轧辊直径、轧辊长度、轧辊密度、轧辊弹性模量以及轧辊泊松比;
具体包括如下参数:工作辊辊身直径、工作辊辊身长度、工作辊密度、工作辊弹性模量、工作辊泊松比;中间辊副辊直径、中间辊副辊长度、中间辊辊身直径、中间辊辊身长度、中间辊密度、中间辊弹性模量、中间辊泊松比;支撑辊辊身直径、支撑辊辊身长度、支撑辊密度、支撑辊弹性模量以及支撑辊泊松比等。工作辊、中间辊和支撑辊的材料属性相同,密度为7800 kg/m3,弹性模量为205Gpa,泊松比为0.3。轧辊的尺寸参数如表1所示。
表1轧辊尺寸参数
步骤2.2:采集轧制工艺参数,包括:摩擦系数、轧制速度、前张力、后张力、压下率、工作辊弯辊力、中间辊弯辊力以及中间辊横移量。轧制工艺参数如表2所示。
表2轧制工艺参数
步骤2.3:采集轧制前的带钢参数包括:带钢宽度和入口横向厚度分布数据;将带钢两侧距边部100mm处设为边部减薄区,剩余部分设为中心区,并在带钢横向上以1.5mm为单位间隔取测点,分别测量中心区和两个边部减薄区每个测点的厚度值并连线获取入口带钢横向厚度分布曲线;采集轧制后的带钢参数包括:出口带钢横向厚度分布数据,轧制前后的带钢宽度视为相等。
步骤3:根据步骤1获取的真实应力-应变曲线和步骤2获取的参数建立关于带钢-轧辊变形耦合分析的横向厚度分布仿真模型,所述步骤3具体为:
步骤3.1:根据获取的轧辊数据建立轧机的三维弹塑性有限元模型;
步骤3.2:根据轧前带钢宽度数据和入口带钢横向厚度分布曲线建立带钢的三维弹塑性有限元模型;
步骤3.3:沿着宽度方向将带钢模型分成32个单元,并将步骤1中获取的32条真实应力-应变曲线写入对应位置的单元中,在仿真求解过程中,每个单元的变形过程均遵循对应的真实应力-应变曲线。
步骤4:制定针对不同板形执行机构的实验计划,利用横向厚度分布仿真模型对带钢轧制过程进行模拟实验,所述步骤4具体为:
步骤4.1:根据获取的轧制工艺参数,制定针对于不同执行机构的实验计划,板形执行机构包括:中间辊横移、工作辊弯辊和中间辊弯辊。
具体实施时,将三种执行机构的调控值各设置五个实验值点,分别进行仿真实验,实验设置条件如表3。
表3 工作辊弯辊力、中间辊弯辊力以及中间辊横移量的设定值
步骤4.2:依据对照实验计划,更改横向厚度分布仿真模型中的中间辊横移量、工作辊弯辊力、中间辊弯辊力的取值,进行模拟仿真。
步骤5:构建板形执行机构的厚度调控功效系数计算模型,提取步骤4中各模拟实验稳定轧制阶段的带钢横向厚度分布数据,在不同工作辊弯辊力调节量下的带钢横向厚度分布变化情况如图3,在不同中间辊弯辊力调节量下的带钢横向厚度分布变化情况如图4,在不同中间辊横移量调节量下的带钢横向厚度分布变化情况如图5,计算厚度值和各板形执行机构的厚度调控功效系数,所述步骤5具体为:
步骤5.1:建立以带钢宽度、长度、厚度为坐标轴的三维坐标系,以带钢的中心位置点作为坐标原点;将带钢假设为连续离散化的纵向纤维条,α点设为距带钢中心宽向位置x α 的离散化的纵向纤维条的离散点编号,设板形执行机构编号为i,其调节量为Δu i 。
步骤5.2:通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量:
通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量:
式中,H(x α )为选取的纵向纤维条的轧前横向厚度,h(x α )为选取的纵向纤维条的轧后横向厚度;Δh α,i (x α )为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量,为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量。
步骤5.3:通过带钢轧制前后横向厚度的相对改变量对厚度调控功效系数进行表示为:
式中,T α,i (x α )为第i种板形执行机构在x α 处的厚度调控功效系数。
步骤5.4:提取步骤4中各模拟实验稳定轧制阶段带钢横向厚度分布数据,根据步骤5.2计算获得位置x α 处的带钢轧制前后横向厚度的相对改变量,并根据步骤5.3计算获得各板形执行机构的厚度调控功效系数及其分布曲线。
获取厚度调控功效系数得到三种板形执行机构在不同调节量下的厚度调控功效系数曲线如图6-8所示。中部与边部的板形差值越大,板形执行机构对带钢横截面形状的调控能力就越强,对带钢板凸度以及边部减薄的调节能力也更强。
步骤6:提取步骤4中各模拟实验稳定轧制阶段的带钢宽度数据和出口带钢横向厚度分布曲线,基于出口带钢横向厚度分布曲线的拟合系数和厚度调控功效系数曲线的拟合系数建立带钢横向厚度分布预测计算方程,输入板形执行机构的调控数值以获取对应横向厚度分布曲线,所述步骤6具体为:
步骤6.1:提取步骤4中各模拟实验稳定轧制阶段带钢宽度数据和出口带钢横向厚度分布曲线,对宽度数据进行归一化处理;并对得到的出口带钢横向厚度分布曲线使用五次多项式进行拟合,得到多项式拟合系数B 1、B 2、B 3、B 4、B 5;五次多项式如下:
其中,y为出口带钢横向厚度分布曲线;THK_Center为截距;B 1、B 2、B 3、B 4、B 5分别为出口带钢横向厚度分布曲线的一次、二次、三次、四次和五次项拟合系数;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1]。
步骤6.2:对步骤5.4中获取的三种板形执行机构的厚度调控功效系数曲线使用六次Legendre正交多项式进行拟合,获取三种板形执行机构的厚度调控功效系数曲线的拟合系数A 0、A 1、A 2、A 4、A 6,采用的六次Legendre正交多项式如下式:
其中,E(x)为带钢横向厚度分布曲线;A 0为常数项、A 1、A 2、A 4、A 6分别为厚度调控功效系数曲线的一次、二次、四次和六次项拟合系数,其绝对值大小表明厚度调控功效系数曲线中对一次、二次、四次和六次板形缺陷的调控分量;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1];e(x)为拟合误差。
步骤6.3:根据步骤2.3中划分的边部减薄区和中心区,以归一化的带钢宽度作为x轴、厚度预测值为y轴建立二维坐标系,建立带钢中心区每坐标点厚度预测值的计算方程为:
其中,H P 为计算得到的每坐标点厚度预测值;
其中,V WRB 为工作辊弯辊力,V IRB 为中间辊弯辊力,V IRS 为中间辊横移量,A 0W 、A 0I 和A 0S 都是常数项,A 1W 、A 2W 、A 4W 和A 6W 分别为工作辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1I 、A 2I 、A 4I 和A 6I 分别为中间辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1S 、A 2S 、A 4S 和A 6S 分别为中间辊横移调控功效系数曲线一次、二次、四次、六次项拟合系数;
步骤6.4:由于轧制过程中发生辊系弹性变形和带钢三维变形,在距两边部一定位置处厚度会发生急剧减小现象,即上文构建的带钢中心区的带钢横向厚度分布每坐标点计算方程并不适用于两边部减薄区,需要对带钢两边部减薄区横向厚度构建单独的每坐标点计算方程,建立带钢两个边部减薄区每坐标点厚度预测值的计算方程为:
其中,为带钢每个边部减薄区基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;/>和/>分别是带钢操作侧边部和驱动侧边部基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;H OS 和H DS 别带钢中心区靠近操作侧最边部点厚度值和驱动侧最边部点厚度值;V e是边部减薄值;f t是边部减薄程度权重;f b是位置系数(每个宽度坐标点相对于最边部位置宽度坐标)。
步骤6.5:对步骤6.1中得到的归一化宽度数据进行反归一化处理后分别代入步骤6.3和6.4的计算方程式,计算得到的每坐标点厚度值并获得带钢整体横向厚度分布曲线。即得到带钢横向厚度分布预测模型,工作辊弯辊力、中间辊弯辊力和中间辊横移量在该模型中均作为输入变量。
具体实施时,还可以对步骤6.5中的带钢横向厚度分布预测模型进行精度验证。验证步骤如下:
首先,将步骤2.2中获得的工作辊弯辊力、中间辊弯辊力、中间辊横移量作为变量输入到模型中,获得计算的横向厚度分布曲线;
然后,将预测结果与步骤5.2中获得的横向厚度分布曲线进行对比,验证结果表明模型计算结果与模拟实验仿真值吻合良好,具有较高预测精度,可以为实际生产中快速调整执行机构调节量以达到良好横向厚度分布提供数据支撑。
以上所述仅为本发明的较佳实施例,并不用以限制本发明的思想,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,包括:
步骤1:轧制前截取部分带钢,并对该段带钢横向不同位置处进行拉伸试验获取相应位置的真实应力-应变曲线;
步骤2:获取轧辊参数、轧制工艺参数以及轧制前后的带钢参数;
步骤3:根据步骤1获取的真实应力-应变曲线和步骤2获取的参数建立关于带钢-轧辊变形耦合分析的横向厚度分布仿真模型;
步骤4:制定针对不同板形执行机构的实验计划,利用横向厚度分布仿真模型对带钢轧制过程进行模拟实验;
步骤5:构建板形执行机构的厚度调控功效系数计算模型,提取步骤4中各模拟实验稳定轧制阶段的带钢横向厚度分布数据,计算各板形执行机构的厚度调控功效系数;
步骤6:提取步骤4中各模拟实验稳定轧制阶段的带钢宽度数据和出口带钢横向厚度分布曲线,基于出口带钢横向厚度分布曲线的拟合系数和厚度调控功效系数曲线的拟合系数建立带钢横向厚度分布预测计算方程,输入板形执行机构的调控数值以获取对应横向厚度分布曲线。
2.如权利要求1所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤1具体为:
轧制前截取部分带钢,沿带钢宽度方向将带钢等间距分成N份带钢试样,对每份带钢试样进行拉伸实验,获取所有带钢试样的N条真实应力-应变曲线。
3.如权利要求2所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,获取真实应力-应变曲线具体为:
步骤1.1:定义真实应力-应变曲线由弹性变形阶段、均匀塑性变形阶段和局部塑性变形阶段组成;
步骤1.2:求出弹性变形阶段和均匀塑性变形阶段的分界点,即屈服点对应的屈服应力σs:
其中,P s是材料开始屈服时的载荷,由试验机载荷刻度盘读出;S 0是试样原始横截面面积;
步骤1.3:找出均匀塑性变形阶段各瞬时的流动应力即真实应力Y和对数应变∈:
其中,P为加载瞬间的载荷,由试验机载荷刻度盘读出;S为加载瞬间试样的横截面面积,由体积不变条件求出:
其中,l 0是试样标距的原始长度;Δl是试样标距的伸长量;同样塑性失稳点的应力和应变也由上式求出;
步骤1.4:到达均匀塑性变形阶段和局部塑性变形阶段的分界点,即塑性失稳点后试样开始出现颈缩,找出断裂时的真实应力Y K1及其对应的对数应变∈ K1:
其中,P K1是试样断裂时的载荷;S K1是试样断裂时的横截面面积;l K1是试样断裂时的标距总长度;
步骤1.5:使用下式对颈缩后的局部塑性变形阶段的曲线进行修正:
其中,Y K2是去除形状硬化后的真实应力;d是试件颈缩处直径;ρ是试件颈缩处外形曲率半径。
4.如权利要求2所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤2具体为:
步骤2.1:采集轧辊参数,包括:轧辊直径、轧辊长度、轧辊密度、轧辊弹性模量以及轧辊泊松比;
步骤2.2:采集轧制工艺参数,包括:摩擦系数、轧制速度、前张力、后张力、压下率、工作辊弯辊力、中间辊弯辊力以及中间辊横移量;
步骤2.3:采集轧制前的带钢参数包括:带钢宽度和入口带钢横向厚度分布数据;将带钢两侧距边部100mm处设为边部减薄区,剩余部分设为中心区,并在带钢横向上以1.5mm为单位间隔取测点,分别测量中心区和两个边部减薄区每个测点的厚度值并连线获取入口带钢横向厚度分布曲线;采集轧制后的带钢参数包括:出口带钢横向厚度分布数据,轧制前后的带钢宽度视为相等。
5.如权利要求4所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤3具体为:
步骤3.1:根据获取的轧辊数据建立轧机的三维弹塑性有限元模型;
步骤3.2:根据轧前带钢宽度数据和入口带钢横向厚度分布曲线建立带钢的三维弹塑性有限元模型;
步骤3.3:沿着宽度方向将带钢模型分成N个单元,并将步骤1中获取的N条真实应力-应变曲线写入对应位置的单元中,在仿真求解过程中,每个单元的变形过程均遵循对应的真实应力-应变曲线。
6.如权利要求1所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤4具体为:
步骤4.1:根据获取的轧制工艺参数,制定针对于不同执行机构的实验计划,板形执行机构包括:中间辊横移、工作辊弯辊和中间辊弯辊;
步骤4.2:依据对照实验计划,更改横向厚度分布仿真模型中的中间辊横移量、工作辊弯辊力、中间辊弯辊力的取值,进行模拟仿真。
7.如权利要求4所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤5具体为:
步骤5.1:建立以带钢宽度、长度、厚度为坐标轴的三维坐标系,以带钢的中心位置点作为坐标原点;将带钢假设为连续离散化的纵向纤维条,α点设为距带钢中心宽向位置x α 的离散化的纵向纤维条的离散点编号,设板形执行机构编号为i,其调节量为Δu i ;
步骤5.2:通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量:
通过下式计算在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量:
式中,H(x α )为选取的纵向纤维条的轧前横向厚度,h(x α )为选取的纵向纤维条的轧后横向厚度;Δh α,i (x α )为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的改变量,为在第i种板形执行机构的调节量Δu i 的作用下在x α 处的带钢轧制前后横向厚度的相对改变量;
步骤5.3:通过带钢轧制前后横向厚度的相对改变量对厚度调控功效系数进行表示为:
式中,T α,i (x α )为第i种板形执行机构在x α 处的厚度调控功效系数;
步骤5.4:提取步骤4中各模拟实验稳定轧制阶段带钢横向厚度分布数据,根据步骤5.2计算获得位置x α 处的带钢轧制前后横向厚度的相对改变量,并根据步骤5.3计算获得各执行机构的厚度调控功效系数及其分布曲线。
8.如权利要求7所述的基于带钢横向强度不均的横向厚度分布预测方法,其特征在于,所述步骤6具体为:
步骤6.1:提取步骤4中各模拟实验稳定轧制阶段带钢宽度数据和出口带钢横向厚度分布曲线,对宽度数据进行归一化处理;并对得到的出口带钢横向厚度分布曲线使用五次多项式进行拟合,得到多项式拟合系数B 1、B 2、B 3、B 4、B 5;五次多项式如下:
其中,y为出口带钢横向厚度分布曲线;THK_Center为截距;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1];
步骤6.2:对步骤5.4中获取的三种板形执行机构的厚度调控功效系数曲线使用六次Legendre正交多项式进行拟合,获取三种板形执行机构的厚度调控功效系数曲线的拟合系数A 0、A 1、A 2、A 4、A 6,采用的六次Legendre正交多项式如下式:
其中,E(x)为带钢横向厚度分布曲线;A 0为常数项、A 1、A 2、A 4、A 6分别为厚度调控功效系数曲线的一次、二次、四次和六次项拟合系数,其绝对值大小表明厚度调控功效系数曲线中对一次、二次、四次和六次板形缺陷的调控分量;x为带钢宽度方向归一化后无量纲坐标,x∈[-1,1];e(x)为拟合误差;
步骤6.3:根据步骤2.3中划分的边部减薄区和中心区,以归一化的带钢宽度作为x轴、厚度预测值为y轴建立二维坐标系,建立带钢中心区每坐标点厚度预测值的计算方程为:
其中,H P 为计算得到的每坐标点厚度预测值;
其中,V WRB 为工作辊弯辊力,V IRB 为中间辊弯辊力,V IRS 为中间辊横移量,A 0W 、A 0I 和A 0S 都是常数项,A 1W 、A 2W 、A 4W 和A 6W 分别为工作辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1I 、A 2I 、A 4I 和A 6I 分别为中间辊弯辊调控功效系数曲线一次、二次、四次、六次项拟合系数;A 1S 、A 2S 、A 4S 和A 6S 分别为中间辊横移调控功效系数曲线一次、二次、四次、六次项拟合系数;
步骤6.4: 建立带钢两个边部减薄区每坐标点厚度预测值的计算方程为:
其中,为带钢每个边部减薄区基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;/>和/>分别是带钢操作侧边部和驱动侧边部基于步骤6.1中得到的归一化宽度数据和横向厚度分布曲线的多项式拟合系数计算得到的每坐标点厚度基础值;H OS 和H DS 分别带钢中心区靠近操作侧最边部点厚度值和驱动侧最边部点厚度值;V e是边部减薄值;f t是边部减薄程度权重;f b是位置系数;
步骤6.5:对步骤6.1中得到的归一化宽度数据进行反归一化处理后分别代入步骤6.3和6.4的计算方程式,计算得到的每坐标点厚度值并获得带钢整体横向厚度分布曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310652406.7A CN116371942B (zh) | 2023-06-05 | 2023-06-05 | 一种基于带钢横向强度不均的横向厚度分布预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310652406.7A CN116371942B (zh) | 2023-06-05 | 2023-06-05 | 一种基于带钢横向强度不均的横向厚度分布预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116371942A true CN116371942A (zh) | 2023-07-04 |
CN116371942B CN116371942B (zh) | 2023-10-20 |
Family
ID=86979126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310652406.7A Active CN116371942B (zh) | 2023-06-05 | 2023-06-05 | 一种基于带钢横向强度不均的横向厚度分布预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116371942B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117019884A (zh) * | 2023-10-08 | 2023-11-10 | 东北大学 | 一种冷连轧各机架出口板形预测可视化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104588417A (zh) * | 2015-01-21 | 2015-05-06 | 燕山大学 | 一种基于机理的带钢厚度横向分布的特征参数识别方法 |
CN110947774A (zh) * | 2019-12-06 | 2020-04-03 | 东北大学 | 一种考虑轧制宽展的板形预测方法 |
CN112685687A (zh) * | 2020-11-13 | 2021-04-20 | 西安交通大学 | 一种基于各向异性的单向拉伸试验数据分析方法 |
CN112916624A (zh) * | 2021-01-29 | 2021-06-08 | 华北电力大学(保定) | 一种ucm轧机板形执行机构调控功效系数获取方法 |
CN113275387A (zh) * | 2021-01-29 | 2021-08-20 | 华北电力大学(保定) | 一种ucm轧机辊系横纵刚度特性曲线的获取方法 |
WO2022057925A1 (en) * | 2020-09-18 | 2022-03-24 | Jiangsu Shagang Group Co., Ltd. | Improving surface quality of twin roll cast and hot rolled thin strip steel |
-
2023
- 2023-06-05 CN CN202310652406.7A patent/CN116371942B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104588417A (zh) * | 2015-01-21 | 2015-05-06 | 燕山大学 | 一种基于机理的带钢厚度横向分布的特征参数识别方法 |
CN110947774A (zh) * | 2019-12-06 | 2020-04-03 | 东北大学 | 一种考虑轧制宽展的板形预测方法 |
WO2022057925A1 (en) * | 2020-09-18 | 2022-03-24 | Jiangsu Shagang Group Co., Ltd. | Improving surface quality of twin roll cast and hot rolled thin strip steel |
CN112685687A (zh) * | 2020-11-13 | 2021-04-20 | 西安交通大学 | 一种基于各向异性的单向拉伸试验数据分析方法 |
CN112916624A (zh) * | 2021-01-29 | 2021-06-08 | 华北电力大学(保定) | 一种ucm轧机板形执行机构调控功效系数获取方法 |
CN113275387A (zh) * | 2021-01-29 | 2021-08-20 | 华北电力大学(保定) | 一种ucm轧机辊系横纵刚度特性曲线的获取方法 |
Non-Patent Citations (1)
Title |
---|
王传杰: "《板材成形性能与塑性失稳理论》", 哈尔滨工业大学出版社, pages: 175 - 178 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117019884A (zh) * | 2023-10-08 | 2023-11-10 | 东北大学 | 一种冷连轧各机架出口板形预测可视化方法 |
CN117019884B (zh) * | 2023-10-08 | 2023-12-29 | 东北大学 | 一种冷连轧各机架出口板形预测可视化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116371942B (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110947774B (zh) | 一种考虑轧制宽展的板形预测方法 | |
CN112949108B (zh) | 热轧高强钢板形缺陷全流程预测方法及图形用户界面装置 | |
CN116371942B (zh) | 一种基于带钢横向强度不均的横向厚度分布预测方法 | |
Wang et al. | Numerical and experimental analysis of strip cross-directional control and flatness prediction for UCM cold rolling mill | |
Behrens et al. | Development of an analytical 3D-simulation model of the levelling process | |
Fan et al. | Experimental study on bearing capacity of stainless steel lipped C section stub columns | |
CN112916624B (zh) | 一种ucm轧机板形执行机构调控功效系数获取方法 | |
CN111487869B (zh) | 一种矫直控制方法及装置 | |
Chen et al. | Determination of formability considering wrinkling defect in first-pass conventional spinning with linear roller path | |
CN112926173B (zh) | 一种热轧高强钢板成形极限图的计算方法 | |
Aerens et al. | Springback prediction and elasticity modulus variation | |
Mohammdi Najafabadi et al. | Effect of forming parameters on edge wrinkling in cold roll forming of wide profiles | |
CN113275387B (zh) | 一种ucm轧机辊系横纵刚度特性曲线的获取方法 | |
Petruška et al. | Novel approach to computational simulation of cross roll straightening of bars | |
Shu-zong et al. | Online calculation model of rolling force for cold rolling mill based on numerical integration | |
CN111611699A (zh) | 一种获取板材残余应力分布特征的智能矫直机 | |
CN116329297B (zh) | 一种基于轧件横向力学性能差异的板形预测方法 | |
CN100377802C (zh) | 通过测量棒校准精轧辊装置的装置和方法 | |
CN106769439A (zh) | 一种管线钢热轧卷板屈服强度的测试方法 | |
Laugwitz et al. | Development of levelling strategies for heavy plates via controlled FE models | |
CN107983772A (zh) | 一种板材冷轧过程中可靠性有限元模型的建立方法 | |
CN111842547B (zh) | 一种具有拉压不对称性金属板材矫直工艺新方法 | |
CN117019884B (zh) | 一种冷连轧各机架出口板形预测可视化方法 | |
CN111069363B (zh) | 一种原位纳米增强高强韧钢的弯曲成形工艺实现方法 | |
Meng et al. | Study on a straightening process by reciprocating bending for metal profiles |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |