Pendule simple discret

Soit un pendule simple, c’est-à-dire un point matériel, M, de masse m, astreint à se déplacer sur un cercle vertical(C), de centre O, de rayon l, dans un champ de pesanteur g uniforme.

C'est donc un cas particulier de mouvement dans un puits de potentiel.

Fait remarquable (et peu connu) : on peut étudier ce problème redoutable de fonctions elliptiques à un niveau élémentaire de géométrie à condition de discrétiser le mouvement : il s'agit du pendule simple discret.

Le cas intégrable

Pour apprivoiser le sujet, décrivons le cas intégrable : celui où le pendule parti du point le plus bas du cercle ( en A) possède la vitesse $V_o = \sqrt{2gL}$ et atteindra le point le plus haut du cercle (en B) au bout d'un temps infini.

Voici la procédure :

Matériel : crayon, règle, compas.

• Construction : Soit (C) la trajectoire de A en B , demi-cercle de diamètre vertical AB = 2L, de centre O.

Soit O₂, situé à la verticale de O; OO₂ = d . Tracer le demi-cercle (C₂) de rayon O₂B = L − d, qui recoupe la verticale en B': AB' = 2d. (Par exemple l = 10cm et d = 1cm ).

Depuis A , mener la tangente à (C₂) qui recoupe (C) en A₂. Depuis A₂, mener la tangente à (C₂) qui recoupe (C) en A₄. Continuer avec soin jusqu'à A₈, voire A₁₀. Depuis B',tracer la tangente horizontale à (C₂) qui coupe (C) en A₁, de cote h = 2d. Depuis A₁ , mener la tangente à (C₂) qui recoupe (C) en A₃ . Depuis A₃ , mener la tangente à (C₂) qui recoupe (C) en A₅. Continuer avec soin jusqu'à A₉.

• Théorème du pendule discret, cas intégrable:

Les points A_n sont séparés par des intervalles de temps égaux. On ne peut jamais atteindre B.

• Vérification expérimentale : la suite polygonale des A_n doit être circonscrite à un cercle (c₁) de centre o₁ ( compris entre O et O₂), tangent en B à (C) et (C₂).

De même, la suite $A_0 \; A_3 \; A_6 \; A_9$ : cercle (C₃) de centre O₃. On est en général assez fier de le tracer! Compléter l'arbelos. Les segments radiaux, issus des A_n, découpent l'arbélos en quadrangles curvilignes ; les griser en deux valeurs de gris  :la figure donne une idée de la dynamique du pendule qui "s'essouffle en montant". Les élèves de 12-15ans sont généralement contents de l'esthétique de leur figure.

Pour les élèves de lycée, on peut aller un peu plus avant :

Considérons par exemple la suite $A_0 \; A_3 \; A_6 \; A_9$ , tangente en T₀ , T₃ , T₆ , T₉ au cercle (C₃):

La module de la vitesse en A₀ est A₀T₀: tracer le vecteur vitesse en A₀.

De même, pour A₃T₃ , A₆T₆ et A₉T₉ ( = A₉T₆ bien sûr !)

Cette propriété est due au fait suivant : en appelant H_n les projections des A_n sur l'axe radical du faisceau de cercles, on a H_nA_n.2d = (A_nT_n)², donc A_nT_n est proportionnel au module de la vitesse. On peut donc construire aisément beaucoup de point du diagramme horaire de $\dfrac{ds}{dt} = v(t)$: la forme caractéristique du soliton apparaît clairement ( cf pendule simple).

Au niveau Deug, un exercice classique est : trouver la démonstration sur laquelle s'appuie cette construction géométrique.

Le cas du tournoiement lent

On suppose la vitesse V₀ très légèrement plus grande.On se doute qu'avec une énergie supérieure de 10 ^{− 37} joules à 2mgL, le pendule va tournoyer, mais sans que vraiment on puisse distinguer expérimentalement avec le cas précédent. Donc le pendule présentera entre l'intervalle très long d'une période, un phase de vitesse rapide (avec une vitesse v(t) quasi-égale à celle du soliton et dont l'intégrale sera 2π).

Expérience

On la fait avec un pendule de Mach : la boule est lancée de B avec une vitesse minime. Une caméra filme le mouvement jusqu'au moment où l'on passe au régime de grandes oscillations à cause de la très faible (mais impossible à éliminer) déperdition d'énergie. On colle informatiquement les photos prises à des temps réguliers : on aura ainsi plusieurs lots. Il est facile de constater que les résultats précédents sont vrais à l'erreur expérimentale près...

Le cas du tournoiement très rapide

Soit toujours (C) le cercle de tournoiement, de diamètre vertical A₀B. Soit H le point de cote OH telle que l'énergie soit mh.OH :

la vitesse en B sera $v(B) = \sqrt{ \left( 2g.OH \right) } \left( 1 - \dfrac{L}{2.OH} \right)$; celle en A, légèrement plus grande $v(A) = \sqrt{2g.oh} \left( 1 + \dfrac{L}{2.OH} \right)$;

