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
|
fid=fopen('C:\EEUG07\ATPDraw\Atp\total_param.txt','w');
for C1=18E-9 : 10E-9 : 28E-9
%Calcul de PCC
Pcc=390E6
Ueff=6.6e3
rapport=12.5;%Ration X/R
Zcc=(Ueff^2)/Pcc;
Xcc=sqrt(Zcc^2/(1+1/rapport^2));
Lcc=Xcc/(100*pi)
Rcc=Xcc/rapport
%Calcul du moteur sous une charge RL avec les paramètres du schéma équivalent
% Rfe=
% Xm=
% R1=
% X1=
% R2=
% X2=
% Zfer=(Rfe*i*Xm)/(Rfe+i*Xm)
% Zrotor=R2+i*X2
% Zstator=R1+i*X1
% Zeq=Zstator+((Zfer*Zrotor)/(Zfer+Zrotor))
%Calcul du moteur sous une charge RL sans les paramètres du schéma équivalent
%data
Pn=2750E3;
cosphin=0.9;
cosphid=0.12;
K=6.5;
U=6600;
%Calculs
Sn=Pn/cosphin;
Sd=Sn*K;
ZM=U^2/Sd;
RM=ZM*cosphid
XM=sqrt(ZM^2-RM^2);
LM=XM/(100*pi)
CM=1/(4*pi^2*14000^2*LM)
%Calcul de la fréquence de résonance en fonction des capas et inductance
%C1=18.44e-9;
CAV=216.7E-9;
C2=CAV+CM;
Ceq=C1*C2/(C1+C2);
L1=3.1E-6;
L2=20.46E-6;
Leq=L1+L2;
F=1/(2*pi*sqrt(Ceq*Leq))
%Calcul de la fréquence des oscillations en fonction de la PCC et des capas
Cosc=C1+C2;
Fosc=1/(2*pi*sqrt(Cosc*Lcc))
%Calcul de la résistance d'amortissement des inductances
K=40;
Ram1=K*2*pi*F*L1
Ram2=K*2*pi*F*L2
Rccam=2*pi*Fosc*Lcc
% %Calcul de la résistance en fonction de la fréquence en tenant compte de l'effet de peau
% %F=120000
Section=120;%en mm²
Rayon=sqrt(Section/pi)*1E-3;
%ro=3E-8;%Pour un cable aluminium
ro=1.8E-8%Pour un cable cuivre
e=1/sqrt((pi*4*pi*1E-7*F)/ro);
Surf=pi*Rayon^2-pi*(Rayon-e)^2;
Res1=ro*20/(Surf*2)%Surf*2 car 2 cables en parralèles
Res2=ro*235/(Surf*2)%Surf*2 car 2 cables en parralèles
%Tracé de la capa equivalente en fonction de deux capacités
%for C1 = 1 : 200,
%Ceq(C1)=C1*176/(C1+176);
%end
%plot(Ceq)
%Mise à l'échelle ATP des capacités et inductances
Lcc=Lcc*1E3;
L1=L1*1E3;
L2=L2*1E3;
LM=LM*1E3;
C1=C1*1E6;
CAV=CAV*1E6;
CM=CM*1E6;
fprintf(fid,'%-2.12f %-8.8f %-8.8f %-2.12f %-2.12f %-8.8f %-8.8f %-2.12f %-8.8f %-8.8f %-2.12f \n\r',Lcc,Rcc,Rccam,C1,L1,Res1,Ram1,L2,Res2,Ram2,CAV);
end
fclose(fid); |
Partager