基于逼近阶跃函数和拉格朗日插值的改进DDA方法

针对非连续变形分析中开合迭代难以收敛的难题,基于块体接触约束状态和块体位移之间的关系,提出了基于逼近阶跃函数和拉格朗日插值的改进DDA方法。采用双曲正切函数来逼近阶跃函数,利用阶跃函数将块体接触约束状态用块体位移来表达,以此来替代开合迭代,避免了开合迭代难以收敛的难题。利用拉格朗日插值原理,推导得到只含有块体位移为未知量的块体系统势能函数,并利用变尺度法来求解总体势能函数的极值以得到块体位移。分别结合滑块模型和地下洞室模型,分析了改进DDA方法的计算精度和计算速度,验证了文中提出的改进DDA方法的正确性和稳定性。研究表明:基于逼近阶跃函数和拉格朗日插值的改进DDA方法具有较高的精度,且相比较传统DDA方法而言,具有更为稳定的和更为强健的计算收敛性。因此,基于逼近阶跃函数和拉格朗日插值的改进DDA方法是一种稳定有效的数值计算方法,为解决非连续变形中开合迭代难以收敛的问题提供了新思路。     

一种改进步长自调的非连续变形分析法

非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径     

非连续变形分析的精度改进方法及研究趋势

由石根华提出的非连续变形分析(DDA)方法理论严密且计算较为高效,在模拟块体滑移、转动和张开等大变形、大位移问题方面具备独特优势,因而在岩土工程领域得到了广泛运用。然而DDA方法在提出初期不可避免地存在一些精度问题,对此石根华及国内外学者就问题产生原因和改进方法开展了大量研究工作。基于部分相关成果的研究学习,对改进原始DDA计算精度的各类方法进行归纳探讨,主要包括:(1)DDA块体内部应力位移场的精度控制;(2)块体接触问题处理方法的改进;(3)人为参数的合理选取;(4)能量耗散机制的考虑;(5)人工边界的改进等方面。对以上方法的改进效果和计算效率进行了简要分析和讨论,在此基础上对DDA的研究热点和发展趋势进行了概述,为该方法的进一步发展和完善提供思路。     

岩石力学联锁咬合岩块结构稳定性分析

首次建立了岩石力学中的多种联锁咬合岩块结构模型,介绍了块体理论的进展,并采用改进的DDA（块体不连续变形分析）法对块体进行了稳定性分析。     

公路隧道局部塌方洞段的围岩稳定性评价

结合某公路隧道工程,以现场地质调查和监测资料为基础,在结构面发育且可能发生大规模塌方的围岩洞段,根据现场实际统计的优势结构面产状,采用关键块体理论、非连续变形方法（DDA）首先确定围岩主要的失稳破坏模式,然后运用连续介质分析程序FLAC,在DDA方法确定的围岩实际变形破坏塌方形态的基础上进行计算分析,评价围岩的整体稳定性及加固措施的有效性。分析表明,围岩破坏的主要模式为受结构面控制的局部块体的失稳,工程开挖后围岩不会发生大变形,这与实际情况基本一致,为安全施工提供了保证。采用非连续变形方法和连续介质计算相结合的分析方法对隧道工程塌方稳定性进行研究,取得了较好的效果,研究方法和经验可为类似工程参考并借鉴。     

数值流形方法框架下散体系统与连续介质共同作用模拟

数值流形方法是包含流形元、有限元及DDA在内的数值方法体系,建立流形元与DDA块体的接触方程,则可实现流形方法框架下的连续介质和散体系统共同作用模拟。针对填石路堤工程,编制了大型数值计算程序,采用块体随机生成、块体粒径控制及块体自然堆积的方法建立散体系统的DDA模型,对路堤的分层铺设、碾压及工后沉降变形等进行模拟分析。通过算例表明,在数值流形方法框架下,采用流形元与DDA共同作用的方法,可以很好地对同时存在连续变形和散体大变形的体系进行计算分析,其对该类问题的模拟更接近分析对象的实际情况,有助于从根本上揭示分析对象变形的细观机制和规律,并能考察更多因素对工程问题的影响。     

基于离散元的强度折减法分析岩质边坡稳定性

