[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN106971026A - 基于全局可加模型的微合金钢力学性能预报方法 - Google Patents

基于全局可加模型的微合金钢力学性能预报方法 Download PDF

Info

Publication number
CN106971026A
CN106971026A CN201710135462.8A CN201710135462A CN106971026A CN 106971026 A CN106971026 A CN 106971026A CN 201710135462 A CN201710135462 A CN 201710135462A CN 106971026 A CN106971026 A CN 106971026A
Authority
CN
China
Prior art keywords
mechanical properties
alloyed steel
micro alloyed
properties forecast
model
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
Application number
CN201710135462.8A
Other languages
English (en)
Other versions
CN106971026B (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.)
Wuhan University of Science and Engineering WUSE
Wuhan University of Science and Technology WHUST
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN201710135462.8A priority Critical patent/CN106971026B/zh
Publication of CN106971026A publication Critical patent/CN106971026A/zh
Priority to US15/856,030 priority patent/US11475317B2/en
Application granted granted Critical
Publication of CN106971026B publication Critical patent/CN106971026B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00Computing arrangements using knowledge-based models
    • G06N5/02Knowledge representation; Symbolic representation
    • G06N5/022Knowledge engineering; Knowledge acquisition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B1/00Metal-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/22Metal-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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B21MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21BROLLING OF METAL
    • B21B1/00Metal-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/22Metal-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/225Metal-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 hot-rolling
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D7/00Modifying the physical properties of iron or steel by deformation
    • C21D7/13Modifying the physical properties of iron or steel by deformation by hot working
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/24Sheet material
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00Computing arrangements using knowledge-based models
    • G06N5/01Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Control Of Metal Rolling (AREA)
  • Heat Treatment Of Steel (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)

Abstract

本发明提供基于全局可加模型的微合金钢力学性能预报方法,包括以下步骤:确定力学性能预报模型的影响因素;计算微合金钢轧制过程中碳氮化物析出的成分及含量;根据广义可加模型理论,将微合金钢力学性能预报模型表达为若干个子模型的可加形式;估计微合金钢力学性能预报模型;以及验证各子模型的可靠性。上述方案获得的微合金钢产品性能预报模型具有计算精度高、适应范围广等优点,可用于新产品设计和钢种成分优化,达到减少物理试验次数,缩短产品研发周期,降低成本的目的。

Description

基于全局可加模型的微合金钢力学性能预报方法
技术领域
本发明涉及冶金工业轧制过程数学模型,特别涉及基于全局可加模型的微合金钢力学性能预报方法。
背景技术
为了满足社会需求和适应市场竞争,现代钢铁企业必须不断改善钢铁产品的组织性能,同时降低生产成本,缩短产品研发周期。这就需要改变以往昂贵而耗时的物理试验方法,借助先进的计算机辅助支撑工具来完成钢铁产品的设计。热轧带钢力学性能预报模型可用于热轧新产品设计和钢种优化,因此一直都是世界各国冶金工作者研究方向之一,与之相关的大量研究工作也一直在开展。
热轧带钢力学性能预报模型建立的理论依据是产品的力学性能决定于其内部组织,而组织又取决于带钢成分及其生产工艺,这两者与力学性能间存在确定而复杂的关系。力学性能预报模型就是要揭示这些参数之间的定量关系,从而根据这些参数来预报产品的力学性能。性能预报问题经过多年的研究,已取得了大量研究成果,国内外众多钢铁企业都试用或使用过性能预报系统,但现有模型在可靠性、实用性与外延性方面还存在不足,难以满足现场用户使用的要求。
从建模思路来看,性能预报问题建模思路大致分为两种:一种是机理建模,基于实验室物理冶金学的实验结果,建立数学模型来预测产品的力 学性能;另一种是数据建模,根据实际生产数据,利用神经网络等智能算法来预测力学性能。以上两种建模方式均存在不足,机理建模多以钢种为单位进行研究,存在不能覆盖众多钢种的局限性,也不便研究成分、工艺对组织性能的影响规律,计算精度不够;而数据建模的缺陷在于:建立的模型存在不确定性,无法保证外推能力(外延性),对训练样本精度高,但难以适应未经过训练的数据样本。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供一种可靠性高、外延性强、实用性广的基于全局可加模型的微合金钢力学性能预报方法。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,具有这样的特征,包括以下步骤:
步骤A,确定微合金钢力学性能预报模型的影响因素;
步骤B,计算微合金钢轧制过程中析出的碳氮化物成分及含量;
步骤C,根据广义可加模型理论,将微合金钢力学性能预报模型表达为若干个子模型的可加形式;
步骤D,估计微合金钢力学性能预报模型;以及
步骤E,验证子模型的可靠性。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,步骤A,具有以下步骤:
步骤A-1,根据冶金机理、数据采集的过程和先验知识将影响因素分为四大类:化学成分、工艺参数、检测参数、异常标记;以及
步骤A-2,分别在每个大类中寻找影响因素,
其中,步骤A-2,包含以下步骤:
步骤A-2-1,将一些公认的影响因素筛选出来;
步骤A-2-2,寻找可能会起显著作用的影响因素。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,步骤B,具有以下步骤:
步骤B-1,基于多元第二相固溶析出理论,开发能定量计算热轧过程中碳氮化物析出成分及含量的热力学模型,并考虑碳氮化物和AlN的相互作用;
步骤B-2,对热力学模型进行求解,采用牛顿拉夫森算法进行数值求解,通过迭代求解过程,确定不同成分下的碳氮化物和AlN的析出开始温度与析出顺序,最终获得微合金钢轧制过程中析出的碳氮化物成分及含量。
基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,步骤C,具有以下步骤:
微合金钢力学性能预报模型的基本形式为
其中,TS代表钢材的力学性能参数;α为截距;Sj(·)为各自变量Xj的三次光滑样条函数,同时也代表力学性能预报模型的各个子模型;p为自变量个数。
基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,步骤D,钢材力学性能预报模型的基本形式确定后,根据从现场收集的热轧过程数据,对微合金钢力学性能预报模型进行估计,
对三次光滑样条函数的估计,采用Back-fitting算法,其具体实现过程如下:
初始化,令m=0:
迭代过程,令m=m+1:
终止:直到小于预定正值常数或循环次数达到指定次数。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,步骤E,当力学性能预报模型建好后,需要根据机理知识与生产数据来验证各子模型的可靠性。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,热力学模型为:
[Ala][Na]=KAlN
KMC、KMN分别为纯碳化物和氮化物的固溶度积;[Ma]、[Ca]和[Na]为固溶的各元素的摩尔分数;LCN=-4260J/mol,交互作用系数;R=8.314J/(K·mol),气体常数;T为绝对温度;M0、C0、N0为钢中微合金元素、碳元素及氮元素的摩尔分数;Al0为钢中Al的摩尔分数;fp为析出的碳氮化铌的摩尔分数;[Ala]为奥氏体中固溶的Al的摩尔分数;KAIN为AlN的固溶度积;fAlN为析出的AlN的摩尔分数。此外,KMC、KMN和KAlN均为二元第二相的平衡固溶度积,表示为:
KMX=[Ma][Xa]=C×10A-B/T
式中:AFe、AM和AX分别为铁、微合金元素和间隙元素(碳元素或氮元素)的相对原子质量;A、B为固溶度积公式的系数。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,子模型包括两类:成分类子模型和工艺类子模型。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,间隙元素为碳元素。
本发明提供的基于全局可加模型的微合金钢力学性能预报方法,还具有这样的特征:其中,间隙元素为氮元素。
发明作用和效果
根据本发明所涉及基于全局可加模型的微合金钢力学性能预报方法, “全局”是指通过收集多个钢种的生产数据来建模,而不仅仅针对某个具体钢种;“可加”是指模型在形式上是可加的,采用广义可加模型的思路,将力学性能预报模型分解为多个子模型。
另外,需要考虑微合金元素碳氮化物析出对力学性能的影响,基于多元第二相固溶析出理论开发热力学模型,对不同温度、成分下的微合金钢轧制过程中碳氮化物析出成分及含量进行定量计算,将计算获得的碳氮化物析出成分、含量与微合金钢实际生产数据结合起来,建立可靠性高、外延性强、实用性广的力学性能预报模型。
上述技术方案获得的微合金钢产品力学性能预报模型具有计算精度高、适应范围广等优点,可用于微合金钢新产品设计和钢种优化,达到减少传统的物理试验次数,缩短产品研发周期,降低研发成本的目的。
附图说明
图1是本发明在实施例中的基于全局可加模型的微合金钢力学性能预报方法的流程图;
图2是本发明在实施例中的步骤A的流程图;
图3是本发明在实施例中的步骤A-2的流程图;
图4是本发明在实施例中的确定变量重要性的流程图;
图5是本发明在实施例中的微合金钢碳氮化合物析出的热力学数值计算流程;
图6(a)是本发明在实施例中的无AlN析出的计算结果图;
图6(b)是本发明在实施例中的有AlN析出的计算结果图;
图7(a)是本发明在实施例中的影响出抗拉强度TS的出炉温度的光滑函数图;
图7(b)是本发明在实施例中的影响出抗拉强度TS的卷取温度的光滑函数图;
图7(c)是本发明在实施例中的影响出抗拉强度TS的终轧厚度的光滑函数图;
图7(d)是本发明在实施例中的影响出抗拉强度TS的剩余碳(质量分数)的光滑函数图;
图7(e)是本发明在实施例中的影响出抗拉强度TS的硅(质量分数)的光滑函数图;
图7(f)是本发明在实施例中的影响出抗拉强度TS的锰(质量分数)的光滑函数图;
图7(g)是本发明在实施例中的影响出抗拉强度TS的碳化铌(质量分数)的光滑函数图;
图7(h)是本发明在实施例中的影响出抗拉强度TS的氮化铌(质量分数)的光滑函数图;
图8(a)是本发明在实施例中的力学性能预报模型的抗拉强度预测值与实测值的对比图;
图8(b)是本发明在实施例中的力学性能预报模型的屈服强度预测值与实测值的对比图。
具体实施方式
以下参照附图及实施例对本发明所涉及的基于全局可加模型的微合金钢力学性能预报方法作详细的描述。
实施例
下面以热连轧含Nb微合金钢为例进行力学性能预报建模。
如图1、图2和图3所示,步骤A:确定微合金钢力学性能预报模型的影响因素,进入步骤B。
对力学性能预报过程参数进行工业大数据分析,寻找合适的影响因素。由于影响因素具有多样性、复杂性,因此影响因素的寻找过程是一个反复校正的过程。通过随机森林、贝叶斯方法等分析方法,同时结合冶金学机理及先验知识,选取带钢中碳、硅、锰、磷、铌等化学成分,以及出炉温度TF、粗轧温度TR、精轧入口温度TFE、终轧温度TFD、卷取温度TC、终轧厚度HFD等轧制工艺参数作为模型的影响因素(即自变量);抗拉强度、屈服强度等力学性能参数作为模型的响应变量。
影响力学性能的因素很多,而且部分变量之间存在交互作用。这就使得寻找影响因素时需要将冶金机理、统计方法与先验知识相结合。
步骤A-1:将影响因素分类:根据冶金机理、数据采集的过程和先验知识将影响因素分为四大类:化学成分、工艺参数、检测参数、异常标记。与检测相关的变量就是用于描述“数据是怎么来的”。如:取样机组、取样时间、取样方向、试样规格、屈服类型等等。当模型研究到一定的深度后,这些因素会成为限制模型精度和可靠度提升的主要因素。进入步骤A-2。
步骤A-2:分别在每个大类中寻找影响因素,具有以下步骤:
步骤A-2-1:将一些公认的影响因素筛选出来。例如:C、Si、Mn、N、Nb、V、Ti等成分,出炉温度、终轧温度、卷取温度、终轧厚度等工艺参数。然后,寻找可能会起显著作用的影响因素。图4为用随机森林判断一个变量重要性的流程图,根据变量的重要性来判断该变量是否起作用。值得注意的是:随机森林寻找的变量只是表示变量对数据具有较好的区分性,但是不能说明该变量对因变量一定有影响作用。由于数据的复杂性,寻找影响因素的过程是一个反复假设、分析、检验与修正的过程。
步骤B:计算微合金钢轧制过程中析出的碳氮化物成分及含量,进入步骤C。
步骤B-1:基于多元第二相固溶析出理论,开发能定量计算热轧过程中碳氮化物析出成分及含量的热力学模型,并考虑碳氮化物和AlN的相互作用。对于Fe-M-C-N合金体系的热力学模型如下:
上述式中:KMC、KMN分别为纯碳化物和氮化物的固溶度积;[Ma]、[Ca]和[Na]为固溶的各元素的摩尔分数;LCN=-4260J/mol,交互作用系数;R=8.314J/(K·mol),气体常数;T为绝对温度;M0、C0、N0为钢中微合金元素、碳元素及氮元素的摩尔分数;fp为析出的碳氮化铌的摩尔分数。
对于存在AlN析出的情况,由于其具有六方晶格结构,不溶于具有NaCl结构的碳氮化合物而独立存在。对于Fe-Al-Nb-C-N热力学平衡系统,需要将式(3)~(5)替换为下列式子:
[Ala][Na]=KAlN (6)
上述式中:[Ala]为奥氏体中固溶的Al的摩尔分数;KAIN为AlN的固溶度积;fAlN为析出的AlN的摩尔分数;Al0为钢中Al的摩尔分数。此外,KMC、KMN和KAlN均为二元第二相的平衡固溶度积,一般可统一表示为:
KMX=[Ma][Xa]=C×10A-B/T (11)
式中:AFe、AM和AX分别为铁、合金元素和间隙元素(碳元素或氮元素)的相对原子质量;A、B为固溶度积公式的系数。
步骤B-2:对上述热力学模型进行求解。由于上述热力学模型是一组非常复杂的非线性方程组,采用牛顿拉夫森算法进行数值求解,通过迭代求解过程,确定不同成分下的碳氮化物和AlN的析出开始温度与析出顺序,最终获得微合金钢轧制过程中析出的碳氮化物成分及含量,计算流程如图5所示。
研究表明带钢中Nb和C、N会在奥氏体相中析出形成碳氮化铌复合化 合物,所以需要采用热力学模型计算出各卷带钢在轧制过程中析出的NbC、NbN含量,用质量分数来表示。
为了示例,图6(a)和图6(b)为其中2卷带钢的计算结果,图6(a)为无AlN析出情况,图6(b)为有AlN析出情况。
模型可自动根据带钢成分、温度计算是否存在AlN的析出。为便于对比,各元素固溶质量分数采用相对于原始成分的比值来表示。
由图6(b)可见,随着温度的降低,固溶的N、Nb元素的质量分数在不断减少,而析出的NbC、NbN的质量分数在不断增加;整个过程中,碳元素减少的速度很慢。通过上述计算,可以得到各卷带钢在奥氏体中析出的NbC、NbN的质量分数,再用于后续的力学性能预报建模研究。
步骤C:根据广义可加模型理论,将微合金钢力学性能预报模型表达为若干个子模型的可加形式,进入步骤D。
首先,简要介绍广义可加模型相关理论。广义可加模型是广义线性模型的非参数化拓展,通过对自变量用非参数函数形式来拟合估计响应变量与自变量间的关系,模型中每一个加性项使用单个光滑函数估计,在每一加性项中可以解释响应变量如何随自变量的变化而变化。其优点是可以处理高维数据中因变量与自变量之间的非线性关系,适于对数据进行探索性分析或寻找响应变量与自变量间是否存在依存关系。广义可加模型因其结构简单、灵活性与稳定性好等优点在许多学科领域得到了广泛应用,其数学表达可写为:
式中:μ为响应变量Y的条件期望,μ=E(Y|X1,X2,...,XP);g(·)为连接函数;α为截距;fj(·)为自变量Xj的任意单变量函数;p为自变量个数。
对于钢材力学性能预报模型而言,连接函数可选择单位连接函数;依据自变量与响应变量之间的散点图,选取各自变量的单变量光滑函数的估计方法。鉴于三次光滑样条具有计算简单、稳定性好、收敛性有保证等优点,选择三次光滑样条估计方法。
综上分析,钢材力学性能预报模型的基本形式可写为:
其中,TS代表钢材的抗拉强度、屈服强度、延伸率等力学性能参数;α为截距;Sj(·)为各自变量Xj的三次光滑样条函数,同时也代表力学性能预报模型的各个子模型;p为自变量个数。
子模型大体包括两类:C、Mn、Si、P、Nb、Ti、V、Mo、Cr及相互作用产生的碳氮化合物等成分子模型;出炉温度、粗轧温度、精轧温度、卷取温度、厚度、规格等工艺子模型。
采用三次光滑样条来估计各自变量的单变量光滑函数,设置模型形式为:
式中:Si(Xi)为各自变量的三次光滑样条函数,主要包括TF、TR、TFE、TFD、TC、HFD、Cs(剩余碳),以及Si、Mn、P、NbC、NbN等影响因素。
步骤D:估计微合金钢力学性能预报模型,进入步骤E。
钢材力学性能预报模型的基本形式确定后,根据从现场收集的大量热 轧过程数据,对微合金钢力学性能预报模型进行估计,也就是通过估计算法来获得各个子模型的参数形式或非参数形式。
对各单变量光滑函数Sj(Xj)的估计,采用Back-fitting算法,其具体实现过程如下:
初始化,令m=0:
迭代过程,令m=m+1:
终止:直到小于预定正值常数或循环次数达到指定次数。
设置Back-Fitting算法的循环次数为10,循环迭代计算结束时,求得α=539.16。下面给出抗拉强度TS主要影响因素的光滑函数图,如图7(a)、7(b)、7(c)、7(d)、7(e)、7(f)、7(g)、7(h)所示,它们代表了各因素对抗拉强度TS的影响规律。
步骤E:验证子模型的可靠性。
当力学性能预报模型建好后,需要根据冶金机理与生产数据来验证各子模型的可靠性,即验证用它获得的规律是否真实。真实的规律往往都具有可重复性。因此,在验证子模型的可靠性时,要紧紧抓住规律的可重复性,在不同的维度下验证它是否可以被数据重现,从而避免一些虚假规律的产生。如果用一个子模型获得的规律在大多数情况下都能够被数据重现或者找不到明显的证据证明这个规律是虚假的,那么就认为这个子模型是 可靠的。反之,需要分析子模型失效的原因,重新修正子模型。
如图7(a)、7(b)、7(c)、7(d)、7(e)、7(f)、7(g)、7 (h)所示:(1)抗拉强度随出炉温度的升高而增大,这是因为加热温度 越高,固溶进奥氏体中的铌质量分数越多,后续降温过程中能析出的碳氮 化物也越多;(2)抗拉强度随卷取温度的升高而降低,这是因为TC越高, 则卷取时的相变温度和析出温度就越高,相变温度增加会减小相变驱动力, 使组织粗化,从而降低强度;同时,析出物在较高的温度下更容易发生长 大,降低析出强化效果;(3)终轧厚度的增加会导致抗拉强度降低,这是 因为带钢厚度越薄,则在层流冷却过程中冷得越快,因而过冷奥氏体的转 变温度降低,驱动力增加,相变后的组织变细;(4)剩余碳质量分数越大, 则抗拉强度越大,这是因为碳元素有很强的固溶强化效果,含碳量增加使 固溶强化作用增加;另外,碳元素明显增强奥氏体的稳定性,使过冷奥氏 体转变温度降低,从而细化组织;(5)与碳元素类似,硅和锰元素也有固 溶强化作用,同时也能增加奥氏体稳定性;(6)氮化铌、碳化铌析出物的 质量分数越大,则抗拉强度越大,这是因为碳氮析出物从变形奥氏体的晶 界处析出,钉扎晶界,阻碍奥氏体再结晶,从而保持了奥氏体的形变效果,为随后的铁素体相变提供更多的形核点,细化铁素体晶粒,起到细晶强化 的作用。此外,图5也表明,碳、硅、锰、碳氮化铌等成分对抗拉强度的影 响程度大于出炉温度、终轧厚度等轧制工艺参数的影响。
根据式(14)及图7(a)、7(b)、7(c)、7(d)、7(e)、7(f)、 7(g)、7(h)中求得的各单变量光滑函数,可得到热轧含Nb微合金钢抗 拉强度的预报模型,其中成分、工艺各因素对应的函数值可根据图8中的样 条函数插值得到。同样地,应用上述方法可得到屈服强度的预报模型。
采用上述模型,计算7740卷某热连轧机组生产的含Nb微合金钢钢(含多钢种)的抗拉强度TS与屈服强度YS,模型预测值与实测值的对比如图8 所示,表1为模型预报的均方根误差和平均误差ME,表2为其中几卷带钢抗拉强度的中间计算结果。
表1模型预报的均方根误差和平均误差ME
表2其中几卷带钢抗拉强度的中间计算结果
可见,得到的子模型均是符合已有经验和机理知识的,且新模型具有较高的预报精度;此外,计算过程无需人工修正,模型预测值偏离实测值均在一定范围内,模型具有较强的适应能力。
实施例的作用与效果
根据本实施例所涉及基于全局可加模型的微合金钢力学性能预报方法,“全局”是指通过收集多个钢种的生产数据来建模,而不仅仅针对某个具体钢种;“可加”是指模型在形式上是可加的,采用广义可加模型的思路,将力学性能预报模型分解为多个子模型。
另外,需要考虑微合金元素碳氮化物析出对力学性能的影响,基于多元第二相固溶析出理论开发热力学模型,对不同温度、成分下的微合金钢轧制过程中碳氮化物析出成分及含量进行定量计算,将计算获得的碳氮化物析出成分、含量与微合金钢实际生产数据结合起来,建立可靠性高、外延性强、实用性广的力学性能预报模型。
上述技术方案获得的微合金钢产品力学性能预报模型具有计算精度高、适应范围广等优点,可用于微合金钢新产品设计和钢种优化,达到减少传统的物理试验次数,缩短产品研发周期,降低研发成本的目的。
上述实施方式为本发明的优选案例,并不用来限制本发明的保护范围。

Claims (10)

1.一种基于全局可加模型的微合金钢力学性能预报方法,其特征在于,包括以下步骤:
步骤A,确定微合金钢力学性能预报模型的影响因素;
步骤B,计算微合金钢轧制过程中析出的碳氮化物成分及含量;
步骤C,根据广义可加模型理论,将所述微合金钢力学性能预报模型表达为若干个子模型的可加形式;
步骤D,估计所述微合金钢力学性能预报模型;以及
步骤E,验证所述子模型的可靠性。
2.根据权利要求1所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述步骤A,具有:
步骤A-1,根据冶金机理、数据采集的过程和先验知识将影响因素分为四大类:化学成分、工艺参数、检测参数、异常标记;以及
步骤A-2,分别在每个大类中寻找所述影响因素,
其中,所述步骤A-2,包含以下步骤:
步骤A-2-1,将一些公认的影响因素筛选出来;
步骤A-2-2,寻找可能会起显著作用的所述影响因素。
3.根据权利要求1所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述步骤B,具有以下步骤:
步骤B-1,基于多元第二相固溶析出理论,开发能定量计算热轧过程中碳氮化物析出成分及含量的热力学模型,并考虑碳氮化物和AlN的相互作用;
步骤B-2,对所述热力学模型进行求解,采用牛顿拉夫森算法进行数值求解,通过迭代求解过程,确定不同成分下的所述碳氮化物和所述AlN的析出开始温度与析出顺序,最终获得微合金钢轧制过程中析出的碳氮化物成分及含量。
4.根据权利要求1所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述步骤C,具有以下步骤:
所述微合金钢力学性能预报模型的基本形式为
T S = α + Σ j = 1 p S j ( X j )
其中,TS代表钢材的力学性能参数;α为截距;Sj(·)为各自变量Xj的三次光滑样条函数,同时也代表力学性能预报模型的各个子模型;p为自变量个数。
5.根据权利要求1所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述步骤D,所述钢材力学性能预报模型的基本形式确定后,根据从现场收集的热轧过程数据,对所述微合金钢力学性能预报模型进行估计,
对三次光滑样条函数的估计,采用Back-fitting算法,其具体实现过程如下:
初始化,令m=0:
α 0 = E ( Y ) = Σ i = 1 n y i n
f 1 m ( X 1 ) = f 2 m ( X 2 ) = ... = f p m ( X p ) = 0
迭代过程,令m=m+1:
R i = Y - α 0 - Σ j = 1 i - 1 f j m ( X j ) - Σ j = i + 1 p f j m - 1 ( X j )
f i m ( X i ) = E ( R i | X i )
终止:直到小于预定正值常数或循环次数达到指定次数。
6.根据权利要求1所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述步骤E,当一个所述力学性能预报模型建好后,需要根据机理知识与生产数据来验证该所述子模型的可靠性。
7.根据权利要求3所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述热力学模型为:
l n yK M C [ M a ] [ C a ] + ( 1 - y ) 2 L C N R T = 0
l n ( 1 - y ) K M N [ M a ] [ N a ] + y 2 L C N R T = 0
[Ala][Na]=KAlN
M 0 = 1 2 f p + ( 1 - f p - f A l N ) [ M a ]
Al 0 = 1 2 f A l N + ( 1 - f p - f A l N ) [ Al a ]
C 0 = y 2 f p + ( 1 - f p - f A l N ) [ C a ]
N 0 = ( 1 - y ) 2 f p + 1 2 f A l N + ( 1 - f p - f A l N ) [ N a ]
KMC、KMN分别为纯碳化物和氮化物的固溶度积;[Ma]、[Ca]和[Na]为固溶的各元素的摩尔分数;LCN=-4260J/mol,交互作用系数;R=8.314J/(K·mol),气体常数;T为绝对温度;M0、C0、N0为钢中微合金元素、碳元素及氮元素的摩尔分数;fp为析出的碳氮化铌的摩尔分数;[Ala]为奥氏体中固溶的Al的摩尔分数;KAIN为所述AlN的固溶度积;fAlN为析出的所述AlN的摩尔分数;Al0为钢中Al的摩尔分数。此外,KMC、KMN和KAlN均为二元第二相的平衡固溶度积,表示为:
KMX=[Ma][Xa]=C×10A-B/T
式中:AFe、AM和AX分别为铁、合金元素和间隙元素的相对原子质量;A、B为固溶度积公式的系数。
8.根据权利要求4所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述子模型包括两类:成分类子模型和工艺类子模型。
9.根据权利要求7所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述间隙元素为碳元素。
10.根据权利要求7所述的基于全局可加模型的微合金钢力学性能预报方法,其特征在于:
其中,所述间隙元素为氮元素。
CN201710135462.8A 2017-03-08 2017-03-08 基于全局可加模型的微合金钢力学性能预报方法 Active CN106971026B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201710135462.8A CN106971026B (zh) 2017-03-08 2017-03-08 基于全局可加模型的微合金钢力学性能预报方法
US15/856,030 US11475317B2 (en) 2017-03-08 2017-12-27 Microalloyed steel mechanical property prediction method based on globally additive model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710135462.8A CN106971026B (zh) 2017-03-08 2017-03-08 基于全局可加模型的微合金钢力学性能预报方法

Publications (2)

Publication Number Publication Date
CN106971026A true CN106971026A (zh) 2017-07-21
CN106971026B CN106971026B (zh) 2020-06-16

Family

ID=59329341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710135462.8A Active CN106971026B (zh) 2017-03-08 2017-03-08 基于全局可加模型的微合金钢力学性能预报方法

Country Status (2)

Country Link
US (1) US11475317B2 (zh)
CN (1) CN106971026B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108182340A (zh) * 2017-12-11 2018-06-19 中国科学院合肥物质科学研究院 一种高固溶性多主元合金的快速筛选方法
CN108428038A (zh) * 2018-01-23 2018-08-21 安徽大学 基于gam模型的土壤锌含量影响因素分析方法
CN109145493A (zh) * 2018-09-11 2019-01-04 张家港宏昌钢板有限公司 开发螺纹钢新品的方法
CN109902872A (zh) * 2019-02-28 2019-06-18 北京首钢股份有限公司 一种性能预测模型与产品研发相结合的方法和装置
CN110797096A (zh) * 2019-10-30 2020-02-14 合肥通用机械研究院有限公司 一种基于深度学习的耐热合金高温力学性能预测方法
CN111931425A (zh) * 2020-09-15 2020-11-13 山东钢铁集团日照有限公司 一种基于rf算法的热轧q355b钢种的性能预测系统
CN117110871A (zh) * 2023-10-13 2023-11-24 北京中航科电测控技术股份有限公司 高功率密度永磁电机用测试台

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020180424A1 (en) 2019-03-04 2020-09-10 Iocurrents, Inc. Data compression and communication using machine learning
CN110472349B (zh) * 2019-08-20 2023-04-07 武汉科技大学 一种基于eemd和深度卷积网络的热轧钢性能预测方法
CN111008738B (zh) * 2019-12-04 2023-05-30 云南锡业集团(控股)有限责任公司研发中心 基于多模态深度学习的Sn-Bi合金延伸率和抗拉强度的预测方法
CN111008477B (zh) * 2019-12-06 2023-03-24 马鞍山钢铁股份有限公司 一种基于冷轧镀锌带钢力学性能调整工艺参数的方法
CA3181997A1 (en) * 2020-06-15 2021-12-23 Jfe Steel Corporation Mechanical property measuring apparatus, mechanical property measuring method, substance manufacturing equipment, substance management method, and substance manufacturing method
US12092609B2 (en) 2020-06-15 2024-09-17 Jfe Steel Corporation Mechanical property measuring apparatus, mechanical property measuring method, substance manufacturing equipment, substance management method, and substance manufacturing method
US12099033B2 (en) 2020-06-15 2024-09-24 Jfe Steel Corporation Mechanical property measuring apparatus, mechanical property measuring method, substance manufacturing equipment, substance management method, and substance manufacturing method
CN112100745B (zh) * 2020-09-15 2023-10-13 东北大学 基于lda理论的汽车大梁钢力学性能预测方法
CN112435234A (zh) * 2020-11-23 2021-03-02 东北大学 基于深度置信神经网络的热连轧带钢头部厚度预测方法
CN112598026B (zh) * 2020-12-07 2023-09-22 东北大学 基于不平衡数据深度学习的连退生产过程故障诊断方法
CN112561182B (zh) * 2020-12-21 2024-02-02 东北大学 一种基于熔体结构解析的熔渣黏度预测方法
EP4124398B1 (de) 2021-07-27 2024-04-10 Primetals Technologies Austria GmbH Verfahren zur bestimmung mechanischer eigenschaften eines walzgutes mit hilfe eines hybriden modells
CN113936744B (zh) * 2021-10-22 2024-03-29 兰州理工大学 热连轧流程生产钒微合金结构钢的在线性能预测方法及系统
CN114386325B (zh) * 2021-12-27 2023-04-28 北京科技大学 一种基于规则寻优的带钢力学性能预报方法
CN114580272B (zh) * 2022-02-16 2024-06-25 昆明贵金属研究所 同时对多元电接触合金导电率和硬度进行优化的设计方法
CN116779062A (zh) * 2023-06-14 2023-09-19 东北大学 基于moea/d智能优化算法的热轧带钢生产工艺动态优化方法
CN116825253B (zh) * 2023-07-03 2023-11-21 东北大学 基于特征选择的热轧带钢力学性能预测模型的建立方法
CN118068788B (zh) * 2024-02-07 2024-09-24 北京科技大学 一种模型与数据联合驱动热连轧过程故障诊断方法和装置
CN118571374A (zh) * 2024-05-29 2024-08-30 上海交通大学 一种tc4钛合金卷带力学性能预测方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102560081A (zh) * 2012-02-27 2012-07-11 宝山钢铁股份有限公司 一种基于带钢力学性能预报模型的加热炉节能控制方法
CN102632082A (zh) * 2011-02-11 2012-08-15 宝山钢铁股份有限公司 基于性能预报模型的热轧带钢力学性能的动态控制方法
CN106326525A (zh) * 2016-08-02 2017-01-11 武汉科技大学 基于广义可加模型的热轧变形抗力预报方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101046682A (zh) * 2007-04-28 2007-10-03 武汉科技大学 一种预测热轧含Nb带钢组织及力学性能的方法
CN104694720B (zh) * 2015-03-31 2017-05-17 北京首钢股份有限公司 一种热轧板卷力学性能预测与判定系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102632082A (zh) * 2011-02-11 2012-08-15 宝山钢铁股份有限公司 基于性能预报模型的热轧带钢力学性能的动态控制方法
CN102560081A (zh) * 2012-02-27 2012-07-11 宝山钢铁股份有限公司 一种基于带钢力学性能预报模型的加热炉节能控制方法
CN106326525A (zh) * 2016-08-02 2017-01-11 武汉科技大学 基于广义可加模型的热轧变形抗力预报方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YONG-HONG ZHAO 等: "Prediction of Mechanical Properties of Hot Rolled Strip by Using Semi-Parametric Single-Index Model", 《JOURNAL OF IRON AND STEEL RESEARCH, INTERNATIONAL》 *
向嵩 等: "低碳钢碳氮析出物的热力学计算", 《北京科技大学学报》 *
唐广波 等: "复合微合金化高强度低合金钢奥氏体相中碳氮化物析出热力学数值模拟", 《金属热处理》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108182340A (zh) * 2017-12-11 2018-06-19 中国科学院合肥物质科学研究院 一种高固溶性多主元合金的快速筛选方法
CN108182340B (zh) * 2017-12-11 2023-02-24 中国科学院合肥物质科学研究院 一种高固溶性多主元合金的快速筛选方法
CN108428038A (zh) * 2018-01-23 2018-08-21 安徽大学 基于gam模型的土壤锌含量影响因素分析方法
CN109145493A (zh) * 2018-09-11 2019-01-04 张家港宏昌钢板有限公司 开发螺纹钢新品的方法
CN109902872A (zh) * 2019-02-28 2019-06-18 北京首钢股份有限公司 一种性能预测模型与产品研发相结合的方法和装置
CN109902872B (zh) * 2019-02-28 2021-10-01 北京首钢股份有限公司 一种性能预测模型与产品研发相结合的方法和装置
CN110797096A (zh) * 2019-10-30 2020-02-14 合肥通用机械研究院有限公司 一种基于深度学习的耐热合金高温力学性能预测方法
CN110797096B (zh) * 2019-10-30 2021-04-27 合肥通用机械研究院有限公司 一种基于深度学习的耐热合金高温力学性能预测方法
CN111931425A (zh) * 2020-09-15 2020-11-13 山东钢铁集团日照有限公司 一种基于rf算法的热轧q355b钢种的性能预测系统
CN117110871A (zh) * 2023-10-13 2023-11-24 北京中航科电测控技术股份有限公司 高功率密度永磁电机用测试台
CN117110871B (zh) * 2023-10-13 2024-05-14 北京中航科电测控技术股份有限公司 高功率密度永磁电机用测试台

Also Published As

Publication number Publication date
CN106971026B (zh) 2020-06-16
US20180260717A1 (en) 2018-09-13
US11475317B2 (en) 2022-10-18

Similar Documents

Publication Publication Date Title
CN106971026A (zh) 基于全局可加模型的微合金钢力学性能预报方法
Sohrabi et al. Deformation-induced martensite in austenitic stainless steels: a review
CA2636898C (en) Apparatus and method for constructing prediction model
Guo et al. Modelling the correlation between processing parameters and properties of maraging steels using artificial neural network
Ray et al. Transformation textures in steels
CN104602830B (zh) 材料组织预测装置、产品制造方法及材料组织预测方法
Sha et al. Phase transformations in maraging steels
Li et al. Machine learning guided automatic recognition of crystal boundaries in bainitic/martensitic alloy and relationship between boundary types and ductile-to-brittle transition behavior
CN105740467A (zh) 一种C-Mn钢工业大数据的挖掘方法
Madej et al. Perceptive comparison of mean and full field dynamic recrystallization models
Łach et al. The evolution of the microstructure in AISI 304L stainless steel during the flat rolling–modeling by frontal cellular automata and verification
CN105980592A (zh) 双相不锈钢
CN107092751A (zh) 基于Bootstrap的变权重模型组合预测方法
CN110442954A (zh) 基于物理冶金学指导下机器学习的超高强不锈钢设计方法
CN104881707A (zh) 一种基于集成模型的烧结能耗预测方法
Hao et al. Interpolation-based anisotropic yield and hardening models
Fattouh Hybrid MCDM Model of ARAS-TOPSIS-GRA for Materials Selection Problem
CN114896770B (zh) 一种低合金高强钢连续冷却组织相体积分数计算方法
CN116779062A (zh) 基于moea/d智能优化算法的热轧带钢生产工艺动态优化方法
CN107922991A (zh) 对钢材进行淬火之际的冷却时间的导出方法、钢材的淬火方法以及钢材的淬火回火方法
Colla et al. Model parameters optimisation for an industrial application: A comparison between traditional approaches and genetic algorithms
Uranga et al. Optimization of rolling conditions in Nb microalloyed steel processed by thin slab casting and direct rolling route: processing maps
Bzowski et al. Sensitivity analysis of phase transformation model based on solution of diffusion equation
Sha Ultra high-strength maraging steel
CN110442953A (zh) 基于物理冶金学指导下机器学习的q&p钢的设计方法

Legal Events

Date Code Title Description
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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20170721

Assignee: Hubei Jiuge Energy Technology Co.,Ltd.

Assignor: WUHAN University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2024980011182

Denomination of invention: Prediction method for mechanical properties of microalloyed steel based on global additive model

Granted publication date: 20200616

License type: Common License

Record date: 20240808

EE01 Entry into force of recordation of patent licensing contract