三维离散元改进算法及其在三峡永久船闸高边坡工程中的应用

本文根据面－面接触模型，写出了三维离散元的物理方程及运动方程，对于接触刚度、阻尼系数、时步、计算控制及块体的滑动失稳判据等作了讨论。在NURBM－3D 长方体单元基础上改进了其划分单元的算法，依据岩体中三组真实的节理面划分平行六面体单元，实现了网格的自动剖分，从而可以将单元视为真实的岩块，更客观地体现了节理面对岩体稳定性的影响。运用改进后的三维离散元方法，编制了计算程序，对三峡工程永久船闸高连坡开挖进行了模拟。计算过程包括初始地应力场模拟、岩体开挖过程两个阶段。模拟计算首先给出了为水平、垂直节理面所切割的长方体单元情况下的应力场分析，考证了程序的可行性。将已有的三峡永久船闸实测节理简化为三组节理，完成了初始地应力场的模拟，给出了分四步开挖后的岩体位移场和应力场，计算出岩体开挖过程中产生滑动块体的位置，其结果可以定性地表明：沿节理面的块体滑动是边坡失稳的主要模式。编制的后处理程序提供了计算数据库与绘图软件Origin的接口，可以画出任意剖面或节理面上的应力等值线、位移矢量图及滑动块体分布图，由此可以直观判断岩体不稳定区域和高应力区范围。编制的计算程序在PC机上计算了10万块体单元，用于模拟大型工程，其计算结果与有限元及类似的连续介质模型相比更为合理，为用离散元法进一步研究可变形块体及岩体的蠕变提供了一种切实可行的工具。将离散元法应用于大型工程三维计算是可行的，但在本构模型建立、岩体破坏参数选择等方面均有等于更深入的研究，其结果也有待于更多工程的验证。

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.