La Physique du Coup Franc : L’Effet Magnus au Service du But

Problématique

Comment la modélisation vectorielle des forces aérodynamiques et la résolution numérique d’équations différentielles permettent-elles d’expliquer la trajectoire incurvée d’un coup franc de football (effet Magnus) et de déterminer les paramètres optimaux pour contourner un mur défensif ?

Introduction

Le football est souvent perçu comme un simple affrontement physique et technique de 22 joueurs sur un rectangle vert. Pourtant, certains gestes techniques frôlent l’œuvre d’art et défient notre intuition visuelle. C’est le cas du coup franc direct. Lorsqu’un joueur comme Lionel Messi ou David Beckham contourne un mur de défenseurs pour loger le ballon en pleine lucarne, la balle semble suivre une trajectoire magique, s’incurvant brusquement en plein vol.

Pour le spectateur, c’est du génie ou du talent pur. Pour le physicien et le mathématicien, c’est une application spectaculaire de la mécanique vectorielle et de la dynamique des fluides. La trajectoire du ballon, loin d’être une simple parabole dictée par la seule gravité, est profondément modifiée par son interaction avec l’air environnant lorsqu’il est mis en rotation. Ce phénomène physique porte un nom : l’effet Magnus.

Dès lors, nous pouvons formuler la problématique suivante : Comment la modélisation vectorielle des forces aérodynamiques et la résolution numérique d’équations différentielles permettent-elles d’expliquer la trajectoire incurvée d’un coup franc de football et de déterminer les paramètres optimaux pour contourner un mur défensif ?

Pour y répondre, nous établirons dans un premier temps le modèle vectoriel du mouvement du ballon à l’aide de la deuxième loi de Newton sous forme d’équations différentielles. Dans un deuxième temps, nous expliquerons l’origine physique de l’effet Magnus à travers le principe de Bernoulli. Enfin, nous résoudrons numériquement ces équations pour analyser un cas concret de coup franc historique.

I. Modélisation vectorielle et équations différentielles du mouvement

A. Paramètres du système et référentiel

Pour étudier le mouvement du ballon de football, nous définissons les éléments suivants :

Le système : Le ballon, assimilé à un point matériel $G$ situé en son centre de gravité, de masse $m \approx 430\text{ g}$ et de rayon $R \approx 11\text{ cm}$.
Le référentiel : Terrestre, considéré comme galiléen sur la courte durée du vol (environ 1 à 2 secondes).
Le repère de coordonnées : Un repère orthonormé cartésien $(O, \vec{i}, \vec{j}, \vec{k})$ où :
- $O$ est le point de départ du ballon au sol.
- $\vec{i}$ est un vecteur unitaire horizontal orienté vers la ligne de but adverse.
- $\vec{j}$ est un vecteur unitaire horizontal orienté vers la gauche (perpendiculaire à $\vec{i}$).
- $\vec{k}$ est un vecteur unitaire vertical dirigé vers le haut.

B. Bilan des forces appliquées

Trois forces s’exercent sur le ballon en mouvement dans l’air :

Le poids $\vec{P}$ : Force gravitationnelle dirigée vers le bas. $$\vec{P} = -m \cdot g \cdot \vec{k}$$
La traînée aérodynamique $\vec{f}_d$ : Force de frottement exercée par l’air, opposée au vecteur vitesse $\vec{v}$. Pour un ballon de football à grande vitesse, cette force est proportionnelle au carré de la vitesse (régime turbulent) : $$\vec{f}_d = -\frac{1}{2} C_x \cdot \rho \cdot S \cdot v \cdot \vec{v}$$ Où $C_x$ est le coefficient de traînée du ballon, $\rho \approx 1{,}2\text{ kg}\cdot\text{m}^{-3}$ la masse volumique de l’air, $S = \pi R^2$ la surface frontale du ballon, et $v = |\vec{v}|$ sa vitesse instantanée.
La force de Magnus $\vec{F}_M$ : Force de portance latérale générée par la rotation du ballon sur lui-même à une vitesse angulaire $\vec{\omega}$. Si l’on applique un effet de rotation purement vertical (axe $\vec{k}$, c’est-à-dire un “effet brossé” latéral), le vecteur rotation s’écrit $\vec{\omega} = \omega \cdot \vec{k}$. La force de Magnus est perpendiculaire à la fois au vecteur vitesse et au vecteur rotation : $$\vec{F}_M = \frac{1}{2} C_L \cdot \rho \cdot S \cdot v^2 \cdot \vec{n}$$ Où $C_L$ est le coefficient de portance (qui dépend de la vitesse de rotation) et $\vec{n}$ est un vecteur unitaire horizontal perpendiculaire à la vitesse de la balle, orienté vers la zone de basse pression.

