X ENS Modélisation PSI 2006Sujet

Pas encore noté

Statique et dynamique des poutres (flexion, Euler-Bernoulli)
Équations aux dérivées partielles et analyse modale
Produits scalaires et orthogonalité de modes propres
Réduction des matrices symétriques
Systèmes différentiels linéaires du premier ordre
Stabilité, commandabilité et observabilité (automatique)

Explorer les annales

BanqueX-ENS Année2006 FilièrePSI MatièrePhysique Modelisation

Annales X-ENS PSI Physique Modelisation X-ENS PSI 2006 PSI Physique Modelisation Physique Modelisation 2006

Téléchargements

Télécharger le sujet (PDF)

Lire en ligne LaTeX Word

Corrigé : pas encore disponible
Rapport du jury : non disponible

Présentation du sujet

Contrôle des vibrations d'une poutre élancée : modélisation, analyse modale, discrétisation, commandabilité et observabilité

Le sujet aborde le contrôle des vibrations d'une structure élancée, modélisée par une poutre. La première partie établit le modèle de poutre d'Euler-Bernoulli à partir des équations d'équilibre et de la loi de Hooke. La deuxième étudie les vibrations libres et forcées de la poutre par analyse modale, avec les notions de modes propres et de résonance. La troisième propose deux discrétisations de la poutre (ressorts de traction ou de torsion) sous forme matricielle.

1Partie 1 : étude du modèle de poutreÉtablir la cinématique, les équations d'équilibre et la relation entre contraintes et efforts intérieurs d'une poutre en flexion, puis obtenir le modèle de poutre d'Euler-Bernoulli en statique et en dynamique.
2Partie 2 : étude d'une poutre en vibrationDéterminer par analyse modale les modes propres et pulsations propres d'une poutre d'Euler-Bernoulli encastrée, étudier leurs propriétés d'orthogonalité, puis utiliser ces modes pour résoudre le problème transitoire et le problème de vibrations forcées avec résonance.
3Partie 3 : discrétisation du problèmeDiscrétiser la poutre par deux modèles à ressorts (traction puis torsion), établir les équations de récurrence, étudier la convergence vers le modèle continu et écrire le problème sous forme matricielle avec ses modes propres.
4Partie 4 : observabilité et contrôlabilitéIntroduire la notion d'état, étudier la stabilité, la commandabilité et l'observabilité d'un système linéaire, puis construire un bouclage stabilisant et un observateur asymptotique pour la poutre discrétisée.

Ces sujets peuvent vous intéresser

Pas encore de corrigé pour ce sujet : voici des sujets proches corrigés.

Lecture du sujet en ligne

L'énoncé complet, avec les formules et les figures, sans ouvrir le PDF.

Sujet en PDF

MODÉLISATION EN SCIENCES PHYSIQUES ET SCIENCES DE L'INGÉNIEUR

Durée: 5 HEURES

Aucun document n'est autorisé.

Abstract

L'usage de calculatrices électroniques de poche à alimentation autonome, non imprimantes et sans document d'accompagnement, est autorisé, une seule à la fois étant admise sur la table ou le poste de travail, et aucun n'échange n'est autorisé entre les candidats.

Si , au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il le signale sur sa copie et poursuit sa composition en expliquant les raisons des initiatives qu'il est amené à prendre.

Le sujet comporte 16 pages

Ce sujet se propose de donner une première approche du problème du contrôle de vibration pour les structures élancées. Il s'agit, par exemple, d'éviter les résonances destructives pour un pont en construction dont une extrémité du tablier pend dans le vide. Pour cette première étude, on ramènera le problème à l'étude d'une poutre. La première partie aborde la modélisation d'un tel milieu, la seconde l'étude vibratoire, la troisième propose une discrétisation du problème de manière à aborder dans la quatrième partie la détection et le contrôle de l'état de la structure. Malgré la progression logique précédemment exposée, ces parties peuvent être traitées indépendamment.

Figure 1: Géométrie d'une poutre au repos

Notations générales

Soit une fonction $f(x, t)$ de la variable d'espace $x$ et de la variable de temps $t$ , on note $f˙$ la dérivée partielle par rapport au temps de $f : f˙ = (\partialf)/(\partialt)$
Pour une matrice (ou un vecteur) $M$ , on note $M^t$ la transposée de cette matrice.

1 Étude du modèle de poutre

