Les stratifiés composites sont de plus en plus utilisés dans les pièces de structures des aéronefs ce qui fait émerger de nouvelles problématiques comme celle des Impacts de Faible Energie (IFE). En effet, bien qu’ils possèdent des propriétés rapportées à leur masse très intéressantes ces matériaux peuvent être vulnérables aux petits chocs. Or, compte tenu des nombreux paramètres influents lors d’un tel impact (énergie, vitesse, stratification...), les essais actuellement majoritairement privilégiés à l’échelle industrielle sont long et coûteux. Ainsi, l’apport de la simulation numérique pourrait être d’une grande aide pour les constructeurs. La pratique du « virtual testing », en particulier, permettrait d’aller dans cette direction ce qui aurait pour effet de rationaliser les campagnes d’essais et les coûts financier qui en découlent. Cependant, elle peine à être mise en place ici car le temps CPU nécessaire pour la simulation fine des ndommagements induits par les IFE est trop important avec les méthodes actuelles. Partant de ce constat, ce travail a consisté à tirer avantageusement partie de la localisation spatiale et temporelle des délaminages, fissurations matricielles et ruptures de fibres qui peuvent apparaître pendant l’impact pour diminuer le coût de calcul. Ainsi une méthode multiéchelle en espace et en temps a été mise en place. Elle consiste à découper la structure impactée en deux zones. L’une est située autour du point d’impact, elle contient l’ensemble des non-régularités du problème (contact, loi adoucissante, modèle de zone cohésive). Elle est traitée avec le code de dynamique explicite Europlexus. L’autre correspond à la partie complémentaire. Le problème mécanique y est beaucoup plus régulier et il est traité avec le code de dynamique implicite Zset/Zébulon. Un couplage peu intrusif basé sur la méthode GC est donc réalisé entre ces deux codes. Il permet d’utiliser une modélisation adaptée dans chacune des deux régions ce qui permet en particulier d’utiliser des pas de temps différents. Un rapport supérieur à 1000 peut ainsi être obtenu entre celui du code explicite fixé par la condition de stabilité et celui utilisé dans la partie complémentaire. Un gain de temps CPU significatif confirmé par la simulation d’un impact réalisé sur un panneau composite raidi est ainsi obtenu. Il est également montré que la répartition implicite/explicite peut évoluer au cours du calcul. Pour cela un mécanisme de bascule a été mis en place. Il permet ainsi de faire transiter la résolution d’une partie de la structure initialement traitée dans le code Zebulon dans Europlexus. Un gain de temps supplémentaire est alors obtenu grâce à cette méthode sur le même cas d’application. / The composite laminates are increasingly used in aircraft structural parts which lead to new issues such as the Low Energy Impacts (LEI). Indeed, although they have well mechanical properties relative to their mass, small shocks may be very harmfull for laminates. Controlling such situations is essential for manufacturers that why lot of testing campaigns are currently performed. Yet, they are time consuming and expensive considering the many influential parameters (energy, speed, layup...). Numerical simulations of this phenomenon by practicing the so called “virtual testing” process could be really helpfull to rationalize testing campaigns in order to save money. Yet, this practice remain currently hard to do at the industrial scale due to the excessive CPU time required for fine simulation of damages induced by the LEI. Based on this observation, this work has consisted in taking advantage of the spatial and temporal location of delamination, matrix cracking and fiber breakage that can occur during impact in order to reduce the computational cost. Thus, a space and time multiscale method has been put in place. The impacted structure is split into two areas. One is located around the impacted point, it contains all the non-regularities of the problem (contact, softening law, cohesive zone model). This domain is treated with the explicit dynamics code Europlexus. The other one corresponds to the complementary part. The mechanical problem is much more regular and it is treated with the implicit dynamics code Zset / Zebulon. A low intrusive coupling based on the GC method is carried out between these two codes. It allows to use an adapted model in both regions different time step are in particular used. A time step ratio upper to 1000 can be reach between the one of the explicit code set by the stability condition and the one used in the complementary part. As a results, significant CPU time is saved. This is confirmed by the simulation of a stiffened composite panel impacted. It is also shown that the implicit / explicit allocation can change over the calculation. To do that, a switch mechanism has been established. It thus makes it possible to transit the resolution of a portion of the structure initially solved in the code Zebulon to Europlexus. As a results, further gain is obtained.