C. Application de la deuxième loi de Newton

D’après la relation fondamentale de la dynamique (deuxième loi de Newton) : $$\sum \vec{F} = \vec{P} + \vec{f}_d + \vec{F}_M = m \cdot \vec{a}$$

En projetant cette équation vectorielle sur les trois axes du repère, nous obtenons le système d’équations différentielles couplées suivant, régissant la vitesse du ballon : $$\frac{dv_x}{dt} = -\frac{1}{2m} C_x \cdot \rho \cdot S \cdot v \cdot v_x$$ $$\frac{dv_y}{dt} = \frac{1}{2m} C_L \cdot \rho \cdot S \cdot v \cdot v_x$$ $$\frac{dv_z}{dt} = -g - \frac{1}{2m} C_x \cdot \rho \cdot S \cdot v \cdot v_z$$

Où la vitesse instantanée est $v = \sqrt{v_x^2 + v_y^2 + v_z^2}$. La non-linéarité de ces équations différentielles (due au terme en $v$) rend leur résolution analytique impossible. Nous devrons utiliser une méthode de résolution numérique.

II. L’origine physique de l’effet Magnus : Bernoulli et la pression

Pourquoi la rotation du ballon engendre-t-elle cette force perpendiculaire $\vec{F}_M$ ? Pour le comprendre, il faut s’intéresser à la mécanique des fluides et à l’interaction entre la surface rugueuse du ballon et l’air.

A. La couche limite et la différence de vitesse relative

Lorsqu’un ballon avance dans l’air à une vitesse $v$, l’air s’écoule autour de lui de manière symétrique s’il ne tourne pas. Mais si le tireur frappe le ballon de l’extérieur du pied pour lui imprimer une rotation rapide sur lui-même (sens antihoraire vu du dessus, soit $\vec{\omega} = \omega \cdot \vec{k}$ avec $\omega > 0$) :

La surface rugueuse du ballon (les coutures) entraîne par friction une mince couche d’air dans sa rotation : c’est la couche limite.
Du côté gauche du ballon (dans le sens du déplacement), la surface du ballon tourne dans le même sens que l’écoulement de l’air. Les deux vitesses s’additionnent, l’air s’écoule donc plus rapidement par rapport au ballon.
Du côté droit du ballon, la surface tourne en sens inverse de l’écoulement de l’air. Les vitesses se soustraient, l’air s’écoule donc plus lentement.

B. Application du théorème de Bernoulli

Le théorème de Bernoulli appliqué à un écoulement d’air stationnaire et horizontal relie la pression $P$ et la vitesse $v_{\text{air}}$ d’un fluide le long d’une ligne de courant : $$P + \frac{1}{2} \rho \cdot v_{\text{air}}^2 = \text{constante}$$

Ce théorème montre qu’une augmentation de la vitesse du fluide s’accompagne nécessairement d’une diminution de sa pression :

Du côté gauche, la vitesse de l’air est maximale $\implies$ La pression $P_{\text{gauche}}$ diminue.
Du côté droit, la vitesse de l’air est minimale $\implies$ La pression $P_{\text{droite}}$ augmente.

Cette asymétrie crée une différence de pression $\Delta P = P_{\text{droite}} - P_{\text{gauche}} > 0$. Le ballon est alors littéralement “aspiré” vers la gauche, de la zone de haute pression vers la zone de basse pression. C’est cette force d’aspiration transversale qui constitue la force de Magnus et courbe sa trajectoire.

III. Résolution numérique et application à un coup franc réel

A. La méthode d’Euler comme outil de résolution

Pour visualiser la trajectoire, nous devons intégrer le système d’équations différentielles établi en partie I. Nous utilisons la méthode d’Euler, qui permet de calculer de proche en proche les positions et vitesses du ballon à des intervalles de temps très courts $\Delta t$ (par exemple $\Delta t = 0{,}01\text{ s}$) :

$v_x(t + \Delta t) \approx v_x(t) + a_x(t) \cdot \Delta t$
$x(t + \Delta t) \approx x(t) + v_x(t) \cdot \Delta t$

B. Exemple numérique chiffré

Analysons un coup franc tiré à une distance de $25\text{ mètres}$ face à un mur défensif situé à la distance réglementaire de $9{,}15\text{ mètres}$.

