Accurate modeling of contact using cubic splines

In this paper, we develop and implement a new method for the accurate representation of contact surfaces. This approach overcomes the difficulties arising from the use of traditional node-to-linear surface contact algorithms. In our proposed method, contact surfaces were modeled accurately using C¹-continuous cubic splines, which interpolate the finite element nodes. In this case, the unit normal vectors are defined uniquely at any point on the contact surfaces. These splines preserve the local deformation of the nodes on each flexible contact surface. Consequently, a consistent linearization of the kinematic contact constraints, based on the spline interpolation, was derived. Moreover, the gap between two contact surfaces was modeled accurately using an efficient surface-to-surface contact search algorithm. Since the continuity of the splines is not affected by the number of nodes, accurate stress distribution can be obtained with less finite elements at the contact surface than that using the traditional linear discretization of the contact surface. Two numerical examples are used to illustrate the advantages of the proposed representation. They show a significant improvement in accuracy compared to traditional piecewise element-based surface interpolation. This approach overcomes the problem of mismatch in a finite element mesh. This is very useful, since most realistic engineering problems involve contact areas that are not known a priori.