c’est-à-dire que cette fois le diagramme de vitesse sera extrêmement plat, avec une vitesse moyenne très légèrement inférieure à $\sqrt{2g.OH}$  : comme on connaît v(s), on sait résoudre ce cas par la méthode du diagramme horaire.

Le tournoiement : construction d'Euler

Rappelons la relation d'Euler (voir complément en discussion)dans un triangle quelconque de cercle circonscrit (C), de centre O, de rayon R , et de cercle inscrit, de centre I, de rayon r:

R² = (OI)² + 2Rr

On va donner un cas facile à tracer d'un tournoiement assez rapide, correspondant à une épure calculée, selon la règle précédente:

Tracer le cercle (C) de rayon 100 (cm).Ce sera la trajectoire.

Tracer le cercle (C₂) de rayon 32 de centre O₂ tel que OO₂ = 60.

La tangente horizontale haute correspond à deux points A₂ et A₄ ( la corde $A_2A_4 = \sqrt{6144} \simeq 78,4$ ). La tangente basse correspond à deux points A₁ et A₇ ( la corde $A_1A_8 \simeq 196,7$ ).

Les points $A \; A_1 \; A_2 \; B \; A_4 \; A_5 \; A$ sont atteints à des dates entières. On peut tracer le cercle inscrit à cet hexagone (dit de Euler-Poncelet), soit (C₁) de centre O₁. Soient les points de tangence $T_0 \; T_1 \; T_2$ et leurs symétriques. On obtient les 4 vitesses : en A ,le segment AT₀ , en A₁, le segment A₁T₁ , en A₂, le segment AT₂, et en B le segment BT₂.

Soit H(M) la projection de M sur l'axe radical (D) du faisceau de cercles à points limites (dits de Poncelet en France, dits de Landen ailleurs): la vitesse en M est proportionnelle à $\sqrt{HM}$ conformément à la règle de Baliani-Torricelli (le proto-théorème de l'énergie cinétique, cinquante ans avant Leibniz!): ici OH₀ = 104,8cm . La puissance de H₀ par rapport au faisceau est 983,04cm², correspondant aux deux points de Poncelet L' de cote 136.15 cm et L (point de Landen, ici), de cote : 73.45 cm. Si la construction est correcte, alors A₁A₂ et A₄A₅ s'intersectent en L' ; A₁A₄ et A₂A₅ en L.

Vérification graphique : Proposons la vérification suivante: prendre par convention un pas de temps 0.5 unité. Construire le point sur le cercle M₃ tel que BM₃ = 0,5BT₃. Construire aussi précisément que possible l'hexagone circonscrit M₃M₄M₅M₀M₁M₂ : on aura le plaisir de le voir se refermer en M₃. Tracer les arcs AM₀, A₁M₁, A₂M₂, BM₃ , A₄M₄, A₅M₅ : on devrait avoir le plaisir de voir A₁M₁ < A₅M₅, légèrement ( on peut comprendre : dans le même temps M monte de A₁ en M₁ mais descend de A₂ en M₂ ), et surtout les 3 diagonales se coupent en L, évidemment.

Bien sûr, on peut obtenir les vitesses en M_n, par la même méthode que pour les A_n, avec la troublante satisfaction de voir que le compas ramène de T₀, T₁ , T₂ , T₃, T₄ , T₅ à T₀ ! et bien sûr comme l'arc M₀M₁ est "très vertical", un rapport EXACT des vitesses $\dfrac{M_1T_0}{M_0T_0}$, bien égal à $\sqrt{ \dfrac{H_1M_1}{H_0M_0}}$.

On complètera à loisir ...

Le tournoiement : construction de Landen-Poncelet

Au lieu de jouer avec un triangle, on joue avec un trapèze. La règle des quadrangles de Poncelet est plus compliquée que celle d'Euler pour les triangles, on peut la retrouver dans l'article polygone de Poncelet : le rayon r du cercle inscrit de centre I, avec IO = d est tel que : $r \sqrt{2R^2 + 2d^2} = R^2 - d^2$.

Mais on peut "jouer simple", en acceptant la coupe diamantaire suivante :

Construire le cercle (C), de centre O, de rayon 100, de diamètre horizontal A₁A₇, et le diamètre vertical AB, où l'on portera OM = 127,2.

Le triangle MA₁A₇ recoupe (C) en A₃ et A₅ . Finir la construction du quadrangle orthocentrique en traçant le point L ( point de Landen).

Le milieu HO de ML définit l'axe radical horizontal des cercles à points limites M et L.

Construire le cercle (C₂) inscrit dans le trapèze A₁A₃A₅A₇.

Tracer l'horizontale passant par L qui coupe (C) en A₂ et A₆ ( on vérifiera A₂A₆ = 120). On constatera avec joie que AA₂BA₆ est circonscrit à (C₂).

On parachèvera par la construction du cercle (C1) tangent en 8 points à l'octogone AA₁A₂A₃BA₅A₆A₇, ce qui permet la détermination des vitesses.

Les huit points obtenus sont atteints à dates entières.

On peut prendre le même plaisir ( jouer du compas) qu'avec l'hexagone d'Euler, mais les constructions sont plus délicates ; ainsi en est-il de l'octogone de tangente horizontale basse : on obtient alors 16 points isochrones avec leurs vitesses.

