Abstract
Standard higher order finite elements often perform unsatisfactory in contact problems. The major difficulties are caused by uneven distribution of nodal forces resulting in oscillating contact pressures. The paper presents a new approach that eliminates this drawback. The weight functions are chosen in such a way that even distributions of nodal forces are obtained. It is achieved by applying piece-wise linear functions. Two new 2D isoparametric quadratic elements are derived: 6-node triangle and 8-node quadrilateral, and tested in many examples. The new elements have unsymmetric stiffness matrices, but the provided examples show their good performance in contact problems.
References
[1] ABAQUS User's Manual, Version 5.7. Hibbitt, Karlsson & Sorensen, Inc, 1997.
[2] K. J. Bathe. Finite Element Procedures in Engineering Analysis. Prentice Hall, pp. 341- 379, 1982.
[3] R. Buczkowski and M. Kleiber. Elasto-plastic interface model for 3D-frictional orthotropic contact problems. Int. J. Numer. Methods Eng., 40: 599- 619, 1997.
[4] R. D. Cook. Concepts and Applications of Finite Element analysis. John Wiley & Sons, pp. 181- 183, 311-315, 1981.
[5] A. Jarzębowski, J. Maciejewski. Some remarks on the applicability of rectangular elements to plane strain boundary value problems. Computer Assisted Mechanics and Engineering Sciences, 5: 299-309, 1998.
[2] K. J. Bathe. Finite Element Procedures in Engineering Analysis. Prentice Hall, pp. 341- 379, 1982.
[3] R. Buczkowski and M. Kleiber. Elasto-plastic interface model for 3D-frictional orthotropic contact problems. Int. J. Numer. Methods Eng., 40: 599- 619, 1997.
[4] R. D. Cook. Concepts and Applications of Finite Element analysis. John Wiley & Sons, pp. 181- 183, 311-315, 1981.
[5] A. Jarzębowski, J. Maciejewski. Some remarks on the applicability of rectangular elements to plane strain boundary value problems. Computer Assisted Mechanics and Engineering Sciences, 5: 299-309, 1998.
