Sur le comportement singulier dans un problème de l'élastostatique linéaire plane

Informations fondamentales

• ID de l'article: 2409.07954
• Titre: On singular behaviour in a plane linear elastostatics problem
• Auteurs: H. Gimperlein, M. Grinfeld, R. J. Knops, M. Slemrod
• Classification: math.AP (Mathématiques - Équations aux dérivées partielles)
• Date de publication: 25 novembre 2024 (version v2)
• Lien de l'article: https://arxiv.org/abs/2409.07954

Résumé

Cet article construit un champ vectoriel analogue à celui introduit respectivement par Artstein et Dafermos, provenant des vecteurs tangents à une famille de cercles non concentriques à un paramètre, croissante de manière monotone, qui sont tous tangents à l'origine. Ces cercles définissent et remplissent une région en forme de lentille Ω, dont les frontières extérieure et intérieure sont respectivement les cercles maximal et minimal. La cuspidale double à l'origine produit une singularité géométrique, au point de laquelle le champ vectoriel est indéterminé et présente un comportement limite non unique. Par la méthode semi-inverse impliquant la fonction de contrainte d'Airy, on montre que ce champ vectoriel correspond au champ de déplacement d'un problème d'élastostatique plane linéaire compressible non homogène isotrope sur Ω, avec des rotations rigides relatives sur la frontière, pouvant potentiellement causer une rupture ou une déchirure à l'origine. L'article trouve une série de solutions où non seulement les paramètres de Lamé satisfont la condition d'ellipticité forte, mais le comportement limite non unique du déplacement est également préservé.

Contexte et motivation de la recherche

Formulation du problème

