首页/文章/ 详情

EAAI | 变凸度平整轧制板形控制的一体化建模与模态分层控制策略

3月前浏览713


文章导读

板形是轧制产品质量的关键指标,配备可变凸度(VC)辊的平整机可实现高精度板形控制。本研究建立了融合 VC 辊物理模型与分层多目标优化的新方法,综合塑性变形、轧辊弹性变形及辊系力学行为,采用粒子群算法优化模型参数,并基于 NSGA‑III 算法设计智能分层板形控制策略:先通过弯辊力快速粗调,再结合液压凸度精细调控。经 2503 卷实际生产验证,该方法使平均板形偏差降低 61.4%,标准差降低 56.3%,合格率达 99.4%,为高精度板形控制提供了理论与工程依据。

文章封面 部分截图

     


部分内容 图文解读

如图 1 所示,平整后的带钢若存在边浪或中浪缺陷,将带钢沿宽度方向划分单元并展平后,各纵向单元的长度会出现明显不均,这是板形不良的直接体现,也是本文控制策略所要解决的核心问题。

     

图 1 平整轧制后典型带钢板形缺陷及各单元扁平化长度:(a) 边浪缺陷;(b) 中浪缺陷;(c) 单元划分与扁平化后的边浪缺陷;(d) 单元划分与扁平化后的中浪缺陷。

图 2 展示了 VC 辊的工作原理与平整机布置:VC 辊依靠内部油压使辊套产生弹性变形,实现辊形凸度的在线连续调节,辊套与芯轴采用热装过盈配合以保证密封与扭矩传递;在四辊平整机中,VC 辊作为支承辊使用,配合内外弯辊机构与张力控制系统,通过液压胀形与弯辊力协同作用,实现对带钢板形的高精度调控。

     

图 2  VC 辊单元划分与热装变形示意图:(a) VC 辊单元划分;(b) VC 辊热装变形示意图。初始状态下,芯轴为蓝色虚线,辊套为紫色虚线;热装后配合状态为红色虚线。

图 3 展示了平整机工作辊与支承辊(VC 辊辊套)的辊形及工作辊载荷情况:工作辊受辊间力、轧制力和内外弯辊力作用,轧制时呈横向不对称,引入转角 β 模拟非对称板形,辊套与工作辊接触长度均划分单元且中心单元位置重合;工作辊采用余弦凸度辊形,支承辊中间为平辊形、端部为四次曲线拟合肩台倒角,基于辊间接触力学建立变形协调方程,结合相关公式可求解辊系挠度、压扁系数等关键参数,为板形预测提供支撑。

     

图 3  辊系单元划分与辊形示意图:(a) 辊系单元划分;(b) 支承辊辊形;(c) 工作辊辊形。工作辊的单元划分基于支承辊的单元划分并与之对应。

图 4 展示了板形预测与模型系数优化的完整流程,以构建的板形预测模型为核心,采用粒子群优化(PSO)算法对预测模型中的关键输入参数进行辨识与优化,通过不断缩小模型计算结果与现场实际板形数据的偏差,提升预测精度与鲁棒性。

     

图 4  板形预测与系数优化流程图。板形预测模型为粒子群优化(PSO)的核心,通过减小模型预测值与实际值的偏差实现系数优化。

图 5 展示了模型各输入参数归一化后的数据分布情况,研究采用四分位数法(1.5 倍四分位距)对 6551 组现场生产数据进行异常值剔除,未额外做降噪和平滑处理,以保留真实板形高频特征;从图中可以看出,不同工艺参数的数据密度与聚集区间各不相同,但绝大多数数据都落在合理范围内,直观证明了本研究所用输入数据的分布具有代表性与合理性。

     

图 5  各输入参数的数据分布。

图 6中,NSGA‑Ⅲ 多目标优化中采用分层迭代控制策略:前 100 代以内外弯辊力为主、VC 辊油压固定在 25 MPa 进行快速粗调;100 代后若板形偏差仍大于 3 IU 阈值,则将油压加入决策变量,与弯辊力协同精细优化,充分发挥不同执行器的调节优势,使控制过程更贴合工业现场逻辑。

     

图 6  分层控制策略流程图。

图 7 展示了 VC 辊板形控制系统的一体化应用框架,整个系统分为现场层、数据预处理层、自动化层三级结构,实现信号采集、数据处理、PSO 模型标定与 NSGA‑Ⅲ 分层优化控制,最终通过 PLC 驱动液压执行机构完成闭环控制。

     

图 7  VC 辊板形控制系统整体应用框架。

经 NSGA‑Ⅲ 优化得到Pareto最优解,实现板形分布与板形值双目标平衡,并选出三种典型方案:侧重分布最优、侧重数值最优及综合均衡解。图 8 验证优化后板形明显改善,中浪缺陷得到有效抑制。系统优化耗时仅 185 ms,满足现场实时控制要求。2503 组工业数据表明,相比传统 PID,板形精度、稳定性与合格率均大幅提升,控制效果更优。

     

图 8 Pareto最优解集板形曲线验证:(a) 最优解集与控制前板形对比;(b) 最优解集与目标板形对比。最终解兼顾了板形分布偏差 X₁与板形值偏差 X₂,是最接近目标板形的优化结果。


