CN115165182A - 一种高效、高空间分辨率的残余应力测量方法 - Google Patents

一种高效、高空间分辨率的残余应力测量方法 Download PDF

Info

Publication number
CN115165182A
CN115165182A CN202210661052.8A CN202210661052A CN115165182A CN 115165182 A CN115165182 A CN 115165182A CN 202210661052 A CN202210661052 A CN 202210661052A CN 115165182 A CN115165182 A CN 115165182A
Authority
CN
China
Prior art keywords
strain
ring
residual stress
depth
diameter
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.)
Pending
Application number
CN202210661052.8A
Other languages
English (en)
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.)
Yangtze River Delta Advanced Materials Research Institute
Original Assignee
Yangtze River Delta Advanced Materials Research Institute
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 Yangtze River Delta Advanced Materials Research Institute filed Critical Yangtze River Delta Advanced Materials Research Institute
Priority to CN202210661052.8A priority Critical patent/CN115165182A/zh
Publication of CN115165182A publication Critical patent/CN115165182A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0047—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00—Technologies related to metal processing
    • Y02P10/20—Recycling

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Toxicology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种高效、高空间分辨率的残余应力测量方法，属于材料分析测试技术领域。本发明解决了现有常用材料残余应力测试方法很难实现微米尺度上的空间定位，阻碍了材料中“微观组织‑残余应力”关系构筑的技术问题。本发明方法：在待测区沉积散斑点以帮助识别；在环切位置内侧沉积保护环以校准和防反沉积；用聚焦离子束逐层切削环槽，并使用电子成像；用有限元软件计算切削带来的理论应变释放；使用数字图像关联技术计算实际切削过程中的应变释放；对比理论和实际应变释放，计算测试位置的残余应力。本发明具有高效、高空间分辨率、适用于大部分固体材料的特点，可测单晶、非晶、薄膜等X射线难测样品。

Description

