Моделирование распределения плотности тока в сложном неоднородном проводнике. Часть 1 Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Гришенцев А.Ю.

В статье рассмотрена математическая основа моделирования распределения плотности тока в сложном, неоднородном проводнике, имеющем произвольную форму поперечного сечения.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Гришенцев А.Ю.

Текст научной работы на тему «Моделирование распределения плотности тока в сложном неоднородном проводнике. Часть 1»

МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ ПЛОТНОСТИ ТОКА В СЛОЖНОМ НЕОДНОРОДНОМ ПРОВОДНИКЕ.

Научный руководитель - доктор технических наук, профессор К.Г. Коротков

В статье рассмотрена математическая основа моделирования распределения плотности тока в сложном, неоднородном проводнике, имеющем произвольную форму поперечного сечения.

Получивший в последнее время широкое распространение во всем мире процесс газоразрядной визуализации (ГРВ) [1], непосредственно связан с протеканием токов высокой частоты через исследуемые объекты. При этом наблюдается неравномерное распределение плотности тока по сечению, называемое поверхностным эффектом. Моделирование поверхностного эффекта в однородных проводниках подробно рассмотрено в различной литературе по электротехнике. В данной статье построена математическая модель поверхностного эффекта для случая неоднородных проводников. В этом случае вклад в неравномерность распределения тока по сечению вносит неоднородность среды, например, биологических объектов, исследуемых на ГРВ оборудовании.

Моделирование электропроводящих свойств, сложных биологических тканей, рассмотрено в ряде современной литературы (например: [2-4] и др.). К сожалению не одна из предложенных моделей не позволяет оценить плотность распределения тока по биологическим тканям в процессе ГРВ, ИПЧ исследований.

Автором статьи предложен возможный вариант оценки плотности распределения тока, по биологическим тканям применимый к ГРВ и ИПЧ измерениям. Моделирование производилось на основе законов классической электродинамики ([5-9]), с учетом особенностей протекания токов по биологическим тканям ([3, 4, 10-14]).

Модель плотности распределения высокочастотного тока по сечению во время исследования на ГРВ оборудовании позволяет:

• рассмотреть различия распределения плотности тока в зависимости от электрических параметров исследуемого объекта;

• определить части объекта, по которым протекает превалирующая часть тока;

• предсказать тенденцию изменения картины ГРВ-граммы при изменении электрических параметров исследуемого объекта;

• предсказать тенденцию изменения интегральной суммы тока для приборов-измерителей, интегральной суммы токов поверхностных частот (ИПЧ), при изменении электрических параметров исследуемого объекта;

• с высокой точностью оценить значения протекающих по различным тканям токов (при исследовании человека) с целью выяснения безопасности обследования на ГРВ, ИПЧ оборудовании.

Построение математической модели

Рассмотрим картину магнитного поля в плоскости поперечного сечения проводника с током, выполненного из однородного материала (рис. 1). Представим этот провод в виде совокупности нитей, параллельных его оси. Чем ближе нить расположена к оси проводника, тем с большим числом магнитных линий она сцеплена.

Рис. 1. Магнитное поле в плоскости поперечного сечения проводника с током

При периодическом изменении тока изменяется магнитное поле, и в нитях проводника наводится э.д.с., противодействующие изменениям тока. Это противодействие тем значительнее, чем больше э.д.с. (чем больше магнитных линий сцеплено с нитью), т.е. чем ближе нить проводника расположена к оси провода. В результате плотность тока в различных точках поперечного сечения получается не одинаковой: наибольшая на периферии и наименьшая на его оси.

Рассмотренное явление концентрации переменного тока в поверхностном слое проводника называют поверхностным эффектом. Резкость проявления его возрастает с увеличением частоты /, диаметра проводника ё, магнитной проницаемости ц и удельной проводимости у материала проводника. Следствием поверхностного эффекта является некоторое уменьшение собственной индуктивности проводника ввиду ослабления магнитного поля во внутренней части проводника. В предельном случае, когда ток концентрируется на поверхности проводника в бесконечно тонком слое, магнитное поле внутри проводника отсутствует.

Другим следствием поверхностного эффекта является возрастание тепловых потерь при одинаковых значениях переменного и постоянного тока (равенстве значения постоянного тока и действующего значения переменного тока): тепловые потери больше при переменном токе. Поэтому сопротивление проводника переменному току (активное сопротивление) выше, чем сопротивление проводника постоянному току.