1. Problème central: Étudier le comportement singulier du champ de déplacement dans les problèmes d'élastostatique plane avec singularités géométriques. Plus précisément, explorer l'indétermination et le comportement limite non unique du champ de déplacement au point singulier lorsque la frontière subit des rotations rigides relatives sur une région en forme de lentille (avec une cuspidale double à l'origine).
2. Importance du problème:
   • D'un point de vue mathématique, ce problème concerne l'existence, l'unicité et la régularité des solutions d'équations aux dérivées partielles sur des domaines avec singularités géométriques
   • D'un point de vue physique, le comportement singulier au point singulier peut correspondre à des phénomènes importants tels que la rupture, la déchirure ou la défaillance fragile du matériau
   • Cette recherche prolonge les travaux d'Artstein et Dafermos sur les solutions non uniques, fournissant une nouvelle interprétation dans le contexte de la mécanique des milieux continus
3. Limitations des méthodes existantes:
   • Artstein a considéré une variante de ce champ vectoriel dans les systèmes de contrôle, mais il n'existe pas de contrôle stable lisse
   • Dafermos l'a utilisé comme exemple de problème de valeur initiale pour les équations différentielles ordinaires d'évolution, montrant des solutions non uniques non lisses
   • Les travaux antérieurs (référence 3) ont prouvé que ce champ vectoriel constitue les trajectoires lagrangiennes des particules fluides dans un écoulement compressible stable
   • Cependant, le champ vectoriel n'a pas été systématiquement étudié dans le cadre de l'élastostatique
4. Motivation de cet article:
   • Adopter une approche différente d'Artstein et Dafermos pour construire le champ vectoriel (basée sur les tangentes à la famille de cercles)
   • Interpréter le comportement singulier dans le contexte de l'élastostatique plane linéaire
   • Trouver des paramètres de matériau non homogène satisfaisant la condition d'ellipticité forte
   • Approfondir la compréhension de la nature du point singulier par le calcul de la force totale, du moment total et de l'énergie de déformation

Contributions principales

1. Construction géométrique du champ vectoriel: Propose une méthode de construction du champ vectoriel basée sur les tangentes à une famille de cercles non concentriques à un paramètre, méthode qui peut être naturellement généralisée et qui montre clairement l'origine de la singularité géométrique
2. Interprétation en élastostatique: Par la méthode semi-inverse et la fonction de contrainte d'Airy, on prouve que ce champ vectoriel correspond au champ de déplacement d'un corps élastique isotrope linéaire compressible non homogène, en équilibre sous l'action de forces volumiques nulles
3. Condition d'ellipticité forte: Trouve une série de solutions dont les paramètres de Lamé satisfont la condition d'ellipticité forte $\mu(\lambda + 2\mu) > 0$, tout en préservant le comportement limite non unique du déplacement
4. Analyse de la singularité:
   • Prouve l'indétermination du champ vectoriel à l'origine et le comportement limite non unique (Proposition 4.1)
   • Calcule la force totale (nulle) et le moment total (proportionnel à $R^2-1$) (Proposition 6.1)
   • Analyse la singularité de l'énergie de déformation ($O(r^{-1})$ quand $k>0$, sans singularité quand $k=0$)
5. Interprétation physique: Interprète les rotations rigides relatives de la frontière et le comportement singulier à l'origine comme des phénomènes possibles de rupture, de déchirure ou de défaillance fragile

Détails de la méthode

Définition de la tâche

Étudier un problème aux limites d'élastostatique linéaire plane sur une région en forme de lentille $\Omega \subset \mathbb{R}^2$:

• Entrée: Région géométrique (définie par une famille de cercles non concentriques), conditions aux limites (rotations rigides)
• Sortie: Champ de déplacement $\mathbf{u}(x_1, x_2)$, champ de contrainte $\sigma_{\alpha\beta}$, paramètres de Lamé non homogènes $\lambda(x), \mu(x)$
• Contraintes: Équation d'équilibre $\sigma_{\alpha\beta,\beta} = 0$, relations constitutives, condition d'ellipticité forte

Construction géométrique et champ vectoriel

1. Définition de la famille de cercles Considérer une famille de cercles non concentriques à un paramètre: $\psi_c(x) := c^{-2}\{(x_1-c)^2 + x_2^2\} = 1$ où $1 \leq c \leq R$, les centres sont en $(c,0)$, les rayons sont $c$, tous les cercles sont tangents à l'origine $(0,0)$.

2. Champ vectoriel tangent Le vecteur unitaire tangent au point $x$ sur le cercle: $t_1(x_1,x_2) = \frac{x_2}{c}, \quad t_2(x_1,x_2) = \frac{c-x_1}{c}$

3. Champ vectoriel de déplacement En posant $\mathbf{u} = c\mathbf{t}$, on obtient: $u_1(x_1,x_2) = x_2$$u_2(x_1,x_2) = \frac{x_2^2-x_1^2}{2x_1} = c-x_1$

4. Définition de la région

• Frontière extérieure: $\psi_R(x) = 1$ (rayon $R$)
• Frontière intérieure: $\psi_1(x) = 1$ (rayon $1$)
• La région $\Omega$ est en forme de lentille, avec une cuspidale double symétrique à l'origine

Analyse des propriétés du champ vectoriel

1. Composantes de déformation$e_{11} = 0, \quad e_{22} = \frac{x_2}{x_1}, \quad e_{12} = \frac{x_1^2-x_2^2}{4x_1^2}$

2. Représentation en coordonnées polaires En coordonnées polaires $(r,\theta)$, l'équation du cercle est $r = 2c\cos\theta$, les composantes radiale et tangentielle du champ vectoriel: $u_r(c,\phi) = 0, \quad u_\phi(c,\phi) = -c$ montrent que le champ vectoriel est tangent au cercle, de magnitude constante $c$.

3. Décomposition en mouvement rigide$\mathbf{u}(x_1,x_2) = c\mathbf{e}_2 + (x_2\mathbf{e}_1 - x_1\mathbf{e}_2)$ s'exprime comme: translation de $c$ dans la direction $x_2$ + rotation rigide unitaire dans le sens horaire autour de l'origine.

4. Singularité (Proposition 4.1)

• À l'origine, $u_2$ donné par $u_2 = (x_2^2-x_1^2)/(2x_1)$ est indéterminé
• En tendant vers l'origine le long du cercle $\psi_c=1$: $\lim_{x_1\to 0} u_2 = c$
• La valeur limite dépend du cercle le long duquel on tend, donc elle est non unique
• Saut entre différents cercles: $[u_2] = c_1 - c_2$

Méthode semi-inverse et fonction de contrainte d'Airy

1. Fonction de contrainte d'Airy La solution générale de l'équation d'équilibre $\sigma_{\alpha\beta,\beta}=0$ peut être exprimée par la fonction d'Airy $\chi(x_1,x_2)$: $\sigma_{11} = -\chi_{,22}, \quad \sigma_{12} = \chi_{,12}, \quad \sigma_{22} = -\chi_{,11}$

2. Équation de contrôle En définissant le paramètre: $\Lambda(x_1,x_2) := \frac{e_{11}-e_{22}}{e_{12}} = \frac{\sigma_{11}-\sigma_{22}}{\sigma_{12}}$ on obtient une équation hyperbolique indépendante des paramètres de Lamé: $\chi_{,11} - \chi_{,22} - \Lambda(x_1,x_2)\chi_{,12} = 0$

3. Résolution par variables complexes En introduisant $z = x_1 + ix_2$, $\bar{z} = x_1 - ix_2$, l'équation devient: $z^2\chi_{,zz} + \bar{z}^2\chi_{,\bar{z}\bar{z}} = 0$

Dans la classe spéciale $z^2\chi_{,zz} = J(z,\bar{z}) = -\bar{z}^2\chi_{,\bar{z}\bar{z}}$ où: $J(z,\bar{z}) = z\bar{z} + k$$k$ est une constante non négative.

4. Solution de la fonction d'Airy$\chi(z,\bar{z}) = (z\bar{z}-k)\log\frac{z}{\bar{z}}$ En coordonnées polaires: $\chi(r,\theta) = (r^2-k)2i\theta$

Contrainte, déformation et paramètres matériaux

Dérivés de la fonction d'Airy:

Composantes de contrainte: $\sigma_{11} = -2\left[2\theta + \frac{r^2+k}{r^2}\sin 2\theta\right]$$\sigma_{22} = 2\left[-2\theta + \frac{r^2+k}{r^2}\sin 2\theta\right]$$\sigma_{12} = 2\frac{r^2+k}{r^2}\cos 2\theta$

Paramètres de Lamé: $\mu = 4\frac{r^2+k}{r^2}\cos^2\theta$$\lambda = -4\theta\cot\theta - 4\frac{r^2+k}{r^2}\cos^2\theta$

Condition d'ellipticité forte: $\lambda + 2\mu = 4\cos\theta\left(\cos\theta - \frac{\theta}{\sin\theta}\right) + \frac{k}{c^2}$

Quand $k=0$, $\lambda+2\mu < 0$, ne satisfait pas l'ellipticité forte. Quand $k > 0$ et suffisamment grand, on peut assurer $\lambda+2\mu > 0$ et $\mu > 0$, satisfaisant l'ellipticité forte.

Points techniques innovants

1. Systématicité de la construction géométrique: Construire le champ vectoriel à partir des tangentes à la famille de cercles, méthode plus intuitive géométriquement et plus généralisable que celle d'Artstein et Dafermos
2. Technique des variables complexes: Utiliser astucieusement les variables complexes pour transformer l'équation aux dérivées partielles hyperbolique en une forme plus facile à résoudre
3. Contrôle paramétrique: En introduisant le paramètre $k$, on peut contrôler l'ellipticité forte des paramètres de Lamé tout en préservant la singularité du déplacement
4. Matériau non homogène: Prouve que ce problème ne peut être réalisé que dans le cadre de matériaux non homogènes; les matériaux homogènes ne peuvent pas satisfaire l'équation d'équilibre

Configuration expérimentale

Cet article est une recherche théorique mathématique pure, ne comportant pas d'expériences numériques ou de données expérimentales, la recherche étant menée principalement par des déductions mathématiques rigoureuses et des analyses.

Outils d'analyse

1. Calcul analytique: Expressions explicites des déformations, contraintes, paramètres de Lamé
2. Analyse des limites: Étude du comportement limite du champ vectoriel à l'origine
3. Calcul intégral: Calcul de la force totale, du moment total, de l'énergie de déformation
4. Analyse asymptotique: Comportement asymptotique des quantités physiques quand $a\to 0$

Cas de calcul

L'article calcule en détail les quantités physiques suivantes:

1. Paramètres géométriques (Tableau 1):

• Coordonnées des points d'intersection du cercle $\psi_c=1$ avec la sphère $B(0,a)$ et valeurs des fonctions trigonométriques
• Paramétrisation des frontières pour les calculs intégraux ultérieurs

2. Forces de traction aux limites: $F_1 = -4\theta\cos 2\theta$$F_2 = -4\theta\sin 2\theta + 2\left(1+\frac{k}{r^2}\right)$

Résultats expérimentaux

Résultats principaux

1. Force totale (Proposition 6.1) Par intégration sur la région perforée $\Omega_a(c_2,c_1)$ et en prenant $a\to 0$:

• Force totale dans la direction $x_1$: $T_1 = 0$ (car $F_1$ est une fonction impaire de $\theta$)
• Force totale dans la direction $x_2$: $T_2 = 0$

2. Moment total (Proposition 6.1)$\Gamma = 4\pi(R^2-1)$ représente la rotation relative des frontières intérieure et extérieure produisant un moment non nul, proportionnel à la différence des carrés des rayons.

3. Énergie de déformation Sur la région $\Omega_a(1,R)$: $\lim_{a\to 0}|E_a(1,R) - \frac{2k}{a}(R-1)| = 0$

• Quand $k > 0$: l'énergie de déformation $E_a = O(a^{-1})$, singularité à l'origine
• Quand $k = 0$: l'énergie de déformation sans singularité (mais les paramètres de Lamé ne satisfont pas l'ellipticité forte)