一种高效、高空间分辨率的残余应力测量方法
技术领域
本发明涉及的是一种高空间分辨率的残余应力测量方法，具体地说，是一种使用聚焦离子束和数字图像关联技术的残余应力测量方法，属于材料分析测试技术领域。
背景技术
残余应力指无外载力或力矩时存在于材料内部的应力，其水平和分布会显著影响部件的服役性能，因此残余应力检测在评估和预测部件服役性能时至关重要。目前在工业和研究中对于材料残余应力的测试存在多种方法，包括以钻孔法、切割法为代表的有损方法和以超声法、电磁法、衍射法为代表的无损方法。前者一般通过机械加工部分或者完全破坏样品，采用应变片测量释放的弹性应变以计算材料中的应力。此类方法受限于机械加工精度和应变片尺寸，空间分辨率较差。无损方法通常要求无应力状态的样品作为测试的参考基准，且对样品类型有一定程度的限制，例如电磁法仅适用于铁磁性材料、衍射法仅适用于晶体材料等。随着近年来的发展，中子衍射的空间分辨率已经达到毫米尺度，同步辐射 X射线的高亮度和高准直性进一步将衍射法的空间分辨率推进至微米量级，但这些方法通常需要使用大型科学设备，严重限制了企业和科研单位使用的灵活性。实验室X射线的空间分辨率依赖相应孔径的准直器，提高空间分辨率往往以牺牲有效亮度、测试实验延长为代价。此外，一个共性问题在于现有常用材料残余应力测试方法很难实现微米尺度上的空间定位，阻碍了材料中“微观组织-残余应力”关系构筑，是高性能材料研发与设计亟需突破的关键问题之一。
因此迫切的需要开发一种可以高效、高空间分辨率的探测材料残余应力的方法，从而节约成本，提高生产效率，对理解材料的失效机制，改进部件设计具有重大意义。
发明内容
本发明针对如上所述的现有技术的不足和缺陷，提供了一种使用聚焦离子束切削环芯，并结合数字图像关联技术的残余应力测量方法，适用于高效、高空间分辨率的残余应力的测量，且适用于大部分固体材料，包括但不限于一般晶体和非晶材料、非铁磁性材料、单晶等。
为了解决上述技术问题，本发明采取了以下的技术方案：
本发明的一种高效、高空间分辨率的残余应力测量方法是通过以下技术方案实现的：(1)制标记点：在待测区沉积散斑点以帮助识别；(2)制保护环：在环切位置内侧沉积保护环以校准和防反沉积；(3)环切并成像：用聚焦离子束逐层切削环槽，并使用电子成像；(4)计算理论应变释放：用有限元软件计算切削带来的理论应变释放；(5)计算实际应变释放：使用数字图像关联技术计算实际切削过程中的应变释放；(6)应力计算：对比理论和实际应变释放，计算测试位置的残余应力。具体包括实验部分、数值计算部分和后处理部分：
A.实验部分
A.1准备待测样品(固体材料)，将其表面处理至没有划痕确保制备过程中不引入和释放应力，然后置于配有聚焦离子束的扫描电子显微镜中的样品台；
A.2使用电子束或离子束在待测区沉积指定占空比的、随机分布的散斑；
A.3在待切削环槽内环(环芯圆柱外环)使用电子束或离子束沉积一个直径等于环芯直径D、宽度大幅小于直径的Pt圆环；
A.4使用离子束进行深度为Δh，且Δh＜＜D的环切；
A.5使用电子束和离子束成像，并保存采取了倾转角度补偿的扫描电子显微图片，补偿角度与样品台相对于电子束的倾转角度一致；
A.6交替重复步骤A.4和A.5，直至环槽的总深度(h)与环芯直的比例h/D，暨相对深度，达到目标值；通常h/D大于1/2后圆柱内应变释放完全，环芯表面应变趋于稳定，继续切削不会带来变化，结束切削。
B.有限元计算部分
B.1在有限元计算软件中建立完整和聚焦离子束逐层切除环槽后的材料模型，每道切除的深度无需和实验过程一致；
B.2对模型施加任意已知等双轴应力，经平衡计算得到不同的环槽相对深度h/D下的环芯表面应变分布；
B.3取不同有效直径(有效直径Deff＝μD，μ为直径计算有效系数)对应的圆内应变分布，并分别计算随着环槽加深，应变释放的演化；
B.4根据应变释放得到应变释放关于相对深度h/D的函数SRμ(h/D)，将应变释放函数关于完整模型的应变进行归一，可以得到不受材料属性影响、受计算有效系数影响的通用归一化应变释放函数SRμ,norm.(h/D)；
B.5通过应变释放函数得到体平均残余应力释放系数γ；
B.6对归一化应变释放函数SRμ,norm.(h/D)求导，得到指定有效直径下的应变响应函数 SIμ(h/D)；
B.7通过实验测得的应变相应函数与步骤B.6中的理论应变响应函数在相同深度区间内的比值得到该深度区间内残余应力对环芯表面应变的影响。
C.后处理部分
C.1使用数字图像关联技术对扫描电子显微图像进行计算，得到环芯表面随相对深度增加而演化的应变场；
C.2取指定有效直径并沿指定方向计算相应圆内平均应变释放，关于相对深度作图并拟合，得到实验应变释放曲线；
C.3将实验应变释放曲线的释放平台，暨相对深度趋于无限时的应变释放带入步骤B.5中的理论关系进行计算，计算圆柱体积内该方向的平均残余应力，过程中保持有效直径一致；C.4根据实验应变释放曲线计算实验应变响应函数，与步骤B.6中的理论应变响应函数进行对比，可以得到任意两张扫描电子显微图像(SEM图像)对应深度区间内的平均残余应力，过程中保持有效直径一致；
C.5改变数字图像关联技术计算的应变方向并重复步骤C.1至C.4，得到不同方向的应力分量。
进一步地限定，在步骤A.1中，表面处理包括但不仅限于电解抛光和机械抛光等常见表面处理方法。
在上述步骤中，散斑点和保护环的沉积方式包括但不仅限于采用电子束或离子束沉积。
进一步地限定，在步骤A.2中，散斑点的占空比为15％～85％。
进一步地限定，在步骤A.2和A.3中，散斑点和保护环的厚度为0.05μm及以上。
进一步地限定，在步骤A.2和A.3中，沉沉积散斑点和保护环覆盖的范围是0.5μm～100μm。
进一步地限定，在步骤A.4中，聚焦离子束切削的内环直径范围为0.5μm～95μm，外环的直径范围为0.5μm～100μm。
进一步地限定，在步骤A.4中，切削步进值Δh为25nm～2μm。
进一步地限定，在步骤B.4中，应变释放函数SRμ(h/D)的计算方法为：
Figure BDA0003690881190000041
其中，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，δε为应变释放，r为环芯圆柱表面圆心为原点的圆柱极坐标系的径向坐标，ε(r,h/D)为圆柱表面内任意点的局部应变。
进一步地限定，在步骤B.5中，体平均残余应力释放系数γ的计算方法为：
Figure BDA0003690881190000042
其中，γ为体平均残余应力释放系数，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，σapp为数值计算时施加的已知等双轴应力。
进一步地限定，在步骤B.6中，指定有效直径下的应变响应函数SIμ(h/D)的计算方法为：
Figure BDA0003690881190000043
其中，h为环槽总深度，D为环芯圆柱直径，SRμ为计算有效系数为μ时应变释放关于相对深度h/D的函数。
与现有技术相比，本发明具有以下有益效果：
本发明所提供的应力测量技术具有高效、高空间分辨率、适用于大部分固体材料的特点，可测单晶、非晶、薄膜等X射线难测样品；可以应用于一般的晶体材料、强织构材料、单晶合金、焊接样品、增材产品和半导体材料等多种固体合金和其它固体上，无需要精确评估无应变晶格参数，且操作简单，数据处理耗时短，尤其对于企业来说，本发明可以快速且精准的得到部件的应力结果，从而可以有效地调整工艺，提高生产效率；针对研究目的，本发明方法可以灵活选择空间分辨率，精确得到不同方向的应力分量，且适用于绝大多数固体材料。
为了能够更进一步了解本发明的特征及技术内容，请参阅以下有关本发明详细说明与附图，然而所附的附图仅提供参考和说明之用，并非用来对本发明加以限制。
附图说明
图1是散斑点示例(a)和图环切过程(b至d)的扫描电子显微图像；
图2是仿真模型示意图；
图3是计算得到应变释放函；
图4是实验测得环芯表面应变释放结果例图。
具体实施方式
下面结合附图详细说明本发明的具体实施。以下实施方式仅用于说明本方法，而非对本发明的限制。说明中的方向性术语仅基于附图所示方位关系，意在简化和清晰化描述，不指示或暗示特定方位的必须性，故不能理解为对本发明的限制。说明中的具体参数和工艺，包括但不限于使用的商用软件、样品尺寸、工作距离、表面处理方法、测式区尺寸、电流、电压、放大倍数、切削深度等因具体设备而异，故只用于帮助描述说明本实施例，不能理解为对本发明的限制。
实施例1：
本实施例一种高效、高空间分辨率的残余应力测量方法是通过以下技术方案实现的：是通过以下技术方案实现的，具体包括实验部分、数值计算部分和后处理部分：
A.实验部分
A.1准备待测固体材料，将其表面处理至没有划痕并置于配有聚焦离子束的扫描电子显微镜中，确保制备过程中不引入和释放应力；
A.2使用离子束在待测区沉积指定占空比的、随机分布的散斑；
A.3在待切削环槽内环使用离子束沉积一个直径等于环芯直径D、宽度大幅小于直径的Pt圆环；
A.4将样品和样品台倾转至样品表面正对离子束，使用离子束进行深度为Δh，且Δh＜＜ D的环切；
A.5使用电子束和离子束成像，并保存采取了倾转角度补偿的扫描电子显微图片，补偿角度与步骤A.4中的倾转角度一致；
A.6交替重复步骤A.4和A.5，直至环槽的总深度(h)与环芯直的比例h/D，暨相对深度，达到目标值；通常h/D大于1/2后圆柱内应变释放完全，环芯表面应变趋于稳定，继续切削不会带来变化。
B.有限元计算部分
B.1在有限元计算软件中建立完整和聚焦离子束逐层切除环槽后的材料模型，每道切除的深度无需和实验过程一致；
B.2对模型施加任意已知等双轴应力，经平衡计算得到不同的环槽相对深度h/D下的环芯表面应变分布；
B.3取不同有效直径(有效直径Deff＝μD，μ为直径计算有效系数)对应的圆内应变分布，并分别计算随着环槽加深，应变释放的演化；
B.4根据应变释放得到应变释放函数，其计算方法分别为：
Figure BDA0003690881190000061
其中，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，δε为应变释放，r为以环芯表面圆心为原点的圆柱极坐标系的径向坐标，ε(r,h/D)为圆柱顶面内任意点的局部应变；将应变释放函数关于完整模型的应变进行归一，可以得到不受材料属性影响、受计算有效系数影响的通用归一化应变释放函数SRμ,norm.(h/D)；
B.5通过应变释放函数得到体平均残余应力释放系数，其计算方法为：
Figure BDA0003690881190000062
其中，γ为体平均残余应力释放系数，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，σapp为数值计算时施加的已知等双轴应力；
B.6对归一化应变释放函数SRμ,norm.(h/D)求导，得到指定有效直径下的应变响应函数，其计算方法为：
Figure BDA0003690881190000063
B.7通过实验测得的应变相应函数与步骤B.6中的理论应变响应函数在相同深度区间内的比值得到该深度区间内残余应力对环芯表面应变的影响。
C.后处理部分
C.1使用数字图像关联技术对扫描电子显微图像进行计算，得到环芯表面随相对深度增加而演化的应变场；
C.2取指定有效直径并沿指定方向计算相应圆内平均应变释放，关于相对深度作图并拟合，得到实验应变释放曲线；
C.3将实验应变释放曲线的释放平台，暨相对深度趋于无限时的应变释放带入步骤 B.5中的理论关系进行计算，计算圆柱体积内该方向的平均残余应力，过程中保持有效直径一致,具体计算方法为：
Figure BDA0003690881190000071
其中，
Figure BDA0003690881190000072
为环芯内沿l方向的平均残余应力，SRμ,exp(h/D)为计算有效系数为μ时的实验应变释放函数，h为环槽总深度，D为环芯圆柱直径；
C.4根据实验应变释放曲线计算实验应变响应函数，与步骤B.6中的理论应变响应函数进行对比，可以得到任意两张扫描电子显微图像对应深度区间内的平均残余应力。过程中保持有效直径一致；
C.5改变数字图像关联技术计算的应变方向并重复步骤C.1至C.4，得到不同方向的应力分量。
在步骤A.1中，表面处理采用电解抛光。
在步骤A.2中，散斑点的占空比为15％～85％。
在步骤A.2和A.3中，散斑点和保护环的厚度为0.05μm及以上。
在步骤A.2和A.3中，沉沉积散斑点和保护环覆盖的范围是0.5μm～100μm。
在步骤A.4中，聚焦离子束切削的内环直径范围为0.5μm～95μm，外环的直径范围为0.5μm～100μm。
在步骤A.4中，切削步进值Δh为25nm～2μm。
在步骤B.4中，应变释放函数SRμ(h/D)的计算方法为：
Figure BDA0003690881190000073
其中，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，δε为应变释放，r为环芯圆柱表面圆心为原点的圆柱极坐标系的径向坐标，ε(r,h/D)为圆柱表面内任意点的局部应变。
在步骤B.5中，体平均残余应力释放系数γ的计算方法为：
Figure BDA0003690881190000074
其中，γ为体平均残余应力释放系数，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，σapp为数值计算时施加的已知等双轴应力。
在步骤B.6中，指定有效直径下的应变响应函数SIμ(h/D)的计算方法为：
Figure BDA0003690881190000081
其中，h为环槽总深度，D为环芯圆柱直径，SRμ为计算有效系数为μ时应变释放关于相对深度h/D的函数。
实施例2：
图1至4为本发明基于配有聚焦离子束的扫描电子显微镜、有限元计算和数字图像关联技术测量残余应力的实施例。本实施例的方法适用于固体部件的残余应力检测，特别是不能造成显著损伤，而实验室无损方法又难以测定的样品，如单晶、非晶、薄膜、强织构材料等。以下以测量单晶镍基高温合金发动机叶片的残余应力为例进行说明。
具体地，本实施例例基于配有聚焦离子束的扫描电子显微镜、计算和数字图像关联技术测量残余应力的方法分为实验、有限元计算和后处理三部分，包含以下步骤：
A.实验部分
A.1将尺寸为10mm×10mm×2mm的单晶镍基高温合金进行电解抛光使待测表面平整；
A.2样品用银胶固定于配有FIB的电子显微镜试样台上，并抽真空；
A.3同时打开离子束；在电子束下采用5kV的电压和0.8nA的束流对单晶合金对焦并调节至8000倍数下，将工作距离调整为4mm，插入Pt针并使用离子束在待测区(本实施例例中为10μm×10μm)沉积0.8μm厚，分布如图1a所示的随机分布的Pt散斑和保护环；
A.4沉积结束后，将样品台倾转52°，在30kV，41pA的离子束下找到已沉积斑点的位置后将离子束流调整至1.2nA；本实施例中环槽的尺寸为：内直径＝10μm，外直径＝15μm；
A.5以从外到内的切削方式，在放大8000倍(与步骤A.3中一致)下调整环切位置。在本实施例中对准指目标环切内圆与沉积的Pt保护环外圈重叠；
A.6如图1b至d所示，使用离子束切除深度为Δh(Δh＜＜D)的环槽，并在每次切削后立刻使用5kV，0.8nA的电子束成像，图像做52°倾斜校正(与步骤A.4中一致)并保存；
A.7交替重复步骤A.5和A.6；在本实施例中，切削和成像的循环持续至总槽深h＝D＝10μm，且Δh并非定值，具体过程为前4μm(图1a至c)的深度每次环切Δh＝0.5μm(图1b)，后6μm(图1c至d)每次环切Δh＝1.0μm；SEM图像分辨率为1536x1094。
B.仿真与计算部分
在了解材料参数的情况下，如本实施例中的单晶镍基高温合金，通过有限元仿真得到已知应力场时环切带来的应变释放。该结果在后处理中将作为计算参照。
B.1获得材料参数；
B.2如图2所示，通过计算机对完整和环切至不同相对槽深(h/D)的样品进行建模。在本实施例中取值为h/D＝0,0.025,0.05,0.075,0.1,0.15,0.2,0.3,0.5,0.6,0.8,1，采用了2D模型并关于圆柱中轴做轴对称以提高计算效率，有限元软件为ANSYS；
B.3施加任意载荷；本例中使用单位压载荷；
B.4通过有限元方法得到完整模型以及不同h/D时的体应变场，并取圆柱表面的应力场；通过对比不同h/D与完整模型的圆柱表面应力场，得到切割到不同相对深度时的应变释放分布；
B.5取不同有效直径，并计算对应圆内的有效应变释放。本例中有效直径选取圆柱直径D的1，0.8和0.6(有效直径Deff＝μD，直径计算有效系数μ＝1，0.8或0.6)，暨只采用圆柱表面一个直径为圆柱总直径1、0.8和0.6倍的同心圆内的数据进行后续运算；
B.6将应变释放关于h/D作图并进行拟合，得到应变释放关于相对深度的函数，计算方法为：
Figure BDA0003690881190000091
其中SRμ(h/D)为有效系数为μ时应变释放关于相对深度h/D的函数，r为环芯圆柱表面圆心为原点的圆柱极坐标系的径向坐标，如图2所示，ε(r,h/D)为圆柱表面的局部应变；将应变释放函数关于完整模型的应变进行归一，可以得到如图3所示的不受材料属性影响的通用归一化应变释放函数SRμ,norm.(h/D)；
B.7应变释放函数在相对深度大于一定值后将趋近一固定值，如图3所示(归一前后函数趋势一致)；该极限只与切割前圆柱体积内的平均残余应力相关，且与残余应力水平成正比；通过施加已知载荷时的最终释放和原点连线算出体平均残余应力释放系数，具体计算方法为：
Figure BDA0003690881190000092
其中σapp为数值计算时施加的已知应力，在本实施例中为单位压应力；
B.8对归一化应变释放函数关于相对深度进行求导，得到应变响应函数，具体计算方法为：
Figure BDA0003690881190000101
C.后处理部分
C.1使用数字图像关联技术技术对步骤A.6和A.7中保存的扫描电子显微图像进行计算，得到环芯圆柱表面随相对深度h/D增加而演化的应变场；在本例中的散斑和图像分辨率条件下，使用了70x70像素的subset和2像素的步长进行计算,数字图像关联技术软件为VIC-2D；
C.2取0.8和0.6(与步骤B.5中一致)的有效直径并得到相应有效范围内，沿计算方向的平均应变释放，将不同相对深度时的应变释放作图，可以得到实验应变释放曲线，如图4所示；
C.3将实验最终释放带入步骤B.7中通过有限元计算得到的理论体平均残余应力释放关系，计算圆柱体积内该方向的平均残余应力，具体计算方法为：
Figure BDA0003690881190000102
其中
Figure BDA0003690881190000103
为环芯内沿l方向的平均残余应力，SRμ,exp(h/D)为计算有效系数为μ时的实验应变释放函数；
C.4对实验应变释放进行拟合并求导，得到实验应变响应函数；与计算得到的应变响应函数进行对比，计算任意两张扫描电子显微图像对应深度区间内的平均残余应力；如第 n张扫描电子显微图像中槽深为0.5μm，相对深度h/D为0.05，第n+1张扫描电子显微图像中槽深为1μm，相对深度h/D为0.1。使用第n+1张的应变释放减去第n张的应变释放，其与理论响应函数SIμ从相对深度从0.05到0.1之间变化量的比值即为样品内深度从0.5 到0.1μm之间的平均残余应力。
C.5改变数字图像关联技术计算的应变方向并重复步骤C.1至C.4，得到不同方向的应力分量，可重构体积内或指定深度范围内的残余应力张量。
实施例3：本实施例与实施例3不同在于，采用电子束替换离子束。其它步骤和参数与实施例相同。
实施例4：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.09μm厚。其它步骤和参数与实施例相同。
实施例5：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.10μm厚。其它步骤和参数与实施例相同。
实施例6：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.11μm厚。其它步骤和参数与实施例相同。
实施例7：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.12μm厚。其它步骤和参数与实施例相同。
实施例8：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.13μm厚。其它步骤和参数与实施例相同。
实施例9：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.14μm厚。其它步骤和参数与实施例相同。
实施例10：本实施例与实施例3不同在于，插入Pt针并使用离子束在待测区沉积0.15μm厚。其它步骤和参数与实施例相同。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (10)

