Cours Syst Asservis
Cours Syst Asservis
1.1 AUTOMATISATION :
Ensemble des procédés visant à réduire ou à supprimer l’intervention humaine dans les processus de
production
La régulation automatique, actuellement rebaptisée «automatique» est noyée dans les techniques
modernes de commande- robotique, productique etc., en raison surtout de l’apparition de l’électronique, puis
vers les années 60 du microprocesseurs et donc de l’informatique.
Mais il est utile de souligner que les vieilles techniques de régulation classiques restent encore très utilisées
dans l'industrie et elles ont encore de beaux jours devant elles car, la théorie en automatique avance bien
plus vite que l'application et cela, parce que les moyens informatiques sont plus «performants» que la
connaissance du système à traiter c’est à dire le modèle.
Il est aussi intéressant de noter qu’aujourd’hui, les mécaniciens souhaitent parrainer l’automatique car, « la
robotique c’est l’automatique » ,disent-ils, et les informaticiens ont les mêmes ambitions car l’informatique
industrielle est leur apanage. Et l’automatique dans tout ça ? Mais cette question, d’actualité d’ailleurs, est
sans doute la conséquence des transformations des sciences de l’ingénieur subies grâce (ou à cause) de
l’informatique.
1.3 HISTORIQUE :
2 Définition
2.1 Systèmes
2.2 Signal
Grandeur physique générée par un appareil ou traduite par un capteur (température débit etc.)
Quantité susceptibles de changer de valeur, elle est associée à la grandeur physique qu’elle représente à un
instant donné dans un système SI exemple : T1(11) = 100 °C (à t=11 sec., T= 100 °C).
On distingue :
• Signal d’entrée : indépendant du système, il se décompose en commandable et non commandable
(perturbations)
• Signal de sortie : dépendant du système et de signal d’entrée. On distingue sortie observable et non
observable
4.1 Définition:
Système à boucle ouverte : système dont on a aucune information sur la grandeur à commander.
Exemple: Le réglage de la température du four est assuré par une personne extérieure (de la salle de contrôle), il n'a
donc aucune information sur la grandeur à régler.
Régulateur
Comparateur
erreur
Actionneur +
+ correcteur
Consigne -
Mesure
ou erreur
+ O
- u
e
r
r
e
u
r
CNA : convertisseur Numérique Analogique
CAN : convertisseur Analogique Numérique
On peut trouver aussi :
I/P : Convertisseur electro pneumatique
P/I : Convertisseur pneumo électrique
8.1 Asservissement
Un système asservi est un SYSTÈME DIT SUIVEUR , c’est la consigne qui varie .
EXEMPLES : une machine outil qui doit usiner une pièce selon un profil donné, un missile qui poursuit une cible .
● Déterniner pour ces deux exemples, la consigne Xc , l'écart ou erreur Xc, la Mesure de la sortie (observation
retournée) Xm , la grandeur de sortie Y .
Xc Xc Xm Y
Machine outil
Missile en
poursuite
8.2 Régulation
Dans ce cas, la consigne est fixée et le système doit compenser l’effet des perturbations.
À TITRE D’EXEMPLE: le réglage de la température dans un four, de la pression dans un réacteur, le niveau
d’eau dans un réservoir.
Xc Xc Xm Y
Régulation de T°
Régulation du
niveau d'eau
Échelon de position
SORTIE Stable Précision
Entrée (Consigne)
S t
t(s)
S t Erreur statique
Stabilité :
On dit qu'un système est stable si pour une entrée constante, la sortie reste constante
quelles que soient les perturbations.
● Les systèmes, dont les courbes de réponse à une entrée de consigne échelon de position (réponse
indicielle) ci-dessus sont numérotées sont :
a) instable pour les N° :
b) Stable pour les N°:
c) Rapidité
● Quel est le système le plus rapide ? Justfifier . De même pour le plus lent ?
● Tr est le retard apporté à la température T(t) compte tenu de la vitesse de l'eau et de la longueur du
tuyau .
● Tc Température de consigne souhaité par la personne , Tm la température mesurée par les capteurs de
la peau du corps .
Schéma Fonctionnel:
● Cependant, l'opérateur agit dans le bon sens et en modérant son action sur la vanne (oscillation de
diminuant), dès qu'il prend en compte correctement le retard Tr apporté aux conséquences de son réglage par
la vitesse limitée de l'eau et la longueur du tuyau .
● Un autre opérateur plus nerveux pourrait provoquer un comportement instable à ce système manuel dont il fait
partie, en tant que régulateur, capteur et actionneur .
Le régulateur proportionnel fournit une sortie U proportionnelle à l'écart ∑ , entre la consigne C,
et la mesure M .
BP est l'inverse Dans cette
de la pente de U(M) bande:
Dans la zone où U = Gr . ∑
U(M) = Gr . ∑.
BP = 1 / Gr
L'écart ∑ dépasse la
bande proportionnelle.
U est maximal.
Saturation Aop....
RADAR
● Définition :
Un système physique est dit linéaire si les différentes grandeurs de sortie sont liées aux entrées par des
équations différentielles linéaires à coefficients constants.
Propriété :
- si on change e(t) en k e(t), s(t) devient k s(t)
- si à e1(t) correspond s1(t) et à e2(t) correspond s2(t),
alors à k1e1(t)+k2e2(t) correspond k1s1(t) + k2s2(t).
Ces remarques relativement triviales jouent cependant un rôle important lorsqu'on veut connaître la réponse
d'un système à une entrée non simple (créneaux, dents de scie etc...)
● un système réel pourra être modélisé par un système linéaire s'il a approximativement les propriétés
ci-dessus .
e(t) = ; s(t ) =
e(t) = ; s(t ) =
Avec E(p) = L [e(t)] et S(p) = L [s(t)] ; les transformations de Laplace de e(t) et s(t)
● La transformée de Laplace permet de faire intervenir les conditions initiales e(t=0) et s(t=0) par
l'intermédiaire du théorème de la dérivation
Si les conditions initiales sont nulles (conditions de Heaviside) on peut exprimer la Fonction de transfert: T(p)
● Exemple pour un système linéaire du premier ordre régi par l'équation : . ds(t) / dt + s(t) = a. e(t)
S p Hp
Montrer que T p= =
Ep 1K p. Hp
S p Hp
T BF p= =
E p 1Hp K(p) = 1 .
Rp
T BO p= rapport des transformées de
Err p
Laplace du signal retourné et de la transformée du
signal d'erreur .
Rp
Montrer que T BO p= =K p. Hp
Err p
● Rq:
L' Échelon de position donne la RÉPONSE INDICIELLE L' Échelon de VITESSE OU RAMPE
Unitaire si A = 1
Échelon de position
SORTIE Stable Précision
Entrée (Consigne)
S t
t(s)
S t Erreur statique
Les modèles les plus utilisés sont les modèles des premiers et second ordre .
L' erreur statique n'a d'intérêt, pour le premier ordre, que pour un système bouclé .
En effet, s t=0 si K0=1 et s t =−9 pour K 0=10 alors que le système
fonctionne parfaitement dans les 2 cas et est même plus rapide dans le second .
H0 p
● Montrer que HF p= ;
1H0 p
KF K0 0
● En déduire que : HF p= avec K F = et F=
1 F . p 1K0 1 K0
● EVALUATION DE LA RAPIDITÉ DU SYSTÈME:
0
a) Montrer que le temps de montée à 95% de la la valeur Finale VF est t 5%= 3. F=F=3.
1 K 0
−t
1 1
On rappelle que : e(t) = 1. u(t) ---> E p= et que --> 1−e
p p.1. p
et que pour ce type de courbes exponentielles, les intervalles de temps se calcule comme suit :
s ∞ −s t1
t2 – t1=. ln
s ∞−s t2
Evaluation de l’erreur statique du Système en BF(retour unitaire)pour une entrée en Échelon unitaire:
soit, S t =
lim
p −>0
A.1−
KF
1 F .p
=A 1−K F = A
1
1K o
Système du 1er Ordre Représentation graphique de s(t) dans les cas B.O. et B.F. :
Ces tracés sont obtenus pour un échelon unitaire avec K0 = 2 et 0 =1s.
Dans le cas de la boucle fermée le retour est considérée comme unitaire.
● Conclusion:
1. L'ordre du système n'est pas modifié par le bouclage. Le système bouclé à retour unitaire reste du premier
ordre.
2. Le bouclage a diminué le gain statique du système : KF K0
3. Le sustème bouclé est plus rapide.On améliore sa rapidité en augmentant le gain Ko .
4. Le système en BF est plus précis , l'erreur statique est plus faible . Cependant, l'erreur statique du système du
1er Ordre en BO, n'a pas vraiment de sens, car si Ko = 1, l'erreur devient nulle .
Pour le système du 1er Ordre en BF(retour unitaire), plus Ko est fort plus l'erreur S est faible .
Exemple : Circuit RC .
Montrer que :
vs t =E 1−e
−
t
0
1. Fonction de transfert isochrone : On remplace p par j. dans la fonction de transfert Ho(p):
K0
H0 j. =
2.a0 j. ²
1 j.
n0 n0 ²
p0 2.
● La pulsation de coupure n0= peut être obtenue par mesure de p0=
1−a 0 ² TP
Tp, pseudo période d'oscillation lors d'un essai de réponse à un échelon unitaire .
Les coefficients d'amortissement vont de a0=0,1 , le plus oscillant, à a0=1 le plus amorti, non oscillant.
Diagrammes de Gain(dB) et de phase (°) , en bas , pour a0=0,1 (Bleu) et a0=0,4 (rouge)
● Le gain se lit sur le diagramme du haut pour tendant vers 0
● La pulsation de cassure n est à l'intersection des deux asymptotes, et aussi pour la phase de –90°
● La valeur du coefficient d'amortissement a0 se trouve par des abaques.
1−a 0 ²
2. Réponse indicielle(échelon unité) : e(t) = u(t) => E(p) = 1/p ; avec =arccos a0
T p 1
.sinn0 . 1−a 0 ².t]
−a . .t
Donc, Sp= et st =K 0. [1− .e 0 n0
p 1−a 0 ²
du C d²u C
u CR.C. L.C . = ue
dt dt Application numérique : R = 1 k ; L = 10 mH ; C = 10nF ,
2.a0 dst 1 d²st
st . . =K 0 .et , et montrer que : n0= 1 R
et a0= .
C
; Ko = 1.
n0 dt n0 dt² L .C 2 L
S p
Calculer la fonction de transfert : T p=
Ep
Montrer que S(p) pour une excitation en échelon de 5V , peut de mettre sous la forme:
5
Sp=
p 1
2.a0
n0
p
p²
n0 ²
● Calculer numériquement : an0 , K0 et n0 .
● Montrer, à l'aide des transformées inverses de Laplace en Annexe 2, que s(t) , réponse à un échelon e(t) de 5V
peut se mettre sous la forme :
1
.sin n0 . 1−a0 ². t ]
−a . .t
s t =5 .[1− .e 0 n0
=arccos a0
1−a 0 ²
● Faire tracer à la calculatrice, au tableur ou avec un logiciel de simulation (scilab...) , s(t) , la réponse à un
échelon e(t) de 5V .
+ ACTIONNEUR PROCESSUS
-
CAN TRANSMETTEUR
1. Action qualitative
2. Équation de récurrence :
Pour réaliser un filtre passe-bas FIR( RIF en Anglais) simple (qui atténue les fréquences élevées du signal), il suffit
d'effectuer la moyenne des valeurs entre l'échantillon présent et l'échantillon précédent, tel qu'illustré sur le schéma-bloc
ci-dessous et tel qu'exprimé par l'équation :
Yn
Xn
Equation aux différences (pas de récurrence ici) : y[n] = 0.5 x[n] + 0.5 x[n-1]
Intuitivement, on peut comprendre l'effet passe-bas de cette opération puisqu'en effectuant la moyenne, on atténue des
variations brusques du signal, ce qui conduit à un lissage du signal d'entrée.
Simulation du filtre par matlab(c) ou scilab(logiciel libre de l'INRIA en GPL) :
script matlab :
a=[0.5 0.5]
[h,w]=freqz(a,1,512,8000);
plot(w,abs(h));
Ce filtre à une fréquence de coupure de 2KHz pour une fréquence d’échantillonnage de 8kHz(la fréquence utilisé pour
les lignes téléphoniques).
1. La transformée de Laplace d'une fonction f(t) est unique . f(t) <--> F(p) unique
2. Linéarité: f(t) = a f1(t) + b. f2(t) = a F1(p) + b F2 (p) .
3. Produit: ATTENTION ! f1(t) . f2 (t) ≠ F1(p) . F2(p)
4. Dérivation :
df t +
f(t) --> F(p) Alors f ' t = --> p.Fp−f 0
dt
d²f t + +
de même: f ' 't = --> p [ p.Fp−f 0 ]−f ' 0
dt²
+
Fp g 0
5. Intégration: L[∫ f t dt ]= en effet,
p p
+
dg t Fp g 0
f(t) --> F(p) Alors si f t = −−−> Fp=p.G p−g 0+ d ' où G p=
dt p p
Remarques:
+ +
● si les conditions initiales sont nulles, (conditions d'Heaviside), f(0 ) = 0 et f'(0 ) = 0 , alors :
a) dériver dans le domaine temporel (t) revient à multiplier par p dans le domaine symbolique (p)
b) intégrer dans le domaine temporel (t) revient à diviser par p dans le domaine symbolique (p)
● Ainsi, les équations différentielles, se transforme en un systèmes d'équations algébriques, bien plus
facile à résoudre .
EXEMPLE: . ds(t) / dt + s(t) = a. e(t) donne . p.S(p) + S(p) = a. E(p) ; avec e(0) = 0 et s(0) = 0
− To . p
8. Théorème du retard: g t=f t−T 0−−>G p=e . Fp
Décomposition en éléments de première espèce (dans ℂ , ce sont toujours des éléments de première espèce):
La fraction est dite irréductible si l'on a simplifié , haut et bas les termes identiques .
P x R x
● =E x Si D° P < D° Q => le quotient E(x) = 0, si D° P = D° Q => E(x) = Cte = C ,
Q x Q x
si D° P > D° Q => E(x) = polynôme de D° = D° P - D° Q
P p A1 A2 A3 B1
● =Ep par exemple pour Q(p) dont les racines ou
Q p p−p1 p−p 1 p−p 1 p−p2
3 2
A3 A4
T p=0 si = 0 , p0 pôle double .
p−p0 ² p−p0
● Détermination des Ai :
A 1=p−p 1.T p pour p=p1 ; A 2=p−p 2.T p pour p=p2
A 3=p−p0². T p pour p=p 0 ; A 4 = lim p−p 0.T p
p ->∞
10. Décomposition en éléments simples de F(p) pour retrouver l'originale f(t), Transformée inverse.
D'après: http://www.cem2.univ-montp2.fr/~moreau/cours/laplace_c.pdf
du t
(t) est la distribution de Dirac ou fonction impulsion . (0) -> ∞ , (t) = 0 pour t ≠ 0 ; t=
dt
[1 −
−
t
e ]. u t
1
−
t
e . u t
1
p.1 . p 1. p
Annexe 3.
Table de Transformée en z et quelques méthodes ,