TpAsserv 2015
TpAsserv 2015
TpAsserv 2015
travaux pratique
ASSERVISSEMENT LINEAIRE
Ralise par:
Mohamedamine tahiri
Tarik mlilas
Ahmed mounsef
introduction
TP1 : Simulation Des Systmes Linaires
TP2 : Correction dun systme de troisime ordre laide
des correcteurs avance de phase, retard de phase et
combine
TP 3 :
Rgulation de niveau
I.
introduction :
MATLAB
Simulink : cest lextension graphique de MATLAB permettant de travailler avec des schma
en blocs, pour modliser et simuler des systmes ;
Blocksets : ce sont des collections de blocs SIMULINK dvelopps pour des domaines
dapplication spcifiques (DSP BLOCKSET, POWER SYSTEM BLOCKSET, etc.)
Toolboxes : ( botes outils) ce sont des collections de fichiers M dvelopps pour des
domaines dapplication spcifiques
I.
Y (s)
E (s ) =
0,5
1+0,1 s
un
Y (s)
E (s )
0,5
1+2 s+ s 2
b) - Pour = 0,15 :
La valeur finale stabilise de y(t) est : yp() = 0,5
D1 p =
D1p% =
y max y ( )
.100
y ( )
0,810,5
.100=62
0,5
On a
Y (s)
E (s )
0,5
2
1+0,3 s+ s
Y ( s) =
0,5
E (s )
2
1+0,3 s+s
0,5
1
2
1+0,3 s+ s S
0,5
Y ( s) =
S ( SS 1 ) ( SS 2)
S 1,2=0,15 j 0,99
avec
0,15 t
y (t ) volue entre
0,15 t
e
e
1+
2
10,15 ] et 10,152 ]
0,5
0,5
Donc : y th ( )=0,5
D1% = 100 e
1 2
0,15
On a =0,15
D1% = 100 e
10,152
On a = 0,15 et 0 = 1
Donc
III.
T P=
2
0 12
T P=
2
10,152
Tpth = 6,35 s
b) pour = 0,15 et K = 1 :
Le premier dpassement D1p% : D 1 p =
D1 p =
y max y ( )
.100
y ( )
0,560,33
.100=69,70
0,33
D 1 th =100 e
1 2
0,15
D1 th =100 e
10,152
=62,08
La pseudo-priode Tpth :
T Pth =
2
2
=
=5,2 s
2
0 1 1,22 10,152
pour = 0,15 et K = 2 :
D1 p =
D1 p =
y max y ( )
.100
y ( )
0,850,5
.100=70
0,5
D1 th =100 e
1 2
0,15
D 1 th =100 e
10,152
=62,08
La pseudo-priode Tpth :
T Pth =
2
2
=
=4,55 s
2
0 1 1,14 10,152
pour = 0,15 et K = 3:
D1 p =
D1 p =
y max y ( )
.100
y ( )
1,050,6
.100=75
0,6
D 1 th =100 e
1 2
0,15
D1 th =100 e
10,152
=62,08
La pseudo-priode Tpth :
T Pth =
2
2
=
=4,02 s
2
0 1 1,58 10,152
Thoriquement :
On a
H ( s )=
Y (s )
0,5
= 2
E (s ) s +0,3 s+1,5
=E ( s )Y ( s )=E ( s ) 1
s0
0,5
E(s)
s + 0,3 s+ 1,5
Y ( s)=
1
E( s)
s + 0,3 s+ 2
Y ( s)=
1,5
E(s)
s + 0,3 s+ 2,5
0,5
0,5
1
= 1 2
s +0,3 s +1,5
s + 0,3 s+1,5 s
2
( )=lim s (s )=lim 1
s 0
][
Y ( s)=
0,5
=0,67
s + 0,3 s+ 1,5
2
H ( s )=
Y (s )
1
= 2
E (s ) s +0,3 s+2
=E ( s )Y ( s )=E ( s ) 1
s0
1
1
1
= 1 2
s +0,3 s +2
s + 0,3 s+2 s
2
( )=lim s (s )=lim 1
s 0
][
1
=0,5
s + 0,3 s+ 2
2
H ( s )=
Y (s )
1,5
= 2
E (s ) s +0,3 s+2,5
=E ( s )Y ( s )=E ( s ) 1
s0
1,5
1,5
1
= 1 2
s +0,3 s +2,5
s + 0,3 s+ 2,5 s
( )=lim s (s )=lim 1
s 0
][
1,5
=0,4
s + 0,3 s+ 2,5
2
H ( s )=
K
s +3 s +2 S + K
3
s3
s2 3
s
K 0
6K
3
s0 K
Pour que le systme reste stable, il faut que les lments de la premire
colonne
Conditions de stabilit :
K >0
6K
>0
3
0< K <6
11
Pour K = 3 :
Pour K = 4 :
Pour K = 4 : () = 0
d) - Conclusion :
Lorsque K prend une valeur entre 0 et 6 ( 0< K <6 le systme reste stable et lerreur
statique est nulle.
=2,
=3 et
G(s)
9
( 1+ a s ) ( 1+ b s ) (1+ c s)
=1.5
A. Partie thorique :
1) - Calcule de la pulsation de rsonnance wr :
La fonction de transfert en BF de systme non corrig est la
suivante :
13
avec
G(s)
H(s) = 1+G( s)
Alors
H(s) =
H(s) =
H ( jw)=
( 1+ a s ) ( 1+ b s ) ( 1+ c s ) +9
( 1013.5 w ) +(6.5 w9 w )
2 2
3 2
9
( 1+ a s ) ( 1+ b s ) (1+ c s)
9
1+
( 1+ a s ) ( 1+ b s ) (1+ c s)
9
2
9 s +13.5 s +6.5 s+ 10
9
D(w)
14
dD ( w )
2
2
2
=0 27 ( 1013.5 w ) + ( 6.527 w )( 6.59 w )=0
dw
4
2
Donc 243 w +130 w 227.75=0
130+ 238273
=0.74
2243
Alors
Wr = 0.86 rad/s
Calcule du facteur de rsonnance M :
M=
H ( j wr )
9
= 7.386.
Donc :
Do:
Alors :
Si :
(BF) on a :
Et
Si :
(HF) on a :
Et
( G ( jw ) ) =0 :
donc :
R(G( jw))
( G ( jw ) ) =5,7
-1
-2
-3
-4
-5
-6
-7
-4
-2
Mg=20 log
avec
10
=0.85 rad /s
Mg=0.24 db
donc :
Pour w 0
G( jw)dB 20 log ( 9)
arg ( G ( jw )) 0
Pour w
G( jw)dB
arg ( G ( jw ) )
3
2
Cherchons la pulsation
w1
w 1=0.8606 rad /s
arctg(G(j
w1
)) = 180.91
Cherchons la pulsation w
arctg
C'est--dire
w9 w 3
6.5 w9 w 3
=
=0 6.5 w9 w 3=0
( 6.5113.5
w2 )
113.5 w 2
w =0.85 rad / s
1+s
1+s
Tracer le
courbe
secondes
Les marges de phase et de gain sont positives, donc le system est stable.
Correcteur retard de phase :
C ( s )=K
1+s
1+s
La marge de gain est positive et de phase ngative donc le systme pour ce cas
est instable.
Tracer le diagramme de Nichols
On obtient la courbe suivante :
Pour
Pour ce cas on voit que le systme devient plus stable que le cas prcdent donc le
correcteur retard de phase rend le systme plus performent que ce avance de phase.
2) - Correcteur avance et retard de phase :
C ( s )=
( 1+ 2 s ) ( 1+ 3 s )
( 1+ 1 s ) ( 1+ 4 s )
Les marges de phase et de gain sont positives pour les deux graphes, donc le system est
stable.
Pour
TP 3 :
Rgulation de niveau
1) - Etude du processus :
a) - On ouvre brusquement la vanne dentre que lon maintient ouverte .Si le dbit est
considr comme un chelon alors la variation du volume dans ce cas sera linaire
car :
Qe =
On a
Or
dV
V ( t )=Qe t +Cste
dt
b) - Lorsque la cuve est remplie la vanne de sortie est ouverte. Lorsque le niveau H
est stabilis au niveau h0 (niveau dquilibre). Le dbit Qe exactement le dbit de
sortie Qs (Qe0 =dbit dquilibre).
Etablir le modle linaire :
dh
On a: Qe =q e +Qe 0= A dt +Qs
Qe =q e +Qe 0= A
Q s=
h0 +h = h0
R
Q e 0=
dh
+Qs
dt
1+
1+
h
h0
h
h
1+
h0
2 h0
h 0 q =A dh + h0 (1+
R
dt
h
1)
2 h0
dh
1
+
h=qe
dt 2 R h 0
dh
1
+
h=qe
dt 2 R h0
1
1
H ( s )=Qe ( s ) H ( s ) As+
=Qe ( s )
2 R h0
2 R h0
H (s)
=
Qe ( s )
1
As+
1
2 R h0
2 R h0
H (s)
K
=
=
1+
s
Q e ( s ) 1+ 2 RA h0 s
K=2 R h0 et =2 RA h 0
h
dh
q e dt= Adh qe ( t0 )= A ( h c 0 ) q e =A c
dt
t
3) Partie simulation :
Schma fonctionnelle du systme:
1re cas :
3 1
Rsultat :
3 1
2eme cas : pour h c =0.5 mq e =0.667 m s
Rsultat :
Au niveau de la prcision on voit quelle est parfaite pour les deux figures, et
concernant la rapidit cest la mme.
4) - Commande en boucle ferme :
1) - Schma fonctionnel :
0.15 A 1
s ( 1+ 0.1 ) (1+2 S)
2.1 0.15 A 1
2.10.03 A 1
0
2.1
0.15 A 1
A 1 >0 et A 1<
2.1
=70
0.03
On a
arg ( F ( jw ) ) =
arctg
2.1 w
Alors
F ( s )=
3.32
3.32
=
3
s ( 1+ 0.1 ) (1+2 S) 0.2 s +2.1 s2 + s
M G=10 dB
M =15
La rponse indicielle
t r=20 s
t m=2.010.555=1.445 s
La sortie h(t) :
Alors
h ( t ) =0.0003 e10 t 2.7958 e0.5 t +2.79541.3944 t+
0.6640 2
t
2
Lerreur est :
( t )=hc h ( t )=0.0003 e10t +2.7958 e0.5 t2.7954 +1.5944 t
0.6640 2
t
2