1.一种高效、高空间分辨率的残余应力测量方法，其特征在于，包括以下步骤：
A.实验部分
A.1对待测样品表面处理至没有划痕，然后置于配有聚焦离子束的扫描电子显微镜内的样品台上；
A.2用电子束或离子束在待测区沉积指定占空比的、随机分布的散斑；
A.3在待切削环槽内环用电子束或离子束沉积一个直径等于环芯直径D、宽度大幅小于直径的圆环；
A.4使用离子束进行深度为Δh，且Δh＜＜D的环切；
A.5使用电子束和离子束束成像，并保存采取了倾转角度补偿的扫描电子显微图片，补偿角度与样品台倾转角度一致；
A.6交替重复上述A.4与A.5步骤直至环槽的总深度h与环芯直径D的比值，暨相对深度h/D达到目标值，通常h/D大于1/2后圆柱内应变释放完全，环芯表面应变趋于稳定，切削结束；
B.有限元计算部分
B.1在有限元计算软件中建立完整和聚焦离子束逐层切除环槽后的材料模型，每道切除的深度无需和实验过程一致；
B.2对模型施加任意已知等双轴应力，经平衡计算得到不同的环槽相对深度h/D下的环芯表面应变分布；
B.3取不同有效直径对应的圆内应变分布，并分别计算随着环槽加深，应变释放的演化；
B.4根据应变释放得到应变释放函数SRμ(h/D)，并关于完整模型内的均匀应变进行归一得到不受材料属性影响、受计算有效系数影响的通用归一化应变释放函数SRμ,norm.(h/D)；
B.5通过应变释放函数得到体平均残余应力释放系数γ；
B.6对归一化应变释放函数SRμ,norm.(h/D)求导，得到指定有效直径下的应变响应函数SIμ(h/D)；
B.7通过实验测得的应变相应函数与上述响应函数在相同深度区间内的比值得到该深度区间内残余应力对环芯表面应变的影响；
C.后处理部分
C.1使用数字图像关联技术对扫描电子显微图像进行计算，得到环芯圆柱表面随环槽相对深度增加而演化的应变场；
C.2取指定有效直径并沿指定方向计算相应圆内平均应变释放，关于相对深度作图，得到实验应变释放曲线；
C.3将实验应变释放曲线的释放平台，暨相对深度趋于无限时的应变释放带入步骤B.5中的理论关系进行计算，计算圆柱体积内该方向的平均残余应力，过程中保持有效直径一致；
C.4根据实验应变释放曲线计算实验应变响应函数，与步骤B.6中的理论应变响应函数进行对比，以得到任意两张SEM图像对应深度区间内的平均残余应力，过程中保持有效直径一致；
C.5改变DIC计算的应变方向并重复步骤C.1至C.4，得到不同方向的应力分量。
2.根据权利要求1所述方法，其特征在于，在步骤A.1中，表面处理包括但不限于电解抛光或机械抛光常见表面处理方法。
3.根据权利要求1所述法，其特征在于，在步骤A.2中，散斑点的占空比为15％～85％。
4.根据权利要求1所述方法，其特征在于，在步骤A.2和A.3中，散斑点和保护环的厚度为0.05μm及以上。
5.根据权利要求1所述方法，其特征在于，在步骤A.2和A.3中，沉积散斑点和保护环覆盖的范围是0.5μm～100μm。
6.根据权利要求1所述方法，其特征在于，在步骤A.4中，聚焦离子束切削的内环直径范围为0.5μm～95μm，外环的直径范围为0.5μm～100μm。
7.根据权利要求1所述方法，其特征在于，在步骤A.4中，切削步进值Δh为25nm～2μm。
8.根据权利要求1所述方法，其特征在于，在步骤B.4中，应变释放函数SRμ(h/D)的计算方法为：
Figure FDA0003690881180000021
其中SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，δε为应变释放，r为环芯圆柱表面圆心为原点的圆柱极坐标系的径向坐标，ε(r，h/D)为圆柱表面内任意点的局部应变。
9.根据权利要求1所述方法，其特征在于，在步骤B.5中，体平均残余应力释放系数γ的计算方法为：
Figure FDA0003690881180000031
其中γ为体平均残余应力释放系数，SRμ(h/D)为计算有效系数为μ时应变释放关于相对深度h/D的函数，h为环槽总深度，D为环芯圆柱直径，σapp为数值计算时施加的已知等双轴应力。
10.根据权利要求1所述方法，其特征在于，在步骤B.6中，指定有效直径下的应变响应函数SIμ(h/D)的计算方法为：
Figure FDA0003690881180000032
CN202210661052.8A 2022-06-13 2022-06-13 一种高效、高空间分辨率的残余应力测量方法 Pending CN115165182A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210661052.8A CN115165182A (zh) 2022-06-13 2022-06-13 一种高效、高空间分辨率的残余应力测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210661052.8A CN115165182A (zh) 2022-06-13 2022-06-13 一种高效、高空间分辨率的残余应力测量方法

