CN114757047B - 一种轴承钢材料m50合金的多尺度建模计算方法 - Google Patents

一种轴承钢材料m50合金的多尺度建模计算方法 Download PDF

Info

Publication number
CN114757047B
CN114757047B CN202210461687.3A CN202210461687A CN114757047B CN 114757047 B CN114757047 B CN 114757047B CN 202210461687 A CN202210461687 A CN 202210461687A CN 114757047 B CN114757047 B CN 114757047B
Authority
CN
China
Prior art keywords
alloy
phase
sample
steel material
bearing steel
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
CN202210461687.3A
Other languages
English (en)
Other versions
CN114757047A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202210461687.3A priority Critical patent/CN114757047B/zh
Publication of CN114757047A publication Critical patent/CN114757047A/zh
Application granted granted Critical
Publication of CN114757047B publication Critical patent/CN114757047B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

本发明公开了一种轴承钢材料M50合金的多尺度建模计算方法,属于多相合金多尺度建模领域;包括以下步骤:划定试样观测区域;在所述观测区域对试样进行物相种类、分布鉴定,元素组成和晶体结构分析,得到介观元素质量配比和不同物相的元素组成;根据物相种类和晶体结构建立单一物相的原子尺度模型;通过介观元素质量配比、物相元素组成建立等式,求解原子模型中每种物相的数量和铁基体适当尺寸;在铁基体中根据所得到的物相比例随机插入几种物相组织,并优化结构,得到轴承钢材料M50合金的微‑介观结构模型;可以基于模型在介观尺度上模拟合金的各项性能,仿真观察合金在热处理以及损伤时在原子尺度上的结构转变过程以及M50合金的损伤机理。

Description

