Использование метода конечных элементов для расчёта парциального распределения ёмкости микрополосковой линии Текст научной статьи по специальности «Физика»

УДК 519.632.4:621.3.01

М.Г. Рубанович, Д.В. Вагин, В.А. Хрусталев, П.Г. Богомолов, А.А. Столяренко

Использование метода конечных элементов для расчёта парциального распределения ёмкости микрополосковой линии

Рассмотрены два метода расчёта парциального распределения ёмкости в поперечном сечении микрополосковой линии: 1) непосредственное решение двумерной краевой задачи методом конечных элементов; 2) предварительное применение конформного отображения расчётной области. В обоих случаях для нахождения потенциала электрического поля решается эллиптическая краевая задача с оператором Лапласа. В первом случае используются прямоугольные конечные элементы с билинейными базисными функциями. Во втором случае из-за наличия в расчётной области сложной границы между воздухом и диэлектриком используются треугольные конечные элементы с линейными базисными функциями. По полученному распределению потенциала вычисляются значения напряженности электрического поля вблизи поверхности микрополосковой линии. Далее по теореме Гаусса определяется распределение заряда по поверхности микрополосковой линии. Разработанные алгоритмы позволяют определять парциальную ёмкость на любом дискретном интервале поперечного сечения несимметричной полосковой линии, расположенной на диэлектрической подложке. Для второго рассмотренного метода сделана оценка погрешности расчёта суммарной ёмкости. Такая же оценка сделана для расчёта суммарной ёмкости по предложенным авторами формулам. Ключевые слова: микрополосковая линия, ёмкость, двухмерная задача электростатики, эффективная диэлектрическая проницаемость, волновое сопротивление, краевая задача, метод конечных элементов, конформно отображённая область.

При проектировании радиотехнических устройств необходимо вычислять реактивные параметры микрополосковой линии передачи. Это позволяет рассчитывать её рабочие параметры. Для расчёта ёмкости на единицу длины микрополосковой линии существует ряд методов [1, 2]. Также существует ряд других специальных методов расчёта ёмкости [3]. Целью данной работы является качественная оценка метода расчёта парциального распределения ёмкости на основе конформного отображения на примере задачи расчёта ёмкости микрополосковой линии.

Для проверки полученных результатов в работе использован метод конечных элементов, являющийся методом точного численного моделирования [4]. Также приводится сравнение с результатами расчёта суммарной ёмкости для микрополосковой линии по формулам, предложенным авторами.

Математическая модель. Вычислить значение емкости металлического полоска можно, определив распределение потенциала электрического поля и(х,у) , такого, что Е = -gradu в области О

(ввиду симметрии расчетной области рассматривается только ее половина), схематично изображенной на рис. 1. Граница 55 является плоскостью симметрии, поэтому на ней должно быть задано однородное второе краевое условие на потенциал и . Границы 53 и 54 предполагаются достаточно удаленными, поэтому для них также будет справедливым однородное второе краевое условие для потенциала и . Граница является поверхностью металлического полоска.

Соответствующая дифференциально-краевая задача будет иметь вид

Для решения этой краевой задачи воспользуемся методом конечных элементов (МКЭ) с прямоугольными конечными элементами [4]. Выпишем для уравнения (1) с краевыми условиями (2) вариационную постановку в форме Галеркина. Для этого умножим левую и правую части уравнения (1) на пробные функции

-div(sgradM) = 0,

(1)

(2)

| (- d/v(sgradu)) и'о <0 = 0, и е Н^ м>о е Н0

(3)

О

где Н^т - это пространство функций, имеющих суммируемые с квадратом производные и удовле-творяющих заданным в (2) краевым условиям первого рода, а Но - пространство функций, имеющих суммируемые с квадратом производные и равные нулю на границах 51 и 52 • Уравнение (3) можно переписать в виде

. ЯЬ/ 1 1

(4)

Jsgradu• gradwоdО+|е—м^dS = 0, и еН^ м'о еНО .

О

5,

^=о дп

ди=о

дп

5 и = о

Рис. 1. Исходная расчётная область

где 5 = У . Если выбрать достаточно уда -

I=1

ленные границы 53 и 54 , то на них будет

ди

справедливо краевое условие — = о, таким