Publications (1)

Publication Number Publication Date
CN115165182A true CN115165182A (zh) 2022-10-11

Family

ID=83485517

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210661052.8A Pending CN115165182A (zh) 2022-06-13 2022-06-13 一种高效、高空间分辨率的残余应力测量方法

Country Status (1)

Country Link
CN (1) CN115165182A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117371171A (zh) * 2023-08-31 2024-01-09 湖北江城实验室科技服务有限公司 密封环可靠性的评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103411714A (zh) * 2013-08-08 2013-11-27 南京航空航天大学 一种基于局部剥层和系数修正的残余应力测量方法
CN106599496A (zh) * 2016-12-20 2017-04-26 上海电气电站设备有限公司 基于数值计算的环芯法测量残余应力的方法
CN111751035A (zh) * 2020-06-16 2020-10-09 武汉钢铁有限公司 残余应力分析方法及应用
CN114611354A (zh) * 2022-03-09 2022-06-10 西北工业大学 一种基于纳米压痕试验的电子封装材料残余应力计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103411714A (zh) * 2013-08-08 2013-11-27 南京航空航天大学 一种基于局部剥层和系数修正的残余应力测量方法
CN106599496A (zh) * 2016-12-20 2017-04-26 上海电气电站设备有限公司 基于数值计算的环芯法测量残余应力的方法
CN111751035A (zh) * 2020-06-16 2020-10-09 武汉钢铁有限公司 残余应力分析方法及应用
CN114611354A (zh) * 2022-03-09 2022-06-10 西北工业大学 一种基于纳米压痕试验的电子封装材料残余应力计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
J.G. ZHU等: "Interfacial Residual Stress Analysis of Thermal Spray Coatings by Miniature Ring-Core Cutting Combined with DIC Method", EXPERIMENTAL MECHANICS, vol. 54, pages 127 - 136 *
KORSUNSKY等: "Residual stress evaluation at the micrometer scale: Analysis of thin coatings by FIB milling and digital image correlation", SURFACE AND COATINGS TECHNOLOGY, vol. 205, no. 7, pages 2393 - 2403, XP027533953 *
SEBASTIANI等: "Depth-resolved residual stress analysis of thin coatings by a new FIB-DIC method", MATERIALS SCIENCE AND ENGINEERING:A, vol. 528, no. 27, pages 7901 - 7908 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117371171A (zh) * 2023-08-31 2024-01-09 湖北江城实验室科技服务有限公司 密封环可靠性的评估方法
CN117371171B (zh) * 2023-08-31 2024-05-24 湖北江城实验室科技服务有限公司 密封环可靠性的评估方法