一种轴承钢材料M50合金的多尺度建模计算方法
技术领域
本发明属于轴承钢多尺度模拟领域,具体涉及一种轴承钢材料M50合金的多尺度建模计算方法。
背景技术
航空发动机主轴承工作在高温、高速、乏油、时变冲击载荷环境下,极易出现损伤失效,属于航空发动机中的薄弱环节。从滚动轴承的宏观故障看,有点蚀剥落、打滑蹭伤、摩擦磨损和腐蚀等失效形式,而从材料结构层面看,材料的失效在微观尺度下有原子空位、间隙和置换等点缺陷以及位错等线缺陷;在介观尺度下有晶界、相界和层错等面缺陷,轴承的失效过程涉及到多个尺度下的耦合作用,仅从单一尺度进行模拟或试验分析难以揭示材料的变形及失效机理,因此需要开展多尺度下的分析方法研究,综合考虑多尺度效应的主轴承材料特征,进而揭示轴承典型故障从微观到宏观的演化过程。
M50合金是一种高温轴承钢,主要应用于航空发动机轴承,M50钢的基体为α-Fe,热处理后的合金中存在残余奥氏体、马氏体几种物相,并且添加了Cr、Mo、V、Ni等元素形成碳化物作为增强相存在,形成的碳化物的形态、颗粒大小、分布的均匀性都对钢的性能、使用寿命等都有着十分重要的影响。目前M50合金的现有研究主要集中在其处理办法、强化手段以及其优化其元素配比几个方向,对其微观尺度及多尺度建模方向展开的研究较少。而对于金属的多尺度研究领域,研究对象多为纯金属或元素种类较少的合金,对于元素种类多、掺杂物相复杂的合金研究较少;研究方法多为基于分子动力学各项理论的仿真模拟计算,实验作为最终验证手段,仿真与实验的结合度较低,对于M50合金并不完全适用;因此,现有的轴承钢材料M50合金的多尺度建模技术有待进一步改进。
发明内容
为了解决现有技术中存在的问题,本发明提供一种轴承钢材料M50合金的多尺度建模计算方法,在保证建模时效性的情况下,通过试验手段获得M50合金中的元素、物相、晶体结构信息,在模型中加入更多的元素和物相掺杂,提高所得微观结构的准确性和真实性。
为了实现上述目的,本发明采用的技术方案是:一种轴承钢材料M50合金的多尺度建模计算方法,包括以下步骤:
划定试样观测区域;
在所述观测区域对试样进行物相鉴定,获得试样含有的物相种类和空间点阵类型;
根据物相种类信息,获取所述观测区域试样的物相分布、尺寸以及取向等信息;
对所述观测区域内每种组织进行元素定性定量分析,得到材料介观尺度元素质量配比以及不同物相元素组成;
根据物相种类及其元素组成,使用Materials Studio软件建立每种物相组织的单胞模型,适当扩胞建成单一物相的原子尺度模型;
通过介观尺度元素质量配比和不同物相的元素组成建立等式,来求解原子模型中每种物相的数量和铁基体的适当尺寸即Fe原子数;
根据所求得的物相数量,在铁基体中随机插入几种物相组织,并优化结构,得到轴承钢材料M50合金的微-介观结构模型;
基于所述轴承钢材料M50合金的微-介观结构模型在原子尺度上模拟合金的各项性能,并仿真观察合金在热处理以及损伤时在原子尺度上的结构转变过程以及M50合金的损伤机理。
所述试样通过以下步骤制备:
使用中走丝线切割方法制备材料试样,进行一次淬火和三次回火热处理;
对热处理后试样进行预处理,将试样打磨抛光,完成测试前准备。
适用显微硬度计或纳米压痕仪完成观测区域划分,先手动使用记号笔圈出观测区域,使用记号笔在试样侧面做出标记,作为试样朝向的定位,然后根据划区精度要求和构建模型需要选择仪器,大区域划分使用显微硬度计设定载荷在试样上完成一次硬度测量,小区域划分使用纳米压痕仪,完成观测区域划分。
观测区域试样的物相种类和空间点阵类型鉴定采用XRD方法获取。
划定区域试样的物相分布、尺寸以及取向等信息采用EBSD方法获取。
划定区域内的介观尺度元素质量配比和每种物相的元素定性定量分析采用EPMA方法获取。
在Materials Studio软件建立每种组织的单胞模型时,其中,Fe-C、Fe-Mo、Mo-Cr、V-Mo原子之间的相互作用力均用其对应的meam势函数描述,Fe-Cr原子之间的相互作用力用cdeam势函数描述,Fe-V原子之间的相互作用力用eam/fs势函数描述,C元素与Cr、Mo、V元素和Cr-V元素之间均采用lj/cut势函数进行描述,lj/cut势函数的参数均由通用力场下拟合得来。
获得每种物相数量的同时获得材料微-介观结构模型的铁基体尺寸,同时能够保证M50合金组成元素的质量配比。
在铁基体中根据所得到的物相比例随机插入几种物相组织,并优化结构,得到轴承钢材料M50合金的微-介观结构模型,同时保证结构稳定且与实际观测到的组织一致。
进一步的,观测区域划定采用显微硬度计和纳米压痕仪,大区域划定采用显微硬度计,使用恰当载荷在试样中心或四角完成硬度测量,直接由凹痕确定观测区域或由四角凹痕完成观测区域划定;小区域划定采用纳米压痕仪,在大区域划分基础上进行更精细地划分,便于选择物相观察。
进一步的,采用电子探针显微分析仪可以对划定区域内每种物相进行元素定性定量分析,直接获得奥氏体、马氏体的元素组成,结合文献确定各类碳化物的种类。
与现有技术相比,本发明至少具有以下有益效果:
本发明在保证建模时效性的情况下,通过实验手段获得M50合金中的元素、物相、晶体结构信息,在模型中加入更多的元素和物相掺杂,提高所得微观结构的准确性和真实性;实现了M50合金在微-介观尺度上的模型建立,在原子尺度上实现了多物相的掺杂,同时又可以模拟合金在介观尺度上的各项性能,初步实现了多相合金M50的多尺度建模。
附图说明
图1为本发明对轴承钢材料M50合金多尺度建模的流程示意图。
具体实施方式
下面结合附图对本发明进行详细阐述。
参考图1,一种轴承钢材料M50的多尺度建模计算方法,属于合金多尺度建模领域;所述多尺度建模方法包括以下步骤:
步骤一:对M50合金棒料进行试验样品制备,完成热处理;
步骤二:将热处理后试样进行预处理,完成试验测试前准备;具体通过以下步骤实现:使用180#,400#,600#,800#,1000#,1500#和2000#的砂纸将M50合金试样表面进行逐级打磨,使用丝绒配W7金刚石抛光剂完成粗抛光,使用呢绒配W3.5金刚石抛光剂精抛光,保证表面粗糙度Ra≤0.1μm,之后依次用丙酮和无水乙醇超声清洗,超声清洗时间均为10min,吹干。
步骤三:使用显微硬度计或纳米压痕仪对试样进行观测区域划定;具体通过以下步骤实现:先用记号笔在试样表面中心圈出计划观测区域,在侧面划线作为试样朝向定位。随后用镊子将试样放于工作台上并调整至测量位置,使用40X物镜,转动升降手轮使得在目镜中可以观察到清晰的试样表面,转动调节手轮使得测微计零线对齐。选择载荷0.49N,加载时间30s,测量HK硬度,测量前Clear测微计清零,按下开始按钮,转动台自动切换至压头,完成硬度测试。卸载后,转动台自动切换至物镜,微调手轮,在目镜中确认压痕清晰完整,观测区域划定完成。
步骤四:使用X射线衍射技术(X-ray diffraction,XRD)对材料试样进行物相鉴定(马氏体、残余奥氏体、碳化物等),获得材料试样的物相信息;使用Bruker D8 Advance型多晶衍射仪。采用Cu-Kα辐射,全自动弯(平)晶石墨单色器,管电压35kV,管电流200mA,扫描步长0.02°,扫描速度1°/min,测量2θ角度范围0~-90°。
步骤五:使用电子背散射衍射技术(Electron Backscattered Diffraction,EBSD)获得划定区域试样的物相分布、尺寸以及取向等信息;使用蔡司场发射扫描电镜GeminiSEM 500和Oxford公司的EBSD设备,具体通过以下步骤实现:戴一次性手套,按标记朝向将试样装台,用导电胶带将试样固定在样品台上,将样品台放于观察室,合上封闭门。按观察室真空键,抽真空过程2分30秒,达到要求时,真空指示灯绿,可以进行试验。首先使用SEM观察,寻找显微硬度计标记点,拍摄出其在记号笔圈出区域的空间位置以方便后续试验寻找。根据XRD获得的物相种类对M50合金的物相组织进行寻找,确定第二相形状、尺寸、数量及取向信息。
步骤六:使用电子探针显微分析仪(Electron probe micro-analyzer,EPMA)对划定区域内每种物相进行元素定性定量分析,面扫获得试样介观尺度元素质量配比,点扫直接获得奥氏体、马氏体的元素组成,同时结合文献确定各类碳化物的种类和元素组成;
步骤七:使用Materials Studio软件建立每种物相的单胞模型;通过EBSD、EPMA获得M50合金的物相种类和元素组成,结合XRD获得物相空间点阵类型,使用MaterialsStudio软件建立每种物相的单胞模型,适当扩胞建成单一物相的原子尺度模型;
步骤八:通过介观尺度元素组成比例和不同物相的元素组成建立等式求得合金中每种物相的数量和M50合金微观建模适当尺寸即Fe原子数值;
步骤九:在铁基体中根据所得到的物相数据比例随机插入几种物相组织,得到轴承钢材料M50合金的微-介观结构模型。
步骤十:本发明所建立的轴承钢材料M50合金的微-介观结构模型,实现了在原子尺度上掺杂不同物相的介观结构模型建立,进一步可以在LAMMPS中模拟合金在介观尺度上的各项性能,仿真观察合金在热处理以及损伤时在原子尺度上的结构转变过程,研究M50合金的损伤机理,为主轴承故障检测研究提供理论支持。其中,Fe-C、Fe-Mo、Mo-Cr、V-Mo原子之间的相互作用力均用其对应的meam势函数描述,Fe-Cr原子之间的相互作用力用cdeam势函数描述,Fe-V原子之间的相互作用力用eam/fs势函数描述,C元素与Cr、Mo、V元素和Cr-V元素之间均采用lj/cut势函数进行描述,lj/cut势函数的参数均由通用力场下拟合得来。

