# Discontinuous Deformation Analysis

**Discontinuous Deformation Analysis** (DDA) is a type of discrete element method (DEM) originally proposed by Shi^{[1]} in 1988. DDA is somewhat similar to the finite element method for solving stress-displacement problems, but accounts for the interaction of independent particles (blocks) along discontinuities in fractured and jointed rock masses. DDA is typically formulated as a work-energy method, and can be derived using the principle of Minimum Potential Energy^{[1]} or by using Hamilton's principle. Once the equations of motion are discretized, a step-wise linear time marching scheme in the Newmark family is used for the solution of the equations of motion. The relation between adjacent blocks is governed by equations of contact interpenetration and accounts for friction. DDA adopts a stepwise approach to solve for the large displacements that accompany discontinuous movements between blocks. The blocks are said to be "simply deformable". Since the method accounts for the inertial forces of the blocks' mass, it can be used to solve the full dynamic problem of block motion.

## DDA vs DEM

Although DDA and DEM are similar in the sense that they both simulate the behavior of interacting discrete bodies, they are quite different theoretically. While DDA is a displacement method, DEM is a force method. While DDA uses displacement as variables in an implicit formulation with opening-closing iterations within each time step to achieve equilibrium of the blocks under constrains of the contact, DEM employs an explicit, time marching scheme to solve the equations of motion directly (Cundall and Hart^{[2]}). The system of equation in DDA is derived from minimizing the total potential energy of the system being analyzed. This guarantee that equilibrium is satisfied at all times and that energy consumption is natural since it is due to frictional forces. In DEM, unbalanced forces drive the solution process, and damping is used to dissipate energy. If a quasi-static solution is desired in which the intermediate steps are not of interest, the type of damping and the type of relaxation scheme can be selected in DEM to obtain the most efficient solution method (Cundall^{[3]}). The application of damping in DEM for quasi-static problem is somewhat analogues to the setting to zero of the initial velocities of the block in the static analysis of DDA. In dynamic problem, however, the amount and type of damping in DEM, which are very difficult to qualify experimentally, has to be selected very carefully to as not to damp out real vibrations. On the other hand, the energy consumption in DDA is due to the frictional resistance at contact. By passing the velocities of the blocks at the end of a time step to the next time step, DDA gives real dynamic solution with correct energy consumption.^{[1]} By using an energy approach, DDA does not require an artificial damping term to dissipate energy as in DEM, and can easily incorporate other mechanisms for energy loss.

## Strengths and Limitations of the DDA method

DDA has several strengths recommending it for use in slope stability problems in jointed rock masses, which are balanced by serious limitations be accounted when DDA is used for larger scale, faster moving problems.

### Strengths

- Very good for problems with small characteristic as time marching scheme provides necessary numerical damping to control resonance interactions within and between particles.
- Step-wise linear implicit time marching allows so-called
*quasi-static*solutions, where step-wise velocities are never used. Quasi-static analysis is useful for examining slow, or creeping failures.

### Limitations

- The most serious limitation of the DDA method is the reduction of numerical damping which occurs as the characteristic length of a problem grows. The numerically, damping is a function of . Typically,

the stiffness doesn't vary over 1 or 2 orders of magnitude, while the mass is a function of the cube of the characteristic length.

## Modification and improvement of the DDA method

Various modifications to the original DDA formulation have been reported in the rock mechanics literature. In the original DDA formulation a first order polynomial displacement function was assumed, so the stresses and strains within a block in the model were constant. This approximation precludes the application of this algorithm to problems with significant stress variations within the block. However, in cases where the displacement inside the block is high and cannot be ignored, the blocks can be divided by mesh. An example of this approach is the research by Chang et al.^{[4]} and Jing^{[5]} who resolved this problem by adding finite element meshes in the two-dimensional domain of the blocks so that stress variations within the blocks can be allowed.

Higher order DDA method for two-dimensional problems has been developed in both theory and computer codes by researchers like Koo and Chern,^{[6]} Ma et al.^{[7]} and Hsiung.^{[8]} Additionally, Grayeli and Mortazavi^{[9]} implemented finite element mesh based on six node triangular element into the original DDA blocks in order to account for stress variations within the blocks. The DDA contact model which was originally based on penalty method was improved by adopting the Lagrange type approach reported by Lin et al.^{[10]}

Since a blocky system is a highly non-linear system due to non-linearity within blocks and between blocks, Chang et al.^{[4]} implemented a material non-linearity model to DDA using strain hardening curves. Ma^{[11]} developed a non-linear contact model for analysis of slope progressive failure including strain softening using the stress and strain curve.

Recent progress in DDA algorithm is reported by Kim et al.^{[12]} and Jing et al.^{[13]} which considers coupling of fluid flow in fractures. The hydro-mechanical coupling across rock fracture surfaces is also taken into account. The program computes water pressure and seepage throughout the rock mass of interest. In its original formulation, a rock bolt was modeled as a line spring connecting two adjacent blocks. Later, Te-Chin Ke^{[14]} suggested an improved bolt model, followed by the rudimentary formulation of lateral constraint of rock bolting.