将通用离散元UDEC与强度折减法结合,对含多结构面的岩质边坡的稳定性进行了分析。通过对节理岩质边坡的UDEC模型中的可变形块体和节理单元的强度参数进行折减,使模型不能再达到平衡状态,此时的折减系数就是边坡的安全系数,另外,由对应的边坡块体的速度矢量可以确定滑动面和边坡的破坏形态。通过与传统的条分法的结果比较,表明基于UDEC的强度折减法是一种可靠、有效的方法,为复杂节理岩质边坡的滑动面确定与安全系数计算开辟了新的途径。     

节理岩体爆破的DDA方法模拟

在非连续变形分析（DDA）方法中,通过跟踪炮孔扩张和炮孔周边裂隙的发展贯通求得爆腔的即时体积,进而根据爆生压力状态方程计算爆腔即时压力,并将爆生压力载荷作用到主炮孔内壁和贯通裂隙面上,实现了爆生产物作用下节理岩体爆破的DDA方法模拟。采用DDA方法模拟了节理岩体中的水平柱状炮孔抛掷爆破问题,得到了爆腔的体积扩张和压力衰减时间曲线,模拟很好的再现了岩石爆破过程中的炮孔扩张、岩体破坏、块体抛掷和爆堆形成过程。     

基于DDA方法一种流-固耦合模型的建立及裂隙体渗流场分析和应用

以非连续变形分析方法（DDA）为基础并采用稳态流体计算方法将二者结合进行裂隙岩体流-固耦合分析。利用DDA方法生成裂隙岩体模型,在此基础上采用矩阵搜索等方法形成新的裂隙水通网络模型。采用稳态迭代算法和立方定律求得裂隙水压力,并把裂隙水压力作为线载荷施加到块体边界,在DDA算法中每个迭代步完成后更新裂隙开度和水压值,与DDA算法结合研究裂隙水与块体之间相互作用关系。利用以上裂隙岩体流-固耦合计算方法研究了某水封油库开挖和运行过程洞室围岩流量和密封性,为该工程预测水封效果提供了有益的主要依据,也是国内首次采用DDA方法做大型工程的流-固耦合模型分析。     

混合多位移模式的非连续变形分析法研究

传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。     

Boundary Element based Discontinuous Deformation Analysis

A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd.     

节理对爆炸波传播影响的数值研究

采用加入无反射边界条件的DDA程序,研究了节理面对应力波传播的影响。结果表明,节理面能阻碍波的传播,有利于波的衰减,节理面越多,波的反射越强,而波的透射越弱。模拟了一个现场爆炸试验,研究爆炸产生的应力波在节理岩体中传播、衰减的规律,模拟结果与现场试验结果比较吻合。研究表明,DDA方法可以模拟节理面对应力波传播的阻碍作用,用它来模拟爆炸波在节理岩体中的传播是适用的。     

Discontinuous deformation analysis of super section tunnel surrounding rock stability based on joint distribution simulation

Based on binocular stereo photography measurement technology, the actual distribution information of joints and fractures in the tunnel face were obtained via image processing and feature extraction, and the preliminary evaluation of surrounding rock stability of Laohushan tunnel was conducted according to surrounding rock classification method. Since all available surrounding rock classification approaches didn’t consider the influences of the size effect of tunnel excavation span and unfavorable geologic bodies such as weak-fracture zones, an improved discontinuous deformation analysis (DDARF) method was adopted to conduct a numerical simulation of the deformation and fracture processes of the surrounding rock. Using the traveling wave method, triangle DDA blocks were automatically generated in the calculation zone, and the block boundaries were divided into real joints and virtual joints. Based on the real joint information obtained via aforementioned photographic measurement, the real joints in the tunnel face were dynamically modified, in order to achieve the simulation of joint distribution. The results revealed that the fracture evolution regularity, deformation failure mechanism, and block dropping phenomenon that the DDARF joint simulation model calculated are in good agreement with actual conditions, while those obtained based on conventional models present differences from field situation. Additionally, focusing on the localised rockfall phenomenon of Laohushan tunnel, the crack extension rate was introduced to conduct a quantitative comparison of the rock crack evolution process with or without anchor supporting. The research results offer practical guidance for field construction and anchorage support scheme optimization.     

