1.2. CONTAC12 - 2D Point-to-Point Contact

Matrix or VectorShape Functions Integration Points
Stiffness MatrixNone (nodes may be coincident)None
Load TypeDistribution
Element TemperatureNone - average used for material property evaluation
Nodal TemperatureNone - average used for material property evaluation

1.2.1. Element Matrices

CONTAC12 may have one of three conditions if the elastic Coulomb friction option (KEYOPT(1) = 0) is used: closed and stuck, closed and sliding, or open. The following matrices are derived assuming that θ is input as 0.0.

  1. Closed and stuck. This occurs if:

    (1–21)

    where:

    μ = coefficient of friction (input as MU on TB command with Lab = FRIC or MP command)
    Fn = normal force across gap
    Fs = sliding force parallel to gap

    The normal force is:

    (1–22)

    where:

    kn = normal stiffness (input as KN on R command
    un,I = displacement of node I in normal direction
    un,J = displacement of node J in normal direction
    d = distance between nodes

    The sliding force is:

    (1–23)

    where:

    ks = sticking stiffness (input as KS on R command)
    us,I = displacement of node I in sliding direction
    us,J = displacement of node J in sliding direction
    uo = distance that nodes I and J have slid with respect to each other

    The resulting element stiffness matrix (in element coordinates) is:

    (1–24)

    and the Newton-Raphson load vector (in element coordinates) is:

    (1–25)

  2. Closed and sliding. This occurs if:

    (1–26)

    In this case, the element stiffness matrix (in element coordinates) is:

    (1–27)

    and the Newton-Raphson load vector is the same as in Equation 1–25. If the unsymmetric option is chosen (NROPT,UNSYM), then the stiffness matrix includes the coupling between the normal and sliding directions; which for STAT = 2 is:

    (1–28)

  3. Open - When there is no contact between nodes I and J. There is no stiffness matrix or load vector.

Figure 1.2: Force-Deflection Relations for Standard Case shows the force-deflection relationships for this element. It may be seen in these figures that the element is nonlinear and therefore needs to be solved iteratively. Further, since energy lost in the slider cannot be recovered, the load needs to be applied gradually.

Figure 1.2: Force-Deflection Relations for Standard Case

Force-Deflection Relations for Standard Case

1.2.2. Orientation of the Element

The element is normally oriented based on θ (input as THETA on R command). If KEYOPT(2) = 1, however, θ is not used. Rather, the first iteration has θ equal to zero, and all subsequent iterations have the orientation of the element based on the displacements of the previous iteration. In no case does the element use its nodal coordinates.

1.2.3. Rigid Coulomb Friction

If the user knows that a gap element will be in sliding status for the life of the problem, and that the relative displacement of the two nodes will be monotonically increasing, the rigid Coulomb friction option (KEYOPT(1) = 1) can be used to avoid convergence problems. This option removes the stiffness in the sliding direction, as shown in Figure 1.3: Force-Deflection Relations for Rigid Coulomb Option. It should be noted that if the relative displacement does not increase monotonically, the convergence characteristics of KEYOPT(1) = 1 will be worse than for KEYOPT(1) = 0.

Figure 1.3: Force-Deflection Relations for Rigid Coulomb Option

Force-Deflection Relations for Rigid Coulomb Option