一种基于原位SEM的裂纹尖端张开角度测量方法与流程

本发明涉及断裂力学研究,尤其涉及一种基于原位sem的裂纹尖端张开角度（ctoa）测量方法。

背景技术

弹塑性断裂力学及其测试方法是评估延性材料断裂特性的重要研究手段。在工程实际中，对于飞机蒙皮、石油管道等薄壁、薄板结构，其断裂失效工况接近于低约束条件。对于这类问题，测定临界裂纹尖端张开角（cracktipopeningangle，ctoa）被认为是一种非常适合用于评估该类材料与构件失效条件的断裂评估准则，并得到广泛的应用。裂纹尖端张开角在裂纹稳态扩展阶段达到临界值且保持恒定，该临界值可以用于表征材料的稳定裂纹扩展韧性，即断裂韧度。断裂韧度反映材料的抗裂纹扩展能力，断裂韧度越大，韧度越高。基于ctoa的理论判据为：ctoa≤ctoac，当ctoa达到临界起始值ctoac时，裂纹开始稳定扩展。

目前，ctoa的试验测试方法主要有直读法、3d显微形貌法、δ5间接测量法。直读法记录试验过程中的裂纹稳态扩展过程的图像，试验完毕后需测定采集的图像来计算ctoa，通过裂尖和裂尖后的定位点计算裂纹尖端张开角的平均值。3d显微形貌法通过试验后的法向分离测量断裂表面的高度，进而获得两个离散定义的数学表面，两个离散表面经最小二乘拟合出平均斜率进而求得裂纹尖端张开角。δ5间接测量法技术通过测量f-δa曲线与f-δ5曲线，进而得出δ5-δa曲线，ctoa近似等于δ5-δa曲线的斜率。

现有试验测试方法在实际应用中普遍存在一定不足，因此受到不同程度的制约。直读法中，光学测量需准确识别裂纹尖端点，识别精度不够容易出现数据异常离散化；3d显微形貌法是在试验后对试样的裂纹面进行测量，不能实时地观察记录加载过程中ctoa的变化规律；δ5间接测量法需采用标准试样，试验过程中需采用柔度法并借助cod规来测量位移增量，程序比较复杂。

原位sem试验技术是近年来发展起来的一种新型材料力学性能测试与研究手段，该技术为从微观尺度分析与表征材料疲劳与断裂失效行为提供了强有力的研究手段，为从机理层面揭示材料的物理失效过程提供了必要的技术支撑。所谓的原位sem试验技术，即扫描电子显微镜（sem）可实现在原位加载条件下实时观察试样在变形过程中的表面形貌与损伤演化，而且实时记录的sem图像可与试样的宏观应力-应变曲线在时间上建立一一对应关系，便于进行断裂力学相关分析。本发明借助于原位sem测试技术，开发一种基于matlab的裂纹尖端张开角的自动测量和读取方法。该方法具有非接触、多尺度、高精度、操作简单等优点，可实时获取裂纹尖端变形形貌和裂纹尖端张开角，自动化程度高，测量结果准确、可靠，在材料性能研究中具有应用潜力。

技术实现要素：

本发明主要内容为借助于原位sem测试技术获取试样在原位加载下的sem图像，然后利用基于matlab语言的“ctoa计算方法”模型计算获取图像上对应的裂纹尖端张开角度。原位加载过程中，使用扫描电镜可以测量出图像中的裂纹长度a，扫描电镜可以根据一定的加载时间或裂纹增长量δa依次采集多组图像。最后，根据“ctoa计算方法”模型计算得到的ctoa值和由扫描电镜测量出的裂纹长度a,即可获得ctoa-δa的演化趋势。

该方法操作简单，直观，可广泛应用于材料的ctoa测量，可以显著提高效率。同时，使研究者在原位观察材料失效过程的同时，获得可定量表征其断裂韧性的参数，拓宽了原位测试技术的应用范围。

基于原位sem试验技术的裂纹尖端张开角度的测量方法具体包括如下步骤：

步骤(1)试样准备：对试样进行机加工，达到所需的几何尺寸及表面粗糙度。

步骤(2)装夹试样，扫描电镜成像：将试样装夹于原位试验机的试样台上，调整扫描电镜的倍率和对比度，直至可清晰观察到试样表面；

步骤(3)设置试验条件开始试验：设定原位试验机的加载参数，采用拉-拉疲劳载荷预制裂纹。

步骤(4)断裂试验和采集图像：预制裂纹完毕后，采用单向拉伸载荷进行断裂试验，与此同时，扫描电镜每隔一定的加载时间对裂纹尖端进行图像采集。

步骤(5)计算ctoa值和测量裂纹长度a：试验完毕后，用扫描电镜依次测量图像中的裂纹长度a，并将图像依次输入到“ctoa计算方法”模型中计算对应的ctoa值。

步骤(6)获得ctoa-δa的演化趋势：将计算得到的ctoa和扫描电镜测得的裂纹扩展增量δa绘制在同一个坐标系中，得到ctoa-δa的趋势图。

在步骤(1)中，试样采用单边缺口拉伸试样（singleedgenotchtensile）,即sent试样，缺口采用中走丝线切割加工。依次使用1000#和2000#的砂纸对试样打磨，随后用直径为3μm和0.5μm的金刚石抛光液抛光，清洗，以光学显微镜（约400x）下观察不到试样表面划痕为标准。

在步骤(2)中，调整扫描电镜的视野清晰度，以1000x下能清晰观察试样表面为标准。观测区域需包含裂尖，且包含裂尖后距离裂尖不小于1.5mm的区域。

在步骤(5)中，“ctoa计算方法”模型是基于matlab语言编写的程序代码，通过对图像中的灰度值进行分析比较，求出两条裂纹线上若干个离散点的位置，对两条裂纹线上的点分别进行最小二乘拟合，可认为两条拟合直线之间夹角的大小即为ctoa值。

本发明与现有技术相比，具有以下优点及突出性效果：

(1)提供了一种原位观察裂纹扩展行为的试验手段。本发明中涉及的原位sem试验技术，可以在介观尺度下实时观测材料的变形行为和裂纹扩展行为。有助于研究材料的微观变形行为和断裂机理。

(2)本发明中涉及的sem成像系统，相对于光学显微镜具有更高的分辨率和清晰度。

(3)具有非接触测量的优点，无需对试样表面进行特殊处理，过程中无复杂的试验操作步骤，后期处理方便，精度较高。试验过程中，可同步进行裂纹扩展速率研究和裂纹扩展路径分析等其他研究工作。

附图说明

图1为“ctoa计算方法”的原理流程图。

图2为“ctoa计算方法”的原理示意图。

图3为实例中sent试样尺寸示意图。

图4为实例中扫描电镜采集的其中一张图像。

图5为实例中“ctoa计算方法”的输出图。

图6为实例中“ctoa计算方法”的输出图的局部放大图。

图7为实例中通过数据处理得到的ctoa-δa的演化趋势图。

具体实施方式

本发明基于原位sem测试技术获取试样在原位加载下的sem图像，然后利用基于matlab语言的“ctoa计算方法”模型计算获取图像上对应的裂尖端张开角度。使用扫描电镜测量出图像中的裂纹扩展量δa，再结合由“ctoa计算方法”模型计算得到的ctoa值，即可获取ctoa-δa的演化趋势。

图1为本发明的“ctoa计算方法”原理流程图，流程图中对“ctoa计算方法”按程序处理的先后顺序进行简要概述。

图2为本发明的“ctoa计算方法”原理示意图。“ctoa计算方法”的原理：导入裂纹尖端图像；调用matlab中的im2bw函数，以设定的阈值将图像二值化，使得图像中的像素点只有两种灰度等级，灰度值大于阈值的像素点为1，灰度值小于阈值的像素点为0；选择裂尖点a和裂尖后1.5mm处点b，在y轴方向上，a→b区域即参与计算的区域；l是以a为起点的y轴方向线段，终点位于点b对应的x轴向线上，在线段l上均布n个点li，i=1~n；以li为起点，沿x轴正方向对像素点进行逐个扫描，当扫描到某个像素点的灰度等级和li的灰度等级不同时，该像素点即存在于裂纹线上的像素点，记录此像素点；随后从li+1点处开始扫描，以此类推，直到i=n；同理，以li为起点，沿x轴负方向扫描，直到i=n；至此，记录了以裂尖为起点，每条裂纹线上的n个像素点。分别将每条裂纹线上的像素点按最小二乘法进行直线拟合，计算两条直线间的夹角，即对应的ctoa值。

