Dual finite-element analysis using second-order cone programming for structures including contact

Chadi El Boustani (1) , Jérémy Bleyer (1) , Mathieu Arquier (2) , Mohammed-Khalil Ferradi (2) , Karam Sab (1)

1 NAVIER UMR 8205 - Laboratoire Navier
2 Strains

Computation of elastic structures in contact is performed by means of a dual analysis combining displacement-based and equilibrium-based finite elements. Contact conditions are formulated in the framework of second-order cone programmaing (SOCP) and an efficient interior point method (IPM) algorithm is presented for the resolution of the associated optimization problems. The dual approach allows the user to assess the quality of convergence and to efficiently calculate a discretization error estimator which includes a contact error term. An efficient remeshing scheme, based on the local contributions of the elements to the global error, can then be used to efficiently improve the solution accuracy. The whole process is illustrated on some examples and applied to a typical steel assembly. Its efficiency, in particular concerning the IPM solver, is demonstrated in comparison with the industrial finite element code Abaqus.

Mots clés: contact elastostatics; equilibrium finite elements; error estimator; second order cone programing; interior point method

Domaines: Physique [physics]; Mécanique [physics]; Matériaux et structures en mécanique [physics.class-ph]

Chadi El Boustani, Jérémy Bleyer, Mathieu Arquier, Mohammed-Khalil Ferradi, Karam Sab. Dual finite-element analysis using second-order cone programming for structures including contact. Engineering Structures, 2020, 208, 20 p. ⟨10.1016/j.engstruct.2019.109892⟩. ⟨hal-02428540⟩