Linéarisation équivalente — Wikipédia

De nombreux problèmes de dynamique basés sur la deuxième loi de Newton aboutissent à une équation du deuxième ordre que l'on peut mettre sous la forme

${\displaystyle M{\ddot {x}}+g(x,{\dot {x}})=f(t)}$

On la supposera symétrique : l'excitation f(t) est centrée tandis que la fonction g est impaire par rapport à ses deux paramètres position et vitesse.

Ce type d'équation ne possède, sauf exception, de solution explicite que lorsqu'elle est linéaire (voir Équation différentielle linéaire d'ordre deux) :

${\displaystyle M{\ddot {x}}+B{\dot {x}}+Kx=f(t)}$

Dans le cadre de l'hypothèse des petits mouvements il est possible d'obtenir une telle équation en négligeant simplement tous les termes d'ordre supérieur au premier. Si cette hypothèse est trop grossière, la linéarisation équivalente ou linéarisation optimale donne une meilleure approximation.

Principe[modifier | modifier le code]

La raideur K et le coefficient d'amortissement linéaire B étant a priori inconnus, l'idée consiste à exprimer l'erreur comme la différence entre les termes des deux équations ${\displaystyle \varepsilon (t)=B{\dot {x}}+Kx-g(x,{\dot {x}})}$,

On minimise alors le carré de sa moyenne quadratique ${\displaystyle {\overline {\varepsilon (t)^{2}}}}$ en annulant ses dérivées par rapport à B et K, ce qui fournit deux équations en B et K.

Ces équations permettent d'exprimer les deux inconnues en fonction des caractéristiques de la réponse. Par ailleurs, en utilisant l'équation linéarisée, le calcul de cette réponse en fonction de B et K permet d'obtenir finalement une équation algébrique qui donne une valeur généralement approchée de la solution.

Cette méthode s'utilise dans deux cas particuliers, celui de l'excitation sinusoïdale et celui de l'excitation par un processus de Gauss. Le calcul des moyennes est évidemment différent mais le principe de la méthode peut s'exposer dans les deux cas en représentant une moyenne par la grandeur surlignée, indépendamment du calcul exact de cette moyenne.

Équations en B et K[modifier | modifier le code]

C'est une méthode des moindres carrés tout à fait classique : on annule par rapport à B et K les dérivées de l'expression de l'erreur pour obtenir les deux équations :

${\displaystyle B{\overline {{\dot {x}}^{2}}}+K{\overline {x{\dot {x}}}}-{\overline {{\dot {x}}g}}=0}$

${\displaystyle B{\overline {x{\dot {x}}}}+K{\overline {x^{2}}}-{\overline {xg}}=0}$

${\displaystyle {\overline {x{\dot {x}}}}}$ représente la moitié de la dérivée par rapport au temps de ${\displaystyle {\overline {x^{2}}}}$ qui est nulle du fait de la symétrie supposée au début.

Les équations qui donnent les valeurs des coefficients en fonction de la réponse se réduisent donc à

${\displaystyle B={\overline {{\dot {x}}g}}/{\overline {{\dot {x}}^{2}}}\quad K={\overline {xg}}/{\overline {x^{2}}}}$

Pour finir de résoudre le problème il faut calculer la réponse de l'équation linéarisée en fonction des coefficients mais auparavant il est bon d'expliciter les forces couramment rencontrées.

Cas des forces de rappel et d'amortissement les plus courantes[modifier | modifier le code]

Une force de rappel s'exprime plus souvent par un développement limité que par une fonction mathématique plus élaborée. La fonction étant impaire, elle s'écrit donc :

${\displaystyle F_{R}=\sum _{i=0}^{n}{k_{2i+1}x^{2i+1}}}$

d'où

${\displaystyle K=\sum _{i=0}^{n}{{k_{2i+1}{\overline {x^{2i+2}}}} \over {\overline {x^{2}}}}}$

En ce qui concerne l'amortissement, il existe trois sortes d'amortissements classiques (voir Systèmes oscillants à un degré de liberté), le frottement sec, l'amortissement linéaire et l'amortissement quadratique. Pour la description de la méthode il est commode de les regrouper dans un seul terme bien qu'il soit peu vraisemblable qu'ils agissent simultanément. On obtient ainsi :

${\displaystyle F_{A}=b_{0}\ {\textrm {sgn}}\ {\dot {x}}+b_{1}{\dot {x}}+b_{2}|{\dot {x}}|{\dot {x}}}$

Le coefficient d'amortissement équivalent s'écrit

${\displaystyle B={{b_{0}{\overline {|{\dot {x}}|}}+b_{1}{\overline {|{\dot {x}}|^{2}}}+b_{2}{\overline {|{\dot {x}}|^{3}}}} \over {{\overline {|{\dot {x}}|}}^{2}}}}$

On peut remarquer que la raideur équivalente ne faisait également intervenir que les valeurs absolues.

Excitation sinusoïdale[modifier | modifier le code]

À l'équation

${\displaystyle M{\ddot {x}}+g(x,{\dot {x}})=F\cos \omega t}$

on substitue l'équation linéaire

${\displaystyle M{\ddot {x}}+B{\dot {x}}+Kx=F\cos \omega t}$

dont la solution est de la forme

${\displaystyle x=X\cos(\omega t+\varphi )\quad {\dot {x}}=-\omega X\sin(\omega t+\varphi )}$

Calcul de B et K en fonction de la réponse[modifier | modifier le code]

En reportant ces dernières expressions dans les formules qui définissent B et K et en tenant compte des symétries, le calcul de toutes les quantités se ramène à celui d'intégrales de la forme

${\displaystyle I_{m}=\int _{0}^{\pi \over 2}\cos ^{m}x\,dx}$

c'est-à-dire de fonctions bêta.

Tous calculs faits, on obtient la raideur

${\displaystyle K=\sum _{i=0}^{n}{(2i+1)k_{2i+1}X^{2i}/2^{2i}}}$

et le coefficient d'amortissement

${\displaystyle B={4 \over \pi }{b_{0} \over {\omega X}}+b_{1}+{8 \over {3\pi }}b_{2}\omega X}$

Solution du problème[modifier | modifier le code]

L'amplitude X de la réponse est donnée en fonction de F et de ω par l'équation

${\displaystyle ((K(X)-M\omega ^{2})^{2}+B(\omega X)^{2}\omega ^{2}X^{2}=F^{2}\,}$

Excitation gaussienne[modifier | modifier le code]

Dans ce cas, le second membre est un processus aléatoire F(t) représentant une infinité de signaux. La solution de l'équation est elle-même un processus X(t) constitué par une infinité de signaux :

${\displaystyle M{\ddot {X}}+g(X,{\dot {X}})=F(t)}$

Si le processus F(t) est gaussien, l'équation linéarisée

${\displaystyle M{\ddot {X}}+B{\dot {X}}+KX=F(t)}$

a pour solution un processus également gaussien de densité de probabilité

${\displaystyle p_{X{\dot {X}}}(x,{\dot {x}})={\frac {1}{2\pi \sigma _{X}\sigma _{\dot {X}}}}e^{-{\frac {x^{2}}{2\sigma _{X}^{2}}}-{\frac {{\dot {x}}^{2}}{2\sigma _{{\dot {X}}^{2}}}}}}$

Des calculs effectués en utilisant l'ergodicité pour remplacer les moyennes temporelles par les moyennes statistiques conduisent à des coefficients de raideur et d'amortissement équivalents que l'on peut calculer comme précédemment. La fin du calcul ressemble à celle du cas sinusoïdal mais nécessite en général une série d'itérations.

Si la non linéarité ne porte que sur la force de rappel et si l'excitation peut être assimilée à un bruit blanc, l'équation de Fokker-Planck donne une solution exacte, ce qui est généralement impossible avec l'excitation sinusoïdale. Dans le cas d'un amortissement non linéaire il n'est plus possible de trouver une solution exacte mais seulement une approximation meilleure que celle de la linéarisation.

Exemple[modifier | modifier le code]

L'équation différentielle

${\displaystyle {\ddot {x}}+{\dot {x}}+0.1x^{3}=F\cos(0.1t)}$

a été résolue en utilisant

• une linéarisation sommaire négligeant le terme cubique,
• la linéarisation équivalente,
• la technique numérique de Newmark.

Dans ce cas (une plus forte non-linéarité donnerait évidemment des résultats plus discutables), les résultats de la linéarisation se comparent très raisonnablement avec ceux de la simulation.

Bibliographie[modifier | modifier le code]

• (en) (en) Y. K. Lin, Probabilistic Theory of Structural Dynamics, New York, Robert E. Krieger Publishing Company, juillet 1976, 368 p. (ISBN 0882753770) pour le cas gaussien

• Portail de l'analyse
• Portail de la physique