首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 671 毫秒
1.
马永政  蔡可键  郑宏 《岩土力学》2016,37(3):867-874
传统的非连续变形分析法(DDA)法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。  相似文献   

2.
林兴超  汪小刚  王玉杰  李旭  韩鑫晔 《岩土力学》2011,32(10):3065-3070
数值流形法是至少包含流形法(NMM)、有限元法(FEM)和非连续变形分析(DDA)的数值方法体系。将数值流形法中物理单元与数学单元完全重合,去掉接触理论,流形元能够回归到有限元,将通过简单的板压缩数值试验验证这一点。在以前的数值流形法法中,质量守恒问题一直被忽视,物理单元的质量会随着单元体积改变,计算结果存在一定的误差。通过改变计算过程中单元密度实现计算过程中的"质量守恒",完善了现有数值流形法的理论基础。  相似文献   

3.
Updating the block configuration on the basis of additive decomposition and its linearized expression of the displacement increment leads to the low calculation accuracy of the original discontinuous deformation analysis (DDA) and false volume expansion. In this study, the displacement expressions of a small deformation, a large rotation, and the corresponding velocity and acceleration terms on the basis of the initial configuration are presented using multiplicative decomposition. With the use of the principle of virtual work, the stiffness matrix, mass matrix, and force vector of blocks are obtained. Compared with the original DDA, each of the block deformation parameters has obvious physical meaning as a parameter of mechanics, which can be obtained by adding the incremental deformation components of each time step directly without co-ordinate transformation. Moreover, the proposed modification automatically considers the block deformation produced by centrifugal and Coriolis forces. The analysis of some typical numerical examples have verified the accuracy of the strain and stress calculated by the proposed method, and the current configuration is updated by the total displacements, which completely overcomes the false volume expansion and provides reasonable linear strains.  相似文献   

4.
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.  相似文献   

5.
李小凯  郑宏 《岩土力学》2014,35(6):1787-1794
非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。  相似文献   

6.
A modified three‐dimensional discontinuous deformation analysis (3D‐DDA) method is derived using four‐noded tetrahedral elements to improve the accuracy of current 3D‐DDA algorithm in practical applications. The analysis program for the modified 3D‐DDA method is developed in a C++ environment and its accuracy is illustrated through comparisons with several analytical solutions that are available for selected problems. The predicted solutions for these problems using the modified 3D‐DDA approach all show satisfactory agreement with the corresponding analytical results. Results presented in this paper demonstrate that the modified 3D‐DDA method with discontinuous modeling capabilities offers a useful computational tool to determine stresses and deformations in practical problems involving fissured elastic media with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
赵强  焦玉勇  张秀丽  谢壁婷  王龙  黄刚海 《岩土力学》2019,40(11):4515-4522
非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。  相似文献   

8.
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.  相似文献   

9.
Summary  Numerical Manifold Method (NMM) is an extension of Discontinuous Deformation Analysis which contains and combines FEM and joint/block oriented DDA in a unified form. At present, application of NMM is mainly limited to two-dimensional analysis which is not realistic for many underground excavation works. In the present paper, three-dimensional NMM based on tetrahedron and hexahedron elements is proposed. The corresponding three-dimensional cover contact detection algorithms which have some fundamental properties different from the corresponding two-dimensional analysis are also developed. Based on the present formulations, engineering problems with arbitrarily oriented discontinuities can be modeled and realistic modeling of construction works becomes possible. Author’s address: Y. M. Cheng, Department of Civil and Structural Engineering, Hong Kong Polytechnic University, Room TV707, Hung Hom Kowloon, Hong Kong, P.R. China  相似文献   

10.
As a hybrid method, the nodal‐based discontinuous deformation analysis (NDDA) greatly improves the stress accuracy within each DDA block by coupling a well‐defined finite element mesh inside the DDA block; at the same time, the NDDA inherits the unique block kinematics of the standard DDA method. Each finite element mesh line inside the DDA block is treated as a potential crack, which enables the transformation of the block material from continuum to discontinuum through the tensile and shear fracturing mechanism. This paper introduces a double minimization procedure into the NDDA method to further improve the accuracy of the stresses evaluated at the finite element mesh lines and thus to obtain a more realistic fracture model. Three numerical examples are employed to demonstrate the improved stress accuracy by the implemented double minimization procedure and the accuracy and capability of the enhanced NDDA method in capturing brittle fracturing process. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

