Испытания схем увязки САУТ-ЦМ и МПЦ в лаборатории
Дипломная работа - Компьютеры, программирование
Другие дипломы по предмету Компьютеры, программирование
?ы увязки на ЭВМ.
Расчет надежности с помощью ЭВМ производился по следующей программе, написанной на языке программирования Pascal:
Program Nadeg;V, I, Z: longint;, LL, II:array [1..5] of real;:array [1..200] of real;, S0, T, P: real;:string;:text;(lst, result.txt); rewrite(lst);:=0;[1]:=12; LL[1]:=6.00e-8;[2]:=5; LL[2]:=1.32e-7;[3]:=3; LL[3]:=2.50e-7;[4]:=2; LL[4]:=2.50e-7;[5]:=44; LL[5]:=8.40e-9;(коэффициент нагрузки K=); readln(K);
write (время наработки V=); readln(V);
str (k, s1); writeln (lst, коэффициент нагрузки K=+s1);(V, s1); writeln (lst, время наработки V=+s1);(lst, I=1-резистор; I=2-конд-р; I=3-транзистор;,
I=4-индуктивность; I=5-пайка);
writeln (lst, введем количество элементов N(I));(lst, введем интенсивность отказов L(I));
for I:=1 to 5 do begin(lst, I=, I, N(I)=, NN[I]:5:0, L(I)=, LL[I]:10);[I]:=NN[I]*LL[I];:=S0+II[I];;; writeln;(lst, время безотказной работы);:=1/(K*S0);(lst, T=, T:1:1);
writeln (lst, вероятность и время наработки (час));
for Z:=0 to (V DIV 1000) do begin:=exp (-K*S0*Z*1000);(lst, P=, P:1:6, Z=, Z*1000);;
close(lst);.
6. Расчёт параметров кабеля iБ
6.1 Математическая модель расчета первичных параметров кабелей iБ
Для разработки схемы увязки устройств САУЦ-ЦМ и МПЦ необходимо рассчитать параметры кабеля iБ на частотах до 30 кГц.
Расчет первичных параметров кабелей iБ удобнее рассчитывать с помощью ЭВМ. Для данных расчетов будем использовать TURBO РACKAL 7.0. Программа будет использовать следующие формулы:
R=R0+RПЭ+RБЛ+RМ+RЗ,
R0=23,3 Ом/км.
Удельное сопротивление земли:
RЗ=1/?з ? r2;
где ?з - удельная проводимость земли =10-4;
r - радиус жилы кабеля =0,6 мм.
Сопротивления за счет эффекта близости:
,
где d - диаметр проводника (жилы кабеля) =1,2 мм,
а - расстояние между жилами,
G(kr), H(kr) - табличные значения функций Бесселя,
kr - коэффициент, зависящий от материала проводника и для меди, равный .
Сопротивление за счет поверхностного эффекта:
RПЭ=?2f 10 - 4, Ом /км.
Километрическая индуктивность может быть определена по следующим формулам:
где
Емкость жилы кабеля относительно земли можно найти по формуле:
,
где ? - относительная диэлектрическая проницаемость оболочки кабеля, для полиэтилена =2,39
h-высота подвески,
Километрическая проводимость может быть представлена в следующем виде:
G= ? C tg?,
где tg? - тангенс потерь в диэлектрике покрытия, который для полиэтиленовых покрытий может быть принят 0,0004 для частот до 10кГц.
Рассмотрим также изменение вторичных параметров в зависимости от частоты. Для этого воспользуемся приближенными формулами, полученными при условиях R<<?L G<<?C, откуда:
Коэффициент затухания:
, дБ/ км.
Коэффициент фазы:
, Рад/км.
Волновое сопротивление:
, Ом.
Скорость распространения:
, км/с.
6.2 Программа расчета параметров
Данная программа производит расчет первичных и вторичных параметров по приведенным выше формулам. Расчет осуществляется в цикле: от 1 до 30 кГц. Результаты расчетов выводятся в файл в виде таблицы.
Program parametr;
Var
Rp, Rm, Lp, Lm, Cp, Cm, Gp, Gm, ap, am, bp, bm, kr, Fkr, Gkr, Hkr, Zvp, Zvm,
Vp, Vm, tp, tm, sigma, Rpe, Rbe1, Rbe2, Ro, Rz, Qkr:real;:integer;:text;(zz, A:\otvet.pas);:=0;(zz);(zz, f , Rp , Rm , Lp , Lm , Cp ,
Cm , Gp , Gm , ap , am , bp ,
bm , Zvp , Zvm , Vp , Vm , tp ,
tm );f:=1 to 30 do:=0.0126*sqrt (1000*f);:=2*sqrt(2)/kr;f=1 then Fkr:=0.000325; If f=1 then Gkr:=0.000975; If f=1 then Hkr:=0.04;f=2 then Fkr:=0.000326; If f=2 then Gkr:=0.000976; If f=2 then Hkr:=0.042;f=3 then Fkr:=0.00039; If f=3 then Gkr:=0.01000; If f=3 then Hkr:=0.049;f=4 then Fkr:=0.0048; If f=4 then Gkr:=0.01480; If f=4 then Hkr:=0.050;f=5 then Fkr:=0.005; If f=5 then Gkr:=0.015; If f=5 then Hkr:=0.051;f=6 then Fkr:=0.00519; If f=6 then Gkr:=0.01519; If f=6 then Hkr:=0.053;f=7 then Fkr:=0.0052; If f=7 then Gkr:=0.0152; If f=7 then Hkr:=0.057;f=8 then Fkr:=0.006; If f=8 then Gkr:=0.02; If f=8 then Hkr:=0.059;f=9 then Fkr:=0.01; If f=9 then Gkr:=0.040; If f=9 then Hkr:=0.065;f=10 then Fkr:=0.02; If f=10 then Gkr:=0.05; If f=10 then Hkr:=0.075;f=11 then Fkr:=0.021; If f=11 then Gkr:=0.06; If f=11 then Hkr:=0.08;f=12 then Fkr:=0.023; If f=12 then Gkr:=0.061; If f=12 then Hkr:=0.083;f=13 then Fkr:=0.025; If f=13 then Gkr:=0.065; If f=13 then Hkr:=0.089;f=14 then Fkr:=0.0258; If f=14 then Gkr:=0.0691; If f=14 thenHkr:=0.092;f=15 then Fkr:=0.03; If f=15 then Gkr:=0.07; If f=15 then Hkr:=0.095;f=16 then Fkr:=0.04; If f=16 then Gkr:=0.09; If f=16 then Hkr:=0.98;f=17 then Fkr:=0.046; If f=17 then Gkr:=0.095; If f=17 then Hkr:=0.1;f=18 then Fkr:=0.05; If f=18 then Gkr:=0.1; If f=18 then Hkr:=0.11;f=19 then Fkr:=0.07; If f=19 then Gkr:=0.14; If f=19 then Hkr:=0.14;f=20 then Fkr:=0.071; If f=20 then Gkr:=0.141; If f=20 then Hkr:=0.141;f=21 then Fkr:=0.074; If f=21 then Gkr:=0.15; If f=21 then Hkr:=0.149;f=22 then Fkr:=0.075; If f=22 then Gkr:=0.16; If f=22 then Hkr:=0.15;f=23 then Fkr:=0.078; If f=23 then Gkr:=0.165; If f=23 then Hkr:=0.16;f=24 then Fkr:=0.078; If f=24 then Gkr:=0.17; If f=24 then Hkr:=0.161;f=25 then Fkr:=0.0782; If f=25 then Gkr:=0.1724; If f=25 then Hkr:=0.169;f=26 then Fkr:=0.0785; If f=26 then Gkr:=0.173; If f=26 then Hkr:=0.17;f=27 then Fkr:=0.079; If f=27 then Gkr:=0.174; If f=27 then Hkr:=0.18;f=28 then Fkr:=0.085; If f=28 then Gkr:=0.18; If f=28 then Hkr:=0.19;f=29 then Fkr:=0.09; If f=29 then Gkr:=0.2; If f=29 then Hkr:=0.2;f=30 then Fkr:=0.1; If f=30 then Gkr:=0.21; If f=30 then Hkr:=0.21;:=Gkr*sqr (1.2/4.34)/(1-Hkr*sqr (1.2/4.34));:=Gkr*sqr (1.2/0.7)/(1-Hkr*sqr (1.2/0.7));:=31.5;:=sqr(pi)*f*0.1;:=1/(0.01*sqr (0.6)*pi);:=Ro+Rbe1+Rpe;:=Ro+Rbe2+Rpe+Rz;:=1/sqrt (pi*1000*f*0.01);:=(4*ln((4.34-0.6)/0.6)+Qkr)*0.0001;:=(1.336-0.1*ln (1000*f*sigma))*0.001;:=0.000001*1.25/(36*ln (4.34/0.6));:=0.000000001*55.5/(23.9*ln (0.7/0.6));:=2*pi*1000*f*Cp*0.0004;:=2*pi*1000*f*Cm*0.0004;:=(Rp/2*sqrt (Cp/Lp)+Gp/2*sqrt (Lp/Cp));:=(Rm/2*sqrt (Cm/Lm)+Gm/2*sqrt (Lm/Cm));:=2*pi*1000*f*sqrt (Lp*Cp);:=2*pi*1000*f*sqrt (Lm*Cm);:=sqrt (Lp/Cp);:=sqrt (Lm/Cm);:=2*pi*1000*f/bp;:=2*pi*1000*f/bm;:=2/Vp;:=2/Vm;(zz, , f, , Rp:2:3, , Rm:2:3, , Lp:2:9, , Lm:2:9, , Cp:1:11, , Cm:1:11, , Gp:1:8,
, Gm:2:8, , ap:2:4, , am:2:2, , bp:2:3, , bm:2:3, , Zvp:4:2, , Zvm:4:2,
, Vp:5:3, , Vm:5:3, , tp:2:6, , tm:2:7);;
close(zz);
end.
.3 Расчёт параметров кабеля методом Х.Х и К.З