On the use of an advanced Kirchhoff rod model to study mooring lines

Cette étude présente et valide un modèle avancé de poutre de Kirchhoff, enrichi d'une fonction barrière pour le contact avec le fond marin et de charges hydrodynamiques, démontrant sa capacité à prédire avec précision le comportement statique et dynamique des lignes d'ancrage d'éoliennes flottantes, notamment en révélant des régimes de transition entre traînée et masse ajoutée ainsi qu'un couplage fort entre les dynamiques axiales et de flexion sous sollicitations tangentielles.

L'essentiel

🌊 Le Grand Défi : Tenir les Éoliennes en Mer

Imaginez une éolienne géante flottant au milieu de l'océan. Elle ne repose pas sur le fond, elle flotte ! Pour qu'elle ne parte pas à la dérive avec le vent et les courants, on l'attache au fond de la mer avec de grosses chaînes ou des câbles, un peu comme un chien en laisse. C'est ce qu'on appelle un système d'ancrage.

Le problème ? L'océan est imprévisible. Les vagues poussent, le vent tire, et le fond de la mer est dur. Si ces câbles cassent ou bougent trop, l'éolienne tombe. Les ingénieurs doivent donc simuler comment ces câbles se comportent avant même de les construire.

🧵 L'Analogie du "Fil de Couture Magique"

Jusqu'à présent, pour simuler ces câbles, les ingénieurs utilisaient des modèles un peu "brouillons". Ils imaginaient le câble comme une série de perles enfilées sur un fil (un modèle à masses ponctuelles). C'est rapide à calculer, mais ça manque de précision quand le câble se tord ou se plie de manière complexe.

Dans cet article, les chercheurs (Bruno, Hoa, et leurs collègues) proposent une nouvelle méthode. Ils utilisent ce qu'on appelle un modèle de "bâton de Kirchhoff".

• L'analogie : Imaginez que vous tenez un fil de couture très fin et flexible. Si vous le laissez tomber, il forme une courbe naturelle (une chaînette). Si vous le secouez, il vibre. Ce modèle mathématique traite le câble non pas comme une suite de perles, mais comme un fil continu, lisse et élastique, capable de se courber, de s'étirer et de réagir à la pression de l'eau avec une grande finesse. C'est comme passer d'un dessin fait avec des points à un dessin au trait continu et fluide.

🚧 Le Mur Invisible (Le Fond de la Mer)

Un défi majeur est que le câble touche souvent le fond de la mer. Il ne doit pas passer à travers le sable, mais glisser dessus.

• La métaphore : Imaginez que vous essayez de faire glisser un tapis sur un sol en bois. Si vous poussez trop fort, il reste coincé. Les chercheurs ont ajouté une "fonction barrière" à leur modèle. C'est comme un champ de force invisible juste au-dessus du fond de la mer. Plus le câble s'approche du sol, plus cette force invisible le repousse doucement, l'empêchant de traverser le sable, tout en lui permettant de glisser le long. C'est une astuce mathématique pour simuler le contact sans que le câble ne "traverse" le sol par erreur.

🌊 La Danse avec l'Eau

L'eau n'est pas vide. Quand le câble bouge, il pousse l'eau (comme quand vous marchez dans une piscine). L'eau résiste et pousse en retour.

• L'explication simple : Le modèle calcule deux choses principales :
  1. La traînée (Drag) : C'est la résistance de l'eau quand le câble bouge lentement (comme quand on nage doucement).
  2. La masse ajoutée (Added Mass) : C'est l'inertie. Si vous essayez de bouger très vite un objet dans l'eau, il semble plus lourd car vous devez aussi accélérer l'eau autour de lui.
  • La découverte clé : Les chercheurs ont vu que selon la vitesse à laquelle on secoue le câble (la fréquence), la nature de la résistance change. À basse vitesse, c'est la "traînée" qui domine (comme freiner une voiture). À haute vitesse, c'est la "masse ajoutée" qui domine (comme essayer de courir dans l'eau).

🧪 Les Tests : Est-ce que ça marche ?

Pour vérifier si leur "fil magique" est fiable, ils l'ont mis à l'épreuve de trois façons :

1. Le test de la chaînette : Ils ont laissé un câble tomber entre deux points. Le résultat correspondait parfaitement à la forme mathématique classique attendue.
2. Le test du "pouls" : Ils ont secoué l'extrémité du câble (là où il est attaché à l'éolienne) de manière rythmée. Ils ont observé comment les vibrations se propageaient. Ils ont découvert que si on secoue le câble dans le sens de sa longueur, cela crée des mouvements complexes et imprévisibles, un peu comme un fouet qui se détend.
3. Le test du géant : Ils ont comparé leur modèle avec un logiciel célèbre utilisé par l'industrie (OpenFAST) pour une vraie éolienne de 15 mégawatts (la taille d'un immeuble de 10 étages !).
   • Résultat : Leur modèle a donné des résultats presque identiques à ceux du logiciel industriel, mais avec une précision supérieure sur la forme du câble.

