Rigid Body Contact (TYPE6)

It works like TYPE3 interface except that the total interface force is a user-defined function of the maximum penetration. The input and computational algorithms are the same as for TYPE3 interfaces. This interface is used extensively in vehicle occupant simulations; example, knee bolsters.

Limitations

Some of the main limitations for this interface type are:

Surface 1 must be part of one and only one rigid body.
Surface 2 must be part of one and only one rigid body.
The interface stiffness (user-defined function) can reduce the time step.

Other limitations are the same as for TYPE3 interfaces.

Interface Stiffness

When two surfaces contact, a massless stiffness is introduced to reduce the penetration's nodes of the other surface into the surface.

The nature of the stiffness depends on the type of interface and the elements involved.

If a segment is a shell as well as the face of brick element, the shell stiffness is used.

Interface Friction

TYPE6 interface allows sliding between contact surfaces. Coulomb friction between the surfaces is modeled. The input card requires a friction coefficient. No value (default) defines zero friction between the surfaces. The friction computation on a surface is the same as for TYPE3 interface (refer to Contact Processing).

Interface Gap

Refer to Contact Detection for TYPE3 interface.

Time Step Calculation

The stable time step used for time integration equations is computed by:(1)

Δ t = 0.1 \sqrt{\frac{M}{K}}

$Δ t = 0.1 \sqrt{\frac{M}{K}}$

Where,

$M$ $M$: min(M rigid body 1, M rigid body 2)
$K$ $K$: Tangent of user force function

Time step

Δ t

$Δ t$ is affected by the actual stiffness derived from function

f

$f$ :(2)

Δ t \leq \sqrt{\frac{M}{f^{'} (p)}}

$Δ t \leq \sqrt{\frac{M}{f^{'} (p)}}$

The function $f$ $f$ refers to a function number given in input and must be user-defined.

Contact Force

(3)

F_{N_{i}} = \frac{- {\vec{p}}_{i} + {\vec{q}}_{N}_{i}}{‖ \sum {\vec{p}}_{j} - \sum {\vec{q}}_{k} ‖} f (M A X (‖ {\vec{p}}_{k} ‖, ‖ {\vec{q}}_{j} ‖))

$F_{N_{i}} = \frac{- {\vec{p}}_{i} + {\vec{q}}_{N}_{i}}{‖ \sum {\vec{p}}_{j} - \sum {\vec{q}}_{k} ‖} f (M A X (‖ {\vec{p}}_{k} ‖, ‖ {\vec{q}}_{j} ‖))$

${\vec{q}}_{N_{i}}$ ${\vec{q}}_{N_{i}}$ is the contribution to node $N_{i}$ $N_{i}$ of vector ${\vec{q}}_{i}$ ${\vec{q}}_{i}$ distributed on the segment penetrated by node $Q_{i}$ $Q_{i}$ .