1.1 Hypothèses, vocabulaire et notations

Géométrie

Une poutre est un milieu solide dont une dimension est très supérieure aux deux autres (on parle de structure élancée). On se place ici dans le cadre simple où, au repos, la poutre est un prisme droit à section rectangulaire. On munit la poutre d'un repère naturel cartésien (

O, x⃗, y⃗, z⃗

) où

O

est le centre de la face de gauche, on note (voir figure 1):

$L$ longueur de la poutre ( $x \in [0, L]$ );
$h$ demi-hauteur de la poutre $(y \in [ - h, h])$ ;
$b$ demi-profondeur de la poutre $(z \in [ - b, b])$ ;
$S$ aire des sections;
$S_x$ ensemble des points d'abscisse $x$ au repos ("section droite");
$G_x$ centre de $S_x$ ;
$Γ = {G_x, x ∈ [0, L]}$ "ligne moyenne".

Inertie

On suppose la poutre homogène sur chaque section et on note :

$ρ$ la masse linéique de la poutre ( $kg/m$ );
$J_z$ moment d'inertie linéique de la section $S_x$ autour de ( $G_x, z⃗$ ) (kg.m).

Cinématique (déformations)

On s'intéresse à un problème de flexion plane: la ligne moyenne se déforme dans le plan (

x⃗, y⃗

), et on suppose que les allongements selon

x⃗

sont négligeables (voir figure 2). L'hypothèse fondamentale est que chaque section

S_x

reste plane au cours de la transformation. On introduit alors pour chaque instant

t

(voir figure 3 ) :

$v(x, t)$ le déplacement de $G_x$ selon $y⃗$ (également appelé flèche);
$θ(x, t)$ la rotation selon $z⃗$ de la direction normale à la ligne moyenne par rapport à la verticale ;
$γ(x, t)$ la rotation selon $z⃗$ d'une section droite par rapport à la direction normale à la ligne moyenne.

Figure 2: Déformation d'une poutre

Figure 3: Déformation d'une poutre - schématisation et paramétrage

Efforts intérieurs

Les efforts intérieurs ou efforts de cohésion représentent les efforts qu'une partie de la matière exerce sur le reste de la poutre : si on suppose que l'on coupe la poutre selon une section, les efforts de cohésion en cette
section sont les efforts qu'il faudrait appliquer sur la partie gauche de la poutre pour que son état ne soit pas modifié. Étant donnée la cinématique de sections rigides, on peut modéliser les efforts de cohésion sur chaque section par un torseur

ζ(x, t).ζ(x, t)

représente donc les efforts appliqués sur la partie gauche de la section

S_x

par la partie droite. Les éléments de réduction de

ζ(x, t)

sont une résultante

T(x, t)y⃗

(effort tranchant) et un moment en

G_x

(moment fléchissant)

M(x, t)z⃗

(voir figure 4).

ζ(x, t) = {T(x, t)y⃗; M(x, t)z⃗}_(G_x)

Figure 4: Schématisation des efforts intérieurs sur la section

S_x

Efforts extérieurs

Pour solliciter la poutre en flexion, on considère un effort linéique réparti

p(x, t)y⃗

, appliqué sur l'ensemble de la poutre, et un éventuel effort concentré

F(t)y⃗

appliqué sur l'extrémité droite de la poutre (point

G_L

), voir figure 5. Notamment

p(x, t)

peut servir à modéliser le poids de la poutre.

Figure 5: Schématisation des efforts extérieurs répartis et concentrés

Petites perturbations

L'hypothèse des petites perturbations consiste à supposer les déplacements et les déformations suffisamment petits pour pouvoir raisonner sur la configuration au repos pour obtenir les équations d'équilibre, et pour pouvoir
assimiler la dérivation le long de la ligne moyenne à une dérivation par rapport à la variable

x

. En conséquence, on peut systématiquement considérer les angles

θ

γ

"petits" et linéariser les fonctions trigonométriques.

1.2 Préliminaires

Question 1.1 Donner les expressions de

S

J_z

en fonction de

h, b

ρ

.
Question 1.2 Donner la relation différentielle entre

θ(x, t)

v(x, t)

