🔢 mathematics

A 2.5D NURBS-Trace Infinite-Element Method for Moving-Load Wave Propagation and Soil--Structure Interaction in Semi-Infinite Ground

Cet article présente une méthode d'éléments infinis par traçage de NURBS en 2,5D (NBIEM) qui combine l'analyse isogéométrique pour le champ proche avec des fonctions NURBS à produit tensoriel et des fonctions radiales exponentielles pour le champ lointain afin de modéliser efficacement et précisément la propagation d'ondes sous charge mobile et l'interaction sol-structure dans des milieux géotechniques semi-infinis et hétérogènes sans nécessiter de quadrature radiale ou de frontières artificielles.

Auteurs originaux : Yanhui Zhong, Hao Hong, Bei Zhang, Quansheng Zang, Hussein Rappel, Stéphane P. A. Bordas

Publié 2026-07-15

📖 8 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Yanhui Zhong, Hao Hong, Bei Zhang, Quansheng Zang, Hussein Rappel, Stéphane P. A. Bordas

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez que vous essayez d'écouter un train passant à toute allure dans un champ, mais que ce champ est en réalité un océan de terre infini qui s'étend à l'infini dans toutes les directions. Si vous essayiez de construire un modèle informatique de tout cet océan infini pour voir comment le sol tremble, votre ordinateur risquerait probablement de fondre. C'est comme essayer de compter chaque grain de sable sur une plage juste pour voir comment un seul pas fait onduler le sable.

Ce document présente une nouvelle astuce ingénieuse appelée la méthode des éléments infinis NURBS-Trace 2.5D (ou NBIEM pour plus de simplicité) pour résoudre ce problème d'« océan infini » sans avoir besoin d'un supercalculateur de la taille d'une ville.

Le tour de magie : la 2.5D

Parlons d'abord de la partie « 2.5D ». Habituellement, pour modéliser un train circulant sur une voie, vous avez besoin d'un modèle 3D complet. Mais si le sol et la voie se ressemblent tout au long de la ligne (comme un long tunnel droit ou une route droite), vous n'avez pas besoin de modéliser tout le monde en 3D. Vous n'avez besoin de modéliser qu'une seule tranche du sol (la section transversale) et d'utiliser un « sortilège mathématique » (une transformée de Fourier) pour imaginer comment cette tranche se répète au fur et à mesure que le train passe. Cela transforme un puzzle 3D massif en une série de puzzles 2D plus petits et plus gérables.

Le problème : l'extrémité infinie

Même avec cette tranche 2D, il y a un bémol. Le sol ne s'arrête pas ; il continue à l'infini. Dans les anciens modèles informatiques, les scientifiques coupaient simplement le sol à une certaine distance et faisaient comme s'il s'arrêtait là. Mais c'est comme couper l'extrémité d'un trampoline ; les ondes rebondiraient sur votre faux bord coupé et fausseraient les résultats.

Pour correr cela, ils ont utilisé des « éléments infinis ». Considérez-les comme des extensions spéciales et magiques attachées au bord de votre modèle, qui savent comment laisser les ondes s'échapper dans le vide sans rebondir.

La nouvelle solution : NURBS-Trace

La grande innovation des auteurs est la façon dont ils connectent le « champ proche » (le sol directement sous le train) à ce « champ infini » (le reste du monde).

L'ancienne méthode : Imaginez que le champ proche est fait de briques Lego et que le champ infini est fait de pâte à modeler. Pour les connecter, vous devez écraser la pâte à modeler pour qu'elle s'adapte aux Lego, ce qui crée souvent des vides ou des connexions instables.
La nouvelle méthode (NBIEM) : Les auteurs utilisent une technique appelée Analyse Isogéométrique. Au lieu de Lego et de pâte à modeler, ils utilisent le même matériau lisse et flexible (des courbes mathématiques appelées NURBS) pour le champ proche et le champ infini.

Ils appellent cela une méthode « NURBS-trace ». C'est comme avoir un ruban unique et sans couture qui coule de la voie ferrée vers l'infini. Comme le ruban est du même matériau des deux côtés, ils n'ont pas besoin de les forcer à s'emboîter. La connexion est parfaite, fluide et mathématiquement exacte.

Le raccourci « Sans Quadrature »

Voici la partie vraiment géniale. Habituellement, lorsqu'on calcule comment ces ondes infinies se comportent, les ordinateurs doivent effectuer des millions de petites étapes mathématiques répétitives (appelées quadrature numérique) pour deviner la réponse. C'est lent et cela peut devenir désordonné.