💡 Pourquoi c'est important ?

Ce travail est comme passer d'une carte dessinée à la main à une carte satellite ultra-précise.

• Pour les ingénieurs : Cela permet de mieux comprendre comment les câbles vont réagir aux tempêtes ou aux mouvements violents de l'éolienne.
• Pour l'avenir : Cela aide à construire des éoliennes offshore plus sûres, plus durables et moins chères, car on peut optimiser la longueur et la solidité des câbles sans avoir à en construire trop ou trop peu.

En résumé, ces chercheurs ont créé un simulateur ultra-réaliste pour les câbles d'ancrage, capable de gérer le contact avec le fond de la mer et les forces complexes de l'océan, offrant ainsi une nouvelle vision pour l'énergie éolienne de demain.

Résumé technique

1. Problématique et Contexte

Le développement des éoliennes flottantes offshore (FOWT) impose des défis majeurs dans la conception et l'analyse des systèmes d'amarrage. Ces systèmes sont soumis à des charges environnementales complexes (vent, vagues, courants) et à des interactions avec le fond marin.
Les méthodes d'analyse actuelles présentent des limites :

• Approches statiques/quasi-statiques : Souvent basées sur des équations de caténaire élastique ou des modèles à paramètres concentrés (LPM), elles négligent les forces hydrodynamiques et inertielles, conduisant à une sous-estimation des charges.
• Analyses dynamiques : Bien que les modèles LPM (comme MoorDyn dans OpenFAST) soient efficaces, ils peuvent manquer de précision pour capturer les phénomènes non linéaires complexes, les interactions couplées (axial-flexion) et les effets de raideur géométrique.
• Modèles distribués (DPM) : Les modèles de poutres distribuées sont plus précis mais souvent coûteux en calcul ou complexes à mettre en œuvre, notamment pour gérer le contact avec le fond marin et les forces fluides non conservatives.

L'objectif de cet article est de développer et valider un modèle numérique avancé, nommé ARMoor (Advanced Rod Model for Mooring Lines), capable de simuler avec précision le comportement statique et dynamique des lignes d'amarrage, y compris les interactions avec le fond marin et le fluide environnant.

2. Méthodologie et Modélisation

Le modèle proposé est une extension du modèle de poutre de Kirchhoff non linéaire (sans cisaillement et sans torsion), initialement présenté dans [35], enrichi de deux fonctionnalités clés :

A. Formulation de la Poutre

• Hypothèses : La poutre est transversalement isotrope, initialement droite, et ne subit ni cisaillement ni torsion.
• Équations du mouvement : Déduites du principe de Hamilton, elles incluent l'énergie cinétique et l'énergie potentielle (élastique).
• Discrétisation spatiale : Utilisation de l'Analyse Isogéométrique (IGA) avec des fonctions B-splines, permettant une continuité élevée et une approximation précise de la géométrie sans nécessiter de variables de directeur supplémentaires.
• Intégration temporelle : Un schéma implicite d'ordre 2 (combinaison des règles du point milieu et trapézoïdale) est utilisé pour préserver le moment linéaire et angulaire, et approximer la conservation de l'énergie.

B. Interaction avec le Fond Marin (Fonction Barrière)

Pour modéliser le contact unilatéral avec le fond marin, les auteurs intègrent une fonction barrière pénalisée dans le Lagrangien du système.

