Розрахунок потужності двигуна механізмів переміщення
3. Розрахунок за допомогою системи matlab 6.5.1
3.1. Порядок виконуваних дій
Вводимо м-файл
format bank
disp('7. ВИБІР ЕЛЕКТРИЧНОЇ МЕРЕЖІ КРАНУ')
disp('Введіть параметри крану')
Lb=input('Переміщення візка:\n Lb=');
Lm=input('Переміщення моста:\n Lm=');
G=input('Вантажність крану:\n G=');
G0=input('Маса вантажозахоплюючого пристрою:\n G0=');
Gb=input('Маса візка:\n Gb=');
Gm=input('Маса моста:\n Gm=');
tn=input('Час роботи механізму підйому:\n tn=');
tb=input('Час роботи механізму переміщення візка:\n tb=');
tm=input('Час роботи механізму переміщення моста:\n tm=');
tc=input('Час циклу:\n tc=');
disp('Введіть параметри двигунів крану')
Pcmn=input('Розрахункова потужність двигуна механізму підйому:\n Pcmn=');
Pn=input('Потужність двигуна механізму підйому:\n Pn=');
cn=input('cosfi двигуна механізму підйому:\n cn=');
hn=input('ККД двигуна механізму підйому:\n hn=');
Pcmb=input('Розрахункова потужність двигуна механізму переміщення візка:\n Pcmb=');
Pb=input('Потужність двигуна механізму переміщення візка:\n Pb=');
cb=input('cosfi двигуна механізму переміщення візка:\n cb=');
hb=input('ККД двигуна механізму переміщення візка:\n hb=');
Pcmm=input('Розрахункова потужність двигуна механізму переміщення візка:\n Pcmm=');
Pm=input('Потужність двигуна механізму переміщення моста:\n Pm=');
cm=input('cosfi двигуна механізму переміщення моста:\n cm=');
hm=input('ККД двигуна механізму переміщення моста:\n hm=');
Pn100=Pn*(100/40)^(1/2)/hn
Pb100=Pb*(100/40)^(1/2)/hb
Pm100=Pm*(100/40)^(1/2)/hm
tgn=(1-cn^2)^(1/2)/cn
tgb=(1-cb^2)^(1/2)/cb
tgm=(1-cm^2)^(1/2)/cm
Qn100=Pn100*tgn
Qb100=Pb100*tgb
Qm100=Pm100*tgm
P100=Pn100+Pb100+Pm100
Q100=Qn100+Qb100+Qm100
Sn100=(Pn100^2+Qn100^2)^(1/2)
Sb100=(Pb100^2+Qb100^2)^(1/2)
Sm100=(Pm100^2+Qm100^2)^(1/2)
S100=(P100^2+Q100^2)^(1/2)
disp('7.3. Розрахунок низьковольтної мережі')
Ip1=Sn100/3^(1/2)/0.4
Ip2=Sb100/3^(1/2)/0.4
Ip3=Ip1+Ip2
Ip4=Sm100/3^(1/2)/0.4
Ip5=Ip3+Ip4
disp('7.3.2. Вибір і розрахунок кабельної лінії КЛ1')
disp('7.3.2.1. Струм у кабельній лінії')
IKL1=Ip5
LKL1=0.1
xKL1=0.07
IKL2=Ip3;
LKL2=(Lb+3);
xKL2=0.07;
In=Pn/3^(1/2)/0.4/cn/hn;
Inyck=3*In;
Pcym=Pn+Pb+Pm;
Pcnow=Pcym/hn;
nef=Pcym^2/(Pn^2+Pb^2+Pm^2);
Imax=((Pcnow*0.5)^2+(Pcnow*1.5)^2)^(1/2)/3^0.5/0.4;
Im=Inyck+Imax-In;
Lmr=(Lm+1);
xmr=0.164;
rmr=0.315;
Umr=3^0.5*Im*Lmr*1e-3*(rmr*0.5+xmr*0.87);
Ucym=39;
while (Ucym>38)
if IKL1<65
rKL1=3.15;
elseif IKL1<90
rKL1=1.97;
elseif IKL1<115
rKL1=1.26;
elseif IKL1<135
rKL1=0.91;
elseif IKL1<165
rKL1=0.63;
elseif IKL1<200
rKL1=0.45;
else
rKL1=0.33;
end
UKL1=3^0.5*Ip5*LKL1*(rKL1*0.75+xKL1*0.661);
if IKL2<65
rKL2=3.15;
elseif IKL2<90
rKL2=1.97;
elseif IKL2<115
rKL2=1.26;
elseif IKL2<135
rKL2=0.91;
elseif IKL2<165
rKL2=0.63;
elseif IKL2<200
rKL2=0.45;
else
rKL2=0.33;
end
UKL2=3^0.5*Ip3*LKL2*1e-3*(rKL2*0.75+xKL2*0.661);
Ucym=12.42+UKL1+UKL2+Umr;
if (Ucym<38)
disp('Cумарні втрати напруги меньші за припустимі.')
else
IKL1=IKL1+20;
IKL2=IKL2+5;
end
end
disp('Кабельна лінія КЛ1')
if IKL1<65
disp('АВВГ 3x10+1x6, Iдоп = 65 А')
rKL1=3.15
elseif IKL1<90
disp('АВВГ 3x16+1x10, Iдоп = 90 А')
rKL1=1.97
elseif IKL1<115
disp('АВВГ 3x25+1x16, Iдоп = 115 А')
rKL1=1.26
elseif IKL1<135
disp('АВВГ 3x35+1x16, Iдоп = 135 А')
rKL1=0.91
elseif IKL1<165
disp('АВВГ 3x50+1x25, Iдоп = 165 А')
rKL1=0.63
elseif IKL1<200
disp('АВВГ 3x70+1x35, Iдоп = 200 А')
rKL1=0.45
else
disp('АВВГ 3x95+1x50, Iдоп = 240 А')
rKL1=0.33
end
UKL1
disp('7.3.3. Вибір і розрахунок кабельної лінії КЛ2')
disp('7.3.3.1. Струм у кабельній лінії')
Ip3
disp('Рекомендована довжина кабельної лінії КЛ2')
LKL2=(Lb+3)
xKL2=0.07
disp('Кабельна лінія КЛ2')
if IKL2<65
disp('АВВГ 3x10+1x6, Iдоп = 65 А')
rKL2=3.15
elseif IKL2<90
disp('АВВГ 3x16+1x10, Iдоп = 90 А')
rKL2=1.97
elseif IKL2<115
disp('АВВГ 3x25+1x16, Iдоп = 115 А')
rKL2=1.26
elseif IKL2<135
disp('АВВГ 3x35+1x16, Iдоп = 135 А')
rKL2=0.91
elseif IKL2<165
disp('АВВГ 3x50+1x25, Iдоп = 165 А')
rKL2=0.63
elseif IKL2<200
disp('АВВГ 3x70+1x35, Iдоп = 200 А')
rKL2=0.45
else
disp('АВВГ 3x95+1x50, Iдоп = 240 А')
rKL2=0.33
end
UKL2
disp('7.5. Вибір тролея.')
In=Pn/3^(1/2)/0.4/cn/hn
Inyck=3*In
Pcym=Pn+Pb+Pm
Pcnow=Pcym/hn
nef=Pcym^2/(Pn^2+Pb^2+Pm^2)
Imax=((Pcnow*0.5)^2+(Pcnow*1.5)^2)^(1/2)/3^0.5/0.4
disp('7.5.7. Тимчасовий пусковий струм під час пуску найпотужнішого двигуна')
Im=Inyck+Imax-In
disp('7.5.9. Втрати напруги у тролейному шинопроводі')
disp('Рекомендована довжина тролея')
Lmr=(Lm+1)
xmr=0.164
rmr=0.315
Umr=3^0.5*Im*Lmr*1e-3*(rmr*0.5+xmr*0.87)
disp('7.6. Сумарні втрати напруги в низьковольтній мережі')
Ucym=12.42+UKL1+UKL2+Umr
disp('де')
UKL1
UKL2
Umr
disp('7.1. Визначимо коефіцієнти завантаження відповідних двигунів при номінальному навантаженні')
K1=Pcmn/Pn
K2=Pcmb/Pb
K3=Pcmm/Pm
disp('7.2. Визначимо коефіцієнти завантаження відповідних двигунів без вантажу')
K31=K1*0.3
K32=K2*(G0+Gb)/(G+Gb)
K33=K3*(G0+Gm)/(G+Gm)
hb1=kz(K1)
hb2=kz(K2)
hb3=kz(K3)
hb31=kz(K31)
hb32=kz(K32)
hb33=kz(K33)
hf1=hb1*hn
hf31=hb31*hn
hf2=hb2*hb
hf32=hb32*hb
hf3=hb3*hm
hf33=hb33*hm
Wg=(Pn*K1*(tn+tn)/hf1+Pn*K31*(tn+tn)/hf31+Pb*K2*tb/hf2+Pb*K32*tb/hf32+Pm*K3*tm/hf3+Pm*K33*tm/hf33)/tc
disp('7.3. Кількість енергії, що споживається за рік')
Wp=Wg*5000
clear