摘要
Contact analysis is recognized as being the most challenging problem in computational mechanics,because the functional system of contact problems is nonlinear and non-smooth and the convergence and accuracy of contact algorithms are difficult to guarantee.In the traditional finite element method(FEM)-based contact analysis[1,2],the contact body is spatially discretized,and the contact boundary is described using a low-order Lagrange interpolation polynomial.
Contact analysis is recognized as being the most challenging problem in computational mechanics, because the functional system of contact problems is nonlinear and non-smooth and the convergence and accuracy of contact algorithms are difficult to guarantee. In the traditional finite element method (FEM)-based contact analysis [1,2], the contact body is spatially discretized, and the contact boundary is described using a low-order Lagrange interpolation polynomial. While low-order Lagrange interpolation polynomial works well for planar and low-curvature surfaces [3], artificial oscillations in the contact force and numerical difficulties in the nonlinear solution algorithm occur when a large curvature or finite deformation is encountered [2].
基金
supported by the National Natural Science Foundation of China(Grant No.51421064)
the National Key Research and Development Plan of China(Grant No.2016YFB0201001)