CN104792808B - 一种金属基复合材料微观残余应力的检测方法 - Google Patents
一种金属基复合材料微观残余应力的检测方法 Download PDFInfo
- Publication number
- CN104792808B CN104792808B CN201410717568.5A CN201410717568A CN104792808B CN 104792808 B CN104792808 B CN 104792808B CN 201410717568 A CN201410717568 A CN 201410717568A CN 104792808 B CN104792808 B CN 104792808B
- Authority
- CN
- China
- Prior art keywords
- situ
- value
- detection
- sub
- residual stress
- 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.)
- Expired - Fee Related
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种金属基复合材料微观残余应力的检测方法,属于复合材料微观表征技术领域。所述方法以EBSD理论为基础,结合材料力学理论,建立微观区域残余应力与原位取向差参数之间的关系,先获取各检测点的位置数据、欧拉角数据和原位取向差数据;然后建立X‑Y坐标系,表征整体检测区域原位取向差分布函数;最后再采用公式σ=(cosθ‑1)E表征整体检测区域的残余应力分布函数,从而通过对EBSD参数的数据处理,最终定量表征金属基复合材料的微观残余应力分布。
Description
技术领域
本发明涉及一种金属基复合材料微观残余应力的检测方法,特别涉及一种采用电子背散射衍射(Electron Backscattered Diffraction,简称EBSD)表征金属基复合材料微观残余应力的方法,属于复合材料微观表征技术领域。
背景技术
对比金属材料,金属基复合材料具有低密度、高比强度、高比模量、低热膨胀系数等特点,因此成为各类多芯片组件和大电流功率模块理想的基底材料,以及各种质量小、强度高的零配件(如航空零件)的理想材料。由于热膨胀系数差异等因素,金属基复合材料在制备过程中形成微观残余应力,特别是材料内不同相界面处的微观残余应力是影响金属基复合材料力学性能的主要因素,如何有效表征界面处微观残余应力的是提高金属基复合材料力学性能的关键。
目前表征材料残余应力的方法包括钻孔法、环芯法、XRD表征法和拉曼表征法:钻孔法、环芯法仅满足宏观残余应力定性表征要求,无法满足结构复杂的多相材料的研究需求;XRD表征法的分辨尺寸大于EBSD分辨尺寸,且X射线穿深大于EBSD,导致该方法无法用于微米级区域残余应力研究;拉曼表征法则受到拉曼效应适用范围的影响,对大量无拉曼效应的材料,特别是无拉曼效应的金属材料无法进行表征。
目前EBSD技术主要用于表征单一金属材料的取向分布,特别集中于经轧制后金属织构取向表征方面,但目前EBSD无法用于定量表征微观残余应力,其主要原因在于:取向表征结果仅用于直观定性描述织构分布,即通过表征区域内具有特定取向的晶粒所占百分比来描述织构分布情况,尚未建立起取向结果与残余应力的关系,故无法将EBSD相比于其他检测手段的优点用于标识微观区域的残余应力。
发明内容
针对现有应力表征技术无法满足金属基复合材料等具有复杂结构的材料的残余应力检测需求的缺陷,本发明的目的在于提供一种金属基复合材料微观残余应力的表征方法,该方法以EBSD理论为基础,建立EBSD参数与残余应力的关系,从而通过对EBSD参数的数据处理,最终定量表征金属基复合材料的微观残余应力分布。
本发明所表征的材料为金属基复合材料,即材料包含金属相与非金属相,包括颗粒增强型金属基复合材料、纤维增强型金属基复合材料、三维网络结构增强型金属基复合材料等,不限于以上所述。
本发明的目的由以下技术方案实现:
一种金属基复合材料微观残余应力的检测方法,所述方法理论如下:
设待检测试样表面用于表征残余应力的一个晶粒为中心晶粒,其晶体坐标系(x,y,z)为参考坐标系,设主应力为σ1,σ2,σ3,其中σ3为垂直于待检测试样表面的方向,由于待检测试样表面为平面应力状态,所以σ3=0,σ1和σ2在待检测试样表面内,设待检测试样表面为x’y’平面,所以晶轴为z’轴,形成理论坐标系(x’,y’,z’);采用电子背散射衍射EBSD技术中的角轴对法表征晶体取向,则r=(0,0,1),θ=LMn,其中,r为旋转轴,θ为原位取向差,代表旋转角,参考坐标系沿旋转轴r旋转θ角度后与理论坐标系重合;LMn(LocalMisorientation)为检测区域内编号为n的中心晶粒的原位取向差,n代表中心晶粒在检测区域内的编号;
设施加载荷的方向与参考方向的角度差为φ,沿施加载荷的方向某段材料变形前长度为D0φ,该变形段起点与终点分别为(x1,y1),(x2,y2);变形后该段长度为D1φ,起点与终点分别为(x11,y11),(x21,y21),设检测方向与施加载荷的方向的角度差为即
则有,
根据应变定义,在φ方向上的应变εφ为:
设检测材料为各项同性材料,即不同方向弹性模量E相同,则应力为:
设即检测方向与施加载荷方向的角度差为0,则
σ=(cosθ-1)E (6)
其中σ为待检测区域微观残余应力,E为待检测相的弹性模量;
在以上理论的基础上,所述方法步骤如下:
第一步:选用包含EBSD系统的场发射电镜,在待检测试样表面收集各检测点的位置数据、欧拉角数据,并利用EBSD系统推导出在0.5°~5°区间内的原位取向差分布数据;
第二步:检测结束后,采用X-Y坐标系表征检测结果,设检测方向为X轴,沿X轴划分若干等长度子区域,每个子区域在X轴上的长度为0.5~5μm,
针对每个子区域,对所有检测点在X轴的取值采用唯一的X值表征,所有检测点在X轴的取值求平均值,即获得唯一的X值;
每个子区域的原位取向差分布数据表示为(fi,i),i为子区域内任意检测点的原位取向差值,i取0.5°~5°;fi是原位取向差等于i值的检测点数量占该子区域内所有检测点数量的百分比;
对各子区域内所有检测点的原位取向差值采用唯一Y值表征,表征方法如下:
Y值即公式(6)中的θ值;使各个子区域拥有唯一的X值和Y值以表征该子区域的原位取向差,集合各子区域的X值和Y值,拟合得出用于表征整体检测区域原位取向差分布函数;
第三步:采用公式σ=(cosθ-1)E将原位取向差值转化为微观残余应力值,从而将整体检测区域的原位取向差分布函数转化为表征整体检测区域的残余应力分布函数。
有益效果
本发明与现有技术相比,通过EBSD检测所得的原位取向差参数为已知量,以EBSD理论为基础,结合材料力学理论,建立微观区域残余应力与原位取向差参数之间的关系,从而通过EBSD检测定量表征金属基复合材料微观区域的残余应力分布,特别适用于无法采用现有方法表征微观区域残余应力分布的金属基复合材料。本发明适用范围不限于金属基复合材料,同样适用于多相导电材料的微观区域残余应力分布的研究。
附图说明
图1为实施例1的总体表征区域Local Misorientation示意图。
图2为实施例1的各子区域的X值、Y值原始数据。
图3为实施例1整体检测区域的残余应力分布函数示意图。
图4为实施例2的总体表征区域Local Misorientation示意图。
图5为实施例2的各子区域的X值、Y值原始数据。
图6为实施例2整体检测区域的残余应力分布函数示意图。
具体实施方式
本发明所表征的材料为金属基复合材料,即材料包含金属相与非金属相,包括颗粒增强型金属基复合材料、纤维增强型金属基复合材料、三维网络结构增强型金属基复合材料等,不限于以上所述。
本发明的目的由以下技术方案实现:
一种金属基复合材料微观残余应力的检测方法,所述方法理论如下:
设待检测试样表面用于表征残余应力的一个晶粒为中心晶粒,其晶体坐标系(x,y,z)为参考坐标系,设主应力为σ1,σ2,σ3,其中σ3为垂直于待检测试样表面的方向,由于待检测试样表面为平面应力状态,所以σ3=0,σ1和σ2在待检测试样表面内,设待检测试样表面为x’y’平面,所以晶轴为z’轴,形成理论坐标系(x’,y’,z’);采用电子背散射衍射EBSD技术中的角轴对法表征晶体取向,则r=(0,0,1),θ=LMn,其中,r为旋转轴,θ为原位取向差,代表旋转角,参考坐标系沿旋转轴r旋转θ角度后与理论坐标系重合;LMn(LocalMisorientation)为检测区域内编号为n的中心晶粒的原位取向差,n代表中心晶粒在检测区域内的编号;
设施加载荷的方向与参考方向的角度差为φ,沿施加载荷的方向某段材料变形前长度为D0φ,该变形段起点与终点分别为(x1,y1),(x2,y2);变形后该段长度为D1φ,起点与终点分别为(x11,y11),(x21,y21),设检测方向与施加载荷的方向的角度差为即
则有,
根据应变定义,在φ方向上的应变εφ为:
设检测材料为各项同性材料,即不同方向弹性模量E相同,则应力为:
设即检测方向与施加载荷方向的角度差为0,则
σ=(cosθ-1)E (6)
其中σ为待检测区域微观残余应力,E为待检测相的弹性模量;
在以上理论的基础上,所述方法步骤如下:
第一步:选用包含EBSD系统的场发射电镜,在待检测试样表面收集各检测点的位置数据、欧拉角数据,并利用EBSD系统推导出在0.5°~5°区间内的原位取向差分布数据;
第二步:检测结束后,采用X-Y坐标系表征检测结果,设检测方向为X轴,沿X轴划分若干等长度子区域,每个子区域在X轴上的长度为0.5~5μm,
针对每个子区域,对所有检测点在X轴的取值采用唯一的X值表征,所有检测点在X轴的取值求平均值,即获得唯一的X值;
每个子区域的原位取向差分布数据表示为(fi,i),i为子区域内任意检测点的原位取向差值,i取0.5°~5°;fi是原位取向差等于i值的检测点数量占该子区域内所有检测点数量的百分比;
对各子区域内所有检测点的原位取向差值采用唯一Y值表征,表征方法如下:
Y值即公式(6)中的θ值;使各个子区域拥有唯一的X值和Y值以表征该子区域的原位取向差,集合各子区域的X值和Y值,拟合得出用于表征整体检测区域原位取向差分布函数;
第三步:采用公式σ=(cosθ-1)E将原位取向差值转化为微观残余应力值,从而将整体检测区域的原位取向差分布函数转化为表征整体检测区域的残余应力分布函数。
下面结合附图和具体实施例来详述本发明,但不限于此。
实施例1
检测材料为三维网络增强SiC3D/Al金属基复合材料,检测流程如下:采用热场发射电镜,配合EBSD系统,检测区域选取SiC/Al界面,推导在0.5°~5°区间内的原位取向差分布数据;如图1所示。选取垂直于SiC/Al界面的方向为检测方向,即X轴方向,划分子区域,子区域宽度为1μm,
针对每个子区域,对所有检测点在X轴的取值采用唯一的X值表征,所有检测点在X轴的取值求平均值,即获得唯一的X值;
每个子区域的原位取向差分布数据表示为(fi,i),i为子区域内任意检测点的原位取向差值,i取0.5°~5°;fi是原位取向差等于i值的检测点数量占该子区域内所有检测点数量的百分比;如图2所示;对各子区域内所有检测点的原位取向差值采用唯一Y值表征,表征方法如下:
集合各子区域的X值和Y值,利用最小二乘法进行拟合,建立X和Y值的函数关系,如图3所示。通过公式(6)计算,Y值即公式(6)中的θ值,得出结论:最大残余应力值σ=(cosθ-1)E=22.96MPa。
实施例2
检测材料为颗粒增强SiCp/Cu金属基复合材料,检测流程如下:采用冷场发射电镜,配合EBSD系统,检测区域选取SiC/Al界面,推导在0.5°~5°区间内的原位取向差分布数据;如图4所示。选取垂直于SiC/Cu界面的方向为检测方向,即X轴方向,划分子区域,子区域宽度为1μm,
针对每个子区域,对所有检测点在X轴的取值采用唯一的X值表征,所有检测点在X轴的取值求平均值,即获得唯一的X值;
每个子区域的原位取向差分布数据表示为(fi,i),i为子区域内任意检测点的原位取向差值,i取0.5°~5°;fi是原位取向差等于i值的检测点数量占该子区域内所有检测点数量的百分比;如图5所示;对各子区域内所有检测点的原位取向差值采用唯一Y值表征,表征方法如下:
集合各子区域的X值和Y值,利用最小二乘法进行拟合,建立X和Y值的函数关系,如图6所示。通过公式(6)计算,Y值即公式(6)中的θ值,得出结论:最大残余应力值σ=(cosθ-1)E=192.12MPa。
本发明包括但不限于以上实施例,凡是在本发明精神的原则之下进行的任何等同替换或局部改进,都将视为在本发明的保护范围之内。
Claims (3)
1.一种金属基复合材料微观残余应力的检测方法,其特征在于,
设待检测试样表面用于表征残余应力的一个晶粒为中心晶粒,其晶体坐标系(x,y,z)为参考坐标系,设主应力为σ1,σ2,σ3,其中σ3为垂直于待检测试样表面的方向,由于待检测试样表面为平面应力状态,所以σ3=0,σ1和σ2在待检测试样表面内,设待检测试样表面为x’y’平面,所以晶轴为z’轴,形成理论坐标系(x’,y’,z’);采用电子背散射衍射EBSD技术中的角轴对法表征晶体取向,则r=(0,0,1),θ=LMn,其中,r为旋转轴,θ为原位取向差,代表旋转角,参考坐标系沿旋转轴r旋转θ角度后与理论坐标系重合;LMn为检测区域内编号为n的中心晶粒的原位取向差,n代表中心晶粒在检测区域内的编号;
设施加载荷的方向与参考方向的角度差为φ,沿施加载荷的方向某段材料变形前长度为D0φ,该变形段起点与终点分别为(x1,y1),(x2,y2);变形后该段长度为D1φ,起点与终点分别为(x11,y11),(x21,y21),设检测方向与施加载荷的方向的角度差为即
则有,
根据应变定义,在φ方向上的应变εφ为:
设检测材料为各项同性材料,即不同方向弹性模量E相同,则应力为:
设即检测方向与施加载荷方向的角度差为0,则
σ=(cosθ-1)E (6)
其中σ为待检测区域微观残余应力,E为待检测相的弹性模量;
利用公式(6)的检测方法步骤如下:
第一步:选用包含EBSD系统的场发射电镜,在待检测试样表面收集各检测点的位置数据、欧拉角数据,并利用EBSD系统推导出在0.5°~5°区间内的原位取向差分布数据;
第二步:检测结束后,采用X-Y坐标系表征检测结果,设检测方向为X轴,沿X轴划分若干等长度子区域,每个子区域在X轴上的长度为0.5~5μm,
针对每个子区域,对所有检测点在X轴的取值采用唯一的X值表征,所有检测点在X轴的取值求平均值,即获得唯一的X值;
每个子区域的原位取向差分布数据表示为(fi,i),i为子区域内任意检测点的原位取向差值,i取0.5°~5°;fi是原位取向差等于i值的检测点数量占该子区域内所有检测点数量的百分比;
对各子区域内所有检测点的原位取向差值采用唯一Y值表征,表征方法如下:
Y值即公式(6)中的θ值;使各个子区域拥有唯一的X值和Y值以表征该子区域的原位取向差,集合各子区域的X值和Y值,拟合得出用于表征整体检测区域原位取向差分布函数;
第三步:采用公式σ=(cosθ-1)E将原位取向差值转化为微观残余应力值,从而将整体检测区域的原位取向差分布函数转化为表征整体检测区域的残余应力分布函数。
2.根据权利要求1所述的一种金属基复合材料微观残余应力的检测方法,其特征在于,所述金属基复合材料由金属相和非金属相组成。
3.根据权利要求2所述的一种金属基复合材料微观残余应力的检测方法,其特征在于,所述金属基复合材料为颗粒增强型金属基复合材料、纤维增强型金属基复合材料或三维网络结构增强型金属基复合材料。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410717568.5A CN104792808B (zh) | 2014-12-01 | 2014-12-01 | 一种金属基复合材料微观残余应力的检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410717568.5A CN104792808B (zh) | 2014-12-01 | 2014-12-01 | 一种金属基复合材料微观残余应力的检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104792808A CN104792808A (zh) | 2015-07-22 |
| CN104792808B true CN104792808B (zh) | 2017-06-06 |
Family
ID=53557791
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410717568.5A Expired - Fee Related CN104792808B (zh) | 2014-12-01 | 2014-12-01 | 一种金属基复合材料微观残余应力的检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104792808B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105136361A (zh) * | 2015-09-07 | 2015-12-09 | 哈尔滨工业大学 | 一种用x射线衍射测定立方单晶体材料应力的方法 |
| CN106546481B (zh) * | 2016-09-28 | 2019-01-11 | 同济大学 | 类岩石材料力学特性的测试方法 |
| CN107014539B (zh) * | 2017-06-05 | 2020-05-22 | 中国人民解放军装甲兵工程学院 | 一种多相复合涂层残余应力检测方法 |
| CN109870257B (zh) * | 2017-12-04 | 2020-12-18 | 有研工程技术研究院有限公司 | 一种板材厚度方向淬火残余应力分布预测方法 |
| JP2019124481A (ja) * | 2018-01-12 | 2019-07-25 | 株式会社神戸製鋼所 | 残留応力測定方法 |
| CN110333219B (zh) * | 2019-07-15 | 2020-07-07 | 天津大学 | 应力分量的解耦检测方法与装置 |
| CN112326084B (zh) * | 2020-10-14 | 2021-09-14 | 北京科技大学 | 一种利用x射线测量含织构材料残余应力的方法 |
| CN112100248B (zh) * | 2020-11-18 | 2021-02-26 | 树根互联技术有限公司 | 基于智能切割边界的应力数据统计方法、装置及电子设备 |
| KR102366284B1 (ko) * | 2020-12-28 | 2022-02-23 | 현대제철 주식회사 | 핫 스탬핑 부품 및 그 제조방법 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1793872A (zh) * | 2005-12-29 | 2006-06-28 | 哈尔滨工业大学 | 微小区域残余应力的无损检测方法 |
| CN102735703A (zh) * | 2012-07-05 | 2012-10-17 | 首钢总公司 | 采用ebsd定量评价钢中残余奥氏体的方法 |
| CN103837557A (zh) * | 2014-01-21 | 2014-06-04 | 首钢总公司 | 采用ebsd表征热轧钢板截面氧化铁皮微观结构的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013194246A (ja) * | 2012-03-15 | 2013-09-30 | Mitsubishi Shindoh Co Ltd | 残留応力の少ないリードフレーム用Cu−Cr−Sn系銅合金板 |
-
2014
- 2014-12-01 CN CN201410717568.5A patent/CN104792808B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1793872A (zh) * | 2005-12-29 | 2006-06-28 | 哈尔滨工业大学 | 微小区域残余应力的无损检测方法 |
| CN102735703A (zh) * | 2012-07-05 | 2012-10-17 | 首钢总公司 | 采用ebsd定量评价钢中残余奥氏体的方法 |
| CN103837557A (zh) * | 2014-01-21 | 2014-06-04 | 首钢总公司 | 采用ebsd表征热轧钢板截面氧化铁皮微观结构的方法 |
Non-Patent Citations (3)
| Title |
|---|
| Mark D. Vaudin等.Comparison of nanoscale measurements of strain and stress using.《Applied Physics Letters》.2008, * |
| 染佰峰等.颗粒增强复合材料残余应力.《材料导报》.2009, * |
| 黄亚敏等.基于电子背散射衍射(EBSD)技术的材料微区.《电子显微学报》.2010, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104792808A (zh) | 2015-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104792808A (zh) | 一种金属基复合材料微观残余应力的检测方法 | |
| An et al. | VULCAN: A “hammer” for high-temperature materials research | |
| Melitz et al. | Kelvin probe force microscopy and its application | |
| Hara et al. | Sintering analysis of sub-micron-sized nickel powders: Kinetic Monte Carlo simulation verified by FIB–SEM reconstruction | |
| Lee et al. | Cold-temperature deformation of nano-sized tungsten and niobium as revealed by in-situ nano-mechanical experiments | |
| LI et al. | A method to calculate the dislocation density of a TWIP steel based on neutron diffraction and synchrotron X-ray diffraction | |
| Balzar et al. | Analyzing microstructure by Rietveld refinement | |
| CN108844816A (zh) | 一种基于十字形试件双轴拉伸试验的屈服准则评估方法 | |
| Okoro et al. | Synchrotron-based measurement of the impact of thermal cycling on the evolution of stresses in Cu through-silicon vias | |
| CN104458373A (zh) | 一种适用于ebsd检测的金属基复合材料试样的制备方法 | |
| Byeon et al. | Diffusion Along Dislocations Mitigates Self‐Limiting Na Diffusion in Crystalline Sn | |
| Oni et al. | Large area strain analysis using scanning transmission electron microscopy across multiple images | |
| Mao et al. | Bending fracture behavior of freestanding (Gd0. 9Yb0. 1) 2Zr2O7 coatings by using digital image correlation and FEM simulation with 3D geometrical reconstruction | |
| Li et al. | Fabrication of heat-resistant grids and their application to deformation measurements using a sampling moiré method | |
| Liu et al. | Quartet structure generation set algorithm based 3D reconstruction on porous structures of sintered copper joints for power electronics packaging | |
| Kumar et al. | Microstructural, microhardness and electrical conductivity analysis of AD31T alloy processed by friction stir processing | |
| Zhao et al. | Experimental study of residual stress distribution in interfacial micro-region of SiCf/Ti17 composites via micro-slotting method | |
| Deka et al. | Localized strain evolution and accumulation at the notch tip of an aluminum single crystal under cyclic loading | |
| Zhang et al. | In-situ 3D deformation and damage evolution of C/SiC composites under high-temperature tensile loading: μCT and digital volume correlation analysis | |
| CN107014539B (zh) | 一种多相复合涂层残余应力检测方法 | |
| Nakajo et al. | Characterization of cracks and their effects on the effective transport pathways in Ni-YSZ anodes after reoxidation using X-ray nanotomography | |
| Maniatty et al. | Effect of microstructure on electromigration-induced stress | |
| CN112050978B (zh) | 一种机匣用Ti2AlNb合金的X射线残余应力测试方法 | |
| CN110501366A (zh) | 一种低维功能复合材料温度相关等效电学性能的预测方法 | |
| Huang et al. | A comparative study of the mechanical properties of sintered submicron and micron Cu |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170606 |