1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
| %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%DEBUT DE LA FONCTION PRINCIPALE%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function calcul_Icc0
%%%%initialisatio de x%%%%
x=0;
%%%%saisie des données%%%%
%Alternateur1
Una1=15.75e3;
sna1=167e6;
X1=0.12;
%transforamteur1
Unt1=247e3;
Ucc1=0.14;
snt1=180e6;
%ligne
X0=0.4;
L=150;
%transformateur2
Unt2=225e3;
Ucc2=0.115;
snt2=120e6;
%alternateur2
Una2=99e3;
sna2=500e6;
X2=0.14;
%calcul des relatives monophasés
E1=(Una1/sqrt(3))*(Unt1/Una1);
E2=(Una2*(Unt2/Una2))/sqrt(3);
% Création de l'objet Figure
handles(1)=figure('Position', [100 100 700 500]);
'color', [0.8 0.5 0.2]; ... % couleur RGB du fond
set(handle(1),'name','Calcul des courants de Court circuit'),
% puis créer l'axe en précisant sa position et dimension en pixels
axes('Units', 'pixels', 'Position', [80 500 700 200]);%axis off,
title('Trançon défectueux du réseau'),
%insertion des boutons poussoirs dans la fenêtre
handles(2)=uicontrol('Style', 'pushbutton', ... % de type Bouton pour afficher le trançon defectueux du réseau
'String', 'Réseau', ... % nom du bouton
'Position', [10 10 90 25], ... % position dans la fenetre et dimension dans la fenêtre f
'Callback',@reseau);
handles(3)=uicontrol('Style', 'pushbutton', ... % de type Bouton pour afficher le schéma électrique équivalent
'String', 'Défaut3~', ... % nom du bouton
'Position', [110 10 90 25], ... % position dans la fenetre et dimension dans la fenêtre f
'Callback', @defauttriphase);
handles(4)=uicontrol('Style', 'pushbutton', ... % de type Bouton pour afficher le schéma électrique équivalent
'String', 'Défaut2~', ... % nom du bouton
'Position', [310 10 90 25], ... % position dans la fenetre et dimension dans la fenêtre f
'Callback', @defautbiphase);
handles(5)=uicontrol('Style', 'pushbutton', ... % de type Bouton pour afficher le schéma électrique équivalent
'String', 'Défaut1~', ... % nom du bouton
'Position', [410 10 90 25], ... % position dans la fenetre et dimension dans la fenêtre f
'Callback', @defauttriphase);
%insertion de zone de texte dans la fenêtre
handles(6)=uicontrol('Style', 'text', ... % de type Bouton message pour demander là où on veut créer le défaut
'String', 'Emplacement du défaut en pourcentage en %', ... % le texte affiché dans le bouton
'Position', [30 175 280 25], ... % position dans la fenetre et dimension
'Callback','0');
%insertion d'éditeur de texte dans la fenêtre
handles(7)=uicontrol('Style', 'edit', ... % de type Bouton qui servira pour ecrire le pourcentage du défaut par raport à la ligne
'String', '',...
'Position',[330 175 90 25],...
'Callback',@defaut);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%DEBUT DE LA FONCTION IMBRIQUEE defaut%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function reseau(obj,event)
h=imread('res.jpg');
imshow(h);
end
function defaut(obj,event)
%calcul des impédances directes disp('Calcul des impédances directes ')
%calcul de Xalt1
Xd1=X1*(Unt1^2)/sna1;
%calcul de Xt1
Xd2=Ucc1*((Unt1^2)/snt1);
%calcul de Xligne
set(handle(7),'string',x)
Xd3=X0*L*x/100;
Xd4=X0*L*(1-(x/100));
%calcul de Xt2
Xd5=Ucc2*(Unt2^2)/snt2;
%calcul de Xalt2
Xd6=X2*(((Una2^2)/sna2)*((Unt2/Una2)^2));
%calcul des impédances inverses disp('les impédances inverses ')
%calcul de Xalt1
%on prend Zi=0,2Zd
Xi1=0.2*Xd1;
%calcul de Xt1
% Zi=Zd pour les transformateurs
Xi2=Xd2;
%calcul de Xligne
% Zi=Zd pour les lignes
Xi3=Xd3;
Xi4=Xd4;
%calcul de Xt2
% Zi=Zd pour les transformateurs
Xi5=Xd5;
%calcul de Xalt2
%on prend Zi=0.2Zd
Xi6=0.2*Xd6;
%calcul des impédances homopolaires disp('Calcul des impédances homopolaires ')
%calcul de Xalt1
%on prend Zo=0,05Zd
Xo1=0.05*Xd1;
%calcul de Xt1 homopolaire
% Zo= infini, pour le cas où les neutres des deux transformateurs sont isolés
%Xo2==Xo5==inf
% Zo=Zd, pour N lié à la terre d'1 transformateur
Xo2=Xd2;
Xo5=Xd5;
%calcul de Xligne
% on prend pour les lignes Zo=2Zd
Xo3=2*Xd3;
Xo4=2*Xd4;
%calcul de Xalt2
%on prend Zo=0.07Zd
Xo6=0.07*Xd6;
function calcul_icc0(obj,event) %#ok<*DEFNU>
function defauttriphase(obj,event)
%calcul des courants de défaut par la méthode de superposition
%Icc en défaut triphasé
Zd1=Xd1+Xd2+Xd3;
Zd2=Xd4+Xd5+Xd6;
Icc3=(E1/Zd1)+(E2/Zd2)
end
%Icc en défaut biphasé
function defautbiphase(obj,event)
Zi1=Xi1+Xi2+Xi3;
Zi2=Xi4+Xi5+Xi6;
Icc2=Unt1/(2*(Zd1+Zi1)) + Unt2/(2*(Zd2+Zi2))
end
%Icc en défaut monophasé
function defautmonophase(obj,event)
Icc11=0,
disp('Cas de défaut monophasé avec neutre directement lié à la terre ')
Zo1=Xo1+Xo2+Xo3;
Zo2=Xo4+Xo5+Xo6;
Icc12=Unt1*sqrt(3)/(Zd1+Zi1+Zo1) + Unt2*sqrt(3)/(Zd2+Zi2+Zo2),
disp('Cas de défaut monophasé avec neutre lié à la terre via une impédance Zn ')
Zn=40;
Zo11=Zo1+3*Zn;
Zo22=Zo2+3*Zn;
Icc13=Unt1*sqrt(3)/(Zd1+Zi1+Zo11) + Unt2*sqrt(3)/(Zd2+Zi2+Zo22),
end
end
end
end |