11.
Over the last decade, researchers in the discontinuous deformation analysis (DDA) community have dedicated a great deal of effort to document the accuracy of the method by performing validation studies. This paper contains a summary of more than 100 published and unpublished validation studies which comprise the body of DDA validation information to which the authors have access. The studies are grouped into three general categories: (a) validation with respect to analytical solutions, (b) validation with respect to results of other numerical techniques, and (c) validation with respect to laboratory and field data. Three general techniques for validation are described: qualitative assessment visually examining runtime behaviour of simulations, semi‐quantitative assessment comparing numerical results of simulations, and quantitative where numerical simulation results are evaluated in detail with respect to similar analytical, laboratory or field results. We find that for many of the problems addressed by the papers in this review, DDA performs more than adequately for engineering analysis. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
刘君  黄盛铨  孔宪京 《岩土力学》2006,27(Z2):286-290
在自行研制的非连续变形分析(discontinuous deformation analysis,DDA)程序中实现了自动强度折减法以模拟边坡的稳定性和安全系数。通过锦屏高边坡和某重力坝深层滑动的计算分析以及与刚体极限平衡法和有限元方法的比较研究可以看出,强度折减DDA法可以较高精度求出岩石高边坡的安全系数,但求得的重力坝深层抗滑稳定安全系数要低于刚体极限平衡法的结果。基于强度折减技术的DDA与FEM的耦合方法,可以较准确地求出重力坝深层抗滑稳定安全系数。  相似文献   

13.
Nodal-based three-dimensional discontinuous deformation analysis (3-D DDA)   总被引:2,自引:0,他引:2  
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.  相似文献   

14.
刘永茜  杨军 《岩土力学》2011,32(8):2544-2548
非连续变形分析(DDA)是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径  相似文献   

15.
张友良  刘登学  刘高敏 《岩土力学》2016,37(8):2404-2410
在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。  相似文献   

16.
The coupled numerical manifold method (NMM) and discontinuous deformation analysis (DDA) are enhanced to simulate deformations of continuous soil and discontinuous masonry structures. An elasto-plastic NMM-DDA is formulated that incorporates elasto-plastic constitutive laws into incremental forms of the equation of motion. A node-based uniform strain element is applied to avoid volumetric locking, which often occurs in conventional NMM-DDA. The proposed method is applied to three fundamental boundary value problems: a beam bending problem, a bearing capacity problem of a footing, and a bearing capacity problem of a masonry structure. The method is verified through comparisons with conventional solutions.  相似文献   

17.
包承纲 《岩土力学》2011,32(Z2):1-9
介绍了在非饱和土、数值分析、测试技术和土工合成材料工程应用方面的成就。关于非饱和土方面,指出了加强非饱和土基本试验研究工作的紧迫性;关于数值分析方面,认为流形元法为连续和非连续变形的力学分析提供了一个统一的方法,应加强在岩体力学和土力学中的开发与应用;关于岩土力学的测试技术,需密切注意近年高新技术在岩土领域的引进和推广;对于土工合成材料,应加强机制的研究,注意在应用中掌握有关原理,最后对岩土工程研究中当前存在的问题和今后发展作了简单的评述  相似文献   

18.
江巍  郑宏 《岩土力学》2007,28(12):2603-2606
非连续变形分析方法(DDA) 自创立以来,已在分析不连续体的大变形、大位移的问题中显示了其优越性。但DDA方法需要人为设置一些控制参数,通过一个简单的滑块模型,分别采用不同的时步大小、不同的接触弹簧刚度、并比较相对误差,分析了人为设定的参数对DDA计算结果的影响,指出仅在合理的取值范围内,大时步和大的接触弹簧刚度才能产生较好的计算结果;同时对造成这种影响的原因做了一定的探讨,并建议了参数取值的上下限原则。  相似文献   

19.
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.  相似文献   

20.
Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号