基于岩体结构面随机强度的非连续变形分析方法技术

本发明专利技术提供了一种基于岩体结构面随机强度的非连续变形分析方法。首先，对岩体中的结构面进行分类，通过相关力学试验获取不同类型的结构面对应的强度参数，包括凝聚力、摩擦角及单轴抗拉强度的均值和标准差；然后，根据各类型结构面强度参数的均值和标准差，利用坐标变换法生成相应结构面的随机强度；最后，基于每一条结构面的随机强度对岩体进行非连续变形分析模拟，包括依据随机强度计算滑动摩擦力，以及依据随机强度进行接触判断、求解开闭状态等。依本发明专利技术的方法，结构面强度的取值更接近岩体的实际情况，非连续变形分析的计算结果更为准确、真实。

Discontinuous deformation analysis method based on random strength of rock mass discontinuities

The present invention provides a discontinuous deformation analysis method based on random strength of rock mass structure plane. First of all, the classification of structural plane of rock mass, the mechanical test for different types of structures should face strength parameters, including mean and standard cohesion, friction angle and the uniaxial tensile strength difference; then, according to the mean and standard strength of various types of structural plane parameters, the corresponding structure of the stochastic intensity generated by the coordinate transformation method; finally, a structure of the random strength of rock mass by non continuous deformation based on simulation analysis, including based on random strength calculation of sliding friction force, and according to the random strength contact judgment, solving the opening and closing state etc.. According to the method of the invention, the strength of the structural surface is closer to the actual situation of the rock mass, and the calculation result of the discontinuous deformation analysis is more accurate and true.

【技术实现步骤摘要】
基于岩体结构面随机强度的非连续变形分析方法
本专利技术涉及一种基于岩体结构面随机强度的非连续变形分析方法，属于水利水电工程

