共查询到20条相似文献，搜索用时 281 毫秒
1.
《Geomechanics and Geoengineering》2013,8(3):221-236
This paper presents a new point-to-face contact algorithm for contacts between two polyhedrons with planar boundaries. A new discrete numerical method called three-dimensional discontinuous deformation analysis (3-D DDA) is used and formulations of normal contact submatrices based on the proposed algorithm are derived. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. This approach does not need to use an iterative algorithm in each time step to obtain the contact plane, unlike the ‘Common-Plane’ method applied in the existing 3-D DDA. In the present 3-D DDA method, block contact constraints are enforced using the penalty method. This approach is quite simple, but may lead to inaccuracies that may be large for small values of the penalty number. The penalty method also creates block contact overlap, which violates the physical constraints of the problem. These limitations are overcome by using the augmented Lagrangian method that is used for normal contacts in this research. This point-to-face contact model has been programmed and some illustrative examples are provided to demonstrate the new contact rule between two blocks. A comparison between results obtained by using the augmented Lagrangian method and the penalty method is presented as well. 相似文献
2.
S. Amir Reza Beyabanaki Roozbeh Geraili Mikola Kianoosh Hatami 《Computers and Geotechnics》2008,35(3):346-356
This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC++. Examples are provided to demonstrate the new contact rule between two blocks. 相似文献
3.
In this paper, a new disk-based DDA formulation is presented. In the original disk-based DDA, disks are considered to be rigid and the penalty method is used to enforce disk contact constraints. In order to improve the accuracy of the disk-based DDA, new formulations of stiffness and force matrices for non-rigid disks using a new efficient contact model are presented in this paper. Blocks are considered deformable without need to do more computations for contact detection. In the proposed contact model, disk–disk and disk–boundary contacts are transformed into the form of point-to-line contacts and normal spring, shear spring and frictional force sub-matrices are derived by vector analysis. The penalty method is quite simple to implement, but has some major disadvantages. In the presented contact model, not only the simplicity of the penalty method is retained but also the limitations are overcome by using the augmented Lagrangian method. Moreover, unlike the contact model used in the original disk-based DDA, reference line can be obtained directly by using only coordinates of disk centers and their radii, and no more computations are needed. The validity and capability of the new disk-based DDA formulation are demonstrated by several illustrative examples. 相似文献
4.
S. Amir Reza Beyabanaki Ahmad Jafari S. Omid Reza Biabanaki M. Ronald Yeung 《Computers and Geotechnics》2009
This paper presents a new numerical model that can add a finite element mesh into each block of the three-dimensional discontinuous deformation analysis (3-D DDA), originally developed by Gen-hua Shi. The main objectives of this research are to enhance DDA block’s deformability. Formulations of stiffness and force matrices in 3-D DDA with conventional Trilinear (8-node) and Serendipity (20-node) hexahedral isoparametric finite elements meshed block system due to elastic stress, initial stress, point load, body force, displacement constraints, inertia force, normal and shear contact forces are derived in detail for program coding. The program code for the Trilinear and Serendipity hexahedron elements have been developed, and it has been applied to some examples to show the advantages achieved when finite element is associated with 3-D DDA to handle problems under large displacements and deformations. Results calculated for the same models by use of the original 3-D DDA are far from the theoretical solutions while the results of new numerical model are quite good in agreement with theoretical solutions; however, for the Trilinear elements, more number of elements are needed. 相似文献
5.
Accurate estimation of rockfall trajectory and motion behaviors is essential for rockfall risk assessment and the design and performance evaluation of preventive structures. Numerical simulation using discontinuous deformation analysis (DDA) is effective and helpful in rockfall analysis. Up to now, there have been many reports on application of two-dimensional (2-D) DDA programs. In this paper, the major advantages of rockfall analysis using 2-D and extensions to three-dimensional (3-D) analysis are presented. A practical 3-D DDA code is demonstrated to be capable of simulating free falling, rolling, sliding, and bouncing with high accuracy. Because rockfall trajectories and motion behaviors can be described as combinations of these four types, this demonstration indicates that the implemented code is capable of providing reliable rockfall analysis. Finally, specific tests are conducted to compare 2-D and 3-D DDA rockfall analysis in predicting trajectory and dynamic behavior. The results indicate that 3-D DDA simulations are more appropriate for rough tree-laden inclined slopes in providing detailed spatial distribution, whereas 2-D DDA simulations have better efficiency for slopes dominated by valleys and ravines. These results can help in selecting the appropriate DDA simulation for rockfall analysis. 相似文献
6.
传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。 相似文献
7.
G.H. Shi 《Geomechanics and Geoengineering》2014,9(2):80-96
Continuous computation and limit equilibrium computation are the two independent computations for practical rock engineering. For global stability analysis, limit equilibrium is still the fundamental method. For any numerical method, reaching limit equilibrium requires large displacements, discontinuous contacts, precise friction law, multistep computation and stabilised time-step dynamic computation. Therefore three convergences are unavoidable: convergence of equilibrium equations, convergence of open-close iterations for contacts and convergence of the contact forces of dynamic computations. This paper utilises mainly two dimensional discontinuous deformation analysis (DDA) and an available simple version of three dimensional DDA. The applications show DDA has the ability to reach the limit equilibrium of block systems. 相似文献
8.
In this paper, a lattice-type model to simulate the micro-mechanical behaviour of particulate/granular media is presented. In this numerical model, a particulate assembly is simulated as a lattice/truss. Nodes located at contacts between a particle and its neighbours are linked by bars to each other. Each particle is represented by a lattice within its microstructure and particle interact through load transfer at the nodes. Constraints are prescribed at the nodes to describe active, deactivated and reactivated contacts. When a particulate assembly develops into a mechanism (deformation with zero incremental load), further deformation is simulated through a framework that describes the kinematics of the particles (sliding, rolling and rotation of particles). This framework is formed by introducing nodes at the particle centroids and linking them with bars. Bars-linking particles with a non-sliding contact are assigned large stiffnesses relative to bars linking particles with a sliding contact. Numerical tests are conducted on two-dimensional assemblies of disks, arranged as very loose and very dense packing under simple shear loading conditions. The results concord with the results of numerical tests conducted using the discrete element method and with photoelastic experiments. Additionally, the model is applied to study the effects of initial imperfections caused by particles with low elastic modulus. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
9.
通过定义距离势函数,提出一种适用于空间任意凸多面体单元的三维距离势函数离散单元法。该方法采用归一化的计算方式,将势函数表征为接触体间的距离函数,并基于此建立接触力计算方程,明确了势函数的物理意义,使接触力计算更加合理,无需对各种可能的接触形式进行差异化处理。新方法克服了原有势函数物理意义不明、接触力计算受单元形式影响等重要缺陷,并突破了四面体单元的限制,可采用空间任意凸多面体单元。通过若干算例说明新方法的正确性和有效性。数值模拟结果表明新方法能够很好处理空间任意多面体单元复杂接触变换过程,准确处理复杂非连续介质的运动过程。 相似文献
10.
R. Resende L. N. Lamas J. V. Lemos R. Calçada 《Rock Mechanics and Rock Engineering》2010,43(6):741-761
The goal of this paper is to simulate the interaction of stress waves and rock fractures in a particle micromechanical model. Stress waves travelling in fractured rock masses are slowed down and attenuated by natural heterogeneities, voids, microcracks and, above all, by faults and fractures. Considerable laboratory and theoretical investigation have uncovered the major aspects of this phenomenon, but models that cover the core mechanisms of the wave propagation in rock masses are necessary to investigate aspects of wave–fracture interaction, which are not completely clear, and in the future simulate full-scale real problems. The micromechanical model is based on the particle discrete element model that reproduces rock through a densely packed non-structured assembly of 2D disks with point contacts. The model of a hard rock core is developed and an irregular rock joint is generated at mid-height. A new contact constitutive model is applied to the particles in the joint walls. Numerical static joint compression tests are performed and a typical hyperbolic stress–displacement curve is obtained. Conditions for good quality wave transmission through non-jointed unorganized particulate media are determined, hybrid static–dynamic boundary conditions are established and plane waves are emitted into the compressed joint. The transmitted and reflected waves are extracted and analysed. Joint dynamic stiffness calculated according to the hypotheses of the Displacement Discontinuity Theory shows to increase with the static joint compression until the joint is completely closed. Still in its early stages of application, this rock micromechanical model enables the joint behaviour under static and dynamic loading to be analysed in detail. Its advantages are the reproduction of the real mechanics of contact creation, evolution and destruction and the possibility of visualizing in detail the joint geometry changes, which is hard to accomplish in the laboratory. 相似文献
11.
A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis 总被引：15，自引：0，他引：15
Summary The key to three-dimensional discontinuous deformation analysis (3D DDA) is a rigorous contact theory that governs the interaction of many three-dimensional blocks. This theory must provide algorithms to judge contact types and locations and the appropriate state of each contact, which can be open, sliding or locked. This paper presents a point-to-face contact model, which forms a part of the contact theory, to be used in 3D DDA. Normal spring, shear spring and frictional force submatrices are derived by vector analysis and the penalty method. Also given are the open-close iteration criteria and operations performed for different changes in contact state. Sliding at a contact can occur in any direction parallel to the contact face, as opposed to one of two directions in two-dimensional DDA. This point-to-face contact model has been implemented into a 3D DDA computer program, and numerical results from several test cases demonstrate the validity of the model and the capability of the program. 相似文献
12.
David M. Doolin 《国际地质力学数值与分析法杂志》2005,29(12):1199-1207
Displacement boundary constraints in discontinuous deformation analysis (DDA) are applied using stiff penalty springs. A co‐ordinate‐free formulation for displacement boundary constraints is presented here for DDA, which unifies previous derivations for points of fixity, and for points constrained to induce or prohibit block motion in specified directions as a function of location or time. Examples for each type of constraint are used to illustrate the behaviour of the algorithm and provide a link with previous formulations for each case. The new, unified formulation has five benefits: (1) simple to express algorithmically; (2) easy to program and verify; (3) penalty values in different directions may be chosen to allow fixed points, lines, curves or planes; (4) formulation works for 2D and 3D; (5) displacement constraint may be a function of time or location or both. Feedback in the algorithm may induce internal resonance in homogeneously deformable discrete elements used in DDA, and resonance in block‐to‐block contact interactions. Consequently, high mass problems with insufficient damping may suffer from excessive ‘vibrational hammering’, inducing physically implausible behaviour such as elastic rebound. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
13.
A New Viscous Boundary Condition in the Two-Dimensional Discontinuous Deformation Analysis Method for Wave Propagation Problems 总被引：1，自引：1，他引：0
Viscous boundaries are widely used in numerical simulations of wave propagation problems in rock mechanics and rock engineering. By using such boundaries, reflected waves from artificial boundaries can be eliminated; therefore, an infinite domain can be modeled as a finite domain more effectively and with a much greater accuracy. Little progress has been made, thus far, with the implementation and verification of a viscous boundary in the numerical, discrete element, discontinuous deformation analysis (DDA) method. We present in this paper a new viscous boundary condition for DDA with a higher absorbing efficiency in comparison to previously published solutions. The theoretical derivation of the new viscous boundary condition for DDA is presented in detail, starting from first principles. The accuracy of the new boundary condition is verified using a series of numerical benchmark tests. We show that the new viscous boundary condition works well with both P waves as well as S waves. 相似文献
14.
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. 相似文献
15.
任意形状块体的接触模拟是发展三维不连续变形分析方法（3D-DDA）的关键问题。提出了三维离散块体单元的 边-边接触模型,并给出了边-边接触进入面的判定方法和嵌入准则。基于矢量分析和罚函数法,推导了接触子矩阵的表达式,算例分析结果表明了该模型的有效性。 相似文献
16.
A robust contact theory can be regarded as key to three dimensional discontinuous deformation analyses (3D-DDA). Not only must this theory provide an efficient algorithm to judge the type and location of contacts but also it must be able to present comprehensive formulations for every kind of contact (open, sliding and locked contact). There are six types of contact in three dimensional discontinuous deformation analyses (vertex-to-vertex, vertex-to-edge, vertex-to-face, edge-to-edge, edge-to-face and face-to-face) that can be converted to vertex-to-face and edge-to-edge contacts. This paper presents a new model of edge-to-edge contact to three dimensional discontinuous deformation analyses (3D-DDA). This new model considers both kinds of edge-to-edge contact (cross-over and parallel edge-to-edge contact) and presents a criterion for inter-penetration. Sub matrices of normal and shear spring and friction force are derived by geometrical analysis and penalty method. This new model is implemented in a 3D-DDA computer programme, and the example results demonstrate the validity of the model. 相似文献
17.
This paper describes a three-dimensional random network model to evaluate the thermal conductivity of particulate materials. The model is applied to numerical assemblies of poly-dispersed spheres generated using the discrete element method (DEM). The grain size distribution of Ottawa 20–30 sand is modeled using a logistic function in the DEM assemblies to closely reproduce the gradation of physical specimens. The packing density and inter-particle contact areas controlled by confining stress are explored as variables to underscore the effects of micro- and macro-scales on the effective thermal conductivity in particulate materials. It is assumed that skeletal structure of 3D granular system consists of the web of particle bodies interconnected by thermal resistor at contacts. The inter-particle contact condition (e.g., the degree of particle separation or overlap) and the particle radii determine the thermal conductance between adjacent particles. The Gauss–Seidel method allows evaluation of the evolution of temperature variation in the linear system. Laboratory measurements of thermal conductivity of Ottawa 20–30 sand corroborate the calculated results using the proposed network model. The model is extended to explore the evolution of thermal conduction depending on the nucleation habits of secondary solid phase as an anomalous material in the pore space. The proposed network model highlights that the coordination number, packing density and the inter-particle contact condition are integrated together to dominate the heat transfer characteristics in particulate materials, and allows fundamental understanding of particle-scale mechanism in macro-scale manifestation. 相似文献
18.
非连续变形分析（DDA）是一种隐式求解的动力学计算方法,且采用在块体界面加减刚硬弹簧的方式来满足块体界面无张拉和无嵌入的接触准则,其中时间步长和弹簧刚度两个物理量的取值直接影响DDA的计算结果。基于对DDA时间步和弹簧刚度在程序运行过程中的调整策略和块体接触的简化力学概念模型,研究了惯性力在DDA收敛求解中的作用过程。采用数值模拟试验对自由落体和斜面单滑块模型在3种力学状态下的相关力学问题进行了数值模拟研究,通过对自由落体运动的模拟,研究了时间步长单一因素对计算结果的影响规律,并初步确定了时间步长的合理取值区间。在此基础上,采用斜面单滑块模型,研究了时间步长和弹簧刚度对计算结果的共同影响,确定了不同时间步长条件下弹簧刚度的合理取值区间。研究成果表明,合适的时间步长和弹簧刚度的取值组合构成一个单连通参数取值域,当时间步和弹簧刚度的取值组合位于此“域”范围内时,DDA的计算结果是合理的。 相似文献
19.
Lu Zheng Guangqi Chen Yange Li Yingbin Zhang Kiyonobu Kasama 《Geomechanics and Geoengineering》2014,9(2):142-152
Rockfall is the most frequent major hazard in mountainous areas. For hazard assessment and further countermeasure design, realistic and accurate prediction of rockfall trajectory is an important requirement. Thus, a modeling method to represent both geometrical parameters of slope and falling rock mass is required. This study, suggests taking the advantages of discontinues deformation analysis (DDA) and geographical information system (GIS). In this study, after developing a three dimensional (3D) DDA program, firstly a special element named contact face element (CFE) was introduced into 3D DDA; secondly, effectively modeling tools with GIS support were developed. The implementation of CFE also improves the efficiency of both the contact searching and solution process. Then a simple impact model was devised to compare the 3D DDA implemented directly with a sliding model with theoretical analysis to verify the reliability of the modified 3D DDA program and investigate the parameter settings. Finally, simulations concerning rock shapes and multi-rocks were carried out to show the applicable functions and advantages of the newly developed rockfall analysis code. It has been shown that the newly developed 3D DDA program with GIS support is applicable and effective. 相似文献
20.
In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献