Interaction sol-structure non-linéaire en analyse sismique

The present work addresses a computational methodology to solve dynamic problems coupling time and Laplace domain discretizations within a domain decomposition approach. In particular, the proposed methodology aims at meeting the industrial need of performing more accurate seismic risk assessments by accounting for three-dimensional dynamic soil-structure interaction (DSSI) in nonlinear analysis. Two subdomains are considered in this problem. On the one hand, the linear and unbounded domain of soil which is modelled by an impedance operator computed in the Laplace domain using a Boundary Element (BE) method; and, on the other hand, the superstructure which refers not only to the structure and its foundations but also to a region of soil that possibly exhibits nonlinear behaviour. The latter subdomain is formulated in the time domain and discretized using a Finite Element (FE) method. In this framework, the DSSI forces are expressed as a time convolution integral whose kernel is the inverse Laplace transform of the soil impedance matrix. In order to evaluate this convolution in the time domain by means of the soil impedance matrix (available in the Laplace domain), a Convolution Quadrature-based approach called the Hybrid Laplace-Time domain Approach (HLTA), is thus introduced. Its numerical stability when coupled to Newmark time integration schemes is subsequently investigated through several numerical examples of DSSI applications in linear and nonlinear analyses. The HLTA is finally tested on a more complex numerical model, closer to that of an industrial seismic application, and good results are obtained when compared to the reference solutions.Ce travail détaille une approche de calcul pour la résolution de problèmes dynamiques qui combinent des discrétisations en temps et dans le domaine de Laplace reposant sur une technique de sous-structuration. En particulier, la méthode développée cherche à remplir le besoin industriel de réaliser des calculs dynamiques tridimensionnels pour le risque sismique en prenant en compte des effets non-linéaires d'interaction sol-structure (ISS). Deux sous-domaines sont considérés dans ce problème. D'une part, le domaine de sol linéaire et non-borné qui est modélisé par une impédance de bord discrétisée dans le domaine de Laplace au moyen d'une méthode d'éléments de frontière ; et, de l'autre part, la superstructure qui fait référence pas seulement à la structure et sa fondation mais aussi, éventuellement, à une partie du sol présentant un comportement non-linéaire. Ce dernier sous-domaine est formulé dans le domaine temporel et discrétisé avec la méthode des éléments finis (FE). Dans ce cadre, les forces liées à l'ISS s'écrivent sous la forme d'une intégrale de convolution en temps dont le noyau est la transformée de Laplace inverse de la matrice d'impédance de sol. Pour pouvoir évaluer cette convolution dans le domaine temporel à partir d'une impédance de sol définie dans le domaine de Laplace, une approche basée sur des Quadratures de Convolution (QC) est présentée : la méthode hybride Laplace-Temps (L-T). La stabilité numérique de son couplage avec un schéma d'intégration de type Newmark est ensuite étudiée sur plusieurs modèles d'ISS en dynamique linéaire et non-linéaire. Finalement, la méthode L-T est testée sur un modèle numérique plus complexe, proche d'une application sismique de caractère industriel, et des résultats satisfaisants sont obtenus par rapport aux solutions de référence