Степень неравномерности распределения тока по сечению проводника благодаря поверхностному эффекту в каждом конкретном случае можно оценить, сравнивая линейные размеры поперечного сечения провода с длиной волны электромагнитных колебаний в проводнике X или с эквивалентной глубиной проникновения электромагнитной волны ([8], стр. 6):

где ю - круговая частота, ц - абсолютная магнитная проницаемость вещества, у -удельная объемная электропроводность вещества.

На распределение переменного тока в проводе также оказывают влияние токи соседних проводников. Это явление называют эффектом близости ([6, 8]).

Следует отметить, что в поверхностном слое плотность тока не только убывает по величине, но и изменяется по фазе. На поверхности проводника плотность тока опережает по фазе полный ток на 45°: плотность тока на поверхности пропорциональна комплексной напряженности электрического поля Е, полный ток пропорционален комплексной напряженности магнитного поля Н, Е опережает Н по фазе на 45°. По мере перехода к более глубоким слоям плотность тока запаздывает относительно поверхностного тока, и на глубине, равной половине длины волны, ток изменяет направление на противоположное. Однако этот обратный ток не играет существенной роли - в поверхностном слое толщиной в половину длины волны Х/2 протекает свыше 95% общего тока [6].

Если проводники имеют сложную форму поперечного сечения, то из-за эффектов поверхностного и близости по отдельным частям проводников протекают токи различной плотности.

Допустим, что имеется проводник длиной I, произвольного поперечного сечения рис. 2 с током:

где 1т - амплитуда тока [А]; ю = 2п/- угловая частота тока [рад/сек]; /- частота тока [Гц]; I - текущее время [с]. Разобьем поперечное сечение проводника на элементарные площадки. В проводнике площадки могут отличаться по размерам одна от другой. Однако, для удобства расчетов, желательно иметь одинаковые по размерам площадки. Диаметр площадки рекомендуется выбирать примерно вдвое меньше, чем глубина проникновения ё. Будем рассматривать образованные разбиением проводники как элементарные со своим током ¡к, где принимает значения 0-„. Заменим получившийся контур эквивалентной электрической схемой рис. 3, состоящей из („+1) ветвей. Если найти токи в каждой ветви, а затем плотность тока в каждом элементарном проводнике, то задача по определению поля плотности тока будет решена.

Рис. 2. Проводник произвольного поперечного сечения

Рис. 3. Эквивалентная схема

На эквивалентной схеме (рис. 3) обозначены: I - мгновенное значение суммарного тока [А]; ¡0, ¡1. 1„ - мгновенные значения токов в элементарных проводниках [А]; Ь0,0, Ьц. Ь„,„ - собственные индуктивности элементарных проводников [Гн]; Я0, Я],. Я„ -активные сопротивления элементарных проводников [Ом].

Для расчета собственных индуктивностей ветвей с 0-й по п-ю рассмотрим произвольный к-й элемент (рис. 4). Координаты центра сечения элементарного проводника хк, ук; размеры сечения Ьк, ск [м]. При моделировании удобно проводить разбиение на элементарные проводники квадратного сечения, поэтому Ьк=ск. Если размеры сечения элементарного проводника составляют половину или меньше от глубины проникновения 3, формула (1) (случай низкой частоты), можно считать плотность тока в пределах сечения единичного проводника постоянной. В этом случае индуктивность и активное сопротивление единичного проводника можно находить, как при постоянном токе. В случае высокой частоты, когда размеры сечения больше половины глубины проникновения, надо учитывать неравномерность распределения тока по сечению. В дальнейшей реализации модели есть возможность разбиения проводника на равные элементарные площадки квадратного (со стороной ребра Ьк) и кругового (диаметр ёк) сечений, поэтому будем рассматривать квадратное и круговое сечения. В приведенных ниже формулах: ук - удельная проводимость, 5к - глубина проникновения для к-го элементарного проводника.

Рис. 4. Произвольный элемент проводника

В случае низкой частоты имеем:

для проводников квадратного сечения собственную индуктивность будем вычислять по следующей формуле (Ьк< 25) ([8], с. 101):