Homogenization of rod-like metamaterials as a special Cosserat rod

Cet article présente un schéma d'homogénéisation variationnel fondé sur la théorie géométriquement exacte des poutres spéciales de Cosserat pour dériver la réponse constitutive non linéaire et la rigidité de métamatériaux de type poutre assemblés périodiquement, validé par des exemples numériques allant de réseaux simples à des structures complexes auxétiques et de muscles artificiels.

L'essentiel

Imaginez que vous avez un long tube flexible, non pas fait d'un seul matériau solide, mais d'un motif complexe et répétitif de minuscules bâtonnets connectés entre eux, comme une échelle microscopique ou une clôture en maillons de chaîne enroulée en cylindre. C'est ce que les auteurs appellent un « matériau métamorphe de type tige ».

Le problème qu'ils abordent est le suivant : si vous voulez savoir comment ce long tube entier se plie, s'étire ou se tord, vous ne pouvez pas simplement regarder un seul minuscule bâtonnet. Vous devez examiner comment l'ensemble du réseau de milliers de bâtonnets interagit. Simuler chaque bâtonnet individuel pour un long tube revient à essayer de compter chaque grain de sable sur une plage pour comprendre comment la plage bouge sous l'effet du vent : cela demande trop de puissance informatique.

Les auteurs proposent un raccourci ingénieux, une « recette » pour prédire le comportement de tout le tube en n'étudiant qu'un tout petit morceau représentatif de celui-ci. Voici comment ils procèdent, expliqué avec des analogies simples :

1. Le « Zoom Magique » (Homogénéisation)

Imaginez le matériau métamorphe comme un immense motif de papier peint répétitif. Au lieu d'analyser tout le mur, vous regardez simplement un seul carré du papier peint (appelé RVE, ou Élément de Volume Représentatif).

L'astuce des auteurs consiste à supposer que si vous étirez ou tordez tout le long tube, ce petit carré s'étire ou se tord également, mais d'une manière très spécifique et hélicoïdale. Ils appellent cela une déformation « hélicoïdale ». Imaginez que vous tirez sur un ressort : les spires ne s'allongent pas seulement ; elles tournent également légèrement. Les auteurs ont réalisé qu'en forçant ce petit carré à imiter exactement ce mouvement en spirale, ils pouvaient déterminer comment tout le long tube réagirait sans avoir à simuler l'ensemble.

2. Les Tiges « Parfaitement Flexibles »

La plupart des modèles informatiques traitent les tiges comme rigides et immuables, comme une règle en acier. Mais dans la réalité, surtout avec ces matériaux métamorphes microscopiques, les tiges peuvent se plier, s'étirer et se cisailer (glisser latéralement) simultanément, même lorsque la déformation est énorme.

Les auteurs utilisent un modèle mathématique spécial appelé « tige de Cosserat spéciale ».

• Analogie : Imaginez un morceau de spaghetti cuit. Il peut se plier, s'étirer un peu et se tordre. Maintenant, imaginez que ce spaghetti est fait d'un matériau capable de faire toutes ces choses parfaitement et précisément, même si vous le pliez en cercle ou l'étirez jusqu'à doubler sa longueur. C'est ce que fait leur modèle. Il ne se contente pas d'approximer ; il capture la géométrie exacte du pliage et de la torsion.

3. Les Règles de la « Piste de Danse » (Conditions aux Limites)

Pour faire en sorte que le petit carré se comporte comme s'il faisait partie d'un tube géant et répétitif, les auteurs ont dû inventer un ensemble de règles régissant la façon dont les bords de ce carré communiquent entre eux.