(on rappelle que

θ(x, t)≪1

1.3 Obtention des équations d'équilibre

Question 1.3 En appliquant le principe fondamental de la dynamique à un tronçon élémentaire soumis à un chargement linéique (voir figure 6), trouver deux équations aux dérivées partielles reliant

v(x, t), θ(x, t), γ(x, t), T(x, t)

M(x, t)

p(x, t)

Figure 6: Équilibre sous chargement réparti

Question 1.4 En considérant l'équilibre d'un tronçon élémentaire à l'extrémité droite de la poutre (voir figure 7), déterminer la valeur de l'effort tranchant et du moment fléchissant en

G_L : T(L, t)

M(L, t)

Figure 7: Équilibre sous chargement concentré à l'extrémité

1.4 Analyse locale des efforts

L'objectif de ce passage est d'expliciter les champs

T(x, t)

M(x, t)

définis sur chaque section par des contributions élémentaires réparties sur chaque section. On définit ainsi en chaque point

Q_x(y, z)

de la section

S_x

(on a

G_x Q_x(y, z)^(\to-) = yy⃗ + zz⃗

) un effort surfacique élémentaire (on parle de vecteur contrainte, voir figure 8):

d^2 f⃗_x(y, z, t) = (σ_x(y, z, t)x⃗ + τ_x(y, z, t)y⃗)d^2 S

ζ(x, t)

est le torseur associé à la répartition des efforts

d^2 f⃗_x(y, z, t)

sur la section

S_x

Figure 8: Vecteur contrainte

Question 1.5 On suppose que la contrainte normale est indépendante de z et varie linéairement en

y

σ_x(y, z, t) = yf_x(t)

. Montrer que la relation entre

σ_x(y, z, t)

M(x, t)

se met sous la forme

σ_x(y, z, t) = - y(M(x, t))/(I_z)

où

I_z

est un coefficient géométrique à déterminer.

Question 1.6 Des raisons de continuité et de symétrie conduisent à envisager la contrainte tangentielle sous la forme

τ_x(y, z, t) = (1 - (y/h)^2)g_x(t)

. Montrer que la relation entre

τ_x(y, z, t)

T(x, t)

se met sous la forme

τ_x(y, z, t) = (1 - (y/h)^2)(T(x, t))/(S~)

où

S~

est un coefficient géométrique à déterminer.

1.5 Relation avec les déformations

On donne la loi de l'élasticité linéaire (loi de Hooke) qui relie l'extension à la contrainte normale :

(dL)/(dx)(x, y, t) = (σ_x(y, t))/E

où

dL

traduit l'allongement d'un brin de longueur originale

dx

(voir figure 9 ).

E

est appelé module de Young.

Figure 9: Mise en évidence de l'allongement

Question 1.7 Trouver la relation entre

dL, dθ

dγ

. En déduire la relation entre

M(x, t), θ(x, t)

γ(x, t)

On admet la relation suivante sur les efforts tranchants,

μ

est appelé module de cisaillement.

T(x, t) = μS~γ(x, t)

1.6 Application

On considère le problème, indépendant du temps, d'une poutre encastrée sur sa face gauche (

v(0) = 0

θ(0) = 0

) soumise à son propre poids (figure 10).

Figure 10: Poutre pesante encastrée

Question 1.8 Déterminer les champs

T(x), M(x), γ(x), θ(x)

v(x)

.
Question 1.9 Évaluer

γ/θ(x)

, montrer que pour une configuration élancée (rapport

L/h

"grand", par exemple supérieur à 20),

γ(x)

devient négligeable devant

θ(x)

dès que l'on s'éloigne de l'encastrement (ie

x

pas trop près de 0 ). Pour les ordres de grandeur, on peut prendre

E \sim 2μ

1.7 Poutres d'Euler-Bernoulli

Ces poutres sont celles pour lesquelles on choisit de négliger

γ

. On a donc, dans les équations précédentes,

γ = 0.T

est alors obtenu par les équations d'équilibre (questions 1.3 et 1.4) et non plus par la relation (4).

On considère, toujours dans le cas statique (pas de dépendance en temps), une poutre d'Euler-Bernoulli non-pesante encastrée sur sa face gauche

(v(0) = 0

θ(0) = 0

) et soumise à un effort

F

à son extrémité droite (figure 11).

Figure 11: Poutre encastrée chargée à son extrémité

Question 1.10 Déterminer le champ v. Que vaut la flèche en bout de poutre

v(L)

Question 1.11 Placer le résultat de la poutre pesante obtenu précédemment dans le cadre des hypothèses actuelles, puis déterminer pour quel effort

