Une équation différentielle classique en Mécanique - Exercice 1
45 min
70
Il existe des équations différentielles "modèles" qui interviennent dans de très nombreux domaines. Cela signifie qu'il existe des comportements similaires dans de nombreux domaines des Sciences. En voici un exemple.
Question 1
On note par t le temps, qui est le paramètre d'évolution en Physique classique. On a alors t⩾0. On considère l'équation différentielle suivante, qui modélise le comportement oscillatoire d'une quantité f qui caractérise un système physique réel : dt2d2f(t)+2λdtdf(t)+ω02f(t)=Acos(ωt) Dans cette équation différentielle linéaire, du second ordre, à coefficients constants, on a : ∙ le terme λ caractérise l'amortissement temporel de f, et est un réel strictement positif ; ∙ le terme ω0 caractérise l'oscillation intrinsèque (on dit aussi propre) de f, et est un réel strictement positif ; ∙ le terme A caractérise l'amplitude de l'excitation externe du système physique, et est un réel strictement positif ; ∙ le terme ω caractérise l'oscillation externe du système (on dit aussi que c'est celle d'un opérateur extérieur), et est également un réel strictement positif. Pour des raisons purement d'ordre physique, on admettra que laconditionsuivanteesttoujoursreˊaliseˊe lors de cet exercice : λ<ω0 Les deux {\color{blue}{\textbf{conditions initiales}}}, qui définissent le problème de \textit{Cauchy} considéré, sont les suivantes : ⎩⎨⎧f(t=0)dtdf(t=0)==00 On adoptera la notation simplificatrice suivante : Ω=ω02−λ2>0
Déterminer le discriminant Δ associée à l'équation caractéristique de l'équation différentielle, et précisez son signe.
Correction
Le discriminant Δ, associée à l'équation caractéristique de l'équation différentielle, est : Δ=(2λ)2−4×1×ω02=4λ2−4ω02=4(λ2−ω02) Soit encore, avec l'hypothèse 0<λ<ω0, l'expression suivante : Δ=4(λ2−ω02)<0
Question 2
En déduire les expressions des racines de cette équation caractéristique.
Correction
Les expressions des deux racines r1 et r2, complexes conjuguées l'une de l'autre, de cette équation caractéristique sont (avec i2=−1) : ⎩⎨⎧r1r2==2−2λ+i−4(λ2−ω02)2−2λ−i−4(λ2−ω02) En simplifiant : ⎩⎨⎧r1r2==−λ+iω02−λ2−λ−iω02−λ2 Soit : ⎩⎨⎧r1r2==−λ+iΩ−λ−iΩ
Question 3
En déduire l'expression de la solution homogène, notée fH(t), de l'équation différentielle proposée.
Correction
L'expression de la solution homogène, notée fH(t), de l'équation différentielle proposée est : fH(t)=e−λt(Acos(Ωt)+Bsin(Ωt))avec:(A;B)∈R2
Question 4
La solution particulière recherchée, notée fP(t), de l'équation différentielle proposée, sera choisie de la forme suivante : fP(t)=αcos(ωt)+βsin(ωt) Avec α et β qui sont deux nombres réels. Déterminer l'expression de fP(t) associée à notre problème, c'est-à-dire les expressions des deux coefficients réels α et β.
Correction
La solution particulière recherchée, notée fP(t), de l'équation différentielle proposée, sera choisie de la forme suivante : fP(t)=αcos(ωt)+βsin(ωt) Avec α et β qui sont deux nombres réels. Ainsi : dtdfP(t)=fP′(t)=−αωsin(ωt)+βωcos(ωt) Et aussi : dt2d2fP(t)=fP′′(t)=(fP′(t))′=−αω2cos(ωt)−βω2sin(ωt) Donc, en remplaçant dans l'équation différentielle originelle, on obtient : −αω2cos(ωt)−βω2sin(ωt)+2λ(−αωsin(ωt)+βωcos(ωt))+ω02(αcos(ωt)+βsin(ωt))=Acos(ωt) Soit encore : (−αω2+2λβω+αω02)cos(ωt)+(−βω2−2λαω+βω02)sin(ωt)=Acos(ωt)+0sin(ωt) D'où : (α(ω02−ω2)+2λβω)cos(ωt)+(β(ω02−ω2)−2λαω)sin(ωt)=<br/>Acos(ωt)+0sin(ωt) On obtient alors le système suivant : ⎩⎨⎧α(ω02−ω2)+2λβωβ(ω02−ω2)−2λαω==A0 Ainsi, on en déduit que : α=2λω(ω02−ω2)β On a alors : 2λω(ω02−ω2)β(ω02−ω2)+2λβω=A Soit encore : β(2λω(ω02−ω2)2+2λω)=A Ce qui nous donne : β(2λω(ω02−ω2)2+4λ2ω2)=A Ainsi : β=(ω02−ω2)2+4λ2ω22λωA On en déduit alors que : α=2λω(ω02−ω2)(ω02−ω2)2+4λ2ω22λωA En simplifiant: α=(ω02−ω2)2+4λ2ω2(ω02−ω2)A Dès lors, l'expression de fP(t) devient : fP(t)=(ω02−ω2)2+4λ2ω2(ω02−ω2)Acos(ωt)+(ω02−ω2)2+4λ2ω22λωAsin(ωt) En factorisant : fP(t)=(ω02−ω2)2+4λ2ω2A((ω02−ω2)cos(ωt)+2λωsin(ωt))
Question 5
Donner alors l'expression de la solution mathématique globale f(t) ainsi recherchée.
Correction
L'expression de la solution mathématique globale f(t) est donc donnée par : f(t)=fH(t)+fP(t) A savoir, avec (A;B)∈R2 : f(t)=e−λt(Acos(Ωt)+Bsin(Ωt))+(ω02−ω2)2+4λ2ω2A((ω02−ω2)cos(ωt)+2λωsin(ωt))
Question 6
A l'aide des deux conditionsinitiales proposées, déterminer l'expression de la solution physique f(t) qui décrit le comportement oscillant réel de f.
Correction
Les deux conditions initiales proposées sont : ⎩⎨⎧f(t=0)dtdf(t=0)==00 Or, on a : f(t=0)=e−λ0(Acos(Ω0)+Bsin(Ω0))+(ω02−ω2)2+4λ2ω2A((ω02−ω2)cos(ω0)+2λωsin(ω0)) Soit encore : f(t=0)=A+(ω02−ω2)2+4λ2ω2A(ω02−ω2) La première des deux conditions initiales nous donne : A+(ω02−ω2)2+4λ2ω2A(ω02−ω2)=0 Ce qui nous conduit à l'expression de A : A=−(ω02−ω2)2+4λ2ω2A(ω02−ω2) Maintenant, déterminons l'expression de la dérivée première temporelle f′. On a alors : f′(t)=+−λe−λt(Acos(Ωt)+Bsin(Ωt))+Ωe−λt(−Asin(Ωt)+Bcos(Ωt))(ω02−ω2)2+4λ2ω2Aω(−(ω02−ω2)sin(ωt)+2λωcos(ωt)) Posons maintenant t=0, on obtient alors : f′(t=0)=+−λe−λ0(Acos(Ω0)+Bsin(Ω0))+Ωe−λ0(−Asin(Ω0)+Bcos(Ω0))(ω02−ω2)2+4λ2ω2Aω(−(ω02−ω2)sin(ω0)+2λωcos(ω0)) Soit encore : f′(t=0)=−λA+ΩB+(ω02−ω2)2+4λ2ω2Aω(2λω) Mais également : f′(t=0)=−λA+ΩB+(ω02−ω2)2+4λ2ω22λω2A Ainsi, la seconde condition initiale nous donne : −λA+ΩB+(ω02−ω2)2+4λ2ω22λω2A=0 Ce qui nous permet d'écrire que : ΩB=λA−(ω02−ω2)2+4λ2ω22λω2A Soit : ΩB=−(ω02−ω2)2+4λ2ω2Aλ(ω02−ω2)−(ω02−ω2)2+4λ2ω22λω2A Donc : ΩB=−(ω02−ω2)2+4λ2ω2Aλ((ω02−ω2)+2ω2) Ce qui nous donne : ΩB=−(ω02−ω2)2+4λ2ω2Aλ(ω02−ω2+2ω2) On aboutit alors à : ΩB=−(ω02−ω2)2+4λ2ω2Aλ(ω02+ω2) Finalement : B=−Ω((ω02−ω2)2+4λ2ω2)Aλ(ω02+ω2) Donc, la solution f(t) qui satisfait à notre problème de Cauchy est donc donnée par : f(t)=+e−λt(−(ω02−ω2)2+4λ2ω2A(ω02−ω2)cos(Ωt)−Ω((ω02−ω2)2+4λ2ω2)Aλ(ω02+ω2)sin(Ωt))(ω02−ω2)2+4λ2ω2A((ω02−ω2)cos(ωt)<br/>+2λωsin(ωt)) En factorisant, on obtient : f(t)=+−(ω02−ω2)2+4λ2ω2Ae−λt((ω02−ω2)cos(Ωt)−Ωλ(ω02+ω2)sin(Ωt))(ω02−ω2)2+4λ2ω2A((ω02−ω2)cos(ωt)<br/>+2λωsin(ωt)) Finalement, on trouve que : f(t)=+(ω02−ω2)2+4λ2ω2A(−e−λt(<br/>(ω02−ω2)cos(Ωt)−Ωλ(ω02+ω2)sin(Ωt))(ω02−ω2)cos(ωt)+2λωsin(ωt))
Question 7
On note par T une quantité qui à les dimension d'un temps, et qui caractérise le comportement oscillatoire de f. Cette quantité T s'appelle la " pseudo-peˊriode de f. Est-il possible, et pourquoi, de définir T par la relation T=Ω2π ?
Correction
Il est impossible de définir T par la relation T=Ω2π car la solution trouvée précédemment présente des termes trigonométriques d'arc Ωt mais aussi d'arc ωt. Ces deux arcs sont, à priori, différents, d'où l'impossibilité de la formule proposée pour T.
Question 8
Parmi les trois courbes proposées sur le graphique ci-dessous, pour t∈[0;8], laquelle est selon vous celle qui correspond à un comportement physique qui peut être modélisé par les conditions de notre exercice ? Indiquer clairement sa couleur : rouge, bleue ou vert. Bien évidemment, votre choix sera accompagné d'une justification claire, argumentée, et sans ambiguïté.
Correction
Seule la solution graphique de couleur bleue bleue satisfait aux deux conditions initiales. Donc la bonne solution est la courbe de couleur bleue.