技术介绍
在岩体工程的稳定性分析中，各种分析方法日趋成熟，经典的极限类分析方法和近代数值类分析方法得到的结果基本一致，并能相互印证，方法本身不再是岩体稳定性分析的制约因素。影响分析结果的主要因素是岩体的强度参数。对于强度参数的取值，传统上采用确定性方法，通过地质调查和室内试验对岩体中的结构面进行分类定值，对每种类型的结构面给出一组确定的强度参数。在数值模拟时，对属于同一类型的结构面全部采用同一组强度参数进行计算。这种取值方法简单明了易于理解，但与实际情况不符，得到的分析结果存在偏差。实际上，岩体结构面的力学性能受到含水率，夹层性态，产状，间距等多种因素的影响，其强度参数在空间分布上具有较大的离散性和随机性。即便对属于同一类型的结构面来说，每条结构面的强度也并不是定值，而是服从某种随机分布规律。一般认为，结构面的强度参数服从正态分布，强度参数围绕某一均值左右浮动。在数值模拟时，应当根据同一类型结构面的多个样本的统计规律，对每条结构面的强度参数进行随机取值，依据随机强度进行计算分析，从而取得更加接近实际的结果。非连续变形分析(DiscontinuousDeformationAnalysis简称DDA)方法是研究包括岩体在内的非连续介质的有力工具。该方法自20世纪80年代提出以来，受到了岩土工程界学者的广泛关注，陆续在国内外大中型岩体工程中得到了应用，其有效性已经被大量验证。DDA方法的研究对象是离散的块体系统，其优点在于，能够充分考虑块体与块体间的相互作用，通过在块体界面上增、减接触弹簧，能正确模拟块体与块体间接触面的闭合，张开和剪切滑移。但是，现有的DDA方法仍然采用结构面的定值强度进行计算分析，虽然简单明了，但与岩体的实际情况不符，计算结果往往存在偏差，不能很好的反映实际。
技术实现思路
鉴于上述原因，本专利技术的目的在于提供一种基于岩体结构面随机强度的非连续变形分析方法，对岩体中的每一条结构面生成随机强度，基于每一条结构面的随机强度进行DDA模拟，这种强度参数的取值方法更接近岩体的实际情况，计算及分析结果更为准确、真实。为实现上述目的，本专利技术采用以下技术方案：一种基于岩体结构面随机强度的非连续变形分析方法，包括：获取岩体中不同类型结构面的强度参数对应的均值和标准差；根据各类型结构面的强度参数对应的均值和标准差，生成各类型结构面对应的每条结构面的随机强度；基于每条结构面的随机强度进行非连续变形分析模拟。进一步的，基于结构面的随机强度进行非连续变形分析的方法为：根据第j条结构面的随机强度按照公式(1)计算该结构面上的摩擦力F，其中，p为法向接触刚度，d为法向嵌入距离，s为摩擦力的方向函数，为第j条结构面的随机摩擦角；然后，根据计算得到的摩擦力生成等效荷载列阵{F}，并集成到总体平衡方程进行求解。基于结构面的随机强度进行非连续变形分析的方法为：根据结构面的随机强度进行接触判断，设Nj为第j条结构面上的法向接触力，Tj为第j条结构面上的切向接触力，为第j条结构面上的切向位移矢量，为摩擦力矢量，||表示两个矢量指向同一方向，为第j条结构面的随机强度，则：A、当前的接触状态为张开时，若Nj＞ftj，则接触状态保持不变；B、当前的接触状态为张开时，若Nj＜0，且则将接触状态更新为滑动；C、当前的接触状态为张开时，若Nj＜0，且则将接触状态更新为锁定；D、当前的接触状态为滑动时，若Nj＞ftj，则将接触状态更新为张开；E、当前的接触状态为滑动时，若Nj＜0且则接触状态保持不变；F、当前的接触状态为滑动时，若Nj＜0且刚将接触状态更新为锁定；G、当前的接触状态为锁定时，若Nj＞ftj，则将接触状态更新为张开；H、当前的接触状态为锁定时，若Nj＜0，且则将接触状态更新为滑动；I、当前的接触状态为锁定时，若Nj＜0，且则接触状态保持不变。通过力学试验获取结构面的凝聚力c的均值与标准差σc，摩擦角的均值与标准差单轴抗拉强度ft的均值和标准差σf。生成结构面的随机强度的方法是：对于第i种类型的第j条结构面，利用随机函数生成在(0，1)区间内均匀分布的随机数；利用坐标变换法对均匀分布的随机数进行处理，生成服从标准正态分布的随机数；根据第i种类型结构面对应的强度参数的均值和标准差，按照公式(2)生成第j条结构面的随机凝聚力随机摩擦角随机单轴抗拉强度其中，x1为服从正态分布的随机数。利用随机函数生成在(0，1)区间内均匀分布的随机数r1、r2，利用坐标变换法将随机数r1、r2处理生成服从标准正态分布的随机数x1、x2，交替使用随机数x1、x2生成结构面的随机强度。本专利技术的优点是：本专利技术的基于岩体结构面随机强度的非连续变形分析方法，对岩体中的每一条结构面生成服从正态分布的随机强度，基于每一条结构面的随机强度进行非连续变形分析，包括依据随机强度计算摩擦力，依据随机强度进行接触判断、求解开闭状态等。这种结构面强度的取值方法更接近岩体的实际情况，计算及分析结果更为准确、真实。附图说明图1是包含所有岩体结构面的DDA分析模型示意图。图2A是岩体结构面的凝聚力正态分布示意图。图2B是岩体结构面的摩擦角正态分布示意图。图2C是岩体结构面的单轴抗拉强度正态分布示意图。图3是基于岩体结构面随机强度进行DDA分析的方法流程图。具体实施方式以下结合附图和实施例对本专利技术作进一步详细的描述。本专利技术公开的基于岩体结构面随机强度的非连续变形分析方法，包括以下步骤：S1：获取岩体中不同类型结构面的强度参数对应的均值和标准差；岩体中共包含m条结构面，被分为n种类型，根据室内相关力学试验的结果，通过回归分析，获得岩体中n种类型结构面的强度参数对应的均值和标准差；具体包括：S11：对岩体中第i种类型的结构面进行多组抗剪强度试验，得到多组剪断面上的剪应力τ和正应力σ的样本值(τ1，σ1)，(τ2，σ2)……(τk，σk)；对得到的多组样本值进行回归分析，得到第i种类型结构面的凝聚力c的均值与标准差σc，摩擦角的均值与标准差同理，对上述类型的结构面进行多组单轴抗拉强度试验，对得到的多组试验数据样本值进行回归分析，得到第i种类型结构面的单轴抗拉强度ft的均值和标准差σf。S12：重复步骤S11，获得岩体中n种类型结构面对应的n组强度参数分别对应的均值和标准差，即：其中，i＝1，2，3……n。S2：根据岩体的地质剖面信息、几何信息生成包含所有结构面的DDA分析模型；模型如图1所示，根据岩体的地质剖面信息、几何信息生成包含所有结构面的DDA分析模型，属于现有技术，本专利技术不再详细描述。S3：将DDA计算控制参数及岩体中不同类型结构面的强度参数对应的均值和标准差载入DDA分析模型；该步骤可参见图3，DDA计算控制参数包括时间步长、位移比、接触刚度、动力系数等参数，其为DDA方法中的公知技术，本专利技术不再详细描述。S4：根据各类型结构面的强度参数对应的均值和标准差，生成每一条结构面的随机强度；以属于第i种类型的第j条结构面为例，说明利用第i种类型结构面的强度参数对应的均值和标准差生成第j条结构面的随机强度的过程，具体包括：S41：利用随机函数生成均匀分本文档来自技高网...

