exponenta event banner

Решатель проводов

wire solver

Ниже приведено сокращенное описание выведения уравнений, решаемых решателем проводов. Для получения более подробной информации см. [1].

Для проводников PEC (Perfect electric conducting) граничные условия требуют, чтобы тангенциальная составляющая общего электрического поля на поверхности провода была равна нулю:

(E + Ei)  tan = 0 (1)

Можно также выразить уравнение в терминах плотности поверхностного тока, J→s, на проводе PEC и его пространственной производной. Это делается путем интегрирования эффекта поверхностной плотности тока с использованием функции Грина, g (|r→−r→'|). Функция Грина абстрактно представляет влияние бесконечно малого источника поверхностно-токовой плотности, расположенного в r→ 'на электрическое и магнитное поля в точке наблюдения r→. Функция Грина для распространения свободного пространства:

g (|r→−r→'|) = exp (−jk'r→−r→'|)      4π'r→−r→'| (2)

где k = startα0pci0 - волновое число, а λ - угловая частота. В общем, для решения электромагнитной (ЭМ) проблемы и получения ЭМ полей в среде необходимы следующие шаги. Для точного представления ЭМ-проблемы на поверхности провода принимаются во внимание как продольные, так и поперечные компоненты интегрированного тока и требуется, чтобы величина (1) поддерживалась по всей поверхности провода. Дискретизируйте Eq. (1) над конечным множеством точек или базисных функций, чтобы получить множество уравнений. Дискретизируйте поверхностную плотность тока в конечное множество базисных функций, чтобы получить множество неизвестных. Решите набор уравнений неизвестных, чтобы получить дискретное приближение к поверхностным токам-плотностям на проводе. Эти плотности тока позволяют вычислять электромагнитные поля в любой требуемой точке. Для тонких проводов указанные выше шаги являются исчерпывающими и упрощают функцию Green для:

g (ra) = exp (jkra )    4āra (3)

где ra = r2 + a2 - приблизительное среднее расстояния |r→−r→'|, r - расстояние между бесконечно малым поверхностным источником плотности тока и осью провода, a - радиус провода. Кроме того, плотность тока на проводе заменяется I=2πa'J→s|| тока провода. Используя вышеприведенное приближение, запишите Eq. (1) как:

(−jωμ0∫s [e (s) g (Ра) +1k2dI (s) ds ∇→ g (Ра)]     ds+E→i), tan=0 (4)

где s - продольное положение вдоль провода, а e→ (ы) - единичный вектор, представляющий ориентацию провода в этом положении. Дискретизируйте провод в конечное число сегментов s (m), где m =1....N, а затем дополнительно дискретизировать Eq. (4) в конечное множество уравнений, наложенных в точках, соответствующих p =1....N для получения:

∑m=1N∫0hm[e→pe→mIm (sm) +e→pk2dIm (sm) dsm∇]g (ra)      dsm=e→pE→ijωμ0 (5)

где Im (sm) - распределение тока вдоль сегмента m, аппроксимируемое многочленом вида:

Im (sm) =∑i=0nmIm,i (sm hm2hm) i, 0≤sm≤hm, m =  1 , . . .., N (6)

где hm - длина сегмента и nm - выбранная степень полинома (которая отличается от сегмента к сегменту). Подставляя Eq. (6) в (5), набор уравнений может быть записан в матричном формате как:

[Zp, mi] [Imi] = [ Vex ,  p ] (7)

где [Zp, mi] - матрица импеданса, [Imi] - вектор неизвестных коэффициентов, представляющих ток на проводе, и [Vexp, p] - вектор возбуждения. Более подробная информация о расчете элементов [Zp, mi] приведена в [ 1].

Ссылки

[1] Попович, Б.С., М.Б. Драгович и А.Р. Джорджевич. Анализ и синтез исследований проводных антенн Пресса, 1982