DDA程序二次开发及工程应用

在石根华先生DDA程序的基础上,中国水利水电科学研究院作了二次程序开发,实现了Windows用户界面,丰富了前后处理的功能,计算中也可以模拟边坡的开挖过程。介绍了程序的新功能、新特点以及程序在向家坝重力坝抗滑稳定性研究中的应用。     

陡倾断层上下盘开挖引起地表变形的数值模拟分析

考虑地下开挖区与断层相对位置关系,采用简化数值模型对陡倾断层上下盘开挖引起的地表变形特征进行了数值模拟。结果表明,当断层位于地下开挖引起地表变形的压缩区且开挖区位于下盘时,地表出现不连续变形的可能性小;当断层位于地下开挖引起地表变形的压缩区且开挖区位于上盘时,随下向开挖进行地表断层处裂缝有减小趋势;而当断层位于地下开挖引起的地表变形拉张区时,断层出露处水平拉张明显,表现为拉张裂缝,随下向开挖深度增加、规模增大,不论开挖区位于上盘还是下盘,极可能导致地表出现正断层式的错动．     

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

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

基于DDA方法的膨胀土边坡稳定性研究

非连续变形分析方法（DDA）是近年来发展起来的数值计算方法,其以新颖的思想、先进的数值技术,得到了学术界的认可和广泛关注。针对宁明高速公路膨胀土路堑边坡破坏的特点,建立地质模型和计算模型,编制了二维DDA程序,对DDA理论模型进行了推广和发展,将非饱和土的本构模型引入程序。对非连续变形分析方法中的控制参数设置进行了探讨,应用DDA程序模拟了膨胀土边坡破坏特征。程序计算结果与工程实际情况比较吻合,验证了程序的合理性和有效性,分析成果可为设计、施工提供参考依据。     

水布垭马崖高边坡岩体地下开挖三维数值模拟研究

清江水布垭水电站是清江上游拟建的一个特大型水电工程,马崖高边坡位于水布垭峡谷出口右侧,总体地质结构具上硬下软特征,电站导流洞及尾水洞将从马崖高陡边坡坡脚穿出,地下开挖主要在软岩层中。应用FINAL及ROCKY程序对水布垭马崖高陡边坡岩体地下开挖进行数值模拟,研究了边坡岩体内进行地下洞室开挖后边坡力学状态的改变及洞室岩体（含岩柱）的稳定性,基于削坡是水布垭马崖高陡边坡的重要整治措施之一,研究中进行了各种分步开挖工况下的削坡有效性定量研究,重要工况还同时进行了VI度地震作用下的动力作用分析。     

冰-结构动力相互作用过程的DDA模拟

应用非连续变形分析(DDA)方法, 模拟了海冰与直立结构的动力相互作用过程的模拟.该方法的优点是能将冰对结构的碰撞及冰的破碎过程可视化, 同时发生碰撞后结构在每一时间步的响应都可计算出来.数值算例表明, DDA方法进行海冰破碎过程的模拟, 比有限元或边界元等方法更为合理和接近实际.     

彭水工程地下厂房开挖及支护的数值模拟

根据彭水工程地下厂房开挖设计，采用弹塑性有限元法进行中，尾部方案的数值模拟计算，对开挖洞群的应力位移状况，各洞之间的相互影响及稳定性等问题进行了对比分析与评价，并提出了该地地下厂房支护加固设计的主要原则。     

连拱隧道围岩变形和破坏形态的DDA分析

连拱隧道围岩一般为由节理结构面相互切割的非连续岩体,利用非连续性分析方法研究这类围岩的变形和破坏形态可以更好地反映工程实际情况。采用非连续变形分析方法 DDA,对金鸡山连拱隧道围岩的变形和破坏过程进行了模拟,将整个变形破坏过程划分为了3个阶段：中墙上方岩体变形、地表下沉、滑移面产生阶段、隧道左右洞两侧滑移带（或滑移面）的形成阶段以及中墙顶部块体失稳、隧道上方岩体快速塌落阶段。研究了金鸡山隧道浅埋围岩和深埋围岩的变形破坏特征,其中边墙部位、靠近中墙的内侧拱顶或拱肩部位、外侧拱肩部位会首先受到变形破坏。     

基于PSO-MLSSVR的土体参数反演方法在深基坑工程中的应用

为了确保基坑工程安全,常常会采用数值模拟的方法预测支护结构的位移,其中岩土体力学参数的选取对于结果的影响最大.本文使用了一种粒子群(PSO)算法结合多输出最小二乘支持向量回归机(MLSSVR)的基坑土体参数位移反分析法,以深圳某深基坑的支护桩顶水平位移监测数据为依据,基于正交设计生成具有代表性的土体参数组合,通过有限元...     

非开挖HDD地表变形数值模拟及DCRP识别技术研究

采用非开挖水平定向钻进技术近地表施工时,由于地层及钻井液压力等原因易产生地表变形。结合美国梅萨市一处通信管道施工实例,对0-0.9 MPa钻井液压力作用下的地表变形程度进行了数值分析,得到地表最大隆起量为26 mm。采用DCRP技术及水准仪方法对该变形区域钻井液压力作用前后的地表进行位移检测,得到地表最大隆起量分别为28 mm和30 mm,表明数值模拟分析结果同现场检测结果比较吻合。同时,在现场所布置的26个检测点中,两种方法检测到的隆起量最大相差2.7 mm,表明其精度比较接近。此外,由于DCRP技术属于非接触式三维检测,受施工现场条件及人员操作限制少,故具有较为理想的应用前景。     

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

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

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

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

GPU-based parallel computation for discontinuous deformation analysis (DDA) method and its application to modelling earthquake-induced landslide

Graphic Processing Unit (GPU), as a computing device, has upgraded from single-subject graphical processors to multi-core processors with tremendous computational horsepower. This paper proposes to accelerate the DDA using parallel Jacobi Preconditioned Conjugate Gradient (JPCG) technique on GPUs. Based on the results of two numerical examples, the calculation accuracies of the DDA with serial and parallel solvers are validated, and we found that the DDA with parallel solvers exhibits a much higher execution efficiency. The movement process of Daguangbao landslide triggered by the Wenchuan earthquake is replicated and the modeled deposit pattern coincides well with the actual topography after earthquake.     

Applying disk-based discontinuous deformation analysis (DDA) to simulate Donghekou landslide triggered by the Wenchuan earthquake

In this paper, for the first time, disk-based discontinuous deformation analysis (DDA) is applied to simulate a real landslide triggered by an earthquake. For this purpose, the kinematic behaviour of the Donghekou landslide triggered by the Wenchuan earthquake is simulated and the results obtained using disk-based DDA are compared with actual data. The comparisons show that there is a good agreement between the results obtained using disk-based DDA and observed data. The simulation results provided an understanding of the failure behaviour and temporal evolution of the landslide. This study shows that disk-based DDA is a practical numerical tool that can be used to simulate the post-failure behaviour of landslides triggered by an earthquake.     

节理岩质边坡变形的DDA模拟

对原非连续变形分析方法（DDA方法）中边界约束方法进行了推广，将其应用于节理岩质边坡变形规律的数值模拟，并利用模型实验结果验证了推广方法的正确性。     

三维激光扫描技术及其在岩土工程中的应用

地面三维激光扫描技术是一项较新的数据收集与分析技术,它可以快速扫描目标,获得高密度,高精度的目标空间坐标,配合后处理软件,进行各种测绘、分析、仿真工作。介绍了三维激光扫描技术的原理,对目前的主要激光扫描硬件进行了比较,最后介绍了三维激光扫描技术在岩土工程应用的主要两个方面：边坡变形监测与岩石裂隙面勘察。相比常规监测和地质调查技术,三维激光扫描技术具有高精度、速度快、节省人工的优势,具有较大的应用潜力。     

自适应非线性有限元的多重网格法求解及应用

描述了非线性自适应有限元的多重网格法求解过程。 通过算例说明了方法的有效性以及自适应多重网格有限元用于开挖施工过程模拟分析的优越性。     

物质点强度折减法及其在边坡中的应用

物质点法适用于模拟连续介质大变形,如边坡失稳全过程。在物质点法中应用强度折减法,用于边坡稳定性评价。与极限平衡法相比,二者安全系数计算值、滑动面位置结果基本一致;与有限元强度折减法相比,物质点法失稳评价标准的物理意义明确。利用物质点法大变形计算优势,评价边坡失稳后的破坏后果,通过算例说明其评价不同安全系数下的滑坡堆积形态及滑移距离的能力,尤其是评价滑坡对临近建筑物的影响程度的能力。物质点强度折减法可用于边坡稳定性评价及边坡破坏后果评价。