共查询到20条相似文献，搜索用时 15 毫秒
1.
原有非连续变形分析（DDA）采用一阶近似后的位移增量表达式更新块体构形,推导相关子矩阵,且对不同时步计算出的应变增量直接叠加,当模拟的块体发生大转动时往往会产生较大误差。为考虑块体转动与变形的耦合作用,引入先变形、后转动的块体位移增量表达式。重新推导了惯性力子矩阵,将块体转动时的离心力与科氏力加到荷载矩阵中。计算时对应变分量及其相关变量进行坐标变换与修正,并采用新引入的位移增量表达式计算块体顶点位移,进行后接触修正与更新块体构形。数值算例表明,改进后的程序能够消除转动带来的误差,自动考虑了块体转动时离心力和科氏力引起的变形,应变计算精度更高。改进方法克服了块体体积自由膨胀、应变场畸变等问题,给出了合理的块体应变。 相似文献
2.
Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
3.
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. 相似文献
4.
非连续变形分析方法（DDA）是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。 相似文献
5.
An extensive examination of the discontinuous deformation analysis (DDA) in block dynamic sliding modeling is carried out in this paper. Theoretical solutions for a single block sliding on an arbitrarily inclined plane by applying the horizontal/vertical seismic loadings to the sliding block as acceleration time histories or to the base as constraint displacement time histories are derived. As compared with the theoretical solutions, for a single block sliding, the DDA predicts the sliding displacements and block interaction forces accurately under various base incline angles and friction angles under both the harmonic loadings and a real seismic loading. The vertical seismic component may influence the block sliding displacements to different extent, and the DDA can capture these phenomena successfully and give accurate results. For the calculation of the single block relative sliding, both the theoretical and the DDA solutions indicate that applying the seismic accelerations as constraint displacement time histories (derived by integrating the seismic accelerations twice) to the base is equivalent to applying the seismic accelerations as volume forces to the sliding block in the opposite directions. The DDA modeling also demonstrates that this conclusion still stands for the case of multi‐block sliding. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。 相似文献
7.
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. 相似文献
8.
This study explores a new form of Discontinuous Deformation Analysis (DDA) method, which uses mesh free displacement functions instead of linear polynomial ones adopted in the original DDA, hence it can effectively describe complex block displacements and deformations. Moreover, the capability of describing a block’s nonlinear mechanical behaviors, i.e., plasticity and fracture, are developed under classical fracture and increment plasticity mechanics theories. With consideration of computation efficiency and convenience, the Sibson natural neighbor interpolation technique for block plasticity analysis and the enriched Moving Least-Squares (MLS) approximation for block fracture analysis are employed, respectively. Numerical results show the applicability of the proposed mesh free DDA method. 相似文献
9.
非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。 相似文献
10.
《Geomechanics and Geoengineering》2013,8(2):147-153
With high-order numerical manifold method (NMM) or high-order discontinuous deformation analysis method (DDA), computational accuracy of structure deformation can be improved greatly. However, poor accuracy is obtained and even computation is not convergent while treating large deformation problems, due to inaccurate or incorrect high-order initial stress formulae. Based on 2D triangular mathematical meshes and polynomial cover functions in high-order NMM, exact formulae for high-order initial stresses are deduced to depict configuration change of structures under large deformations. The formulae are expressed in polynomial forms so as to be used in simplex integrations. The approach is also extended to high-order DDA. Comparing with analytical solutions, accurate results for large deformation of a cantilever beam prove the validity of the proposed formulae. 相似文献
11.
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. 相似文献
12.
定量研究了非连续变形分析(DDA)方法中的黏性阻尼与数值阻尼。首先,基于Newmark直接积分法,推导了块体系统的运动方程。其次,通过动力学中的黏性阻尼理论建立了DDA中动力系数、时间步长与黏性阻尼比的关系式,探讨了DDA中的常加速度积分方案的数值阻尼分区及阻尼比计算方法,进而得到两种阻尼共同作用时的阻尼比表达式,并分析了频域内阻尼比的分布情况。最后,以谐振激励下的块体振动为例,通过对比不同阻尼作用下块体位移的DDA计算值与理论解,验证了本文提出的阻尼比计算公式的正确性。研究表明：黏性阻尼对低频的衰减作用明显,数值阻尼则可以很快地消除高频干扰,而二者共同作用下可降低阻尼的频率相关性。该研究成果为DDA的振动、波动等动力计算的阻尼取值提供了理论依据。 相似文献
13.
《Geomechanics and Geoengineering》2013,8(4):293-302
This paper presents a comparative study of two methods, Sarma's method and the discontinuous deformation analysis (DDA), for rock slope stability analysis. The comparison concerns the stability analysis of two classic rock slopes. The study shows that the DDA, which accounts for the block kinematics, provides a very different factor of safety as compared with Sarma's method. More realistic reaction forces around each rock block can be obtained by the DDA, including the thrust forces between rock blocks and the forces between the base and the blocks. The DDA's result shows two possible directions for the relative movement between two contiguous blocks at the initiation of slope failure. It also indicates that the limit equilibrium condition may not occur along the interfaces of rock blocks at the initiation of slope failure. The determination of realistic interaction forces around each block will be very important in rock slope stability analysis if nonlinear failure criteria are considered. 相似文献
14.
非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径 相似文献
15.
In the original formulation of DDA by Shi, a linear displacement function term is used. This has the limitations of uncontrolled block and stress distortion due to rigid body rotation. In the present paper, the authors propose a new iterative method which can avoid the distortion due to the rotation even when the rotation or number of time‐step is large. Furthermore, the authors propose a simple internal discretization scheme which is applicable for both concave and convex polygon which is particularly important for a large block. The stress and strain distribution with a large block can be obtained with ease under this scheme. The numerical examples as shown have demonstrated the advantages of the present proposal in DDA analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
16.
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. 相似文献
17.
非连续变形分析（DDA）是一种隐式求解的动力学计算方法,且采用在块体界面加减刚硬弹簧的方式来满足块体界面无张拉和无嵌入的接触准则,其中时间步长和弹簧刚度两个物理量的取值直接影响DDA的计算结果。基于对DDA时间步和弹簧刚度在程序运行过程中的调整策略和块体接触的简化力学概念模型,研究了惯性力在DDA收敛求解中的作用过程。采用数值模拟试验对自由落体和斜面单滑块模型在3种力学状态下的相关力学问题进行了数值模拟研究,通过对自由落体运动的模拟,研究了时间步长单一因素对计算结果的影响规律,并初步确定了时间步长的合理取值区间。在此基础上,采用斜面单滑块模型,研究了时间步长和弹簧刚度对计算结果的共同影响,确定了不同时间步长条件下弹簧刚度的合理取值区间。研究成果表明,合适的时间步长和弹簧刚度的取值组合构成一个单连通参数取值域,当时间步和弹簧刚度的取值组合位于此“域”范围内时,DDA的计算结果是合理的。 相似文献
18.
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. 相似文献
19.
《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. 相似文献
20.
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. 相似文献