Comment intégrer l'expression polynomiale sur un élément 3D à 4 nœuds?

12

Je souhaite intégrer une expression polynomiale sur un élément à 4 nœuds en 3D. Plusieurs livres sur FEA couvrent le cas où l'intégration est effectuée sur un élément arbitraire à 4 non-plats. La procédure habituelle dans ce cas est de trouver la matrice de Jacobi et d'utiliser son déterminant pour changer la base d'intégration en celle normalisée dans laquelle j'ai les limites d'intégration plus simples [-1; 1] et la technique de quadrature de Gauss-Legendre est facilement utilisée.

En d'autres termes, est réduit à la forme deSf(x,y) dxdy1111f~(e,n) |det(J)|dedn

Mais dans le cas 2D, je change l'élément arbitraire plat en un élément plat mais bien carré 2 par 2.

L'élément 3D à 4 nœuds n'est pas plat en général, mais je suppose qu'il peut toujours être mappé avec le système de coordonnées 2D qui est en quelque sorte lié au système de coordonnées cartésiennes. Je ne peux pas comprendre comment exprimer {x, y, z} en termes de {e, n} et quelle serait la taille de la matrice de Jacobi dans ce cas (elle est censée être carrée).

Domaines 2D et 3D

danny_23
la source

Réponses:

8

Vous intégrez une fonction sur une variété à 2 dimensions intégrée dans ; Les livres d'analyse sur les variétés (comme le livre accessible de Munkres ou les livres de Lee sur les variétés) sont utiles pour discuter de la théorie définissant ce type d'intégrale.R3

Supposons que soit une fonction à valeur réelle définie sur la variété , qui est votre élément 3D à 4 nœuds.fM

Vous souhaitez calculer:

MfdS.

Supposons que est une fonction qui mappe à . alorsφ[1,1]2M

MfdS=[1,1]2f(φ(x,y))(det(DφT(x,y)Dφ(x,y)))1/2dxdy

(J'ai utilisé cet ensemble de notes pour rafraîchir ma mémoire.) Ci-dessus, est la matrice jacobienne de et est sa transposition.DφφDφT

Une fois que vous pouvez écrire l'intégrale sur , vous pouvez alors utiliser des méthodes numériques pour l'évaluer.[1,1]2

Certains commentaires:

  • Je suis presque sûr que votre élément 3D à 4 nœuds est un collecteur. Si tel est le cas, la fonction existe (par définition), est continue par morceaux (pour les variétés topologiques) et est inversible. C'est à vous de trouver une fonction avec ces propriétés.φ
  • L'argument ci-dessus suppose que est une variété lisse, ce qui implique qu'il existe un qui est continuellement différenciable. Dans votre cas, l'élément que vous décrivez peut ne pas être différenciable en continu. Si cela est vrai, vous pourriez probablement toujours partitionner votre manifold en deux manifolds lisses, et l'argument ci-dessus tient toujours. Encore une fois, vous devez trouver satisfaisant les propriétés d'invertibilité et de différentiabilité continue.Mφφ
Geoff Oxberry
la source
Merci beaucoup. Le livre que je lis ne couvre que le cas où une matrice Jacobi carrée (2 par 2) est impliquée pour garder les choses simples. L'expression ci-dessus si je l'ai bien fait permet d'utiliser des matrices Jacobi de taille arbitraire (2 par 3). Malheureusement, je reçois toujours pour le moment mais c'est beaucoup mieux que ce que j'avais auparavant. Je vais créer un autre fil sur le bon choix de la fonction de mappage. Merci encore. det(DφT(x,y)Dφ(x,y))=0
danny_23
3
Votre matrice jacobienne doit être 3 par 2, donc doit être une matrice 2 par 2. DφDφTDφ
Geoff Oxberry
2
Geoff, c'est exact. Je mets simple formule générale , plus un exemple travaillé ici: theoretical-physics.net/dev/src/math/integration.html
Ondřej Čertík