• Une contrainte unilatérale $C(\phi) > 0$ est imposée pour empêcher la poutre de pénétrer le fond.
• Différentes fonctions barrières (logarithmique, réciproque, quadratique, etc.) sont étudiées.
• Une étude paramétrique sur le paramètre de pénalité $\mu$ révèle que la fonction barrière réciproque est la plus robuste pour les simulations d'amarrage, permettant de trouver l'équilibre statique via une relaxation dynamique (dissipation d'énergie par le fluide) plutôt que par une résolution statique directe instable.

C. Interaction Fluide-Structure

Un modèle de fluide simplifié mais complet est intégré pour calculer les forces non conservatives par unité de longueur :

• Masse ajoutée : Force d'inertie due à l'accélération relative du fluide.
• Traînée tangentielle et normale : Forces de frottement dépendant du carré de la vitesse relative.
• Poussée d'Archimède : Force hydrostatique verticale.
• Le modèle prend en compte des profils de courant dépendant de la profondeur.

3. Contributions Clés

1. Développement d'ARMoor : Une formulation novatrice combinant la poutre de Kirchhoff non linéaire, l'IGA, une fonction barrière pour le contact sol/poutre et un modèle hydrodynamique complet.
2. Analyse des fonctions barrière : Une étude comparative approfondie des différentes fonctions de pénalité pour la simulation du contact avec le fond marin, identifiant les conditions de stabilité numérique et l'impact du paramètre de pénalité.
3. Découverte de couplages dynamiques : Mise en évidence d'un couplage fort entre la dynamique axiale et la flexion sous des forces tangentes au point d'amarrage (fairlead), contrairement aux forces normales où les dynamiques restent découplées.
4. Transition de régime hydrodynamique : Identification d'une transition du régime dominé par la traînée (basses fréquences) au régime dominé par la masse ajoutée (hautes fréquences) sous chargement pulsatoire.

4. Résultats Numériques et Validation

Trois cas d'étude principaux ont été utilisés pour valider le modèle :

Cas 1 : Comparaison avec la caténaire élastique et contact sol

• Configuration statique : Comparaison avec la solution analytique de la caténaire élastique 3D. Le modèle ARMoor montre un excellent accord géométrique et sur les réactions aux appuis. La prise en compte de la rigidité en flexion rend la solution légèrement plus rigide que la caténaire pure.
• Contact avec le fond : Validation contre des solutions de référence ([6]) pour une ligne reposant sur le fond. Les positions de touchdown et les tensions au point d'amarrage correspondent avec une erreur inférieure à 4 %, confirmant la capacité du modèle à gérer le contact complexe.

Cas 2 : Forces pulsatoires (Dynamique)

• ROD-C (Ligne suspendue) : Validation de la conservation de l'énergie sous chargement harmonique.
• ROD-D (Ligne d'amarrage) : Application de forces pulsatoires tangentes et normales au point d'amarrage.
  • Forces normales : Le comportement suit la théorie des oscillateurs harmoniques amortis (amplitude décroissante avec la fréquence). Transition claire d'un régime dominé par la traînée à basse fréquence vers un régime dominé par la masse ajoutée à haute fréquence.
  • Forces tangentes : Mise en évidence d'un couplage fort entre les modes axiaux et de flexion. L'analyse des portraits de phase et des cartes de Poincaré suggère des comportements non linéaires complexes, voire une transition vers le chaos (mouvement périodique d'ordre supérieur), absents dans les modèles linéaires classiques.

Cas 3 : Plateforme marine UMaine VolturnUS-S (Éolienne 15 MW)

• Comparaison directe avec OpenFAST (utilisant le modèle LPM MoorDyn) pour une éolienne flottante sous chargement de vent par paliers.
• Résultats : ARMoor reproduit avec une grande précision la configuration statique et la distribution de tension (différence maximale de ~0,73 % sur la position du fairlead).
• Dynamique : Sous chargement dynamique, les positions des fairleads et les tensions prédites par ARMoor concordent très bien avec OpenFAST (écarts moyens de l'ordre de 10-40 cm sur une longueur de 850 m). Les écarts mineurs sont attribués à la nature couplée aéro-hydro-élastique d'OpenFAST (mouvements de la plateforme) contre l'approche forcée d'ARMoor.

5. Limitations et Perspectives

• Limitations actuelles : Le modèle suppose une section transversale circulaire et un comportement élastique linéaire. Il ne gère pas encore les matériaux composites ou les comportements viscoplastiques. Le coût computationnel est supérieur à celui des modèles LPM (OpenFAST/OrcaFlex).
• Perspectives : Intégration de la viscoplasticité pour mieux modéliser les câbles réels et réduction du coût de calcul pour des simulations couplées temps réel.

6. Signification et Conclusion

Cet article démontre que l'utilisation d'un modèle de poutre de Kirchhoff non linéaire avancé, enrichi par des fonctions barrière et des modèles hydrodynamiques, offre une alternative robuste et précise aux méthodes traditionnelles pour l'analyse des lignes d'amarrage offshore.

La principale valeur ajoutée réside dans la capacité du modèle à capturer des phénomènes non linéaires complexes (couplage axial-flexion, transitions de régime hydrodynamique, comportements quasi-chaotiques sous chargement tangent) que les modèles simplifiés négligent. Bien que plus coûteux en calcul, ARMoor constitue un outil prometteur pour la conception de systèmes d'amarrage dans des conditions environnementales extrêmes et pour l'analyse de la fatigue, comblant ainsi un vide entre les modèles quasi-statiques rapides et les simulations CFD coûteuses.