дп

образом, окончательно получаем уравнение Jsgradu • gradw0dО = о, и е Н^, мо е Но. (5) О

Будем искать приближенное решение и^ потенциала и в виде линейной комбинации

финитных базисных функций уj с весами qj . п

и = Е qj уj, j=1

(6)

где п - число узлов конечноэлементной сетки.

Выбирая в качестве пробных функций мо функции у,, / й , где - набор индексов, соответствующих узлам сетки, в которых заданы краевые условия первого рода, и подставляя (6) в (5),

получаем систему линейных алгебраических уравнений (СЛАУ) вида:

( \

тЪ<

j=1\О

Jsgradуj • gradу^О

qj = 0, , й №

/

Ъ<0

qj = 0, j е N qj = 1, ]■ е

(7)

где N 0 и N 1 - это наборы индексов, соответствующих узлам сетки, расположенным на границах 51 и 52 соответственно, NЪd = ^^ УNЪdl .

Парциальное распределение емкости или суммарная ёмкость С могут быть вычислены по полученному в результате решения задачи (1), (2) распределению электрического поля. Например, для определения суммарной ёмкости МПЛ надо проинтегрировать нормальную составляющую вектора напряженности электрического поля по поверхности полоска (по границе 52):

С = \&Еп<5 . (8)

52

Математическая модель для решения задачи в конформно отображенной области. После конформного отображения расчетная область примет вид, схематично изображенный на рис. 2, где 5>1 - 5>4 - границы конформно отображенной расчетной области.

В этом случае краевые условия для потенциала электрического поля и = и(х,у) на соответствующих границах такие же, как и в исходной расчетной области. Для данной области удобнее использовать треугольные конечные элементы. В остальном метод численного расчета распределения потенциала и полностью совпадает с описанным выше алгоритмом расчёта в исходной области (рис. 1). Однако для вычисления емкости по описанному в работе методу интегрировать нормаль-

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

Результаты численного моделирования. Ниже приведены результаты расчетов для двух моделей. Расчеты выполнялись на нескольких вложенных сетках. Погрешность результатов вычислений не превышает 1%.

В первой модели ширина полоска равна 1 мм, толщина полоска - 6 мкм, толщина диэлектрика - 1 мм, а значение диэлектрической проницаемости диэлектрика е = 9,6ео. Во второй модели ширина полоска равна 6 мм, толщина полоска - 6 мкм, толщина диэлектрика - 4 мм, а значение диэлектрической проницаемости диэлектрика е = 6,6е0.

На рис. 3 показано распределение удельно -

го значения заряда Q , Кл/(В-м ) по нижней

Рис. 2. Часть конформно отображенной расчетной области вблизи раздела воздуха и диэлектрика

(со стороны диэлектрика) и по верхней (со стороны воздуха) поверхностям металлического полоска для обеих моделей.

На рис. 4 показано распределение удельного значения заряда Q' 52 в конформно отображенной расчетной области.

Кл/(В-м2) на границах 51 и

Кл

10-

10-

10-

-10-

-10-

Часть проводящей плоскости, расположенной непосредственно^ под полоском

3 2,5 2 1,5 1 0,5 Рис. 3. Распределение удельного значения заряда по поверхности полоска и на медном основании под полоском

4 3 2 1 0 -1

-6 -4 -2 0 2 4 Рис. 4. Распределение удельного значения

заряда на границах ¿1 и 52

Сравнение точности значений суммарной ёмкости, полученных методом конечных элементов, с соответствующими ёмкостями, полученными по предложенным формулам. Значение еэфф определяется отношением емкостей Сд к Сбд [5]:

еэфф :

:Сд

Сбд

(9)

где Сбд - погонная ёмкость МПЛ без диэлектрика; Сд - погонная ёмкость МПЛ с диэлектриком. Расчётное еэфф определяется по формуле [6]:

0

X , мм

ег +1 ег-1 10-Я. 2

Зэфф =——+—— (1 +-) 2, (10)

эфф 2 2 ' ^фф

где Я - толщина диэлектрика; ^фф - эффективная ширина МПЛ. Эффективная ширина МПЛ Wэфф по формуле [5]:

Яэфф = Я + (П)

1 + 1п(2Я)

8

(11)

