1.19. CONTA176 - 3D Line-to-Line Contact

Matrix or VectorShape Function Integration Point Locations
Stiffness MatrixW = C1 + C2x + C3x2 None

1.19.1. Other Applicable Sections

The CONTA176 description is the same as for CONTA174 - 3D 8-Node Surface-to-Surface Contact except that it is a 3D line contact element.

1.19.2. Contact Kinematics

Three different scenarios can be modeled by CONTA176:

Use KEYOPT(3) = 0 for the first two scenarios (internal contact and parallel beams). In both cases, the contact condition is only checked at contact nodes.

Use KEYOPT(3) = 1 for the third scenario (beams that cross). In this case, the contact condition is checked along the entire length of the beams. The beams with circular cross-sections are assumed to come in contact in a point-wise manner.

Figure 1.14: Beam Sliding Inside a Hollow Beam

Beam Sliding Inside a Hollow Beam

Figure 1.15: Parallel Beams in Contact

Parallel Beams in Contact

Figure 1.16: Crossing Beams in Contact

Crossing Beams in Contact

Contact is detected when two circular beams touch or overlap each other. The non-penetration condition for beams with a circular cross-section can be defined as follows.

For internal contact:

(1–128)

and for external contact:

(1–129)

where:

g = gap distance
rc and rt = radii of the cross-sections of the beam on the contact and target sides, respectively.
d = minimal distance between the two beam centerlines (also determines the contact normal direction).

Contact occurs for negative values of g.

1.19.3. Contact Models

The contact model can be either contact force- based (KEYOPT(3) = 0 or 1) or contact traction-based (KEYOPT(3) = 2 or 3). For the contact traction-based model, the program determines the area (based on the underlying beam element length and the contact radius, R2) associated with the contact node.

1.19.4. Contact Forces

In order to satisfy contact compatibility, forces are developed in a direction normal (n-direction) to the target that will tend to reduce the penetration to an acceptable numerical level. In addition to normal contact forces, friction forces are developed in directions that are tangent to the target plane.

(1–130)

where:

Fn = normal contact force
Kn = contact normal stiffness (input as FKN on R command)
un = contact gap size

(1–131)

where:

FT = tangential contact force
KT = tangential contact stiffness (input as FKT on R command)
uT = contact slip distance

For orthotropic friction, μeq is computed using the expression:

(1–132)

where:

μeq = equivalent coefficient of friction for orthotropic friction
μ1, μ2 = coefficients of friction in first and second principal directions (input as MU1 and MU2 using TB command with Lab = FRIC)