F

on obtient la même flèche

v(L)

en bout de poutre.

1.8 Poutre d'Euler-Bernoulli en dynamique

Question 1.12 Mettre les équations obtenues aux questions 1.3, 1.4 et 1.7 dans le cadre des hypothèses d'Euler-Bernoulli (conserver les termes dynamiques). Se ramener à des équations différentielles faisant intervenir les inconnues

v

M

(penser à donner également les conditions aux limites).

Question 1.13 Montrer que, sous les mêmes hypothèses d'élancement, le terme d'inertie en rotation (terme en

J_z

) peut être négligé devant le terme d'inertie en translation (terme en

ρ

) (pour les calculs d'ordres de grandeur, on peut utiliser

(\partialv)/(\partialx) \sim v/L

) .

2 Étude d'une poutre en vibration

On admet dans cette partie qu'une poutre droite d'Euler Bernoulli encastrée sur sa face gauche et soumise à un effort de flexion réparti

p(x, t)

et un effort concentré

F(t)

en son extrémité droite, au repos à l'instant

t = 0

est gouvernée par un système différentiel de la forme :

ρ(\partial^2 v)/(\partialt^2)(x, t) + EI_z(\partial^4 v)/(\partialx^4)(x, t), = p(x, t); v(0, t), = 0; (\partialv)/(\partialx)(0, t), = 0; (\partial^2 v)/(\partialx^2)(L, t), = 0; (\partial^3 v)/(\partialx^3)(L, t), = - (F(t))/(EI_z); v(x, 0) = 0, et (\partialv)/(\partialt)(x, 0) = 0

2.1 Analyse modale - méthode directe

On s'intéresse pour l'instant à la recherche de solution du système (

5, 6, 7, 8, 9

) homogène (

p(x, t) = 0, F(t) = 0

), oscillant dans le temps à la pulsation

ω

, autrement dit en notation complexe

v(x, t) = e^(iωt)V(x)

(où

i

est solution de

i^2 = - 1, ω

est un réel positif).

Question 2.1 Quelle équation différentielle vérifie

V(x)

? Préciser les conditions aux limites.

Question 2.2 Montrer que pour obtenir des solutions en

V

non triviales (ie non nulles)

ω

doit vérifier une relation de la forme :

cos(λL)cosh(λL) + 1 = 0

avec

λ^4 = ω^2 ρ/(EI_z)

Question 2.3 Montrer schématiquement que cette relation admet une infinité dénombrable de solutions notées

λ_k

(desquelles on déduit

ω_k

). Exprimer les

V_k

associés sous la forme d'une fonction de (

λ_k, x

), montrer qu'ils sont définis à une constante multiplicative près.

Pour la suite on conservera la notation

V_k

sans chercher à employer leur expression analytique.

2.2 Étude des modes propres

On définit l'espace de fonctions

H

H = {u ∈ C^2([0, L]), u(0) = 0, (du)/(dx)(0) = 0}

On définit les formes bilinéaires suivantes pour

v

w

des fonctions de

H

(on note

v^'

la dérivée de

v

selon

x

) :

Forme bilinéaire de masse : $M(v, w) = \int_0^L ρvwdx$
Forme bilinéaire de rigidité : $K(v, w) = \int_0^L EI_z v^('')w^('')dx$

Question 2.4 Vérifier que

M

K

sont des produits scalaires sur

H

.
Question 2.5 Montrer que les modes propres (

V_i

) vérifient la relation suivante pour tout champ

u \in H

- ω_i^2 M(V_i, u) + K(V_i, u) = 0

Question 2.6 En déduire que pour

i \neq j

, sachant que

ω_i \neq ω_j

, les modes

V_i

V_j

sont orthogonaux au sens de

M

K

On introduit le quotient de Rayleigh

R

u \in H, R(u) = (K(u, u))/(M(u, u))

Question 2.7 Que vaut le quotient de Rayleigh d'un mode propre?
Par la suite on choisit de normer les

(V_k)

par la masse :

M(V_k, V_k) = 1

.
Question 2.8 Montrer que les modes propres rendent stationnaire le quotient de Rayleigh ou autrement dit montrer que la différentielle de

R

est nulle en

V_k

(on pourra s'intéresser à un champ légèrement "écarté" par rapport à un mode propre (

u = V_k + δ

) et étudier les variations de