4. Ellipticité forte

• $k=0$: $\mu > 0$ mais $\lambda+2\mu < 0$, ne satisfait pas l'ellipticité forte
• $k > 0$ suffisamment grand: $\mu > 0$ et $\lambda+2\mu > 0$, satisfait l'ellipticité forte

Découvertes importantes

1. Préservation de la non-unicité: Même dans le cas où l'ellipticité forte est satisfaite, le comportement limite non unique du déplacement à l'origine persiste, ce qui est une caractéristique essentielle de la singularité géométrique
2. Relation entre moment et énergie:
   • La force totale nulle indique l'équilibre global du système
   • Le moment non nul reflète la rotation relative des frontières
   • L'existence/absence de singularité de l'énergie de déformation est liée à l'ellipticité forte des paramètres matériaux
3. Interprétation physique de la singularité: Le comportement singulier à l'origine peut correspondre à:
   • Rupture ou déchirure du matériau
   • Défaillance fragile
   • Dislocations ou autres défauts
4. Isotropie le long des cercles: Sur chaque cercle $\psi_c=1$, le champ vectoriel est isotrope ($e_{rr}+e_{\theta\theta}=0$)

Travaux connexes

Origines du champ vectoriel

1. Artstein (1983): Introduit une variante de ce champ vectoriel dans les systèmes de contrôle, étudiant les cas où il n'existe pas de contrôle stable lisse
2. Dafermos (2012): Utilise ce champ vectoriel comme exemple d'équation différentielle ordinaire d'évolution, montrant des solutions non uniques non lisses, et propose un critère de "taux d'entropie" pour sélectionner une solution unique
3. Travaux antérieurs des auteurs (2024): Prouve que ce champ vectoriel constitue les trajectoires lagrangiennes des particules fluides dans un écoulement compressible stable, les trajectoires pouvant être non uniques

