A finite element algorithm for incremental analysis of large three-dimensional frictional contact problems of linear elasticity

1993 ◽  
Vol 46 (4) ◽  
pp. 669-677 ◽  
Author(s):  
G. Zboiński
1999 ◽  
Vol 66 (2) ◽  
pp. 460-467 ◽  
Author(s):  
S. H. Ju ◽  
R. E. Rowlands

A three-dimensional contact element based on the penalty function method has been developed for contact frictional problems with sticking, sliding, and separation modes infinite element analysis. A major advantage of this contact element is that its stiffness matrix is symmetric, even for frictional contact problems which have extensive sliding. As with other conventional finite elements, such as beam and continuum elements, this new contact element can be added to an existing finite element program without having to modify the main finite element analysis program. One is therefore able to easily implement the element into existing nonlinear finite element analysis codes for static, dynamic, and inelastic analyses. This element, which contains one contact node and four target nodes, can be used to analyze node-to-surface contact problems including those where the contact node slides along one or several target surfaces.


1995 ◽  
Vol 30 (3) ◽  
pp. 195-203 ◽  
Author(s):  
M P Whelan ◽  
J E Mottershead ◽  
P D Edwards ◽  
E G Little

This study deals with the development of a specialized finite element algorithm suitable for a non-linear contact stress analysis of a model of a plastic tibial plateau of a typical unicondylar knee prosthesis. The principle feature of the contact algorithm is the use of Lagrange multiplier methods for the application of displacement constraints to surface Gauss points of a contacting body to prevent mesh overlap. This allows the effective modelling of two- and three-dimensional contact problems, with or without friction, using higher order elements. Through the selection of suitable benchmark tests, the performance and accuracy of the algorithm was assessed prior to the model analysis. Good agreement was obtained between the finite element results for the contact model and existing theoretical and experimental data. It was found that the Hertzian theory failed to accurately predict localized stresses at the contact interface when the indenter was much stiffer than the model.


Author(s):  
W. Habchi ◽  
J. Issa

This paper presents a reduced full-system finite element solution of isothermal elastohydrodynamic (EHD) line contact problems. The proposed model is based on a full-system finite element resolution of the EHL equations: Reynolds, linear elasticity and load balance. A reduced model is proposed for the linear elasticity problem. For this, three different techniques are tested: the classical “Modal reduction” and “Ritz-vector” methods and a novel “EHL-basis” method. The reduction order in the first two appears to be insufficient and a large number of degrees of freedom is required in order to attain an acceptable solution. On the other hand, the “EHL-basis” method shows up to be much more efficient, requiring only a few degrees of freedom to compose the elastic deformation of the solid components. In addition, a comparison with the full model shows an order of magnitude cpu time gain with errors of the order of only 1‰ for the central and minimum film thicknesses.


Sign in / Sign up

Export Citation Format

Share Document