R

autour de

V_k

2.3 Utilisation des modes propres pour une étude transitoire

On admet que n'importe quel champ de

H

peut s'écrire comme la combinaison linéaire de modes propres (éventuellement en nombre infini). On s'intéresse maintenant au problème non-homogène et on choisit de chercher une solution

v(x, t)

sous la forme suivante (le nombre de modes considérés, noté

m

, est un entier fixé a priori) :

v(x, t) = \sum_(i = 1)^m α_i(t)V_i(x)

Question 2.9 Par le même raisonnement qu'à la question 2.5 montrer que les fonctions

α_i

vérifient des équations différentielles découplées.

Question 2.10 Donner l'expression des

α_i

pour le cas d'une poutre au repos pour

t < 0

puis soumise à un effort concentré constant en bout de poutre.

2.4 Utilisation des modes propres pour une étude de vibrations forcées

On considère maintenant que la poutre est soumise uniquement à un effort linéique

p(x, t)

et que la dépendance en temps de

p

est une sinusoïde de pulsation

ω

p(x, t) = sin(ωt)P(x)

Question 2.11 Donner l'équation différentielle vérifiée par chaque coefficient

α_j

. Donner l'expression de la solution (sans chercher à déterminer les constantes d'intégration).

Question 2.12 Que se passe

t

-il quand

ω

tend vers une pulsation propre (

ω \to ω_k

)? Quel est le nom de ce phénomène, quelles peuvent en être les conséquences?

3 Discrétisation du problème

Dans le cas général d'une structure complexe, il est rare d'obtenir un modèle continu exploitable, on applique alors une démarche de discrétisation. De manière à contourner le problème de la discrétisation des équations aux dérivées partielles, on choisit ici de proposer deux modèles discrets distincts pour représenter, de manière approchée, la poutre. On étudie donc une poutre décomposée de deux manières :

par une succession de ressorts de traction verticaux de raideur $Δk$ reliant des segments rigides (appelés "éléments") horizontaux de longueur $ΔL$ (voir figure 12); le dernier élément est chargé par un effort concentré $F$ ,
par une succession de ressorts de torsion de raideur $Δr$ reliant des segments rigides (appelés "éléments") de longueur $ΔL$ (voir figure 13); le dernier élément est chargé par un moment de flexion $M_L = FL$ .

Figure 12: Première discrétisation : ressorts de traction

Figure 13: Deuxième discrétisation : ressorts de torsion

On note pour

i \in [ [0, N] ]x_i = iΔL

(avec

x_N = L

). Le segment

]x_(i - 1), x_i]

forme le

i^(eme)

élément. On appelle

v_i(t)

l'ordonnée du point d'abscisse

x_i

θ_i

l'angle formé par le

i^(eme)

élément de poutre par rapport à l'horizontale. On néglige les actions de la pesanteur.

3.1 Formalisme séquentiel

Question 3.1 Pour la première discrétisation, écrire l'équilibre d'un élément de poutre, en déduire une récurrence entre les

v_i

. Écrire l'équilibre du dernier élément.

Question 3.2 Dans le cadre statique (accélération négligeable), déterminer les

v_i(i \in [ [0, N] ])

Question 3.3 Pour la seconde discrétisation, calculer l'équilibre d'un élément de poutre, en déduire une récurrence entre les (

θ_i

). Écrire l'équilibre du dernier élément.

Question 3.4 Dans le cadre statique, déterminer les

(θ_i), i \in [ [1, N] ]

. On suppose

\foralli \in [ [1, N] ]θ_i≪1

, relier les (

θ_i

) aux (

v_i

) et trouver l'expression des

(v_i), i \in [ [0, N] ]

Les modèles discrets sont fondamentaux pour l'étude de modèles donnés par des équations aux dérivées partielles. Un critère essentiel pour les évaluer est leur limite quand on raffine la discrétisation.

Question 3.5 On s'intéresse au passage à la limite des deux discrétisations pour des problèmes statiques :