Fondements de la théorie de l'élasticité

4. Manuels classiques:
   • Gurtin (1972): Théorie de l'élasticité linéaire
   • Landau & Lifshitz (1970): Théorie de l'élasticité
   • Love (1927): Traité sur la théorie mathématique de l'élasticité
5. Ellipticité forte: Knops & Payne (1971) sur les théorèmes d'unicité en élasticité

Théorie des singularités

6. Grisvard (1985), Kondratiev (1967): Problèmes aux limites d'équations elliptiques sur des domaines non lisses, fournissant une base théorique pour comprendre les singularités géométriques

Contributions uniques de cet article

• Innovation méthodologique: Méthode de construction géométrique basée sur les tangentes à la famille de cercles
• Approfondissement théorique: Étude systématique de la singularité dans le cadre de l'élastostatique
• Contrôle paramétrique: Réalisation de l'équilibre entre l'ellipticité forte et la singularité par le paramètre $k$
• Analyse complète: Calculs complets de la force, du moment et de l'énergie

Conclusions et discussion

Conclusions principales

1. Existence: Construction d'une classe de champs de déplacement pour des corps élastiques isotropes compressibles non homogènes linéaires, satisfaisant l'équation d'équilibre sur la région en forme de lentille
2. Caractéristiques de la singularité:
   • Le déplacement est indéterminé à l'origine et possède une limite non unique
   • La valeur limite dépend du chemin d'approche (le long de quel cercle)
   • C'est une conséquence directe de la singularité géométrique (cuspidale double)
