2019/04/02 by Y. F. Zhang, J. H. Yue, Junhong Yue +4
Computer Science · Engineering · Mathematics · #Applied mathematics #Boundary element method #Boundary knot method #Contact Mechanics and Variational Inequalities #Discretization #Engineering #Extended finite element method #Finite element method #Galerkin method #Mathematical analysis #Mathematics #Mechanical stress and fatigue analysis #Mixed finite element method #Numerical analysis #Numerical methods in engineering #Singularity #Smoothed finite element method #Stiffness #Structural engineering
paper · doi:10.1142/s0219876219400127
openalex publication_date 2019/04/02 · crossref created 2019/04/02 · crossref issued 2019/06/14 · crossref published 2019/06/14 · crossref published-online 2019/06/14 · crossref published-print 2020/06/01 · crossref deposited 2020/12/05 · openalex created_date 2025/10/10 · crossref indexed 2026/08/04 · openalex updated_date 2026/08/04
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.