CN112380690A - 一种亚稳复相金属材料流变应力模型的构建方法 - Google Patents

一种亚稳复相金属材料流变应力模型的构建方法 Download PDF

Info

Publication number
CN112380690A
CN112380690A CN202011255111.9A CN202011255111A CN112380690A CN 112380690 A CN112380690 A CN 112380690A CN 202011255111 A CN202011255111 A CN 202011255111A CN 112380690 A CN112380690 A CN 112380690A
Authority
CN
China
Prior art keywords
model
stress
metastable
metal material
phase
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
CN202011255111.9A
Other languages
English (en)
Other versions
CN112380690B (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.)
Anhui University of Technology AHUT
Original Assignee
Anhui University of Technology AHUT
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 Anhui University of Technology AHUT filed Critical Anhui University of Technology AHUT
Priority to CN202011255111.9A priority Critical patent/CN112380690B/zh
Publication of CN112380690A publication Critical patent/CN112380690A/zh
Application granted granted Critical
Publication of CN112380690B publication Critical patent/CN112380690B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

本发明公开了一种亚稳复相金属材料流变应力模型的构建方法，属于金属材料分析技术领域，包括以下步骤：S1：构建组成相体积分数物理模型；S2：构建组成相应力物理模型；S3：构建流变应力模型；S4：计算得到流变应力模型。本发明针对含亚稳相金属材料，在建立了具有物理冶金意义、体现组成相对其流变应力贡献的应力‑应变模型的同时，还便于深入区分理解组成相在亚稳金属材料加工硬化过程中的作用和贡献；另外采用本发明的模型构建方法及模型可编写成标准计算软件，进一步提高获取相关信息的效率，值得被推广使用。

Description

