Modèle multi-échelle de la fatigue des lignes d’ancrage câblées pour l’éolien offshore flottant

The main function of mooring systems of floating offshore wind turbines is to ensure station keeping. The mooring lines can be composed of chains, wire ropes, synthetic ropes, or even a combination of them. In this thesis we focus on wire ropes, whose advantage over chain is to sustain high tension at a lower weight. Their design must consider the successive tension and bending loading induced by the floater movement for various wind and waves conditions.The thesis purpose is to develop a new numerical model, dedicated to the prediction of fatigue damage in mooring wire ropes of a floating wind turbine. In particular it has to simulate the relative movements between the wires when the rope is bent. Results from free-bending fatigue tests in the literature show the importance of these effects, since the first rupture is localized near the neutral plane, where fretting is more important. This phenomenon affecting the fatigue life is not considered by fatigue criteria of current offshore standards, which are related to tension-tension loading.It is worth noting that the use of a detailed model of wire rope in a fatigue design procedure represents a real challenge. The high number of contact interactions to be modeled, which are several thousands per meter of rope, and the large amount of loading cases make this type of computations extremely time-consuming.The loading used in the developed local model of wire rope is obtained from global computations performed with a dedicated multiphysics software (Deeplines). This software allows to simulate the environmental conditions (wind, waves, current) applied on the whole structural system.Some preliminary computations showed that the nonlinear bending behavior of the wire rope, linked to the wire contact interactions, does not significantly affect the output of the global model. This observation justifies the use of a top-down scheme, with a prior computation of the global scale.The global scale tension and curvature are then uniformly imposed on the central wire of the local model. The continuity of the rope is represented by periodic conditions which link the end sections to points within the model, at the same circumferential locations. The wires are modeled by beam elements. The outputs at the local scale are the stress resultants on the wires, and the contact forces and relative displacements at contact locations.Small sliding between the wires has been observed from first numerical analysis, for a representative loading case. Therefore, in order to reduce the computational cost of the wire rope model, a new node-to-node contact element has been developed, dedicated to the modeling of contact between non-parallel beams with circular cross section. It assumes fixed contact pairing and finite rotations. Numerical benchmarks and experimental tests on wire ropes show the improvement with results closer to a reference surface-to-surface model, when compared to standard algorithm for the simulation of contact between beams. Moreover, the new model reduces significantly the CPU cost and is also more robust, which is crucial for fatigue life estimates.The outputs of the local scale model are then used to obtain the complete 3D stress state by means of analytical solutions of contact between solids with cylindrical shape. Finally, a multiaxial fatigue criterion is applied in order to assess the safety of the system.La fonction principale des systèmes d'ancrage des éoliennes offshore flottantes est de limiter les mouvements du support. Les lignes d'ancrage qui les composent sont typiquement constituées de chaînes, de câbles aciers, de câbles synthétiques ou d'une combinaison de ces composants.Dans cette thèse, on se concentre sur les câbles en acier qui permettent de réduire le poids et d'augmenter la résistance en tension par rapport aux chaînes. Leur dimensionnement dépend des chargements en tension et flexion, liés aux mouvements du flotteur sous l'action de la mer et du vent.L'objectif de la thèse est le développement d'un nouveau modèle numérique pour prédire la durée de vie en fatigue des câbles d'ancrage d'une éolienne offshore flottante. Il doit notamment simuler les glissements relatifs entre les fils au cours d'une flexion du câble. Des résultats d'essais de tension-flexion de la littérature ont en effet montré que la première rupture est localisée près du plan neutre de flexion, où ces déplacements relatifs sont les plus grands. Cet effet majeur sur la durée de vie du câble n'est pas pris en compte par les lois de fatigue en tension-tension des normes de design offshore actuelles.Il faut aussi remarquer que l'utilisation d'un modèle détaillé de câble dans une démarche de dimensionnement à la fatigue représente un vrai défi. Le nombre élevé d'interactions de contact à modéliser, de l'ordre de plusieurs milliers par mètre de câble, et le grand nombre de cas de chargement rendent ce type de calculs très coûteux.Les chargements qui sont utilisés dans le modèle local de câble sont issus de calculs globaux réalisés à l'aide d'un logiciel multiphysique (Deeplines). Ce logiciel permet de simuler les conditions environnementales (vent, houle, courant) appliquées sur l'ensemble de la structure offshore.Nous montrons que le comportement non linéaire en flexion du câble, lié aux interactions de contact entre les fils, n'influence pas significativement les résultats du modèle global. Cette observation justifie une démarche de type descendante, les calculs globaux pouvant être réalisés en première étape. Les évolutions temporelles des tensions et courbures globales sont appliquées uniformément sur le fil central du modèle local du câble. La continuité du câble est représentée par des conditions de périodicité reliant les sections de bord à des points internes du modèle situés sur la même position circonférentielle. Les fils sont modélisés par des éléments poutres. On obtient les contraintes généralisées sur les fils, les forces de contact et les glissements relatifs. Des premières analyses ont montré que les déplacements relatifs entre les fils restent petits dans notre cadre d'application. Afin de réduire le coût calcul, nous avons développé un nouvel élément de contact entre poutres non parallèles, avec un appariement fixe de contact, dans l'hypothèse de petits glissements mais en grands déplacements et grandes rotations. Des tests numériques montrent l'amélioration obtenue, avec un résultat plus proche d'un modèle de référence qui considère un contact surfacique. De plus, le nouveau modèle réduit significativement le coût calcul et se montre plus robuste en convergence, ce qui s'avère crucial pour un calcul de fatigue. Les sorties du modèle local sont ensuite utilisées pour prédire un état de contrainte 3D, en exploitant des solutions analytiques de contact entre corps cylindriques. Finalement, un critère de fatigue multiaxial de la littérature est appliqué pour évaluer le risque en dommage