Les auteurs ont trouvé un moyen d'écrire la réponse pour la partie infinie en utilisant une formule à forme fermée. Pensez-y de cette façon : au lieu de compter chaque pas d'un coureur pour atteindre la ligne d'arrivée, ils ont trouvé une équation unique qui indique exactement combien de temps dure la course. Cela signifie que leur ordinateur n'a pas besoin de découper le sol infini en petits morceaux ou de faire des millions de suppositions. Il lui suffit de brancher les chiffres dans la formule, et boum — la réponse est là. Cela rend la méthode incroyablement rapide.

Ce qu'ils ont testé (et ce qu'ils n'ont pas testé)

L'équipe n'a pas seulement imaginé cela ; elle l'a passé au crible lors d'une série de simulations :

Le test de vitesse : Ils ont simulé des trains se déplaçant à différentes vitesses :
- Sub-Rayleigh : Plus lent que la vitesse de vibration naturelle du sol (comme une voiture roulant sur une route cahoteuse).
- Super-shear : Plus rapide que les ondes de cisaillement mais plus lent que les ondes de compression (un peu comme le début d'un bang supersonique).
- Super-compressional : Plus rapide que toutes les ondes du sol (un véritable bang supersonique dans le sol).
- Le résultat : Dans tous ces cas, leur nouvelle méthode correspondait presque exactement aux solutions mathématiques « parfaites » connues, même pour les contraintes et les déformations (combien le sol est compressé ou étiré).
La comparaison : Ils ont comparé leur nouvelle méthode à deux autres méthodes populaires :
- FEM/IEM : La méthode standard, plus ancienne.
- IGA/SBIGA : Une méthode plus récente et plus complexe.
- Le résultat : Leur nouvelle méthode était beaucoup plus rapide que la méthode standard et beaucoup plus rapide que la méthode plus récente et complexe (qui a mis 168 secondes pour résoudre un problème que leur méthode a résolu en 0,023 seconde !). Elle était également plus précise que la méthode standard, surtout à des fréquences plus élevées.
Scénarios réels : Ils ont simulé :
- Sol stratifié : Comme un gâteau avec différentes saveurs de terre.
- Tunnels enterrés : Comme la ligne 1 du métro de Zhengzhou. Ils ont vérifié à quelle profondeur le tunnel doit être pour empêcher les vibrations d'atteindre la surface. Ils ont trouvé qu'une fois le tunnel à environ 12 à 14 mètres de profondeur, les vibrations au-dessus diminuent considérablement, et creuser plus profondément apporte de moins en moins de bénéfices.
- Couches tendres : Ils ont testé ce qui se passe s'il y a une couche de sol molle et spongieuse sous la route. Ils ont trouvé que traiter le sol jusqu'à une profondeur d'environ 1,2 à 1,5 mètre est ce qui aide le plus ; aller plus profond n'aide pas beaucoup plus.

Ce qu'ils excluent explicitement

Le papier est très clair sur ce que cette méthode n'est pas :

Ce n'est pas une méthode qui nécessite de connaître la réponse « parfaite » à l'avance pour ajuster vos paramètres. Ils ont créé un flux de travail « sans référence » qui fonctionne sans avoir besoin d'une solution analytique pour vérifier.
Ce n'est pas une méthode qui repose sur les « Couches Parfaitement Adaptées » (PML), qui sont des couches absorbantes artificielles épaisses nécessitant un réglage minutieux. Leur méthode utilise les éléments infinis directement attachés à la limite.
Ce n'est pas une méthode suggérant que la version « calibrée » (qui utilise des mathématiques supplémentaires complexes pour ajuster l'atténuation) est toujours nécessaire. Ils ont constaté que la version plus simple « all-S » (utilisant un seul type d'information d'onde) fonctionne très bien pour la plupart des besoins d'ingénierie, et que la calibration complexe n'est qu'un bonus pour des cas spécifiques de basse fréquence.

À quel point sont-ils sûrs ?

Les auteurs sont très confiants dans leurs simulations. Ils ont prouvé via leurs propres modèles informatiques que :

La méthode fonctionne pour différentes vitesses de train.
Elle est plus rapide et plus précise que les méthodes existantes.
Elle gère bien les formes complexes (comme les tunnels courbes et les sols stratifiés).

