Analyse des systèmes dynamiques
Analyse des systèmes dynamiques
1 - Pour un tel système, peut-on considérer que si Com(t) = constante alors (t) est aussi une constante ? Si non modifier le
schéma pour faire apparaître les causes d’une variations de (t).
2 - On considère que le bloc se décompose en trois blocs en série représentant successivement, le variateur, le radiateur
(élément chauffant seul) et la pièce. Tracer le nouveau schéma bloc en faisant apparaître les grandeurs intermédiaires entre ces blocs,
puis positionner les causes de la question 1 sur ce nouveau schéma.
3 - Le système de maintien est modifié en une régulation de température par l’installation d’un thermostat d’ambiance.
Donner le nouveau schéma bloc. Quel est l’intérêt fondamental de cette installation par rapport à celle décrite à la question 1.
Exercice n°3
1 - Déterminer la transformée de Laplace d’un créneau : x(t) = 0 pour t < 0 et T0 t
Exercice n°4
1 - Calculer la transformée de Laplace de la fonction définie par : x(t) = 0 pour t < 0
x(t) = X0 = V0 T0 pour T0 t
2 - Déterminer et tracer la réponse à l’entrée définie ci dessus, d’un système du premier ordre
K
( fonction de transfert H(p) = )
1+p
3 - Retrouver le résultat en utilisant le principe de superposition.
Exercice n°5 : Système du premier ordre
Un système d’entrée e(t) et de sortie s(t) nulles à l’origine est régi par l’équation différentielle
d e(t) d s(t)
3 e(t) + 2 = s(t) + 2
dt dt
1 - Déterminer la fonction de transfert du système.
2 - L’entrée e(t) est un échelon d’amplitude 0,5. Déterminer la sortie s(t) correspondante.
6 - Déterminer et représenter graphiquement le signal v(t) de transformée de Laplace V(p) qui est envoyé au système linéaire
de fonction de transfert G(p) dans la configuration avec boucle si x(t) est un échelon d’amplitude X0.
1 - Quelle est la valeur de la sortie au bout d’un temps ‘infini’ : y() = lim y(t) ?
t
2 - En raison du vieillissement, K diminue de 10%. De quelle proportion change y() ?
On insère ce système du premier ordre dans un
X(p) V(p) Y(p)
système asservi représenté ci-contre. La chaîne d’action est + A G(p)
-
composée d’un gain pur A et de G(p). Le retour est unitaire.
Y(p)
3 - Déterminer la fonction de transfert H(p) =
X(p)
4 - Quelle est la valeur du gain statique KH de H(p) ( rappel : KH = H(0) )
5 - Déterminer la valeur à donner au gain A pour que H(p) ait le même gain statique que G(p).
6 - Dans ces conditions, quel est l’effet de la diminution de 10% de K, sur la valeur y() de la réponse à x(t) = U0.u(t).
7 - Comparer à la question 2 et conclure.
Exercice n°11 ( Mise à part la question 1, attendre le cours de physique correspondant pour traiter cet exercice )
Le montage électrique ci-dessous est utilisé dans un asservissement. On e(t)
a T0
considère que l’entrée de ce montage est la tension e(t) et la sortie la tension s(t).
4 - Eliminer les variables inutiles pour faire apparaître la fonction de e(t) R1 = 1 s(t)
S(p)
transfert du montage G(p) = R2 = 1
E(p)
5 - Calculer et tracer la réponse s(t) à l’entrée e(t).
Exercice n°12 : D’après concours X-ENS Cachan 2005 Z(p)
X(p) W(p) Y(p)
+
+ K + A + B + C + D
- - - -
E
F
G
Z(p)
Déterminer les expressions H1 , H2 , H3 et H4 fonctions des
X(p) W(p) Y(p)
expressions K, A, B, C, D, E, F et G réalisant l’équivalence entre + H1 H2 ++ H3
-
le schéma blocs ci-dessus et le schéma blocs ci-contre.
H4
Performance en poursuite : Pour la question suivante, la route est plane, la perturbation Fp(t) est nulle. On agit sur l’entrée de
consigne yc(t), afin de faire varier la « garde au sol », c'est-à-dire la distance entre le plancher du véhicule et la route.
Y(p)
3 - Déterminer la fonction de transfert H1(p) = en l’absence de perturbation FP(p).
Yc(p)
Performance en régulation : Pour la question suivante, la consigne est nulle yc(t), la route impose une perturbation FP(t) non nulle.
Y(p)
4 - Déterminer la fonction de transfert H2(p) = en l’absence de consigne Yc(p).
FP(p)
Corrigé
Exercice n°1
.
1- En prenant la transformée de Laplace de l’équation Sg y(t) + Kd y(t) = Kd x(t),
Y(p) 1 Sg
on obtient Sg p Y(p) + Kd Y(p) = Kd X(p) soit = H(p) = avec τ =
X(p) 1+τp Kd
X0 1 X0
2- Si l’entrée est x(t) = X0 u(t), alors X(p) = d’où Y(p) = H(p) X(p) =
p 1+τp p
1 X0
On peut déjà calculer lim y(t) = lim p Y(p) = lim p = X0 → valeur finale X0
t→∞ p→0 p→0 1 + τp p
. 1 X0 X0
lim y(t) = lim p2 Y(p) = lim p2 = → tangente à l’origine atteint X0 en τ
t→0 p→∞ p→∞ 1 + τ p p τ
1 X0 A1 A2
Pour avoir y(t) pour tout t, il faut décomposer Y(p), c'est-à-dire écrire Y(p) = = +
1+τp p p 1+τp
En calculant p Y(p) en p = 0, on obtient A1 = X0
-1
En calculant ( 1 + τ p ) Y(p) en p = , on obtient A2 = - τ X0
τ
X0 τ X0 t
d’où Y(p) = - par transformée de Laplace inverse, on a y(t) = X0 [ 1 - exp( - ) ] u(t)
p 1+τp τ
V0 1 V0
3- Si l’entrée est x(t) = V0 t u(t), alors X(p) = d’où Y(p) = H(p) X(p) =
p2 1+τp p
2
1 V0
On peut déjà calculer lim y(t) = lim p Y(p) = lim p 2 =∞
t→∞ p→0 1+τp p
p→0
. 1 V0
lim y(t) = lim p2 Y(p) = lim p2 2 = V0 → direction asymptotique
t→∞ p→0 p→0 1+τp p
. 1 V0
lim y(t) = lim p2 Y(p) = lim p2 2 = 0 → tangente à l’origine horizontale
t→0 p→∞ p→∞ 1 + τp p
1 V0 A1 A A*
Pour avoir y(t) pour tout t, il faut décomposer Y(p), c'est-à-dire écrire Y(p) = 2 = + + 2
1+τp p 1+τp p p
-1
En calculant ( 1 + τ p ) Y(p) en p = , on obtient A1 = τ2 V0
τ
En calculant p2 Y(p) en p = 0, on obtient A* = V0
A1
En calculant p Y(p) en p = ∞, on obtient 0= +A soit A = - τ V0
τ
V0 τ V0 τ V0 2
t
d’où Y(p) = - + soit y(t) = V [ t - τ + τ exp( - ) ] u(t)
p2 p 1+τp τ
4- Comme les conditions initiales ne sont pas toutes nulles, on repart de l’équation différentielle
.
Sg y(t) + Kd y(t) = Kd x(t) dont la transformée de Laplace donne Sg ( pY(p) - y0 ) + Kd Y(p) = 0
Sg y0 τ y0 t
soit Y(p) = = par transformée de Laplace inverse, on a y(t) = y0 exp( - ) u(t)
Sg p + Kd 1 + τ p τ
X0 ω 1 X0 ω
5- Si l’entrée est x(t) = X0 sin( ωt ) u(t), alors X(p) = d’où Y(p) = H(p) X(p) =
p2+ ω2 1 + τ p p2 + ω2
.
On peut déjà calculer lim y(t) = lim p2 Y(p) = 0 → tangente à l’origine horizontale
t→0 p→∞
1 X0 ω A Bp C
Pour avoir y(t) pour tout t, il faut décomposer Y(p), c'est-à-dire écrire Y(p) = = + +
1 + τ p p2 + ω2 1 + τ p p2 + ω2 p2 + ω2
-1 X0 τ2 ω
En calculant ( 1 + τ p ) Y(p) en p = , on obtient A=
τ 1 + τ2 ω2
X0 ω
En calculant ( p2 + ω2 ) Y(p) en p = i ω , on obtient = B i ω+ C
1+τiω
X0 ω ( 1 - τ i ω ) X0 ω - X0 τ ω
soit = B i ω+ C qui conduit à C= et B=
1+τ ω 2 2
1+τ ω2 2
1 + τ2 ω2
X0 τ2 ω 1 X0 ω 1 X0 τ ω p
On obtient donc Y(p) = + -
1 + τ2 ω2 1 + τ p 1 + τ2 ω2 p2 + ω2 1 + τ2 ω2 p2 + ω2
X0 τ ω t X0 X0 τ ω
soit y(t) = exp( - ) u(t) + sin( ωt ) u(t) - cos( ωt ) u(t)
1 + τ2 ω2 τ 1 + τ2 ω2 1 + τ2 ω2
X0 τ ω t X0 1 τω
y(t) = exp( - ) u(t) + [ sin( ωt ) - cos( ωt ) ] u(t)
1 + τ2 ω2 τ 1 + τ2 ω2 1 + τ2 ω2 1 + τ2 ω2
X0 τ ω t X0
finalement y(t) = exp( - ) u(t) + sin ( ω t + ϕ ) u(t) avec tan ( ϕ ) = - τ ω
1 + τ2 ω2 τ 1 + τ2 ω2
6- Le système étant linéaire, on peut raisonner par superposition des cas 4 et 5, on obtient le cas 6.
réponse réponse
libre forcée
t X0 τ ω t X0
y(t) = [ y0 exp(- )+ exp(- ) + sin ( ω t + ϕ ) ] u(t) avec tan ( ϕ ) = - τ ω
τ 1 + τ2 ω2 τ 1 + τ2 ω2
réponse réponse
transitoire permanente
3- L’intérêt fondamental de cette installation, par rapport à celle décrite question 1, est la possibilité de s’opposer aux
perturbations. On peut remarquer que le bouclage sur le comparateur modifie la valeur informationnelle de l’entrée qui devient la
consigne de température. Ce système Perturbations
fonctionne comme un asservissement si la Cons θ I Q θ
Variateur Radiateur Pièce
consigne de température est plus élevée +
-
que la température extérieure, c’est à dire
si le système fonctionne en chauffage et Mes θ
Thermostat
non en climatisation.
Exercice n°3
1- Le signal x(t) demandé, en créneau d’amplitude X0 et de durée T0,
peut être considéré comme la somme de deux échelons d’amplitude X0 et X0
décalés de T0. L’un est X0 u(t), l’autre - X0 u(t-T0), T0
soit : x(t) = X0 [u(t) - u(t-T0)] t
X0 τ X0
donc Y(p) = K [ 1 - exp(-T0 p) ] [ - ]
p 1+τp
τ T0 τ+T0
1 τ exp(-T0 p) τ exp(-T0 p) t
soit Y(p) = K X0 [ - - + ]
p 1+τp p 1+τp
t t-T
soit par transformée inverse y(t) = K X0 [ u(t) - u(t-T0) - exp(- τ ) u(t) + exp(- τ 0 ) u(t-T0) ]
t
3- La réponse à X0 u(t) est K X0 [ u(t) - exp(- τ ) u(t) ]
t-T
La réponse à X0 u(t-T0) est K X0 [ u(t-T0) - exp(- τ 0 ) u(t-T0) ]
t t-T
La réponse à X0 u(t) - X0 u(t-T0) est K X0 [ u(t) - u(t-T0) - exp(- τ ) u(t) + exp(- τ 0 ) u(t-T0) ]
Exercice n°4
1- Le signal x(t) demandé, en rampe d’amplitude X0 pendant une
X0
durée T0, puis à valeur constante X0 peut être considéré comme la différence
de deux rampes d’amplitude X0 pendant une durée T0, et décalées de T0.
X0 X0
L’une est t u(t) , l’autre (t-t0) u(t-T0), T0 t
T0 T0
X0
soit : x(t) = [ t u(t) - (t-T0) u(t-T0) ]
T0 X0
Par transformée de Laplace, on en déduit
X0
L [x(t)] = X(p) = [ 1 - exp(-T0 p) ]
T0 p2
T0 t
K
2- Si le système est de fonction de transfert H(p) = , alors la sortie y(t) est de transformée de Laplace Y(p) avec :
1+τp
K X0 K X0 1
Y(p) = H(p) X(p) = [ 1 - exp(-T0 p) ] = [ 1 - exp(-T0 p) ] 2
1+τp T0 p2 T0 p (1 + τ p)
1 1 τ τ2
On sait que = 2 - + (voir exercice n°1, question 3)
p (1 + τ p) p p
2
1+τp
K X0 1 τ τ2
donc Y(p) = [ 1 - exp(-T0 p) ] [ 2 - + ]
T0 p p 1+τp
K X0 1 τ τ2 exp(-T0 p) τ exp(-T0 p) τ2 exp(-T0 p)
soit Y(p) = [ 2 - + - + - ]
T0 p p 1+τp p2 p 1+τp
K X0
soit par transformée inverse y(t) =
T0
[ t u(t) - τ u(t) + τ exp(- τt ) u(t) - (t-T0) u(t-T0) + τ u(t-T0) - τ exp(-
t-T0
τ
) u(t-T0) ]
X0 KX0
3- La réponse à t u(t)
T0
K X0 t
est
T0
[ t u(t) - τ u(t) + τ exp(- τ ) u(t) ]
X0
La réponse à (t-T0) u(t-T0)
T0 τ T0
K X0 t-T0 t
est
T0
[ (t-T0) u(t-T0) - τ u(t-T0) + τ exp(- τ ) u(t-T0) ]
X0
La réponse à [ t u(t) - (t-T0) u(t-T0) ]
T0
K X0
est
T0
[ t u(t) - τ u(t) + τ exp(- τt ) u(t) - (t-T0) u(t-T0) + τ u(t-T0) - τ exp(- t-Tτ 0 ) u(t-T0) ]
Exercice n°5 : Système du premier ordre
d e(t) d s(t)
1- 3 e(t) + 2
= s(t) + 2 donne en transformée de Laplace 3 E(p) + 2 p E(p) = S(p) + 2 p S(p)
dt dt
S(p) 3 + 2 p
soit la fonction de transfert H(p) = =
E(p) 1 + 2 p
0,5 3 + 2 p 0,5
2- L’entrée e(t) = 0,5 u(t) a une transformée de Laplace d’où la transformée de la sortie S(p) =
p 1+2p p
3 + 2 p 0,5 a b
Pour pouvoir prendre la transformée inverse de S(p), on décompose en éléments simples soit S(p) = = +
1+2p p p 1+2p
1,5 2
p S(p) en p = 0 conduit à 1,5 = a ( 1 + 2 p ) S(p) en p = -0,5 conduit à -2 = b soit S(p) = -
p 1+2p
donc s(t) = [ 1,5 - exp(-0,5 t ) ] u(t)
2- Le gain est 0,15 donc s(0-) = 2 . 0,15 = 0,3 La réponse à u(t) est 0,15 [ 1 - exp(-20 t) ] u(t)
Par superposition, la réponse à l’essai de lâcher 2 - u(t) est 0,3 - 0,15 [ 1 - exp(-20 t) ] u(t)
Z(p)
2-
2- Si K est multiplié par 0,9, alors y(∞) = K U0 est multiplié par 0,9 aussi.
A G(p) AK
3- H(p) = =
1 + A G(p) 1+AK+τp
AK
4- Le gain statique de H(p) est KH = H(0) =
1+AK
AK 1
5- Le gain statique de H(p) est K, soit à résoudre K= ⇒ A = = 5
1+AK 1-K
A K U0 5 . 0,8 U0 5 . 0,8 . 0,9 U0
6- Si K est multiplié par 0,9, alors y(∞) = passe de = 0,8 U0 à = 0,783 U0
1+AK 1 + 5 . 0,8 1 + 5 . 0,8 . 0,9
0,783
Ceci représente une multiplication par = 0,98 de la réponse finale y(∞).
0,8
7- Le vieillissement de K est moins sensible avec la structure en boucle.
Exercice n°11
1
1- e(t) = a [ t u(t) - 2 (t -T0) u(t -T0) + (t -2T0) u(t -2T0)] et E(p) = [ 1 - 2 exp(-T0p) + exp(-2T0 p) ]
p2
d
2- i1 R1 = e-s i2 = C (e - s) ( i1 + i2 ) R2 = s
dt
3- I1 R1 = E - S I2 = C p (E - S) (I1 + I2) R2 = S
R2 1 + C R1 p 1 1+p
4- G(p) = =
R1 + R2 C R1 R2 2 1
1+ p 1+ p
R1 + R2 2
1 1
5- s(t) =
2
[ t u(t) - 2 (t -T0) u(t -T0) + (t -2T0) u(t -2T0) ] +
2
[ u(t) - 2 u(t -T0) + u(t -2T) ]
1
-
2
[ exp(-2t) u(t) - 2 exp[-2(t -T0)] u(t -T0) + exp[-2(t -2T0)] u(t -2T0) ]
Exercice n°12 : D’après concours X-ENSCachan 2005
Le schéma initial de l’énoncé est équivalent à
Z(p)
X(p) W(p) Y(p)
C +
+ KA + + B + + D
- - - -
E
F
GA
Z(p)
X(p) W(p) Y(p)
C +
+ KA + B + + D
soit - - -
E
GA+F
Z(p)
X(p) W(p) Y(p)
C +
+ KAB + + D
- -
soit
B(GA+F)+E
Z(p)
X(p) W(p) Y(p)
C +
+ + KAB + D
- -
soit
B(GA+F)+E
KAB
Z(p)
X(p) W(p) Y(p)
C +
+ KAB + D
-
soit
G F E
1+ + +
K KA KAB
G F E
Finalement : H1 = K A B H2 = C H3 = D H4 = 1 + + +
K KA KAB
Remarque : Si dans cet exercice, on commence par supprimer la boucle interne pour transformer
Z(p) Z(p)
W(p) Y(p)
C + D W’(p) Y(p)
+ + C + D
- en + 1+CDE
On obtient un schéma dans lequel la grandeur W’(p) qui entre dans le bloc C n’est pas la grandeur voulue W(p).
Exercice n° 13 : Suspension hydraulique
d2 y(t) d y(t)
1- m = Kr [ x(t) - y(t) ] - f + FP(t) → m p2 Y(p) = Kr ( X(p) -Y(p) ) - f p Y(p) + FP
dt2 dt
d x(t)
q(t) = S → Q(p) = S p X(p)
dt
q(t) = Kd [ yc(t) - y(t) ] → Q(p) = Kd ( Yc(p) - Y(p) )
fp
3- FP(p) = 0
- par le calcul : m p2 Y(p) = Kr ( X(p) -Y(p) ) - f p Y(p)
Q(p) = S p X(p)
Q(p) = Kd ( Yc(p) - Y(p) )
Kr Kr
⇒ ( m p2 + f p + Kr ) Y(p) = Kr X(p) = Q(p) = K ( Yc(p) - Y(p) )
Sp Sp d
⇒ ( S m p3 + S f p2 + S Kr p + Kr Kd ) Y(p) = Kr Kd Yc(p)
Y(p) Kr Kd 1
⇒ H1(p) = = =
Yc(p) S m p3 + S f p2 + S Kr p + Kr Kd S Sf Sm 3
1+ p+ p2 + p
Kd Kd Kr Kd Kr
1
m p2 1
- par le schéma bloc en introduisant W(p) en amont du triple comparateur : Y(p) = W(p) = W(p)
fp f p + m p2
1+ 2
mp
Kr
f p + m p2 Kr
⇒ Y(p) = X(p) = X(p)
Kr Kr + f p + m p 2
1+ 2
fp+mp
Kd Kr
S p Kr + f p + m p 2 Kr Kd
⇒ Y(p) = Yc(p) = Y (p)
Kd Kr S m p 3 + S f p 2 + S Kr p + Kr Kd c
1+ 2
S p Kr + f p + m p
4- Yc(p) = 0
- par le calcul : m p2 Y(p) = Kr ( X(p) -Y(p) ) - f p Y(p) + FP
Q(p) = S p X(p)
Q(p) = - Kd Y(p)
Kr Kr
⇒ ( m p2 + f p + Kr ) Y(p) = Kr X(p) + FP = Q(p) + FP = - K Y(p) + FP
Sp Sp d
⇒ ( S m p3 + S f p2 + S Kr p + Kr Kd ) Y(p) = S p FP(p)
Y(p) Sp
⇒ H2(p) = =
FP(p) S m p3 + S f p2 + S Kr p + Kr Kd