蒙特卡罗法在Matlab环境.docVIP

作者：朱川曲 张道兵 朱海燕 【摘要】 应用工程结构可靠性理论，建立煤巷锚杆支护结构可靠性模型，针对其结构稳定性的极限状态功能函数高度非线性的特征，采用蒙特卡罗法在Matlab环境下直接产生服从各相应概率分布函数的随机变量数组计算其结构可靠度，编程过程简化，计算速度快，精度高，且不受极限状态方程非线性、随机变量非正态的限制，开辟了煤巷锚杆支护结构可靠度计算的新途径。建立了支护参数与可靠度的关系，指出了提高煤巷锚杆支护结构可靠度的方法和措施，为煤巷锚杆支护参数设计和优化提供了科学依据。【关键词】 煤巷；锚杆支护结构；可靠性；蒙特卡罗法；Matlab 0 引 言 煤巷锚杆支护的工程实践表明，巷道支护失效主要是因为锚杆支护参数选择不合理造成的。究其原因，主要有两个方面：一是巷道围岩的力学特性存在很大差异，许多参数为随机变量，以往用确定性模型进行锚杆支护参数选择，容易出现误差；二是在进行煤巷锚杆支护参数选择时，往往把载荷效应作为一个定值，再给一定的安全系数。巷道支护的工程实践表明，其载荷效应具有不确定性，即随机性。当载荷效应的变异性较大时，很可能使支护结构达不到设计的安全水平，从而处于不安全状态。因此，这就需要在充分考虑参数随机性的可靠性理论指导下进行巷道支护设计。 然而，煤巷锚杆支护结构的极限状态功能函数是高度非线性的，如果采用求导数的可靠度计算方法，极限状态功能函数变得很难处理，且计算结果存在一定的误差。而蒙特卡罗法在目前可靠度计算中，被认为是一种相对精确的方法，且它不受极限状态方程非线性、随机变量非正态的限制。如果采用Fortran，Basic或C语言等编制计算程序实现蒙特卡罗法计算结构可靠度，其编程复杂，效率低下，而且还需要系统掌握语言的算法和编程技巧，对广大工程技术人员来说，就有相当的难度。而在Matlab中采用蒙特卡罗直接抽样法计算结构可靠度，能很好地解决上述问题。 因此，笔者根据煤巷围岩稳定性极限状态方程，利用Matlab的强大数值计算功能，建立可靠度计算模型，实现在Matlab中采用蒙特卡罗直接抽样法计算结构可靠度，并对煤巷锚杆支护结构的可靠性进行系统的分析。 1 蒙特卡罗法计算工程结构可靠度与Matlab实现 1.1 基本原理 蒙特卡罗法的基本原理为：根据大数定理，设x1，x1，…，xn是n个独立的随机变量，若它们来自同一母体，有相同的分布，且具有相同的有限均值和方差，分别用μ和σ2表示，则对于任意ε0有  另外，若随机事件A发生的概率为P(A)，在n次独立试验中，事件A发生的频数为m，频率为W(A)=m/n，则对于任意ε0有  蒙特卡罗法的基本思路是先对影响结构可靠度的随机变量进行大量随机抽样，然后把所抽样值一组一组地代入结构功能函数式，确定结构失效与否，最后从中求得结构的可靠度。  1.2 随机数产生 Matlab实现蒙特卡罗法的最突出的优点是它拥有功能强大的随机数发生器指令。通常在用蒙特卡罗直接抽样法时，必须解  以求得服从相应分布类型的随机变量x，再代入功能函数求解。而Matlab提供了23种随机变量分布类型的随机数发生器，如正态分布、对数正态分布、泊松分布、威布尔分布等，可以简单方便地产生工程实际中常出现的变量分布的随机数，直接代入功能函数，省去了可能会带来很大麻烦的求分布函数反函数这一步，极大地提高了效率。  1.3 求解步骤 设功能函数为g(x1，x2，…，xn)，xi(i=1，2，…，n)为服从各相应概率分布函数的随机变量。首先用相应随机数发生器指令产生m×n的随机变量数组xi(k，l)(i=1，2，…，n；k=1，2，…，m；l=1，2，…，n)，然后将各数组中元素一一对应代入功能函数，最后得到功能函数结果数组。统计结果数组中不小于0的元素个数j，可得可靠度：  Matlab实现蒙特卡罗法计算可靠度程序框图如图1所示。 2 煤巷锚杆支护结构的可靠度 2.1 顶板锚杆支护结构的可靠度 巷道顶板的破坏类型主要有两种：一是沿拱座处(巷道两顶角上方)发生剪切破坏；二是巷道顶板强度较低时，在拱座处和巷道直接顶因受挤压而超过岩石的残余强度时发生压缩破坏。  图1 Matlab实现蒙特卡罗法计算可靠度程序框图 顶板不发生剪切破坏时的极限状态方程为  顶板不发生剪切破坏时的功能函数为  式中， N——顶角锚杆处于冒落范围以外的锚固力； W—一潜在冒落块体的重量； H——冒落块高度； σt——潜在冒落范围以内最大水平应力的稳定值；