Intégrales en physique/Intégrales multiples

Intégrale double

On peut généraliser l’idée précédente à des problèmes à plusieurs dimensions. Par exemple, au lieu de considérer les contributions de points présents sur une ligne, on peut considérer la contribution de points d'une surface.

Exemple tiré de la mécanique

Une pièce d'une machine est en contact avec le plan (Oxy) du sol. La surface de contact est appelée Σ.

Lorsque la machine est en marche, elle exerce sur le sol en tout point M de Σ une pression p(M)=p(x,y) dans la direction $\overrightarrow{z}$, où (x,y) sont les coordonnées de M.

Sur une surface élémentaire d'ordre 2 ${\mathrm{d}}^{2}S=\mathrm{d}x\mathrm{d}y$ au voisinage de M, la force (infinitésimale d'ordre 2) exercée par la pièce sur le sol vaut alors ${\mathrm{d}}^2\overrightarrow{F}(x,y)=-p(x,y)\overrightarrow{z}{\mathrm{d}}^{2}S$.

Pour trouver la force totale s'exerçant sur le sol, il faut sommer les contributions de tous les points de Σ. Comme cette somme se fait dans deux dimensions : x et y, il faut intégrer deux fois:

• Pour un y donné, on somme les contributions de tous les points M(x,y) tels que M soit dans Σ (intégrale suivant x). On obtient la force qui s'exerce sur la surface infinitésimale d'ordre 1 dS, de largeur dy et de longueur x₂(y)-x₁(y) : $\mathrm{d}\overrightarrow{F}(y)=-\mathrm{d}y\underset{x_1(y)}{\overset{x_2(y)}{\int }}p(x,y)\mathrm{d}x\overrightarrow{z}$

• On fait ensuite l'intégrale sur y pour sommer les contributions de toutes les bandes horizontales. On obtient au final $\overrightarrow{F}=\underset{y_{min}}{\overset{y_{max}}{\int }}\mathrm{d}\overrightarrow{F}(y)=-\underset{y_{min}}{\overset{y_{max}}{\int }}\underset{x_1(y)}{\overset{x_2(y)}{\int }}p(x,y)\mathrm{d}x\mathrm{d}y\overrightarrow{z}$.

Cette dernière expression rend compte de la somme des contributions de tous les points de la surface Σ, de dimension 2. Pour condenser l'écriture, on utilise une intégrale double pour noter $\overrightarrow{F}=\underset{\mathrm{\Sigma }}{\iint }-p(x,y)\mathrm{d}x\mathrm{d}y\overrightarrow{z}$.

On peut bien entendu choisir de faire l'opération dans l'autre sens : intégrer suivant y à x choisi, puis sur x.

Exemple tiré de la statique des fluides

Considérons le barrage de Bimont. Il retient l'eau en grande quantité.

