Small-sliding interaction between bodies

Although four internal element types are used to model the various types of small-sliding contact interactions supported by Abaqus/Standard, all four formulations are based on the notion that a given slave node always interacts with the same subset of master surface nodes. This nodal subset is initially determined by the Abaqusanalysis input file processor from the undeformed model definition, thus avoiding the need to “track” the slave node during the course of the analysis. This set of nearest-neighbor nodes to the point on the master surface closest to the slave node is used to parametrize a contact plane with which the slave node will interact during the analysis. This concept is illustrated next for the case of a two-dimensional slave node interacting with a first-order master surface. This formulation can be generalized to second order as well as three-dimensional situations, but this generalization will not be discussed here.

Consider the contact interaction of three nodes—102, 103, and 104—on the slave surface with a master surface made up of first-order element faces described by nodes 1, 2, and 3, as shown in Figure 1.

Figure 1. Slave nodes interacting with a two-dimensional master surface.

Before initiating the search for the nodal subset of the master surface nodes that will interact with each node on the slave surface, unit normal vectors are computed for all the nodes on the master surface. For example, the unit normal vector ${\mathbf{N}}_2$ is computed by averaging the unit normal vectors of segments 1–2 and 2–3. The user can also specify the normal vector for each node on the master surface. Additional unit normal vectors are computed for each segment a distance $\epsilon L$ from each end of the segment, where $\epsilon$ is a fraction and $L$ is the length of the segment; e.g., ${\mathbf{N}}_{\epsilon }$. Currently the value of $\epsilon$ is set to 0.5, and the user cannot change this value. The unit normals computed are then used to define a smooth varying normal vector, $\mathbf{N}\left(\mathbf{X}\right)$, at any point, $\mathbf{X}$, on the master surface.

An “anchor” point on the master surface, ${\mathbf{X}}_0$, is computed for each slave node so that the vector formed by the slave node and ${\mathbf{X}}_0$ coincides with the normal vector $\mathbf{N}\left({\mathbf{X}}_0\right)$. Suppose that a search for the anchor point, ${\mathbf{X}}_0$, of slave node 103 reveals that ${\mathbf{X}}_0$ is on segment 1–2. Then, we find that

where ${\mathbf{X}}_1$ and ${\mathbf{X}}_2$ are the coordinates of nodes 1 and 2, respectively, and $u_0$ is calculated so that ${\mathbf{X}}_0-{\mathbf{X}}_{103}$ coincides with $\mathbf{N}\left({\mathbf{X}}_0\right)$. Moreover, the contact plane tangent direction, ${\mathbf{v}}_0$, at ${\mathbf{X}}_0$ is chosen so that it is perpendicular to $\mathbf{N}\left({\mathbf{X}}_0\right)$; i.e.,

where $\mathbf{T}$ is a (constant) rotation matrix.

The small-sliding contact constraint is achieved by requiring that slave node 103 interact with the tangent plane whose current anchor point coordinates are, at any time, given by

where $N_1\left(u_0\right)=1-u_0$ and $N_2\left(u_0\right)=u_0$, and whose current tangent direction is given by

where $N_1^u\left(u_0\right)=-1$ and $N_2^u\left(u_0\right)=1$. Since the above expressions for the point ${\mathbf{x}}_0$ and the vector $\mathbf{v}$ resulted from barycentric (affine) combinations of the points ${\mathbf{x}}_1$ and ${\mathbf{x}}_2$—that is,

the contact plane will be mapped properly under affine transformations such as translation, scaling (stretching), and rotation.

Next, suppose that a search for the anchor point of slave node 104 reveals that the anchor point is coincident to the master node 2. In this case the anchor point is chosen to be ${\mathbf{X}}_2$, or in terms of the coordinates of the three master nodes 1, 2, and 3,

where $N_1\left(u_0\right)=0$, $N_2\left(u_0\right)=1$, and $N_3\left(u_0\right)=0$. The contact tangent direction at ${\mathbf{X}}_2$ is simply

However, we want to express ${\mathbf{v}}_0$ in terms of the coordinates of the three nodes 1, 2, and 3 to be able to track the evolution of the tangent plane. To this end, we solve for the $N_i^u$ from the equation