为了更好地理解本发明，下面结合实例对本发明作进一步的详细说明。本发明以单边缺口拉伸试样为例，对其加载疲劳载荷以预制裂纹，随后进行拉伸断裂试验。断裂试验过程中，扫描电镜观测试样表面裂纹尖端区域并采集裂纹尖端区域的图像。然后利用基于matlab语言的“ctoa计算方法”模型计算获取图像上对应的裂尖端张开角度，再结合扫描电镜测量出的裂纹长度a，获得ctoa-δa的演化趋势图。具体实施方式如下：

试样准备。本例选用材料为工程中常用的316l不锈钢，该材料具有很好的强度和延展性。将316l不锈钢进行机加工，试样为长度50mm，宽度10mm，厚度1mm，缺口长度1mm的单边缺口拉伸试样，试样尺寸示意图见图3。将加工好的试样用1000#和2000#的砂纸打磨，随后将试样置于抛光机上进行抛光，依次使用直径3μm和0.5μm的金刚石抛光液，磨抛过后用超声波清洗机清洗试样，然后烘干。

装夹试样，调节扫描电镜使清晰成像。将试样装夹在原位试验机的夹头上，将试样连同夹头整体置于扫描电镜真空箱中，使得试验过程处于扫描电镜的观察视场中。调整扫描电镜的焦距、对比度、像散，直至能在1000x下清晰观察到试样表面。

设置试验条件开始试验。使用拉-拉疲劳载荷对试样进行预制裂纹，应力比r=0.1，σmax=0.6σy(σmax为最大名义应力，σy为屈服强度)，最终预制裂纹长度为2mm，此时，裂纹总长度为3mm。

断裂试验和图像采集。预制裂纹完毕后，对试样施加轴向拉伸载荷进行断裂试验，试验采用位移控制，应变率为10^-4/s。每40s暂停试验，扫描电镜采集裂纹尖端区域的清晰图像，倍率为14x-100x，直到试样发生断裂破坏。

计算ctoa值和测量裂纹扩展量δa。试验完毕后，用扫描电镜依次测量出图像中的裂纹扩展量δa，将图像依次输入到“ctoa计算方法”模型中计算对应的ctoa值。本例中选取其中一张图像，见图4，将此图的处理过程予以展示。图5为“ctoa计算方法”的输出图，该图像对应的ctoa值为23.91°。图6为输出图5的局部放大图，可以看到通过程序识别的裂纹面上的若干个离散点，以及由这些离散点拟合成的直线。由此，获得每幅图像对应的裂纹扩展量δa和裂纹尖端张开角ctoa。进而得到裂纹稳定扩展过程中裂纹扩展量δa与相应裂纹尖端张开角ctoa的关系，而ctoa趋于稳定状态时的下限值即被认为是裂纹尖端张开角的临界值ctoac。

获取ctoa-δa的演化趋势。将计算得到的ctoa和扫描电镜测得的裂纹扩展增量δa绘制在同一个坐标系中，可以得到ctoa-δa的趋势图，如图7所示。由图7可以得到临界裂纹尖端张开角ctoac，该临界值ctoac可作为该材料（316l不锈钢）的一种断裂韧性，用于其在低约束工况下的断裂失效评估。

本发明提供了一种原位观察裂纹扩展行为的试验手段，可以在介观尺度下实时观测材料的变形行为和裂纹扩展行为。有助于研究材料的微观变形行为和断裂机理。

本发明中通过sem成像系统采集的裂纹尖端图像，具有高的分辨率和清晰度。

本发明实现了ctoa的非接触测量，无需对试样表面进行特殊处理，过程中无复杂的试验操作步骤，后期处理方便，精度较高。试验过程中，可同步进行裂纹扩展速率研究和裂纹扩展路径分析等其他研究工作。