Données du ballon : masse $m = 0{,}43\text{ kg}$, rayon $R = 0{,}11\text{ m} \implies S \approx 0{,}038\text{ m}^2$.
Conditions initiales :
- Vitesse initiale $v_0 = 25\text{ m/s}$ (soit $90\text{ km/h}$) projetée vers l’avant.
- Vitesse de rotation $N = 10\text{ rotations/seconde} \implies \omega = 2\pi N \approx 62{,}8\text{ rad/s}$.
- Pour cette rotation, le coefficient de portance estimé est $C_L \approx 0{,}15$.

Calculons la valeur de la force de Magnus au départ du coup franc : $$F_M = \frac{1}{2} C_L \cdot \rho \cdot S \cdot v_0^2 = 0{,}5 \times 0{,}15 \times 1{,}2 \times 0{,}038 \times 25^2 \approx 2{,}14\text{ N}$$

Comparons cette force latérale au poids du ballon : $$P = m \cdot g = 0{,}43 \times 9{,}81 \approx 4{,}22\text{ N}$$

La force de Magnus représente environ $50\ %$ du poids du ballon. C’est une force considérable, capable d’infléchir fortement la trajectoire.

C. Calcul de la déviation latérale

À cette vitesse, le temps de vol pour parcourir les $25\text{ mètres}$ est d’environ $t \approx 1{,}0\text{ seconde}$. L’accélération latérale moyenne $a_y$ subie par le ballon est : $$a_y = \frac{F_M}{m} = \frac{2{,}14}{0{,}43} \approx 5\text{ m/s}^2$$

La déviation latérale $y(t)$ après une seconde de vol s’obtient par double intégration : $$y(1) \approx \frac{1}{2} \cdot a_y \cdot t^2 = 0{,}5 \times 5 \times 1^2 \approx 2{,}5\text{ mètres}$$

Une déviation de $2{,}5\text{ mètres}$ vers la gauche est largement suffisante pour contourner le mur défensif (généralement large de $2$ mètres) et revenir se loger dans le cadre du but, rendant le tir totalement imparable pour le gardien de but resté au centre.

Conclusion

En conclusion, la trajectoire courbe d’un coup franc n’a rien de magique : elle s’explique de manière déterministe par la physique et les mathématiques. La deuxième loi de Newton modélise l’influence conjointe de la gravité, de la traînée de l’air et de la force de Magnus. Cette dernière trouve son origine physique dans la différence de vitesse d’écoulement de l’air induite par la rotation de la balle, créant une différence de pression modélisée par le théorème de Bernoulli.

Bien qu’un joueur de football n’effectue pas consciemment des calculs différentiels ou de double intégration sur le terrain, ses muscles et son cerveau ont intégré ces lois physiques à force d’entraînement. En ajustant l’angle d’impact, la force de frappe et la zone de contact pour maximiser la rotation $\omega$, le footballeur résout, par le geste, un problème complexe de mécanique des fluides pour tromper la vigilance de ses adversaires.

Questions Potentielles du Jury

1. Comment modéliseriez-vous le tracé de cette trajectoire en Python en utilisant la méthode d’Euler ?

Réponse : Pour modéliser la trajectoire sous Python, on définit une boucle temporelle qui actualise à chaque itération la position, la vitesse et l’accélération du ballon. Voici un exemple simple d’algorithme :

import numpy as np

# Constantes
g = 9.81; rho = 1.2; m = 0.43; R = 0.11; S = np.pi * R**2
Cx = 0.25; CL = 0.15; dt = 0.001

# Initialisation des vecteurs position et vitesse
r = np.array([0.0, 0.0, 0.0]) # x, y, z
v = np.array([25.0, 0.0, 3.0]) # v_x, v_y, v_z (légèrement vers le haut)

trajectory = []

for t in range(1500): # 1.5 seconde de simulation
    v_norm = np.linalg.norm(v)
    # Calcul des forces
    P = np.array([0, 0, -m*g])
    Fd = -0.5 * Cx * rho * S * v_norm * v
    # Force de Magnus (perpendiculaire à v dans le plan horizontal si rotation selon z)
    FM = np.array([-v[1], v[0], 0]) 
    FM = FM / np.linalg.norm(FM) * (0.5 * CL * rho * S * v_norm**2)
    
    # Deuxième loi de Newton
    a = (P + Fd + FM) / m
    
    # Méthode d'Euler
    v = v + a * dt
    r = r + v * dt
    trajectory.append(r.copy())

Cet algorithme permet d’obtenir les coordonnées $x(t)$, $y(t)$, $z(t)$ du ballon pas à pas et de les afficher sous forme de courbe 3D.