где 8 - толщина полоска.

Яэфф не рассчитывалась, так как при 8 = (5-6) мкм Я отличается от Яэфф в пятом знаке.

Я

Точность расчёта по формуле (10) ±2% [2] в широком диапазоне ег и — .

Я

Волновое сопротивление МПЛ с диэлектриком между полосковой линией и основанием определялось по формуле [5]

2в = 7—I-=—(Т^' (12)

Сбд 'с ^еэфф с \ Сбд ' Сд

где с - скорость света.

При этом справедливы следующие соотношения (13)—(15) [5]:

2в (13)

у8эфф

Я 8Я Я

при —<1, 20 = 601п(-+-), (14)

Я Я 4Я

Я ъ1 2 _120п__(15)

при — > 1, 20 = -ф-я-Я7, (15)

Я ^ + 2,42 - 0,44Я+(1- Я )6

Я Я Я

где 20 - волновое сопротивление МПЛ с воздухом между полосковой линией и основанием.

Приведённые выше выражения (9) и (12) являются системой уравнений относительно Сбд и Сд. Если разрешить систему уравнений относительно Сбд и Сд, то имеем [5]

Сбд = -1-, (16)

с - 20

Сд = ^ . (17)

с - ¿0

Рассчитаем емкостные параметры для микрополосковой линии передачи, выполненной на основе брокерита с относительной диэлектрической проницаемостью ег = 6,6 размером 6x4 мм2 (6 мм шириной, 4 мм толщиной) и для микрополосковой линии передачи, выполненной на основе поликора с относительной диэлектрической проницаемостью ег = 9,6 размером 1x1 мм2 (1 мм шириной, 1 мм толщиной), полученных конформным отображением и МКЭ. Определим суммарную ёмкость при использовании метода конформного отображения (МКО) для микрополосковой линии (МПЛ) размером 6x4 мм2. Суммарная ёмкость без диэлектрика равна Сбд = 30,44 пФ/м. При расчёте суммарной ёмкости МКЭ для МПЛ размером 6x4 мм2 без диэлектрика ёмкость равна Сбд = 32,08 пФ/м. Разница при расчёте суммарной ёмкости составляет 1,64 пФ/м. Суммарная ёмкость при использовании метода конформного отображения для МПЛ размером 1x1 мм2 без диэлектрика Сбд = 24,5 пФ/м. При

использовании метода конечных элементов для МПЛ 1x1 мм2 без диэлектрика Сбд = 26,43 пФ/м.

Разница при расчёте суммарной ёмкости составляет 1,93 пФ/м.

Причина такой разницы в том, что расчёт МКО не учитывает взаимодействие зарядов на краях МПЛ, так как рассчитывается распределение зарядов на одном торце и прилегающих поверхностях сверху и снизу МПЛ, в то время как при расчёте МКЭ взаимодействие зарядов на краях учитывается

автоматически. Таким образом, непосредственный расчёт МКЭ (см. рис. 1) дал более точное распределение емкости и ее суммарное значение.

Рассчитаем ёмкость в поперечном сечении МПЛ, если МПЛ расположена на диэлектрике, с относительной диэлектрической проницаемостью ег = 6,6 для МПЛ размером 6x4 мм2. Расчёт проводился конформным отображением одного торца с прилегающими отрезками верха и низа МПЛ, отображёнными в прямую линию (см. рис. 2) с последующим расчётом образовавшейся области методом конечных элементов (МКО с МКЭ). Кривая граница - это граница между воздухом и диэлектриком. Также эта задача (см. рис. 1) рассчитывалась непосредственно методом конечных элементов. Суммарная ёмкость с диэлектриком при использовании МКО с МКЭ Сд = 151,63 пФ/м. При

использовании МКЭ Сд = 151,24 пФ/м для МПЛ размером 6x4 мм2 с диэлектриком под МПЛ с

ег = 6,6 . Разница при расчёте суммарной ёмкости составляет 0,26%.

При использовании МКО с МКЭ для МПЛ размером 1x1 мм2 с диэлектриком под МПЛ с ег = 9,6 суммарная ёмкость равна Сд = 169,55 пФ/м. При использовании МКЭ суммарная ёмкость

равна Сд = 170,07 пФ/м. Разница полученных значений составляет 0,31%. Столь малое отличие при

расчёте суммарной ёмкости методом МКО с МКЭ и методом МКЭ объясняется тем, что увеличению ёмкости на краях поперечного сечения МПЛ при расчёте МКЭ соответствует примерно равное уменьшение ёмкости на его торцах по сравнению с расчётом МКО с МКЭ.

Для оценки точности предложенных формул (16) и (17) были рассчитаны значения Сбд и Сд для МПЛ с относительной диэлектрической проницаемостью ег = 6,6 размером 6x4 мм2 и для МПЛ с относительной диэлектрической проницаемостью ег = 9,6 размером 1x1 мм2, которые были сопоставлены с аналогичными значениями, полученными МКЭ. При размере МПЛ 6x4 мм2, ег = 6,6, эффективная диэлектрическая проницаемость, рассчитанная по формуле (10), еэфф = 4,81. Определённая МКЭ суммарная погонная ёмкость для МПЛ с воздухом между полосковой линией и основанием Сбд1 = 32,077пФ/м. Рассчитав по формуле (15) волновое сопротивление для МПЛ с воздухом между полосковой линией и основанием, получим 20 =103,936 Ом. Определённая по формуле (16) суммарная погонная ёмкость МПЛ с воздухом между полосковой линией и основанием Сбд2 =32,072 пФ/м. Относительная разница между Сбд1 и Сбд2 равняется 0,02%. Результаты расчёта по МКЭ и по выражению (16) мало отличаются. По формуле (17) определим суммарную погонную ёмкость МПЛ с диэлектриком между полосковой линией и основанием. Ёмкость Сд2 =

= 154,3 пФ/м. Определённая МКЭ суммарная погонная ёмкость для МПЛ с диэлектриком между полосковой линией и основанием Сд1 = 151,24 пФ/м. Относительная разница между Сд2 и Сд1 равняется 2%.

При размере МПЛ 1x1 мм2, е г = 9,6 эффективная диэлектрическая проницаемость, рассчитанная по формуле (10), еэфф = 6,596. Определённая МКЭ суммарная погонная ёмкость для МПЛ с воздухом между полосковой линией и основанием Сд1 = 26,43 пФ/м. Рассчитав по формуле (14)

волновое сопротивление для МПЛ с воздухом между полосковой линией и основанием, получим 20 =126,61 Ом. Определённая по формуле (16) суммарная погонная ёмкость МПЛ с воздухом между полосковой линией и основанием Сбд2 = 26,32 пФ/м. Относительная разница между Сд1 и Сбд2

равняется 0,4%. Результаты расчёта по МКЭ и по выражению (16) тоже мало отличаются. По формуле (17) определим суммарную погонную ёмкость МПЛ с диэлектриком между полосковой линией и основанием. Ёмкость Сд2 = 173,6 пФ/м. Определённая МКЭ суммарная погонная ёмкость для

МПЛ с диэлектриком между полосковой линией и основанием Сд1 = 170,08 пФ/м. Относительная разница между Сд2 и Сд1 равняется 2%.

Меньшая относительная погрешность определения ёмкости Сбд объясняется тем, что в формуле (16) в знаменателе стоит 20, относительная погрешность определения которой около 1% [5]. Формула (17), определяющая Сд в числителе, имеет еэфф, относительная погрешность которой [5]

еэфф

2%. Относительная погрешность частного от деления —равняется сумме относительных по-

¿0

грешностей еэфф и ¿о [7]. Таким образом, относительная погрешность при расчете Сд равняется 3%.

Выводы. Расчёт ёмкости методом МКО с МКЭ не учитывает взаимодействие заряда на краях микрополосковой линии. Поэтому если необходимо получать значение ёмкости с погрешностью, не превышающей 3-5%, нужно решать задачу напрямую без конформного отображения. При более мягких требованиях к ошибке определения суммарной ёмкости удобно использовать предложенные формулы (16)-(17). Для качественной оценки полученных значений расчёта суммарных ёмкостей Сбд и Сд для микрополосковой линии передачи, выполненной на основе брокерита с относительной диэлектрической проницаемостью ег = 6,6 размером 6x4 мм2, и для микрополосковой линии передачи, выполненной на основе поликора с относительной диэлектрической проницаемостью ег = 9,6 размером 1x1 мм2, были выполнены расчёты методом конечных элементов без конформного отображения.

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

Работа выполнена при финансовой поддержке Министерства образования и науки РФ в рамках базовой части государственного задания, тема проекта: Теория построения и практическая реализация широкополосных микроволновых устройств, предназначенных для использования в измерительном оборудовании цифровых систем телевидения и телекоммуникаций, включая вычислительные системы.

Литература

1. Ковалев И.С. Теория и расчёт полосковых волноводов / И.С. Ковалев. - Минск: Наука и техника, 1967. - 233 с.

2. Кухаркин Е.С. Машинные методы расчета в инженерной электрофизике / Е.С. Кухаркин, Б.В. Сестрорецкий. - М.: Моск. энерг. ин- т, 1986. - 68 с.

3. Иоссель Ю.Я. Расчет электрической емкости / Ю.Я. Иоссель, Э.С. Кочанов, М.Г. Струнский. -Ленингр. отд.: Энергия, 1969. - 240 с.

4. Соловейчик Ю. Г. Метод конечных элементов для решения скалярных и векторных задач / Ю.Г. Соловейчик, М.Э. Рояк, М.Г. Персова. - Новосибирск: НГТУ, 2007. - 895 с.

5. Красов В.Г. Толстоплёночная технология в СВЧ-микроэлектронике / В.Г. Красов, Г.Б. Петрау-скас, Ю.С Чернозубов. - М.: Радио и связь, 1985. - 168 с.

6. Справочник по расчёту и конструированию СВЧ-полосковых устройств / С.И. Бахарев, В.И. Вольман, Ю.Н. Либ и др.; под ред. В.И. Вольмана. - М.: Радио и связь, 1982. - 328 с.

7. Поршнев С.В. Вычислительная математика: курс лекций. - СПб.: БХВ-Петербург, 2004. - 302 с.

Рубанович Михаил Григорьевич

Канд. техн. наук, доцент каф. общей физики Новосибирского гос. технического университета (НГТУ)

Тел.: 8 (383) 346-08-34

Эл. почта: mihail-nstu@mail.ru

Вагин Денис Владимирович

Канд. техн. наук, доцент каф. прикладной математики НГТУ Тел.: 8 (383) 346-27-76 Эл. почта: vdv_wk@mail.ru

Хрусталев Владимир Александрович

Д-р техн. наук, профессор каф. электронных приборов НГТУ

Тел.: 8 (383) 346-06-64

Эл. почта: VA_Khructalev@ngs.ru

Богомолов Павел Геннадьевич

Аспирант каф. общей физики НГТУ Тел.: 8 (383) 210-02-81 Эл. почта: rk_6-41@mail.ru

Столяренко Алексей Андреевич

Магистрант каф. радиоприёмных и радиопередающих устройств НГТУ

Тел.: 8 (383) 346-08-34

Эл. почта: mihail-nstu@mail.ru

Rubanovich M.G., Vagin D.V, Khrustalyov V. A., Bogomolov P.G., Stolyarenko A.A.

Using of the finite element method for calculating the partial capacity distribution in the microstrip line

The paper considers two methods for calculating the partial capacity distribution in the cross section of the microstrip line: 1) direct solution of two-dimensional boundary-value problem by the finite element method; 2) calculation with preliminary conformal mapping of the computational domain. In both cases, elliptic boundary-value problem with Laplace operator is solved to find the electric field potential. In the first case, rectangular finite elements with bilinear basis functions are used. In the second case, triangular finite elements with linear basis functions are used because of the complex boundary between air and dielectric in the computational domain. Upon receipt of the potential distribution, values of the electric field near the surface of the microstrip line are calculated. Then we define the charge distribution on the microstrip line surface by means of Gauss theorem. The developed algorithms are used to determine the partial capacity on any discrete interval of the asymmetric stripline cross-section disposed on a dielectric substrate. For the second method we calculated error estimation of the total capacity. The same estimation is made for the total capacity calculation by the formulas, proposed by the authors.

Keywords: microstrip, capacity, two dimensional electrostatic problem, effective permittivity, wave resistance, boundary problem, finite element method, conformal area mapping.