Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 84
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 105
Scalable Algorithms for Contact Problems with Additional Nonlinearities J. Dobiáš^{1}, S. Pták^{1}, Z. Dostál^{2} and V. Vondrák^{2}
^{1}Institute of Thermomechanics, Prague, Czech Republic
J. DobiáÂš, S. Pták, Z. Dostál, V. Vondrák, "Scalable Algorithms for Contact Problems with Additional Nonlinearities", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Fifth International Conference on Engineering Computational Technology", CivilComp Press, Stirlingshire, UK, Paper 105, 2006. doi:10.4203/ccp.84.105
Keywords: contact problem, domain decomposition, numerical scalability, geometric nonlinearity, material nonlinearity, finite element method.
Summary
Contact modelling is still a challenging problem of nonlinear
computational mechanics. The complexity of such problems is related
to the a priori unknown contact interface and contact tractions.
Their evaluation have to be part of the solution. In addition, the
solution across the contact interface is nonsmooth.
In 1991 Farhat and Roux came up with a novel domain decomposition method called the finite element tearing and interconnecting (FETI) method [1]. This method belongs to the class of nonoverlapping totally disconnected spatial decompositions. Its key concept stems from the idea that the spatial subdomains, into which the domain is partitioned, are 'glued' by Lagrange multipliers, or forces in this context. After eliminating the primal variables, which are displacements in the displacement based analysis, the original problem is reduced to a small, relatively well conditioned, typically equality constrained quadratic programming problem that is solved iteratively. The CPU time that is necessary for both the elimination and iterations can be reduced nearly proportionally to the number of the processors, so that the algorithm exhibits parallel scalability. This method has proved to be one of the most successful algorithms for parallel solution of problems governed by elliptic partial differential equations. Observing that the equality constraints may be used to define so called 'natural coarse grid', Farhat, Mandel and Roux [2] modified the basic FETI algorithm so that they were able to prove its numerical scalability, i.e. asymptotically linear complexity. If the FETI method is applied to the contact problems, the same methodology can be used to prescribe conditions of nonpenetration between bodies. We shall obtain a new minimization problem with additional nonnegativity constraints which replace more complex general nonpenetration conditions [3]. It turned out that the scalability of the FETI methods may be preserved even for the solution of contact problems [3,5]. In this paper we are concerned with application of one new variant of the FETI domain decomposition method, called TFETI (Total FETI) method, to the solution of contact problems. Both compatibility between adjacent subdomains and Dirichlet boundary conditions are enforced by Lagrange multipliers acting along the boundaries or mutual interfaces. Such an approach is very advantageous especially from the computational point of view, because the stiffness matrices of all subdomains exhibit the same defect which is, in addition, known beforehand. See [4] for discussion of this topic. We describe theoretical foundation of the TFETI algorithm and its implementation into the inner loop of the code which treats the material and geometric nonlinear effects in the outer loop. To solve a contact problem by the FETI or TFETI method, we use modified proportioning with the reduced gradient projection (MPRGP) algorithm. More details can be found in [6]. Numerical experiments were carried out with our inhouse general purpose finite element package PMD (package for machine design) [7] and they include solutions of contact problems both with and without geometric and material nonlinear effects. To demonstrate the ability of the method to solve contact problems, we computed the problem of contact of two cylindric bodies with parallel axes. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
