Поле и емкость двухпроводной линии с учетом влияния Земли

Поле и емкость двухпроводной линии с учетом влияния Земли Рис. 8. Дано: над плоской поверхностью Земли подвешены горизонтально два цилиндрических провода с параллельными осями (рис. 8). h1 – высота подвеса 1-го провода; h2 – высота подвеса 2-го провода;    R – радиусы проводов; d – расстояние между нормальными проекциями осей проводов на поверхности Земли. По условию задачи требуется: вывести уравнения, связывающие между собой линейные плотности зарядов на проводах  и потенциалы проводов. Определить параметры этих уравнений: потенциальные и емкостные коэффициенты, частичные емкости и рабочую емкость линии, если d, h1 и h2>>R. Для решения поставленной задачи можно воспользоваться методом изображений. Распределение поля над поверхностью Земли не изменится, если Землю убрать, а под поверхностью Земли расположить на глубинах h1 и h2 провода с линейной плотностью заряда . После такого преобразования можно считать, что в системе действует электростатическое поле двух пар параллельных разноименно заряженных осей (рис. 9). Рис. 9. Поскольку d, h и h >> R, смещением электрических осей относительно геометрических осей можно пренебречь. Используя принцип наложения, выразим потенциалы проводов через линейные плотности зарядов   Из этих уравнений видно, что потенциалы проводов являются линейными комбинациями линейных плотностей зарядов или Коэффициентыназываются потенциальными коэффициентами единицы длины проводов. ,– это собственные потенциальные коэффициенты проводов, , – это взаимные потенциальные коэффициенты. Как видно, матрица симметричная, значит, для линии выполняется принцип взаимности. Из системы уравнений (1) выразим τ1 и τ2. Коэффициентыназываются емкостными коэффициентами на единицу длины линии и измеряются в Ф/м. Собственные потенциальные и емкостные коэффициенты всегда положительны. Взаимные потенциальные коэффициенты положительны, а взаимные емкостные коэффициенты всегда отрицательны. Систему уравнений (2) можно записать иначе Коэффициенты Сij называют частичными емкостями на единицу длины. Если провода линии не связаны с Землей и питаются от незаземленного источника ЭДС, то суммарный заряд линии равен нулю, т.е Вычтем второе уравнение из первого и получим Отношение линейной плотности заряда провода к напряжению называют в данном случае рабочей ёмкостью линии на единицу длины Можно изобразить эквивалентную схему системы заряженных проводников линии (рис. 10). Рис. 10. Анализируя эту схему, можно получить другое выражение для рабочей емкости линии Cраб = C12 + C11C22/(C11 + C22)                                                  (4) Можно доказать, что выражения (3) и (4) тождественны. Ниже представлен текст вычислительного сценария расчёта потенциальных и ёмкостных коэффициентов, а также частичных ёмкостей многопроводной воздушной линии с учётом влияния земли. % ElStatLin - Расчёт потенциальных и емкостных коэффициентов многопроводной линии. % Смещение электрических осей относительно геометрических не учитываетсq. % Входные параметры: % x - горизонтальные координаты подвеса проводов; % y - вериткальные координаты подвеса проводов; % D - диаметры всех проводов % Все эти переменные - строковые матрицы % Выходные параметры: % al - потенциальные коэффициенты проводов; % be - ёмкостные коэффициенты проводов; % c  - частичные ёмкости проводов   eps0=8.85e-12;  % Абcолютнаq диэлектрическаq проницаемость вакуума, Ф/м rp=sqrt((repmat(x,length(x),1)-repmat(x,length(x),1).').^2+...  (repmat(y,length(y),1)-repmat(y,length(y),1).').^2)+diag(D/2); rm=sqrt((repmat(x,length(x),1)-repmat(x,length(x),1).').^2+...  (repmat(y,length(y),1)+repmat(y,length(y),1).').^2); al=log(rm./rp)/eps0/2/pi             % Потенциальные коэффициенты, м/Ф be=inv(al)                           % Ёмкостные коэффициенты, Ф/м c=diag(sum(be))+diag(diag(be))-be    % Частичные ёмкости, Ф/м