3. Ellipticité forte: En choisissant le paramètre $k$ de manière appropriée, on peut assurer que les paramètres de Lamé satisfont la condition d'ellipticité forte tout en préservant le comportement singulier
4. Interprétation physique:
   • La rotation rigide relative de la frontière produit un moment non nul
   • Une rupture, une déchirure ou une défaillance fragile peut se produire à l'origine
   • La singularité de l'énergie de déformation dépend des paramètres matériaux

Limitations

1. Restriction bidimensionnelle: L'analyse actuelle est limitée aux problèmes plans; la généralisation en trois dimensions nécessite l'utilisation de fonctions de Maxwell-Morera
2. Géométrie spéciale: La région en forme de lentille et la famille de cercles sont des choix spéciaux; bien que la méthode soit généralisable, les calculs spécifiques sont complexes
3. Théorie linéaire: Adoption de la théorie de l'élasticité linéaire; les cas de grandes déformations ne sont pas considérés
4. Statique: Seul l'équilibre statique est considéré; le comportement dynamique n'est pas abordé
5. Analyse théorique: Absence de vérification numérique ou de comparaison expérimentale

Directions futures

L'article propose plusieurs directions de recherche dans la section 8:

1. Compréhension approfondie de la singularité: Par analogie avec l'étude de l'équation de Laplace dans les régions en forme de secteur (Grisvard, Kondratiev), récupérer l'unicité par des conditions de croissance appropriées
2. Généralisation géométrique:
   • Considérer des familles de cercles ne se touchant pas à l'origine ($d_2 > 0$)
   • Étudier d'autres familles de courbes (ellipses, paraboles, etc.)
3. Extension du domaine:
   • Prolonger le champ vectoriel à $\Omega(0,1)$ (rotation rigide)
   • Prolonger à la région extérieure $\Omega(R,\infty)$
4. Problème de petit trou: Remplacer la cuspidale par un petit cercle $c_2=\epsilon$, étudier la région avec trou
5. Problème du demi-plan: Prendre $R\to\infty$, étudier la distribution de singularité sur le demi-plan
6. Champ vectoriel général: Appliquer la méthode à un champ vectoriel plan conservatif arbitraire
7. Fonction de contrainte différente: Explorer d'autres solutions de l'équation (5.1)
8. Décomposition de Helmholtz: Exprimer le champ vectoriel comme $\mathbf{u}=\nabla\Psi+\nabla\times\boldsymbol{\psi}$, lier à la théorie des dislocations
9. Généralisation en trois dimensions: Utiliser les fonctions de Maxwell-Morera pour généraliser en trois dimensions
10. Contrôle relaxé: Considérer $\mathbf{u}=cw(x)\mathbf{t}$, où $w(x)$ est une fonction scalaire arbitraire

Évaluation approfondie

Points forts

1. Rigueur mathématique