一种亚稳复相金属材料流变应力模型的构建方法
技术领域
本发明涉及金属材料分析技术领域，具体涉及一种亚稳复相金属材料流变应力模型的构建方法。
背景技术
先进高强度金属材料既是发展高新技术的基础也是制造环保型设备的关键，在汽车、石油化工等行业有着广泛需求。虽然均匀结构、单相超细晶金属材料具有超高强度，但伴随的有限加工硬化导致的极小均匀塑性变形能力限制了其作为结构材料的用途。
近年来的研究结果表明，将显微结构设计为高强度基体+亚稳相构成的复相组织具有同时提高金属材料强度和塑性的巨大潜力。含亚稳相的超细晶金属材料具有良好塑性变形能力的原因在于亚稳相具有的TRIP效应(相变诱导塑性)可使材料具有优异的加工硬化能力。因此，建立体现含亚稳相金属材料中各组成相对流变应力或加工硬化率贡献的模型将有利于深入理解含亚稳相金属材料的塑性变形机理，以便基于塑性变形机理开发相关调控方法进一步改善其塑性变形能力。
目前，所建立的相关模型多针对单个应变诱导相变过程，但是实际上近期开发的含亚稳相金属材料存在多个应变诱导相变过程。例如，亚稳超细晶奥氏体不锈钢会同时发生奥氏体向体心立方马氏体、奥氏体向六方马氏体和六方马氏体向体心立方马氏体的应变诱导相变过程。现有基于单个应变诱导相变的流变应力模型显然无法准确描述具有多个应变诱导相变过程的金属材料的塑性变形行为。因此，有必要在现有明确物理含义模型的基础上针对含多个应变诱导相变过程的金属材料，建立体现各组成相对流变应力或加工硬化率贡献的模型以便深入；理解这类金属材料的塑性变形机理。为此，提出一种亚稳复相金属材料流变应力模型的构建方法。
发明内容
本发明所要解决的技术问题在于：如何建立体现各组成相对流变应力或加工硬化率贡献的模型以便深入理解这类金属材料的塑性变形机理，提供了一种亚稳复相金属材料流变应力模型的构建方法。
本发明是通过以下技术方案解决上述技术问题的，本发明包括以下步骤：
S1：使用Olson-Cohen模型或其他具有物理意义的模型构建亚稳金属材料在塑性变形过程中产生的新相P1，P2，…Pi的体积分数VP1，VP2…VPi随应变(ε)变化的关系；
步骤S1的具体过程如下：
S11：分别将P1作为一个整体，(P1+P2)作为一个整体，…(P1+P2+…+Pi)作为一个整体，使用Olson-Cohen模型VP＝1-(1-V0)exp{-β[1-exp(-αε)]m}或其他具有物理意义的模型建立它们的体积分数VP1，VP1+P2，…VP1+P2+…+Pi随应变ε变化的关系，如下所示；
VP1＝1-(1-V1 0)exp{-β1[1-exp(-α1ε)]m1} (1-1)
VP1+P2＝1-(1-V2 0)exp{-β2[1-exp(-α2ε)]m2} (1-2)
…
VP1+P2+…+Pi＝1-(1-Vi 0)exp{-βi[1-exp(-αiε)]mi} (1-i)
S12：那么，亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi的体积分数VP0，VP1，VP2…VPi可表示为；
VP0＝1-VP1+P2+…+Pi (2-0)
VP1＝VP1 (2-1)
VP2＝VP1+P2-VP1 (2-2)
…
VPi＝VP1+P2+…+Pi-VP1+P2+…+Pi-1 (2-i)
S13：将(P1+P2+…+Pi)作为一个整体、(P2+P3+…+Pi)作为一个整体、…(Pi-1+Pi)作为一个整体和Pi的体积分数VP1+P2+…+Pi、VP2+P3+…+Pi、…VPi-1+Pi和VPi表示如下：
VP1+P2+…+Pi＝VP1+P2+…+Pi (3-1)
VP2+P3+…+Pi＝VP1+P2+…+Pi-VP1 (3-2)
…
VPi-1+Pi＝VP1+P2+…+Pi-VP1+P2+…+Pi-2 (3-i-1)
VPi＝VPi (3-i)
S2：构建亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi的应力随应变ε变化的物理模型；
步骤S2的具体过程如下：
S21：使用Kocks-Mecking模型(式4)或其他具有物理意义的模型描述只有亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi单独构成的材料的拉伸应力FP0，FP1，FP2，…FPi随应变ε变化的关系，如下所示：
Figure BDA0002772868880000021
Figure BDA0002772868880000022
Figure BDA0002772868880000031
…
Figure BDA0002772868880000032
S22：使用Kocks-Mecking模型分别描述由(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi构成的材料的拉伸应力FP1+P2+…+Pi、FP2+P3+…+Pi、…FPi随应变ε变化的关系，如下所示：
Figure BDA0002772868880000033
Figure BDA0002772868880000034
…
FPi＝FPi (5-i)
S3：将亚稳金属材料看作由基体相P0和(P1+P2+…Pi)整体构成，将P1，P2，…Pi组成的材料看作由P1和(P2+P3+…Pi)整体构成，将P2，P3，…Pi组成的材料看作由P2和(P3+P4+…Pi)整体构成，…将Pi-2，Pi-1，…Pi组成的材料看作由Pi-2和(Pi-1+Pi)整体构成，使用修正Gladman混合定则(式6)构建体现各组成相对亚稳金属材料流变应力贡献的模型，如下所示：
Figure BDA0002772868880000035
Figure BDA0002772868880000036
Figure BDA0002772868880000037
…
Figure BDA0002772868880000038
Figure BDA0002772868880000039
其中：
Figure BDA00027728688800000310
Figure BDA00027728688800000311
…
Figure BDA0002772868880000041
S4：将步骤S1和步骤S2中构建的组成相体积分数物理模型和组成相应力物理模型嵌套在步骤S3中构建的流变应力模型，便得到了体现亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对其流变应力(S)贡献的应力-应变模型，如下所示：
Figure BDA0002772868880000042
步骤S2的具体过程如下：
S41：基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对亚稳金属材料流变应力S的贡献S0，S1，S2，…Si可如下式表示；
Figure BDA0002772868880000043
Figure BDA0002772868880000044
Figure BDA0002772868880000045
…
Figure BDA0002772868880000046
Si＝S-S1-S2-...Si-1 (9-i)
S42：对式9-0至式9-i求导便得到了基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对亚稳金属材料加工硬化率的贡献(dS/dε)P0，(dS/dε)P1，(dS/dε)P2，…(dS/dε)Pi。
本发明相比现有技术具有以下优点：该亚稳复相金属材料流变应力模型的构建方法，针对含亚稳相金属材料，在建立了具有物理冶金意义、体现组成相对其流变应力贡献的应力-应变模型的同时，还便于深入区分理解组成相在亚稳金属材料加工硬化过程中的作用和贡献；另外采用本发明的模型构建方法及模型可编写成标准计算软件，进一步提高获取相关信息的效率，值得被推广使用。
附图说明
图1是本发明实施例一中亚稳复相金属材料流变应力模型构建方法的流程示意图；
图2是本发明实施例二中实测亚稳超细晶奥氏体不锈钢的各组成相体积分数和拟合结果；
图3是本发明实施例二中实测各相对亚稳超细晶奥氏体不锈钢的流变应力贡献和拟合结果；
图4是本发明实施例二中实测各相对亚稳超细晶奥氏体不锈钢的加工硬化率贡献和拟合结果。
具体实施方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例一
如图1所示，本实施例提供一种技术方案：一种亚稳复相金属材料流变应力模型的构建方法，包括以下步骤：
S1：使用Olson-Cohen模型或其他具有物理意义的模型构建亚稳金属材料在塑性变形过程中产生的新相P1，P2，…Pi的体积分数VP1，VP2…VPi随应变(ε)变化的关系；
步骤S1的具体过程如下：
S11：分别将P1作为一个整体，(P1+P2)作为一个整体，…(P1+P2+…+Pi)作为一个整体，使用Olson-Cohen模型VP＝1-(1-V0)exp{-β[1-exp(-αε)]m}或其他具有物理意义的模型建立它们的体积分数VP1，VP1+P2，…VP1+P2+…+Pi随应变ε变化的关系，如下所示；
VP1＝1-(1-V1 0)exp{-β1[1-exp(-α1ε)]m1} (1-1)
VP1+P2＝1-(1-V2 0)exp{-β2[1-exp(-α2ε)]m2} (1-2)
…
VP1+P2+…+Pi＝1-(1-Vi 0)exp{-βi[1-exp(-αiε)]mi} (1-i)
S12：那么，亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi的体积分数VP0，VP1，VP2…VPi可表示为；
VP0＝1-VP1+P2+…+Pi (2-0)
VP1＝VP1 (2-1)
VP2＝VP1+P2-VP1 (2-2)
…
VPi＝VP1+P2+…+Pi-VP1+P2+…+Pi-1 (2-i)
S13：将(P1+P2+…+Pi)作为一个整体、(P2+P3+…+Pi)作为一个整体、…(Pi-1+Pi)作为一个整体和Pi的体积分数VP1+P2+…+Pi、VP2+P3+…+Pi、…VPi-1+Pi和VPi表示如下：
VP1+P2+…+Pi＝VP1+P2+…+Pi (3-1)
VP2+P3+…+Pi＝VP1+P2+…+Pi-VP1 (3-2)
…
VPi-1+Pi＝VP1+P2+…+Pi-VP1+P2+…+Pi-2 (3-i-1)
VPi＝VPi (3-i)
S2：构建亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi的应力随应变ε变化的物理模型；
步骤S2的具体过程如下：
S21：使用Kocks-Mecking模型(式4)或其他具有物理意义的模型描述只有亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi单独构成的材料的拉伸应力FP0，FP1，FP2，…FPi随应变ε变化的关系，如下所示：
Figure BDA0002772868880000061
Figure BDA0002772868880000062
Figure BDA0002772868880000063
…
Figure BDA0002772868880000064
S22：使用Kocks-Mecking模型分别描述由(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi构成的材料的拉伸应力FP1+P2+…+Pi、FP2+P3+…+Pi、…FPi随应变ε变化的关系，如下所示：
Figure BDA0002772868880000065
Figure BDA0002772868880000066
…
FPi＝FPi (5-i)
S3：将亚稳金属材料看作由基体相P0和(P1+P2+…Pi)整体构成，将P1，P2，…Pi组成的材料看作由P1和(P2+P3+…Pi)整体构成，将P2，P3，…Pi组成的材料看作由P2和(P3+P4+…Pi)整体构成，…将Pi-2，Pi-1，…Pi组成的材料看作由Pi-2和(Pi-1+Pi)整体构成，使用修正Gladman混合定则(式6)构建体现各组成相对亚稳金属材料流变应力贡献的模型，如下所示：
Figure BDA0002772868880000071
Figure BDA0002772868880000072
Figure BDA0002772868880000073
…
Figure BDA0002772868880000074
Figure BDA0002772868880000075
其中：
Figure BDA0002772868880000076
Figure BDA0002772868880000077
…
Figure BDA0002772868880000078
S4：将步骤S1和步骤S2中构建的组成相体积分数物理模型和组成相应力物理模型嵌套在步骤S3中构建的流变应力模型，便得到了体现亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对其流变应力(S)贡献的应力-应变模型，如下所示：
Figure BDA0002772868880000079
步骤S2的具体过程如下：
S41：基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对亚稳金属材料流变应力S的贡献S0，S1，S2，…Si可如下式表示；
Figure BDA00027728688800000710
Figure BDA00027728688800000711
Figure BDA0002772868880000081
…
Figure BDA0002772868880000082
Si＝S-S1-S2-...Si-1 (9-i)
S42：对式9-0至式9-i求导便得到了基体相P0和在塑性变形过程中产生的新相P1，P2，…Pi对亚稳金属材料加工硬化率的贡献(dS/dε)P0，(dS/dε)P1，(dS/dε)P2，…(dS/dε)Pi。
实施例二
下述实施例中亚稳超细晶奥氏体不锈钢，其按重量百分比，合金元素包括0.010～0.080％的C、≤2.00％的Mn、≤0.035％的P、≤0.020％的S、≤0.75％的Si、8.00～10.50％的Ni、18.00～20.00％的Cr、≤0.10％的N，余量为Fe。
S1、使用Olson-Cohen模型构建亚稳超细晶奥氏体不锈钢在塑性变形过程中产生的体心立方马氏体(P1)和六方马氏体(P2)的体积分数VP1，VP2随ε变化的关系；
S11、分别将P1作为一个整体，(P1+P2)作为一个整体，使用Olson-Cohen模型VP＝1-(1-V0)exp{-β[1-exp(-αε)]m}建立它们的体积分数VP1，VP1+P2随应ε变化的关系，如下所示；
VP1＝1-(1-V1 0)exp{-β1[1-exp(-α1ε)]m1} (10-1)
VP1+P2＝1-(1-V2 0)exp{-β2[1-exp(-α2ε)]m2} (10-2)
S12、那么，亚稳超细晶奥氏体不锈钢的基体相P0和在塑性变形过程中产生的新相P1，P2的体积分数VP0，VP1，VP2可表示为；
VP0＝1-VP1+P2 (11-0)
VP1＝VP1 (11-1)
VP2＝VP1+P2-VP1 (11-2)
S13：将(P1+P2)作为一个整体和P2作为一个整体的体积分数VP1+P2和VP2表示如下：
VP1+P2＝VP1+P2 (12-1)
VP2＝VP2 (12-2)
S2、构建亚稳超细晶奥氏体不锈钢基体相P0和在塑性变形过程中产生的新相P1，P2的应力随ε变化的物理模型；
S21、使用Kocks-Mecking模型描述只有基体相P0和在塑性变形过程中产生的新相P1，P2单独构成的材料的拉伸应力FP0，FP1，FP2随ε变化的关系；
Figure BDA0002772868880000091
Figure BDA0002772868880000092
Figure BDA0002772868880000093
S22、使用Kocks-Mecking模型描述由(P1+P2)构成的材料的拉伸应力FP1+P2，随ε变化的关系；
Figure BDA0002772868880000094
S3、将亚稳超细晶奥氏体不锈钢看作由基体相P0和(P1+P2)整体构成，将P1，P2组成的材料由P1和P2构成，使用修正Gladman混合定则构建体现各组成相对亚稳超细晶奥氏体不锈钢流变应力(S)贡献的模型；
Figure BDA0002772868880000095
Figure BDA0002772868880000096
其中，
Figure BDA0002772868880000097
Figure BDA0002772868880000098
S4、将S1步骤和S2步骤中构建的组成相体积分数物理模型和组成相应力物理模型嵌套在S3步骤中构建的流变应力模型，便得到了体现亚稳超细晶奥氏体不锈钢基体相P0和在塑性变形过程中产生的新相P1，P2对其流变应力(S)贡献的应力-应变模型；
Figure BDA0002772868880000099
S41、基体相P0和在塑性变形过程中产生的新相P1，P2对亚稳超细晶奥氏体不锈钢流变应力(S)的贡献可如下式表示；
Figure BDA00027728688800000910
Figure BDA00027728688800000911
S2＝S-S1-S0 (18-2)
S42、对式18-0至式18-2求导便得到了基体相P0和在塑性变形过程中产生的新相P1，P2对亚稳超细晶奥氏体不锈钢加工硬化率的贡献(dS/dε)P0，(dS/dε)P1和(dS/dε)P2。
为验证所建模型的合理性：
将亚稳超细晶奥氏体不锈钢进行基于同步辐射的高能X射线原位拉伸试验，以实测的各相体积分数为纵坐标轴、ε为横坐标轴，使用S11步骤中式10-1和式10-2拟合，拟合参数见表1，结果与实测数据符合的很好，如图2所示(□，○，△和▽代表实测数据，实线代表拟合结果)；
表1 Olson-Cohen模型拟合实测数据获得的参数表
Figure BDA0002772868880000101
将亚稳超细晶奥氏体不锈钢进行基于同步辐射的高能X射线原位拉伸试验，以实测的各相对其流变应力贡献为纵坐标轴、ε为横坐标轴，采用P0和P1对应的F0、A、G、b实际物理参数(表2)以及表1的参数作为输入参数，使用S41步骤构建的模型进行拟合，拟合参数见表3，结果与实测数据符合的很好，如图3所示(□，○，△，▽和◇代表实测数据，实线代表拟合结果)；
表2基体相P0和新相P1对应的F0、A、G、b实际物理参数表
Figure BDA0002772868880000102
表3所建应力-应变模型的拟合参数表
Figure BDA0002772868880000103
更进一步地，以实测的各相对流变应力贡献的导数为纵坐标轴、ε为横坐标轴，使用S41步骤构建的模型的导数形式拟合，结果与实测数据也符合的很好，如图4所示(□，○，△和■代表实测数据，实线代表拟合结果)；
综上所述，上述实施例的亚稳复相金属材料流变应力模型的构建方法，针对含亚稳相金属材料，在建立了具有物理冶金意义、体现组成相对其流变应力贡献的应力-应变模型的同时，还便于深入区分理解组成相在亚稳金属材料加工硬化过程中的作用和贡献；另外采用本发明的模型构建方法及模型可编写成标准计算软件，进一步提高获取相关信息的效率，值得被推广使用。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (6)

1.一种亚稳复相金属材料流变应力模型的构建方法，其特征在于，包括以下步骤：
S1：构建组成相体积分数物理模型
构建亚稳金属材料在塑性变形过程中产生的所有新相的体积分数随应变变化的物理模型；
S2：构建组成相应力物理模型
构建亚稳金属材料在塑性变形过程中组成相的应力随应变变化的物理模型；
S3：构建流变应力模型
构建包含亚稳金属材料在塑性变形过程中组成相的体积分数和应力在内的流变应力模型；
S4：计算得到流变应力模型
将步骤S1构建的组成相体积分数物理模型和步骤S2中构建的组成相应力物理模型嵌套在步骤S3构建的流变应力模型中，得到描述多个应变诱导相变过程的亚稳金属材料流变应力模型。
2.根据权利要求1所述的一种亚稳复相金属材料流变应力模型的构建方法，其特征在于：所述步骤S1具体包括以下步骤：
S11：分别将P1作为一个整体、(P1+P2)作为一个整体、…(P1+P2+…+Pi)作为一个整体，使用Olson-Cohen模型VP＝1-(1-V0)exp{-β[1-exp(-αε)]m}建立它们的体积分数VP1、VP1+P2、…VP1+P2+…+Pi随应ε变化的关系，如下所示；
VP1＝1-(1-V1 0)exp{-β1[1-exp(-α1ε)]m1}
VP1+P2＝1-(1-V2 0)exp{-β2[1-exp(-α2ε)]m2}
…
VP1+P2+…+Pi＝1-(1-Vi 0)exp{-βi[1-exp(-αiε)]mi}
其中，P1、P2、…Pi分别为亚稳金属材料在塑性变形过程中产生的新相；V1 0、V2 0、…Vi 0分别为P1、P2、…(P1+P2+…+Pi)对应的初始体积分数；β1、β2、…βi分别为P1、P2、…(P1+P2+…+Pi)对应Olson-Cohen模型的第一系数；α1、α2、…αi分别为P1、P2、…(P1+P2+…+Pi)对应Olson-Cohen模型的第二系数；m1、m2、…mi分别为P1、P2、…(P1+P2+…+Pi)对应Olson-Cohen模型的第三系数；
S12：将亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1、P2、…Pi的体积分数VP0、VP1、VP2…VPi表示如下：
VP0＝1-VP1+P2+…+Pi
VP1＝VP1
VP2＝VP1+P2-VP1
…
VPi＝VP1+P2+…+Pi-VP1+P2+…+Pi-1；
S13：将(P1+P2+…+Pi)作为一个整体、(P2+P3+…+Pi)作为一个整体、…(Pi-1+Pi)作为一个整体和Pi的体积分数VP1+P2+…+Pi、VP2+P3+…+Pi、…VPi-1+Pi和VPi表示如下：
VP1+P2+…+Pi＝VP1+P2+…+Pi
VP2+P3+…+Pi＝VP1+P2+…+Pi-VP1
…
VPi-1+Pi＝VP1+P2+…+Pi-VP1+P2+…+Pi-2
VPi＝VP。i
3.根据权利要求2所述的一种亚稳复相金属材料流变应力模型的构建方法，其特征在于：所述步骤S2具体包括以下步骤：
S21：使用Kocks-Mecking模型描述只有亚稳金属材料基体相P0和在塑性变形过程中产生的新相P1、P2、…Pi单独构成的材料的拉伸应力FP0、FP1、FP2、…FPi随应变ε变化的关系，如下所示：
Figure FDA0002772868870000021
Figure FDA0002772868870000022
Figure FDA0002772868870000023
…
Figure FDA0002772868870000024
其中，F0为点阵摩擦力；A、G和b分别是与Taylor因子有关的常数、切变模量和伯氏矢量；ε0为塑性变形开始时的应变；ρ0 P0、ρ0 P1、ρ0 P2、…ρ0 Pi分别为P0、P1、P2、…Pi对应的初始位错密度；k1 P0、k1 P1、k1 P2、…k1 Pi分别为P0、P1、P2、…Pi对应的位错增殖速率；k2 P0、k2 P1、k2 P2、…k2 Pi分别为P0、P1、P2、…Pi对应的位错湮灭速率；
S22：使用Kocks-Mecking模型分别描述由(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi构成的材料的拉伸应力FP1+P2+…+Pi、FP2+P3+…+Pi、…FPi随应变ε变化的关系，如下所示：
Figure FDA0002772868870000025
Figure FDA0002772868870000031
…
FPi＝FPi
其中，ρ0 P1+P2+…+Pi、ρ0 P2+P3+…+Pi、…ρ0 Pi分别为(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi对应的初始位错密度；k1 P1+P2+…+Pi、k1 P2+P3+…+Pi、…k1 Pi分别为(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi对应的位错增殖速率；k2 P1+P2+…+Pi、k2 P2+P3+…+Pi、…k2 Pi分别为(P1+P2+…+Pi)、(P2+P3+…+Pi)、…Pi对应的位错湮灭速率。
4.根据权利要求3所述的一种亚稳复相金属材料流变应力模型的构建方法，其特征在于，所述步骤S3的过程具体为：将亚稳金属材料看作由基体相P0和(P1+P2+…Pi)整体构成，将P1，P2，…Pi组成的材料看作由P1和(P2+P3+…Pi)整体构成，将P2、P3、…Pi组成的材料看作由P2和(P3+P4+…Pi)整体构成，以此类推，将Pi-2、Pi-1、…Pi组成的材料看作由Pi-2和(Pi-1+Pi)整体构成，使用修正Gladman混合定则构建体现各组成相对亚稳金属材料流变应力贡献的模型，如下所示：
Figure FDA0002772868870000032
Figure FDA0002772868870000033
Figure FDA0002772868870000034
…
Figure FDA0002772868870000035
Figure FDA0002772868870000036
其中，S表示流变应力；n1、n2、n3、...ni-1、ni分别是(P0+P1+P2+…+Pi)、(P1+P2+…+Pi)、(P2+P3+…+Pi)、…(Pi-2+Pi-1+Pi)、(Pi-1+Pi)对应的Gladman混合定则系数；
Figure FDA0002772868870000037
Figure FDA0002772868870000038
Figure FDA0002772868870000039
其中，a1、a2、...ai分别是(P0+P1+P2+…+Pi)、(P1+P2+…+Pi)、…Pi-1+Pi对应的第一拟合系数；l1、l2、...li分别是(P0+P1+P2+…+Pi)、(P1+P2+…+Pi)、…Pi-1+Pi对应的第二拟合系数。
5.根据权利要求4所述的一种亚稳复相金属材料流变应力模型的构建方法，其特征在于：在所述步骤S4中，描述多个应变诱导相变过程的亚稳金属材料流变应力模型如下所示：
Figure FDA0002772868870000041
6.根据权利要求5所述的一种亚稳复相金属材料流变应力模型的构建方法，其特征在于：所述步骤S4具体包括以下步骤：
S41：将基体相P0和在塑性变形过程中产生的新相P1、P2、…Pi对亚稳金属材料流变应力S的贡献S0、S1、S2、…Si表示如下：
Figure FDA0002772868870000042
Figure FDA0002772868870000043
Figure FDA0002772868870000044
…
Figure FDA0002772868870000045
Si＝S-S1-S2-...Si-1；
S42：对步骤S41中的各式求导数，得到基体相P0和在塑性变形过程中产生的新相P1、P2、…Pi对亚稳金属材料加工硬化率的贡献(dS/dε)P0、(dS/dε)P1、(dS/dε)P2、…(dS/dε)Pi。
CN202011255111.9A 2020-11-11 2020-11-11 一种亚稳复相金属材料流变应力模型的构建方法 Active CN112380690B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011255111.9A CN112380690B (zh) 2020-11-11 2020-11-11 一种亚稳复相金属材料流变应力模型的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011255111.9A CN112380690B (zh) 2020-11-11 2020-11-11 一种亚稳复相金属材料流变应力模型的构建方法

Publications (2)

Publication Number Publication Date
CN112380690A true CN112380690A (zh) 2021-02-19
CN112380690B CN112380690B (zh) 2022-07-26

Family

ID=74582738

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011255111.9A Active CN112380690B (zh) 2020-11-11 2020-11-11 一种亚稳复相金属材料流变应力模型的构建方法

Country Status (1)

Country Link
CN (1) CN112380690B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB771932A (en) * 1953-10-06 1957-04-10 Commw Of Australia Method of and apparatus for indicating the plastic deformation of materials
US20100143181A1 (en) * 2008-07-18 2010-06-10 Northwestern University High strength austenitic trip steel
CN106383928A (zh) * 2016-08-30 2017-02-08 西华大学 一种粉末烧结材料高温塑性流变致密化本构模型构建方法
US20170101702A1 (en) * 2015-10-07 2017-04-13 Seoul National University R&Db Foundation Metallic glass composites with controllable work-hardening capacity
CN107991179A (zh) * 2017-11-03 2018-05-04 合肥通用机械研究院 一种测定应变诱导马氏体相变动力学曲线的方法
CN110442954A (zh) * 2019-07-31 2019-11-12 东北大学 基于物理冶金学指导下机器学习的超高强不锈钢设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB771932A (en) * 1953-10-06 1957-04-10 Commw Of Australia Method of and apparatus for indicating the plastic deformation of materials
US20100143181A1 (en) * 2008-07-18 2010-06-10 Northwestern University High strength austenitic trip steel
US20170101702A1 (en) * 2015-10-07 2017-04-13 Seoul National University R&Db Foundation Metallic glass composites with controllable work-hardening capacity
KR20170041594A (ko) * 2015-10-07 2017-04-17 서울대학교산학협력단 가공경화능 제어 비정질 합금 기지 복합재의 제조 방법 및 그에 따라 제조된 복합재
CN106383928A (zh) * 2016-08-30 2017-02-08 西华大学 一种粉末烧结材料高温塑性流变致密化本构模型构建方法
CN107991179A (zh) * 2017-11-03 2018-05-04 合肥通用机械研究院 一种测定应变诱导马氏体相变动力学曲线的方法
CN110442954A (zh) * 2019-07-31 2019-11-12 东北大学 基于物理冶金学指导下机器学习的超高强不锈钢设计方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
张玉兵等: "金属材料非平衡相变的热/动力学协同", 《中国材料进展》 *
强洪夫等: "一种新型亚稳β钛合金的热变形本构模型", 《稀有金属》 *
朱玉亮等: "描述稳定面心立方金属真应力与真应变关系的数学模型", 《钢铁研究学报》 *
李志杰等: "中碳钢温变形的力学特性", 《北京理工大学学报》 *
郭志凯: "高锰TRIP钢的组织控制与力学行为研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑(月刊)》 *

Also Published As

Publication number Publication date
CN112380690B (zh) 2022-07-26

Similar Documents

Publication Publication Date Title
Bleck et al. A comparative study of the forming-limit diagram models for sheet steels
Roberts et al. Dynamic recrystallization kinetics
Barraclough et al. Static recrystallization and restoration after hot deformation of type 304 stainless steel
Iwamoto et al. Computational simulation of the dependence of the austenitic grain size on the deformation behavior of TRIP steels
Sieurin et al. Modelling solid solution hardening in stainless steels
CN108645706B (zh) 一种通过硬化强度和抗拉强度预测金属材料疲劳强度的方法
Cai et al. Strain hardening behavior of high performance FBDP, TRIP and TWIP steels
Zhu et al. Constitutive modelling of hot deformation behaviour of metallic materials
Li et al. Stress relaxation in tensile deformation of 304 stainless steel
Yang et al. Effect of cooling rate and austenite deformation on hardness and microstructure of 960MPa high strength steel
CN113094885B (zh) 含缺陷结构高熵合金强度预测方法
CN112380690B (zh) 一种亚稳复相金属材料流变应力模型的构建方法
Kumar et al. Numerical evaluation of JR curve using small punch test data
Wu et al. Composition-processing-property correlation mining of Nb–Ti microalloyed steel based on industrial data
Donnini et al. Evaluation of microstructure quality in ductile irons based on tensile behaviour analysis
Majta et al. Modeling of ferrite structure after deformation in the two-phase region
Parnian et al. Effect of drawing strain on development of martensitic transformation and mechanical properties in AISI 304L stainless steel wire
Song et al. Effect of cooling rate and coiling temperature on microstructure and precipitation behavior of a 700 MPa weathering steel
Marrapu Effect of localization criteria and yield criteria in predicting the forming limit diagram (FLD) of DP590 steel sheets
Krauer et al. Forming limit prediction of metastable materials with temperature and strain induced martensite transformation
Han et al. The Influence of Chemical Element on Properties of Deformed Steel Bar
Stough et al. The Stress Corrosion and Pyrophoric Behavior of Titanium and Titanium Alloys
Szewczyk Investigation of the tensile deformation and fracture of a dual-phase steel
Liu et al. Strain-Induced Precipitation Behavior and Microstructure Evolution of Ti-V-Mo Complex Microalloyed Steel
Zhang et al. Dynamic and Postdeformation Recrystallization of Nuclear-Grade 316LN Stainless Steel

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