Université Kasdi Merbah Ouargla 2022/2023
Faculté des Sciences Appliquées
Département De Génie Électrique
Master 02 Génie Électrique Module: TP Stabilité et dynamique des réseaux électriques
Travail pratique no. 3
ANALYSE DE STABILITÉ TRANSITOIRE ET PETITE PERTURPATION– SYSTEME
MACHINE A BUS INFINITE
1. OBJECTIF
-Se familiariser avec divers aspects de l'analyse de la stabilité transitoires et de petites perturbations du
système SMIB (Single-Machine-Infinite Bus)
- Etudier et simuler à l'aide d'un logiciel de simulation l'analyse de la stabilité transitoires et des petites
perturbations du système SMIB (Single-Machine-Infinite Bus)
2. THÉORIE: Stabilité
Le problème de stabilité concerne le comportement du système électrique lorsqu'il est soumis à des
perturbations et est classé en petit problème de stabilité du signal si les perturbations sont petites et en
problème de stabilité transitoire lorsque les perturbations sont importantes.
La stabilité transitoire: Lorsqu'un système électro-énergétique est en régime permanent, la charge plus la
perte de transmission est égale à la production dans le système. La machine synchrone fonctionne à une
vitesse de rotor synchrone et la fréquence, la tension, le courant et la puissance du système sont stables.
FORMULE
Pm puissance mécanique générée en p.u.
Pe puissance électrique générée en p.u.
E = f.e.m Générateur en p.u.
V = La tension de bus infini en p.u.
X1 =Réactance avant défaut en p.u.
X2 =Réactance pendant défaut en p.u.
X3 =Réactance après défaut en p.u.
H = Constant d’inertie de générateur en sec.
f = Fréquence du systeme en Hz
Dt = Intervale du temps
tc = le temps d’élimination du défaut en sec
tf = Temps final de equation swing en sec
Puissance réactive Qe = sin(cos-1(p.f))
S*
Courant du stator It =
Et*
Pe - jQe
=
Et*
Tension derrière une condition transitoire
E1 = Et + j Xd1It
Tension du bus infini
EB = Et - j( X3 + Xtr )It
X1 X2
où, X3 =
X1 + X2
1
Angle entre E et EB
δo = ∠ E1 - ∠ EB
1) Condition avant défaut:
X1 X2
1
X = j Xd + jXtr +
X1 + X2
1
Dr. LAROUCI Benyekhlef
E1 x EB
Puissance: Pe = sinδo
X
Pe * X
δo = sin -1
E1 * EB
2) Condition Pendant du défaut: Pe = PEii = 0
Trouvez X à partir du circuit équivalent pendant la condition de défaut
3) Condition après défaut:
Trouvez X à partir du circuit équivalent pendant la condition l’application du défaut
E1 x EB
Power Pe = sinδo
X
δmax = π - δo
Pm
Pe =
sinδmax
a) Angle Critique après élimination du défaut
δmax - δo ) + P3maxcosδ
Pm(δ δmax - P2maxcosδ
δo
Cosδcr =
P3max - P2max
Pe1max= E*V/X1 : puissance électrique max générée avant du défaut en p.u
Pe2max=E*V/X2 : puissance électrique max générée pendant du défaut en p.u
Pe3max= E*V/X3 : puissance électrique max générée après élimination du défaut en p.u
a) Temps Critique après élimination du défaut
2H (δcr - δo)
tcr = Secs
√ πfo Pm
3. Illustration:
Un générateur synchrone à 60 Hz ayant une constante d’inertie H = 5MJ / MVA et une réactance transitoire
à axe direct Xd ’= 0,3p.u. est connecté à un bus infini via un circuit purement réactif comme indiqué sur la
figure 1. Les réactances sont représentées sur le schéma dans une base de système commune. Le générateur
délivre une puissance réactive Pe = 0,8p.u. et Q = 0,074 p.u. au bus infini à une tension de 1p.u. Un défaut
triphasé se produit au milieu d'une ligne et est éliminé en isolant le circuit défectueux simultanément aux
deux extrémités, comme indiqué sur la figure 4.1. Le défaut est effacé en 0,3 seconde. La solution
numérique est obtenue pendant 1,0 seconde
dans Matlab7.0. La courbe d'oscillation est utilisée pour déterminer la stabilité du système et le temps de
compensation critique est déterminé. La simulation a été répétée et les swing plots obtenus à l'aide de
SIMULINK.
4.1 CODE MATLAB POUR LA CONCEPTION DE LA STABILITÉ DES ÉTATS TRANSITOIRES
Pm = 0.80; E = 1.17; V = 1.0;
X1 = 0.65; X2 = 1.80; X3 = 0.8;
H = 5.0; f = 60; tc = 0.3; le temps d’élimination du défaut 0.3 sec.
2
Dr. LAROUCI Benyekhlef
a) Un défaut triphasé temporaire se produit à l'extrémité émettrice de la ligne au point F. Lorsque le
défaut est éliminé, les deux lignes sont intactes. Déterminez l'angle d'élimination critique et le temps
d'élimination du défaut critique :
Pe1max= E*V/X1 = 1,80= Pm= P1max * sinδmax=1,80*sinδmax
Pe2max=E*V/X2 = Pemax = 0,650 Pm= P2max * sinδmax=0,65*sinδmax
Pe3max= E*V/X3 = 1,4625 Pm= P3max * sinδmax=1,4625*sinδmax
Pemax=Pe2max
ck=pi*f/H
d0 =asin(Pm/Pe1max);
dmax = pi-asin(Pm/Pe3max);
cosdc = (Pm*(dmax-d0)+Pe3max*cos(dmax)-Pe2max*cos(d0))/(Pe3max-Pe2max);
dc = acos(cosdc);
Criticalclearingangle = dc*180/pi
Criticalclearingtime=2*(sqrt(2*H*( dc- d0)/(2*pi*f*Pm)))
b) Faite, vérifiez et interpréter le résultat à l'aide du programme Simulink MATLAB :
Pm=0.8 Radian à degree
X'2 X2 X1
1 1 delta
pi*60/5 180/pi
s s
Step
Subtract Integrator Integrator1 delta
Gain Gain1
Aavant le défaut Pe1
1.8*sin(u)
Pm=0.8 Radi an à degree
X'2 X2 X1
1 1 del ta
pi *60/5 180/pi
s s
Step
Subtract Integrator Integrator1 del ta
Gai n Gain1
Pendant l e défaut Pe2
0.65*si n(u)
Pm=0.8 Radi an à degree
X'2 X2 X1
1 1 del ta
pi *60/5 180/pi
s s
Step
Subtract Integrator Integrator1 del ta
Gai n Gain1
Apres le défaut Pe2
1.4625*sin(u)
3
Dr. LAROUCI Benyekhlef
Partie II. Programme MATLAB pour résoudre l'équation d’oscillation à l'aide de la méthode point
par point
Objectif : Résoudre l'équation d’oscillation (swing) du problème donné en utilisant la méthode point par
point et écrire un programme MATLAB pour vérifier le résultat.
PROBLÈME:
Un générateur de 20 MVA, 50 Hz délivre 18 MW sur une ligne à double circuit vers un bus infini. Le
générateur a un KE de 2,52MJ/MVA à la vitesse nominale. La réactance transitoire du générateur est
Xd=0.35p.u. Chaque ligne de transport a R=0 et une réactance de 0.2pu sur 20 MVA Base. |E|=1,1 p.u et
tension de bus infinie V=1,0. Un court-circuit triphasé se produit au milieu de l'une des lignes de transport.
Tracer les courbes d'oscillation avec défaut éliminé par ouverture simultanée des disjoncteurs aux deux
extrémités de la ligne à 6,25 cycles après l'apparition du défaut. Tracez également la courbe d'oscillation
sur la période de 0,5 s si le défaut persiste.
Résolvez théoriquement l'équation du swing par la méthode point par point et vérifiez à l'aide du
programme MATLAB.
Commenter sur les résultats de la stabilité du système.
PROGRAMME MATLAB :
%Defining the function swing
function[time ang]=swing(tc)
k=0;v=1;E=1.1;pm=0.9;T=0.5;delT=0.05;ddelta=0;time(1)=0;ang(1)=21.64; xdf=1.25;xaf=0.55;t=0;
delta=21.64*pi/180;i=2; m=2.52/(180*50);
while t<T
if t<tc
x=xdf;
else x=xaf;
end
pmax=(E*v)/x;
pa=pm-pmax*sin(delta);
ddelta=ddelta+(delT^2*(pa/m));
delta=(delta*180/pi+ddelta)*(pi/180);
deltadeg=delta*180/pi;
t=t+delT;
time(i)=t;
ang(i)=deltadeg;
i=i+1;
end
end
Program 2: Main program that is dependent on swing.m
clc
clear all
close all
for i=1 : 2
tc=input('enter the value of clearing time:\n');
[time,ang]=swing(tc)
t(:,1)=time;
a(:,i)=ang;
end
plot(t,a(:,1),'*-',t,a(:,2),'d-')
axis([0 0.5 0 inf])
t,a
Enter the value of clearing time as 0.25 sec, and 5 sec
- Cnclusion:
4
Dr. LAROUCI Benyekhlef
Expt 7 . - ANALYSE DE LA STABILITÉ TRANSITOIRE DE LA MACHINE DE BUS
INFINI (SMIB)
OBJECTIF:
Écrire un programme MATLAB pour analyser la stabilité transitoire d'une e machine à bus infinie,
en supposant le modèle classique du générateur (tension constante derrière réactance transitoire)
THÉORIE:
La stabilité du système électrique peut être définie au sens large comme la propriété d'un système
électrique qui permet de rester dans un état d'équilibre de fonctionnement dans des conditions
normales de fonctionnement et de regagner un état d'équilibre acceptable après avoir été soumis à
une perturbation
La stabilité du système d'alimentation peut être classée en gros comme (i) la stabilité de l'angle du
rotor et (ii) la tension stabilité. La stabilité de l'angle du rotor est la capacité des machines
synchrones interconnectées d'une puissance système de rester en synchronisme.
La stabilité de l'angle du rotor peut être classée en stabilité transitoire et stabilité de petit signal
selon le type de perturbation. La stabilité transitoire est l'étude de la stabilité de l'angle du rotor
d'un système à la suite de perturbations importantes.
La stabilité des petits signaux (ou des petites perturbations) est la capacité du système électrique à
maintenir synchronisme sous de petites perturbations. Les perturbations sont considérées comme
suffisamment faibles pour la linéarisation des équations du système doit être admissible aux fins de
l'analyse. L'instabilité qui peut résultat peut être de deux formes.
I. Augmentation constante de l'angle du rotor en raison d'un manque de couple de synchronisation
suffisant.
II. Oscillations du rotor d'amplitude croissante dues à un manque de couple d'amortissement
suffisant. Là
sont quatre modes d'oscillations provoquant une petite instabilité du signal dans un système
d'alimentation. Elles sont:
• Les modes locaux ou les modes système de la machine sont associés au basculement des unités à
un
centrale par rapport au reste du système électrique. La gamme de fréquence de
l'oscillation est (0,8 à 2) Hz.
• Les modes inter zones sont associés au basculement de plusieurs machines dans une partie du
système contre des machines dans d'autres parties. La gamme de fréquences pour les modes
interzones est (0,2 à
0,8) Hz.
• Les modes de contrôle sont associés aux groupes électrogènes et à d'autres contrôles.
• Les modes de torsion sont associés à la rotation du système d'arbre turbine-alternateur
Composants
Exemple numérique :
Le système de test considéré pour l'analyse de la stabilité transitoire est le bus infini à une seule
machine
5
Dr. LAROUCI Benyekhlef
Fig. 1. Schéma unifilaire du système SMIB
La figure montre la représentation du système applicable à la centrale thermique se composent
de quatre unités 555MVA, 24kV, 60Hz.
Les réactances du réseau indiquées dans la figure sont exprimées par unité sur une base de
2220MVA, 24kV (référence côté basse du transformateur élévateur). Les résistances sont
supposées être négligeables.
L'état du système après défaut par unité sur la base 2220MVA, 24kV est comme suit
P=0.9 Q=0.3 (surexcité ) Et=1.0 36° EB=0.995∟0°
Les générateurs doivent être modélisés comme un seul générateur équivalent représenté par le
modèle classique avec les paramètres suivants exprimés en PU de 2220 MVA 24kV base
Xd=0.3 H=3.5MW-s/MVA
Programme:
p=0.9;q=0.3;et=1.0;anget=0.628;eb=0.995;angeb=0;xd=.3*i;f=60;xtr=.15*i;x1=.5*i;x2=.93*i;
h=input('enter the value of inertia:');
fault=input('enter the line no:');
if fault==1
xt=xd+xtr+x2;
else
xt=xd+xtr+x1;
end
it=conj((p+q*i)/et);
edas=et+(xd*it);
angedas=angle(edas);
delta=anget+angeb+angedas;
ks=(edas*eb*cos(delta))/xt;
kd=input('enter the value of kd:');
A=[-kd/(2*h) -abs(ks)/(2*h); (2*pi*f) 0 ]
lamda=eig(A)
[v d]=eig(A)
l=inv(v)
p=abs(v.*1)
s2=det(A);
wn=sqrt(s2);
f=sqrt(s2)/(2*pi)
s1=-(-kd/(2*h))+0;
zita=s1/(2*wn)
wd=wn*(sqrt(1-(zita*zita)));
wdh=wd/(2*pi)
%for plots using zero inputresponse
th=acos(zita);
t=0:0.01:3; dd=anget/sqrt(1-zita^2)*exp(-zita*wn*t).*sin(wd*t+th); d=(delta+dd)*180/pi;
%relative to synchronously revolving field
dw=-wn*anget/sqrt(1-zita^2)*exp(-zita*wn*t).*sin(wd*t);
ff=f+dw;
6
Dr. LAROUCI Benyekhlef
subplot(2,1,1),plot(t,d),grid
xlabel('t sec'),ylabel('delta degre')
subplot(2,1,2),plot(t,ff),grid
xlabel('t sec'),ylabel('rotor speed ')
Résultat: Cas-1 : Après l'élimination du défaut sur la ligne-1
Constant KD Valeurs propres ωn (rad/sec) rapport ωn
inertie, H amortissement (Hertz)
0
3.5 10
-10
70 10
-10
Cas-1 : Après l'élimination du défaut sur la ligne-2
Constant KD Valeurs propres ωn (rad/sec) rapport ωn
inertie, H amortissement (Hertz)
0
3.5 10
-10
70 10
-10
RÉSULTATS:
Une augmentation du coefficient de couple d'amortissement KD augmente le taux
d'amortissement, alors que une augmentation de la constante d'inertie diminue à la fois ωn et ?
7
Dr. LAROUCI Benyekhlef