期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 257, 期 -, 页码 333-351出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2013.10.005
关键词
Contact mechanics; Active set strategy; Conjugate gradients method; Preconditioners; Toeplitz matrices; Fast Fourier transform
资金
- VORtech BV (Mark Roest)
This paper presents our new solver BCCG+FAI for solving elastic normal contact problems. This is a comprehensible approach that is based on the Conjugate Gradients (CG) algorithm and that uses FFTs. A first novel aspect is the definition of the FFT-based Approximate Inverse preconditioner. The underlying idea is that the inverse matrix can be approximated well using a Toeplitz or block-Toeplitz form, which can be computed using the FFT of the original matrix elements. This preconditioner makes the total number of CG iterations effectively constant in 2D and very slowly increasing in 3D problems. A second novelty is how we deal with a prescribed total force. This uses a deflation technique in such a way that CGs convergence and finite termination properties are maintained. Numerical results show that this solver is more effective than existing CG-based strategies, such that it can compete with Multi-Grid strategies over a much larger problem range. In our opinion it could be the new method of choice because of its simple structure and elegant theory, and because robust performance is achieved independently of any problem specific parameters. (C) 2013 Elsevier Inc. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据