Claims (9)

1.一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,包括以下步骤:
划定试样观测区域;
在所述观测区域对试样进行物相鉴定,获得试样含有的物相种类和空间点阵类型;
根据物相种类信息,获取所述观测区域试样的物相分布、尺寸以及取向等信息;
对所述观测区域内每种组织进行元素定性定量分析,得到材料介观尺度元素质量配比以及不同物相元素组成;
根据物相种类及其元素组成,使用Materials Studio软件建立每种物相组织的单胞模型,适当扩胞建成单一物相的原子尺度模型;
通过介观尺度元素质量配比和不同物相的元素组成建立等式,来求解原子模型中每种物相的数量和铁基体的适当尺寸即Fe原子数;
根据所求得的物相数量,在铁基体中随机插入几种物相组织,并优化结构,得到轴承钢材料M50合金的微-介观结构模型;
基于所述轴承钢材料M50合金的微-介观结构模型在原子尺度上模拟合金的各项性能,并仿真观察合金在热处理以及损伤时在原子尺度上的结构转变过程以及M50合金的损伤机理。
2.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,所述试样通过以下步骤制备:
使用中走丝线切割方法制备材料试样,进行一次淬火和三次回火热处理;
对热处理后试样进行预处理,将试样打磨抛光,完成测试前准备。
3.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,适用显微硬度计或纳米压痕仪完成观测区域划分,先手动使用记号笔圈出观测区域,使用记号笔在试样侧面做出标记,作为试样朝向的定位,然后根据划区精度要求和构建模型需要选择仪器,大区域划分使用显微硬度计设定载荷在试样上完成一次硬度测量,小区域划分使用纳米压痕仪,完成观测区域划分。
4.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,观测区域试样的物相种类和空间点阵类型鉴定采用XRD方法获取。
5.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,划定区域试样的物相分布、尺寸以及取向等信息采用EBSD方法获取。
6.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,划定区域内的介观尺度元素质量配比和每种物相的元素定性定量分析采用EPMA方法获取。
7.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,在Materials Studio软件建立每种组织的单胞模型时,其中,Fe-C、Fe-Mo、Mo-Cr、V-Mo原子之间的相互作用力均用其对应的meam势函数描述,Fe-Cr原子之间的相互作用力用cdeam势函数描述,Fe-V原子之间的相互作用力用eam/fs势函数描述,C元素与Cr、Mo、V元素和Cr-V元素之间均采用lj/cut势函数进行描述,lj/cut势函数的参数均由通用力场下拟合得来。
8.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,获得每种物相数量的同时获得材料微-介观结构模型的铁基体尺寸,同时能够保证M50合金组成元素的质量配比。
9.根据权利要求1所述一种轴承钢材料M50合金的多尺度建模计算方法,其特征在于,在铁基体中根据所得到的物相比例随机插入几种物相组织,并优化结构,得到轴承钢材料M50合金的微-介观结构模型,同时保证结构稳定且与实际观测到的组织一致。
CN202210461687.3A 2022-04-28 2022-04-28 一种轴承钢材料m50合金的多尺度建模计算方法 Active CN114757047B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210461687.3A CN114757047B (zh) 2022-04-28 2022-04-28 一种轴承钢材料m50合金的多尺度建模计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210461687.3A CN114757047B (zh) 2022-04-28 2022-04-28 一种轴承钢材料m50合金的多尺度建模计算方法