• Le Problème : Si vous coupez un morceau d'un escalier en colimaçon, le bord supérieur ne s'aligne pas parfaitement avec le bord inférieur.
• La Solution : Ils ont créé une « condition aux limites hélicoïdale ». Imaginez que le côté gauche de votre petit carré tient la main du côté droit, mais que le côté droit est légèrement tourné et décalé, tout comme les marches d'un escalier en colimaçon.
• L'Innovation : Les méthodes précédentes ne pouvaient gérer que des mouvements petits et doux. La nouvelle règle des auteurs fonctionne même si le tube est tordu en forme de bretzel ou étiré jusqu'à devenir fin comme un fil. Elle est « géométriquement exacte », ce qui signifie qu'elle ne perd jamais en précision, quelle que soit l'extravagance de la forme.

4. Les « Articulations » et la « Colle »

À l'intérieur de ce petit carré, les tiges sont connectées par des articulations.

• Articulations Rigides : Certaines articulations agissent comme une colle ultra-forte ; les tiges ne peuvent pas bouger les unes par rapport aux autres au point de connexion.
• Les Mathématiques : Les auteurs ont mis en place un système où l'ordinateur calcule la meilleure position de chaque tige dans ce petit carré, en veillant à ce que les articulations restent connectées et que les règles de « l'escalier en colimaçon » soient respectées, tout en utilisant le minimum d'énergie possible.

5. Ce Qu'ils Ont Découvert (Les Résultats)

Une fois qu'ils ont résolu les mathématiques pour le petit morceau, ils ont pu prédire comment tout le tube se comporterait. Ils ont testé cela avec différentes formes :

• Les Formes en Croix et en Carré : Ils ont commencé par des formes simples (comme un signe plus ou un carré fait de tiges) pour prouver que leurs mathématiques fonctionnaient. Ils ont constaté que si les minuscules tiges sont épaisses et courtes, il est très important de savoir si elles peuvent s'étirer ou se cisailer. Si elles sont très fines et longues, les anciennes mathématiques, plus simples, fonctionnent bien.
• Les Tiges Hélicoïdales (Ressorts) : Ils ont examiné un carré composé de tiges déjà courbées comme des ressorts (hélices).
  • L'Étirement en Forme de « J » : Lorsqu'ils ont tiré sur ce matériau, il était d'abord mou (comme un ressort qui se déroule) mais devenait très rigide en se redressant. Cela crée une courbe en forme de « J ». C'est exactement ainsi que se comportent les tissus biologiques (comme les muscles), raison pour laquelle les auteurs mentionnent qu'ils pourraient être utilisés pour des muscles artificiels.
  • L'Adoucissement au Pliage : Lorsqu'ils l'ont plié, le matériau est devenu plus mou à mesure qu'ils le pliaient davantage. Cela s'est produit parce que la tige-ressort de connexion a commencé à se tordre hors du plan, agissant comme une charnière.
