共查询到20条相似文献，搜索用时 962 毫秒
1.
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. 相似文献
2.
非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。 相似文献
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.
非连续变形分析（DDA）是一种隐式求解的动力学计算方法,且采用在块体界面加减刚硬弹簧的方式来满足块体界面无张拉和无嵌入的接触准则,其中时间步长和弹簧刚度两个物理量的取值直接影响DDA的计算结果。基于对DDA时间步和弹簧刚度在程序运行过程中的调整策略和块体接触的简化力学概念模型,研究了惯性力在DDA收敛求解中的作用过程。采用数值模拟试验对自由落体和斜面单滑块模型在3种力学状态下的相关力学问题进行了数值模拟研究,通过对自由落体运动的模拟,研究了时间步长单一因素对计算结果的影响规律,并初步确定了时间步长的合理取值区间。在此基础上,采用斜面单滑块模型,研究了时间步长和弹簧刚度对计算结果的共同影响,确定了不同时间步长条件下弹簧刚度的合理取值区间。研究成果表明,合适的时间步长和弹簧刚度的取值组合构成一个单连通参数取值域,当时间步和弹簧刚度的取值组合位于此“域”范围内时,DDA的计算结果是合理的。 相似文献
6.
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. 相似文献
7.
基于DDA方法一种流-固耦合模型的建立及裂隙体渗流场分析和应用 总被引：1，自引：0，他引：1
以非连续变形分析方法（DDA）为基础并采用稳态流体计算方法将二者结合进行裂隙岩体流-固耦合分析。利用DDA方法生成裂隙岩体模型,在此基础上采用矩阵搜索等方法形成新的裂隙水通网络模型。采用稳态迭代算法和立方定律求得裂隙水压力,并把裂隙水压力作为线载荷施加到块体边界,在DDA算法中每个迭代步完成后更新裂隙开度和水压值,与DDA算法结合研究裂隙水与块体之间相互作用关系。利用以上裂隙岩体流-固耦合计算方法研究了某水封油库开挖和运行过程洞室围岩流量和密封性,为该工程预测水封效果提供了有益的主要依据,也是国内首次采用DDA方法做大型工程的流-固耦合模型分析。 相似文献
8.
非连续变形分析(DDA)方法对大规模工程问题的数值模拟耗时太长,其中线性方程组求解耗时可占总计算时间的70%以上,因此,高效的线性方程组解法是重要研究课题。首先,阐述了适用于DDA方法的基于块的行压缩法和基于试验-误差迭代格式的非0位置记录;然后,针对DDA的子矩阵技术,将块雅可比迭代法 (BJ)、预处理的块共轭梯度法 (PCG,包括Jacobi-PCG、SSOR-PCG) 引入DDA方法,重点研究了线性方程组求解过程中的关键运算;最后,通过两个洞室开挖算例,分析了各线性方程组求解算法在DDA中的计算效率。研究表明：与迭代法相比,直解法无法满足大规模工程计算需要;BJ迭代法与块超松弛迭代法(BSOR)的效率差别不大,但明显不如PCG迭代法。因此,建议采用PCG迭代法求解DDA线性方程组,特别是SSOR-PCG值得推广;如果开展并行计算研究,Jacobi-PCG是较好的选择,当刚度矩阵惯性优势明显时,BJ迭代法同样有效。 相似文献
9.
The discontinuous deformation analysis (DDA) is a discontinuum‐based method, which employs a penalty method to represent the contact between blocks. The penalty method is easy to be implemented in the program, but the contact constraint is only approximately satisfied. Penetrations between contacting blocks are unavoidable even if the penalty value is very large. To improve the contact precision in the DDA, an augmented Lagrangian method is introduced, which can make use of advantages of both the Lagrangian multiplier method and the penalty method. This paper provides a detailed implementation of the augmented Lagrangian method in the DDA program and compares it with the standard DDA on the computational efficiency and contact precision. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
Conventional discontinuous deformation analysis (DDA) results in a change of block volume, which is known as free expansion, during rotation calculations because of the use of a linear displacement function to simulate the behavior of a block with a rigid body and elastic behaviors. This study demonstrates that the linear displacement function also generates unsolved elastic distortion, especially when the block undergoes large rotation in each calculation step. The distortion disturbs the contact judgment in the open–close iteration and update calculations of vertex coordinates, stresses, velocities, etc. at the end of each calculation step. A new procedure follows the flow chart of the original DDA, but it adopts additional codes for the coordinate-transformation calculations in vertex coordinate, stress, and velocity updates. When the vertex coordinates are updated, vertex displacements caused by strains are calculated before involving the block-rotation term in the displacement function to mitigate the elastic distortion. In addition, new codes compile formulas to transform stresses and velocities with block rotation. The new DDA ensures the correctness of rotating elastic calculations to solve practical falling rock problems with a large rotational angle in each calculation step. 相似文献
11.
传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。 相似文献
12.
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. 相似文献
13.
Disk clusters are developed to represent the shape of granular materials more precisely (compared to circular particles) and to minimise excessive rolling. Investigating the behaviour of dynamic disk-based discontinuous deformation analysis (DDA) with disk clusters is very important to evaluate the applicability of disk-based DDA to dynamic problems in geomechanics. In this paper, the accuracy of disk-based DDA under dynamic conditions is studied by a comprehensive sensitivity analysis. The results obtained by disk-based DDA are compared with the analytical solutions of a disk cluster on an incline subjected to gravitational force only, and three different accelerations of increasing complexity with sinusoidal input functions as well as gravitational load. In this research, the effects of time step size and interface friction angles on the results are studied. Overall, most of the error for both velocity and displacement occurs at the beginning of the solution. With increasing friction angle, the initial perturbation of the solution increases in the case of sliding under gravitational force only, and decreases in the case of sliding under dynamic loads. This study shows that disk-based DDA predicts accurately the velocities and displacements derived with respect to the frictional resistance offered by the inclines. 相似文献
14.
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. 相似文献
15.
《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. 相似文献
16.
非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径 相似文献
17.
非连续变形分析方法（DDA）是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。 相似文献
18.
In modelling particulate media, such as soils, using 3-D DDA, there are four types of contacts: sphere to sphere, sphere to boundary face, sphere to boundary edge, and sphere to boundary corner contacts. The first two were studied by the authors in a previous work (Beyabanaki and Bagtzoglou 2012). In this paper, we present a new contact model for sphere-boundary edge and sphere-boundary corner contacts in sphere-based 3-D DDA. The model includes a new algorithm to search for contacts, detect the contact types and calculate contact points. Moreover, formulas for contact sub-matrices are derived. The proposed contact model has been implemented into a sphere-based 3-D DDA program and three test cases are studied in order to verify the workability of the new contact model. The numerical results obtained demonstrate the capability of the model to deal with sphere-boundary interaction in particulate media. 相似文献
19.
Lanru Jing 《Engineering Geology》1998,49(3-4):371-381
Continuum-based numerical methods have played a leading role in the numerical solution of problems in rock mechanics and engineering geology. However, for fractured rocks, a continuum assumption often leads to difficult parameters to define and over-simplified geometry to be realistic. In such case, discrete representations of fractures and individual blocks must be adopted. In this paper, a newly emerged member in the family of discrete element methods (DEM), the discontinuous deformation analysis (DDA), is presented, including its variational principle, governing equations, solution techniques and contact representation and detection algorithms. Its relative advantages and shortcomings are compared with the explicit distinct element method and the finite element method. An example of the analysis of tunnel stability is provided to demonstrate the capability of this new method. 相似文献
20.
Boundary setting method for the seismic dynamic response analysis of engineering rock mass structures using the discontinuous deformation analysis method 
下载免费PDF全文
下载免费PDF全文 Large deformations and discontinuous problems can be calculated using the discontinuous deformation analysis (DDA) method by solving time steps, and this method is suitable for simulating the seismic dynamic response of engineering rock mass structures. However, the boundary setting must be carefully analyzed. In this paper, four boundary settings for the DDA method are investigated. First, the contributions to the DDA equations for nonreflecting boundaries (including the viscous boundary and the viscoelastic boundary) are deduced based on the Newmark method. Second, a free‐field boundary is introduced in the DDA method with boundary grid generation and coupling calculation algorithms to accurately simulate external source wave motion, such as earthquakes. Third, seismic input boundary treatments are intensively examined, and the force input method is introduced based on nonreflecting boundaries. Finally, the static‐dynamic unified boundary is implemented to ensure consistent boundary transformation. The boundary setting method in the DDA method is discussed, and the suggested treatments are used to analyze the seismic dynamic response of underground caverns. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献