TD Regul 2005
TD Regul 2005
TD Regul 2005
1 Transformée de Laplace
1.1 Calculs de transformées inverses
Retrouver le signal temporel à partir de la transformée de Laplace Y (p) :
1.
p2 + p + 9
Y (p) =
p(p2 + 9)
2.
4p + 12
Y (p) = 2
p + 8p + 16
3.
−3p2 + 6p − 24
Y (p) =
(p − 1)(p − 2)2
4.
e−p (p2 + 6p + 13)
Y (p) =
p(p + 3)2
5.
5p2 − 15p − 11
Y (p) =
(p + 1)2 (p − 2)2
∗
Pour toutes questions, e-mail : {msekher,eric.magarotto}@greyc.ismra.fr
1.2 Résolution d'équations diérentielles par la transformée
de Laplace
Résoudre les équations diérentielles suivantes :
1.
ÿ − 4ẏ + 8y = 0
ẏ(0) = 3
y(0) = 7
2. ½
ẏ − 2y = sin(t)
y(0) = 1
3.
ÿ + y = t
ẏ(0) = 2
y(0) = 1
4.
ÿ + 3ẏ − 4y = e−2t
ẏ(0) = 1
y(0) = 1
T1 T2
2 Modélisation
Illustration de la notion de fonction de transfert avec résistance et ressort.
Qp
On note :
S : la section du bac (en m2 ) ;
h(t) : la hauteur du bac (en m) ;
Qe (t) : le débit d'alimentation (en m3 /s) ;
Qp (t) : le débit de vidange (en m3 /s).
Nous cherchons à modéliser le phénomène de remplissage et de vidange du bac.
Questions :
1. Ecrire le bilan de masse (ou conservation de la matière) faisant intervenir h,
Qe et Qp et S .
2. En première approximation, nous pouvons écrire Qp = αh. Que devient la
relation précédente ?
3. En utilisant la transformée de Laplace, donner la fonction de transfert du
système. Quel est l'ordre du système ? Tracer les pôles dans le plan complexe.
Q, Ce(t)
V
Q, Cs(t)
Fig. 2 Réacteur
On note :
Ce (t) : la concentration initiale du produit (mole/m3 ) ;
Cs (t) : la concentration nale du produit (mole/m3 ) ;
Q : le débit du produit (en m3 /s) ;
V : le volume du réacteur (en m3 ).
Questions :
1. Ecrire le bilan de masse faisant intervenir Q, Ce et Cs et V .
2. Notons F (p) la fonction de transfert liant Cs (p) à Ce (p), c'est à dire Cs (p) =
F (p)Ce (p). En utilisant la transformée de Laplace, exprimer F (p).
3. Quel est l'ordre du système ? Tracer les pôles dans le plan complexe.
2.3 Amortisseur
Un amortisseur peut être considéré comme un système masse-ressort décrit dans
la gure 3. On cherche à modéliser ce système.
On note :
K la constante de raideur du ressort (en N/m) ;
f la constante de frottement visqueux du vérin (en N s/m) ;
u(t) la force exercée sur la masse (en N ) : c'est l'entrée du système ;
y(t) le déplacement engendré par u(t) (en m) : c'est la sortie du système.
Questions :
1. Appliquer la relation fondamentale de la dynamique (pomme des forces =
masse × accélération).
2. Exprimer la fonction de transfert du système.
3. Quel est l'ordre du système ? Tracer les pôles dans le plan complexe.
TD n◦ 2 : schémas-blocs ∗
I.U.T. Génie Chimique & Génie des Procédés - 2ème année
1 Schéma-bloc physique
Pour l'exemple qui suit, identier tous les éléments typiques d'une boucle (à savoir le
processus, le capteur, l'actionneur et le correcteur). Déterminer aussi quelle est la consigne.
Exemple : prenons le cas d'un réfrigérateur : ce système consiste à refroidir l'air grâce à
un moteur. Un thermomètre indique la température, celle-ci est comparée à la température
désirée. Un circuit électrique met en route le moteur si la température est supérieure à celle
demandée ; sinon le circuit électrique arrête le moteur.
2 Réduction de schémas-blocs
Après simplication du schéma-bloc, calculer la fonction de transfert F (p) liant X(p)
à Y (p), c'est-à-dire Y (p) = F (p)X(p).
∗
pour toutes questions, e-mails : {msekher, eric.magarotto}@greyc.ismra.fr
1.
H(p) ¾
2.
X(p) +- 10 1 Y (p)
p+1 p -
− 6
2
+ ?+
¾
3.
X(p) +- 1 2 1 Y (p)
K p + 10 p+3 p -
− 6
0.1
+ ?+ + ?+
¾ ¾
4.
H1 (p)
−
?
X(p) +
- -
+ - G1 (p)
+
- - G2 (p) - G3 (p)
Y (p)
-
−6 −6
H2 (p)
B
D
OP N
∗
X X W X X X W
\ f _ \
Amplitude :; 9
OP N
xy w
Q
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
@
0
R
r r
Z g
^
/ . -
1 9 ]
? Y TU S
<
b
=
G =
G
F2 (p)
Q
K@
z I
R
A
I @ ? >
8 D
B@
K
x w
TU S
`
K
B@
B c
=
0.5
@
QV :;
s D F
u =
{| z
>
1 0
B@ A
=
:
A
A^
QV 3 =
?
>
<
=
1
}x
B
B@
H
I
>
D
B B
3d
D
{| A
B@ B A
@
/
tr
X 02
K@
A
F
B A
?
C <
z . BZ A
e
C <
B
1.5
;
w
|
<B;
<
=
D
B
<
51
34
D
F D
B@ >
B
;
F E
/
H F E
2
=
D
B
B
A
@
B A
51
<
A
K < <
I D
?
34
H
2.5
G ?
>
>
B
K >
I
@
?
@
Time (sec)
H
H :;
Step Response
[ =
:; 6
J
@
?
H
1
7
:; K
G
3
=
◦
51
H D
B
<B
:;
58
;
D
<B @
:; =
=
3.5
;
= 6
B@
K
<B
1 B
7
o mn l
G
>
A
5
:@
D
G
I
>
> `
B
4
0
58
=
B =
A :;
@
@
G
51
I
=
B G
A
D >
p W
B A
;
B
:;
4.5
B A
I 5 X ;
95%
B
:; =
t
K;
0
<
? H
<
q
H
B@
J
51
D
;
@
5
:;
=
X D
W
>
D B
@
D
B
I
=
;
Amplitude
<B
D
B
;
;
>
B
"
= I
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
2.2
2.4
2.6
2.8
3
G
D
0
h
D
B
G
: <
>
F B@
1
F
$% #
B =
2
<
B A B
<K
&
τ
B;
u(t) = 1
3
L =
!
;
"
>
@ K
>
4
B
D
D
M B
C <
B '
5
>
B
j K
<
2
D )* (
6
F E
C <
7
X
temps (sec.)
<
D
réponse indicielle
∗
8
F (p) =
F E
9
<
3
10
11
1 + 10p
.
F1 (p)
X
u(t) = 2 , t ≥ 0
12
Z
13
IKG
14
k
B
B@
[ f
B
Amplitude Amplitude Amplitude
K
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
100
150
200
250
300
350
400
450
50
<
Z
0
0
g
F3 (p) =
7
B A
K
Z
2
2
B
1
K
1+1p
u <
@
3
BH
B
4
4
G @
:@
<B
6
6
5/
H
F4 (p) =
;
< D
B
:;
8
8
34
D
0.5
B A
1+3p
>
Time (sec)
Time (sec)
Time (sec)
@
Step Response
Step Response
10
10
Step Response
<
10
51
:;
6
12
12
12
1
7
C <
<B
58
14
14
<
14
<
`
F E
K;
B
16
16
16
G ?
51
>
k
B
18
18
18
:;
:;
G d
A
:; =
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
>
0
0
0
IF
>
B
2
2
2
4
@
D
D
K
; 8
I
4
4
4
>
B
<
B@
3
2
M 1
6
6
6
>
A
8
8
8
<
0
=
B
K@
Time (sec)
Time (sec)
Time (sec)
c
Step Response
Step Response
Step Response
10
10
10
7
F1 (p) =
<
C <
1
12
12
12
0.5
<
1+1p
u
D
F E
14
14
14
X
<
16
16
16
?
Y
F2 (p) =
D
B@
F
18
18
18
1
s
1+3p
u
7
.
3
5/
34
51
58
51
3
31
1
g
>
A
f
=
?
?
D
D
B
<
K@
<
<B
<
:;
<B
<
K;
<
:;
<
<
<
<
<
B@
B@
B
C
X
0.5 1
F1 (p) = F2 (p) = F3 (p) =
g
>
>
>
J
=
`
E
D
D
Z
K@
:@
<
:;
<
<
:;
:;
<
B@
C
M
1+1p 1+3p
u
1 0.5
1+1p F4 (p) = 1+3p
u
0.9
120
0.8
100 0.7
0.6
80
Amplitude
Amplitude
0.5
60
0.4
0.3
40
0.2
20
0.1
0 0
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Time (sec) Time (sec)
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
Amplitude
Amplitude
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Time (sec) Time (sec)
0.9 0.9
0.8 0.8
0.7 0.7
0.6 0.6
Amplitude
Amplitude
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0
0 2 4 6 8 10 12 14 16 18 0 2 4 6 8 10 12 14 16 18
Time (sec) Time (sec)
_
T A
G
B
+
H @ ?
cd b cd b O
GH
∗
I
K
^ ^ @
[ n [ m [ [ n [ m [
u(t)
r
@A ?
L
B@
C
{| z
e e
H EZ I
V
s
f f L
H H H H
w H
L
q - ,
v hi g hi g
k k k k
E O
D
H
: B@C
B@ H
I I I I
} L K @
EZ
H
H G M E H
{ z @
C
U U
t≥0
./
a
H O
K
C C Y
G
D
ej ej
D U
~ }
TC
< H
D
B@ GH F E
C
H
G
G t
O
a O K
A
K
H
H H
3
B
0
K K
K
I
H H
C H I
U K
{
C
U UZ
E .-
K U
C
C O
~ E
G U H
U D
C <=
T
L
O I@
U
H L C
M
G 2 1
}
U W
GH
C
T
A
H
A
.
AO
@ D
I C O
L
J
z
C
G
@ D
W
H
B@ D
5
C
O J
O
y(∞)
u(∞)
L F
K
H
@
H
L
34
4 T
D
H K
K M
l L K F
G K M
O
H
E
H
K H
2 K
O
H
F
3
O
U
D E
B@
GH
ξ
U K
-
A - C
X D
G
K
@
φ = arcsin
H
H
K M
E
A
X D
M
tm
U
E
G A
H
M
M E
G
D
x
I
T G
5
GH M
C
I
o
64
Y
M
T
D
³p
L
<
T
K
◦
B@
CH
G D
H A B
T
@ K
Y O
H C
"
G
ω0
H
B@ T
E C
H / G
N
7-
O
ω0 ξ
U
V
O
B
C
U
H
1 − ξ2
D
E
T
U
G
<
/
I
K
O
C
C
E [
O I ./
´
@
D
M
O@C
@ I
A P
C
A
y(t) = AK 1 − p
T 3
U
K
[ GH
B@
J
O
H
K
8
D
F (p) =
1
G
B@
C P
L F
UZ K
H
34
F (p) =
K
D
L
H
B
D
K
H
H 2
D
H
O
B@
1−ξ
U
]
U B@
2
GH
p2
T
H
I
G
@ D
K C
G
U
9/
^
E
C
H
Q
B
E
H
^
H
BTGH
D
G 2
k
E
K
X
H 9
U H O@
H
G
U I O
!
O@C
U a
Kω02
G :
U
H H
O
E
/ ;
O
O
T
M K
U
P
V
R E #$ "
4 1
ξ
B@
C
H
H
U
B@C
H
+ 2ξω0 p + ω02
T S
ξ=0
.
0.25p2 + 0.25p + 1
E
A
%
G 1
.
³ p
C
l
H P
H
A O
0
T
!
G U
H
U
C K
.
O
u
C
y(t)
34
U
X
B@ I
P
&
H
A
B@
EZ
I
2
E
l
I@ EZ
() '
<=
H
K
B
O D
GH
e−ξω0 t sin ω0 1 − ξ 2 t + φ
B@
K
.
E
G
´
H
C
U
B
5
#
∗
E
4
C
V
CH
D
J G
L
3
L K F
0<ξ<1
O@C
P
5
[ K
<
I@
\
G
A
C /
C
B
H
M
I
H X
<
/
B@
B@
E >
¢
^ ^
[ m [ [ [ n [ m [ [
cd b cd b
A
X X X
G G G G G
£ w w H w w H
@ ª A A A A A
r
L -
B
1 TC
I I I
e Z Z Z e
H H H H H H
L K L K
f f
M E
O O O O O
3
CO
hi g hi g
_ _ _
G G G G G
T
H
U
D
5
U U U
C
¨
B 6 TC TC
C
D
C
I
O O O H
H O
I
@
X X X
U
7 UZ
E
ej ej
K H U H U H U
G
H H
ξ = 0.25
L
O
D
B
8 OG OG OG
O
9/
H
C C C T
B@
I I I
M
8
OT U E E
EZ
K
U
M M M GH G
T G M
./
T T T T
I A
O 9 @
K K K
M
U
D D D 4
1
B@
CH
D
I U T
D
GH
H
H
I L
B@C B@C B@C
B@ B@ B@
3 H H H K
- B
I
N E
O@
B@
CH
C
C C C
E
tp
< ^ ^ ^
N
@ D
C
5
K
U x
O
A
3 G
@ @ @ L
D D D
S O
[
H
K = 0.5
G
Tp
/ D
E
O A O A O A
H
T
D < M
x
I
_
E
I I I
E E E
H
H T
H H H
G
D D D
P P P 3
D
H
;
G
u
O
A B
U
U D D D H U
M
H H H D
U
B@
C CH
7- U
7
N
U U U
N
H
U U U
S
5% tr
B B@C EZ
G H I
GA
.4
54
C T x
I GH
l @
D
[
EZ EZ EZ
K=2
I@
[ A
ω0 = 4
[ C
ξ = 0.6
A
1
H
H
L
3
Y
I@ I@ I@
E
D
H H
3
H
X D
C
L K C C C H
k
H
k
J
:
I ./
D
I
H
H
T G
D
H
I
C C C G
I
3 4
L F
H T K
U
U
A
J U J J
H
H
U U
H
U H
U U M
H
GH
U
C
L F L F L F
K K K 4
C
6
B
C
M E
C
M
H H H
T
2
7- K
3
U
C
K
E
A
K E
U
K C
/
H
U CH
J O
@
-
U
I
H A
U
m P L
O
E E E
K
[ GH
O
3
O
H
8
E
D
E
H A A A
D
@
8
B@
GH M
<
o
GH M @
3 o
GH M
t
@ C
H
D
A E
GH
/
B@ B@ @ B@ G
A W
GH S
¤ @
H
K
H B
I
G
GH
- 8 E E E
© <
2 /
E
H
G G G
UZ
B B l B
l [
D%
l U
GH GH GH
34
CB
O A
E E E
U E
9
U
GH
B@
4 1
G
H H H K
AO
H
a O a O a O
A A A
C
\ A
:
H
H H H
G K M
1
-
U U
C C C
C
U U U
AO
L L L
5
S O
M
K
G
T X D
H
K
G K M
A
L L L
<
H H H C
¡ ¡ ¡
I
B@C
O
H
3
D
C C C
K
T
G A A A
G
H
X D
7
H
L
K
¥
C
I
T G G G
M
AO AO AO
C
H
H
D
3
T
O O O
L E
¦
4 E
GH
B@ @
H
B
G K M G K M G K M
§
H
C C C C
K
K
D D D
H
U
K
[ K K K
X D X D X D
H H H
C C C
I I I
U
ω0 = 2
O
. .
. -
9
9 9
9 I I I- 3
W
0
s l
.
P P
. /
-
. 2 .
D -
6 6 /
-
6
6 6
∗
5
5
z /
} | 9 9
I O /
9 0
G G G 6
4.
~ C C
1 1
g g
r r
F E 1 2 -. ,
v
q
/
4.
9 / -
9 9 89
1 1
/ /
4
1 1
3
9
Y X Y X Y X
I
7 /
5
/ 7 /
9. -
2
4 9. 4
8
. .
. 1 7
7 9
84
/
-
2 2
F 6 6
4 7
4
4 5
4.
4
4
Z Z Z
0
4 4
C
Q
/
. .
C 2
4 4. 4 . .
C C
F1 (p)
F1 (p)
F1 (p)
1 0
89
3
5
9 . 5 -5
1 2
D D J
6 5
7
F (p) =
F (p) =
4.
.
\ \ \
5
h
/ /
8
4
/
4
Y Y Y 2
/ /
6
/ .
4
- / - / 1
4. 3
6
.
F E F E
6
84 /
M
.
-
2
4 4
- C
9
4 .
I
^ ^ ^
/
- 0 0
1 1
9 7
5
/ 5 C
1
4
a
1
. . . 1 2 1 2
5 /
6 D
M 4
4
5 5
0
/
/ I
- 8 .
. / /
` _ n _m `a _
F2 (p)
F2 (p)
4 9.
2
G . 2
a b
C C
6 /
9 9 8
3 3
/
a
V
4
-
E
F
5 .
4 4 4 .
1
.
1 0 6
4
a
4.
F 4
1
/
x C C
b
4
. .
1 2
p−10
C
784
6 6
9
. D
D a /
-
e 9 G
/ /
c
/
-
.
I
t a
9
1 789
1 F
/
/
3
8
5
/
8
b
F E
F E u 4. 8
4 84 T
C 6
a
- / 1 -
6 C
1
.
6
O 4
84 4.
C 1
/
.
d ^ b 4
6
78
9 C
9
8
/
F (p)
F (p)
/
6
0
1
G /
5 9
Y G G
.
6
4
2 9
. 4
. 8 78
4
d c
4.
/
5
3
h a
0 .
C
.
◦
4
8
G
/
2 8
6
F
D G
C d c
C
4. 84
Z /
3
G
F2 (p)
{
`a
.
C
/ 84 4
-
9 1
/
9
4
2U
6
;< :
{ `
F E
h - /
9
4
8
]
4.
6
4
/
C
^ 4
4
1 6
4
- / 6
C
V 0
.
]
9
5 2
I
6
n
F
^ > =
d ^
6
9 /
. 9
h
/
1
6
5
_
4
0
1 G
4
.
h
.
/
- C
1
2
d ^
6 4
2
`
7
F (p) =
?@
6
6
I 2 P
/ 2 <
9
G I
4
9
0 /
n
1
I
6
3 h
2 G `
9
/
4 1
o
4
e Y
F (p) =
6
5
C
j I
C 4
CU
z 6 >
A
6
{ 9
C
j
Q
7
I
Y
6 e
4. 4
9
8
/
Y
1 7 @
f \ Z B
4 Z 7
8 / f
C 8
F1 (p) + F2 (p)
4
I
.
J 1 2
6
.
- 7
. C
5
I
k
k
4
J 1 2
-
`a
f \ O
4.
Z C !
\ 5 84
"
6 Q
4
/
8 /
.
C
6
4
- 9
8 /
ap + b
4 .
.
1
`a
-
4
8
d c -
.
.
D
h
8
L K
8
6 4 .
9
/ 4
$% #
ap2 + bp + c
4 6 -
6
` J
C
/
.
4
F E
2 h
/ 5 1
8 D
4.
8
h
-
8
6 6 /
.
I8
&
i /
1
f I8
1 I
F E
6
i /
p 1 9
4.
4
4
/
/ - /
!
5
"
/
h
84 E
- C
.
4
.
f
8 8
F1 (p)
84
E
/
J 1
F1 (p)F2 (p)
9
_ -
G
5
. G
5
/
6
1
`
'
4
M 4
4. /
.
9
4
4. 5
.
C /
q
C 8 H
.
2
2
)* (
O
1
4
4.
.
9
6
1
9
/
/
5
k
4 -
5 k
S
∗
89 / F
/
-
6
.
F2 (p)
8
G 4.
4.
6 +
1
/
-
4
/
/
0
4. 2
/
7
4 Q
2
4
Q
.
P 4 .
i
0
1
I-
6
1 2 C
/
9 I-
4 /
9 /
I
8
5 C . 8
1
8
-
C
C F 5
3
4
6 4 3
3
. 9
i y
/
/ -
7 I
.
.
V
5
/
8 / 8
8 L
F 8
3
F
1
. 4
6
1
9
I
1
6
.
Q
/ 4 D
4
G
/
J 1
F 3 6
G
6
4
/ G
.
6
8
5 C 8
N
4.
- F E
4 4
- - -
6
2 2 2
w V
4.
ª
q q
4 4
9
3 / 3
/ /
9 9
-
2
8 8
g g
P P
F F
1 1
6 6
Y X
I
C = 10
/ /
1 a « 9 1
9 4 4
4 4. 4
I
C
. .
.
I
.
2 2
^
6 6
C
. .
Z
Y
8 8
. .
4 C 6
3 1
/ [\]
5 1
1 9 /
/
1
8
D D a
1 1
Z \
C
4 4
/ 84 /
E
84 . .
E
[\
]
6
I 6 I
1
F E F E
a
4 4
4 4
\ .
. .
6 - Y
C F
Z
.
F
8
1
-
D D
4 4
6
1
C C
- /
1 1 / Z /
p
9
4
F E F E
9 \ I 6 9
- 4 4 e
2 /
h h
0
1
4. p 1 2 4.
6-
9
6 6
.
6-
n _m
0 /
5
- /
e
I I
. 9 . 9
3 ¨ ¦§ ¥
/ /
1 1
/ /
/
6
/
4
n _m /
I
¤
2
8 G
a
2
6 6
e
X(p) -
2
/ w /
- -
4 4
9
4
9
C C
T
5
0
4 4
5
a
9 0 0
C
8
a
e
6
2 r 2
0
5
b
8
8 / { 8 /
- -
Y
. .
©
4 5
I
a
F (p) =
-
7- 2 7-
F (p) =
4. 8 b L4 8
6 6
-
2
c
4 4 d
L4
Y
I
- -
h h
-
-
M 2
.
8 8
4
L4 `
7 C 7
d c
C=2
/ /
i i
7 7
1 2 1 2 \
- / -
C
O -
2
C
`
I- I-
4
a 3
|
4
-
2
0 e
/
2 9 9
/ /
1 £ 1
8 L5
Y
5
F
Y
G
£
1
L5
84
6 1 6
.
x 4
5 5
4 6
Z
]
84
6
L
C 5 C
]
[
.
[ -
`
2 `
1 O
2 1
p2 + 2p + 2
4.
I8 I8
9 1 6 1
-
/ - /
2
4
¡ 6 f
- F (p)
84 84
E E
y /
Z 6
-
2 i
p4 + 4p3 + 8p2 + 8s + 3
0
p4 + 2p3 + 2p2 − 2p − 3
5
.
4 4
- /
C
4.
9
C C
[
4
[ n
6
9 9
n
/
p
5 5
5 / k k
p
4
9
8
/
1
]
t
/ F /
Y (p)-
Q Q
4
G G
6
6 Z¢
3
Z¢
4.
`
/
-
I
- /
0 6
1
1 2
- / 4
6
5
G
T
5
2
C
2
8
{
Le but de la régulation est de « rajouter » un autre système (le correcteur) de fonction de
transfert C(p) de façon à ce que le système global (réfrigérateur + correcteur) se comporte
« bien ». Il y a deux façon de rajouter le correcteur :
TD n◦6 : Erreur et précision en régime permanent∗ 1. en boucle ouverte (le correcteur et le système sont en série) ;
2. en boucle fermée.
I.U.T. Génie Chimique & Génie des Procédés - 2ème année d
d Boucle fermée ?
H(p)
Boucle ouverte ?
H(p) +
1 Intérêt de la boucle fermée r -
+ ² - C(p) u - G(p) -?
+ -y
+ -
6
r - C(p) u - G(p) -?
+ -y
Dans le TD précédent, nous avons vu que la boucle fermée est intéressante car elle permet
entre autre de stabiliser un système instable. En même temps, il est important de prendre Dans les deux schémas, r est la consigne. Le correcteur est donc chargé de trouver l’entrée
des précautions car la boucle fermée peut aussi déstabiliser un système stable. Rappelons u du réfrigérateur pour que ce dernier maintienne une température de 3 ◦ C (la consigne) à
que la stabilité est une propriété essentielle qu’il tient de toujours vérifier avant toute autre l’intérieur sans l’intervention d’un être humain. De la même manière, il est chargé de trouver
chose. Nous allons voir dans cet exercice un autre intérêt de la boucle fermée. l’entrée u de façon à ce que la température à l’intérieur du réfrigérateur reste constante même
Considérons un réfrigérateur que l’on modélise par la fonction de transfert suivante : si la porte reste ouverte par mégarde. On veut donc que y vale 3 lorsque r vaut 3 même
2 lorsque d n’est pas nulle (en régime permanent) : le système global se comporte bien.
G(p) = .
1 + 1000p
1.1 Boucle ouverte
Cette modélisation est approximative, mais c’est pour donner une idée et pouvoir faire un
exemple simple. Examinons d’abord le système : 1. Cas de la poursuite : déterminer la fonction de transfert entre y et r (d = 0). Pour
simplifier, nous prendrons C(p) une constante. Quelle valeur doit-on choisir pour que
u - G(p) -y la sortie tende vers la consigne lorsque celle-ci est constante ?
2. Cas de la régulation : déterminer la fonction de transfert entre y et d (r = 0). Est-il
possible de changer cette fonction de transfert par le choix de C(p) ? Qu’est ce que cela
L’entrée u est une tension qui fait tourner un moteur (ce qui refroidit l’air à l’intérieur du veut-il dire ? Que se passe-t-il si la porte reste ouverte ?
réfrigérateur), alors que la sortie y représente la température intérieure du réfrigérateur.
N’oublions pas que lorsque l’on ouvre le réfrigérateur pour prendre quelque chose, de l’air
1.2 Boucle fermée
rentre à l’intérieur, modifiant la température intérieure du réfrigérateur. C’est ce que l’on
appelle une perturbation. Ceci peut se présenter comme dans le schéma suivant : 1. Cas de la poursuite : déterminer la fonction de transfert entre y et r (d = 0). Pour
simplifier, nous prendrons C(p) une constante. Quelle valeur doit-on choisir pour que
d la sortie tende vers la consigne lorsque celle-ci est constante ? Pouvait-on s’en douter ?
? Est-ce possible ?
H(p) 2. Cas de la régulation : déterminer la fonction de transfert entre y et d (r = 0). Est-il
possible de changer cette fonction de transfert par le choix de C(p) ? Qu’est ce que cela
+
u - G(p) -?
+ -y veut-il dire ? Que se passe-t-il si la porte reste ouverte ?
Prenons maintenant pour C(p) la fonction de transfert :
où d représente la température dans la pièce où se trouve le réfrigérateur. La fonction de 1
transfert H(p) représente la dynamique d’échange de chaleur entre l’intérieur et l’extérieur C(p) =
Ti p
du réfrigérateur : l’échange n’est pas instantané. Nous prendrons un exemple simple pour
H(p) : avec Ti strictement positif. La sortie tend-t-elle vers la consigne (en régime permanent) ?
1 Pouvait-on s’en douter ? L’effet de la perturbation est-elle rejeté (en régime permanent) ?
H(p) = .
1 + 1000p Pouvait-on s’en douter ? Que se passe-t-il si la porte reste ouverte ?
∗
pour toutes questions, e-mails : {msekher, eric.magarotto}@greyc.ensicaen.fr
1 2
2 Choix du correcteur pour la garantir précision et rapidité
+
X(p) - ²(p)- U (p)
- Y (p)-
C(p) G(p)
Soit le système : -6
+
X(p) - ²(p)- U (p)
- Y (p)-
C G(p)
-6
Y (p) 1
avec G(p) = U (p) = p2 +9p−10
.
Application numérique : A = 2, τ = 4. 2. Déterminez la condition de stabilité sur P pour obtenir un système stable en boucle
fermée.
3. Mettre F2 (p) sous forme standard et exprimez ses paramètres caractéristiques en fonc-
3 Choix du correcteur pour la précision
tion de P .
Soit le système : 4. Calculez la valeur limite de P pour que la réponse indicielle reste en régime apériodique.
On fixera cette valeur obtenue à la valeur entière la plus proche.
+
X(p) - ²(p)- U (p)
- Y (p)- 5. Déterminez la meilleure erreur de position possible en réponse à un échelon unitaire et
C(p) G(p)
-6 garantissant un régime apériodique. Exprimez cette valeur en pourcentage.
10
avec C(p) qui possède la fonction de transfert C(p) = Kp et G(p) = (p+1)(p+10) . Déterminer
K tel que le système en boucle fermée soit stable et l’erreur de vitesse inférieure à 10%.
3 4
> >
LS <
∗
<B I @P < r IUDR 453
S
DR> hi D=B @P
tuv s DL F HD= ?C >D D> kjl D FA
a F?dd LBR ` L= 67
tu w SD F M
m
◦
E F mn 89 A@Q 89:
Qe (t)
g I<= FF? >R CI N D>>E
xyzu w D HD G oi 9 F? => ;<
Qe
LB e< H FG
xy w = VD <= F? I @H gR> PD RD Da !
t|y} { "#
3
L@ DR F L> `< A@QV B<=> A?@
Qe
m /s
IB ?X <=U
~ a W R QI <H T CPb R< F CI Da > DL F?
Qp (t)
H< >D FGb \ Z[ Y #
m2
CD N @D T> CD B CD D FA W EHD FG
G(p) =
M] DF
h
~
> A?@C LSD < CI LS> D @ )
A@Q a W
h(t)
!
G(p) =
~ F? IF ??Q L bI R ^I
H(p)
C
Qe (p)
M <=T ?
> @D DU ? LB CI JK "#
=
~ IB
0, 5
U @D L CI bI LS FA? DLB F + *#
1 + 4p
.
<R F T ?@D N D M
1
α
=K ?@I LD F =R X< D= D @P
<L> LI
Qp (t) = αh(t)
IBL > OC . ,-
1 + Sα p
D BL a= KD _I
∗
= VD @ =< F? <L>
= K PI
/
h(t)
IE<=> <B F?
R QI +*
Qp
>D L= c R>D ?J L= @P
m3 /s
~ BI=> F @D IU<=>
0
2
N B W = VD F @D D
m
~ B XD UL L FA? R QI
D=<=> e a W
g F => @I <LD>B f LLD F bI => @I LSD #! 1&
~ !2
Qe (t)
LB B
α
CD CD= <=B CP #
D F CD CD Md@D X< D
S
A@ NN p A@
cm
N @ HT C
D @ F ?Q Mq F ?Q a F?
LR Q<
>L L= T> R< RD U @D >L > @D U HI= TL C > qD
D> F? o i ?
C
BDH CD X DR IB A??C 76 8
B CD LDB F N D> 6 BH ?X< =R> @< BI F @D <BL ? o 9i
BIT D ?QC DA CI RD> L<R @D D=B I F IU DR f SI= R< CI N L 6 BT ACI 6i
>D F =R X< >E = A?C R F ¥ on kl9 B< D F <=T ? <B B< F? @ HH >D OC > ?e imm <=> mm 94
T ?@Da H ?@D @ D>D C IB ? D 9 DR BD CP B> C< 949l D 9l
IB =L ? =< F? HD FG D ¦ <L F 6 DL F DL <LB LS ? D> CD L CI bI =I = F @< LS 8k9 DHT F 8k9
Qe (t) = 0, 02 m3 /s
R I= K= ? LD F N @D RT>
X(p) -
= P <=> DD C KI= R< CI B> @D LDL F bI I > F??C @P >L B < <
h(t)
UI @D CD
BDT @ <BDR <B F? LB JK <HH HH
U (p)
C
@D a R> @D DA
BTI @D R< F CI F R fD D A > DT> DL F
< B
BT A D= ? D @ L= @D
C
¦ F PD DHT> F R< F??X ?@D L CI @D HH =< CI <B>D F? L bI B AA?@ <=> FMd <L FF RA
F (p)
??Q DIL @P
cm
N HH I= D a I @P
= bD C<
Y (p)-
== q< B c@ BT ACI ? <= F? D=L F DLB F N D= a = FA?C > A @ DHT F <= B Q<
DBL DB ? @D BD= e D CD >D @D
CI H<= <=> HDT> G @D N D C D @P TTI <L> BDT A? <BD TB L= @ IHT @P
=D =DB L Q< <=>
R BD F D < TIT F LR ¢Q< L> ¡< ILB DD C IU<=> L> C D C D=T @ R A
DD F C D @CI
CI
@D ??Q D AC CD D A L @D <L = bD C<
B
L F ?C
?C D @D
h(t)
F A
h(t)
C = 2
UI L Q< BH CD XD
LB JK = BH XD L> F? N D C DBDLB CP > QD <H N R<H IH @P
0, 02
Fig. 1 système
L'objectif de la régulation est d'imposer une trajectoire particulière au signal de sortie y(t).
Au lieu d'analyser le comportement de la sortie à partir des équations diérentielles il est plus
facile d'utiliser les transformée de Laplace (voir TD1, partie Transformée de Laplace).
L'idée ici pour maîtriser la sortie y(t) est d'utiliser une boucle de régulation (gure 2). Le bloc
C(p) est appelé le régulateur. Il génère le signal u(t) qui va permettre d'obtenir le signal y(t)
désiré.
D'après le TD5 on ne peut choisir C(p) n'importe comment, sous peine de rendre la boucle
instable. De plus si on veut que la sortie y(t) suive la signal de consigne x(t) il faut placer plus ou
moins d'intégrateurs dans la boucle de régulation (cf TD6). Quant à la rapidité, les trois actions
(proportionnel, intégrateur et dérivateur) peuvent être utiles.
L'étude de quelques procédés (cf TD1, partie Modélissation) a permis de montrer que le com-
portement de système parfois relativement complexe peut être analysé à partir de modèle du
premier ou second ordre.
L'objectif de ce TD est de présenter 2 régulateurs possibles : P et PI. Le comportement de ces
régulateurs est présenté sur des systèmes du premier et second ordre.
A
G(p) =
1 + Tp
∗
pour toutes questions, e-mails : {msekher,eric.magarotto}@greyc.ensicaen.fr
1
+ ε(t)- C(p) u(t)- G(p)
x(t) - - y(t)
-6
C(p) = K.
1. Calculer la fonction de transfert entre X(p) et ε(p). Montrer qu'il s'agit d'un premier ordre
et donner sa constante de temps T 0 .
2. Déterminer l'erreur statique de la boucle lorsque l'entrée x(t) est un échelon. Théoriquement
pour quelle valeur de K cette erreur est elle nulle ?
3. Pour la valeur de K précédente il peut apparaitre quelques problèmes. Pour le comprendre
calculer la fonction de transfert entre X(p) et U (p). Déterminer ensuite l'expression tempo-
relle de u(t) lorsque x(t) est un échelon et tracer l'allure de u(t). Conclusion ?
4. Application : A = 2, T = 3. On suppose que u(t) ne peut dépaser une valeur de 10. Lorsque
l'entrée x(t) est un échelon unitaire, en déduire la valeur maximum de K . Calculer alors les
valeurs de T 0 et de l'erreur statique.
1. Calculer la fonction de transfert entre X(p) et ε(p). Montrer qu'il s'agit d'un second ordre
et donner les expressions de ξ et ωo en fonction des diérents paramètres de la boucle : A,
T , K et Ti .
2. Discuter de la stabilité de la boucle si tous les coecients A, T , K et Ti sont positifs.
3. Déterminer l'erreur statique de la boucle lorsque l'entrée x(t) est un échelon.
4. D'après la question précédente, la présence d'un intégrateur sur C(p) assure une erreur
statique nulle pour une entrée en échelon. Maintenant il s'agit de voir l'inuence d'une
perturbation d(t) appliquée sur le signal de sortie y(t) (gure 3).
2
d(t)
Calculer la fonction de transfert entre D(p) et Y (p). Lorsque cette perturbation est un
échelon, déterminer l'inuence sur y(t) à l'inni, c'est à dire déterminer limt→∞ y(t).
5. Application : A = 2, T = 3 et K = 1. On veut ξ = 1. En déduire Ti . Déterminer l'expression
de y(t). Calculer le maximum de y(t) et l'instant de ce maximum.
C(p) = K
T (p)
2. Donner l'expression de la fonction de transfert en boucle fermée F (p) = Tref (p)
du système.
3. Calculer la plage de variation de K pour que le système soit stable.
4. On désire régler notre système avec un facteur d'amortissement de0, 5 . Calculer l'expression
de K en fonction de a.
5. On suppose maintenant que a = 20, calculer l'erreur statique en % si on envoie un échelon
de température de 10◦ (application numérique).
3
6. On souhaite améliorer cet asservissement en supprimant l'erreur de position. Dans ce but,
on remplace le correcteur proportionnel par un correcteur PI de fonction de transfertC(p)
à !
1
C(p) = K 1 +
Ti p
T (p)
Donner l'expression de la nouvelle fonction de transfert en boucle fermée F (p) = Tref (p)
du
système.
7. On suppose K > 0. Calculer la plage de variation de Ti pour que le système soit stable.
8. calculer l'erreur statique pour ce type de régulateur si on envoie un échelon sur Tref .
4
T.D. n◦9 : Exemples de régulation∗
I.U.T. Génie Chimique 2ème année
1 Amortisseur
On considère un amortisseur. D'après le T.D. sur la modélisation (T.D. 1), on sait qu'il s'agit
d'un système du second ordre. Sa fonction de transfert est donc du type :
Aω02
G(p) =
p2 + 2ξω0 p + ω02
où A est le gain statique, ξ le facteur d'amortissement et ω0 la pulsation propre. Mais on ne
connaît pas les valeurs de A, ξ et ω0 . La seule information disponible est la réponse indicielle de
ce système donnée par la gure 1.
Step Response
5
4.5
3.5
3
Amplitude
2.5
1.5
0.5
0
0 2 4 6 8 10 12 14 16 18 20
Time (sec)
1
X(p) +
- - - G(p) - Y (p)
P
-6
3. A partir de maintenant, on choisira la valeur de P pour que le système bouclé soit stable.
Déterminer alors les valeurs du gain statique A0 , du facteur d'amortissement ξ 0 et de la
pulsation propre ω00 du système bouclé F (p).
4. Suivant les valeurs du gain P , les pôles de F (p) sont soit réelles, soit complexes. Calculer les
pôles de F (p). Tracer l'évolution dans le plan complexe de la position des pôles en fonction
du gain de correction P .
5. On veut ξ 0 < 1, déterminer alors la contrainte sur P .
6. On veut que le premier dépassement de y(t) pour une entrée en échelon de x(t) soit inférieur
à 20%, déterminer alors la contrainte sur P .
+? +?
ϕ(p) - G0 (p) - - - Te (p)
+ +
Montrer que G0 (p), G00 (p) et G000 (p) sont tous les trois stables : on notera K1 , K2 et K3 les
gains statiques et τ la constante de temps.
2
2. Pour l'instant, on considère la boucle de régulation de la gure 4 : les perturbations liées à
T1 et T2 sont négligées (T1 (t) = T2 (t) = 0). Pour concevoir un correcteur C(p) et ajuster ses
+ ε(p) ϕ(p)
Tc (p) - - C(p) - G0 (p) - Te (p)
-6
paramètres, il faut d'une part connaître le modèle du procédé et d'autre part spécier les
performances : la boucle devra se comporter comme un système du premier ordre, l'entrée
et la sortie de ce système bouclé sont représentées gure 5. En observant la gure 5, en
Tc (t) Te (t)
6 6
50% 50% .........................
....±
...
...
...
- t .. - t
τo
Fig. 5 tracé de Tc et Te
déduire l'expression de la fonction de transfert F (p) entre Tc (p) et Te (p). Déterminer l'erreur
statique du système en boucle fermée.
3. Calculer l'expression de F (p) en fonction de C(p) et de G0 (p). En déduire l'expression de
C(p) en fonction de G0 (p) et F (p).
4. Montrer que C(p) est un régulateur PI de fonction de transfert :
µ ¶
1
C(s) = K 1 + .
Ti s
Calculer alors les valeurs de K et de Ti en fonction de τ , τo et K1 .
5. On désire que le système réponde à l'échelon aussi vite en boucle fermée qu'en boucle ouverte
c'est-à-dire τ = τo . Dans ces conditions calculer K et Ti .
6. Calculer la fonction de transfert H(p) entre Tc (p) et ϕ(p). Donner l'expression de ϕ(t) en
fonction de Tc (t).
7. Les valeurs maximales de la consigne Tc (t), de la mesure Te (t) et de la puissance ϕ(t) sont
de 100%. Quelle est alors la valeur maximale de consigne Vc admissible pour démarrer
l'installation sans saturer ni la commande, c'est à dire la puissance ϕ(t), ni la sortie, c'est-
à-dire Te (t).
3
T1 (p)
- G00 (p)
+ ε(p) ϕ(p) +?
Tc (p) - - C(p) - G0 (p) - - Te (p)
-6 +