Study on elastic P-wave propagation law in unfavorable geologic structures with discontinuous deformation analysis method

The propagation of elastic P-wave in unfavorable geologic structures is a complex process which contains reflex, transmission, and refraction. In this paper, the discontinuous deformation analysis (DDA) method is introduced, and the viscous boundary condition and stress wave input way are applied to enable the DDA method to simulate the seismic dynamic response of underground rock mass. The improved DDA program is used to analyze the propagation laws of elastic P-wave in joints and weak interlayer. The results show that (1) for a range of joint stiffness, the difference of transmission coefficient of joint surface between DDA calculation and theoretical solution is small; (2) the lower the joint stiffness, the more obvious the attenuation of the elastic P-wave is; (3) varied joint spacing and joint number play different roles on the attenuation of elastic P-wave; (4) the DDA calculation result is consistent with the theoretical solution for the transmission coefficient of interface; (5) the thickness and elastic modulus of weak interlayer will affect the interference effect of multiple transmission waves; (6) with the increasing weak interlayer thickness H (or a), the attenuation of the elastic P-wave is more obvious; and (7) the lower the wave impedance ratio of weak interlayer, the more obvious the attenuation of the elastic P-wave is.     

强度折减DDA方法及其在稳定分析中的应用

在自行研制的非连续变形分析（discontinuous deformation analysis,DDA）程序中实现了自动强度折减法以模拟边坡的稳定性和安全系数。通过锦屏高边坡和某重力坝深层滑动的计算分析以及与刚体极限平衡法和有限元方法的比较研究可以看出,强度折减DDA法可以较高精度求出岩石高边坡的安全系数,但求得的重力坝深层抗滑稳定安全系数要低于刚体极限平衡法的结果。基于强度折减技术的DDA与FEM的耦合方法,可以较准确地求出重力坝深层抗滑稳定安全系数。     

水电站地下洞室群分步开挖的非连续变形分析

利用大岗山水电站的地质勘探资料建立不同贯穿节理组合下地下洞室群的DDA模型,分析了洞室群开挖的破坏模式,并简单阐述其破坏机制。基于Monte-Carlo方法模拟节理的随机分布,结合大岗山水电站的开挖设计方案,对地下洞室群在分步开挖过程中的位移及应力演化进行研究。研究表明：洞室群开挖完成后,洞周块体的位移均指向洞内,尾调室上下游与主厂房上游的变形较大;洞周块体的大主应力均平行于洞室边界线,尾调室的应力释放比较明显,其边墙块体发生了滑移;洞室群在开挖过程中互相影响,其中主变室受主厂房和尾调室的影响比较明显。研究结果指出了洞室群开挖过程中需要关注和加固的具体部位,为工程支护措施的优化提供了依据。     

Modified predictor-corrector solution approach for efficient discontinuous deformation analysis of jointed rock masses

The high computational costs associated with the implicit formulation of discontinuous deformation analysis (DDA) have been one of the major obstacles for its implementation to engineering problems involving jointed rock masses with large numbers of blocks. In this paper, the Newmark-based predictor-corrector solution (NPC) approach was modified to improve the performance of the original DDA solution module in modeling discontinuous problems. The equation of motion for a discrete block system is first established with emphasis on the consideration of contact constraints. A family of modified Newmark-based predictor-corrector integration (MNPC) scheme is then proposed and implemented into a unified analysis framework. Comparisons are made between the proposed approach and the widely used constant acceleration (CA) integration approach and central difference (CD) approach, regarding the stability and numerical damping features for a single-degree-of-freedom model, where the implications of the proposed approach on open-close iteration are also discussed. The validity of the proposed approach is verified by several benchmarking examples, and it is then applied to two typical problems with different numbers of blocks. The results show that the original CA approach in DDA is efficient for the simulation of quasi-static deformation of jointed rock masses, while the proposed MNPC approach leads to improved computational efficiency for dynamic analysis of large-scale jointed rock masses. The MNPC approach therefore provides an additional option for efficient DDA of jointed rock masses.