Solution of 3D contact shape optimization problems with Coulomb friction based on TFETI

The present paper deals with the numerical solution of 3D shape optimization problems in frictional contact mechanics. Mathematical modelling of the Coulomb friction problem leads to an implicit variational inequality which can be written as a fixed point problem. Furthermore, it is known that the discretized problem is uniquely solvable for small coefficients of friction. Since the considered problem is nonsmooth, we exploit the generalized Mordukhovich’s differential calculus to compute the needed subgradient information. The state problem is solved using successive approximations combined with the Total FETI (TFETI) method. The latter is based on tearing the bodies into “floating” subdomains, discretization by finite elements, and solving the resulting quadratic programming problem by augmented Lagrangians. The presented numerical experiments demonstrate our method’s power and the importance of the proper modelling of 3D frictional contact problems. The state problem solution and the sensitivity analysis process were implemented in parallel.