В качестве примера составим математическую модель, позволяющую исследовать электромагнитные переходные процессы, возникающие в схеме при ударе молнии в линию в непосредственной близости от кабельной вставки, в конце которой подключен ненагруженный трансформатор. На рис.3 представлены исходная и расчетная схемы замещения, содержащая три узла и связывающие эти узлы отрезки воздушной и кабельной линий. В схеме приняты следующие исходные данные:
e(t) - воздействующий импульс напряжения конечной длительности,
l1= 1500 м и l2= 600 м - соответственно, длины воздушной и кабельной линии,
n1= 300 м/мкс и n2= 150 м/мкс - скорости распространения электромагнитной волны по воздушной и кабельной линии,
Z = ZВ1= 400 Ом и ZВ2 = 40 Ом - волновые сопротивления ВЛ и КЛ,
Стр = 2000 пФ - входная емкость трансформатора на частоте грозового импульса напряжения.
Исходные уравнения для узлов записываются в виде:
узел 1:
V21(t) = u1(t) + i21(t) ZВ1 ,
u1(t) = e(t) – i(t) Z , (4)
i(t) = – i21(t),
W12(t) = 2u1(t) – V21(t).
а)
б)
Рис. 3. Исходная (а) и расчетная (б) схемы замещения для анализа
волновых процессов в схеме с распределенными параметрами
узел 2:
V12(t) = u2(t) + i12(t) ZВ1 ,
V32(t) = u2(t) + i32(t) ZВ2 ,
i12(t) = – i32(t), (5)
W21(t) = 2×u2(t) – V12(t) ,
W23(t) = 2×u2(t) – V32(t) .
узел 3:
V23(t) = u3(t) + i23(t) ZВ2 ,
i23(t) = Стр du3(t)/dt , (6)
W32(t) = 2u3(t) – V23(t) .
При неучете потерь в воздушной и кабельной линии приходящие в узлы волны определятся как:
V12(t) = W12(t–t12),
V21(t) = W21(t–t21) ,
V23(t) = W23(t–t23), (7)
V32(t) = W32(t–t32),
где t12 и t23 - времена пробега волны по воздушной и кабельной линии, соответственно.
Для определения u1(t), u2(t) и u3(t) системы уравнений (4)…(6) целесообразно переписать в виде:
u1(t) = [V21(t) Z + e(t) ZВ1]/(Z + ZВ1) ,
u2(t) = [V12(t) ZВ2 + V32(t) ZВ1]/(ZВ1 + ZВ2) ,
du3(t)/dt = [V23(t) – U3(t)]/Cтр/ZВ2 , (8)
W12(t) = 2u1(t) – V21(t) , W21(t) = 2u2(t) – V12(t) ,
W23(t) = 2u2(t) – V32(t) , W32(t) = 2u3(t) – V23(t) .
При использовании для численного интегрирования обыкновенного дифференциального уравнения dy(t)/dt = f [y(t), t], например, неявного метода трапеций, решение уравнения на (k+1)-м шаге интегрирования (t= k×h) запишется в виде
y[(k+1)h] = y(kh) + 0,5×h [f(yk+1,tk+1) + f(yk,tk)] . (9)
Соответственно, напряжение в узле 3 на (k+1)-ом шаге интегрирования можно определить как
u3(k+1) = [u3(k) + (V23(k) +V23(k+1) – u3(k))×h/2/ZВ2/Стр]/(1+ h/2/ZВ2/Cтр). (10)
Таким образом, расчет переходного процесса в исследуемой схеме сводится к решению системы уравнений (8) при известных значениях волн, приходящих в узлы схемы, определению волн, отраженных от узлов и имитации движения обобщенных волн по линиям в соответствии с волновыми уравнениями (7).
Для реализации этой модели может быть использована какая-либо из систем инженерных расчетов: MATLAB, Mathcad, EXCEL и др. Ниже приведены примеры использования для реализации рассмотренного алгоритма систем MATLAB и EXCEL
Пример реализации математической модели в системе MATLAB
Программа расчета переходных процессов в схеме, приведенной на рис.3, при использовании метода бегущих волн имеет вид:
% Исходные данные
tr=40e-6; - время расчета переходного процесса
h=0.1e-6; - шаг расчета
L1=1500;L2=600; - длина ВЛ и КЛ
v1=300.e6;v2=150.e6; - скорости распр. волн по ВЛ и КЛ
Z1=400;Z2=40; - волновые сопротивления ВЛ и КЛ
C=2e-9; - входная емкость трансформатора
U0=1; - амплитуда воздействующей волны
ti=2e-6; - длительность срезанного импульса
U3k=0; - начальное значение U3(t)
nt=tr/h; - количество циклов по времени
% Начальные значения массивов напряжений и времени
UM1=0;UM2=0;UM3=0;tm=0;
% Формирование размерности для массивов волн
N1=round(1+L1/v1/h);
N2=round(1+L2/v2/h);
% Формирование начальных (нулевых) массивов волн
V12=zeros(1,N1);
V21=zeros(1,N1);
V23=zeros(1,N2);
V32=zeros(1,N2);
% Расчет волнового процесса (цикл по времени)
for k=1:nt
t=k*h;
% if(t<=ti)
% e=2*U0-2*U0*(ti-t)/ti;
% else - срезанный импульс
% e=0;
% end
e=2*U0; - прямоугольный импульс
% Сдвиг волн на шаге расчета
for i=N1:-1:2
V12(i)=V12(i-1);
V21(i)=V21(i-1);
end
for i=N2:-1:2
V23(i)=V23(i-1);
V32(i)=V32(i-1);
end
% Вычисление значений напряжения в узлах
U1=(V21(N1)+e)/2;
U2=(V12(N1)*Z2+V32(N2)*Z1)/(Z1+Z2);
U3=(U3k+h*(V23(N2)+V23(N2-1)-U3k)/2/Z2/C)/(1+h/2/Z2/C);
U3k=U3;
% Формирование массивов напряжений и времени
UM1=[UM1 U1];
UM2=[UM2 U2];
UM3=[UM3 U3];
tm=[tm t*1.e6];
% Вычисление отраженных (падающих) волн
V12(1)=2*U1-V21(N1);
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.