CN107983772B - 一种板材冷轧过程中可靠性有限元模型的建立方法 - Google Patents

一种板材冷轧过程中可靠性有限元模型的建立方法 Download PDF

Info

Publication number
CN107983772B
CN107983772B CN201711209282.6A CN201711209282A CN107983772B CN 107983772 B CN107983772 B CN 107983772B CN 201711209282 A CN201711209282 A CN 201711209282A CN 107983772 B CN107983772 B CN 107983772B
Authority
CN
China
Prior art keywords
value
curve
avg
rolling
cold
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.)
Active
Application number
CN201711209282.6A
Other languages
English (en)
Other versions
CN107983772A (zh
Inventor
杨明绪
许开春
殷福星
何继宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hebei University of Technology
Original Assignee
Hebei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN201711209282.6A priority Critical patent/CN107983772B/zh
Publication of CN107983772A publication Critical patent/CN107983772A/zh
Application granted granted Critical
Publication of CN107983772B publication Critical patent/CN107983772B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21B—ROLLING OF METAL
    • B21B1/00—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations
    • B21B1/22—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations for rolling plates, strips, bands or sheets of indefinite length
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21B—ROLLING OF METAL
    • B21B1/00—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations
    • B21B1/38—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations for rolling sheets of limited length, e.g. folded sheets, superimposed sheets, pack rolling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21B—ROLLING OF METAL
    • B21B1/00—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations
    • B21B1/22—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations for rolling plates, strips, bands or sheets of indefinite length
    • B21B2001/221—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations for rolling plates, strips, bands or sheets of indefinite length by cold-rolling
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21B—ROLLING OF METAL
    • B21B1/00—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations
    • B21B1/38—Metal-rolling methods or mills for making semi-finished products of solid or profiled cross-section; Sequence of operations in milling trains; Layout of rolling-mill plant, e.g. grouping of stands; Succession of passes or of sectional pass alternations for rolling sheets of limited length, e.g. folded sheets, superimposed sheets, pack rolling
    • B21B2001/386—Plates

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Mechanical Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Control Of Metal Rolling (AREA)

Abstract

本发明提供了一种板材冷轧过程中可靠性有限元模型的建立方法，属于有限元仿真模型建模领域。本发明所描述的有限元建模方法基于实际的冷轧生产过程，建立板材轧制过程有限元模型并根据实际冷轧过程的轧制力载不同压下率下进行校对，从而建立有针对性的、可以反映实际冷轧过程的有限元仿真模型，达到有限元仿真在板材冷轧过程中的应用的目的。

Description

一种板材冷轧过程中可靠性有限元模型的建立方法
技术领域
本发明属于有限元仿真模型建模技术领域，尤其是涉及一种板材冷轧过程中可靠性有限元模型的建立方法。
背景技术
冷轧钢板被广泛应用于航空、电力、汽车、建筑、产业设备等各个国民生产的关键领域，属于高附加值的钢铁产品。从轧板产品的发展历史看，冷轧带钢是在热轧带钢的基础上发展起来的，并且有逐渐取代热轧带钢的趋势，这是因为冷轧带钢克服了热轧带钢板厚精度差、表面氧化严重、能源消耗大等诸多缺点，其产品性能好、品种多，用途极为广泛。
近年，随着有限元仿真结果在工业中的应用，有限元仿真结果也在一定程度上对冷轧工艺设计提供了一定的参考，但有限元仿真结果与试验的一致性以及其对试验的指导性还存在很多问题，只能以目前的生产经验为基础结合有限元仿真得到的规律性结果为工艺设计在有限范围内提供一定的参考。
所以，为了利用有限元仿真结果达到减少调试时间、节约生产成本的目的，建立可靠性有限元模拟模型是有限元结果指导生产设计及生产过程分析的前提条件。
发明内容
有鉴于此，本发明旨在提出一种板材冷轧过程中可靠性有限元模型的建立方法，以提高有限元仿真结果的准确性。
为达到上述目的，本发明的技术方案是这样实现的：
一种板材冷轧过程中可靠性有限元模型的建立方法，包括如下步骤：
1）将被加工板材进行冷轧试验，获得至少两种压下率下的轧制力曲线，并确定稳定轧制阶段的压下率曲线，得到试验轧制力平均值曲线，根据测量轧制前后板材厚度计算出压下率R1与R2；
2）利用有限元仿真软件DEFORM3D，建立压下率为R1下的冷轧过程的有限元模型，设被加工板材变形抵抗式为σ=Kε n，作为变形抵抗模型，其中，σ为变形过程中的应力，ε为变形过程中的应变，K为变形抗力，n为加工硬化值，K值、n值采用初始值K0、n0；根据冷轧试验的轧制参数统一有限元轧制参数，运行计算程序；
3）当n= n0时，K值设定为以K0为基准值上下浮动一定范围，设定K值浮动间隔，运行压下率R1对应的有限元模型，得到不同K值时的行程-载荷曲线；当K= K0时，n值设定为以n0为基准值上下浮动一定范围，设定n值浮动间隔，运行压下率R1对应的有限元模型，获得不同n值时的行程-载荷曲线；
4）将不同K值下的行程-载荷曲线取稳定轧制阶段并求取轧制力平均值F avg ，做出平均值F avg 随K值变化曲线，记为F avg-K ，将F avg-K 与冷轧试验获得的平均轧制力拟合求出两条曲线交点对应的K1与n1值；
5）将不同n值时的行程-载荷曲线取稳定轧制阶段并求取轧制力平均值Favg，做出平均轧制力F avg 随n值变化曲线，记为将F avg-n 与F avg 拟合求出两条曲线交点对应的K2与n2值，将F avg-n 与F avg 拟合求出两条曲线交点对应的K2与n2值；
6）将K1与n1值、K2与n2值分别带入变形抵抗式为σ=Kε n，得到σ=K1 ε n1与σ=K2 ε n2，将σ=K1 ε n1与σ=K2 ε n2作为变形抵抗式代入压下率为R2的有限元模型中，其他条件不变，运行程序；将有限元求得平均轧制力F avg-K1 、F avg-K2 与试验平均轧制力曲线F avg-R2 比较，相似结果即为被加工板材在试验用轧制上冷轧过程的变形抵抗式，所求得结果更接近于真实值，以此变形抵抗式来校正有限元模型。
进一步的，所述K值、n值的初始值可由被加工板材的真应力-真应变曲线拟合所得。
相对于现有技术，本发明具有以下优势：
本发明基于实际的冷轧生产过程，建立板材轧制过程有限元模型并根据实际冷轧过程的轧制力载不同压下率下进行校对，从而建立有针对性的、可以反映实际冷轧过程的有限元仿真模型，达到有限元仿真在板材冷轧过程中的应用的目的。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例冷轧过程有限元模型（1/2模型）；
图2为本发明实施例当n=0.33时不同k值下的轧制力曲线与试验轧制力曲线（压下率=55%）。
图3为本发明实施例当n=0.33时K值-轧制力曲线（压下率=55%）。
图4为本发明实施例当K=600时 不同n值下的轧制力曲线（压下率=55%）。
图5为本发明实施例当K=600时n值-轧制力曲线（压下率=55%）。
图6为本发明实施例R=55%时K值-n值-Favg曲线图。
图7为本发明实施例变形抵抗σ=620.5ε 0.33与σ=600ε 0.17有限元仿真所得轧制力与试验结果比较（压下率=40%）。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
本实施例所用被加工板材为Q235，通过单轴拉伸试验获得被加工板材的真应力-应变曲线，并将其拟合得到σ=480ε 0.33；
通过Q235板材冷轧试验获得轧制力曲线，轧制参数为：Q235板材宽度为50mm，厚度2.9mm，轧后厚度为1.3mm，压下率为55%，轧制速度为0.05m/s，轧辊直径为110mm。其中，轧制曲线取稳定轧制阶段用于有限元结果对比。
建立可靠性有限元模拟模型，主要包括如下步骤：
（1）如图1所示，利用有限元仿真软件DEFORM3D建立冷轧过程的有限元模型，板材变形抵抗式为σ=Kε n，其中，K为变形抗力，n为加工硬化值，K，n值以拉伸试验拟合结果为参照分别设定为K=480MPa，n=0.33，并在材料库中输入弹性模量、泊松比等材料物理特性参数；按照冷轧试验的轧制参数统一有限元轧制参数，包括轧辊尺寸、轧制前被加工板材厚度、压下率R为55%、轧辊轧制速度等尺寸参数，摩擦系数按照金属间干摩擦系数设定为0.08，设定冷轧模型的边界条件并运行计算程序；
（2）当n=0.33时，K值分别设定为450-750，数值间隔为50，分别运行程序，得到各K值下的轧制力曲线（如图2所示，图中最长的曲线为试验轧制力曲线），取稳定轧制阶段的平均值Favg，做出K值-轧制力曲线（如图3所示），由图3可知，当n=0.33时，K=600时，所得轧制力曲线与试验相吻合；
（3）同样的，当K=600时，n值分别设定为0.2-0.5，数值间隔为0.05，分别运行程序后，得到不同n值下的轧制力曲线（如图4所示），取稳定轧制阶段的平均值Favg，做出n值-轧制力曲线（如图5所示），并结合步骤（2）所得K值-轧制力曲线，得到K值-n值-Favg曲线图（如图6所示）；
（4）K值-平均轧制力曲线拟合结果为式（2）所示：
K = 0.5972F + 4.7615 （2）
n 值-平均轧制力曲线拟合结果为式（3）所示：
n = -0.0035F + 3.7781 （3）
将冷轧试验获得的平均轧制力F=1031kN分别代入式（2）与（3）所得曲线，并求解得到：K约为620.5，n约为0.17；即当n=0.33时K约为620.5，当K=600时n约为0.17，与实验最为接近；最终得到：σ=620.5ε 0.33 ，σ=600ε 0.17；
（5）将步骤（4）所得解作为变形抵抗式代入有限元模型中，将压下率R调整为40%，将通过有限元仿真所得轧制力与同样压下率40%的冷轧试验平均轧制力做对比，（如图7所示）可知当变形抵抗值为σ=620.5ε 0.33时，有限元所得轧制力与试验吻合程度较高。
本发明方法通过有限元结果与试验结果的多次对比，实现了适用于对比试验用四辊可逆冷轧机有限元模型的建立。
同样的材料在使用不同冷轧机时，均可用此方法进行模型的校正，从而建立具有针对性的可靠性有限元仿真模型。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种板材冷轧过程中可靠性有限元模型的建立方法，其特征在于包括如下步骤：
1）将被加工板材进行冷轧试验，获得至少两种压下率下的轧制力曲线，并确定稳定轧制阶段的压下率曲线，得到试验平均轧制力曲线，根据测量轧制前后板材厚度计算出压下率R1与R2；
2）利用有限元仿真软件DEFORM3D，建立压下率为R1下的冷轧过程的有限元模型，设被加工板材变形抵抗式为σ=Kε n，作为变形抵抗模型，其中，σ为变形过程中的应力，ε为变形过程中的应变，K为变形抗力，n为加工硬化值，K值、n值采用初始值K0、n0；根据冷轧试验的轧制参数统一有限元轧制参数，运行计算程序；
3）当n= n0时，K值设定为以K0为基准值上下浮动一定范围，设定K值浮动间隔，运行压下率R1对应的有限元模型，得到不同K值时的行程-载荷曲线；
4）将不同K值下的行程-载荷曲线取稳定轧制阶段并求取轧制力平均值F avg ，做出平均轧制力F avg 随K值变化曲线，记为F avg-K ，将F avg-K 与冷轧试验获得的平均轧制力拟合求出两条线交点对应的K1与n1值；
5）当K= K0时，n值设定为以n0为基准值上下浮动一定范围，设定n值浮动间隔，运行压下率R1对应的有限元模型，获得不同n值时的行程-载荷曲线；
将不同n值时的行程-载荷曲线取稳定轧制阶段并求取轧制力平均值F avg ，做出平均值F avg 随n值变化曲线，记为将F avg-n ，将F avg-n 与冷轧试验获得的平均轧制力拟合求出两条线交点对应的K2与n2值；
6）将K1与n1值、K2与n2值分别带入变形抵抗式为σ=Kε n，得到σ=K1 ε n1与σ=K2 ε n2，将σ=K1 ε n1与σ=K2 ε n2作为变形抵抗式代入压下率为R2的有限元模型中，其他条件不变，运行程序；将有限元求得平均轧制力F avg-K1 、F avg-K2 与试验平均轧制力曲线F avg-R2 比较，相似结果即为被加工板材在试验用轧制上冷轧过程的变形抵抗式，所求得结果更接近于真实值，以此变形抵抗式来校正有限元模型。
2.根据权利要求1所述的一种板材冷轧过程中可靠性有限元模型的建立方法，其特征在于：所述K值、n值的初始值可由被加工板材的真应力-真应变曲线拟合所得。
CN201711209282.6A 2017-11-27 2017-11-27 一种板材冷轧过程中可靠性有限元模型的建立方法 Active CN107983772B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711209282.6A CN107983772B (zh) 2017-11-27 2017-11-27 一种板材冷轧过程中可靠性有限元模型的建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711209282.6A CN107983772B (zh) 2017-11-27 2017-11-27 一种板材冷轧过程中可靠性有限元模型的建立方法

Publications (2)

Publication Number Publication Date
CN107983772A CN107983772A (zh) 2018-05-04
CN107983772B true CN107983772B (zh) 2019-06-04

Family

ID=62033270

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711209282.6A Active CN107983772B (zh) 2017-11-27 2017-11-27 一种板材冷轧过程中可靠性有限元模型的建立方法

Country Status (1)

Country Link
CN (1) CN107983772B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113569443B (zh) * 2021-06-29 2023-08-22 河钢股份有限公司 一种多机架冷连轧板形有限元仿真方法
CN117350105B (zh) * 2023-09-22 2024-03-22 河北工程大学 一种金属压缩实验考虑鼓肚时数据修正与校核的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3882709A (en) * 1972-10-16 1975-05-13 Nippon Steel Corp Method for controlling the profile of workpieces on rolling mills
CN1190220A (zh) * 1996-10-29 1998-08-12 摩根建社公司 优化轧机轧辊库存的***与方法
EP1491268A2 (de) * 2003-06-25 2004-12-29 ABB PATENT GmbH Verfahren zum Kaltwalzen metallischen Bandes
CN104376424A (zh) * 2014-11-27 2015-02-25 东北大学 一种钢铁企业冷轧区多产线钢卷协调调度方法
CN106964654A (zh) * 2016-11-07 2017-07-21 南阳师范学院 一种求解轧制力的建模和仿真方法
CN107066737A (zh) * 2017-04-14 2017-08-18 北京科技大学 一种预测热轧过程板带温度场的二维交替差分方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3882709A (en) * 1972-10-16 1975-05-13 Nippon Steel Corp Method for controlling the profile of workpieces on rolling mills
CN1190220A (zh) * 1996-10-29 1998-08-12 摩根建社公司 优化轧机轧辊库存的***与方法
EP1491268A2 (de) * 2003-06-25 2004-12-29 ABB PATENT GmbH Verfahren zum Kaltwalzen metallischen Bandes
CN104376424A (zh) * 2014-11-27 2015-02-25 东北大学 一种钢铁企业冷轧区多产线钢卷协调调度方法
CN106964654A (zh) * 2016-11-07 2017-07-21 南阳师范学院 一种求解轧制力的建模和仿真方法
CN107066737A (zh) * 2017-04-14 2017-08-18 北京科技大学 一种预测热轧过程板带温度场的二维交替差分方法

Also Published As

Publication number Publication date
CN107983772A (zh) 2018-05-04

Similar Documents

Publication Publication Date Title
Qian et al. 3D coupled macro–microscopic finite element modelling and simulation for combined blank-forging and rolling process of alloy steel large ring
CN107983772B (zh) 一种板材冷轧过程中可靠性有限元模型的建立方法
CN102539315B (zh) 一种快速准确确定金属塑性成形过程摩擦系数的方法
CN101692028B (zh) 金属板料大变形流动应力曲线的测定方法
CN104561848B (zh) 一种蠕变时效成形工艺方法
CN103341503B (zh) 一种自适应凸度变化热轧板形控制模型
CN103344478B (zh) 一种测定薄板反向加载包辛格效应的方法
CN109351785A (zh) 一种轧制力优化方法及装置
CN105243193B (zh) 一种确定徐变试验棱柱体试件抗压强度折算系数的方法
Jiao-Jiao et al. A novel approach to springback control of high-strength steel in cold roll forming
CN112926173B (zh) 一种热轧高强钢板成形极限图的计算方法
CN106802256A (zh) 一种构建金属薄板成形极限应力图的方法
CN101417292B (zh) 依托常规轧制模型控制中低牌号电工钢轧制的方法
Hussain et al. The effect of variation in the curvature of part on the formability in incremental forming: An experimental investigation
CN111950098A (zh) 基于有限元仿真的金属板料冲孔工艺参数优化方法
Morris et al. Some fundamental considerations for the control of residual flatness in tension levelling
CN116136892B (zh) 一种二十辊轧机轧制力的计算方法及***
Li et al. Low cycle fatigue and ratcheting properties of steel 40Cr under stress controlled tests
CN113275387B (zh) 一种ucm轧机辊系横纵刚度特性曲线的获取方法
CN104568611B (zh) 基于dic应变测量***的板料成形能力与变形均化能力评价方法
Chen et al. Research on the constitutive model of cold-formed thick-walled steel under cyclic loading
CN110303070B (zh) 一种双金属复合板矫直工艺新方法
CN106599446A (zh) 一种板壳材料起皱失稳极限图的建立方法及***
CN113204870B (zh) 一种现场原级配堆石料力学参数推测方法
KR100467231B1 (ko) 냉간 조질압연에서의 압연하중 예측방법

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