Moosavi and Grayeli^{[15]} implemented a fully grouted cable bolt element into the DDA algorithm which can be considered another step towards improving the potential of DDA for solving blocky rock masses requiring reinforcements for stability purposes.

Shahami et al^{[16]} tried to analyze the blocks interaction of key-groups, first to study performance of key-groups and then to choose the best method for developing between numerical methods were used fuzzy multi criteria decision maker. At the end the these, discontinuous deformation analysis as compatible numerical method with standards developed was selected and then combining these two methods were considered and then the algorithm and the primary program was prepared in mathematica software.

## Software

### Open source and non-commercial software

- UC Berkeley DDA for Windows (2003–present) discontinuous deformation analysis open source project located in Sourceforge web site.

- : :

## References

- 1 2 3 Shi G.H. Discontinuous deformation analysis: A new numerical model for the statics and dynamics of block systems. University of California, Berkeley. 1988
- ↑ Cundall, P. A., and R. D. Hart. “Numerical Modelling of Discontinua,” in Proceedings of the 1st U.S. Conference on Discrete Element Methods (Golden, Colorado, October 1989), pp. 1-17. G.G.W. Mustoe, M. Henriksen and H-P. Huttelmaier, Eds. Golden, Colorado: CSM Press, 1989.
- ↑ Cundall, P. A. “Distinct Element Models of Rock and Soil Structure,” in Analytical and Computational Methods in Engineering Rock Mechanics, Ch. 4, pp. 129-163. E. T. Brown, Ed. London: George Allen & Unwin, 1987.
- 1 2 CHANG, C. T., MONTEIRO, P., NEMATI, K., & SHYU, K. (1996). Behavior of marble under compression. Journal of materials in civil engineering, 8(3), 157-170.
- ↑ Jing L. Formulation of discontinuous deformation analysis (DDA)—an implicit discrete element model for block systems.
*Int J Eng Geol*1998;49:371–81. - ↑ Koo CY, Chern JC. The development of DDA with third order displacement function. In: Salami MR, Banks D, editors.
*Discontinuous deformation analysis (DDA) and simulations of discontinuous media. 1996.* - ↑ Ma MY, Zaman M, Zhu JH. Discontinuous deformation analysis using the third order displacement function. In: Salami MR, Banks D, editors.
*Discontinuous deformation analysis (DDA) and simulations of discontinuous media. 1996.* - ↑ Hsiung SM. Discontinuous deformation analysis (DDA) with nth order polynomial displacement functions. 38th US rock mechanics symposium, July 7–10, Washington DC; 2001.
- ↑ Roozbeh Grayeli, Ali Mortazavi. Discontinuous deformation analysis with second-order finite element meshed block.
*Int. J. Numer. Anal. Meth. Geomech.*, 2006, 30, 1545-1561 - ↑ Lin CT, Amadei B, Jung J, Dwyer J. Extension of discontinuous deformation analysis for jointed rock masses.
*Int J Rock Mech Min Sci Geomech Abstr*1996;33:671–94. - ↑ Ma MY.
*Development of discontinuous deformation analysis, the first ten years; 1986–1996. 1999. in ICADD-3: Third International Conference on Analysis of Discontinuous Deformation---From theory to practice, pages 17-32. American Rock Mechanics Association.* - ↑ Kim Y, Amadei B, Pan E. Modelling the effect of water, excavation sequence and rock reinforcement with discontinuous deformation analysis.
*Int J Rock Mech Min Sci Geomech Abstr*1999;36:949–70. - ↑ Jing L, Ma Y, Fang Z. Modeling of fluid flow and solid deformation for fractured rocks with discontinuous deformation analysis (DDA) method.
*Int J Rock Mech Min Sci Geomech Abstr*2001; 38:343–55. - ↑ Te-Chin K. Improved modeling of rock bolting in DDA. Computer methods and advances in geomechanics; 1997.
- ↑ Moosavi, M & Grayeli, R; A model for cable bolt-rock mass interaction:integration with discontinuous deformation analysis (DDA) algorithm.
*Int J Rock Mech Min Sci Geomech Abstr*2006;43:661–670. - ↑ Shahami, M.H. & Yarahmadi, A.R. & Grayeli, R; Evaluate the capability of key-groups method development using numerical methods for finding the critical key-group in analysis of stability of rock slopes.
*Iranian Journal of Mining Engineering*

## Additional References

- Shi GH.
*Block system modeling by discontinuous deformation analysis.*Computational Mechanics Publications; 1993. - Shi GH. Discontinuous deformation analysis technical note. First international forum on discontinuous deformation analysis, June 12–14. Berkeley, California; 1996.