Publications (2)

Publication Number Publication Date
CN114757047A CN114757047A (zh) 2022-07-15
CN114757047B true CN114757047B (zh) 2024-04-02

Family

ID=82332417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210461687.3A Active CN114757047B (zh) 2022-04-28 2022-04-28 一种轴承钢材料m50合金的多尺度建模计算方法

Country Status (1)

Country Link
CN (1) CN114757047B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109002662A (zh) * 2018-09-12 2018-12-14 华北水利水电大学 一种多尺度参数耦合的圆弧粗糙表面建模方法
CN109191582A (zh) * 2018-09-12 2019-01-11 华北水利水电大学 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法
AU2020101597A4 (en) * 2020-07-31 2020-09-10 Donghua University A visualization method of ship structural parts based on a multi-scale model
AU2020104167A4 (en) * 2020-12-18 2021-03-04 Ningbo University A Preparation Method for Low-stress Electrode Material Sintering Based on the Mesoscale Model
CN112487584A (zh) * 2020-12-03 2021-03-12 天津工业大学 一种基于动力学的滚动轴承数字孪生建模方法
CN113096749A (zh) * 2021-06-10 2021-07-09 武汉大学深圳研究院 n型共掺杂金刚石半导体材料制备的多尺度耦合仿真方法
CN114038521A (zh) * 2021-11-11 2022-02-11 上海理工大学 一种锡基轴承合金力学性能的跨尺度仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110210103B (zh) * 2019-05-27 2020-11-27 北京工业大学 一种多相复合材料力学行为的多尺度模拟方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109002662A (zh) * 2018-09-12 2018-12-14 华北水利水电大学 一种多尺度参数耦合的圆弧粗糙表面建模方法
CN109191582A (zh) * 2018-09-12 2019-01-11 华北水利水电大学 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法
AU2020101597A4 (en) * 2020-07-31 2020-09-10 Donghua University A visualization method of ship structural parts based on a multi-scale model
CN112487584A (zh) * 2020-12-03 2021-03-12 天津工业大学 一种基于动力学的滚动轴承数字孪生建模方法
AU2020104167A4 (en) * 2020-12-18 2021-03-04 Ningbo University A Preparation Method for Low-stress Electrode Material Sintering Based on the Mesoscale Model
CN113096749A (zh) * 2021-06-10 2021-07-09 武汉大学深圳研究院 n型共掺杂金刚石半导体材料制备的多尺度耦合仿真方法
CN114038521A (zh) * 2021-11-11 2022-02-11 上海理工大学 一种锡基轴承合金力学性能的跨尺度仿真方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
泡沫铝弹性接触静力学行为的多尺度模拟;张元章;尹益辉;;计算力学学报;20081215(S1);全文 *
耐热不锈钢X8CrNi25-21介观尺度材料本构建模与仿真;周天丰;车江涛;鲁陈立;梁志强;王西彬;;北京理工大学学报;20170915(09);全文 *

