[发明专利]一种花岗岩细观裂纹结构的数值表征方法

说明书

本发明为解决目前针对花岗岩数值模拟的研究不能兼顾考虑花岗岩中材料非均匀性、晶内裂纹、晶界裂纹影响的问题，提供了一种花岗岩细观裂纹结构的数值表征方法。包括以下步骤：1)、随机生成晶粒；2)、赋予晶粒材料编号；3)、随机选取若干个晶粒并裁剪其晶界，生成晶界裂纹；4)、利用数值模拟软件给裁剪后的晶粒赋予各晶粒的材料编号对应的材料模型，再通过布尔运算将裁剪后的晶粒组合生成包含晶界裂纹的花岗岩模型；5)、设定的晶内裂纹长度和取向，随机生成晶内裂纹，并分割晶粒；6)、通过草图切割的方式，根据晶内裂纹分割花岗岩模型，并赋予晶内裂纹属性。7)、建立裂纹接触模型，得到花岗岩细观裂纹宽度的数值表征。

技术领域

本发明属于岩土工程声波探测领域，具体涉及一种花岗岩细观裂纹结构的数值表征方法。

背景技术

花岗岩作为一种常见的地质材料，其力学特性对地质工程、岩土爆破工程有着重要影响，而声速是其力学特性的重要体现。

花岗岩宏观上常表现为连续的、致密的岩石，但内部却包含大量微观结构。材料组成上，花岗岩内部含有石英、长石、云母等不同晶粒，在晶粒边界和晶粒内部含有大量微裂纹。不同的花岗岩，其各类晶粒含量不同，晶界裂纹和晶内裂纹的含量以及宽度亦不相同。声波在花岗岩等岩石介质中传播时，遇到裂纹等微结构时会发生透射、反射、绕射等现象，进而对宏观声速、声透射系数产生影响，这也是通过声波对岩石品质进行评价的依据。建立包含内部微裂纹的花岗岩细观数值模型，通过数值模拟的方式对声波在花岗岩中的传播机制进行深入分析，可对基于声波探测的工程质量评价提供依据。

目前的数值模拟研究常将花岗岩当做连续介质或仅考虑材料非均匀性的多晶材料处理，尚缺乏可以同时考虑花岗岩中材料非均匀性、晶内裂纹、晶界裂纹的数值建模方法。

发明内容

本发明的目的是解决目前针对花岗岩数值模拟的研究不能兼顾考虑花岗岩中材料非均匀性、晶内裂纹、晶界裂纹影响的技术问题，而提供一种花岗岩细观裂纹结构的数值表征方法。

为达到上述目的，本发明采用的技术方案为：

一种花岗岩细观裂纹结构的数值表征方法，其特殊之处在于，包括以下步骤：步骤1)、随机生成晶粒

在目标区域内随机布设与待生成的花岗岩晶粒数目相等的n个种子，采用泰森多边形分割法将目标区域划分成n个多边形子区域，每个多边形子区域看作一个晶粒，每个晶粒中包含1个种子，n个晶粒的序列编号分别为S₁,S₂,S₃,…,S_k…,S_n，n≥k≥1，且n为正整数；

步骤2)、赋予晶粒材料编号

将n个晶粒分别赋予材料编号M_i，不同的编号代表不同的晶粒种类；其中，i＝1,2,3,...,m，m为晶粒的种类数，m≥2；

步骤3)、随机选取若干个晶粒并裁剪其晶界，生成宽度为d_deta的晶界裂纹；

步骤4)、利用数值模拟软件前处理模块给步骤3)中裁剪后的晶粒赋予步骤2)中各晶粒的材料编号对应的材料模型，再通过布尔运算将裁剪后的晶粒组合生成包含晶界裂纹的花岗岩模型；

步骤5)、设定晶内裂纹长度和取向，在步骤3)裁剪后的晶粒中随机生成宽度为d_inner的晶内裂纹，利用晶内裂纹线段分割晶粒，并建立晶内裂纹集合；

步骤6)、通过草图切割的方式，根据步骤5)晶内裂纹集合中的晶内裂纹对步骤4)所得花岗岩模型进行分割，并赋予晶内裂纹在数值模拟软件中对应的裂纹属性；

步骤7)、建立裂纹接触模型，得到花岗岩细观裂纹宽度的数值表征：d_outer＝d_inner+d_deta。