subject to the barycentric constraint

The barycentric constraint ensures that the resulting expression for the contact plane tangent direction behaves properly under affine transformations such as translations and rotations.

At each slave node that can come into contact with a master surface we construct a measure of overclosure $h$ (penetration of the node into the master surface) and measures of relative slip $s_i$. These kinematic measures are then used, together with appropriate Lagrange multiplier techniques, to introduce surface interaction theories for contact and friction, as described in Contact pressure definition and Coulomb friction.

In two dimensions the overclosure $h$ along the unit contact normal $\mathbf{n}$ between a slave point ${\mathbf{x}}_{N+1}$ and a master line $\mathbf{p}\left(\xi \right)$, where $\xi$ parametrizes the line, is determined by finding the vector $\left(\mathbf{p}-{\mathbf{x}}_{N+1}\right)$ from the slave node to the line that is perpendicular to the tangent vector $\mathbf{v}$ at $\mathbf{p}$. Mathematically, we express the required condition as

when

Similarly, in three dimensions the overclosure $h$ along the unit contact normal $\mathbf{n}$ between a slave point ${\mathbf{x}}_{N+1}$ and a master plane $\mathbf{p}\left({\xi }_1,{\xi }_2\right)$, where ${\xi }_i$ parametrize the plane, is determined by finding the vector $\left(\mathbf{p}-{\mathbf{x}}_{N+1}\right)$ from the slave node to the plane that is perpendicular to the tangent vectors ${\mathbf{v}}_1$ and ${\mathbf{v}}_2$ at $\mathbf{p}$. Mathematically, we express the required condition as

when

If at a given slave node $h<0$, there is no contact between the surfaces at that node, and no further surface interaction calculations are needed. If $h\ge 0$, the surfaces are in contact. The contact constraint $h=0$ is enforced by introducing a Lagrange multiplier, $\tilde{p}$, whose value provides the contact pressure at the point. To enforce the contact constraint, we need the first variation $\delta h$; and for the Newton iterations, we need the second variation, $d\delta h$. Likewise, if frictional forces are to be transmitted across the contacting surfaces, the first variations of relative slip, $\delta s_i$, and the second variations, $d\delta s_i$, are needed in the formulation. The derivation of some of these quantities is described next for all four small-sliding contact formulations.

For the case of two-dimensional small-sliding deformable contact, a point on the contact line associated with a slave node $n_{N+1}$ is represented by the vector

where, as described previously, the line's anchor point—${\mathbf{x}}_0$—and its tangent vector—$\mathbf{v}$—are functions of the current master node coordinates, ${\mathbf{x}}_1\mathrm{\dots }{\mathbf{x}}_N$. Hence, the vector $\mathbf{v}$ is, in general, a nonunit vector. Linearization of Equation 1 yields

where $\delta {\mathbf{u}}_{N+1}=\delta \left({\mathbf{X}}_{N+1}+{\mathbf{u}}_{N+1}\right)=\delta {\mathbf{x}}_{N+1}$, $\delta {\mathbf{x}}_0=N_i\left(u_0\right)\delta {\mathbf{x}}_i$, and $\delta \mathbf{v}=\mathbf{R}\cdot N_i^u\left(u_0\right)\delta {\mathbf{x}}_i$.

Taking the dot product of Equation 4 with $\mathbf{n}$ results in the following expression for $\delta h$:

Likewise, taking the dot product of Equation 4 with the normalized tangent vector $\mathbf{t}=\mathbf{v}/\parallel \mathbf{v}\parallel$ and setting $h=0$ results in the following expression for the variation in slip:

Suitable expressions for $d\delta h$ and $d\delta s$ can be derived by linearizing Equation 4 and applying the techniques highlighted above. Since the resulting expressions do not provide additional insight into the understanding of this capability, they will not be presented here.

The three-dimensional small-sliding deformable contact formulation is a straightforward generalization of the previous two-dimensional formulation. A point on the contact plane associated with a slave node $n_{N+1}$ is represented by the vector