L'eau retenue exerce des efforts de pression sur une partie Σ de la surface du barrage. Seulement voilà, tous les points de la masse d'eau retenue ne sont pas à la même profondeur (distance à la surface de l'eau) h. La force de pression n'est donc pas la même en tout point de la surface Σ du barrage en contact avec l'eau.

Supposons pour simplifier que le barrage est :

• rectangulaire
• inclus dans le plan (Oyz)
• de hauteur H
• de longueur L

Tous les points de Σ qui sont à la profondeur h sont à la même pression $p(h)=p_{atm}+\rho gh$.

La normale à Σ est $\overrightarrow{x}$ en tout point M de Σ.

• Si on considère une surface infiniment petite d'ordre 2 d²S, de longueur infinitésimale d'ordre 1 dy et de hauteur infiniment petite d'ordre 1 dh située à la profondeur h, la force infinitésimale ${\mathrm{d}}^2{\overrightarrow{F}}_e(y,h)$ qu'exerce l'eau sur d²S vaut ${\mathrm{d}}^2{\overrightarrow{F}}_e(y,h)=p(y,h)\mathrm{d}y\mathrm{d}h\overrightarrow{x}=(p_{atm}+\rho gh)\mathrm{d}y\mathrm{d}h\overrightarrow{x}$.
• Comme p ne dépend pas de y, l'intégration sur y donne directement, à la profondeur h, $\mathrm{d}{\overrightarrow{F}}_e(h)=\underset{y=0}{\overset{y=L}{\int }}p(y,h)\mathrm{d}y\mathrm{d}h\overrightarrow{x}=\underset{y=0}{\overset{y=L}{\int }}(p_{atm}+\rho gh)\mathrm{d}y\mathrm{d}h\overrightarrow{x}=L(p_{atm}+\rho gh)\mathrm{d}h\overrightarrow{x}$.
• Ensuite, sur la surface infiniment petite d'ordre 1 dS, de longueur L et de hauteur infiniment petite dh située à la profondeur h, la force infinitésimale $\mathrm{d}{\overrightarrow{F}}_e(h)$ qu'exerce l'eau sur dS vaut $\mathrm{d}{\overrightarrow{F}}_e(h)=p(h)\mathrm{d}S\overrightarrow{x}=(p_{atm}+\rho gh)L\mathrm{d}h\overrightarrow{x}$.
• De l'autre côté du barrage, l'air exerce une force de pression $\mathrm{d}{\overrightarrow{F}}_a(h)=-p_{atm}\mathrm{d}S\overrightarrow{x}$.
• Au total, la surface dS est soumise à $\mathrm{d}\overrightarrow{F}(h)=\rho gh\mathrm{d}S\overrightarrow{x}$.
• Pour trouver l'effort exercé sur l'intégralité de Σ, il faut « sommer les contributions de toutes les forces de pression s'exerçant aux profondeurs comprises entre 0 et H sur le barrage » :

$\begin{array}{cc}\overrightarrow{F}& =\underset{0}{\overset{H}{\int }}\mathrm{d}\overrightarrow{F}(h)\\ & =\underset{0}{\overset{H}{\int }}\rho gh\mathrm{d}S\overrightarrow{x}\\ & =\underset{0}{\overset{H}{\int }}\rho gh(L\mathrm{d}h)\overrightarrow{x}\\ & =\rho gL\underset{0}{\overset{H}{\int }}h\mathrm{d}h\overrightarrow{x}\\ & =\rho gL\overrightarrow{x}\cdot {\left[\frac{h^2}{2}\right]}_{h=0}^{h=H}\\ \end{array}$

Remarque : En remarquant que p est indépendant de y, on aurait pu commencer directement la résolution en considérant la surface infiniment petite d'ordre 1 dS=L dh, sur laquelle s'exerce $\mathrm{d}{\overrightarrow{F}}_e(h)=p(h)\mathrm{d}S\overrightarrow{x}=(p_{atm}+\rho gh)L\mathrm{d}h\overrightarrow{x}$.

Intégrale triple

On peut poursuivre la généralisation dans l'espace. Voyons sur un exemple tiré de l'électrostatique.

Soit une distribution de charges réparties dans un volume V telle qu'en un point courant M de V, la densité de charge volumique vale $\rho (M)$. Cela signifie que pour un volume élémentaire d³τ autour d'un point M, la charge que porte le volume d³τ vaut d³q = ρ d³τ.

En un point P, le champ électrostatique généré par la charge d³q présente en M vaut ${\mathrm{d}}^3\overrightarrow{E}(P)=\frac{{\mathrm{d}}^{3}q{\overrightarrow{u}}_r}{4\pi {\epsilon }_0{\mathrm{M}\mathrm{P}}^2}=\frac{\rho (M){\mathrm{d}}^3\tau }{4\pi {\epsilon }_0{\mathrm{M}\mathrm{P}}^2}{\overrightarrow{u}}_r$, avec ${\overrightarrow{u}}_r$ vecteur unitaire de même sens et même direction que $\overrightarrow{\mathrm{M}\mathrm{P}}$.

Pour obtenir le champ généré par l'intégralité des charges du volume V, on somme les contributions de tous les volumes élémentaires de V, ce qui se fait à l'aide d'intégrales suivant 3 dimensions : x, y et z.

On note alors cette opération à l'aide d'une intégrale triple sur le volume V : le champ électrostatique créé en P par la distribution vaut $\overrightarrow{E}(P)=\underset{V}{\iiint }\frac{\rho (M){\overrightarrow{u}}_r}{4\pi {\epsilon }_0{\mathrm{M}\mathrm{P}}^2}{\mathrm{d}}^3\tau$ avec ${\overrightarrow{u}}_r$ vecteur unitaire de même sens et même direction que $\overrightarrow{MP}$.

Ordre des différentielles

On remarque, lors de la manipulation des intégrales, la manipulation constante des différentielles d'ordres différents. Prêter attention à l’ordre des différentielles est un bon moyen pour dépister les erreurs de calcul ou de manipulation. Les ordres permettent de vérifier l'« homogénéité » des résultats. Par exemple, un infiniment petit d'ordre 1 ne peut pas être égal à une grandeur macroscopique.

Modèle:Bas de page