Linearisation equivalente

Linearisation equivalente

Linéarisation équivalente

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

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) :

M \ddot{x} + B \dot{x} + K x = 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.

Sommaire

Principe

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 \varepsilon(t) = B \dot{x} + K x - g(x,\dot{x}) ,

On minimise alors le carré de sa moyenne quadratique \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 aléatoire gaussien. 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

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 :

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

\overline{x \dot{x}} représente la moitié de la dérivée par rapport au temps de \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 à

B = \overline{\dot{x} g} / \overline{\dot{x}^2} \quad K = \overline{x g} / \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

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 :

F_R = \sum_{i=0}^n {k_{2i+1} x^{2i+1}}

d'où

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 :

F_A = b_0 {sgn} \dot{x} + b_1 \dot{x} + b_2 |\dot{x}| \dot{x}

Le coefficient d'amortissement équivalent s'écrit

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

À l'équation

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

on substitue l'équation linéaire

M \ddot{x} + B \dot{x} + K x = F \cos \omega t

dont la solution est de la forme

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

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

I_m = \int_{0}^{\pi \over 2} \cos^m x\, dx

c'est-à-dire de fonctions Beta.

Tous calculs faits, on obtient la raideur

K = \sum_{i=0}^n {(2i+1) k_{2i+1} X^{2i} / 2^{2i}}

et le coefficient d'amortissement

B = {4\over\pi} {b_0\over{\omega X}} + b_1 + {8\over{3\pi}} b_2 \omega X

Solution du problème

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

((K(X) - M\omega^2)^2 + B(\omega X)^2 \omega^2 X^2 = F^2\,

Excitation gaussienne

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 :

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

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

M \ddot{X} + B \dot{X} + K X = F(t)

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

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

L'équation différentielle

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

a été résolue en utilisant

Linearisation equivalente.png
  • 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.

  • Portail des mathématiques Portail des mathématiques
  • Portail de la physique Portail de la physique
Ce document provient de « Lin%C3%A9arisation %C3%A9quivalente ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Linearisation equivalente de Wikipédia en français (auteurs)

Игры ⚽ Нужно сделать НИР?

Regardez d'autres dictionnaires:

  • Linéarisation Équivalente — 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 On la supposera symétrique : l excitation f(t) est centrée tandis que la fonction g est… …   Wikipédia en Français

  • Linéarisation équivalente — 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 On la supposera symétrique : l excitation f(t) est centrée tandis que la fonction g est… …   Wikipédia en Français

  • Liste des articles de mathematiques — Projet:Mathématiques/Liste des articles de mathématiques Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou probabilités et statistiques via l un des trois bandeaux suivants  …   Wikipédia en Français

  • Projet:Mathématiques/Liste des articles de mathématiques — Cette page n est plus mise à jour depuis l arrêt de DumZiBoT. Pour demander sa remise en service, faire une requête sur WP:RBOT Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou… …   Wikipédia en Français

  • Equation de Fokker-Planck — Équation de Fokker Planck L équation de Fokker Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d un processus de Markov. A l origine, une forme simplifiée de cette équation a… …   Wikipédia en Français

  • Processus de Gauss — Cette notion se rencontre dans des domaines variés allant des vibrations mécaniques aux vagues de la mer. De même que le théorème de la limite centrale permet de considérer une somme de variables aléatoires indépendantes comme une variable de… …   Wikipédia en Français

  • Processus de gauss — Cette notion se rencontre dans des domaines variés allant des vibrations mécaniques aux vagues de la mer. De même que le théorème de la limite centrale permet de considérer une somme de variables aléatoires indépendantes comme une variable de… …   Wikipédia en Français

  • Équation de Fokker-Planck — L équation de Fokker Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d un processus de Markov. A l origine, une forme simplifiée de cette équation a permis d étudier le… …   Wikipédia en Français

  • Équation de fokker-planck — L équation de Fokker Planck est une équation aux dérivées partielles linéaire que doit satisfaire la densité de probabilité de transition d un processus de Markov. A l origine, une forme simplifiée de cette équation a permis d étudier le… …   Wikipédia en Français

  • Determination des constantes d'equilibre — Détermination des constantes d équilibre Les constantes d équilibre sont évaluées pour quantifier les équilibres chimiques à partir de mesures de concentrations, directes ou indirectes, et mettant en œuvre des techniques numériques. Cet article… …   Wikipédia en Français

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”