where the plane's anchor point—${\mathbf{x}}_0$—and its two tangent vectors—${\mathbf{v}}_1$ and ${\mathbf{v}}_2$—are functions of the current master node coordinates ${\mathbf{x}}_1,\mathrm{\dots },{\mathbf{x}}_N$. Linearization of Equation 2 yields

where $\delta {\mathbf{x}}_0=\mathbf{f}\left(\delta {\mathbf{x}}_1,\mathrm{\dots },\delta {\mathbf{x}}_N\right)$ and $\delta {\mathbf{v}}_i={\mathbf{g}}_i\left(\delta {\mathbf{x}}_1,\mathrm{\dots },\delta {\mathbf{x}}_N\right)$.

Taking the dot product of Equation 7 with $\mathbf{n}$ results in the following expression for $\delta h$:

Likewise, taking the dot product of Equation 7 with ${\mathbf{t}}_1={\mathbf{v}}_1/\parallel {\mathbf{v}}_1\parallel$ and setting $h=0$ results in the following expression for the variation of the first slip component:

Similarly, taking the dot product of Equation 7 with ${\mathbf{t}}_2={\mathbf{v}}_2/\parallel {\mathbf{v}}_2\parallel$ and setting $h=0$ results in the following expression for the variation of the second slip component:

The formulation for two-dimensional small-sliding rigid contact follows from its deformable counterpart by exploiting the fact that the evolution of the contact plane is fully determined by the motion of the rigid body's reference node. Figure 2 shows how the undeformed coordinates ${\mathbf{X}}_0$ of the contact plane's anchor point are related vectorially to the undeformed coordinates of the rigid reference node, ${\mathbf{X}}_{rs}$, and the relative position vector $\mathbf{R}$.We can express this relationship as

Figure 2. Rigid body reference geometry.

Suppose the rigid reference node undergoes a motion described by the displacement vector ${\mathbf{u}}_{rs}$ and the rotation vector ${\varphi }_{rs}{\mathbf{e}}_z$, then the current coordinates of the contact plane's anchor point are given by

where $\mathbf{C}$ is the orthogonal matrix that produces the rotation ${\varphi }_{rs}{\mathbf{e}}_z$ (see Rotation variables) and $\mathbf{r}$ is the current position vector from the rigid reference node to the anchor point. The rotation matrix is also used to obtain the contact plane's current tangent and current normal as follows:

where ${\mathbf{t}}_0$ and ${\mathbf{n}}_0$ are the initial contact tangent and normal at ${\mathbf{X}}_0$, respectively, as shown in Figure 2. By definition, a rigid surface cannot stretch; therefore, a point on the rigid contact line associated with a slave node $n_{N+1}$ is represented by the vector

where—unlike the vector $\mathbf{v}$ in Equation 3—the tangent $\mathbf{t}$ is always a unit vector.

Linearization of Equation 11 and Equation 12 results in the following expressions for the first variations of the anchor coordinates and the contact plane's tangent:

Replacing $\mathbf{v}$ by $\mathbf{t}$ in Equation 5, noting that $\mathbf{t}\cdot \delta \mathbf{t}=0$, and substituting for $\delta {\mathbf{x}}_0$ and $\delta \mathbf{t}$ from Equation 13 results in the following expression for $\delta h$:

Similar manipulation of Equation 6 yields the following expression for $\delta s$:

The three-dimensional small-sliding rigid contact formulation can also be derived from its deformable counterpart by generalizing some of the expressions that were introduced in the previous section. In particular, in this formulation the rigid reference node can undergo an arbitrary finite rotation described by the rotation vector ${\mathit{\varphi }}_{rs}$. Consequently, Equation 13 for the variations of the anchor point coordinates and the contact plane's tangent generalize to

where ${\widehat{\delta \mathit{\theta }}}_{rs}$ is the skew-symmetric matrix associated with the linearized rotation $\delta {\mathit{\theta }}_{rs}$, as explained in Rotation variables.

Replacing ${\mathbf{v}}_i$ by ${\mathbf{t}}_i$ in Equation 8 through Equation 10 and substituting for $\delta {\mathbf{x}}_0$ and $\delta {\mathbf{t}}_i$ from above results in the following expressions for $\delta h$, $\delta s_1$, and $\delta s_2$: