Condition Number and Clustering-Based Efficiency Improvement of Reduced-Order Solvers for Contact Problems Using Lagrange Multipliers

Abstract

This paper focuses on reduced-order modeling for contact mechanics problems treated by Lagrange multipliers. The high nonlinearity of the dual solutions lead to poor classical data compression. A hyper-reduction approach based on a reduced integration domain (RID) is considered. The dual reduced basis is the restriction to the RID of the full-order dual basis, which ensures the hyper-reduced model to respect the non-linearity constraints. However, the verification of the solvability condition, associated with the well-posedness of the solution, may induce an extension of the primal reduced basis without guaranteeing accurate dual forces. We highlight the strong link between the condition number of the projected contact rigidity matrix and the precision of the dual reduced solutions. Two efficient strategies of enrichment of the primal POD reduced basis are then introduced. However, for large parametric variation of the contact zone, the reachable dual precision may remain limited. A clustering strategy on the parametric space is then proposed in order to deal with piece-wise low-rank approximations. On each cluster, a local accurate hyper-reduced model is built thanks to the enrichment strategies. The overall solution is then deeply improved while preserving an interesting compression of both primal and dual bases.