2. Pourquoi le théorème de Bernoulli utilisé ici s’applique-t-il sur l’air, et quelles sont les limites de cette simplification ?

Réponse : Le théorème de Bernoulli s’applique sous l’hypothèse d’un écoulement fluide parfait (viscosité nulle), incompressible (densité constante) et stationnaire.

Dans notre cas, l’air est assimilé à un fluide parfait pour simplifier la modélisation à l’échelle globale.
Cependant, en réalité, l’air a une viscosité non nulle. C’est précisément cette viscosité qui permet à l’air “d’adhérer” à la surface du ballon et de former la couche limite en rotation. Sans viscosité, il n’y aurait aucun entraînement de l’air, donc pas d’effet Magnus (ce qu’on appelle le paradoxe de d’Alembert). Le modèle utilisant Bernoulli est donc une approximation commode pour expliquer qualitativement la différence de pression, mais les calculs rigoureux de traînée et de portance doivent intégrer les équations de Navier-Stokes.

3. La force de frottement de l’air est modélisée par une force quadratique (proportionnelle à $v^2$). Pourquoi ne pas utiliser une force linéaire (proportionnelle à $v$) ?

Réponse : Le choix entre un modèle de frottement linéaire ou quadratique dépend d’un nombre sans dimension appelé le nombre de Reynolds ($Re$), qui compare les forces d’inertie aux forces de viscosité : $$Re = \frac{\rho \cdot v \cdot d}{\eta}$$ Où $d = 2R$ est le diamètre du ballon et $\eta \approx 1{,}8 \times 10^{-5}\text{ Pa}\cdot\text{s}$ la viscosité dynamique de l’air.

Pour un ballon de football à $v = 25\text{ m/s}$ : $$Re = \frac{1{,}2 \times 25 \times 0{,}22}{1{,}8 \times 10^{-5}} \approx 3{,}7 \times 10^5$$ Un nombre de Reynolds aussi élevé ($Re \gg 1000$) indique un écoulement hautement turbulent autour du ballon où les forces d’inertie de l’air dominent largement. Dans ce régime, la traînée est proportionnelle au carré de la vitesse ($v^2$). Le frottement linéaire est réservé aux écoulements laminaires et lents (comme une petite bille tombant dans l’huile, où $Re < 1$).

4. Qu’est-ce que la crise de traînée du ballon de football et quel est son impact sur le coup franc de Roberto Carlos en 1997 ?

Réponse : La crise de traînée est un phénomène où le coefficient de traînée $C_x$ d’une sphère chute brutalement à partir d’un certain nombre de Reynolds critique (autour de $Re \approx 2 \times 10^5$). Lors du célèbre coup franc de Roberto Carlos en 1997 (frappe surpuissante à $137\text{ km/h}$, soit $38\text{ m/s}$) :

Au début de la trajectoire, la vitesse est très élevée, le ballon est au-delà du Reynolds critique. La couche limite est turbulente, ce qui réduit le coefficient de traînée $C_x$. Le ballon ralentit peu et la force de Magnus est relativement faible par rapport à l’énergie cinétique. La trajectoire est presque rectiligne.
En cours de vol, les frottements ralentissent le ballon. Dès qu’il repasse sous le Reynolds critique, l’écoulement redevient laminaire autour de la sphère, le $C_x$ remonte brusquement et le coefficient de portance $C_L$ augmente nettement.
Le ballon subit alors une décélération rapide vers l’avant et une déviation latérale de Magnus beaucoup plus violente à la fin du tir, donnant cette impression de courbe brisée ou de “banane” qui a surpris le gardien de but français Fabien Barthez.

5. En quoi ce sujet s’inscrit-il dans votre projet d’orientation ?

Réponse : Ce sujet illustre parfaitement l’intérêt de la physique théorique et de la programmation informatique appliquées à l’aérodynamique et à la modélisation industrielle. Je souhaite m’orienter vers des études d’ingénieur en aéronautique, en mécanique des fluides ou en simulation numérique (par exemple via une classe préparatoire MPSI/PSI ou une école d’ingénieurs comme l’ISAE-SUPAERO or l’ENSEEIHT). La modélisation de trajectoires d’objets en rotation dans un flux d’air est au cœur de la conception des pales d’éoliennes, des ailes d’avions ou de la pénétration dans l’air des carrosseries automobiles. Ce projet me permet d’aborder des concepts fondamentaux d’intégration numérique et de mécanique des fluides que j’aurai à manipuler au cours de ma formation supérieure.

Attention : Ce sujet n'est qu'un exemple !

La Physique du Coup Franc : L'Effet Magnus au Service du But