2025《Matlab建立随机节理网络的案例分析》1800字.docxVIP

Matlab建立随机节理网络的案例分析

1岩体节理分布特征

结合上文对于岩体节理概率分布的研究和现场地质勘查数据，坡体节理面迹长和倾角的几何参数分布规律如图5-1至图5-3，经过函数曲线拟合并根据皮尔森相关系数进行检验。

（1）J1倾角与迹长概率模型

图5-1J1倾角与迹长概率分布图

通过软件进行拟合，可以发现J1倾角基本符合正态分布规律，J1迹长基本符合对数正态分布规律，其概率密度分布函数分别为：

f

f

J1倾角和J1迹长的拟合曲线的皮尔森相关系数的平方分别为0.9340和0.9636，其X与Y相关性均较好。

（2）J2倾角与迹长概率模型

图5-2J2倾角与迹长概率分布图

通过软件进行拟合，可以发现J2倾角基本符合对数正态分布规律，J2迹长基本符合正态分布规律，其概率密度分布函数分别为：

f

f

J2倾角和J2迹长的拟合曲线的皮尔森相关系数的平方分别为0.9631和0.9455，其X与Y相关性均较好。

（3）J3倾角与迹长概率模型

图5-3J3倾角与迹长概率分布图

通过软件进行拟合，可以发现J3倾角基本符合对数正态分布规律，J3迹长基本符合正态分布规律，其概率密度分布函数分别为：

f

f

J3倾角拟合曲线的皮尔森相关系数的平方为0.9223，其X与Y相关性较好。J3迹长拟合曲线的皮尔森相关系数的平方为0.0357，其X与Y相关性较差，其概率密度分布函数为fx

根据上述曲线拟合得出J1、J2和J3的倾角和迹长概率分布类型及其均值和标准差，见表5-2。

表5-2蒙特卡洛模拟岩体结构面几何特征参数

组数

倾角/°

迹长/m

间距/m

均值

标准差

分布类型

均值

标准差

分布类型

均值

1

75.18

5.95

正态

11.32

3.68

对数正态

1.69

2

32.78

6.30

对数正态

4.55

1.38

正态

0.77

3

49.92

7.07

对数正态

6.94

2.87

平均

2.26

2随机节理边坡模型在Matlab中的实现

根据第三章测线法实地勘察节理数据并经过概率模型拟合得知，上木江坪边坡的节理裂隙是按照一定概率模型随机分布的，为了实现二维随机节理网络的模拟，在Matlab编程中做出以下几个设定：

（1）本文根据UDEC软件的模型要求，将节理裂隙模拟为二维平面，以线段来模拟节理裂隙，以倾角来表示节理裂隙产状。

（2）节理线段由以下3个数据描述：节理线段中点坐标x0,y

x=

（3）模拟节理裂隙的中点分布在生成域内服从均匀分布，从而推算出节理裂隙的条数为节理生成域的面积除以节理平均间距。

为了使生成节理程序和模型能在UDEC程序中运行，需要生成每条节理的起始和终点坐标，以满足“crack”命令的条件。上木江坪边坡具有3组优势节理，将该程序进行3次模拟，即可得到n组优势节理组成的随机节理网络。根据5.2统计并计算各节理参数得出的概率分布型，生成节理中点坐标、迹长和倾角的随机数。由此即可得知每条节理在生成域中的起点和终点坐标，由此可在生成域内描绘出每一条节理。

运用Matlab软件生成随机节理边坡模型的过程如下：

第一步，生成域参数设置编程：

输入生成生成域的长和宽，确定生成域的面积；

第二步，节理几何参数模拟编程：

输入某组节理间距均值，取其均值，并且假设节理中心点呈均匀分布，以确定生成域内的节理条数；接下来输入节理迹长和倾角参数，首先输入迹长和倾角的概率分布类型（正态分布、对数正态分布、负指数分布、均匀分布），再输入该概率分布类型下的均值和标准差，然后得到某组节理在此概率分布类型下生成的随机节理，即节理线段的起点和终点坐标；

第三步，节理坐标提取及矩阵生成

通过蒙特卡洛随机模拟得到某组节理中心点服从均匀分布的坐标，将得到的中心点坐标编辑成节理中心点矩阵，然后将第三步生成的节理几何参数即迹长和倾角赋予节理中心点矩阵，即生成了一个含有每条节理起始和终点坐标的随机节理矩阵；

第四步，绘图并保存图像：

将每组节理按照以上步骤输入参数并计算，输入完三组节理后，按得到的坐标矩阵绘图，得到随机节理模型，如图5-4所示；

图5-4Matlab随机结构面模拟

第五步，形成边坡模型：

根据边坡高程数据生成边坡模型，将矩形生成域按照边坡拟建铁路隧道入口处剖面切割，剖面线ab如图5-5（a）所示，边坡模型高60m，宽60m。将上述节理数据输入UDEC程序，由于UDEC程序计算时认为边坡沿着贯通的裂隙破坏，所以运用“jdelete”删除多余部分节理，使得只有贯通岩体的裂隙保留在模型中，而未贯通裂隙的区域保留完整岩块的形式，如图5-6（b）所示。

（a）模型剖面线（b）随机节理边坡模型

图5-6某拟建铁路隧道洞口随机节理边坡模型