Rapport Asservissement KOBEWO KCH
Rapport Asservissement KOBEWO KCH
Rapport Asservissement KOBEWO KCH
Paix-Travail-Patrie Peace-Work-Fatherland
********** **********
MINISTERE DE L’ENSEIGNEMENT SUPERIEUR MINISTRY OF HIGHER EDUCATION
********** **********
MINISTERE DES POSTES ET DES TELECOMMUNICATIONS MINISTRY OF POSTS AND TELECOMMUNICATIONS
********** **********
ECOLE NATIONALE SUPERIEURE DES NATIONAL ADVANCED SCHOOL OF
POSTES, TELECOMMUNICATIONS ET DES TECHNOLOGIES POSTS, TELECOMMUNICATIONS AND
DE L’INFORMATION ET DE LA COMMUNICATION INFORMATION AND COMMUNICATION TECHNOLOGIES
********** **********
GROUPE 3
ITT3 ALT RT
Rédigé par :
KOBEWO KOBEWO CHRISTIAN-HOPE
MANGA FERDINAND
ATOUBA TOICK
Sous la Supervision de :
Mme ABENA AUDE
2020-2021
2
I. OBJECTIFS
Modéliser et simuler les opérations liées aux systèmes asservis à l’aide du logiciel MATLAB.
Nous aurons ainsi l’occasion de mettre en pratique et d’expérimenter les notions des enseignements lors des séances de cours
magistraux.
Ces notions incluent les généralités sur les asservissements et régulations. Les réponses temporelles et fréquentielles à un système
d’ordre un et à un système d’ordre deux.
Un système régulé maintenir une ou plusieurs grandeurs physiques de sorties à une valeur constante ou de leur faire suivre un cycle
d’évolution et ceci indépendamment des entrées perturbatrices.
2
b) Système asservi
Un système est asservi si et seulement s’il comprend un dispositif qui va forcer les signaux de sortie à suivre au mieux les consignes.
Exemple : le système (bicyclette- cycliste)
Un système asservi doit permettre de faire suivre à une ou plusieurs grandeurs de sortie une loi évolutive et non fixée à l’avance.
II. Modélisation
En asservissement, l’ensemble des composants-système est considéré comme SISO ou MIMO, continus, causaux, linéaires et invariants.
Leurs description est faite à l’aide d’équations différentielles linéaires à coefficients constants (EDLCC).
Objectif principal de la modélisation = maîtriser la complexité
• Modéliser = abstraire la réalité pour mieux comprendre le système à réaliser
• Le modèle doit être relié au monde réel
La transformée en Z n'est rien d'autre qu'un cas particulier de la transformée de Laplace, bien connue de tous ingénieur. En réalité, il s'agit de
la version discrète de la transformée de Laplace. Elle est donc utilisée dans le monde du numérique, notamment en traitement du signal (pour
designer des filtres numériques par exemple) ou en automatique (pour asservir numériquement un système).
2
IV. Schémas Fonctionnels
Le schéma fonctionnel, appelé aussi schéma bloc, schéma de principe ou en anglais block
diagram, est la représentation graphique simplifiée d'un procédé relativement complexe
impliquant plusieurs unités ou étapes.
La définition des fonctions de transfert de schéma bloc et l’ensemble d’un procédé perme
construire un schéma fonctionnel qui forme l’étude d’analyse et de synthèse de la boucle.
2
ii. Système continu
Un système continu dans lequel entre deux valeurs de variables, quelles qu'elles soient, on peut toujours supposer une valeur intermédiaire.
Remarque : Dans de nombreuses techniques de transmission ou de traitement d’information, on utilise des signaux de type « discret » appelés
aussi quelquefois signaux « échantillonnés ».
On appelle fonction de transfert du système le rapport de la transformée de Laplace à la sortie sur la transformée de Laplace à l’entrée, ainsi :
2
Par identification :
Le gain statique K = 1 et la constante de temps
𝜏= 𝐿/𝑅.
d2 y (t ) dy ( t )
+2 m wo 2 + wo2 y ( t )=wo 2 x ( t )
dt ² dt
𝑤0 = Pulsation propre du système non amorti
m = Coefficient d’amortissement
En appliquant la transformée de Laplace, on obtient la fonction de transfert ci-dessous :
N (p) w20
G( p )= = 2 2
D( P) p + 2m w0 P+w 0
2
IX. Performance des Systèmes Asservis
Nous avons trois critères fondamentaux pour juger la performance d’un système asservi :
1. Stabilité
2. Rapidité
3. Précision
Etudier la stabilité d’un SA revient à vérifier que l’équation caractéristique D(p)=0 n’a pas de racine à partie réelle positive.
Il existe plusieurs méthodes. Nous pouvons citer entre autres :
2
Analyse des pôles par résolution de l’équation D(p)=0
Critère de Routh-Hurwitz
Critère de Nyquist
Méthode de Evans ou lieu des racines
Un système asservi est asymptotiquement stable tel que sa fonction de transfert H(p) n'ait que des pôles à partie réelle
négative.
Un système asservi est stable si chacune des racines du polynôme caractéristique possède une partie réelle strictement
négative.
CRITERE DE ROUTH
2
TP2 : Etude des Systèmes du Premier Ordre : Cas des fonctions continues
A. Objectif
étudier la réponse d’un système de premier ordre à un échelon unité, tracer les lieux de de Bode, Nyquist et Black-Nichols.
B. Méthode
Soit un système de premier ordre avec un gain statique k.
Identifier la forme générale des systèmes linéaires du premier ordre
Identifier la fonction de transfert liée à un système linéaire en fonction du gain statique et de la constante de
temps t
User des commandes Matlab afin de répondre aux questions posées.
1. Fonctions de transfert
Code h1 :
Code h2 :
2
Code h3 :
2
Réponse indicielle de h2 :
2
La courbe de la fonction h3 à un temps de réponse plus faible du système.
4. Temps de stabilisation
2
Montage pour une réponse à échelon unité
2
Code h1 :
Figure 1 :
2
Code h2 :
2
Figure 2 :
8. THÉORÈME
I. Bode
Code :
Figure 1 :
II. Nyquist
Code :
2
Figure :
2
III. Nichols
Code :
2
Figure :
2
9. INTERPRÉTATION DES COURBES
A l’aide des diagrammes obtenues nous pouvons donc dire que le système ayant pour fonction de transfert h1 dispose du
plus faible temps de réponse (donc est le système le plus rapide), que le système ayant pour fonction de transfert h2
quant à lui présente un temps de réponse moyen et enfin le système h3 présente un grand temps de réponse (système le
plus lent).
2
TP3 : Etude des systèmes du premier ordre : Cas des fonctions discrètes
I. Objectif
On considère un système échantillonné, tracer la réponse d’un système de premier ordre à un échelon unité, étudier la réponse temporelle à
un échelon d’une fonction de transfert dans le cas numérique.
II. Procédure
a) Identifier la forme générale des systèmes linéaires d’ordre 1
b) Utiliser les commandes Matlab pour une réponse aux différentes questions posées.
Travail à faire
III.
A. Réponses indicielles
Fonction de transfert échantillonnée h1
Code h1 :
Figure 1 :
2
Fonction de transfert échantillonnée h2
Code h2 :
Figure 2 :
Figure 3 :
2
Fonction de transfert échantillonnée h4
Code h4 :
Figure 4 :
2
Fonction de transfert échantillonnée h5
Code h5 :
Figure 5 :
2
Fonction de transfert échantillonnée h6
Code h6 :
Figure 6 :
2
B. Simulink
1. Montage sur Simulink
2. Réponses indicielles
Réponse indicielle de h1 :
2
Réponse indicielle de h2 :
Réponse indicielle de h3 :
2
Réponse indicielle de h6 :
2
−1 b 1 z−1
C. Etudes des réponses temporelles à un échelon unité de la fonction H ( z )=
1+a1 z−1
Figure 1 :
2
2. Fonction de transfert h2
Code h2 :
2
Figure 2 :
3. Fonction de transfert h3
Code h3 :
2
Figure 3 :
4. Fonction de transfert h4
Code h4 :
Figure 4 :
2
E. Représentation des réponses sur le même schéma
Code :
2
Figure :
2
La valeur de a1=-0 ,2 donne d’un système de premier ordre continu un échantillonnage rapide .
La valeur de a1=-0 ,9 donne un échantillonnage lent d’un système de premier ordre continu.
2
Figure obtenue :
I. Objectifs
Étudier la réponse système de second ordre à
un échelon unité, tracer les lieux de Bode,
Nyquist et Nichols.
2
II. Procédure
Soit un système linéaire de second ordre de gain k, de pulsation propre w0 et son coefficient d’amortissement m.
a. Identifier la forme générale d’un système de second ordre
b. Identifier la fonction de transfert d’un système linéaire de second ordre en fonction de son gain, de sa pulsation propre et de son
coefficient d’amortissement
c. Utiliser les commandes Matlab pour répondre aux questions posées.
III. Travail à faire
A. Création d’un fichier et traçage des réponses unitaires
a. Réponse de la fonction de transfert h1
IV. Code 1 :
V. Figure 1 :
2
b. Réponse de la fonction de transfert h2
VI. Code 2 :
2
c. Réponse de la fonction de transfert h3
VIII. Code 3 :
2
d. Réponse de la fonction de transfert h4
X. Code 4 :
2
XIII. figure 5 :
2
B. Commentaires sur les courbes
Si m ∈ [0, 1[ alors il y a un dépassement de gain
Si m ∈ [1,∞ [ alors il n’y a plus de dépassement de gain et dans nous n’observons aucune oscillation.
C. Tableau récapitulatif des valeurs du dépassement D%
s (t )−S (∞ )
D %= ×100
s (∞ )
Pour m ∈ ]0, 1[ (est proche de 1) , le temps de stabilisation devient faible alors la réponse du système devient plus rapide.
Pour m ∈ [1, ∞[ (est proche de ∞) , le temps de stabilisation devient grand alors la réponse du système devient plu lente .
D. Tableau récapitulatif et schémas du temps de réponse (tr) du système pour divers cas à 5%
3
tr5% = mw ; n
temps de réponse à 5% = [ 95% , 105% ]
2
XIV. Code h1 :
2
Pour la fonction de transfert h2
XVI. Code h2 :
2
Pour la fonction de transfert h3
XVIII. Code 3 :
2
Pour la fonction de transfert h4
XX. Code 4 :
2
Pour la fonction de transfert h5
XXII. Code h5 :
Réponse indicielle de la
fonction de transfert h1
2
Réponse indicielle de la fonction de transfert h2
2
Réponse indicielle de la fonction de transfert h4
2
Réponse indicielle des 04 fonctions de transfert superposées
2
F. THÉORÈME
Codes :
2
2
2
Diagrammes de Bode, Nyquist et Nicols :
2
G. INTERPRÉTATION DES DIFFÉRENTS DIAGRAMMES
Les diagrammes permettent d’étudier la stabilité, la précision, la rapidité et la marge de stabilité d’un système asservi.
2
XXV. Figure 1 :
Fonction de transfert he2
XXVI. Code 2 :
2
Fonction de transfert he3
XXVIII. Code 3 :
2
Pour la fonction de transfert he4
XXX. Code 4 :
2
Pour la fonction de transfert he5
XXXII. Code 5 :
2
Fonction de transfert he6
XXXIV. Code 6 :
2
Nb : Plus la période augmente, plus il devient continu.
y ( z) 1
H ( z)
= 2
z −0.3 z +0.5
XXXVI. Code :
2
XXXVII. Figure :
2
Coefficient d’amortissement Pulsation propre (Wn) Fonction de transfert du
(m) second ordre en continue
iii. Traçage des réponses indicielles unitaires des deux fonctions de transfert en Z et en S
2
2
2
2
2
2
2
2
iv. Calcul approximatif du temps de réponse, du temps de stabilisation à ±5% et le dépassement
3 S ( Am ) −S ( ∞ )
Tr5% = m× wn et D = s (∞)
Période 0.05
d’échantillonnage
Dépassement 45%
Temps de réponse à 0.338s
±5%
Temps de stabilisation 0.624s