Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE EIGHTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Simulation of Orthotropic Frictional Contact with Non-Associated Sliding Rule
G. de Saxcé1, Z.Q. Feng2, M. Hjiaj3 and Z. Mróz4
1Laboratoire de Mécanique de Lille, Villeneuve d'Ascq, France
G. de Saxcé, Z.Q. Feng, M. Hjiaj, Z. Mróz, "Simulation of Orthotropic Frictional Contact with Non-Associated Sliding Rule", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Eighth International Conference on Computational Structures Technology", Civil-Comp Press, Stirlingshire, UK, Paper 229, 2006. doi:10.4203/ccp.83.229
Keywords: contact simulations, anisotropic friction law, non associated law, variational formulations, augmented Lagrangean technique, predictor-corrector scheme.
For many industrial applications, the assumption of isotropic friction is unrealistic because of directional surface machining and finishing operations. For orthotropic frictional contact, Michaowski and Mróz [3,4] pointed out the non-associated nature of the sliding rule. In the present work, their law is derived from a single non-differentiable function, called a bi-potential. On this ground, the time-integration of the contact law takes the form of a predictor-corrector scheme addressing all cases (sticking, sliding, no-contact) thanks to a projection onto the friction cone. Three-dimensional numerical applications show the strong influence of the slip rule.
The asperity model used by Michaowski and Mróz [3,4] to study anisotropic frictional contact phenomenon generates limit friction curves in the contact plane that can be slightly non-convex but very close to ellipses. Here, we consider only a convex friction cone
where is a diagonal matrix with diagonal elements and . The non associated slip rule
is defined by the convex slip potential
where is a diagonal matrix with diagonal elements and . Introducing the sliding non-associativity matrix , it is proved in  that the previous friction law can be put in the following compact form of a differential inclusion
where we use the dual norm and is the indicatory function of , null if and infinite otherwise. Next, we adopt an incremental approach to discretise the quasi-static evolution. Using the concept of bi-potential introduced in  and the implicit integration scheme, we deduce the value of the contact reaction at the end of the step from its value at the begining of the step and the relative displacement increment thanks to one very robust predictor-corrector step
where is a strictly positive number that need to be chosen in a suitable range to ensure convergence.
In this work, we propose a benchmark test to validate the algorithm for a class of non-associated anisotropic friction laws. The test of such frictional contact laws requires a 3D finite element model. The problem under consideration is a deformable elastic cylinder. The lower face is in contact with a rigid horizontal surface with anisotropic friction condition, while the upper face is subjected to a vertical rigid motion.
purchase the full-text of this paper (price £20)