The DDA Method

The Elements

The elements in Shi’s model have arbitrary polyhedral shapes. Each element has a reference point, which is the centre of gravity of the element. Denote the position of the reference point of element p by (x^p, y^p, z^p).

The element can translate and rotate as a rigid body (6 degrees of freedom); and in addition to that, a uniform deformation field (characterized by the usual small strain tensor) belongs to the element (another 6 degrees of freedom). The element has altogether 12 degrees of freedom, and the 12 unknown kinematic characteristics are summarized into the generalized displacement vector of the element:

Assuming small displacements, the translations of an arbitrary (x, y, z) point on the element can uniquely be determined with the help of the u^p generalized displacement vector of the element:

Note that these relations are valid for small incremental rotations only (the issue of large rotations will be discussed at the end of Section 2).

This translation vector is particularly important for those points on the element boundary where a contact is formed with a neighbouring element. The difference between the translations of the two material points belonging to the two elements that form the contact is considered the contact deformation. This relative translation determines the forces transmitted through the contact, according to the contact constitutive relations.

Note that the role of matrix T is similar to that of the transition matrix in FEM; the difference is that in DDA the deformation of the element is also taken into account.