• Le Tube Auxétique : Ils ont modélisé un tube creux qui s'élargit lorsque vous le tirez (comme un nid d'abeilles).
  • Ils ont montré qu'en modifiant l'angle des tiges, vous pouvez régler le tube pour qu'il soit très flexible de côté à côté (bon pour le pliage) mais très rigide contre l'écrasement (bon pour maintenir les vaisseaux sanguins ouverts).
  • Ils ont noté que ces structures peuvent être réglées pour éviter le « raccourcissement » (devenir plus court lors de l'expansion), qui est un problème courant dans les stents cardiovasculaires (tubes en maille utilisés pour maintenir les artères ouvertes).

Résumé

Les auteurs ont construit un « traducteur universel » pour les matériaux métamorphes. Ils ont créé une méthode qui prend un réseau complexe et tridimensionnel de minuscules tiges et le traduit en une description mathématique simple et lisse d'une seule tige. Cela permet aux ingénieurs de concevoir des matériaux complexes et flexibles pour des choses comme des bras robotiques, des muscles artificiels et des stents médicaux en ajustant les motifs internes microscopiques, sachant exactement comment le produit final se pliera et s'étirera, sans avoir besoin de faire tourner une simulation sur un superordinateur pour chaque modification de conception.

Résumé technique

Résumé technique : Homogénéisation de métamatériaux de type tige comme tige de Cosserat spéciale

Énoncé du problème
Les métamatériaux de type tige sont des structures formées par l'assemblage périodique d'unités microstructurales (généralement des réseaux de tiges) dans une seule direction, résultant en une géométrie macroscopique significativement plus longue dans une dimension que dans les deux autres. Ces structures sont utilisées dans des applications allant des stents cardiovasculaires et des muscles artificiels aux actionneurs robotiques. Alors que la simulation numérique directe à l'échelle complète de telles structures hétérogènes est prohibitive en termes de calcul, les techniques d'homogénéisation standard rencontrent souvent des limitations. Plus précisément, les schémas d'homogénéisation du premier ordre sont insuffisants pour capturer les grandes déformations impliquant des gradients de déformation significatifs (tels que la flexion et la torsion) typiques des structures de type tige. De plus, les approches d'homogénéisation tridimensionnelles standard ne parviennent pas à prendre en compte les effets de taille lorsque la séparation des échelles de longueur ne s'applique pas dans les directions transversales. Il existe un besoin d'un cadre capable de modéliser avec précision la réponse constitutive non linéaire de ces métamatériaux sous de grandes déformations arbitraires, en tenant compte du cisaillement, de l'allongement, de la flexion et de la torsion finis aux échelles micro et macro.

Méthodologie
Les auteurs proposent un schéma d'homogénéisation computationnelle qui modélise les métamatériaux de type tige comme une « tige de Cosserat spéciale » à l'échelle macroscopique, tandis que l'unité microstructurale (Volume Élémentaire Représentatif ou VER) est modélisée comme un réseau de tiges de Cosserat géométriquement exactes.

1. Cinématique géométriquement exacte : La tige effective macroscopique et les tiges constitutives microscopiques sont toutes deux modélisées en utilisant la théorie géométriquement exacte de la tige de Cosserat spéciale. Cette formulation prend en compte la flexion, la torsion, le cisaillement et l'allongement axial finis, évitant ainsi les hypothèses de petites déformations des théories classiques des poutres.
2. Extension de la règle d'Hélice Cauchy-Born (HCB) : La méthode étend la règle d'Hélice Cauchy-Born. Elle suppose que le métamatériau est soumis à un champ de déformation uniforme le long de sa longueur d'arc à l'échelle macroscopique. Cette hypothèse permet de réduire le problème structurel complet à celui d'un seul VER.
3. Conditions aux limites périodiques hélicoïdales : Une innovation clé réside dans la dérivation d'une application géométriquement exacte et non linéaire pour appliquer des conditions aux limites périodiques hélicoïdales au VER. Cette application relie les degrés de liberté de translation et de rotation de la frontière droite du VER à sa frontière gauche sur la base des paramètres de déformation macroscopique ($v_0, k_0$). Contrairement aux applications de transition linéaires précédentes, cette formulation reste valable pour des déformations macroscopiques arbitrairement grandes.
4. Minimisation d'énergie contrainte : Le problème du VER est formulé comme une minimisation contrainte de l'énergie élastique. Les contraintes incluent :
   • Contraintes de joints internes : Imposant une connectivité rigide entre les tiges au sein du VER en utilisant une approche maître-esclave.
   • Contraintes de limites hélicoïdales : Imposant la périodicité dérivée de la règle HCB.
   • Contraintes globales : Restreignant la translation et la rotation du corps rigide du VER pour assurer une solution unique et maintenir des paramètres de déformation macroscopique fixes pendant la minimisation.
5. Discrétisation et résolution : L'énergie du VER est discrétisée à l'aide d'« éléments hélicoïdaux », où chaque segment est supposé subir un champ de déformation uniforme, rendant le segment lui-même hélicoïdal. La forme faible non linéaire résultante est résolue à l'aide d'un code interne sur la plateforme deal.II.
6. Dérivation des propriétés macroscopiques : Des expressions analytiques pour la force de contact interne, le moment et le tenseur de rigidité de la tige homogénéisée sont dérivées en différenciant la densité d'énergie du VER minimisée par rapport aux paramètres de déformation macroscopique.

Contributions clés
L'article identifie trois contributions principales novatrices :

1. Exactitude géométrique à double échelle : Les tiges aux échelles micro et macro sont modélisées comme des tiges de Cosserat spéciales géométriquement exactes, capturant le cisaillement et l'allongement finis en plus de la flexion et de la torsion aux deux échelles.
2. Application macro-vers-micro géométriquement exacte : Une application entièrement non linéaire est proposée pour appliquer des conditions aux limites périodiques hélicoïdales aux degrés de liberté de translation et de rotation. Cette application est valable même pour de grandes déformations macroscopiques, répondant ainsi à une limitation des applications de transition linéaires précédentes.
3. Comportement constitutif microscopique arbitraire : Le cadre permet un comportement constitutif arbitraire des tiges constitutives au sein du VER, allant au-delà des hypothèses d'élasticité linéaire.

Résultats et exemples numériques
Les auteurs valident le schéma en utilisant plusieurs exemples numériques avec des complexités de VER variables :

• VER en croix et carrés : Des VER plans simples (formes en croix et carrées) ont été utilisés pour valider les résultats par rapport à des formules analytiques existantes et à la littérature. L'étude a démontré que l'inclusion du cisaillement et de l'extension dans les tiges microscopiques affecte significativement la rigidité macroscopique, en particulier pour les tiges ayant des rapports d'élancement plus élevés. Les résultats correspondaient aux limites analytiques pour les tiges de Kirchhoff (inextensibles/incompressibles).
• VER de tiges hélicoïdales : Un VER carré plus complexe composé de tiges constitutives hélicoïdales a été analysé. Les résultats ont montré une réponse contrainte-allongement macroscopique en forme de J, passant d'un comportement dominé par la flexion à un comportement dominé par l'étirement. L'étude a mis en évidence comment des paramètres tels que le nombre de tours hélicoïdaux et le pas pouvaient ajuster cette réponse. L'analyse de la flexion a révélé un comportement d'adoucissement déclenché par la flexion hors-plan des tiges hélicoïdales de liaison transversale.
• Structures tubulaires auxétiques : Le schéma a été appliqué à des métamatériaux tubulaires auxétiques (pertinents pour les stents). Le modèle a réussi à capturer la dépendance de la rigidité de flexion macroscopique, du coefficient de Poisson (comportement de raccourcissement) et de la rigidité radiale par rapport à l'angle de conception $\theta_a$. Il a également démontré la capture de l'effet Bazier (ovalisation de la section transversale) sous une grande flexion.
• Validation FE2 : Une comparaison entre le schéma d'homogénéisation proposé (approche FE2) et des simulations à l'échelle complète d'un métamatériau à VER en croix a montré une bonne concordance, validant la précision de la méthode même sous de grands moments de torsion où les modèles constitutifs linéaires échouaient.

Importance et affirmations
L'article affirme fournir un cadre robuste pour l'homogénéisation des métamatériaux de type tige qui surmonte les limitations des schémas d'homogénéisation du premier ordre et linéaires. En utilisant une théorie géométriquement exacte de la tige de Cosserat spéciale aux deux échelles et en dérivant une application de transition macro-vers-micro non linéaire, la méthode prédit avec précision la réponse constitutive non linéaire de ces structures sous de grandes déformations.

Les auteurs soulignent que leur approche est générale et peut être étendue à d'autres types de VER (par exemple, solides poreux, structures d'origami) et peut intégrer des phénomènes microscopiques complexes tels que l'élastoplasticité, la viscoélasticité et le couplage multiphysique (thermo-, chemo- ou électro-élasticité). Ce travail sert de fondation pour la conception de métamatériaux sur mesure avec des propriétés macroscopiques spécifiques, telles que des muscles artificiels ou des stents déployables, en ajustant les paramètres microstructuraux. Les auteurs notent modestement que, bien que le travail actuel suppose que le VER est un réseau de tiges, les conditions aux limites dérivées pourraient être adaptées à d'autres représentations microstructurales.