Modélisation de base des réseaux
électriques
• Représentation des réseaux.
• Théorie
Thé i desd graphes
h appliquée
li é aux réseaux
é él t i
électriques,
formation des matrices admittance et impédance d’un RE
• Modification et inversion de la matrice admittance
2
Représentation des Réseaux
3
Exemple de diagramme unifilaire
4
Diagramme des impédances
5
Représentation des Réseaux Electriques
NOEUD ou Jeu
de barre (BUS)
BRANCHES
(Lignes ou
Transformateurs)
RESEAUX
(à vide non
alimenté)
6
Représentation des Réseaux Electriques
Charge Generateur
Reseau
(Alimenté et
Chargé)
7
Représentation des Réseaux Electriques
Modele
approximatif
Modele d’une
Branche
Reseau
(Alimenté et
chargé)
8
Représentation des Réseaux Electriques
On représente le réseau par des impedances (admittance)
L’impedance
L impedance est un nombre complexe zij=rrij+jx
jxij.
On peut negliger la resistance : zij=jxij
z12
1 2
z13
z14
z23
4 3
z34
9
Représentation des Réseaux Electriques
Selon la loi d’ Ohm:
1
I ij (Vi V j )
Courant (amps)
zij
Tension (volts)
Vi Iij
Vj
zij
i j
10
Représentation des Réseaux Electriques
1
I ij (Vi V j )
zij
L’admittance,, yij, est l’inverse de l’impedance,
p , zij:
I ijj yijj (Vi V j )
Vi Iij
Vj
yij
i j
11
Représentation des Réseaux Electriques
On peut représenter le réseau par les admittances yij
Vi : la
l tension
t i au nœud
d i.
i
Ii : le courant injecté dans le nœud i
I1 y12 I2
1 2
V1 V2
y13
y14
y23
V4
V3
4 3
y34 I3
I4
12
Théorie des graphes
Graphe: un dessin avec des points définis appelés noeuds (aussi sommets)
reliés par des éléments appelés segments.
Un graphe fini G = (A, E) est défini par l’ensemble fini A = {a1, a2,..., an}
dont les éléments sont les sommets (noeuds), et par l’ensemble fini
E = {e1,
{e1 e2,...,em}
e2 em} dont les éléments sont les segments.
segments
Un noeud de référence doit être choisi et il est pratique courante
d'assigner à ce noeud le numéro zéro (0)
Graphe orienté: un graph dans lequel on assigne une direction (sens) dans
chaque segment.
13
Exemple
14
Théorie des graphes
Incidence: Un segment (branche) est incident à un sommet
(nœud) si il « touche » ce sommet.
Arbre:
b graphique contenant tous les noeuds, mais pas de maille, il est
équivalent à un graphe à n noeud et n‐1 branches. Les segments
restants sont appelés les liens. LL’ensemble
ensemble des liens forme le co
co‐arbre.
arbre.
15
Représentation matricielle des réseaux
Le réseau peut être décrit par les trois types de matrices suivantes:
1‐ Matrices d'incidence: Elles représentent le lien entre les
branches et les nœuds du réseau avec leur orientation. Les valeurs
de ces matrices sont binaires 1,
1 0,
0 ‐ 1.
1
2‐ Matrices élémentaires (ou primitives): ces matrices décrivent
les composants individuels en tenant compte, le cas échéant, de
leur couplage électromagnétique (capacitif et inductif). Elles sont
de structure diagonale, le couplage des composants est représenté
par des éléments non diagonaux.
3‐ Matrices de transferts: elles décrivent de façon mathématique
le comportement électrique du réseau maillé. maillé Elles sont
essentiellement les matrices impédance ou admittance qui
correspondent aux nœuds du réseau (matrices nodales).
16
1‐ Matrice Incidence nodale
La matrice incidence nodale A décrit la topologie du système en
faisant abstraction de la nature des éléments qui le composent.
P
Pour un graphe h comportant n noeuds d et e segments (branches),
(b h ) la
l
matrice A est formé de n rangées et e colonnes. Ses éléments aij
peuvent prendre les valeurs 0, 0 ‐ 1 ou 1 comme suit :
•aij = 1 si la branche sort du noeud;
•aij = ‐1 si la branche rentre dans le noeud
•aij = 0 si la branche n’est pas connectée au noeud.
Remarques importante:
L matrice
La t i d'incidence
d'i id nodale
d l n'est
' t pas vraiment
i t utile
til dans
d l
le
développement qui suit. On utilise plutôt la matrice incidence
nodale réduite,
réduite qui correspond à la matrice 'incidence
incidence nodale,
nodale
dans laquelle le choix d'un nœud de référence des tensions
conduit à la suppression d'une colonne de la matrice (A)
17
Exemple:
Nœud de
référence 18
Exemple:
Nombre de Branches = 07
Nombre de
neouds = 04
19
2‐ Matrices élémentaires (ou primitives):
Soit v1,v2,v3,v4,v5,v6 et v7 les tensions aux bornes des éléments
du réseau (tensions des branches) et i1, i2, i3, i4, i5,i6,i7 les
courant qui traversent ces éléments.
Z
20
2‐ Matrices élémentaires (ou primitives):
Z estt appelé
lé la
l matrice
t i impédance
i éd d réseau
du é primitif
i itif , c’est
’ t une
matrice à dominance diagonale dont les éléments correspondent
aux impédances de chaque élément du réseau.
réseau Ces impédances sont
les impédances propres Zii et les impédances de couplage entre les
éléments (mutuelle) qui représentent les éléments hors diagonale Zij
De même sachant que:
On peut définir la matrice admittance du réseau primitif
tel que
21
3 ‐Les Matrices Zbus et Y bus
Soit V1, V2, V3 et V4 les tensions aux nœuds (mesurés par
rapport
pp à un nœud de référence nœud 0 ) et I1,I2,I3,I4
, , , les
courants injectés aux nœuds. Les tensions des branches sont
liées aux tension des nœuds par la relation suivante:
22
3‐ Les Matrices Zbus et Y bus
Les courants des branches sont liées aux courants injectés aux
nœuds par les relations suivantes:
23
Les Matrices Zbus et Y bus
Sous forme matricielle
24
Les Matrices Zbus et Y bus
Sous forme matricielle
25
Les Matrices Zbus et Y bus
De la même façon , pour un réseau général de N noeuds, on pourra écrire
Equations
i nodale
d l d
du
réseau
Zbus désigne
Zb dé i lla matrice
i d’impedance
d’i d d b
des barres ((nœuds)
d )
Ybus désigne la matrice d’admittance des barres (nœuds)
Tandis que Vbus et Ibus sont respectivement les matrices tensions et
courants aux N noeuds (jeux de barres) du réseau.
26
Formation de la matrice admittance nodale [Ybus]:
On a trois méthodes pour la détermination de la
matrice [Ybus] qui sont :
• a) Par inversion de la matrice impédance nodale
[Zbus] .
• b) à partir des admittances des éléments du réseau
((matrice p primitive)) et la matrice d’incidence
(topologique)
• c) Méthode d d’inspection
inspection
27
Formation de [Ybus] à partir des admittances des
éléments du réseau et la matrice d’incidence
28
Formation de [Ybus] à partir des admittances des
éléments du réseau et la matrice d’incidence
Dans le cas où le réseau ne présente pas de couplage mutuel
entre ses éléments , la matrice admittance du réseau primitif Y
sera diagonale
di l
29
Formation de [Ybus] à partir des admittances des
éléments du réseau et la matrice d’incidence
Par conséquent , la matrice Ybus aura la forme suivante
30
Formation de [Ybus] à partir des admittances des
éléments du réseau et la matrice d’incidence
31
Formation de [Ybus] à partir des admittances des
éléments du réseau et la matrice d’incidence
32
Formation de [Ybus] par la Methode d’inspection
En Résumé : pour un RE donné la matrice
d’admittance [Ybus] peut être construite comme
suit :
1‐ Les éléments diagonaux pour chaque noeud
sont la somme des admittances connectées en ce
noeud.
2‐ Les éléments non diagonaux sont donnés par
ll'admittance
admittance de la ligne joignant le point i à j
affectée du signe négatif
3 La matrice est symétrique Yij=Yji
3‐
33
Exemple
34
Formation de [Ybus] par la Methode d’inspection
D
Demonstationi
Soit le reseau suivant , selon Kirchoff
I1 I12 I13 I14
I1 y12 I2
1 2
I12 V2
V1
I14 y13
y14 I13
y23
V4
V3
4 3
y34 I3
I4 35
Formation de [Ybus] par la Methode d’inspection
I1 I12 I133 I14
En utilisant la loi d
d’Ohm’s
Ohm s :
I ij yij (Vi V j )
I1 y12 (V1 V2 ) y13 (V1 V 3) y14 (V1 V4 )
I1 V1 ( y12 y13 y14 ) V2 ( y12 ) V 3( y13 ) V4 ( y14 )
36
Pour les autres courants:
I1 y12 I2
1 2
I12 V2
V1
I14 y13
y14 I13
y23
V4
V3
4 3
y34 I3
I4
I 2 V1 ( y21 ) V2 ( y21 y23 y24 ) V 3( y23 ) V4 ( y24 )
I 3 V1 ( y31 ) V2 ( y32 ) V 3( y31 y32 y34 ) V4 ( y34 )
I 4 V1 ( y41 ) V2 ( y42 ) V 3( y41 y42 y43 ) V4 ( y43 )
37
pour les autres noeuds:
I1 V1 ( y12 y13 y14 ) V2 ( y12 ) V 3( y13 ) V4 ( y14 )
I 2 V1 ( y21 ) V2 ( y21 y23 y24 ) V 3( y23 ) V4 ( y24 )
I 3 V1 ( y31 ) V2 ( y32 ) V 3( y31 y32 y34 ) V4 ( y34 )
I 4 V1 ( y41 ) V2 ( y42 ) V 3( y41 y42 y43 ) V4 ( y43 )
Notes:
I1 y12 I2 1. yij=yyji
1 2 2. Si la branche ij
I12 V2 n’existe pas, alors
V1
yij=0.
I14 y13
y14 I13
y23
V4
V3
4 3
y34 I3
I4 38
Sous forme matricielle:
I1 y12 y13 y14 y12 y13 y14 V1
I y y y y y y V
2 21 21 23 24 23 24 2
I3 y31 y32 y31 y32 y34 y34 V3
I4 y41 y42 y43 y41 y42 y43V4
I 1 Y11 Y12 Y13 Y14 V1
I Y Y22 Y23 Y24 V2
2 21
I 3 Y31 Y32 Y33 Y34 V3
I 4 Y41 Y42 Y43 Y44 V4
39
Algorithme de Construction de [Ybus]
Pour un RE donné la matrice d’admittance peut être
construite par en utilisant l’équation suivante
Avec: Yb = jB est l’admittance
’ shunt
B: est la susceptance shunt.
n: nombre
b ded bus
b
40
41
Exemple
42
Algorithme de Construction de [Ybus]
Pour un RE donné la matrice d’admittance peut être
construite par l’addition des éléments un par un.
En séparant le lien correspondant à l’élement 5 qui relie les
nœuds 2 et 3, la matrice Ybus peut être réécrite sous la
forme suivante :
43
Algorithme de Construction de [Ybus]
44
Algorithme de Construction de [Ybus]
On peut montrer que l’effet d’ajouter un élément 5 entre les
nœuds 2 et 3 est:
l’ajout d’une admittance Y55 aux éléments :
Ybus(2,2) et Ybus(3,3)
L’ j
L’ajout d’
d’une admittance
d i – Y55 aux éléments:
élé
Ybus(2,3) et Ybus (3,2).
Pour construire la matrice Ybus, on initialise tous les éléments à
la valeur zero, ensuite on ajoute les éléments du réseau un par
un, pour chaque pas 04 éléments de Ybus sont modifiés (à
démontrer))
45
Modélisation de base des réseaux
électriques
Construction de la matrice Z bus.
Algorithme
g de construction de la matrcie Zbus
La modification de la matrice Zbus implique l'un des 4 cas suivants
Cas 1 : Ajout d'une nouvelle branche d’impedance Zb entre un
nouveau bus (k) et le bus de référence (r)
Cas 2 : Ajout d’une nouvelle branche d'impédance Zb entre un
nouveau bus ((k)) et un bus existant (j)
Cas 3 : Ajout d'une nouvelle branche d’impedance Zb entre un
bus existant (j) et le bus de référence (r)
Cas 4 : Ajout d'une nouvelle branche d’impedance Zb entre
deux bus existants (i) et (j)
Algorithme
g de construction de la matrcie Zbus
Cas 1 : Ajout d'une nouvelle branche d’impedance Zb entre un
nouveau bus (k) et le bus de référence (r)
L'ajout d'un nouveau bus augmentera l'ordre de la matrice Zbus de 1
Algorithme
g de construction de la matrcie Zbus
Cas 2 : Ajout d’une nouvelle branche d'impédance Zb entre un
nouveau bus (k) et un bus existant (j)
Algorithme
g de construction de la matrcie Zbus
Cas 2 : Ajout d’une nouvelle branche d'impédance Zb entre un
nouveau bus (k) et un bus existant (j)
Algorithme
g de construction de la matrcie Zbus
Cas 3 : Ajout d'une nouvelle branche d’impedance Zb entre un bus
existant (j) et le bus de référence (r)
Algorithme
g de construction de la matrcie Zbus
Cas 3 : Ajout d'une nouvelle branche d’impedance Zb entre un bus
existant (j) et le bus de référence (r)
Algorithme
g de construction de la matrcie Zbus
Cas 4 : Ajout d'une nouvelle branche d’impedance Zb entre deux bus
existants (i) et (j)
Algorithme
g de construction de la matrcie Zbus
Cas 4 : Ajout d'une nouvelle branche d’impedance Zb entre deux bus
existants (i) et (j)
Exemple : Formation de Zbus
Etape
p 1 : Ajout
j d'une nouvelle branche d’impedance
p Z1 entre
un nouveau bus (1) et le bus de référence (j)
Etape
p 2 : Ajout
j d'une nouvelle branche d’impedance
p Z2 entre
un nouveau bus (2) et le bus de référence (j)
Etape 3 : Ajout d'une nouvelle branche d’impedance Z4 entre
un nouveau bus (3) et un bus existant (2)
Etape 4 : Ajout d'une nouvelle branche d’impedance Z3 entre
deux bus existants (1) et (2)
Etape 5 : Ajout d'une nouvelle branche d’impedance Z5 entre
deux bus existants (1) et (3)
Devoir maison !!!
Modification de la matrice admittance
( é h d de
(méthode d réduction
éd i de d Kron
K )
Soit le réseau à quatre nœuds définit par l’équation suivante
Si le courant I4 =0,
=0 le nœud 4 peut être supprimé et l’équation
l équation
du réseau devient: (Determiner les nouvelles admittances)
61
Modification de la matrice admittance
( é h d de
(méthode d réduction
éd i de d Kron
K )
L’équation du nœud 4 est de la forme
Remplaçant V4 dans ll’équation
équation du nœud 1 :
62
Modification de la matrice admittance
( é h d de
(méthode d réduction
éd i de d Kron
K )
En comparant avec
63
Modification de la matrice admittance
( é h d de
(méthode d réduction
éd i de d Kron
K )
Conclusion , quand le nœud k est supprimé, les éléments
modifiés de Ybus peuvent être calculer comme suit:
64
Modification de la matrice admittance
( é h d de
(méthode d réduction
éd i de d Kron
K )
Exemple:
Calculer les tension V1, V2, V3 et V4 en utilisant la technique
de réduction du réseau
65
solution:
Eliminant le noeud 1
Eliminant le noeud 3
66