Il est prévu d’augmenter la durée de vie des centrales nucléaires actuellement opérationnelles. Le vieillissement des structures en béton armé telles que les tours aéroréfrigérantes doit être évalué, son incidence sur la capacité portante calculée. Dans le cas de fortes dégradations, le renforcement doit être envisagé, afin d’assurer la pérennité de ces tours face aux sollicitations extrêmes telles que les tempêtes de vent et les séismes. Ce travail vise à quantifier les effets néfastes que peut générer la réduction de section des aciers induite par la corrosion, en particulier sur la capacité portante des tours dans des conditions de sollicitations extrêmes monotones ou cycliques de types vent et séisme. Ces sollicitations sont certainement les plus sévères, entrainant la structure dans le domaine non linéaire, elles sont susceptibles d’induire des endommagements de type fissuration qui dans le cas de sollicitation cycliques peuvent s’avérer néfastes. Des modélisations numériques sont proposées pour déterminer la réponse quasi-statique ou dynamique de la structure, en tenant compte des apparitions de fissures dans le béton et de leur évolution via des lois de comportement appropriées du matériau béton, ainsi que la plastification des aciers. Dans le cas d’une sollicitation sismique, dans le but de comparer les approches de modélisation du séisme et d’évaluer la robustesse des résultats, les réponses dynamiques sont évaluées par trois méthodes différentes de calcul : l’approche dynamique temporelle non linéaire, la méthode spectrale et la méthode modale temporelle. Des études paramétriques portant sur l’amortissement, les combinaisons de charges et les configurations structurales, sont aussi menées. Dans le cas d’une sollicitation de type vent, la technique de renforcement à l’aide de matériaux composite, tel que le tissu de fibres de carbone (TFC) est modélisée. Le comportement de la structure endommagée présentant un taux de corrosion avancée, est évalué dans le régime pré- et post-fissuration, comparativement à la structure intègre. La perte de capacité portante est quantifiée, un renforcement permettant de restaurer l'intégrité et donc d’augmenter la durée de vie de la structure est proposé. / The authorities have planned to increase the lifetime of currently operating nuclear power plants. The ageing of reinforced concrete structures such as cooling towers should be evaluated and its impact on the bearing capacity calculated. In the case of significant damage, the strengthening must be considered to ensure the sustainability of these towers facing the risk of storms and earthquakes becoming more and more frequent. This work aims to quantify the adverse effects that can generate concrete cracks and rebar section loss induced by corrosion, especially on the bearing capacity of nuclear power plant cooling towers under monotonic or cyclic extreme load conditions (wind and earthquake). These loads are certainly the most severe, since they take the structure into the nonlinear domain and can induce or amplify cracking damage. Numerical simulations are proposed to determine the quasi-static or dynamic response of the structure, taking into account appearance of concrete cracks and their evolution via an appropriate material concrete law and rebar's yielding. In the case of a seismic load, the responses are evaluated by three different methods; the nonlinear response history analysis (NLRHA), the response spectrum analysis and the modal response history analysis (MRHA) in order to compare the earthquake modeling approaches and to evaluate the robustness of the results. Parametric studies on damping, load combinations and structural configurations, are also performed. In the case of a wind load, the strengthening technique using composite materials, such as carbon fiber reinforced plastic (CFRP) is modeled. The behavior of the damaged structure with an advanced corrosion rate is estimated in the pre- and post-cracking regime, compared to the undamaged structure. The drop of bearing capacity is quantified, a reinforcement designed is proposed to restore the integrity and thus increase the lifetime of the structure.