0% found this document useful (0 votes)
72 views26 pages

LT TBD Word

This document contains the design calculations for a 3-phase induction motor with the following specifications: 1) Rated power of 1.1-75 kW 2) Number of pole pairs from 2-8 3) Efficiency and power factor inputs from the user 4) Calculates key dimensions such as outer diameter, inner diameter, stack length, number of slots, etc. using standard reference tables. The results are output for the user to select values.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
72 views26 pages

LT TBD Word

This document contains the design calculations for a 3-phase induction motor with the following specifications: 1) Rated power of 1.1-75 kW 2) Number of pole pairs from 2-8 3) Efficiency and power factor inputs from the user 4) Calculates key dimensions such as outer diameter, inner diameter, stack length, number of slots, etc. using standard reference tables. The results are output for the user to select values.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 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