LT TBD Word

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 26

clc;

clear all;
disp(' ||--------------------------------------------------------------------
---------------------------------||');
disp(' || MON HOC : THIET KE THIET BI DIEN
||');
disp(' || ****------*****------****
||');
disp(' || Bai thiet ke: THIET KE DONG CO 3 PHA ROTO LONG SOC
||');
disp(' ||
||');
disp(' || Sinh vien thuc hien: HOANG VU NGOC HUNG
||');
disp(' || Sinh vien lop: KY THUAT DIEN_DIEN TU K37_A
||');
disp(' ||--------------------------------------------------------------------
---------------------------------||');
disp(' ****_______(^.^)_______****'
);
disp(' Nhap cac thong so co ban:');
P=input('- Nhap cong suat cua dong co (KW): P= ');
while(P~=1.1 && P~=1.5 && P~=2.2 && P~=3 && P~=4 && P~=5.5 && P~=7.5 && P~=11
&& P~=15 && P~=18.5 && P~=22 && P~=30 && P~=37 && P~=45 && P~=55 && P~=75 &&
P~=90)
P=input('Ban da nhap sai gia tri P. Vui long nhap lai: P=');
end
pp=input('- Nhap so doi cuc : 2p= ');
while(pp~=2 && pp~=4 && pp~=6 && pp~=8)
pp=input('--> Ban da nhap sai so doi cuc. Vui long nhap lai so doi cuc :
2p= ');
end
hs=input('- Nhap gia tri hieu suat (%): hs= ');
cosfi=input('- Nhap gia tri cosfi : cosfi= ');
ii=input('- Chon cap cach dien: Neu tinh theo cap cach dien B thi nhan PHIM
1. Neu tinh theo cap cach dien F thi nhan PHIM 2 :');
if ii==1
disp(' -> Ban da chon cach tinh theo cap cach dien B');
elseif ii==2
disp(' -> Ban da chon cach tinh theo cap cach dien F');
end
disp(' ');
% Chuong 1: TINH TOAN CAC KICH THUOC CHU YEU CUA DONG CO:
disp(' Chuong 1: TINH TOAN CAC KICH THUOC CHU YEU CUA DONG CO:');
disp(' ------------*****------------');

% 1/ TOC DO DONG BO
ndb=60*50*2/pp;
x=['1/ Toc do dong bo : ndb= ',num2str(ndb),' (vong/phut)'];disp(x);
% 2/ DUONG KINH NGOAI STATOR: (bang 10.1 trang 601)
t=pp/2+1;
if ii==1
load('capB.mat')
for i=1:25
if(capB(i,t)==P)
h=capB(i,1);
end
end
end
if ii==2
load('capF.mat');
for i=1:25

if(capF(i,t)==P)
h=capF(i,1);
end
end
end
x=['2/ Chieu cao tam truc h= ',num2str(h),'(mm) tra theo bang
IV.2/TR.601/602'];disp(x);
load('traDn.mat');
hh=h/10;
for i=1:20
if traDn(i,1)==hh
Dn=traDn(i,2);
end
end
x=['- Theo bang 10.3/230 ta tra duoc duong kinh ngoai tieu chuan:
Dn=',num2str(Dn),' (cm)'];disp(x);
%3/ TRA DUONG KINH TRONG
tk=pp/2;
tt=pp/2+5;
load('traKD.mat');
kDmin=traKD(1,tk);
kDmax=traKD(1,tt);
x=['3/ Theo bang 10.2/230 ta co: kd= ',num2str(kDmin),'-
>',num2str(kDmax),];disp(x);
D1=kDmin*Dn;
D2=kDmax*Dn;
x=['- Gia tri duong kinh trong trong khoang: D = (',num2str(D1),'-
>',num2str(D2),') (cm)'];disp(x);
D=input('--> Nhap gia tri duong kinh trong ma ban muon chon: D= ');
while(D<D1 || D>D2)
D=input('--> Ban da nhap sai. Vui long nhap lai gia tri: D= ');
end
%4/ CONG SUAT TINH TOAN
disp('4/ Cong suat tinh toan: ');
load('traKE.mat');
for i=1:10
if traKE(i,1)<=Dn && traKE(i+1,1)>=Dn
aa=traKE(i,1);
bb=traKE(i+1,1);
cc=traKE(i,t);
dd=traKE(i+1,t);
break;
end
end
kE=(Dn-aa)*(dd-cc)/(bb-aa)+cc;
Ptt=kE*P/(hs*cosfi);
x=['- Gia tri cong suat tinh toan: Ptt= ',num2str(Ptt),'(KW)-> trong do kE
=',num2str(kE),' theo hinh 10.2/tr.231'];disp(x);
%5/ CHIEU DAI TINH TOAN CUA LOI SAT STATO
if pp/2==1
kd=0.96;
elseif pp/2>1
kd=0.92;
end
as=0.64;
ks=1.11;
tt1=pp/2+5;
tt2=tt1+4;
tt3=tt2+4;
if (h<=132)
load('tra_ABmin.mat');
%load('traAB1.mat');
for i=1:9
if Dn>=tra_ABmin(i,1) && Dn<=tra_ABmin(i+1,1)
d1tra=tra_ABmin(i,1);
d2tra=tra_ABmin(i+1,1);

a1min=tra_ABmin(i,t);
a2min=tra_ABmin(i+1,t);
a1max=tra_ABmin(i,tt2);
a2max=tra_ABmin(i+1,tt2);

b1min=tra_ABmin(i,tt1);
b2min=tra_ABmin(i+1,tt1);
b1max=tra_ABmin(i,tt3);
b2max=tra_ABmin(i+1,tt3);
end
end
elseif (h>=160 && h<=250)
%load('traAB2.mat');
load('tra_ABmax.mat');
for i=1:11
if Dn>=tra_ABmax(i,1) && Dn<=tra_ABmax(i+1,1)
d1tra=tra_ABmax(i,1);
d2tra=tra_ABmax(i+1,1);

a1min=tra_ABmax(i,t);
a2min=tra_ABmax(i+1,t);
a1max=tra_ABmax(i,tt2);
a2max=tra_ABmax(i+1,tt2);

b1min=tra_ABmax(i,tt1);
b2min=tra_ABmax(i+1,tt1);
b1max=tra_ABmax(i,tt3);
b2max=tra_ABmax(i+1,tt3);
end
end
end
Amin=(Dn-d1tra)*(a2min-a1min)/(d2tra-d1tra)+a1min;
Amax=(Dn-d1tra)*(a2max-a1max)/(d2tra-d1tra)+a1max;
Bmin=(Dn-d1tra)*(b2min-b1min)/(d2tra-d1tra)+b1min;
Bmax=(Dn-d1tra)*(b2max-b1max)/(d2tra-d1tra)+b1max;
x=['- Gia tri A: A=(Amin-->Amax)= (',num2str(Amin),'-->',num2str(Amax),')
'];disp(x); % NOTE: trang 232/
A=input('--> Nhap gia tri A ma ban muon chon: A= ');
while(A<Amin || A>Amax)
A=input('--> Ban da nhap sai. Vui long nhap lai gia tri: A= ');
end
x=['- Gia tri Bss: Bss=(Bssmin-->Bssmax)= (',num2str(Bmin),'--
>',num2str(Bmax),') '];disp(x);
Bss=input('--> Nhap gia tri Bss ma ban muon chon: Bss= ');
while(Bss<Bmin || Bss>Bmax)
Bss=input('--> Ban da nhap sai. Vui long nhap lai gia tri: Bss= ');
end
%elseif (h>=280)
% load('traAB3.mat');
% for i=1:9
% if Dn>=traAB3(i,1) && Dn<=traAB3(i+1,1)
% aaa=traAB3(i,1);
% bbb=traAB3(i+1,1);
% ccc=traAB3(i,t);
% ddd=traAB3(i+1,t);
% eee=traAB3(i,tt);
% fff=traAB3(i+1,tt);
% end
%end
%end
%A=(Dn-aaa)*(ddd-ccc)/(bbb-aaa)+ccc;
%Bss=(Dn-aaa)*(fff-eee)/(bbb-aaa)+eee;
Ls=round((6.1*10^7*Ptt)/(as*ks*kd*A*Bss*D*D*ndb));
%%Ls=Ls/10;
x=['5/ Chieu dai tinh toan cua loi sat stato: Ls=
',num2str(Ls),'(cm)'];disp(x);
x=['-> Do loi sat ngan nen lam thanh mot khoi -> chieu dai loi sat stato,
roto bang: L1=L2=Ls= ',num2str(Ls),'(cm)'];disp(x);
x=['- Theo hinh 10.3a/tr.233 lay: A= ',num2str(A),'(A/cm), Bs=
',num2str(Bss),'(T)'];disp(x);
%6/ BUOC CUC
to=pi*D/pp;
x=['6/ Buocc cuc: to= ',num2str(to),'(cm)'];disp(x);
%7/ LAP PHUONG AN SO SANH
londa=Ls/(to/10);
x=['7/ Gia tri londa: londa= ',num2str(londa),];disp(x);
%8/ DONG DIEN DINH DIEN MUC.
I1=P*10^3*100/(3*220*hs*cosfi);
I1=round(I1);
I1=I1/100;
x=['8/ Dong dien pha dinh muc: I1= ',num2str(I1),'(A)'];disp(x);
disp(' ');

% CHUONG 2: DAY QUAN, RANH STATO VA KHE HO KHONG KHI


disp(' Chuong 2: DAY QUAN, RANH STATO VA KHE HO KHONG KHI');
disp(' ----------*****----------');
%9/ SO RANH STATOR.
%if P<1
% q1=2;
%elseif (P>1 && P<=22 && ndb==1000) ||(P>1 && P<=22 && ndb==1500 )||(P>1 &&
P<=22 && ndb==750)
% q1=3;
%elseif P>1 && P<7 && ndb==3000
% q1=4;
%elseif P>22 && ndb>=1000
% q1=4;
%elseif P>=7 && ndb==3000
% q1=6;
%end
if P<1
if ndb>=750
q1=2;
end
elseif P>1 && P<7
if ndb<3000
q1=3;
elseif ndb==3000
q1=4;
end
elseif P>7 && P<=22
if ndb<3000
q1=3;
elseif ndb==3000
q1=6;
end
elseif P>22
if ndb<3000
q1=4;
elseif ndb==3000
q1=6;
end
end
x=['- Ta thay cong suat dong co P= ',num2str(P),'(kW) va toc do dong co ndb=
',num2str(ndb),'(vong/phut)--> nen chon q1= ',num2str(q1),];disp(x);
Z1=6*pp*q1/2;
x=['9/ So ranh stator: Z1= ',num2str(Z1),'(ranh)'];disp(x);
%10/ BUOC RANH STATOR
t1= pi*D/Z1;
x=['10/ Buoc ranh stator: t1= ',num2str(t1),'(cm)'];disp(x);
%11/ SO THANH DAN TAC DUNG CUA 1 RANH:
if P<1
a1=1;
n1=1;
elseif P>=1 && P<7
a1=1;
n1=1;
elseif P>=7 && P<= 22
a1=2;
n1=1;
elseif P>22;
a1=4;
n1=2;
end
ur1=A*t1*a1/I1;
ur1=round(ur1);% fix
x=['11/ So thanh dan tac dung cua mot ranh: ur1=
',num2str(ur1),'(vong/phut)'];disp(x);
x=['-> Chon so bo so mach nhanh song song: a1= ',num2str(a1),];disp(x);
x=['-> Chon so bo so so chap n1= ',num2str(n1),];disp(x);
%12/ SO VONG DAY NOI TIEP CUA MOT PHA.
w1=pp*q1*ur1/(2*a1);
w1=round(w1);
x=['12/ So vong day noi tiep cua mot pha: w1=
',num2str(w1),'(vong)'];disp(x);
%13/ TIET DIEN VA DUONG KINH DAY DAN.
load('traAJ1.mat');
load('traAJ2.mat');
load('traAJ3.mat');
if h<=132
for i=1:9
if traAJ1(i,1)<=Dn && traAJ1(i+1,1)>=Dn
a4=traAJ1(i,1);
a5=traAJ1(i+1,1);
b4=traAJ1(i,t);
b5=traAJ1(i+1,t);
end
end
elseif h>=160 && h<= 250
for i=1:9
if traAJ2(i,1)<=Dn && traAJ2(i+1,1)>=Dn
a4=traAJ2(i,1);
a5=traAJ2(i+1,1);
b4=traAJ2(i,t);
b5=traAJ2(i+1,t);
end
end
elseif h>=280 && h<=355
for i=1:9
if traAJ3(i,1)<=Dn && traAJ3(i+1,1)>=Dn
a4=traAJ3(i,1);
a5=traAJ3(i+1,1);
b4=traAJ3(i,t);
b5=traAJ3(i+1,t);
end
end
end
AJ=(Dn-a4)*(b5-b4)/(a5-a4)+b4;
AJ=1800;
Jphay1=AJ/A;
sphay1=(I1/(a1*n1*Jphay1));
load('tratietdien.mat');
for i=1:93
if tratietdien(i,1)<=sphay1 && tratietdien(i+1,1)>=sphay1
smin=tratietdien(i,1);
smax=tratietdien(i+1,1);
dcdmin=tratietdien(i,2);
dcdmax=tratietdien(i+1,2);
dmin=tratietdien(i,3);
dmax=tratietdien(i+1,3);
break;
end
end
x=['13/ Theo hinh 10-4/tr.237 chon tich so: AJ=
',num2str(AJ),'(A2/cm.mm2)'];disp(x);
x=['- Mat do dong dien: Jphay1= ',num2str(Jphay1),'(A/mm2)',];disp(x);
x=['- Tuong ung voi: s1phay= ',num2str(sphay1),'-> Tiet dien day so bo tra
bang VI.1/TR.618 thi s nam trong khoang: (smin-smax)=(',num2str(smin),'-
>',num2str(smax),')(mm2)'];disp(x);
ss=input('--> Neu chon s theo tiet dien smin nhan PHIM 1. Neu chon s theo
tiet dien smax thi nhan PHIM 2: ');
if ss==1
s1=smin;
d=dmin;
dcd=dcdmin;
x=['--> Ban da chon theo thiet dien min voi: s=
',num2str(s1),'(mm2)'];disp(x);
x=['- Theo phuc luc VI, Bang VI.1 chon day dong trang men PETV co duong kinh:
(d/dcd)= (',num2str(dmin),'->',num2str(dcdmin),')'];disp(x);
elseif ss==2
s1=smax;
d=dmax;
dcd=dcdmax;
x=['--> Ban da chon theo thiet dien max voi: s=
',num2str(s1),'(mm2)'];disp(x);
x=['- Theo phuc luc VI, Bang VI.1 chon day dong trang men PETV co duong kinh:
(d/dcd)= (',num2str(dmax),'->',num2str(dcdmax),')'];disp(x);
end
%14/ KIEU DAY QUAN
disp('14/ Chon kieu day quan:');
%yy=input('-> Neu quan dong co theo buoc thieu thi nhan PHIM 1. Neu quan dong
co theo buoc du thi nhan PHIM 2: ');
t=Z1/pp;
if P<7
x=['- Dong co P= ',num2str(P),' < 7(kW)--->nen ta theo kieu day quan 1
lop buoc du:'];disp(x);
giam=0;
y=Z1/pp;
beta=y/t;
x=['- Buoc cuc tinh: t= ',num2str(t),];disp(x);
x=['- Buoc day quan: y= ',num2str(y),];disp(x);
x=['- Voi he so beta: beta= ',num2str(beta),];disp(x);
elseif P>=7 && P<22
x=['- Dong co P= ',num2str(P),' > 7(kW)--->nen ta theo kieu day quan 2
lop buoc thieu:'];disp(x);
giam=1;
y=Z1/pp-1;
beta=y/t;
x=['- Buoc cuc tinh: t= ',num2str(t),];disp(x);
x=['- Buoc day quan: y= ',num2str(y),];disp(x);
x=['- Voi he so beta: beta= ',num2str(beta),];disp(x);
else
x=['- Dong co P= ',num2str(P),' > 7(kW)--->nen ta theo kieu day quan 2
lop buoc thieu:'];disp(x);
giam=2;
y=Z1/pp-2;
beta=y/t;
x=['- Buoc cuc tinh: t= ',num2str(t),];disp(x);
x=['- Buoc day quan: y= ',num2str(y),];disp(x);
x=['- Voi he so beta: beta= ',num2str(beta),];disp(x);
end
%if(yy==1)
%y=Z1/pp-1;
%beta=y/t;
%x=['- Ban da chon quan theo buoc thieu voi: y= ',num2str(y),];disp(x);
%x=['- Voi he so beta: beta= ',num2str(beta),];disp(x);
%elseif yy==2
%y=Z1/pp;
%beta=y/t;
%x=['- Ban da chon quan theo buoc du voi: y= ',num2str(y),];disp(x);
%x=['- Voi he so beta: beta= ',num2str(beta),];disp(x);
%end
%15/ HE SO DAY QUAN
disp('15/ He so day quan:');
ky=sin((beta*pi)/2);
x=['- He so buoc day : ky= ',num2str(ky),];disp(x);
alpha=pp*360/(2*Z1);
kr=sind(q1*alpha/2)/(q1*sind(alpha/2));
kD=ky*kr;
x=['- He so day quan: kD= ',num2str(kD),];disp(x);
x=['-> Trong do: alpha= ',num2str(alpha),'(do)'];disp(x);
%16/ TU THONG KHE HO KHONG KHI
fi=kE*220/(4*ks*kD*50*w1);
x=['16/ Tu thong khe ho khong khi: fi= ',num2str(fi),'(Wb)'];disp(x);
%17/ MAT DO TU THONG KHE HO KHONG KHI
Bs=fi*10^4/(as*to*Ls);
x=['17/ Mat do tu thong kh ho khong khi: Bs= ',num2str(Bs),'(T)'];disp(x);
%18/ SO BO DINH CHIEU RONG RANG VA CHIEU CAO GONG
load('trabz.mat');
load('trabg.mat');
ti=pp/2;
ttt=pp/2+4;
if h>=50 && h<=132
bzmin=trabz(1,ti);
bgmin=trabg(1,ti);
bzmax=trabz(1,ttt);
bgmax=trabg(1,ttt);
elseif h>=160 && h<=250
bzmin=trabz(2,ti);
bgmin=trabg(2,ti);
bzmax=trabz(2,ttt);
bgmax=trabg(2,ttt);
end
disp('18-19/ So bo chieu rong cua rang va chieu cao cua gong:');
x=['- Mat do tu thong Bz tren rang nam trong khoang:Bz= (',num2str(bzmin),'-
>',num2str(bzmax),')'];disp(x);
x=['- Mat do tu thong trong gong stato nam trong khoang:Bg=
(',num2str(bgmin),'->',num2str(bgmax),')'];disp(x);
h41=0.5;
kc=0.95;
isok=false;
for Bz1=bzmax:-0.01:bzmin
if isok==true
break;
end
for BG1=bgmin:0.01:bgmax
hphayg1=fi*10000/(2*BG1*Ls*kc);
bphayz1=Bs*t1/(Bz1*kc);
hr1=0.5*(Dn*10-D*10)-hphayg1*10;
hr1=round(hr1);
d1s=(pi*(D*10+2*h41)-bphayz1*10*Z1)/(Z1-pi);
d1s=round(d1s);
d2s=(pi*(Dn*10-2*hphayg1*10)-bphayz1*10*Z1)/(Z1+pi);
d2s=round(d2s);
h12=hr1-d2s/2-h41;
h12=round(h12);
Sphayr=pi*(d1s*d1s+d2s*d2s)/8+(d1s+d2s)*(h12-d1s/2)/2;
if h>=90 && h<=132
c=0.25;
cphay=0.35;
Scd=(pi*d2s/2+2*h12)*c+pi*d1s*cphay/2;
elseif h>=50 && h<=80
c=0.2;
cphay=0.3;
Scd=(pi*d2s/2+2*h12)*c+pi*d1s*cphay/2;
elseif h>=160 && h<=250
c=0.4;
cphay=0.5;
Scd=(pi*d2s/2+2*h12+(d1s+d2s))*c+pi*d1s*cphay/2;
end
Sr=Sphayr-Scd;
kld=ur1*n1*dcd*dcd/Sr;
if kld>=0.7 && kld<=0.75% && abs(bz1-bphayz1)<0.02
isok = true;
break;
end
end
end
while (kld < 0.7 || kld > 0.75)
disp('He so lap day khong thoa man');
disp(kld);
return;
end;
x=['- So bo chieu rong cua rang: bphayz1= ',num2str(bphayz1),' (cm) -> voi he
so Bz1 tra bang 10.5b/tr241 Bz1= ',num2str(Bz1),'(T)'];disp(x);
x=['- So bo chieu cao cua gong: hphayg1= ',num2str(hphayg1),' (cm) -> voi he
so Bg1 tra bang 10.5a/tr240 Bg1= ',num2str(BG1),'(T)'];disp(x);
disp('20/ Kich thuc ranh va cach dien:');
b41=(dcd+1.5)*10;
b41=round(b41);
b41=b41/10;

x=['- hr1= ',num2str(hr1),'(mm)'];disp(x);


x=['- h12= ',num2str(h12),'(mm)'];disp(x);
x=['- d1s= ',num2str(d1s),'(mm)'];disp(x);
x=['- d2s= ',num2str(d2s),'(mm)'];disp(x);
x=['- b41= ',num2str(b41),'(mm)'];disp(x);
x=['- h41= ',num2str(h41),'(mm)'];disp(x);
x=['- Theo bang VIII-1 o phu luc VIII chieu day cach dien cua ranh la: c=
',num2str(c),'(mm);cua nem la cphay = ',num2str(cphay),'(mm)'];disp(x);
x=['- Dien tich ranh tru nem: Sphayr= ',num2str(Sphayr),'(mm2)'];disp(x);
x=['- Dien tich cach dien ranh: Scd= ',num2str(Scd),'(mm2)'];disp(x);
x=['- Dien tich co ich cua ranh: Sr= ',num2str(Sr),'(mm2)'];disp(x);
x=['- - - - - > He so lap day ranh: kld=
',num2str(kld),];disp(x);
%21/ BE RONG RANG SATTO
bz11=(pi*(D+2*h41*0.1+d1s*0.1))/Z1-(d1s*0.1);
bz12=(pi*(D+2*(h41*0.1+h12*0.1)))/Z1-(d2s*0.1);
bz1=(bz11+bz12)/2;
hg1=(Dn-D)/2-(hr1/10)+(d2s/60);
disp('21/ Be rong rang stator:');
x=['- b1phayz1= ',num2str(bz11),'(cm)'];disp(x);
x=['- b2phayz1= ',num2str(bz12),'(cm)'];disp(x);
x=['--> bz1= (b1phayz1+b2phayz1)/2 = ',num2str(bz1),'(cm)'];disp(x);
%22/ CHIEU CAO GONG STATOR.
x=['22/ Chieu cao gong stato: hg1= ',num2str(hg1),'(cm)'];disp(x);
%23/ KHE HO KHONG KHI
if(P<=20)
if(pp==2)
deltal=0.3+1.5*D*10/1000;
elseif(pp>=4)
deltal=0.25+D*10/1000;
end
elseif (P>20)
deltal=D*10*(1+9/pp)/1200;
end
x=['23/ Khe ho khong khi: delta=', num2str(deltal),'(mm)'];disp(x);
load('tradeta.mat');
j=pp/2+1;
for i=1:14
if(tradeta(i,1)==h)
delta=tradeta(i,j);
break;
end
end
Denta=delta/10;
x=[' Theo nhung may da che tao bang 10.8, ta lay: delta=',
num2str(delta),'(mm) =',num2str(Denta),'(cm)'];disp(x);
disp(' ');

% CHUONG 3: DAY QUAN, RANH VA GONG ROTO


disp(' Chuong 3: DAY QUAN, RANH VA GONG ROTO');
disp(' ----------*****---------');
%24/ SO RANH ROTO CHON THEO BANG 10.6
load('traZ2.mat');
for i=1:13
if traZ2(i,1)==Z1
Z2=traZ2(i,j);
end
end
%%Z2=38;
x=['24/ So ranh roto chon theo bang 10.6/tr246: Z2= ',num2str(Z2),'
(ranh)'];disp(x);
%25/ DUONG KINH NGOAI ROTO
Dphay=D-2*Denta;
x=['25/ Duong kinh ngoai roto: Dphay= ',num2str(Dphay),'(cm)'];disp(x);
%26/ BUOC RANG ROTO
t2=pi*Dphay/Z2;

x=['26/ Buoc rang roto: t2=',num2str(t2),'(cm)'];disp(x);


%27/ SO BO DINH CHIEU RONG RANG
BZ2=Bz1;
bphayZ2=Bs*t2/(BZ2*kc);
x=['27/ So bo dinh chieu rong rang roto: bphayz2= ',num2str(bphayZ2),'(cm);
voi BZ2= ',num2str(BZ2),'(T)'];disp(x);
%28 DUONG KINH TRUC ROTO
Dt=0.3*D;
x=['28/ Duong kinh truc roto: Dt= ',num2str(Dt),'(cm)'];disp(x);
%29/ DONG DIEN TRONG THANH DAN ROTO
load('trakI.mat');
for i=1:11
if trakI(i,1)==cosfi
kI=trakI(i,2);
break;
end
if trakI(i,1)<=cosfi && trakI(i+1,1)>=cosfi
a8=trakI(i,1);
a9=trakI(i+1,1);
b8=trakI(i,2);
b9=trakI(i+1,2);
end
end
kI=(cosfi-a8)*(b9-b8)/(a9-a8)+b8;
Itd=kI*I1*6*w1*kD/Z2;
x=['29/ Dong dien trong thanh dan roto: Itd= ',num2str(Itd),' trong do: kI=
',num2str(kI),' lay theo hinh 10-5/tr.244'];disp(x);
%30/DONG DIEN TRONG VANH NGAN MACH
Iv=Itd/(2*sin(pi*pp/(2*Z2)));
x=['30/ Dong dien trong vanh ngan mach: Iv= ',num2str(Iv),'(A)'];disp(x);
% 31,32,33,34,35,36:
Jv=2.5;
BZ2=Bz1;
BG2=BG1;
isok=false;
if h<=160
for Jtd=2.5:0.1:3.5
if isok==true
break;
end
% for h42=0.5:0.1:1
h42=0.5;
hphayg2=fi*10000/(2*BG1*Ls*kc);
bphayz2=Bs*t2/(BZ2*kc);
hr2=0.5*(Dphay*10-Dt*10)-hphayg2*10;
dr1=round((pi*(Dphay*10-2*h42)-bphayz2*10*Z2)/(Z2+pi));
dr2=round((pi*(Dphay*10-2*hr2)-bphayz2*10*Z2)/(Z2-pi));
Sphaytd=Itd/Jtd;
ka=Sphaytd;
dr=dr1;
b42=1;
Sv=Iv/Jv;% Tiet dien vanh ngan mach
hr12=hr2-dr1/2-dr2/2-h42;
A=1.2*hr2;
%A=ceil(A);
B=Sv/A;
%B=ceil(B);
Dv=D*10-(A+1);
Snm=A*B;% dien tich vanh ngan mach
Sr2=pi*(dr1*dr1+dr2*dr2)/8+0.5*(dr1+dr2)*(hr12);

bZ21=(pi*(Dphay-2*h42*0.1-dr1*0.1)/Z2-dr1*0.1); %
tr.102 sgk
bZ22=(pi*(Dphay-2*0.1*hr2+0.1*dr2)/Z2-dr2*0.1);

bZ213=(bZ21+bZ22)/2;

hg2=((Dphay-Dt)/2)-(hr2/10)+(dr2/60);
if abs(bZ213-bphayz2)<=0.08 && abs(Sr2-Sphaytd)<= 2
isok=true;
break;
end

end

end
if h>=180
isok=false;
for Jtd=2.5:0.1:3.5
if isok==true
break;
end
for hr2=25:1:45
if isok==true
break;
end
for dr=3.5:0.1:7
if isok==true
break;
end
for h42=0.5:0.1:1.5
b42=1.5;
dr1=dr;
dr2=dr;
Sphaytd=Itd/Jtd;
ka=Sphaytd;
Sv=Iv/Jv;% Tiet dien vanh ngan mach
hr12=hr2-dr-h42;
A=1.2*hr2;
A=ceil(A);
B=Sv/A;
B=ceil(B);
Dv=D*10-(A+1);
Snm=A*B;% dien tich vanh ngan mach
Sr2=pi*dr*dr/4+hr12*dr;
bZ213=pi*(Dphay-2*h42*0.1-(4/3)*(hr12+dr)*0.1)/Z2-
(dr*0.1);
hg2=((Dphay-Dt)/2)-(hr2/10)+(dr1/60);
if abs(Sr2-Sphaytd)<=1
isok=true;
break;
end
end
end
end
end
end
abs(Sr2-Sphaytd)
bphayz2=Bs*t2/(BZ2*kc);
while abs(Sr2-Sphaytd)>2
disp(' --> Khong tinh duoc Dien tich ranh (Sr2) roto va Tiet dien thanh
nhom (Sphaytd) gan bang nhau.');
return;
end
%31/ TIET DIEN THANH DAN BANG NHOM
%%bZ21= bZ21
%%Sphaytd= Itd/3;
x=['- So bo chieu rong cua rang: bphayz2= ',num2str(bphayz2),' (cm) -> voi he
so Bz1 tra bang 10.5b/tr241 Bz2= ',num2str(BZ2),'(T)'];disp(x);
x=['31/ Tiet dien thanh nhom : Sphaytd= ',num2str(Sphaytd),'(mm2) -> Trong do
chon J2= ',num2str(Jtd),'(A/mm2)'];disp(x);
%32/
x=['32/ So bo chon mat do dong dien trong vanh ngan mach: Jv=
',num2str(Jv),'(A/mm2)'];disp(x);
x=['- Tiet dien vanh ngan mach: Sv= ',num2str(Sv),' (mm2)'];disp(x);
%33/ KICH THUOC RANH ROTO VA VANH NGAN MAXH
disp('33/ Kich thuoc ranh roto va vanh ngan mach:');
if h<=160
disp('---> Chon ranh roto theo hinh QUA LE');
elseif h>=180
disp('---> Chon ranh roto theo hinh OVAN');
end
x=['- hr2= ',num2str(hr2),'(mm)'];disp(x);
x=['- hr12= ',num2str(hr12),'(mm)'];disp(x);
x=['- dr1= ',num2str(dr1),'(mm)'];disp(x);
x=['- dr2= ',num2str(dr2),'(mm)'];disp(x);
x=['- b42= ',num2str(b42),'(mm)'];disp(x);
x=['- h42= ',num2str(h42),'(mm)'];disp(x);
x=['- Dv=',num2str(Dv),'(mm) Voi:'];disp(x);
x=['+ A =',num2str(A),'(mm)'];disp(x);
x=['+ B = ',num2str(B),'(mm)'];disp(x);
%34/ DIEN TICH RANH ROTO
x=['34/ Dien tich ranh roto: Sr2= ',num2str(Sr2),' (mm2)'];disp(x);
%35/ DIEN TIC VANH NGAN MACH
x=['35/ Dien tich vanh ngan mach: Snm= ',num2str(Snm),' (mm2)'];disp(x);
%36/ BE RONG RANG ROTO O 1/3 CHIEU CAO RANG
x=['36/ Be rong rang roto o 1/3 chieu cao gong: bZ213=
',num2str(bZ213),'(cm)'];disp(x);
%37/ CHIEU CAO GONG ROTO
x=['37/ Chieu cao gong roto: hg2=',num2str(hg2),' (cm)'];disp(x);
%38/ LAM NGHIENG RANH O ROTO
disp('38/ Lam nghieng ranh o roto:');
bn=t1;
x=['- Do nghiheng bang mot buoc ranh stator: bn=t1=
',num2str(bn),'(cm)'];disp(x);
disp(' ');

% CHUONG 4: TINH TOAN MACH TU


disp(' Chuong 4: TINH TOAN MACH TU.');
disp(' --------*****--------' );
%39/ HE SO KHE HO KHONG KHI
v1=(b41/delta)^2/(5+b41/delta);
ks1=t1/(t1-v1*Denta);
v2=(b42/delta)^2/(5+b42/delta);
ks2=t2/(t2-v2*Denta);
kS=ks1*ks2;
x=['39/ He so khe ho khong khi: ks= ',num2str(kS),'voi:'];disp(x);
x=['- ks1= ',num2str(ks1),' trong do: v1= ',num2str(v1),];disp(x);
x=['- ks2= ',num2str(ks2),' trong do: v2= ',num2str(v2),];disp(x);
%40/
disp('40/ Dung thep ky thuat dien can nguoi loai 2212.');
%41/ SUC TU DONG KHE HO KHONG KHI
Fs=fix(1.6*Bs*kS*Denta*10^4);
x=['41/ Suc tu dong khe ho khong khi: Fs= ',num2str(Fs),'(A)'];disp(x);
%42/ MAT DO TU THONG O GONG STATOR.
Bz1=Bs*t1*100/(bz1*kc);
Bz1=round(Bz1);
Bz1=Bz1/100;
x=['42/ Mat do tu thong o gong stator: Bz1= ',num2str(Bz1),' (T)'];disp(x);
%43/
load('bangv_6.mat');
for i=1:200
if bangv_6(i,1)==Bz1
Hz1=bangv_6(i,2);
end
end
x=['43/ Cuong do tu truong tren rang stato, theo bang V-6/tr.608 o phu luc V
co: Hz1= ',num2str(Hz1),'(A/cm)'];disp(x);
%44/ SUC TU DONG TREN RANG STATO.
hphayz1=hr1*0.1-d2s*0.1/3;
Fz1=fix(2*hphayz1*Hz1);
x=['44/ Suc tu dong tren rang stato: Fz1= ',num2str(Fz1),' (A) Trong do:
hphayz1= ',num2str(hphayz1),'(cm)'];disp(x);
%45/ MAT DO TU THONG O RANG ROTO.
t2=fix(t2*100);
t2=t2/100;
Bz2=Bs*t2*100/(bZ213*kc);
Bz2=fix(Bz2);
Bz2=Bz2/100;
x=['44/ Mat do tu thong o rang roto: Bz2= ',num2str(Bz2),' (T)'];disp(x);
%46/
for i=1:200
if bangv_6(i,1)==Bz2
Hz2=bangv_6(i,2);
end
end
x=['46/ Cuong do tu truong tren rang roto, theo bang V-6/tr.608 o phu luc V
co: Hz2= ',num2str(Hz2),'(A/cm)'];disp(x);
%47/ SUC TU DONG TREN RANG ROTO
hphayz2=hr2*0.1-dr2*0.1/3;
Fz2=fix(2*hphayz2*Hz2);
x=['47/ Suc tu dong tren rang roto: Fz2= ',num2str(Fz2),' (A) Trong do:
hphayz2= ',num2str(hphayz2),'(cm)'];disp(x);
%48/ HE SO BAO HOA RANG.
kz=(Fs+Fz1+Fz2)/Fs;
x=['48/ He so bao hoa rang: kz= ',num2str(kz),];disp(x);
%49/ MAT DO TU THONG TREN GONG STATO.
Bg1=fi*10^4*100/(2*hg1*Ls*kc);
Bg1=round(Bg1);
Bg1=Bg1/100;
x=['49/ Mat do tu thong tren gong stato: Bg1= ',num2str(Bg1),' (T)'];disp(x);
%50/
load('bangv_9.mat');
for i=1:170
if bangv_9(i,1)==Bg1
Hg1=bangv_9(i,2);
end
end
x=['50/ Cuong do tu truong o gong stato, theo bang V-9/tr.611 o phu luc V co:
Hg1= ',num2str(Hg1),'(A/cm)'];disp(x);
%51/ CHIEU DAI MACH TU O GONG STATO.
Lg1=pi*(Dn-hg1)/pp;
x=['51/ Chieu dai mach tu o gong stato: Lg1= ',num2str(Lg1),' (cm)'
];disp(x);
%52/ SUC TU DONG O GONG STATO.
Fg1=fix(Lg1*Hg1);
x=['52/ Suc tu dong o gong stato: Fg1= ',num2str(Fg1),' (A)'];disp(x);
%53/ MAT DO TU THONG TREN GONG ROTO.
Bg2=fi*10^4*100/(2*hg2*Ls*kc);
Bg2=round(Bg2);
Bg2=Bg2/100;
x=['53/ Mat do tu thong tren gong roto: Bg2= ',num2str(Bg2),' (T)'];disp(x);
%54/
for i=1:170
if bangv_9(i,1)==Bg2
Hg2=bangv_9(i,2);
end
end
x=['50/ Cuong do tu truong o gong roto, theo bang V-9/tr.611 o phu luc V co:
Hg2= ',num2str(Hg2),'(A/cm)'];disp(x);
%55/ CHIEU DAI MACH TU O GONG ROTO.
Lg2=pi*(Dt+hg2)/pp;
x=['55/ Chieu dai mach tu o gong roto: Lg2= ',num2str(Lg2),' (cm)'];disp(x);
%56/ SUC TU DONG TREN GONG ROTO.
Fg2=fix(Lg2*Hg2);
x=['56/ Suc tu dong tren gong roto: Fg2= ',num2str(Fg2),' (A)'];disp(x);
%57/ TONG SUC TU DONG CUA MACH TU.
F=Fs+Fz1+Fz2+Fg1+Fg2;
x=['57/ Tong suc tu dong cua mach tu: F= ',num2str(F),' (A)'];disp(x);
%58/ HE SO BAO HOA TOAN MACH
ku=F/Fs;
x=['58/ He so bao hoa toan mach: ku= ',num2str(ku),];disp(x);
%59/ DONG DIEN TU HOA.
Iu=fix(10*pp*F/(2*2.7*w1*kD));
Iu=Iu/10;
Iuphantram=Iu*100/I1;
x=['59/ Dong dien tu hoa: Iu= ',num2str(Iu),' (A)'];disp(x);
x=['- Dong dien tu hoa phan tram: Iu%= ',num2str(Iuphantram),'(%)'];disp(x);
disp(' ');

% CHUONG 5: THAM SO DONG CO DIEN O CHE DO DINH MUC


disp(' Chuong 5: THAM SO DONG CO DIEN O CHE DO DINH MUC.');
disp(' -----------*****-----------');
%60/ CHIEU DAI PHAN DAU NOI CUA DAY QUAN STATO.
load('traKD1.mat');
th=pp/2;
Kd1=traKD1(1,th);
ty=pi*(D+hr1*0.1)*y/Z1;
ld1=Kd1*ty+2*1;
x=['60/ Chieu dai phan dau noi cua day quan stato: ld1= ',num2str(ld1),'
(cm)'];disp(x);
x=['- trong do: ty= ',num2str(ty),' ; Kd1= ',num2str(Kd1)];disp(x);
%61/ CHIEU DAI TRUNG BINH NUA VONG DAY CUA DAY QUAN STATO.
ltb= Ls+ld1;
x=['61/ Chieu dai trung binh nua vong day quan stato: ltb= ',num2str(ltb),'
(cm)'];disp(x);
%62/ CHIEU DAI DAY QUAN MOT PHA CUA STATO.
L1=2*ltb*w1*10^-2;
x=['62/ Chieu dai day quan 1 pha cua stato: L1= ',num2str(L1),'
(m)'];disp(x);
%63/ DIEN TRO TAC DUNG CUA DAY QUAN STATO.
r1=L1/(46*n1*a1*s1);
r1sao=r1*I1/220;
x=['63/ Dien tro tac dung cua day quan stato: r1= ',num2str(r1),'
(0m)'];disp(x);
x=['- Tinh theo don vi tuong doi: r1sao= ',num2str(r1sao),];disp(x);
%64/ DIEN TRO TAC DUNG CUA DAY QUAN ROTO
rtd=Ls*10^-2/(23*Sr2);
x=['64/ Dien tro tac dung cua day quan roto: rtd=
',num2str(rtd),'(0m)'];disp(x);
%65/ DIEN TRO VANH NGAN MACH.
rv=pi*Dv*0.1*10^-2/(23*Z2*Snm);
x=['65/ Dien tro vanh ngan mach: rv= ',num2str(rv),' (0m)'];disp(x);
%66/ DIEN TRO ROTO.
DEnta=2*sin(pi*pp/(2*Z2));
r2=rtd+2*rv/(DEnta*DEnta);
x=['66/ Dien tro roto: r2= ',num2str(r2),' (0m)'];disp(x);
x=['- Trong do: DEnta= ',num2str(DEnta),];disp(x);
%67/ HE SO QUY DOI.
gama=4*3*w1*kD*w1*kD/Z2;
x=['67/ He so quy doi: gama= ',num2str(gama),];disp(x);
%68/
r2phay=gama*r2;
r2phaysao=r2phay*I1/220;
x=['68/ Dien tro roto da quy doi: r2phay= ',num2str(r2phay),'(0m)'];disp(x);
x=['- Tinh theo don vi tuong doi: r2phaysao= ',num2str(r2phaysao),];disp(x);
%69/ HE SO TU DAN TAN RANH ROTO.
kpbeta=(1+3*beta)/4;% tr.125
kbeta=1/4+3/4*kpbeta;
h1=hr1-0.1*d2s-2*c-cphay;
h2=-(d1s/2-2*c-cphay);
londar1=h1*kbeta/(3*d1s)+(0.785-b41/(2*d1s)+h2/d1s+h41/b41)*kpbeta;
x=['69/ He so tu dan tan ranh roto: londar1=',num2str(londar1),];disp(x);
x=['- Trong do: beta= ',num2str(beta),' ; kbeta= ',num2str(kbeta),'; kpbeta=
',num2str(kpbeta),' theo cong thuc tr. 125'];disp(x);
x=['- h1= ',num2str(h1),'(mm); h2= ',num2str(h2),' (mm)'];disp(x);
%70/ HE SO TU TAN TAP STATO.
load('bang5_2a.mat');
load('bang5_3.mat');
q2=Z2/(3*pp);
td=Z2*2/pp;
q2dow=floor(q2);
q2up=ceil(q2);
for i=1:10
if bang5_3(i,1)<=td && bang5_3(i+1,1)>=td
tdmin=bang5_3(i,1);
tdmax=bang5_3(i+1,1);
rot11=bang5_3(i,q2dow);
rot12=bang5_3(i+1,q2dow);
rot21=bang5_3(i,q2up);
rot22=bang5_3(i+1,q2up);
end
end
rotmin=(td-tdmin)*(rot12-rot11)/(tdmax-tdmin)+rot11;
rotmax=(td-tdmin)*(rot22-rot21)/(tdmax-tdmin)+rot21;
rot=(q2-q2dow)*(rotmax-rotmin)/(q2up-q2dow)+rotmin;
x=['- TRi so cua he so rot theo bang 5.3/137: rot= ',num2str(rot),];disp(x);
x=['- Trong do: td= ',num2str(td),' ; q2= ',num2str(q2),' ->
q2dow=',num2str(q2dow),' & q2up= ',num2str(q2up),];disp(x);
x=['-> rotmin= ',num2str(rotmin),'; rotmax= ',num2str(rotmax),];disp(x);
q=Z1/(3*pp)+1;
for i=1:11
if bang5_2a(i,1)==giam
xitma1=bang5_2a(i,q);
end
end
xitma1=xitma1/100;
%x=['- tri so detat tra theo bang5.2a/tr.134: detat=
',num2str(detat),];disp(x);
k41=1-0.033*(b41*0.1*0.1*b41)/(t1*Denta);
londat1=0.9*t1*q1*kD*q1*kD*rot*k41*xitma1/(kS*Denta);
x=['70/ He so tu dan tan tap stato: londat1= ',num2str(londat1),];disp(x);
x=['- Trong do: k41= ',num2str(k41),];disp(x);
x=['- TRi so cua he so rot theo bang 5.3/137: rot= ',num2str(rot),];disp(x);
x=['- Voi: td= ',num2str(td),' ; q2= ',num2str(q2),' ->
q2dow=',num2str(q2dow),' & q2up= ',num2str(q2up),];disp(x);
x=['-> rotmin= ',num2str(rotmin),'; rotmax= ',num2str(rotmax),];disp(x);
x=['- Tri so detat tra theo bang5.2a/tr.134: xitma1=
',num2str(xitma1),];disp(x);
%71/ HE SO TU TAN PHAN DAU NOI.
londad1=0.34*q1*(ld1-0.64*beta*to)/Ls;
x=['71/ He so tu tan phan dau noi: londad1= ',num2str(londad1),];disp(x);
%72/ HE SO TU DAN TAN STATO.
Tonglonda1=londar1+londat1+londad1;
x=['72/ He so tu dan tan stato: Tonglonda1= ',num2str(Tonglonda1),];disp(x);
%73/ DIEN KHANG DAY QUAN STATO
x1=0.158*(50/100)*(w1/100)^2*(Ls*2/(pp*q1))*Tonglonda1;
x1sao=x1*I1/220;
x=['73/ Dien khang day quan stato: x1= ',num2str(x1),' (0m)'];disp(x);
x=['- Tinh theo don vi tuong doi: x1sao= ',num2str(x1sao),];disp(x);
%74/ HE SO TU DAN TAN RANH ROTO
h11=hr12+dr2/2-0.1*dr1;

londar2=(h11/(3*dr1))*((1-((pi*dr1*dr1)/(8*Sr2)))^2)+(0.66-
(b42/(2*dr1)))*1+(h42/b42);

x=['74/ He so tu dan tan ranh roto: londar2= ',num2str(londar2),'


(0m)'];disp(x);
%75/ HE SO TU DAN TANTAP ROTO.
load('bang5_2c.mat');
q2=Z2/(3*pp);
if q2<=4
for i=1:8
if bang5_2c(i,1)<=q2 && bang5_2c(i+1,1)>=q2
e7=bang5_2c(i,1);
e8=bang5_2c(i+1,1);
f7=bang5_2c(i,2);
f8=bang5_2c(i+1,2);
end
end
xitma2=((q2-e7)*(f8-f7)/(e8-e7)+f7)/100;
end
if q2>4
xitma2=9.15/(100*q2*q2);
end

londat2=0.9*t2*(q2*ks2)^2*xitma2/(kS*Denta);
ks2=ks2
x=['75/ He so tu dan tan tap roto: londat2= ',num2str(londat2),' Trong do:
xitma2= ',num2str(xitma2),' theo bang 5.2c/tr.136'];disp(x);
%76/ HE SO TU TAN PHAN DAU NOI.
londad2=(2.3*Dv*0.1/(Z2*Ls*DEnta*DEnta))*log10(4.7*Dv*0.1/(A*0.1+B*0.1*2));
x=['76/ He so tu tan tap roto: londad2=',num2str(londad2),];disp(x);
%77/ HE SO TU TAN DO RANH NGHIENG.
londarn=0.5*londat2*(bn/t2)^2;
x=['77/ He so tu tan do ranh nghieng: londarn= ',num2str(londarn),];disp(x);
%78/ HE SO TU TAN RANH ROTO.
Tonglonda2=londar2+londad2+londat2+londarn;
x=['78/ He so tu tan ranh roto: Tonglonda2= ',num2str(Tonglonda2),];disp(x);
%79/ DIEN KHANG TAN DAY QUAN ROTO.
x2=7.9*50*Ls*Tonglonda2*10^-8;
x=['79/ Dien khang tan day quan roto: x2= ',num2str(x2),' (0m)'];disp(x);
%80/ DIEN KHANG ROTO DA QUY DOI.
x2phay=gama*x2;
x2phaysao=x2phay*I1/220;
x=['80/ Dien khang roto da quy doi: x2phay= ',num2str(x2phay),'
(0m)'];disp(x);
x=['- Tinh theo he don vi tuong doi: x2phaysao=
',num2str(x2phaysao),];disp(x);
%81/ DIEN KHANG HO CAM.
x12=(220-Iu*x1)/Iu;
x12phay=x12*I1/220;
x=['81/ Dien khang ho cam: x12= ',num2str(x12),' (0m)']; disp(x);
x=['- Tinh theo don vi tuong doi: x12phay=',num2str(x12phay),];disp(x);
%82/ TINH LAI KE.
ke=(220-Iu*x1)/220;
x=['82/ Tinh lai kE: ke= ',num2str(ke),];disp(x);
x=['--> Ta thay tri so nay khong khac nhieu so voi gia tri ban dau kE=
',num2str(kE),' --> Nen khong can tinh lai'];disp(x);
disp(' ');

% CHUONG 6: TON HAO THEP VA TON HAO CO.


disp(' Chuong 6: TON HAO THEP VA TON HAO CO');
disp(' ---------*****---------');
%83/ TRONG LUONG RANG STATO.
hphayz1=hr1/10-d2s/30;
Gz1=7.8*Z1*bz1*hphayz1*Ls*kc*0.001;
x=['83/ Trong luong rang stato: Gz1= ',num2str(Gz1),' (kg)'];disp(x);
x=['- hphayz1= ',num2str(hphayz1),];disp(x);
%84/ TRONG LUONG GONG TU STATO.
Gg1=7.8*Ls*Lg1*hg1*pp*kc*0.001;
x=['84/ Trong luong gong tu stato: Gg1=',num2str(Gg1),'(kg)'];disp(x);
%85/ TON HAO SAT TRONG LOI SAT STATO.
PFez1=1.8*2.5*Bz1*Bz1*Gz1*0.001; % TRONG RANG
PFeg1=1.6*2.5*Bg1*Bg1*Gg1*0.001; % TRONG GONG
PphayFe=fix(100*(PFez1+PFeg1)); % TRONG CA LOI SAT STATO
PphayFe=PphayFe/100;
disp('85/ Ton hao sat trong loi sat stato:');
x=['- Trong rang: PFez1= ',num2str(PFez1),' (kW)'];disp(x);
disp(Bz1);
x=['- Trong gong: PFeg1= ',num2str(PFeg1),' (kW)'];disp(x);
disp(Bg1);
x=['- Trong ca loi sat stato: PphayFe= ',num2str(PphayFe),' (kW)'];disp(x);
%86/ TON HAO BE MAT TREN RANG ROTO.
k=b41/delta;
load('hinh6_1.mat'); % TR.141
for i=1:9
if hinh6_1(i,1)<=k && hinh6_1(i+1,1)>=k
g1=hinh6_1(i,1);
g2=hinh6_1(i+1,1);
h5=hinh6_1(i,2);
h6=hinh6_1(i+1,2);
end
end
beta0=(k-g1)*(h6-h5)/(g2-g1)+h5;
Bo=beta0*kS*Bs;
pbm=0.5*2*((Z1*ndb/10000)^1.5)*((10*Bo*t1)^2);
Pbm=pp*to*((t2-b42/10)/t2)*Ls*pbm*(10^-7);
x=['86/ Ton hao be mat tren rang roto: Pbm= ',num2str(Pbm),' (kW)'];disp(x);
disp('- Trong do:');
x=['- pbm= ',num2str(pbm),' khi Bo= ',num2str(Bo),' khi tra hinh 6.1/tr.141
duoc: beta0= ',num2str(beta0),];disp(x);
%87/ TON HAO DAP MACH TREN RANG ROTO.
Bdm=(v1*delta/10)*Bz2/(2*t2);
h13=hr2/10-dr2/30
bZ213=bZ213
Gz2=7.8*Z2*(hr2/10-dr2/30)*bZ213*Ls*kc*0.001;
Pdm=0.11*((Z1*ndb/1000*Bdm)^2)*Gz2*0.001;
x=['87/ Ton hao dap mach tren rang roto: Pdm=', num2str(Pdm),'(kW)'];disp(x);
disp('- Trong do:');
x=['- Bdm=', num2str(Bdm),'(T) & Gz2= ',num2str(Gz2),' (kg)'];disp(x);
%88/ TONG TON HAO THEP
PFe=fix(100*(PphayFe+Pdm+Pbm));
PFe=PFe/100;
x=['88/ Tong ton hao thep: PFe= ',num2str(PFe),' (kW)'];disp(x);
%89/ TON HAO CO.
Pco=fix(100*(((ndb/1000)^2)*((Dn/10)^4)*0.001));
Pco=Pco/100;
x=['89/ Ton hao co: Pco=', num2str(Pco),'(kW)'];disp(x);
%90/ TON HAO KHONG TAI
P0=PFe+Pco;
x=['90/ Ton hao khong tai: P0=', num2str(P0),'(kW)'];disp(x);
disp(' ');

% CHUONG 7: DAC TINH LAM VIEC


disp(' Chuong 7: DAC TINH LAM VIEC');
disp(' --------*****--------' );
disp(' Tu cac tinh toan tren ta duoc cac thong so sau:');
x=['- r1=', num2str(r1),'(Ohm)'];disp(x);
x=['- x1=', num2str(x1),'(Ohm)'];disp(x);
x=['- r2phay=', num2str(r2phay),'(Ohm)'];disp(x);
x=['- x2phay=', num2str(x2phay),'(Ohm)'];disp(x);
x=['- x12=', num2str(x12),'(Ohm)'];disp(x);
C1=1+x1/x12;
Idbx=Iu;
Idbr=(PFe*1000+3*Iu*Iu*r1)/(3*220);
E1=220-Iu*x1;
kItt=6*w1*kD/Z2;
Iphay2=Itd/kItt;
sdm=fix((Iphay2*r2phay/E1)*1000);
sdm=sdm/1000

%sdm=sdm/1000;
sm=(r2phay/(x1/C1+x2phay));
%sm=sm/1000;
x=['- C1=', num2str(C1),];disp(x);
x=['- Idbx=', num2str(Idbx),'(A)'];disp(x);
x=['- Idbr=', num2str(Idbr),'(A)'];disp(x);
x=['- E1=', num2str(E1),'(V)'];disp(x);
x=['- kItt=', num2str(kItt),];disp(x);
x=['- Iphay2=', num2str(Iphay2),'(A)'];disp(x);
x=['- sdm=', num2str(sdm),];disp(x);
x=['- sm=', num2str(sm),];disp(x);
% BANG DAT TINH LAM VIEC
bdattinh=zeros(18,6);
sbang=[0.005 0.01 sdm 0.025 0.1 sm];
for i=1:6
s=sbang(1,i);
bdattinh(1,i)=C1*C1*((r1/C1) + (r2phay/s));
%rns
bdattinh(2,i)=C1*C1*(x1/C1 + x2phay);
%xns
bdattinh(3,i)=sqrt(bdattinh(1,i)*bdattinh(1,i) +
bdattinh(2,i)*bdattinh(2,i)); %Zns
bdattinh(4,i)=C1*220/bdattinh(3,i);
%I'2
bdattinh(5,i)=bdattinh(1,i)/bdattinh(3,i);
%cosfi'2
bdattinh(6,i)=bdattinh(2,i)/bdattinh(3,i);
%sinfi'2
bdattinh(7,i)=Idbr + bdattinh(4,i)*bdattinh(5,i)/C1
;%I1r
bdattinh(8,i)=Idbx + bdattinh(4,i)*bdattinh(6,i)/C1;
%I1x
bdattinh(9,i)=sqrt(bdattinh(7,i)*bdattinh(7,i) +
bdattinh(8,i)*bdattinh(8,i)); %I1
bdattinh(10,i)=bdattinh(7,i)/bdattinh(9,i);
%cosfi
bdattinh(11,i)=3*220*bdattinh(7,i)*10.^-3;
%P1
bdattinh(12,i)=3*bdattinh(9,i)*bdattinh(9,i)*r1*10.^-3 ;
%PCu1
bdattinh(13,i)=3*bdattinh(4,i)*bdattinh(4,i)*r2phay*10.^-3;
%Pcu2
bdattinh(14,i)=0.005*bdattinh(11,i);
%Pf
bdattinh(15,i)=P0;
%P0
bdattinh(16,i)=bdattinh(12,i) + bdattinh(13,i) + bdattinh(14,i) +
bdattinh(15,i) ; %Tong P
bdattinh(17,i)=bdattinh(11,i)-bdattinh(16,i);
%P2
bdattinh(18,i)=bdattinh(17,i)*100/bdattinh(11,i);
%hieu suat

end
cnames = ['0.005 0.01 0.015 ', num2str(sdm),' 0.1 ',num2str(sm)];
rnames = ['rns(ohm) xns(ohm) Zns(ohm) Iphay2(A) cosfi2p sinfi2p I1r(A) I1x(A)
I1(A) cosfi P1(kW) Pcu1(kW) Pcu2(kW) Pf(kW) P0(kW) Tong_P(kW) P2(kW)
hieu_suat'];
printmat(bdattinh,'Dac tinh lam viec cua dong co dien khong dong bo roto long
soc',rnames,cnames);
% DO THI DAC TINH LAM VIEC
tong=0;
for sao=0.0001:0.0001:sdm
if sao>0
tong=tong+1;
end
i=tong;
rnsao=(C1^2)*(r1/C1+r2phay/sao);
xnsao=C1*C1*(x1/C1+x2phay);
Znsao=sqrt(rnsao^2+xnsao^2);
Ip2ao=C1*220/Znsao;
cosphip2ao=rnsao/Znsao;
sinphip2ao=xnsao/Znsao;
I1rao=Idbr+Ip2ao*cosphip2ao/C1;
I1xao=Idbx+Ip2ao*sinphip2ao/C1;
I11ao=sqrt(I1rao^2+I1xao^2);
cosphittao=I1rao/I11ao;
P1ao=3*220*I1rao*0.001;
Pcu1ao=3*I11ao*I11ao*r1*0.001;
Pcu2ao=3*Ip2ao*Ip2ao*r2phay*0.001;
Pfao=0.005*P1ao;
TongPao=Pcu1ao+Pcu2ao+Pfao+P0;
P2ao=P1ao-TongPao;
hsttao=P2ao*100/P1ao;
if hsttao<0
hsttao=0;
end
if P2ao<0
P2ao=0;
end
if cosphittao<0
cosphittao=0;
end
if I11ao<0
I11ao=0;
end
hstt(i)=hsttao;
P2(i)=P2ao;
cosphitt(i)=cosphittao*100;
stt(i)=sao*1000;
I11(i)=I11ao;
end
plot(P2,cosphitt,'b',P2,stt,'m');
hold on
plot(P2,hstt,'r',P2,I11,'k');
hold off
grid on
legend('cosphi','s','hs','I1');
xlabel('DO THI DAC TINH LAM VIEC CUA DONG CO KHONG DONG BO ROTO LONG SOC
3(kW) ');
ylabel('y=f(cosphi,s,hs,I1)')
% 91/ BOI SO MOMEN CUC DAI.
Ip2m=bdattinh(4,6);
Ip2dm=bdattinh(4,4);
Idm=bdattinh(9,4);
P2dm=bdattinh(17,4);
cosy=bdattinh(10,4);
hieusuat=bdattinh(18,4);
disp(' --> Cac gia tri dinh muc:');
x=[' Ip2dm= ',num2str(Ip2dm),' (A)'];disp(x);
x=[' Ip2m= ',num2str(Ip2m),'(A)'];disp(x);
x=[' Idm= ',num2str(Idm),' (A)'];disp(x);
x=[' P2dm= ',num2str(P2dm),'( kW)'];disp(x);
x=[' cosy= ',num2str(cosy),];disp(x);
x=[' hieusuat= ',num2str(hieusuat),' (%)'];disp(x);
mmax=(sdm/sm)*(Ip2m/Ip2dm)^2;
x=['91/ Boi so momen cuc dai: mmax= ',num2str(mmax),];disp(x);
% KIEM TRA
load('trammax.mat'); %tr.268
tkt=pp/2+1;
for i=1:19
if trammax(i,1)==P
mkt=trammax(i,tkt);
end
end
x=['--> Tra bang10.10/tr.268 ->Boi so momem cuc dai ung voi cong suat: P=
',num2str(P),'(kW) la mmax = ',num2str(mkt),];disp(x);
if (mmax>=mkt)
disp('--> Gia tri mmax THOA MAN yeu cau');
elseif mmax<mkt
disp('--> Gia tri mmax KHONG THOA MAN yeu cau');
end
disp(' ');

% CHUONG 8: TINH TOAN DAT TINH KHOI DONG


disp(' Chuong 8: TINH TOAN DAT TINH KHOI DONG');
disp(' ---------*****----------');
%92/ THAM SO CUA DONG CO KHI XET DEN HIEU UNG MAT NGOAI VOI s=1;
disp('92/ Tham so cua dong co dien khi xet den hieu ung mat ngoai voi s=1:
');
S=1;
a=ceil(hr2-h42);
exelon=ceil(0.067*a*10);
exelon=exelon/10;
load('hinh10_13.mat');
for i=1:13
if hinh10_13(i,1)<=exelon && hinh10_13(i+1,1)>=exelon
y5=hinh10_13(i,1);
y6=hinh10_13(i+1,1);
x5=hinh10_13(i,2);
x6=hinh10_13(i+1,2);
z5=hinh10_13(i,3);
z6=hinh10_13(i+1,3);
end
end
ci=(exelon-y5)*(x6-x5)/(y6-y5)+x5; % tra ci
phi=(exelon-y5)*(z6-z5)/(y6-y5)+z5; % tra phi
kR=1+phi;
rtde=kR*rtd;
r2e=rtde+(2*rv/(DEnta*DEnta)); % dien tro cua roto khi xet den hieu ung
mat ngoai voi s=1
r2phaye=gama*r2e; % dien tro roto da quy doi.
x=['- exelon= ',num2str(exelon),' Trong do: a= ',num2str(a),'(mm)'];disp(x);
x=['- Theo hinh 10-13/tr.256 , khi exelon= ',num2str(exelon),' tra ra: ci=
',num2str(ci),' ; phi= ',num2str(phi),];disp(x);
x=['- kR= ',num2str(kR),];disp(x);
x=['- rtde= ',num2str(rtde),' (0m)'];disp(x);
x=['- Dien tro toro khi xet den hieu ung mat ngoai voi s=1: r2e=
',num2str(r2e),];disp(x);
x=['- Dien tro roto da quy doi: r2phaye= ',num2str(r2phaye),' (0m)'];disp(x);
londar2e=(h1/(3*dr)*(1-(pi*dr^2/(8*Sr2))+0.66-b42/(2*dr)))*ci+h42/b42 ;
x=['- He so tu dan ranh roto khi xet den hieu ung mat ngoai voi s=1:
londar2e= ',num2str(londar2e),];disp(x);
Tonglonda2e=londar2e+londat2+londad2+londarn;
x=['- Tong he so tu dan roto khi xet den hieu ung mat ngoai voi s=1:
Tonglonda2e= ',num2str(Tonglonda2e),];disp(x);
x2phaye=x2phay*Tonglonda2e/Tonglonda2;
x=['- Dien khang roto khi xet den hieu ung mat ngoai: x2phaye=
',num2str(x2phaye),' (0m)'];disp(x);
disp('- Tong tro ngan mach khi xet den hieu ung mat ngoai voi s=1 :' );
rne=r1+r2phaye;
xne=x1+x2phaye;
zne=sqrt(rne^2+xne^2);
x=['- rne= ',num2str(rne),' (0m)'];disp(x);
x=['- xne= ',num2str(xne),' (0m)'];disp(x);
x=['- zne= ',num2str(zne),' (0m)'];disp(x);
Ine=220/zne;
x=['- Dong dien ngan mach khi chi xet den hieu ung mat ngoai: Ine=
',num2str(Ine),' (A)'];disp(x);
%93/
disp('93/ Tham so cua dong co dien khi xet ca hieu ung mat ngoaiva su bao hoa
cua mach tu tan khi s=1:' );
kbh=1.4;
disp('- So bo chon he so bao hoa: kbh=1.4');
Inbhe=kbh*Ine;
x=['- Dong dien ngan mach khi xet ca hieu ung mat ngoai va su bao hoa cua
mach tu tan: Inbhe= ',num2str(Inbhe),' (A)'];disp(x);
load('hinh10_14.mat'); % TR.259
for i=1:7
if hinh10_14(i,1)==beta
kb=hinh10_14(i,2);
break;
elseif hinh10_14(i,1)<=beta && hinh10_14(i+1,1)>=beta
aa1=hinh10_14(i,1);
aa2=hinh10_14(i+1,1);
bb1=hinh10_14(i,2);
bb2=hinh10_14(i+1,2);
end
end
kb=(beta-aa1)*(bb2-bb1)/(aa2-aa1)+bb1;
Fzbh=0.7*Inbhe*ur1/a1*(kb+ky*kD*(Z1/Z2));
Cbh=0.64+2.5*sqrt(Denta/(t1+t2));
x=['- Suc tu dong trung binh cua mot ranh stato: Fzbh= ',num2str(Fzbh),' ->
Trong do: kb tra theo hinh10-14/tr.259: kb= ',num2str(kb),];disp(x);
x=['- Cbh= ',num2str(Cbh),];disp(x);
Bfid=Fzbh*10^-4/(1.6*Cbh*Denta);
x=['- Bfid= ',num2str(Bfid),' (T)'];disp(x);
load('hinh10_15.mat');%tr.260
for i=1:11
if hinh10_15(i,1)<=Bfid && hinh10_15(i+1,1)>=Bfid
aaa1=hinh10_15(i,1);
aaa2=hinh10_15(i+1,1);
bbb1=hinh10_15(i,2);
bbb2=hinh10_15(i+1,2);
end
end
xd=(Bfid-aaa1)*(bbb2-bbb1)/(aaa2-aaa1)+bbb1;
C11=(t1-b41*0.1)*(1-xd);
%C11=round(C11);
%C11=C11/10;
Dlonda1bh=(h41*0.1+0.58*(dr1/20))/(b41*0.1)*C11/(C11+1.5*0.1*b41);
x=['- C1=',num2str(C11),' Trong do: Theo hinh 10-15/tr260 tra dc: xdeta=
',num2str(xd),];disp(x);
x=['- Dlonda1bh= ',num2str(Dlonda1bh),];disp(x);
londar1bh=londar1- Dlonda1bh;
x=['- He so tu dan tan ranh stato khi xet den bao hoa mach tu tan:
londar1bh= ',num2str(londar1bh),];disp(x);
londat1bh=londat1*xd;
x=['- He so tu tan tap stato khi xet den bao hoa mach tu tan: londat1bh=
',num2str( londat1bh),];disp(x);
Tonglonda1bh=(londar1bh+londat1bh+londad1)*100;
Tonglonda1bh=round(Tonglonda1bh);
Tonglonda1bh=Tonglonda1bh/100;
x=['- Tong he so tu tan tap stato khi xet den bao hoa mach tu tan:
Tonglonda1bh= ',num2str(Tonglonda1bh),];disp(x);
x1bh=x1*Tonglonda1bh/Tonglonda1;
C2=(t2-b42*0.1)*(1-xd);
Dlonda2bh=(h42*0.1)/(b42*0.1)*C2/(C2+0.1*b42);
x=['- Dien khang stato khi xet den bao hoa mach tu tan: x1bh=
',num2str(x1bh),' (0m)'];disp(x);
x=['- C2=',num2str(C2),];disp(x);
x=['- Dlonda2bh= ',num2str(Dlonda2bh),];disp(x);
londar2ebh=londar2e-Dlonda2bh;
x=['- He so tu dan tan ranh roto khi xet den bao hoa mach tu tan va hieu ung
mat ngoai: londar2ebh= ',num2str(londar2ebh),];disp(x);
londat2bh=londat2*xd;
x=['- He so tu tan tap roto khi xet den bao hoa mach tu tan: londat2bh=
',num2str( londat2bh),];disp(x);
londarnbh=londarn*xd;
x=['- He so tu tan do ranh nghieng roto khi xet den bao hoa mach tu tan:
londarnbh= ',num2str(londarnbh),];disp(x);
Tonglonda2ebh=londarnbh+londat2bh+londar2ebh+londad2;
x=['- Tong he so tu tan roto khi xet den bao hoa mach tu tan va hieu ung mat
ngoai: Tonglonda2ebh= ',num2str( Tonglonda2ebh),];disp(x);
x2phayebh=x2phay*Tonglonda2ebh/Tonglonda2;
%x2phayebh=round(x2phayebh);
%x2phayebh=x2phayebh/100;
x=['- Dien khang roto khi xet den hieu ung mat ngoai va su bao hoa cua mach
tu tan: x2phayebh= ',num2str(x2phayebh),' (0m)'];disp(x);
%94/
disp('94/ Cac tham so ngan mach khi xet den hieu ung mat ngoai va su bao hoa
cua mach tu tan:');
rnexelon=r1+r2phaye;
xnebh=x1bh+x2phayebh;
znebh=sqrt(rnexelon^2+xnebh^2);
x=['- rnexelon= ',num2str(rnexelon),' (0m)'];disp(x);
x=['- xnebh= ',num2str(xnebh),' (0m)'];disp(x);
x=['- znebh= ',num2str(znebh),' (0m)'];disp(x);
%95/ DONG DIEN KHOI DONG
Ik=220/znebh;
%Ik=fix(Ik);
x=['95/ Gia tri dong dien khoi dong: Ik= ',num2str(Ik),' (A)'];disp(x);
x=['--> Tri so dong dien khoi dong: Ik= ',num2str(Ik),'(A) khong khac nhieu
so voi gia thiet: Inbhe= ',num2str(Inbhe),' (A) --> nen khong can gia thiet
lai'];disp(x);
%96/ BOI SO DONG DIEN KHOI DONG
ik=Ik/Idm;
x=['96/ Boi so dong dien khoi dong: ik= ',num2str(ik),];disp(x);
load('tra_ik.mat');
tkt=pp/2+1;
for i=1:19
if tra_ik(i,1)==P
IK=tra_ik(i,tkt);
end
end
x=['--> Tra bang10.12/tr.271 ->Boi so ddkd ung voi cong suat: P=
',num2str(P),'(kW) la ik(bang) = ',num2str(IK),];disp(x);
if (ik>=IK)
disp('--> Gia tri ik KHONG THOA MAN yeu cau');
elseif ik<IK
disp('--> Gia tri ik THOA MAN yeu cau');
end
disp('- Dien khang ho cam khi xet den bao hoa:');
x12n=x12*ku;
C2ebh=1+x2phayebh/x12n;
Iphay2k=round(Ik/C2ebh);
x=['- x12n= ',num2str(x12n),];disp(x);
x=['- C2ebh= ',num2str(C2ebh),];disp(x);
x=['- Iphay2k= ',num2str(Iphay2k),];disp(x);
%97/ BOI SO MOMEN KHOI DONG
mk=(Iphay2k/Ip2dm)^2*r2phaye*sdm/r2phay;
x=['97/ Boi so momen khoi dong: mk= ',num2str(mk),];disp(x);
load('tra_mk.mat');
tkt=pp/2+1;
for i=1:19
if tra_mk(i,1)==P
MK=tra_mk(i,tkt);
end
end
x=['--> Tra bang10.11/tr.271 ->Boi so momem khoi dong ung voi cong suat: P=
',num2str(P),'(kW) la mk(bang) = ',num2str(MK),];disp(x);
if (mk>=MK)
disp('--> Gia tri mk THOA MAN yeu cau');
elseif mk<MK
disp('--> Gia tri mk KHONG THOA MAN yeu cau');
end
disp(' ');

disp(' Chuong 9: TRONG LUONG VAT LIEU TAC DUNG VA CHI TIEU SU
DUNG');
disp(' -----------*****-----------');
%106/ TRONG LUONG THEP SILIC CAN CHUAN BI
GFe=((Dn+delta)^2)*Ls*kc*7.8*0.001;
x=['106/ Trong luong thep silic can chuan bi: GFe=',
num2str(GFe),'(kg)'];disp(x);
%107 TRONG LUONG DONG CUA DAY QUAN STATO
GpCu=Z1*ur1*n1*s1*ltb*8.9*0.00001;
GCu=(0.876+0.124*((dcd/d)^2))*GpCu;
disp('107/ Trong luong dong cua day quan stato:');
x=['- Khi khong tinh cach dien: GpCu=', num2str(GpCu),'(kg)'];disp(x);
x=['- Khi tinh cach dien: GCu=', num2str(GCu),'(kg)'];disp(x);
%108/ TRONG LUONG NHOM ROTO
Gtd=Z2*Sr2*Ls*2.6*0.00001;
Gv=2*pi*0.1*Dv*Sv*2.6*0.00001;
GAl=Gtd+Gv;
disp('108/ Trong luong nhom roto:');
x=['- Trong luong nhom o thanh dan: Gtd=', num2str(Gtd),'(kg)'];disp(x);
x=['- Trong luong nhom o vanh ngan mach: Gv=', num2str(Gv),'(kg)'];disp(x);
x=['- Trong luong nhom o roto: GAl=', num2str(GAl),'(kg)'];disp(x);
%109/ CHI TIEU KINH TE VE VAT LIEU TAC DUNG
gFe=GFe/P;
gCu=GCu/P;
gAl=GAl/P;
disp('109/ Chi tieu kinh te ve vat kieu tac dung:');
x=['- Thep ky thuat dien: gFe=', num2str(gFe),'(kg/kW)'];disp(x);
x=['- Dong: gFe=', num2str(gCu),'(kg/kW)'];disp(x);
x=['- Nhom: gAl=', num2str(gAl),'(kg/kW)'];disp(x);
disp(' ');
disp('______________________________________________________*****________THE
END_________*****____________________________________________________________
');

You might also like