一种应变场声发射事件时空强的测量方法
【专利摘要】一种应变场声发射事件时空强的测量方法，步骤为：利用数字图像相关方法，获取材料或结构变形过程中一个表面的最大剪切应变场；将选定的图像分成若干个正方形像素块，根据最大剪切应变场，利用插值方法获得各像素块的最大剪切应变和各像素块中心坐标；确定各破坏像素块及各破坏像素块在测量间隔内释放的弹性应变能、释放弹性应变能的次数；查找相互连通的破坏像素块构成声发射事件；测量各声发射事件在测量间隔内释放的弹性应变能、空间位置、持续时间；确定各声发射事件在各测量间隔内释放的弹性应变能、空间尺度、持续时间的统计规律。优点是：本发明实现了对声发射事件的时空强的全方位测量，扩展了声发射技术的既有功能，应用前景广阔。
【专利说明】
一种应变场声发射事件时空强的测量方法
技术领域
[0001] 本发明属于应变场测量技术领域，特别涉及一种应变场声发射事件时空强的测量 方法。
【背景技术】
[0002] 声发射技术是靠材料在受载条件下发射的弹性波来获取其内部状态和力学特性 的一种实验方法。声发射技术可应用于实验室岩石类材料的微破裂测量、天然地震和矿井 岩爆的监测和预报等多个领域。
[0003] 基于声发射技术的实验结果强烈依赖于实验设备。虽然自70年代以来，国内外已 经普遍采用了数字化全波形声发射测量系统，但由于存在动态范围不足与"死时间"等问 题，使得声发射波形失真，声发射事件大量缺失，从而影响声发射的时空分布、一些统计量 的可靠性以及对材料变形破坏过程的认识。声发射技术对于低频事件一般并不敏感。由于 声发射系统多方面的差异，不同实验室的结果难于直接对比。随着时代的发展和技术的进 步，上述问题已在一定程度上得到了克服，但仍无法根除。和过去相比，人们对材料变形破 坏过程的研究在深度和广度上都有所超越，技术上的局限性严重阻碍了一些研究向前推 进。
[0004] 数字图像相关方法是材料或结构变形破坏过程中位移场和应变场观测的一种重 要手段，其原理是通过比较变形前后两块像素子区的相关程度来实现子区中心点位移和子 区应变的测量，具有实时观测、光路简单、精度高的特点。
[0005] 声发射技术和数字图像相关方法有各自的适用领域，各具优势。前者适用于探测 声发射事件的时空分布规律，后者适于探测位移场和应变场。目前，根据测量获得的应变场 信息挖掘声发射事件的时空强分布规律还未见报道，此规律的获得有助于深刻认识材料的 变形破坏过程。
【发明内容】
[0006] 针对现有技术的不足，本发明提出一种应变场声发射事件时空强的测量方法。该 方法的具体步骤如下：
[0007] 步骤1、利用数字图像相关方法，获取材料或结构变形破坏过程中一个表面的最大 剪切应变场；
[0008] 步骤1.1、利用拍摄设备采集材料或结构变形破坏过程中一个表面的图像，如果材 料或结构表面没有天然纹理或天然纹理质量较差，需要利用涂料在材料或结构表面制作散 斑场，任意相邻两张图像的拍摄时间间隔相同；
[0009] 步骤1.2、选定若干张图像，设置子区尺寸、测点数目、测点位置和测点间距，利用 数字图像相关方法测量材料或结构变形破坏过程中一个表面的应变场，所述测点应布置在 选定的第一张图像上；
[0010] 步骤1.3、利用应变场获得最大剪切应变场；
[0011] 步骤2、将选定的第一张图像分成若干个无间隔、无重叠、成行成列的正方形像素 块，根据最大剪切应变场，利用插值方法获得各像素块的最大剪切应变和各像素块中心坐 标；
[0012] 步骤3、根据图像中各像素块的最大剪切应变和设置的强度参数，确定各破坏像素 块，并分别确定各破坏像素块在测量间隔内释放的弹性应变能和各测量间隔释放弹性应变 能的次数，所述测量间隔包括若干个能量释放计算间隔，所述能量释放计算间隔为任意相 邻两张图像的拍摄时间间隔，所述破坏像素块为发生塑性变形的像素块；
[0013] 步骤3.1、将声发射事件的能量释放计算间隔设置为任意相邻两张图像的拍摄时 间间隔，同时设置测量间隔和测量间隔布置方式，所述测量间隔包括若干个能量释放计算 间隔，所述测量间隔布置方式包括重叠布置、无缝不重叠布置和有缝布置，重叠布置是指两 次相邻的测量间隔包含部分相同的能量释放计算间隔，无缝不重叠布置是指两次相邻的测 量间隔紧邻且包含的能量释放计算间隔均不相同，有缝布置是指两次相邻的测量间隔有间 隙，间隙部分的能量释放计算间隔不属于任何一次测量间隔，且两次相邻的测量间隔包含 的能量释放计算间隔均不相同；
[0014] 步骤3.2、根据图像中各像素块的最大剪切应变和设置的强度参数，确定破坏像素 块，分别测量各破坏像素块在各能量释放计算间隔释放的弹性应变能；
[0015] 步骤3.3、分别对各破坏像素块在各能量释放计算间隔释放的弹性应变能进行求 和，获得各测量间隔内各破坏像素块释放的弹性应变能；
[0016] 步骤3.4、分别测量各破坏像素块在各测量间隔内释放弹性应变能的次数；
[0017] 步骤4、将任一个破坏像素块标记为i，i为声发射事件的标号，取为自然数，遍历与 该像素块具有相同的点或相同的线的所有像素块，若这些像素块中存在破坏像素块，则这 些破坏像素块同样标记为i，同样标记为i的破坏像素块组成一组相互连通的破坏像素块， 构成一个声发射事件，同样标记为i的破坏像素块的数目为声发射事件的空间尺度，以此类 推，找到应变场中所有的声发射事件，进而确定应变场中声发射事件的数目和各声发射事 件的空间尺度，所述一组相互连通的破坏像素块由若干个具有相同的点或者相同的线的破 坏像素块组成，所述各声发射事件的空间尺度为各声发射事件包含的破坏像素块的数目；
[0018] 步骤5、将每组相互连通的破坏像素块在测量间隔内释放的弹性应变能求和，作为 每个声发射事件在测量间隔内释放的弹性应变能；
[0019] 步骤6、根据形心坐标公式，确定各声发射事件所在的空间位置；
[0020] 步骤7、利用各破坏像素块在各测量间隔内释放弹性应变能的次数乘以能量释放 计算间隔，测量各破坏像素块释放能量的持续时间，通过对每组相互连通的破坏像素块释 放能量的持续时间求均值或者求并集，确定各声发射事件在测量间隔内的持续时间；
[0021] 步骤8、在各测量间隔内，确定各声发射事件释放的弹性应变能、各声发射事件的 空间尺度、各声发射事件的持续时间的统计规律；
[0022] 步骤8.1、在各测量间隔内，将各声发射事件释放的弹性应变能按最大值和最小值 为首末端均分为能量释放若干等级、将各声发射事件的空间尺度按最大值和最小值为首末 端均分为空间尺度若干等级、将各声发射事件的持续时间按最大值和最小值为首末端均分 为持续时间若干等级；
[0023] 步骤8.2、将各声发射事件释放的弹性应变能、各声发射事件的空间尺度、各声发 射事件的持续时间分别划归入相应的等级中，并统计能量释放各等级的频次、空间尺度各 等级的频次、持续时间各等级的频次，所述能量释放各等级的频次指能量释放各等级包含 的相应的声发射事件的个数，所述空间尺度各等级的频次指空间尺度各等级包含的相应的 声发射事件的个数，所述持续时间各等级的频次指持续时间各等级包含的相应的声发射事 件的个数；
[0024]步骤8.3、利用能量释放各等级的数据，通过线性回归确定频次-能量释放关系，利 用空间尺度各等级的数据，通过线性回归确定频次-空间尺度关系，利用持续时间各等级的 数据，通过线性回归确定频次-持续时间关系。
[0025] 有益效果：
[0026] 本发明提出了一种应变场声发射事件时空强的测量方法，该方法采用拍摄设备采 集材料或结构变形破坏过程中一个表面的图像，利用数字图像相关方法，获取材料或结构 变形破坏过程中一个表面的应变场，根据图像中各像素块的最大剪切应变和设置的强度参 数确定破坏像素块，通过在各测量间隔内确定图像中相互连通的破坏像素块并计算各相互 连通的破坏像素块释放的弹性应变能，进而确定各声发射事件的空间尺度和各声发射事件 释放的弹性应变能;通过统计各测量间隔内各破坏像素块释放弹性应变能次数，计算各破 坏像素块释放能量的持续时间，进而确定各声发射事件的持续时间；最后利用统计方法确 定频次_能量释放关系、频次_空间尺度关系、频次_持续时间关系，获得弹性应变能、空间尺 度、持续时间三种统计量的演变规律。本发明实现了对声发射事件的时(持续时间）空(空间 尺度)强(能量释放)的全方位测量，拓展了声发射技术的既有功能，应用前景广阔。
【附图说明】
[0027] 图1为本发明的应变场声发射事件时空强的测量方法的流程图；
[0028] 图2为本发明实施例的应变场中声发射事件的分布图；
[0029] 图3为本发明实施例的8个破坏像素块在测量间隔内释放弹性应变能的示意图；
[0030] 图4为本发明实施例的声发射事件空间尺度的分布图；
[0031]图5为本发明实施例的声发射事件能量释放的分布图；
[0032]图6为本发明实施例的声发射事件持续时间的分布图；
[0033]图7为本发明实施例的声发射事件空间尺度的排序图；
[0034]图8为本发明实施例的声发射事件能量释放的排序图；
[0035]图9为本发明实施例的声发射事件持续时间的排序图；
[0036] 图10为本发明实施例的声发射事件频次-空间尺度的统计图；
[0037] 图11为本发明实施例的声发射事件频次-能量释放的统计图；
[0038]图12为本发明实施例的声发射事件频次-持续时间的统计图。
【具体实施方式】
[0039] 下面结合附图对本发明【具体实施方式】做详细说明。一种应变场声发射事件时空强 的测量方法，如图1所示：
[0040] 步骤1、利用数字图像相关方法，获取材料或结构变形破坏过程中一个表面的最大 剪切应变场；
[0041] 步骤1.1、利用拍摄设备采集材料变形破坏过程中一个表面的图像，任意相邻两张 图像的拍摄时间间隔均为〇. 2s;
[0042] 步骤1.2、选定11张图像，设置子区尺寸为21X21像素，测点数目为10X15,测点成 行成列等间距布置，测点间距为50像素，利用数字图像相关方法测量材料或结构变形破坏 过程中一个表面的应变场，所述测点应布置在选定的第一张图像上；
[0043] 步骤1.3、利用应变场获得最大剪切应变场
，其中ymax为最 大剪切应变，ex、ey、y xy分别为x方向正应变，y方向正应变以及剪切应变。
[0044]步骤2、将选定的第一张图像分成无间隔、无重叠、成行成列的20X30个正方形像 素块，根据获得的最大剪切应变场，利用插值方法获得各像素块的最大剪切应变和各像素 块中心坐标。
[0045] 步骤3、根据图像中各像素块的最大剪切应变和设置的强度参数，确定各破坏像素 块，并分别确定各破坏像素块在测量间隔内释放的弹性应变能和各测量间隔释放弹性应变 能的次数；
[0046] 步骤3.1、设置声发射事件的能量释放计算间隔为任意相邻两张图像的拍摄时间 间隔，设置测量间隔为10个能量释放计算间隔，设置测量间隔布置方式为无缝不重叠布置；
[0047] 步骤3.2、根据图像中各像素块的最大剪切应变和设置的强度参数，确定破坏像素 块，所述破坏像素块是发生塑性变形的像素块，若像素块的最大剪切应变大于设置的强度 参数，则认为该像素块发生了塑性变形，即y max多y f，其中y f是设置的强度参数，取为〇. 2， 分别计算各破坏像素块在各能量释放计算间隔释放的弹性应变能；
[0048] -个破坏像素块存储的弹性应变能的计算公式为：
L0050J 其中，E为弹性模量，取为20GPa，ii为泊松比，取为0.15，V为破坏像素块的面积乘以 破坏像素块的厚度，破坏像素块的面积为25X25像素，破坏像素块的厚度取为25像素，〇1和 〇3分别为第1和第3主应力，当 〇1和〇3取为一个能量释放计算间隔开始时的值时，1]=山，当 和〇3取为一个能量释放计算间隔结束时的值时，U = U2, AUiUi-m即为一个能量释放计算 间隔该破坏像素块释放的弹性应变能；
[0051 ] 主应力〇1和〇3可由主应变求得：
[0053] 其中，￡1和￡3分别为第1和第3主应变；
[0054] 主应变￡1和￡3可由应变分量求得：
[0056] 步骤3.3、分别对各破坏像素块在各能量释放计算间隔释放的弹性应变能进行求 和，获得各测量间隔内各破坏像素块释放的弹性应变能；
[0057] 步骤3.4、分别统计各破坏像素块在各测量间隔内释放弹性应变能的次数。
[0058]步骤4、将任一个破坏像素块标记为i，i为声发射事件的标号，取为自然数，遍历与 该像素块具有相同的点或相同的线的所有像素块，若这些像素块中存在破坏像素块，则这 些破坏像素块同样标记为i，同样标记为i的破坏像素块组成一组相互连通的破坏像素块， 构成一个声发射事件，同样标记为i的破坏像素块的数目为声发射事件的空间尺度，以此类 推，找到应变场中所有的声发射事件，如图2所示，灰色像素块为破坏像素块，共84个，其中8 个破坏像素块在测量间隔内的各能量释放计算间隔释放的弹性应变能如图3所示，84个破 坏像素块共构成50个声发射事件，确定应变场中声发射事件的数目和各声发射事件的空间 尺度，如图4和图7所示，图4展示了各声发射事件空间尺度的分布图，声发射事件空间尺度 越小，其对应的圆形面积越小，图7中各声发射事件已按照声发射事件空间尺度的大小重新 标号，空间尺度越小，标号越小，空间尺度为1的声发射事件最多，达到30个，空间尺度为3的 声发射事件最少，为2个。
[0059] 步骤5、将每组相互连通的破坏像素块在测量间隔内释放的弹性应变能求和，作为 每个声发射事件在测量间隔内释放的弹性应变能，如图5和图8所示，图5展示了各声发射事 件能量释放的分布图，声发射事件能量释放越小，其对应的圆形面积越小，图8中各声发射 事件已按照声发射事件释放弹性应变能的大小重新编号，释放弹性应变能越小，编号越小。
[0060] 步骤6、根据形心坐标公式，确定各声发射事件所在的空间位置，
其中x。、y。、z。为一个声发射事件所在的空间位置坐标，Xi、yi、Zi ，- 为一个声发射事件中m个破坏像素块的中心坐标，Si为一个声发射事件中各破坏像素块的 m 面积，m为一个声发射事件包含的破坏像素块的数目，$ = 为一个声发射事件各破坏像 f=l 素块的面积之和。
[0061]步骤7、利用各破坏像素块在各测量间隔内释放弹性应变能的次数乘以能量释放 计算间隔，测量各破坏像素块释放能量的持续时间，通过对每组相互连通的破坏像素块释 放能量的持续时间求均值或者求并集，确定各声发射事件在测量间隔内的持续时间，如图6 和图9所示，图6展示了各声发射事件持续时间的分布图，声发射事件持续时间越短，其对应 的圆形面积越小，图9中各声发射事件以按照声发射事件持续时间的长短重新编号，持续时 间越短，编号越小。
[0062] 步骤8、在各测量间隔内，确定各声发射事件释放的弹性应变能、各声发射事件的 空间尺度、各声发射事件的持续时间的统计规律；
[0063] 步骤8.1、在各测量间隔内，将各声发射事件的空间尺度按最大值和最小值为首末 端均分为4个空间尺度等级、将各声发射事件释放的弹性应变能按最大值和最小值为首末 端均分为5个能量释放等级、将各声发射事件的持续时间按最大值和最小值为首末端均分 为5个持续时间等级；
[0064] 步骤8.2、将各声发射事件释放的弹性应变能、各声发射事件的空间尺度、各声发 射事件的持续时间分别划归入相应的等级中，统计空间尺度各等级的频次，如图10所示，统 计能量释放各等级的频次，如图11所示，统计持续时间各等级的频次，如图12所示；
[0065]步骤8.3、利用能量释放各等级的数据，通过线性回归确定频次-能量释放关系，利 用空间尺度各等级的数据，通过线性回归确定频次-空间尺度关系，利用持续时间各等级的 数据，通过线性回归确定频次-持续时间关系。
【主权项】
1. 一种应变场声发射事件时空强的测量方法，其特征在于:包括如下步骤： 步骤1、利用数字图像相关方法，获取材料或结构变形破坏过程中一个表面的最大剪切 应变场； 步骤1.1、利用拍摄设备采集材料或结构变形破坏过程中一个表面的图像，任意相邻两 张图像的拍摄时间间隔相同； 步骤1.2、选定若干张图像，设置子区尺寸、测点数目、测点位置和测点间距，利用数字 图像相关方法测量材料或结构变形破坏过程中一个表面的应变场； 步骤1.3、利用应变场获得最大剪切应变场； 步骤2、将选定的图像分成若干个无间隔、无重叠、成行成列的正方形像素块，根据最大 剪切应变场，利用插值方法获得各像素块的最大剪切应变和各像素块中心坐标； 步骤3、根据图像中各像素块的最大剪切应变和设置的强度参数，确定各破坏像素块， 并分别确定各破坏像素块在测量间隔内释放的弹性应变能和各测量间隔释放弹性应变能 的次数，所述测量间隔包括若干个能量释放计算间隔，所述能量释放计算间隔为任意相邻 两张图像的拍摄时间间隔，所述破坏像素块为发生塑性变形的像素块； 步骤4、将任一个破坏像素块标记为i，i为声发射事件的标号，取为自然数，遍历与该像 素块具有相同的点或相同的线的所有像素块，若这些像素块中存在破坏像素块，则这些破 坏像素块同样标记为i，同样标记为i的破坏像素块组成一组相互连通的破坏像素块，构成 一个声发射事件，同样标记为i的破坏像素块的数目为声发射事件的空间尺度，以此类推， 找到应变场中所有的声发射事件，进而确定应变场中声发射事件的数目和各声发射事件的 空间尺度，所述一组相互连通的破坏像素块由若干个具有相同的点或者相同的线的破坏像 素块组成，所述各声发射事件的空间尺度为各声发射事件包含的破坏像素块的数目； 步骤5、将每组相互连通的破坏像素块在测量间隔内释放的弹性应变能求和，作为每个 声发射事件在测量间隔内释放的弹性应变能； 步骤6、根据形心坐标公式，确定各声发射事件所在的空间位置； 步骤7、利用各破坏像素块在各测量间隔内释放弹性应变能的次数乘以能量释放计算 间隔，测量各破坏像素块释放能量的持续时间，通过对每组相互连通的破坏像素块释放能 量的持续时间求均值或者求并集，确定各声发射事件在测量间隔内的持续时间； 步骤8、在各测量间隔内，确定各声发射事件释放的弹性应变能、各声发射事件的空间 尺度、各声发射事件的持续时间的统计规律。2. 根据权利要求1所述的一种应变场声发射事件时空强的测量方法，其特征在于:所述 步骤3具体步骤如下： 步骤3.1、将声发射事件的能量释放计算间隔设置为任意相邻两张图像的拍摄时间间 隔，同时设置测量间隔和测量间隔布置方式，所述测量间隔包括若干个能量释放计算间隔， 所述测量间隔布置方式包括重叠布置、无缝不重叠布置和有缝布置，重叠布置是指两次相 邻的测量间隔包含部分相同的能量释放计算间隔，无缝不重叠布置是指两次相邻的测量间 隔紧邻且包含的能量释放计算间隔均不相同，有缝布置是指两次相邻的测量间隔有间隙， 间隙部分的能量释放计算间隔不属于任何一次测量间隔，且两次相邻的测量间隔包含的能 量释放计算间隔均不相同； 步骤3.2、根据图像中各像素块的最大剪切应变和设置的强度参数，确定破坏像素块， 分别测量各破坏像素块在各能量释放计算间隔释放的弹性应变能； 步骤3.3、分别对各破坏像素块在各能量释放计算间隔释放的弹性应变能进行求和，获 得各测量间隔内各破坏像素块释放的弹性应变能； 步骤3.4、分别测量各破坏像素块在各测量间隔内释放弹性应变能的次数。3.根据权利要求1所述的一种应变场声发射事件时空强的测量方法，其特征在于:所述 步骤8具体步骤如下： 步骤8.1、在测量计间隔内，将各声发射事件释放的弹性应变能按最大值和最小值为首 末端均分为能量释放若干等级、将各声发射事件的空间尺度按最大值和最小值为首末端均 分为空间尺度若干等级、将各声发射事件的持续时间按最大值和最小值为首末端均分为持 续时间若干等级； 步骤8.2、将各声发射事件释放的弹性应变能、各声发射事件的空间尺度、各声发射事 件的持续时间分别划归入相应的等级中，并确定能量释放各等级的频次、空间尺度各等级 的频次、持续时间各等级的频次，所述能量释放各等级的频次指能量释放各等级包含的相 应的声发射事件的个数，所述空间尺度各等级的频次指空间尺度各等级包含的相应的声发 射事件的个数，所述持续时间各等级的频次指持续时间各等级包含的相应的声发射事件的 个数； 步骤8.3、利用能量释放各等级的数据，通过线性回归确定频次-能量释放关系，利用空 间尺度各等级的数据，通过线性回归确定频次-空间尺度关系，利用持续时间各等级的数 据，通过线性回归确定频次-持续时间关系。
【文档编号】G01N29/14GK105929027SQ201610206423
【公开日】2016年9月7日
【申请日】2016年3月31日
【发明人】王学滨, 郭翔, 冯威武, 白雪元, 马冰, 齐大雷, 李阳, 武其奡
【申请人】辽宁工程技术大学