Cependant, ils précisent avec prudence qu'il s'agit de simulations. Ils n'ont pas encore comparé leurs résultats à des mesures réelles provenant d'une véritable voie ferrée (bien qu'ils suggèrent que ce serait une excellente étape suivante). Ils notent également que leurs réglages « calibrés » ont été dérivés d'un benchmark spécifique et pourraient nécessiter des ajustements si les conditions du sol changent radicalement.

L'essentiel

Ce papier présente une nouvelle façon super rapide et fluide de modéliser comment le sol tremble lorsqu'un train ou un camion lourd passe. En utilisant un ruban mathématique sans couture pour connecter le sol près de la route au sol infini loin de là, et en utilisant une « formule de raccourci » pour calculer la partie infinie, ils peuvent résoudre ces problèmes beaucoup plus rapidement qu'auparavant. C'est un outil que les ingénieurs peuvent utiliser pour concevoir de meilleures voies et de meilleurs tunnels, garantissant que les vibrations ne dérangent pas les personnes vivant à proximité, le tout sans avoir besoin d'un ordinateur qui coûte plus cher que le PIB d'un petit pays.

Résumé Technique : Une méthode d'éléments infinis par trace de NURBS en 2,5D pour la propagation d'ondes sous charges mobiles

Énoncé du problème
Les vibrations du sol induites par des charges mobiles (ex. : trains à grande vitesse, trafic lourd) constituent une préoccupation critique en géotechnique des transports. La prédiction de ces vibrations nécessite la modélisation de la propagation des ondes dans des domaines de sol semi-infinis et non bornés. Bien que la réponse soit intrinsèquement tridimensionnelle, les systèmes dont la géométrie et les propriétés matérielles sont invariantes le long de la direction de déplacement (l'axe longitudinal $z$ ) peuvent être analysés efficacement à l'aide d'une formulation 2,5D. Cette approche réduit le problème à une séquence de problèmes de sections transversales 2D dans le domaine nombre d'onde–fréquence via une transformation de Fourier. Cependant, modéliser avec précision le domaine extérieur non borné tout en satisfaisant les conditions de rayonnement et en maintenant la cohérence géométrique avec des interfaces complexes (tunnels courbes, sols stratifiés) reste un défi. Les méthodes existantes, telles que le couplage Élément Finis/Élément Infini (FEM/IEM), les Méthodes d'Éléments de Frontière (BEM) et la Méthode des Éléments Finis de Frontière Échelle (SBFEM), font souvent face à des compromis entre coût computationnel, flexibilité géométrique et nécessité de solutions fondamentales ou de couches absorbantes artificielles.

Méthodologie
L'article propose une méthode d'éléments infinis par trace de NURBS en 2,5D (NBIEM). Cette formulation couple un domaine de champ proche borné avec un domaine extérieur semi-infini en utilisant un cadre isogéométrique unifié.

Réduction 2,5D : Les équations du mouvement pour les milieux viscoélastiques linéaires sont transformées dans le domaine fréquence-nombre d'onde. La charge mobile est traitée comme une excitation harmonique où la fréquence temporelle $\omega$ et le nombre d'onde longitudinal $k$ sont couplés via la vitesse de la charge $c_{load}$ ( $k = (\omega - \omega_0)/c_{load}$ ).
Champ proche borné (IGA) : La section transversale du champ proche est discrétisée par l'Analyse Isogéométrique (IGA). La géométrie et le champ de déplacement sont représentés par les mêmes fonctions de base NURBS, permettant une représentation exacte des frontières courbes et une continuité d'ordre élevé.
Extérieur semi-infini (Éléments infinis par trace NURBS) : Le domaine extérieur est modélisé à l'aide d'éléments infinis construits comme des produits tensoriels de :
- La trace NURBS de la frontière (partagée avec la discrétisation du champ proche).
- Des fonctions radiales admissibles ( $P_q$ ) représentant des ondes sortantes ou évanescentes (ex. : décroissance exponentielle).
- Crucialement, le champ proche et le champ lointain partagent le même espace de trace NURBS discret et les mêmes degrés de liberté de contrôle de la frontière. Cela impose la continuité du déplacement directement sans projection, variables de mortier ou opérations de transfert d'interface.
Intégration en forme fermée : Pour des fonctions radiales exponentielles sélectionnées, les moments radiaux requis pour les matrices de rigidité et de masse du champ lointain sont évalués sous forme fermée. Cela élimine le besoin de coupures radiales finies ou de quadrature radiale numérique, garantissant la satisfaction exacte de la condition de rayonnement au sein de l'espace fonctionnel choisi.
Gestion de la condition de rayonnement : La méthode emploie une réalisation « tout-S » où des fonctions radiales informées par les ondes S sont assignées à tous les segments de la frontière. Les paramètres de décroissance sont déterminés par les propriétés matérielles extérieures et les échelles géométriques. Pour assurer la robustesse sans références analytiques, les auteurs introduisent deux indicateurs de sensibilité : $I_\beta$ (sensibilité aux paramètres de décroissance radiale) et $I_D$ (sensibilité au placement de la frontière).

Contributions clés
L'article présente trois contributions principales :

Formulation : Développement d'une formulation d'éléments infinis par trace de NURBS en 2,5D à trois composantes par la méthode de Galerkin. Elle couple de manière unique le domaine IGA borné et l'extérieur semi-infini par une trace NURBS discrète identique, préservant le couplage complet de l'inertie et de la déformation dépendant de la fréquence longitudinale et du nombre d'onde parmi les trois composantes du déplacement.
Efficacité analytique : Évaluation des contributions de rigidité et de masse du champ lointain via des moments radiaux en forme fermée pour des fonctions radiales exponentielles sélectionnées. Cela évite la troncature radiale finie et l'intégration numérique dans la direction radiale.
Caractérisation systématique : Vérification et évaluation complète de la méthode, incluant :
- La validation par rapport à des solutions de demi-espace en forme fermée dans les régimes sub-Rayleigh, super-cisaillement et super-compression.
- Le développement d'indicateurs de sensibilité sans référence ( $I_\beta$ et $I_D$ ) pour filtrer l'adéquation de la frontière.
- Une comparaison directe avec l'approche classique FEM/IEM et IGA/SBIGA concernant la précision, la fidélité de phase et le coût computationnel.

Résultats

Précision : La NBIEM reproduit les fonctions de réponse en fréquence du déplacement et de la contrainte avec une grande précision, correspondant aux solutions analytiques pour les charges mobiles sub-Rayleigh et super-compressibles. Les résultats de contrainte (champs dépendants de la dérivée) confirment la capacité de la méthode à capturer les gradients, bénéficiant de la continuité plus élevée des bases NURBS.
Comparaison avec les alternatives :
- vs. FEM/IEM : La NBIEM atteint une précision de déplacement et de phase comparable ou supérieure avec nettement moins de degrés de liberté. Elle est environ un ordre de grandeur plus rapide en temps de calcul réel pour des maillages de taille modérée à grande en raison de l'élimination de la construction de la rigidité dynamique par échelle.
- vs. IGA/SBIGA : La NBIEM atteint une précision similaire à l'IGA/SBIGA dans la plage de fréquences intermédiaire à élevée, mais à un coût computationnel sensiblement inférieur. Contrairement à l'IGA/SBIGA, elle ne nécessite pas la résolution d'un système de rigidité dynamique par échelle, évitant ainsi les coûts d'assemblage élevés associés.
Applications d'ingénierie : La méthode modélise avec succès :
- Sol stratifié : Gestion des interfaces de matériaux et des contrastes de rigidité.
- Systèmes voie-plateforme : Analyse des configurations de remblai par rapport de déblai avec des géométries multi-patchs, révélant l'amplification dépendante de la vitesse aux interfaces plus profondes.
- Structures enterrées : Simulation de tunnels de métro avec des profondeurs de couverture variables, démontrant comment l'augmentation de la profondeur d'enfouissement atténue la transmission des vibrations vers le sol au-dessus de la voûte du tunnel.
Sensibilité : La réalisation par défaut « tout-S » s'est avérée robuste. Les indicateurs de sensibilité ont identifié efficacement les cas nécessitant un élargissement de la frontière (ex. : basses fréquences ou vitesses super-cisaillement), garantissant la fiabilité de la solution sans avoir recours à des références analytiques.

Signification et affirmations
L'article affirme que la NBIEM fournit un cadre géométriquement cohérent et computationnellement efficace pour la propagation des ondes 2,5D et l'interaction sol-structure dans des domaines semi-infinis. Sa signification réside dans :

Cohérence géométrique : En utilisant des traces NURBS pour les champs proche et lointain, elle gère naturellement les interfaces courbes complexes et les configurations multi-patchs, évitant les erreurs d'approximation géométrique communes aux éléments finis d'ordre bas.
Efficacité computationnelle : Elle offre un équilibre précision-coût favorable, particulièrement pour les analyses répétées dans le domaine fréquentiel, en évitant la construction coûteuse de la rigidité dynamique de la SBFEM et les dépendances aux solutions fondamentales de la BEM.
Robustesse : Les indicateurs de sensibilité proposés permettent un filtrage automatisé et sans référence des réglages de la frontière, rendant la méthode pratique pour les applications d'ingénierie où les solutions analytiques sont indisponibles.

Les auteurs positionnent la NBIEM non pas comme un remplacement de toutes les méthodes existantes, mais comme une alternative distincte et efficace qui comble le fossé entre la flexibilité géométrique de l'analyse isogéométrique et le traitement des domaines non bornés par les éléments infinis, spécifiquement adaptée aux problèmes de charges mobiles en géotechnique.