CN106501105A - 一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 - Google Patents
一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 Download PDFInfo
- Publication number
- CN106501105A CN106501105A CN201611266304.8A CN201611266304A CN106501105A CN 106501105 A CN106501105 A CN 106501105A CN 201611266304 A CN201611266304 A CN 201611266304A CN 106501105 A CN106501105 A CN 106501105A
- Authority
- CN
- China
- Prior art keywords
- mingled
- region
- joint
- crystal
- grain
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/32—Polishing; Etching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
Abstract
本发明涉及一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,对接头各区域截取微小试件,数据处理得到静载力学性能和循环载荷下的疲劳参数;对每个区域截取的微小试件依次进行磨削、抛光与电解;采用混合公式,结合接头各区的循环应力应变曲线,得到各区域夹杂的循环应力应变曲线;把不同的取向属性赋予代表性体积单元中的每个晶粒,建立搅拌摩擦焊接头微结构梯度模型,宏观模型的材料采用非线性运动硬化模型;对模型施加疲劳载荷。观察接头各区应力应变分布情况,重点分析代表性体积单元中多晶粒晶界处、不同取向的晶粒、夹杂附近的应力应变分布情况,以较为明显的局部应力集中区域和应变较大区域为疲劳薄弱区域。
Description
技术领域
本发明涉及微结构的搅拌摩擦焊接头的疲劳载荷下薄弱区域特征确定方法,特指考虑夹杂和晶体取向的焊接接头疲劳载荷下薄弱区域特征确定方法。
背景技术
从微观角度看材料是非均匀、各向异性的,晶粒的形态和晶体学取向会影响材料的力学性能。材料在疲劳载荷下,裂纹常常是从夹杂处或者晶界处萌生,然后扩展为宏观裂纹。研究表明材料的微观结构对疲劳微裂纹的萌生和扩展有重要影响。晶体塑性有限元是将晶体塑性力学本构关系与有限元方法相结合,从微观尺度出发模拟金属的塑性变形过程。
在搅拌摩擦焊接头中存在微观组织结构不同的区域,本发明的方法是考虑接头各区域晶体取向、晶界的作用及夹杂的影响,实现基于微结构的搅拌摩擦焊接头疲劳薄弱区域界定,较好地预测接头疲劳裂纹萌生部位,准确地分析局部微观结构的受力情况。
发明内容
本发明的目的在于能够考虑焊接接头中的夹杂和晶体取向在疲劳损伤过程中的影响,准确地反映焊接接头在疲劳载荷下的薄弱区域特征。
为解决上述技术问题,本发明采用的技术方案一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,该方法的实施流程如下,步骤1)对接头各区域截取微小试件,并对微小试件进行微拉伸实验和疲劳试验,数据处理得到静载力学性能和循环载荷下的疲劳参数;步骤2)对每个区域截取的微小试件依次进行磨削、抛光与电解,用电子背散射衍射(EBSD)技术观察接头各区域的晶体特征,结合金相图片观察夹杂特征,统计出各区域的晶粒平均尺寸、晶体取向分布及夹杂平均尺寸及分布位置情况;步骤3)对各区域夹杂和基体相进行纳米压痕实验测试其弹性模量和硬度,采用混合公式,结合接头各区的循环应力应变曲线,得到各区域夹杂的循环应力应变曲线;步骤4)选定足够表达晶体特征的一个尺寸,统计出这一尺寸下各个区域所包含的夹杂、晶粒个数及晶体取向分布;步骤5)根据接头各个区域的晶体形态和数量利用MATLAB软件生成不同大小、不同数量的voronoi多边形,将图形导入ABAQUS生成代表性体积单元,并参照金相图对代表性体积单元加入夹杂,接头的每个区域都生成相应的代表性体积单元,把晶体塑性本构关系赋予代表性体积单元的基体晶粒,其中的夹杂属性采用非线性运动硬化模型;步骤6)利用MATLAB生成随机晶体取向结合Python脚本语言创建材料库,把不同的取向属性赋予代表性体积单元中的每个晶粒,按比例为晶粒赋予特定织构的取向,其余晶粒赋予随机取向;步骤7)综合硬度和金相微观组织分区建立宏观模型,其形状和尺寸是疲劳试件的有效部分,在接头每个区域插入上述建立的代表性体积单元,建立搅拌摩擦焊接头微结构梯度模型,宏观模型的材料采用非线性运动硬化模型;步骤8)对模型施加疲劳载荷,根据接头各区的应力应变分布判定性能薄弱区域。观察接头各区应力应变分布情况,重点分析代表性体积单元中多晶粒晶界处、不同取向的晶粒、夹杂附近的应力应变分布情况,以较为明显的局部应力集中区域和应变较大区域为疲劳薄弱区域。
本发明所建方法考虑了焊接接头中各区域的微观结构即夹杂和晶体取向等微观特征对焊接接头的力学性能的影响,能够更准确的判断出接头疲劳薄弱区域并更精确地计算局部应力或应变数值,能够更真实的反映早期疲劳损伤过程,便于更准确的进行疲劳寿命预测,有效预防疲劳断裂的发生。
附图说明
图1为一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法流程示意图。
图2为微小试件形状及尺寸。
图3为接头热影响区的EBSD图。
图4为接头模型示意图。
图5为薄弱区域示意图。
具体实施方式
下面结合附图对本发明方法作进一步的描述,如图1所示一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,其步骤为:
步骤1)对接头焊核区、热机影响区、热影响区及母材各区域截取微小试件,试件厚度1mm,如图2,并对微小试件进行微拉伸实验和疲劳试验,数据处理得到静载力学性能和循环载荷下的疲劳参数;步骤2)对每个区域截取的微小试件依次进行磨削、抛光与电解,用电子背散射衍射(EBSD)技术观察接头各个区域的晶体特征,图3为接头热影响区的EBSD图,结合金相图片观察夹杂特征,统计出各区域的晶粒平均尺寸、晶体取向分布及夹杂平均尺寸及分布位置情况;步骤3)对各区域夹杂和基体相进行纳米压痕实验测试其弹性模量和硬度,采用混合公式,结合接头各区的循环应力应变曲线,得到各区域夹杂的循环应力应变曲线;步骤4)选定足够表达晶体特征的一个尺寸,统计出这一尺寸下各个区域所包含的夹杂、晶粒个数及晶体取向分布;步骤5)根据接头各区域的晶体形态和数量利用MATLAB软件生成不同大小、不同数量的voronoi多边形,经过CAD软件对voronoi多边形进行修正后再将图形导入ABAQUS生成代表性体积单元,并参照金相图对代表性体积单元加入夹杂,接头的每个区域都生成相应的代表性体积单元,如图4中的局部放大图为热影响区域的voronoi多边形,其中红色多边形为夹杂,通过子程序UMAT把晶体塑性本构关系赋予代表性体积单元的基体晶粒,夹杂属性采用非线性运动硬化模型;步骤6)利用MATLAB生成随机晶体取向结合Python脚本语言创建材料库,把不同的取向属性赋予代表性体积单元中的每个晶粒,按比例为晶粒赋予特定织构的取向,其余晶粒赋予随机取向;步骤7)综合硬度和金相微观组织分区建立宏观模型,其形状和尺寸是疲劳试件的有效部分,在接头每个区域插入上述建立的代表性体积单元,建立搅拌摩擦焊接头微结构梯度模型,宏观模型的材料采用非线性运动硬化模型,模型示意如图4;步骤8)对模型施加疲劳载荷,根据接头各区的应力应变分布判定性能薄弱区域。观察接头各区应力应变分布情况,重点分析代表性体积单元中多晶粒晶界处、不同取向的晶粒、夹杂附近的应力应变分布情况,以较为明显的局部应力集中区域和应变较大区域为疲劳薄弱区域。图5为模型中微结构薄弱区域示意图,薄弱区域为靠近夹杂的四晶界交界处。
Claims (4)
1.一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,其步骤为:
A)对接头焊核区、热机影响区、热影响区及母材各区域截取微小试件,并对微小试件进行微拉伸实验和疲劳试验,数据处理得到静载力学性能和循环载荷下的疲劳参数;
B)对每个区域截取的微小试件依次进行磨削、抛光与电解,用电子背散射衍射(EBSD)技术观察接头各个区域的晶体特征,结合金相图片观察夹杂特征,统计出各区域的晶粒平均尺寸、晶体取向分布及夹杂平均尺寸及分布位置情况;
C)对各区域夹杂和基体相进行纳米压痕实验测试其弹性模量和硬度,采用混合公式,结合接头各区的循环应力应变曲线,得到各区域夹杂的循环应力应变曲线;
D)选定足够表达晶体特征的一个尺寸,统计出这一尺寸下各个区域所包含的具体夹杂、晶粒个数及晶体取向分布;
E)根据接头各个区域的晶体形态和数量利用MATLAB软件生成不同大小、不同数量的voronoi多边形,经过CAD软件对voronoi多边形进行修正后再将图形导入ABAQUS生成代表性体积单元,并参照金相图对代表性体积单元加入夹杂,接头的每个区域都生成相应的代表性体积单元,把晶体塑性本构关系赋予代表性体积单元中的基体晶粒,其中的夹杂属性采用非线性运动硬化模型;
F)利用MATLAB生成随机晶体取向结合Python脚本语言创建材料库,把不同的取向属性赋予代表性体积单元中的每个晶粒,按比例为晶粒赋予特定织构的取向,其余晶粒赋予随机取向;
G)综合硬度和金相微观组织分区建立宏观模型,其形状和尺寸是疲劳试件的有效部分,在接头每个区域插入上述建立的代表性体积单元,建立搅拌摩擦焊接头微结构梯度模型,宏观模型的材料采用非线性运动硬化模型;
H)对模型施加疲劳载荷,根据接头各区的应力应变分布判定性能薄弱区域;观察接头各区应力应变分布情况,重点分析代表性体积单元中多晶粒晶界处、不同取向的晶粒、夹杂附近的应力应变分布情况,以明显的局部应力集中区域和应变较大区域为疲劳薄弱区域。
2.根据权利要求1所述的一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,其特征在于:所述步骤E)中针对接头各个区域晶粒形态、取向及夹杂分布生成能够反映真实晶粒的图形,并赋予夹杂和基体晶粒的材料属性和本构关系。
3.根据权利要求1所述的一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,其特征在于:所述步骤G)中建立的搅拌摩擦焊接头微结构梯度模型,首先从宏观上对接头进行分区,每个区域的材料属性根据实验得到的真实数据进行输入,其次在接头宏观模型中对每个区域都插入反映微观结构的含夹杂和晶粒特征的代表性体积单元,这样既能从宏观上看到接头各区的应力应变分布,又能从微观出发研究夹杂、晶界和晶体取向对材料应力应变分布的影响。
4.根据权利要求1所述的一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法,其特征在于:所述步骤H)中判定疲劳薄弱区域的方法,综合考虑应力、应变、晶粒特征、晶界形式及夹杂的尺寸与分布因素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611266304.8A CN106501105B (zh) | 2016-12-31 | 2016-12-31 | 一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611266304.8A CN106501105B (zh) | 2016-12-31 | 2016-12-31 | 一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106501105A true CN106501105A (zh) | 2017-03-15 |
CN106501105B CN106501105B (zh) | 2019-02-12 |
Family
ID=58334917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611266304.8A Active CN106501105B (zh) | 2016-12-31 | 2016-12-31 | 一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106501105B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181190A (zh) * | 2017-12-26 | 2018-06-19 | 大连交通大学 | 一种异种材料点焊接头疲劳极限快速预测方法 |
CN109271713A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 考虑晶体微结构力学的齿轮接触疲劳分析方法 |
CN109425630A (zh) * | 2017-08-21 | 2019-03-05 | 深圳市爱能森科技有限公司 | 一种快速寻找四元熔盐体系最低熔点的方法 |
CN109738620A (zh) * | 2018-12-18 | 2019-05-10 | 吉林大学 | 皮质骨径向表面裂纹扩展路径的预测方法 |
CN110118696A (zh) * | 2019-05-20 | 2019-08-13 | 吉林大学 | 纳米压痕结合疲劳载荷诱导Zr基非晶合金纳米晶化的方法 |
CN110222442A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 面心立方材料疲劳过程晶体塑性本构模型建立方法 |
CN110222439A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 基于Abaqus平台疲劳损伤与寿命评估方法 |
CN110231238A (zh) * | 2019-05-30 | 2019-09-13 | 吉林省维尔特隧道装备有限公司 | 掘进机滚刀刀圈截面的梯度硬度检测装置及方法 |
CN110320117A (zh) * | 2019-07-15 | 2019-10-11 | 泰州市金鹰精密铸造有限公司 | 基于单晶叶片蠕变性能的测试方法 |
CN110487788A (zh) * | 2019-07-12 | 2019-11-22 | 中国科学院金属研究所 | 一种单晶高温合金小角度晶界形成倾向性的评定方法 |
CN110907272A (zh) * | 2018-09-14 | 2020-03-24 | 西门子股份公司 | 3d打印元件的寿命预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073162A2 (en) * | 2001-03-07 | 2002-09-19 | Massachusetts Institute Of Technology | Systems and methods for estimation and analysis of mechanical property data |
CN103350279A (zh) * | 2013-07-30 | 2013-10-16 | 上海交通大学 | 搅拌摩擦焊接头弱连接缺陷制备方法 |
CN103970969A (zh) * | 2014-05-30 | 2014-08-06 | 哈尔滨工业大学 | 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法 |
CN104573392A (zh) * | 2015-01-27 | 2015-04-29 | 湖南大学 | 一种焊点疲劳寿命预测方法 |
-
2016
- 2016-12-31 CN CN201611266304.8A patent/CN106501105B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073162A2 (en) * | 2001-03-07 | 2002-09-19 | Massachusetts Institute Of Technology | Systems and methods for estimation and analysis of mechanical property data |
CN103350279A (zh) * | 2013-07-30 | 2013-10-16 | 上海交通大学 | 搅拌摩擦焊接头弱连接缺陷制备方法 |
CN103970969A (zh) * | 2014-05-30 | 2014-08-06 | 哈尔滨工业大学 | 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法 |
CN104573392A (zh) * | 2015-01-27 | 2015-04-29 | 湖南大学 | 一种焊点疲劳寿命预测方法 |
Non-Patent Citations (3)
Title |
---|
M.OHATA, ET AL.: "3D-Simulation of ductile failure in two-phase structural steel with heterogeneous microstructure", 《ENGINEERING FRACTURE MECHANICS》 * |
丛述玲等: "焊接接头细观多晶模型及拉伸应力变化数值模拟", 《兵器材料科学与工程》 * |
孙国芹等: "2219-T6铝合金搅拌摩擦焊接头疲劳性能", 《北京工业大学学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109425630A (zh) * | 2017-08-21 | 2019-03-05 | 深圳市爱能森科技有限公司 | 一种快速寻找四元熔盐体系最低熔点的方法 |
CN108181190A (zh) * | 2017-12-26 | 2018-06-19 | 大连交通大学 | 一种异种材料点焊接头疲劳极限快速预测方法 |
CN110907272B (zh) * | 2018-09-14 | 2022-06-24 | 西门子股份公司 | 3d打印元件的寿命预测方法 |
CN110907272A (zh) * | 2018-09-14 | 2020-03-24 | 西门子股份公司 | 3d打印元件的寿命预测方法 |
CN109271713A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 考虑晶体微结构力学的齿轮接触疲劳分析方法 |
CN109271713B (zh) * | 2018-09-25 | 2023-02-28 | 重庆大学 | 考虑晶体微结构力学的齿轮接触疲劳分析方法 |
CN109738620A (zh) * | 2018-12-18 | 2019-05-10 | 吉林大学 | 皮质骨径向表面裂纹扩展路径的预测方法 |
CN110118696A (zh) * | 2019-05-20 | 2019-08-13 | 吉林大学 | 纳米压痕结合疲劳载荷诱导Zr基非晶合金纳米晶化的方法 |
CN110118696B (zh) * | 2019-05-20 | 2021-07-23 | 吉林大学 | 纳米压痕结合疲劳载荷诱导Zr基非晶合金纳米晶化的方法 |
CN110231238B (zh) * | 2019-05-30 | 2020-04-24 | 吉林省维尔特隧道装备有限公司 | 掘进机滚刀刀圈截面的梯度硬度检测装置及方法 |
CN110231238A (zh) * | 2019-05-30 | 2019-09-13 | 吉林省维尔特隧道装备有限公司 | 掘进机滚刀刀圈截面的梯度硬度检测装置及方法 |
CN110222439A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 基于Abaqus平台疲劳损伤与寿命评估方法 |
CN110222442A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 面心立方材料疲劳过程晶体塑性本构模型建立方法 |
CN110487788A (zh) * | 2019-07-12 | 2019-11-22 | 中国科学院金属研究所 | 一种单晶高温合金小角度晶界形成倾向性的评定方法 |
CN110320117A (zh) * | 2019-07-15 | 2019-10-11 | 泰州市金鹰精密铸造有限公司 | 基于单晶叶片蠕变性能的测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106501105B (zh) | 2019-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106501105A (zh) | 一种考虑微结构夹杂和晶体取向的搅拌摩擦焊接头疲劳薄弱区域的确定方法 | |
Pilchak | A simple model to account for the role of microtexture on fatigue and dwell fatigue lifetimes of titanium alloys | |
CN106383052B (zh) | 一种考虑夹杂的金属疲劳载荷下薄弱区域确定方法 | |
Schweizer et al. | Mechanisms and modelling of fatigue crack growth under combined low and high cycle fatigue loading | |
Guilhem et al. | Investigation of the effect of grain clusters on fatigue crack initiation in polycrystals | |
Prawoto et al. | Computational approach using Johnson–Cook model on dual phase steel | |
Mlikota et al. | Numerical determination of Paris law constants for carbon steel using a two-scale model | |
Vöse et al. | An approach to life prediction for a nickel-base superalloy under isothermal and thermo-mechanical loading conditions | |
Wang et al. | HCF strength estimation of notched Ti–6Al–4V specimens considering the critical distance size effect | |
Nikara et al. | Finite element simulation of the micromechanical changes of the tissue and cells of potato response to impact test during storage by scanning electron microscopy | |
Park et al. | Effect of specimen thickness on fatigue crack growth rate | |
Li et al. | Fatigue cracking and fracture behaviors of coarse-grained copper under cyclic tension–compression and torsion loadings | |
Needleman | Effect of size on necking of dynamically loaded notched bars | |
Subrahmanya Prasad et al. | Numerical simulations of cylindrical void growth in Mg single crystals | |
Sun et al. | Crystal plastic modeling on fatigue properties for aluminum alloy friction stir welded joint | |
Boljanović et al. | Mixed mode crack growth simulation with/without overloads | |
Nyounguè et al. | Fracture characterisation of float glass under static and dynamic loading | |
Mikkola et al. | Mesoscale modelling of crack nucleation from defects in steel | |
Vuherer et al. | Behaviour of small cracks during their propagation from Vickers indentations in coarse-grain steel: an experimental investigation | |
Hammouda et al. | Significance of crack tip plasticity to early notch fatigue crack growth | |
Liu | Numerical analysis for effects of shot peening on fatigue crack growth | |
Radaj et al. | Elastic-plastic fatigue crack growth | |
Zerbst et al. | Fracture mechanics based determination of the fatigue strength of weldments | |
Chakraborty et al. | A unified cohesive zone approach to model the ductile to brittle transition of fracture toughness in reactor pressure vessel steels | |
Cornet et al. | A study of cyclic behaviour of a nickel-based superalloy at elevated temperature using a viscoplastic-damage model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |