论文部分内容阅读
系统深入地研究了三维离散单元法的基本原理和数值解算方法，在改进和推广了现有算法的基础上，自行编制了相应的Ｃ语言程序。采用投影方法求解接触力的作用点和作用面积，提高了接触计算的精度和速度。用矢量方法推导了三维动态松弛离散元法的理论公式，为程序编写提供了方便；同时，对阻尼和时步采取自动调节技术，较好地解决了块体间的侵入问题。首次研究了三维静态松弛的离散单元法，给出了详细的理论推导，并与动态松弛法相耦合，克服了两者的缺点；通过算例的比较，发现静态方法比动态方法具有更高的效率。提出了三维离散元法中地下水和锚杆的模拟方法，认为地下水可能诱发滑坡，锚杆通过悬吊和销接作用对块体产生锚固力。探讨了三维离散元法的数据结构，提出了地下水、锚杆、边界力等数据的组织方法。论述了三维离散元法的程序编制方法和特点，分别开发了动态和静态松弛法的Ｃ语言程序。对海冰撞击钻井平台、杨家槽滑坡和长壁工作面覆岩运动规律进行了计算，从计算后的块体系统变形形态来看，模拟结果与实际情况大致吻合。
The system thoroughly studies the basic principle and numerical solution method of the three-dimensional discrete element method. On the basis of improving and popularizing the existing algorithms, the corresponding C language program is compiled by itself. The projection method is used to solve the point of contact and the area of contact force, which improves the accuracy and speed of contact calculation. The theoretical formulas of the three-dimensional dynamic relaxation discrete element method are deduced by the vector method, which provides convenience for program preparation. At the same time, automatic adjustment techniques for damping and time steps are adopted to solve the problem of intrusion between blocks. For the first time, the discrete element method for three-dimensional static relaxation is studied. The detailed theoretical derivation is given and coupled with the dynamic relaxation method to overcome the shortcomings of the two. By comparison of examples, it is found that the static method has higher efficiency than the dynamic method. . The simulation method of groundwater and anchor rods in the 3D discrete element method is proposed. It is believed that groundwater may induce landslides, and the anchor rods generate anchorage forces through the suspension and pinning effects on the blocks. The data structure of three-dimensional discrete element method is discussed, and the methods for organizing data such as groundwater, anchor rods, and boundary forces are proposed. The programming method and characteristics of three-dimensional discrete element method are discussed. The C language programs of dynamic and static relaxation methods are respectively developed. The law of overlying strata of the sea ice impact drilling platform, Yangjiacao landslide and longwall face was calculated. From the calculated deformation of the block system, the simulation results are generally consistent with the actual conditions.