【技术保护点】
基于岩体结构面随机强度的非连续变形分析方法，其特征在于，包括：获取岩体中不同类型结构面的强度参数对应的均值和标准差；根据各类型结构面的强度参数对应的均值和标准差，生成各类型结构面对应的每条结构面的随机强度；基于每条结构面的随机强度进行非连续变形分析模拟。

【技术特征摘要】
1.基于岩体结构面随机强度的非连续变形分析方法，其特征在于，包括：获取岩体中不同类型结构面的强度参数对应的均值和标准差；根据各类型结构面的强度参数对应的均值和标准差，生成各类型结构面对应的每条结构面的随机强度；基于每条结构面的随机强度进行非连续变形分析模拟。2.根据权利要求1所述的基于岩体结构面随机强度的非连续变形分析方法，其特征在于，基于结构面的随机强度进行非连续变形分析的方法为：根据第j条结构面的随机强度按照公式(1)计算该结构面上的摩擦力F，其中，p为法向接触刚度，d为法向嵌入距离，s为摩擦力的方向函数，为第j条结构面的随机摩擦角；然后，根据计算得到的摩擦力生成等效荷载列阵{F}，并集成到总体平衡方程进行求解。3.根据权利要求2所述的基于岩体结构面随机强度的非连续变形分析方法，其特征在于，基于结构面的随机强度进行非连续变形分析的方法为：根据结构面的随机强度进行接触判断，设Nj为第j条结构面上的法向接触力，Tj为第j条结构面上的切向接触力，为第j条结构面上的切向位移矢量，为摩擦力矢量，||表示两个矢量指向同一方向，为第j条结构面的随机强度，则：A、当前的接触状态为张开时，若Nj＞ftj，则接触状态保持不变；B、当前的接触状态为张开时，若Nj＜0，且则将接触状态更新为滑动；C、当前的接触状态为张开时，若Nj＜0，且则将接触状态更新为锁定；D、当前的接触状态为滑动时，若Nj＞ftj，则将...

【专利技术属性】
技术研发人员：张国新，雷峥琦，程恒，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11