[发明专利]一种基于离散单元法的声发射模拟方法及系统

权利要求书

1.一种基于离散单元法的声发射模拟方法，其特征在于，所述模拟方法包括：

在计算模型表面布置监测点；

监测所述监测点的速度波形；

对所述速度波形进行分析，得到声发射数、能量以及b值；

对所述速度波形进行分析，得到声发射数具体包括：

A1：导入监测点的速度波形；

A2：依次对每个采样数据进行分析；

A3：判断第i个数据是否小于或等于采样总数据数量；

A4：若是，则结束；

A5：若否，则执行下一步；

A6：判断第i个数据是否大于或等于预设的门槛值；

A7：若否，则i＝i+1，返回步骤A2，对下一个数据进行判断；

A8：若是，则执行下一步骤；

A9：判断一次声发射事件是否结束；

A10：若否，则强行切断，该声发射事件结束，执行下一步骤；

A11：若是，则记录本次声发射事件的特征参数，所述特征参数包括到达时间、幅值、信号持续时间、振铃计数；

A12：对下一个声发射信号进行分析，直到所有采样数据被分析完毕；

对所述速度波形进行分析，得到能量具体包括：

B1：判断第i个数据是否小于或等于采样总数据数量；

B2：若否，则结束；

B3：若是，则计算E(i)＝[V_x²(i)+V_y²(i)]；其中，E(i)表示第i个数据的能量，V_x指速度的x方向分量，V_y指速度的y方向分量；

B4：计算其中，M(i)表示第i个数据的震级；

B5：i＝i+1，返回步骤B1，直到所有采样数据分析完毕；

对所述速度波形进行分析，得到b值具体包括：

C1：输入震级范围[M_min,M_max]，并输入M＝M_min；

C2：判断第i个数据是否小于或等于采样总数据数量；

C3：若否，则结束；

C4：若是，则执行下一步骤；

C5：获取M(1)～M(i)中大于或等于M_min的数据个数，计为N1，获取M(1)～M(i)中大于或等于M_min+0.1的数据个数，计为N2，直到获取完M(1)～M(i)中大于或等于M_max的数据个数，记为Nm；其中，M(1)表示第1个数据的震级，M(i)表示第i个数据的震级；

C6：得到m组(M_min～_max)～N_(1～m)数据，根据公式logN＝a-bM进行拟合得b(i)值，i＝i+1，返回步骤C2；其中a，b为常数。

2.一种基于离散单元法的声发射模拟系统，其特征在于，所述模拟系统包括：

监测点布置模块，用于在计算模型表面布置监测点；

速度波形监测模块，用于监测所述监测点的速度波形；

速度波形分析模块，用于对所述速度波形进行分析，得到声发射数、能量以及b值；

所述速度波形分析模块包括：声发射数分析单元、能量分析单元以及b值分析单元；

所述速度波形分析模块的声发射数分析单元具体包括以下步骤：

A1：导入监测点的速度波形；

A2：依次对每个采样数据进行分析；

A3：判断第i个数据是否小于或等于采样总数据数量；

A4：若是，则结束；

A5：若否，则执行下一步；

A6：判断第i个数据是否大于或等于预设的门槛值；

A7：若否，则i＝i+1，返回步骤A2，对下一个数据进行判断；

A8：若是，则执行下一步骤；

A9：判断一次声发射事件是否结束；

A10：若否，则强行切断，该声发射事件结束，执行下一步骤；

A11：若是，则记录声发射特征参数与到达时间、幅值、信号持续时间、振铃计数；

A12：对下一个声发射信号进行分析，直到所有采样数据被分析完毕；

所述能量分析单元具体包括以下步骤：

B1：判断第i个数据是否小于或等于采样总数据数量；

B2：若否，则结束；

B3：若是，则计算E(i)＝[V_x²(i)+V_y²(i)]；

B4：计算

B5：i＝i+1，返回步骤B1，直到所有采样数据分析完毕；

所述b值分析单元具体包括以下步骤：

C1：输入震级范围[M_min,M_max]，并输入M＝M_min；

C2：判断第i个数据是否小于或等于采样总数据数量；

C3：若否，则结束；

C4：若是，则执行下一步骤；

C5：获取M(1)～M(i)中大于或等于M_min的数据个数，计为N1，获取M1～Mi中大于或等于M_min+0.1的数据个数，计为N2，直到获取完M(1)～M(i)中大于或等于M_max的数据个数，记为Nm；其中，M(1)表示第1个数据的震级，M(i)表示第i个数据的震级；

C6：得到m组(M_min～_max)～N_(1～m)数据，根据公式logN＝a-bM进行拟合得b(i)值，i＝i+1，返回步骤C2；其中a，b为常数。