Cette règle expérimentale d'antan s'est perdue ? certes, il fallait savoir manier règle et compas. Mais plus vraisemblablement, la virtuosité des logiciels montre que la construction n'est qu'approchée : OM n'a pas exactement la bonne cote!(cf discussion)

Néanmoins pour la présentation du mouvement à des scolaires, la figure est patente : diminution drastique de vitesse du point B₁ au point B₂ car la corde est quasi verticale. A contrario, passage de B₄ à B₅ en passant par B (alias A₄ ) à vitesse quasi uniforme : le mouvement est ralenti, mais la vitesse initiale en B est grande devant la chute d'altitude : on comprend pourquoi le régime "logarithmique" du cas limite a du mal à s'installer.

Le tournoiement : cas général

Tracer un 2n-polygone de Poncelet dont deux côtés soient horizontaux tangents au cercle (C₁) "adéquat". Tracer alors la 2n-chaîne de Poncelet passant par A et B : les 2n points seront atteints à dates entières. La difficulté évidemment réside dans le mot "adéquat", qui cache la théorie de la duplication des fonctions elliptiques. Si cos(nt) s'exprime aisément à l'aide des polynômes de Tchebycheff à l'aide de cos(t), de même c_n(nt,k) est algébrique en c_n(t) et k², mais la relation est autrement difficile!

Néanmoins, on a compris : le point H₀ de l'axe radical des deux cercles (C) et (C₁) , de cote h est tel que la vitesse V₀ en A est ${V_0}^2 = 2d.h$

Depuis A, on trace la chaîne de Poncelet, qui représentent des points successifs du pendule simple à des temps égaux t₀ : en général,t₀ n'est pas un rationnel fois la période T(h), et la chaîne de Poncelet ne se referme jamais.

Mais évidemment aussi, par continuité, il existe aussi des cercles (C₁) pour lesquels la chaîne de Poncelet se ferme : par contre calculer la distance OO₁ = d et le rayon R₁ correspondant exige de faire le calcul via les fonctions elliptiques ( cf APPELL et LACOUR ou cf GREENHILL).

Le cas des petites oscillations

Le cas pendulaire, on en aura eu l'intuition, est le cas où l'axe radical coupe (C) en C et C' qui seront donc les points d'altitude maximale atteints par le point M dans son oscillation pendulaire. Le point H₀ sera maintenant le milieu de la corde CC': AH₀ = h.

Dans le cas pendulaire, traçons donc le cercle (B) de centre B, de rayon BC. Si on joint AC, évidemment AC est tangente à B , et l'arc AB correspond à une durée K(h), quart de période. La tangente horizontale SS' au cercle (B) correspond forcément aux points de date $\dfrac{K}{2}$, $\dfrac{3K}{2}$ pour S et $\dfrac{5K}{2}$, $\dfrac{7K}{2}$ pour S'.

Puisque $h \ll l$ , on peut approximer arcs et cordes :

AC = l(α) ;et BC = lcos(α), puis $AS = l ( \alpha) \dfrac{\sqrt{2}}{2}$.

Ce qui correspond bien au 1/8 ème de période, pour une oscillation qui est donc quasi sinusoïdale. A déjà été évoqué dans l'article pendule simple, le raisonnement de Torricelli, très proche de ce tracé.

Le cas des très grandes oscillations

Cette fois, on a au contraire $l-h \ll l$. Néanmoins la construction des points S et S' reste identique. Il est patent de reconnaître la différence des longueurs d'arcs huitième de période. Il est aussi facile de reconnaître la vitesse en S par rapport à celle en A.

Par contre , pas de moyen simple d'évaluer la valeur de la période 4K(h) ; et surtout pas moyen d'évaluer les hauteurs h en oscillation ( h < l) et h' en tournoiement ( h' > l ) qui ont les mêmes périodes; alors qu'évidemment la théorie mathématique donne la correspondance entre ces valeurs.

Quelques cas particuliers

Bien sûr , le cas de l'oscillation d'amplitude 90° est auto-similaire à son complément. Le raisonnement de la symétrie de Corinne joue à plein dans ce cas, et K = K'

Le cas particulier le plus connu et utilisé en travaux pratiques est celui du pendule oscillant avec une élongation maximum 150°, et son "complément", l'oscillation de 30°. Alors $K = K' \sqrt{3}$, comme on pourra le vérifier grâce aux vitesses.

Un autre cas particulier est celui de l'élongation 120° et son "complément" de 60°: la figure discrète exacte est très jolie à faire.

Il existe d'autre cas intéressants avec $K = K' \sqrt{n}$; mais cela fait entrer dans la discussion des fonctions elliptiques et de leurs relations algébriques: on pourra consulter le Greenhill (fonctions elliptiques ).

Voir aussi

• Pendule simple
• Démonstration du pendule simple discret
• Théorème de Poncelet sur les polygones inscrits et circonscrits
• Théorème de Steiner sur les colliers d'Apollonius
• Théorème du rebond en zigzag