Similar Documents

Publication Publication Date Title
US6882739B2 (en) Method and apparatus for rapid grain size analysis of polycrystalline materials
Vlieg et al. Relaxations in the missing-row structure of the (1× 2) reconstructed surfaces of Au (110) and Pt (110)
Serrano-Munoz et al. The residual stress in as-built Laser Powder Bed Fusion IN718 alloy as a consequence of the scanning strategy induced microstructure
Dingley et al. Microtexture determination by electron back-scatter diffraction
EP2818852B1 (en) Methods and devices for high throughput crystal structure analysis by electron diffraction
Brückner et al. The influence of the dendritic structure on the γ/γ′-lattice misfit in the single-crystal nickel-base superalloy CMSX-4
Guéninchault et al. Nanox: a miniature mechanical stress rig designed for near-field X-ray diffraction imaging techniques
Slavik et al. Determining fracture facet crystallography using electron backscatter patterns and quantitative tilt fractography
CN115165182A (zh) 一种高效、高空间分辨率的残余应力测量方法
JP2013083574A (ja) 塑性ひずみの評価システムおよび評価方法
Davies et al. Combined application of electron backscatter diffraction and stereo‐photogrammetry in fractography studies
Davidson Uses of electron channelling in studying material deformation
JP2016191562A (ja) 結晶方位の解析方法および解析装置
Klocke et al. EBSD-analysis of flexure hinges surface integrity evolution via wire-EDM main and trim cut technologies
US20220412901A1 (en) Device and method for measuring short-wavelength characteristic x-ray diffraction based on array detection
WO2023101800A1 (en) Method of locating poles on bearing elements
Altinkurt et al. Local strain redistribution in a coarse-grained nickel-based superalloy subjected to shot-peening, fatigue or thermal exposure investigated using synchrotron X-ray Laue microdiffraction
Bettles et al. In situ observation of strain evolution in CP-Ti over multiple length scales
Randle Crystallographic analysis of facets using electron backscatter diffraction
Gussev et al. Analysis of Localized Deformation Processes in Highly Irradiated Austenitic Stainless Steel through In Situ Techniques
Brewer et al. Comparison of diffraction methods for measurement of surface damage in superalloys
Chighizola et al. Intermethod comparison and evaluation of near surface residual stress in aluminum parts subject to various milling parameters
Statnik et al. Polar transformation of 2D X-ray diffraction patterns for 2D strain evaluation
Pai et al. Diffraction-Based Multiscale Residual Strain Measurements
Zurcher et al. Sliding wear resistance and residual stresses of parts repaired by laser metal deposition

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