Une balle tombe sur le sol, s'aplatit un instant et rebondit. La quantité d'énergie qu'elle restitue dépend du matériau : une balle en caoutchouc rebondit presque à la même hauteur, tandis qu'un morceau d'argile s'écrase et reste au sol.
Auparavant, les ingénieurs ne savaient calculer précisément ce rebond que pour deux cas idéaux : lorsque la balle se comporte comme un ressort (la force de compression augmente proportionnellement à l'écrasement) ou comme une sphère de Hertz (la force augmente plus rapidement). Pour la diversité des matériaux réels — mousse, plastique, tissus vivants — il fallait recourir à des approximations.
La nouvelle étude fournit une recette générale. Les auteurs ont remarqué que si l'on redessine le graphique « vitesse de la balle en fonction de la compression » dans des axes particuliers, tout matériau se transforme en un simple ressort avec frottement. Il suffit de choisir correctement le coefficient de frottement dans ce système imaginaire. Le rebond se calcule alors par une formule courte, qui ne fait intervenir que deux grandeurs : la rigidité du matériau et sa « viscosité » (les pertes d'énergie lors de la compression).
L'utilité pratique est directe : les jeux vidéo et les simulateurs d'ingénierie pourront calculer les collisions plus rapidement et plus précisément. Il suffit de mesurer une fois deux constantes du matériau, et la machine prédira fidèlement le rebond à n'importe quelle vitesse.
tel que rédigé par les auteurs, en anglais arXiv:2603.27764
We present an exact phase-space analytical treatment of the power-law damped contact oscillator governed by $$m\ddot{\delta} + \alpha\sqrt{mk_H}\,\delta^{(p-1)/2}\dot{\delta} + k_H\delta^p = 0$$, valid for all force-law exponents $$p \geq 1$$ and all initial impact velocities $$v_0$$. The central result is the transformation $$\delta = Ax^{2/(p+1)}$$, where $$A = [(p+1)/2]^{1/(p+1)}$$, which maps the nonlinear phase-space equation $$v\,dv/d\delta + \dots = 0$$ exactly onto a linear spring-dashpot (LSD) system with effective damping ratio $$\alpha_\text{eff} = \frac{\alpha}{\sqrt{2(p+1)}}$$. The phase portrait $$v(\delta)$$, coefficient of restitution $$e$$, and maximum penetration $$\delta_\text{max}$$ follow in closed form. The physical time-domain solution $$(\delta(t), v(t))$$ is obtained parametrically via a single quadrature, which evaluates analytically for $$p=1$$ and at negligible numerical cost for all other $$p$$. We prove that $$e$$ is exactly independent of $$v_0$$ for all $$p \geq 1$$ and derive the universal calibration formula: $$\alpha = \sqrt{2(p+1)}\cdot\frac{-\ln e}{\sqrt{\pi^2 + \ln^2 e}}$$. This generalises the known results for $$p=1$$ (linear spring-dashpot) and $$p=3/2$$ (Hertz contact, Antypov and Elliott, 2011) to the entire power-law family. A closed-form estimate for the critical timestep of explicit time integration is also derived, exhibiting universal scaling with impact velocity and force-law exponent.
Si l'on comprime une balle en caoutchouc jusqu'à la moitié de son diamètre, la force de compression pour de nombreux matériaux croît non pas linéairement, mais comme le carré de la compression. C'est pourquoi la formule générale utilise un changement de variable inhabituel — comme si l'on regardait le monde à travers une loupe qui redresse les courbes des graphiques.