LAUSR

In the paper, the smoothed finite element method (S-FEM) based on linear triangular elements is used to solve 2D solid contact problems for functionally graded materials. Both conforming and nonconforming contacts algorithms are developed using modified Coulomb friction contact models including tangential strength and normal adhesion. Based on the smoothed Galerkin weak form, the system stiffness matrices are created using the formulation procedures of node-based S-FEM (NS-FEM) and edge-based S-FEM (ES-FEM), and the contact interface equations are discretized by contact point-pairs. Then these discretized system equations are converted into a form of linear complementarity problems (LCPs), which can be further solved efficiently using the Lemke method. The singular value decomposition method is used to deal with the singularity of the stiffness matrices in the procedure constructing the standard LCP, which can greatly improve the stability and accuracy of the numerical results. Numerical examples are presented to investigate the effects of the various parameters of functionally graded materials and comparisons have been made with reference solutions and the standard FEM. The numerical results demonstrate that the strain energy solutions of ES-FEM have higher convergence rate and accuracy compared with that of NS-FEM and FEM for functionally graded materials through the present contact analysis approach.