%% Machine Asynchrone

% v01-v4pdf : code versioning pour suivi des versions A.Hodder
% v02-v5pdf : enlevé la partie freinage et rajouté courant au démarrage
% conversion = MAS2
close all; clear; clc  % on ferme tout et clear tout

%% Préambule
% Le but de cet exercice est multiple :
% Faire le passage des p.u. aux vraies grandeurs pour les paramètres du schéma équivalent.
% Calculer l'impédance équivalente de la machine asynchrone.
% Appliquer l’équivalent de Thévenin pour le calcul du couple et du glissement critique (ceci pourrait évidemment être étendu au calcul du couple).
%% Exercice
% Une machine asynchrone triphasée a les caractéristiques suivantes :
% Stator couplé en étoile
% P_{fer} = 0 (R_{fer} est négligé dans le schéma équivalent)

Pn = 30e3      % [kW] - Puissance active nominale
Unligne_star = 400  % [V] - Tension de ligne nominale en étoile
Inligne_star = 50   % [V] - Courant de ligne nominal en étoile
fn = 50        % [Hz] - Fréquence
p = 2          % [-] - Nombre de paires pôles

rs_pu = 0.02   % [pu] - Résistance statorique en p.u.
xss_pu = 0.12  % [pu] - Réactance de fuites statorique en p.u.
xh_pu = 3      % [pu] - Réactance de champ principal en p.u.
rrp_pu = 0.03  % [pu] - Résistance rotorique en p.u.
xsrp_pu = 0.10 % [pu] - Réactance de fuites rotorique en p.u.

% Questions
% La machine est alimenté sous tension et fréquence nominales.
% 1. Que vaut le courant de ligne à l'instant du démarrage en p.u. ?
% 2. Pour un glissement s = 0.027, calculer
%         2.1 le courant de ligne au stator
%         2.2 la puissance électrique active consommée au stator
%             a. en utilisant le cos\varphi à partir de l’impédance équivalente
%             b. en calculant la puissance apparente au stator
% 3. Calculer le glissement critique sK et le couple de décrochage TK.
%
%% Résolution
% Tension et courant de phase nominaux

Un = % [V]
In = % [A]

% Impédance nominale

Zn = % [Ω]

% Paramètres en vraies grandeurs

Rs = % [Ω] - Résistance statorique
Xss = % [Ω] - Réactance de fuites statorique
Xh = % [Ω] - Réactance de champ principal
Rrp = % [Ω] - Résistance rotorique
Xsrp = % [Ω] - Réactance de fuites rotorique

%
%% 1. Courant de ligne à l'instant du démarrage en p.u.
% glissement à l'instant du démarrage

s = 1 % [-]

% Impédance équivalent complexe en fonction du glissement

Zeq_cmpx_start = % [Ω]

% Norme de l'impédance équivalente complexe

Zeq_start = % [Ω]

% Courant de ligne au démarrage

Istart = % [A]

% Courant de ligne au démarrage en p.u.

istart = % [pu]

%
%% 2. Courant de ligne et puissance active électrique au stator
%% 2.1 Courant de ligne au stator
% Tension de phase pour le cas étudié ici

Uph =  % [V]

% Impédance équivalente pour un glissement de 0.027.
% Glissement

s = 0.027 % [-]

% Impédance équivalent complexe (copy/paste car fait avant)

Zeq_cmpx = % [Ω]

% Norme de l'impédance équivalente complexe (copy/paste car fait avant)

Zeq = % [Ω]

% Courant de ligne (norme)

Iph = % [A]

%
%% 2.2.a  Puissance électrique au stator en utilisant le cos\varphi à partir de l’impédance équivalente
% cos \varphi

cosphi22a = % [-]

% Puissance électrique

Pel22a = % [W]

%
%% 2.2.b Puissance électrique au stator en calculant la puissance apparente
% Courant de ligne (complexe)

Iph_cmpx = % [A]

% Puissance apparente complexe

S_cmpx = % [VA]

% Puissance électrique

Pel22b = % [W]

%
%% 3. Glissement et le couple critique
% Tension de phase complexe

Uph_cmpx = % [V]

% Tension de l'équivalent de Thévenin (complexe)

Ue_cmpx = % [V]

% Tension de l'équivalent de Thévenin (norme)

Ue = % [V]

% Impédance équivalente de de Thévenin

Ze_cmpx = % [Ω]

% Partie réelle et imaginaire de l'impédance équivalente de l'équivalent de Thévenin

Re =  % [Ω]
Xe =  % [Ω]

% Glissement critique

sk = % [-]

% Vitesse angulaire du champ tournant statorique (dans le monde mécanique)

f = % [Hz]
Omega_s = % [rad/s]

% Couple critique

Tk = % [Nm]

%%
%
%
%
%
%
%
%
%
%
%
%