文章总结

本文针对 VC 辊平整机建立了带钢板形一体化建模与智能控制框架,主要结论如下:


01

构建了融合金属塑性变形、轧辊弹性变形及辊系力学耦合的 VC 辊板形高精度预测模型,可准确反映凸度、弯辊与液压耦合作用。


02

 采用粒子群优化(PSO)算法优化模型关键参数,提升了预测精度与工况适应性。


03

基于 NSGA‑Ⅲ 多目标优化提出智能分层板形控制策略:先由内外弯辊力快速粗调,再由液压凸度精准微调。

 现场验证表明,该方法控制效果优异,弯辊与 VC 辊液压胀形协同作用,显著提升板形均匀性与过程稳定性。该框架为 VC 辊平整机高精度板形控制提供了可行路径,其分层协调思路可推广至 CVC、六辊轧机等多执行器轧制系统,工业应用能有效降低后续加工缺陷与废品率,助力高效智能制造升级。本文不足在于未将轧辊热凸度、磨损等长期时变因素纳入动态模型,可能对辊形基准产生小幅影响。


文章信息

主要作者情况:

Ji Zhang(第一作者),吴仁豪(通讯作者),Hyoung Seop Kim(通讯作者),  Zhenhua Bai(通讯作者)

通讯作者单位:

1.日本东北大学材料高等研究所

2.韩国浦项科技大学金属与环保材料研究所

3.燕山大学国家冷轧板带装备及工艺工程技术研究中心

文章标题:

Integrated modeling and modal layered control strategy for flatness

regulation in variable crown temper rolling

期刊名称:

Engineering Applications of Artificial Intelligence

文章链接:

https://doi.org/10.1016/j.engappai.2026.114202

来源:我的博士日记
理论材料控制PLC
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-05-22
最近编辑:3月前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 130文章 153课程 0
点赞
收藏
作者推荐

多尺度晶体塑性模拟文章推荐

文章名称《Concurrent multi-scale crush simulations with a crystal plasticity model》DOI:10.1016/j.tws.2011.12.019在汽车防撞梁、吸能盒和薄壁管结构中,压溃吸能能力直接影响结构安全性。传统有限元分析通常采用各向同性塑性模型,通过宏观应力–应变曲线描述材料响应。但实际金属材料并不是“均匀黑箱”:晶粒取向、滑移系激活、织构演化都会影响局部塑性变形,尤其在薄壁管压溃这类大变形、强局部化问题中,微观结构可能对吸能行为产生重要影响。压溃模型示意图:Najafi、Marin 和 Rais-Rohani 的文章《Concurrent multi-scale crush simulations with a crystal plasticity model》正是针对这一问题展开研究。作者关注的是:在方形薄壁管轴向压溃过程中,材料初始织构及其演化是否会影响整体压溃力、平均吸能能力和局部折叠模式。为解决这一问题,作者提出了一种并发多尺度建模方法:宏观结构层面采用显式有限元模拟方管压溃;每个积分点内部嵌入一个由多个 FCC 晶粒组成的多晶聚集体;晶粒层面采用 Marin 晶体塑性模型描述滑移、硬化和晶格旋转;最后通过 Taylor 型均匀化获得积分点平均应力。这样,宏观有限元计算不再只依赖经验塑性曲线,而是能够实时考虑晶粒取向和织构演化对结构响应的影响。文章中,作者首先通过单元模型分别施加拉伸、压缩和简单剪切,生成不同初始织构;随后将这些织构赋予方管模型,并进行轴向压溃模拟。结果表明,虽然不同织构对整体折叠形貌的影响并不总是非常显著,但对压溃力–位移曲线、平均压溃力和能量吸收能力具有明显影响。尤其是在角部、水平铰线和锥面等局部大塑性区域,晶粒取向会持续演化,形成不同的局部织构模式。文章还指出,拉伸织构和压缩织构在不同压溃模式下表现出不同的吸能优势,这说明“材料制造历史”并不是可以忽略的背景信息,而是可能影响结构服役性能的重要因素。这篇文章对我们的启发在于:晶体塑性并不只能用于单晶拉伸、RVE 或微观变形分析,也可以嵌入显式动力学框架,用于研究真实工程结构中的局部变形、吸能和织构演化。对于高温合金、铝合金薄壁件、微尺度构件等问题,如果材料存在明显织构或晶粒尺度效应,将晶体塑性与结构有限元耦合,能够提供比传统本构更丰富的物理信息。我们可以将我之前推文提到的umat-taylor模型转化为vumat子程序,进一步使用晶体塑性模型模拟大变形结构尺度材料变形行为。案例展示如下:初始模型参考文章的设置(上下两层钢板,中间为薄壳结构):使用通用接触,摩擦系数设置为0.5,共4000个单元,每个单元包含50个具有不同初始取向晶粒。共20万晶粒。边界条件设置为下端钢板固定,上端下压。模拟结果如下:应力分布结果:晶粒1的剪切滑移:晶粒2的剪切滑移:晶粒50的剪切滑移:单元标号5变形结束后的50个欧拉角分布:来源:我的博士日记

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