Also Published As

Publication number Publication date
CN114757047A (zh) 2022-07-15

Similar Documents

Publication Publication Date Title
Guo et al. Effects of sliding amplitude and normal load on the fretting wear behavior of alloy 690 tube exposed to high temperature water
Wang et al. Effect of surface machining on the corrosion behaviour of 316 austenitic stainless steel in simulated PWR water
Kim et al. Evolution of carbides in cold-work tool steels
Lai et al. Effect of temperature on fretting wear behavior and mechanism of alloy 690 in water
Qin et al. Microstructural origin of the outstanding durability of the high nitrogen bearing steel X30CrMoN15-1
Wu et al. Evaluation of pre-shot peening on improvement of carburizing heat treatment of AISI 9310 gear steel
Jackson et al. The influence of microstructure on the ductile to brittle transition and fracture behaviour of HVOF NiCoCrAlY coatings determined via small punch tensile testing
Yuan et al. Tuning the pitting performance of a Cr-13 type martensitic stainless steel by tempering time
Yao et al. Effects of Ni and Cr on cryogenic impact toughness of bainite/martensite multiphase steels
Yang et al. Enhancing stress corrosion cracking resistance of machined surface via surface mechanical grinding treatment for AISI 316 L stainless steel
Kumar et al. Influence of ultrasonic shot peening on microstructure, mechanical, and electrochemical behavior of 316 stainless steel
CN114757047B (zh) 一种轴承钢材料m50合金的多尺度建模计算方法
Brodyanski et al. Microstructural changes caused by friction loading in high manganese TWIP steel and case-hardened 16MnCr5
Zhang et al. Altered microstructure characteristics and enhanced corrosion resistance of UNS S32750 duplex stainless steel via ultrasonic surface rolling process
Harris et al. Multiscale assessment of deformation induced by hydrogen environment-assisted cracking in a peak-aged Ni-Cu superalloy
Cui et al. Stability of retained austenite in martensitic high carbon steels. Part II: Mechanical stability
Yin et al. Influence of multiphase evolution on corrosion resistance of AlxCoCrFeNi alloys determined by transmission electron microscopy
Labonne et al. Precipitation kinetics in a Nb-stabilized ferritic stainless steel
Wanderka et al. Carbides in Co–Re–Cr-based high-temperature alloys
Zhang et al. The effect of bainite type on the evolution of carbide constituent during an accelerated aging in Cr-Mo-V steel
Osei et al. Effect of Al Additions on scale structure and oxidation kinetics of 430-ferritic stainless steel reheated in a combustion atmosphere
Li et al. Study of cold rolling on the transformation mechanism, microstructure, and properties of 304 austenitic stainless steel
Dobrotvorskii et al. Metallographic study of construction materials with atomic force microscopy method
Stygar et al. Physicochemical and mechanical properties of Crofer 22 APU ferritic steel applied in SOFC interconnects
CN111693674A (zh) Co基覆层高温摩擦磨损的组织性能分析方法

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