Vers quelle fonction de $x$ converge $(v_i)$ quand $N \to \infty$ pour la première discrétisation (pour des raisons d'homogénéité, on posera $Δk = (ES)/(ΔL)$ )?
Vers quelle fonction de $x$ converge $(v_i)$ quand $N \to \infty$ pour la deuxième discrétisation (pour des raisons d'homogénéité, on posera $Δr = (EI_z)/(ΔL)$ )?

Question 3.6 Pour chaque discrétisation représenter graphiquement l'allure de la solution pour

N = 2

éléments et la limite obtenue à la question précédente. Commenter (comparer éventuellement à la solution exacte obtenue à la question 1.10).

Question 3.7 Pour un problème statique, on s'intéresse succinctement à un chargement réparti constant par élément (qui se traduit dans l'équilibre d'un élément par un second membre

p_i

pour les ressorts de traction ou

m_i

pour les ressort de torsion). Pour les deux discrétisations, quelle est la régularité de la solution si (

p_i

) (respectivement (

m_i

)) tend vers une fonction continue ? si

(p_i)

(respectivement

(m_i)

) tend vers une fonction continue par morceaux? Plus que des démonstrations mathématiques, donner des interprétations physiques du problème, éventuellement illustrées.

3.2 Formalisme matriciel

L'équilibre dynamique d'une poutre discrétisée peut se mettre sous la forme :

MV¨ + KV = P

où

M

K

sont des matrices carrées

N \times N, V^t = (v_1, \dots, v_N)

P

est un vecteur représentant les efforts extérieurs (puisque

v_0 = 0

, il peut être éliminé).

Question 3.8 Donner l'expression de

M, K

P

pour la poutre discrétisée à l'aide de ressorts de traction, soumise à un chargement

F

à son extrémité droite.

On admet que la matrice

M

est symétrique définie positive et la matrice

K

est symétrique. Si on reprend la définition des modes propres ("solution harmoniques en temps du problème homogène") on trouve que la recherche de modes propres discrets se traduit par un problème aux valeurs propres généralisé :

- ω_i^2 MV_i + KV_i = 0

Question 3.9 Montrer que les hypothèses formulées précédemment sur

M

K

suffisent à l'existence de

N

valeurs et modes propres.

Question 3.10 Prouver que si

ω_i \neq ω_j

alors

V_j

V_i

sont orthogonaux au sens des matrices

M

K

On choisit alors systematiquement de normer les vecteurs propres par la masse (ie

V_i^t MV_i = 1

Pour la suite on prend des matrices de la forme :

M = m(1; ⋱; ⋱; 1), K = k(2, - 1; - 1, ⋱, ⋱; ⋱, ⋱, - 1; - 1, 2), P = F(0; ⋮; 0; 1)

où

m

k

sont des réels positifs homogènes respectivement à une masse et à une raideur (ces matrices ne correspondent pas aux discrétisations précédentes, mais conduisent à des calculs plus "légers"). On pourra noter

Ω = \sqrt(k/m)

.
Question 3.11 Pour une discrétisation à deux éléments, calculer les valeurs et vecteurs propres du système.

4 Notions d'observabilité et de contrôlabilité

On considère le système décrit par l'équation matricielle :

MV¨ + KV = P

où

M

est une matrice symétrique définie positive de dimension

N, K

est une matrice symétrique de dimension

N, M

K

ne dépendent pas du temps,

V

P

sont des vecteurs de dimension

N

dépendant du temps. On admet l'existence de

N

valeurs et vecteurs propres (

ω_i, V_i

) vérifiant

- ω_i^2 MV_i + KV_i = 0

On admet que ces vecteurs propres sont orthogonaux au sens de

M

K

.
Pour les applications on utilisera les matrices de l'équation (18) de dimension

N = 2

(discrétisation en deux éléments).

Cette partie cherche des réponses à deux questions :

Est-ce qu'un protocole permet de connaître l'état du système ?
Est-ce qu'un protocole permet de modifier l'état du système?

4.1 Notion d'état

Par état on entend non seulement la valeur des déplacements mais aussi la vitesse à un instant donné. On introduit donc la variable d'état

Φ = (V/(V˙))

, qui est donc un vecteur de dimension

2N

fonction du temps.
Question 4.1 Montrer que l'évolution de

Φ

est gouvernée par une équation différentielle du premier ordre que l'on écrira sous la forme

Φ˙ = AΦ + P~

Donner les expression de

A

P~

.
Question 4.2 Que devient cette équation quand on l'étudie dans la base des modes propres

V(t) = \sum_(k = 1)^N V_k α_k(t)

? Quelle est l'expression des matrices associées?

4.2 Systèmes linéaires commandés / observés

On considère ici que le chargement

P~

est une commande imposée par l'utilisateur de la forme

P~(t) = Bu(t)

où

u(t)

est une fonction scalaire et

B

un vecteur indépendant du temps. On considère également qu'un appareil de mesure schématisé par un vecteur ligne

D

permet de capter une information scalaire

ψ(t)

sur l'état du système. Autrement dit :

{Φ˙, =, AΦ + Bu; ψ, =, DΦ

Question 4.3 Donner l'expression de

Φ(t)

pour

t ⩾ 0

avec

Φ(0) = Φ_d

donné.

4.2.1 Stabilité

Un système est dit stable, si et seulement si pour toute condition initiale

Φ_d

et pour toute commande

u(t)

bornée, l'état

Φ(t)

reste borné.

Question 4.4 Montrer que le système est stable si et seulement si toutes les valeurs propres de A sont à partie réelle négative. On pourra utiliser la décomposition de l'espace en sous-espaces propres de

A

et étudier l'évolution du système sur chacun de ces sous-espaces.

Question 4.5 Le système de la poutre discrétisée en deux éléments est-il stable? Interpréter le résultat.

On note pour la suite le polynôme caractéristique de la matrice

A

sous la forme

Γ_A(λ) = det(A - λI) = λ^(2N) + a_1 λ^(2N - 1) + \dots + a_(2N)

4.2.2 Commandabilité

Le système (22) est dit commandable si, pour tout couple d'états (

Φ_d, Φ_a

), il existe un temps

T

fini, ainsi qu'une commande

u(t), t \in [0, T]

, qui appliquée au système est telle que (

Φ(0) = Φ_d

) implique (

Φ(T) = Φ_a

L'objectif ici est de démontrer que (22) est commandable si et seulement si la matrice de commandabilité

C = {BAB…A^(2N − 1)B}

est de rang

2N

.
Question 4.6 Montrer que (22) est commandable si et seulement si la matrice

Q_c(T) = \int_0^T e^(sA)BB^t e^(sA^t)ds

est définie pour au moins une valeur de

T

strictement positive.
Question 4.7 Montrer que si

T > 0, Q_c(T)

est définie si et seulement si

C

est de rang

2N

Question 4.8 Quelle est la matrice

B

associée à un effort imposé en bout de poutre? Montrer qu'un tel système est commandable.

De manière à élaborer un contrôleur le plus simplement possible, on se place dans une base adaptée (dite base canonique commandable).

Question 4.9 Montrer que, si le système (22) est commandable, il existe un changement de repère

P_c

tel que :

{A_c, = P_c AP_c^(− 1); B_c, = P_c B

avec

B_c = (0; ⋮; ⋮; 0; 1), et A_c = (0, 1, 0, \dots, 0; ⋮, ⋱, ⋱, ⋱, ⋮; ⋮, ⋱, ⋱, 0; 0, \dots, \dots, 0, 1; - a_(2N), \dots, \dots, - a_2, - a_1)

Pour ce faire, on pourra considérer la dernière ligne

L

de l'inverse de la matrice de commandabilité

C

et étudier des matrices du type:

(L; LA; \dots; LA^(2N - 1))

On appelle bouclage linéaire une commande pour (22) du type

u(t) = - RΦ(t) + v(t)

où le vecteur ligne

R

est appelé vecteur des gains de contre-réaction, et le scalaire

v

commande auxiliaire.

Question 4.10 Montrer, grâce au changement de base

P_c

introduit dans la section précédente, qu'il est toujours possible, par un choix judicieux de

R

, de stabiliser tout système du type (22).

Question 4.11 Proposer un bouclage linéaire

R

de sorte que le système de la poutre discrétisée en deux éléments soit auto-amorti. Autrement dit, chercher un bouclage

R

rendant la partie réelle de toutes les valeurs propres de A strictement négative. On pourra, par exemple, rendre les parties réelles toutes égales à un même nombre strictement négatif.

4.3 Observabilité

L'objectif de l'étude de l'observabilité est de déterminer dans quelle mesure il est possible de reconstituer l'état complet d'un système dont on n'observe qu'une partie.

Kalman a montré que cette notion est duale de la commandabilité. Il n'y a dès lors aucune surprise à ce que les questions suivantes (et leurs réponses) calquent les résultats établis dans la section précédente.
(22) est dit observable si pour tout

Φ_d

, il existe un temps

T

fini tel que si

Φ(0) = Φ_d

, alors la connaissance de

ψ(t)

sur

[0, T]

permet de calculer

Φ_d

L'objectif des deux prochaines questions est de démontrer que (22) est observable si et seulement si la matrice d'observabilité suivante est de rang

2N

O = {D; DA; ⋯; DA^(2N − 1)}}}{{

Question 4.12 Montrer que (22) est observable si et seulement si la matrice

Q_O(T) = \int_0^T e^(sA^t)D^t De^(sA)ds

est définie pour au moins une valeur de

T

strictement positive.
Question 4.13 Démontrer la condition d'observabilité (on pourra invoquer le résultat obtenu à la question 4.7).

Question 4.14 On choisit d'observer le déplacement en bout de poutre. Quelle est la matrice

D

associée? Montrer que le système est observable.

Question 4.15 Si la paire

(A, D)

est observable, montrer qu'il existe une matrice inversible

P_O

, telle que les matrices

A

D

, dans la base induite par

P_O

s'écrivent :

A_O, = P_O AP_O^(- 1); D_O, = DP_O^(- 1)

avec

A_O = (0, \dots, \dots, 0, - a_(2N); 1, ⋱, ⋮, ⋮; 0, ⋱, ⋱, ⋮, ⋮; ⋮, ⋱, ⋱, 0, - a_2; 0, \dots, 0, 1, - a_1) D_O = (0, \dots, 0, 1)

Ayant caractérisé l'observabilité d'un système linéaire, se pose la question de parvenir à construire un système équivalent au système complet, construit uniquement à partir d'observations partielles. Pour cela, on introduit un observateur

Φ^

Φ

, c'est-à-dire une quantité qui converge vers

Φ

quand la variable temps

t

tend vers l'infini.

On appelle observateur asymptotique linéaire un système de la forme :

{Φ^˙, =, AΦ^ + Bu − W(ψ^ − ψ); ψ^, =, DΦ^

tel que

Φ^

tend vers

Φ

lorsque

t \to \infty

, quelles que soient les conditions initiales

Φ(0)

Φ^(0)

Question 4.16 Montrer, en se plaçant dans la base induite par

P_O

, que tout système linéaire commandé-observé admet un observateur asymptotique

W

. On introduira l'erreur de l'observateur :

e = Φ^ - Φ

Question 4.17 Proposer un observateur asymptotique

W

pour la poutre dont on observe le déplacement à l'extrémité.

4.4 Bilan commandabilité / observabilité

Grâce aux résultats démontrés dans les deux sections précédentes, on sait déterminer si un système est commandable et observable. Si tel est le cas, on sait construire un contrôleur du système, ainsi qu'un observateur. L'idée naturelle qui vient, est d'utiliser l'observateur, qui reconstitue l'état du système à partir de mesures partielles, dans le contrôleur du système complet, qui est capable de stabiliser le système vers toute commande prédéfinie.

Cette démarche s'appelle le principe de séparation estimation-commande. On peut montrer que le contrôleur muni de l'observateur, stabilise le système linéaire commandé-observable, quelles que soient les conditions initiales de

Φ

et de

Φ^

Questions fréquentes

4 questions

Sur quels chapitres porte ce sujet de modélisation X-ENS PSI 2006 ?

Il porte sur la mécanique des poutres et leur mise en équation, l'analyse modale des vibrations, la discrétisation matricielle et les notions d'automatique (stabilité, commandabilité, observabilité).

Le sujet demande-t-il des connaissances en automatique ?

Oui, la dernière partie introduit les notions de stabilité, de commandabilité et d'observabilité d'un système linéaire, et construit un contrôleur et un observateur pour la poutre discrétisée.

Les parties de ce sujet sont-elles indépendantes ?

L'énoncé précise que, malgré une progression logique, les quatre parties peuvent être traitées indépendamment.

Ce sujet est-il accessible en première année ?

Non, il mobilise des outils de deuxième année (équations aux dérivées partielles, réduction des matrices symétriques, systèmes différentiels linéaires) appliqués à un problème de mécanique des poutres et d'automatique.

Pas de description pour le moment

Commentaires• X ENS Modélisation PSI 2006

Connectez-vous pour participer aux discussions

Partagez vos avis, posez des questions et échangez avec la communauté

Chargement des commentaires...