LT TBD Word
LT TBD Word
LT TBD Word
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(' ');
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(' ');
londar2=(h11/(3*dr1))*((1-((pi*dr1*dr1)/(8*Sr2)))^2)+(0.66-
(b42/(2*dr1)))*1+(h42/b42);
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(' ');
%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(' ');
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_________*****____________________________________________________________
');