# ContactLink3D

(New page: == General description of the condition == This conditions represents the main part of the algorithm to solve contact problems. The algorithm is based on a surface to surface penalty formu...) |
(→Formulation of the two body contact problem) |
||

Line 8: | Line 8: | ||

[[Image:Structural Appl Contact Figure 2.gif|thumb|right|300px|Slave <math>\Gamma_1</math> and master <math>\Gamma_2</math> contact surface, and gap <math>g=\nu\cdot\left(\boldsymbol{y}_2-\boldsymbol{y}_1\right)</math> between the surfaces]] | [[Image:Structural Appl Contact Figure 2.gif|thumb|right|300px|Slave <math>\Gamma_1</math> and master <math>\Gamma_2</math> contact surface, and gap <math>g=\nu\cdot\left(\boldsymbol{y}_2-\boldsymbol{y}_1\right)</math> between the surfaces]] | ||

[[Image:Structural Appl Contact Figure 3.gif|thumb|right|300px|Three dimensional view on the contact surfaces and the gap between the two surfaces]] | [[Image:Structural Appl Contact Figure 3.gif|thumb|right|300px|Three dimensional view on the contact surfaces and the gap between the two surfaces]] | ||

+ | The contact of two (or more) bodies in three dimension is a constrained continuum mechanical problem. The surface <math>\partial\Omega =\Gamma</math> of the contacting continua <math>\Omega_1</math> and <math>\Omega_2</math> are constructed by the Neumann boundary <math>\Gamma^{\sigma}</math>, the Dirichlet | ||

+ | boundary <math>\Gamma^{u}</math> and the contact boundary <math>\Gamma^c</math>. For the treated area <math>\Omega</math> the following assumptions hold: | ||

+ | |||

+ | *<math>\Omega = \Omega_1 \cup \Omega_2</math> and <math>\Omega_1 \cap \Omega_2=0</math>, | ||

+ | |||

+ | *having the boundaries <math>\partial\Omega=\Gamma^{\sigma}_i\cup\Gamma^{u}_i\cup\Gamma^{c}_i</math> and <math>\Gamma^{\sigma}_i\cap \Gamma^{u}_i=0</math>, <math>\Gamma^{\sigma}_i\cap \Gamma^{c}_i=0</math>, <math>\Gamma^{u}_i\cap \Gamma^{c}_i=0</math> | ||

+ | |||

+ | The problem is constrained by the prescription of a penetration of the two bodies, e.g. of the penetration of the master surface <math>\Gamma^2</math> by the slave surface <math>\Gamma^1</math>. This constrained is expressed by the gap function defined as the scalar projection of the distance vector between the point on the slave surface <math>\mathbf{x}\in \Gamma^1</math> and its closest point projection onto the slave surface <math>\mathbf{y}(\mathbf{x})=\mbox{arg}\;\mbox{min}_{\mathbf{y}\epsilon\Gamma_2^c}\| \mathbf{x}-\mathbf{y}\left(\mathbf{x}\right)\|</math> and the normal vector <math>\nu</math> on the master surface: | ||

+ | *<math>g\left(\mathbf{x}\right)=-\mathbf{\nu}\left(\mathbf{x}-\mathbf{y}\right)</math> | ||

+ | :::<math>=-\mathbf\nu\left(\left[\mathbf{X}+\mathbf{u}_1\right]-\left[\mathbf{Y}-\mathbf{u}_2\right]\right)</math> | ||

+ | :::<math>=g_0-\mathbf\nu\left(\mathbf{u}_1-\mathbf{u}_2\right)</math> | ||

+ | The constrained that defines the contact problem is now stated as: | ||

+ | *<math>g(x)\leq 0</math> | ||

+ | Within the finite element formulation the problem is now solved by use of an penalty treatment of the contact constrained. To solve the problem the Energie of the system must be minimised by fulfilling the penetration constraint. After having introduced a contact stress <math>t_N</math> that acts on the surfaces if they are in contact the constraint can be rewritten in form of Kuhn-Tucker constraints: | ||

+ | *<math>g\leq 0\;\;\;\;\;\;t_N\geq 0\;\;\;\;\;\;g t_N=0\;\;\;\forall\;\mathbf{x}\epsilon\Gamma^c</math> | ||

+ | The penalty method now penalizes violation of the constraint, leading to the following potential: | ||

+ | *<math>\Pi^c=\sum_{i=1}^2\int_{\Gamma_c^{i}}<\epsilon_Ng>^2d\Gamma^i</math> | ||

+ | The minimization of the potential is now been done by its variation: | ||

+ | *<math>\delta W^c=\sum_{i=1}^2\int_{\Gamma_c^{i}}\delta g\epsilon_N<g>d\Gamma^i</math> | ||

+ | Together with the minimisation of the system energy | ||

+ | *<math>\delta W^{int, ext}=\sum_{i=1}^2\delta W_i^{int, ext}</math> | ||

+ | the minimisation of the system energie with fulfillment of the penetration constraint can be written as | ||

+ | *<math>\delta W=\delta W^{int, ext}+\delta W^c=0</math> | ||

+ | This equation states a initial boundary value problem and can be solved by normal discretization methods in time and space. | ||

+ | |||

==References== | ==References== | ||

* T.A. Laursen ''Computational Contact and Impact Mechanics'', Springer, 2002 | * T.A. Laursen ''Computational Contact and Impact Mechanics'', Springer, 2002 |

## Revision as of 12:19, 2 May 2008

## General description of the condition

This conditions represents the main part of the algorithm to solve contact problems. The algorithm is based on a surface to surface penalty formulation using a master and a slave contact surface. It is formulated to solve frictional contact problems in structural mechanics, the algorithm is formulated in a finite deformation context.

This condition and the underlying algorithm for contact problems has not been tested by means of a broader range of benchmark examples so far. It seem tio work quite fine for a number of different examples, but I do not give a guarantee that it is free of errors in its formulation or its implementation. |
---|

## Formulation of the two body contact problem

The contact of two (or more) bodies in three dimension is a constrained continuum mechanical problem. The surface of the contacting continua Ω_{1} and Ω_{2} are constructed by the Neumann boundary Γ^{σ}, the Dirichlet
boundary Γ^{u} and the contact boundary Γ^{c}. For the treated area Ω the following assumptions hold:

- and ,

- having the boundaries and , ,

The problem is constrained by the prescription of a penetration of the two bodies, e.g. of the penetration of the master surface Γ^{2} by the slave surface Γ^{1}. This constrained is expressed by the gap function defined as the scalar projection of the distance vector between the point on the slave surface and its closest point projection onto the slave surface and the normal vector ν on the master surface:

The constrained that defines the contact problem is now stated as:

Within the finite element formulation the problem is now solved by use of an penalty treatment of the contact constrained. To solve the problem the Energie of the system must be minimised by fulfilling the penetration constraint. After having introduced a contact stress *t*_{N} that acts on the surfaces if they are in contact the constraint can be rewritten in form of Kuhn-Tucker constraints:

The penalty method now penalizes violation of the constraint, leading to the following potential:

The minimization of the potential is now been done by its variation:

Together with the minimisation of the system energy

the minimisation of the system energie with fulfillment of the penetration constraint can be written as

- δ
*W*= δ*W*^{int,ext}+ δ*W*^{c}= 0

This equation states a initial boundary value problem and can be solved by normal discretization methods in time and space.

## References

- T.A. Laursen
*Computational Contact and Impact Mechanics*, Springer, 2002