• Toutes les conclusions sont rigoureusement prouvées
• Les calculs sont détaillés et la logique est claire
• Le traitement des processus limites est minutieux (comme l'analyse asymptotique quand $a\to 0$)

2. Innovativité méthodologique

• La méthode de construction géométrique est nouvelle et intuitive
• L'application de la technique des variables complexes est astucieuse
• La combinaison de la méthode semi-inverse et de la fonction de contrainte d'Airy est efficace

3. Profondeur théorique

• Révèle le lien profond entre la singularité géométrique et la non-unicité de la solution
• Réalise l'équilibre entre l'ellipticité forte et la singularité par le paramètre $k$
• Analyse mécanique complète (force, moment, énergie)

4. Perspicacité physique

• Lie la singularité mathématique aux phénomènes physiques (rupture, déchirure)
• La décomposition en mouvement rigide fournit une image physique claire
• Preuve de la propriété d'isotropie

5. Généralisabilité

• La méthode peut s'appliquer à d'autres familles de courbes
• La solution de la fonction d'Airy peut être généralisée (Appendice B)
• Fournit un cadre pour la recherche sur des problèmes connexes

Insuffisances

1. Absence de vérification numérique

• Pas de simulation numérique soutenant les résultats théoriques
• Pas de visualisation de la distribution des contraintes et déformations
• Manque d'analyse quantitative de l'influence du paramètre $k$ sur la solution

2. Comparaison avec expériences physiques

• Pas de données expérimentales validant l'interprétation de rupture/déchirure
• Manque de comparaison avec le comportement réel des matériaux

3. Spécificité géométrique

• La région en forme de lentille et la famille de cercles sont des choix spéciaux
• Le traitement de géométries plus générales n'est pas suffisamment discuté

4. Scénarios d'application limités

• La théorie linéaire limite les applications pratiques
• La signification pratique des problèmes bidimensionnels est relativement limitée

5. Problème d'unicité

• Bien que la non-unicité soit signalée, aucun critère de sélection n'est proposé
• L'application du critère de "taux d'entropie" de Dafermos aux problèmes d'élasticité n'est pas explorée

6. Aspects de rédaction

• Certains endroits ont une notation complexe (comme plusieurs systèmes de coordonnées)
• Manque de diagrammes illustratifs (seulement deux figures simples)
• L'interprétation physique pourrait être plus complète

Influence

1. Contribution théorique

• Fournit une nouvelle perspective pour l'étude de la singularité dans les corps élastiques non homogènes
• Enrichit la compréhension de la non-unicité des solutions sous la condition d'ellipticité forte
• Relie la singularité géométrique, les propriétés des solutions et les phénomènes physiques

2. Valeur méthodologique

• La méthode de construction géométrique peut inspirer la recherche sur des problèmes connexes
• Démonstration de l'application des variables complexes et de la méthode semi-inverse
• Fournit une référence technique pour traiter les singularités géométriques

3. Applications potentielles

• Mécanique de la rupture: Analyse des contraintes et déformations au point singulier
• Science des matériaux: Conception de matériaux non homogènes
• Théorie des dislocations: Lien avec les dislocations de Volterra (mentionné à la section 8)

4. Limitations

• Les applications pratiques directes à court terme peuvent être limitées
• Nécessite un travail numérique et expérimental supplémentaire
• La généralisation en trois dimensions présente des défis

Scénarios applicables

1. Recherche théorique

• Théorie des singularités des équations aux dérivées partielles
• Problèmes de non-unicité dans l'analyse non linéaire
• Étude de l'influence de la singularité géométrique sur les solutions

2. Science des matériaux

• Analyse des contraintes dans les matériaux non homogènes
• Problèmes d'interface dans les matériaux composites
• Conception de matériaux à gradient fonctionnel

3. Mécanique de la rupture

• Singularité des contraintes à la pointe de fissure
• Mécanismes de déchirure et de défaillance fragile
• Interaction entre fissures multiples

4. Physique mathématique

• Fondements mathématiques de la mécanique des milieux continus
• Solutions non uniques des lois de conservation
• Conditions d'entropie et critères de sélection

5. Fins pédagogiques

• Exemple de la méthode semi-inverse dans la théorie de l'élasticité
• Application des variables complexes aux équations aux dérivées partielles
• Cas concret de singularité géométrique

Références

L'article cite 9 références principales:

1. Artstein (1983): Stabilisation par contrôle de relaxation dans les systèmes de contrôle
2. Dafermos (2012): Dissipation maximale dans les équations d'évolution
3. Gimperlein et al. (2024): Non-unicité dans l'écoulement de fluide plan (travaux antérieurs des auteurs)
4. Grisvard (1985): Équations elliptiques sur des domaines non lisses
5. Gurtin (1972): Théorie de l'élasticité linéaire (Encyclopédie de physique)
6. Knops & Payne (1971): Théorèmes d'unicité en élasticité
7. Kondratiev (1967): Problèmes aux limites dans les régions coniques ou avec points singuliers
8. Landau & Lifshitz (1970): Théorie de l'élasticité
9. Love (1927): Traité sur la théorie mathématique de l'élasticité

Évaluation globale

Ceci est un article de haute qualité en physique mathématique, avec des contributions significatives tant en profondeur théorique qu'en innovativité méthodologique. Par une analyse mathématique rigoureuse, les auteurs révèlent le lien profond entre la singularité géométrique, la non-unicité de la solution et les phénomènes physiques. Bien qu'il y ait encore de la place pour l'amélioration en termes de vérification numérique et d'applications pratiques, en tant que recherche théorique, ce travail fournit des outils mathématiques importants et des perspectives physiques pour les domaines connexes. Il est particulièrement remarquable que l'article non seulement résout un problème spécifique, mais propose également une méthode de recherche systématique et plusieurs directions de recherche futures précieuses, avec un potentiel d'influence académique considérable.