This paper presents a fast numerical solver for a nonlinear constrained optimization problem, arising from 3D concentrated frictional shift and rolling contact problems with dry Coulomb friction. The solver combines an active set strategy with a nonlinear conjugate gradient method. One novelty is to consider the tractions of each slip element in a polar coordinate system, using azimuth angles as variables instead of conventional traction variables. The new variables are scaled by the diagonal of the underlying Jacobian. The fast Fourier transform (FFT) technique accelerates all matrix–vector products encountered, exploiting the matrix’ Toeplitz structure. Numerical tests demonstrate a significant reduction of the computational time compared to existing solvers for concentrated contact problems.
Additional Metadata
Keywords Frictional contact problem, Azimuth angle variables, Active set strategy, Nonlinear conjugate gradient, Preconditioner, Fast Fourier transform
THEME Null option (theme 11)
Publisher Elsevier
Stakeholder Unspecified
Persistent URL dx.doi.org/10.1016/j.jcp.2015.02.016
Journal Journal of Computational Physics
Citation
Zhao, J, Vollebregt, E.A.H, & Oosterlee, C.W. (2015). A fast nonlinear conjugate gradient based method for 3D concentrated frictional contact problems . Journal of Computational Physics, (288), 86–100. doi:10.1016/j.jcp.2015.02.016