Фрактальные антенны привлекают внимание благодаря своей уникальной геометрии, основанной на принципе самоподобия. Они позволяют объединить несколько ценных характеристик: уменьшить физические размеры устройства, сохранив при этом резонансную длину проводника, обеспечить работу в нескольких диапазонах и предоставить возможности для точной настройки диаграммы направленности. В статье подробно рассматривается процесс расчёта, моделирования, оптимизации и производства таких антенн, а также обсуждаются возможные компромиссы при их внедрении в реальные системы.
Математика и физика фракталов
Фрактал формируется путём многократного применения простого оператора итерации к исходной форме. На каждом этапе часть элемента заменяется набором более мелких элементов в соответствии с определённым правилом. В результате получается геометрическая структура, включающая элементы различных масштабов.
Фрактальная размерность
Для самоподобных фракталов размерность D можно вычислить по формуле
где N — количество копий, возникающих при масштабировании с множителем r$ а каждая копия имеет линейный масштаб r относительно исходного элемента.
Например:
Физическая интерпретация для антенн
Чем ближе фрактальная размерность к 2 (как у плоской шины), тем более плотно «заполняется» поверхность, что позволяет увеличить электрическую длину проводника при ограниченных размерах. Однако эффективность зависит не только от размерности $D$, но и от конкретной топологии, включая соединения, протяжённость контура и связность элементов.
Электрическая длина и распределение тока
Фрактальная топология создаёт множество путей для тока разной длины и ориентации. Резонансные характеристики фрагмента определяются его эффективной электрической длиной $L_{eff}$ и окружающей средой (подложка, соседние элементы). Наличие различных значений $L_{eff}$ приводит к появлению нескольких пиков в характеристике $S11$.
Популярные фрактальные формы
Различные фрактальные формы обладают уникальными характеристиками, преимуществами и недостатками, что определяет сферы их применения:
Кривая Коха** увеличивает периметр при ограниченной площади, что делает её подходящей для создания компактных узкополосных элементов. Однако рост периметра может привести к появлению острых углов и высокой плотности тока, что увеличивает потери. Применяется в встраиваемых компактных антеннах.
Треугольник Серпинского** благодаря удалению центральных частей создаёт множество масштабов, обеспечивая естественную многодиапазонность и простоту геометрии. Однако требует тщательной подгонки размеров и согласования. Используется в многодиапазонных планарных антеннах.
Кривые Хилберта и Пеано** обеспечивают плотное заполнение площади, максимально увеличивая длину проводника в заданных размерах. Однако их изготовление при высоком числе итераций может быть сложным, а потери — высокими. Применяются в случаях, когда критически важна миниатюризация.
Множественные «змеи» Минковски** с зигзагами и шагами на разных масштабах позволяют гибко управлять плотностью и обладают хорошей управляемостью. Однако требуют тщательного подбора параметров. Используются в портативных устройствах.
Каскадные и люнеточные структуры**, сочетающие фракталы с классическими элементами, позволяют достичь комбинированных свойств, но требуют более сложного анализа. Применяются в антеннах с заданной диаграммой направленности.
Миниатюризация и резонансные характеристики
Физическая резонансная длина проводника определяется его электрической, а не только геометрической длиной. Размещение сложного контура с большим периметром в ограниченной площади позволяет «удлинить» путь тока без увеличения занимаемого пространства, что способствует миниатюризации.
Оценочные правила
Для классического полуволнового диполя в свободном пространстве $L_{\lambda/2} \approx \lambda/2$. Например, для частоты 2,4 ГГц длина волны $\lambda \approx 125$ мм, а $\lambda/2 \approx 62,5$ мм. Применение фрактальных структур (одна-две итерации) часто позволяет сократить физические размеры антенны на 20–70 % по сравнению с базовой геометрией. Однако уменьшение габаритов может сопровождаться ростом реактивности, сужением рабочей полосы и увеличением потерь в проводнике или диэлектрике.
Расчёт эффективной длины
Эффективная длина $\ell_{eff}$ сложного пути приблизительно пропорциональна сумме эффективных сегментов, учитывая распределение тока и фазовые изменения из-за ветвления. Универсального аналитического выражения для расчёта не существует, поэтому применяются электромагнитные симуляции и экстраполяция на основе экспериментальных данных.
Многодиапазонность и широкополосность
Благодаря самоподобию фрактальные структуры содержат подэлементы разной длины $L_i$, каждый из которых может иметь свою резонансную частоту $f_i \approx \frac{v}{2L_i}$ для симметричных контуров или $f_i \approx \frac{v}{4L_i}$ для моноопорных/монопольных паттернов, где $v$ — скорость распространения волны в среде. При плотном распределении масштабов резонансы могут перекрываться, формируя более широкую рабочую полосу. Однако её ширина часто ограничена реактивностью и добротностью структуры.
Методика проектирования
Проектирование фрактальной антенны начинается с определения исходных требований:
— рабочие частоты или диапазоны;
— максимально допустимая площадь монтажа;
— требуемый коэффициент усиления, диаграмма направленности, поляризация;
— ограничения по материалам, толщине печатной платы, допустимым потерям;
— условия установки (встроенная или внешняя антенна) и влияние корпуса.
Затем выбирается базовая фрактальная форма с учётом поставленных задач:
— для многодиапазонной работы подойдёт треугольник Серпинского;
— для компактности и простоты — кривая Коха;
— для максимального увеличения длины проводника — кривые Хилберта или Пеано. Определение параметров и количество итераций
Рекомендации: целесообразно начинать процесс с одной-двух итераций, поскольку при большем их числе (более трёх) возрастает сложность и возможны потери. Коэффициент масштаба следует подбирать так, чтобы ключевой крупный элемент соответствовал необходимому уровню для резонанса на низкой частоте — как правило, это половина или четверть длины волны.
Шаг 3: выбор метода питания
В зависимости от выбранной схемы можно использовать различные подходы:
— при применении коаксиального вывода может потребоваться согласующий трансформатор;
— в случае микрополосковой подпитки необходимо учитывать параметры подложки, ширину линии и особенности перехода на контур;
— для достижения необходимого согласования можно применить конденсаторные или индуктивные схемы.
Шаг 4: первичное моделирование
Для оценки ключевых характеристик, таких как коэффициент отражения (S11), диаграммы направленности и распределение тока, следует воспользоваться симуляторами электромагнитных полей (EM-симуляторами) в двух- или трёхмерном режиме.
Шаг 5: оптимизация
Сначала проводятся параметрические расчёты, после чего осуществляется автоматическая оптимизация с применением различных методов, таких как генетические алгоритмы, метод роя частиц (PSO) или градиентные методы. В процессе оптимизации следует учитывать такие параметры, как масштабный коэффициент, угол ветвления, ширина линий, расстояние до заземляющего слоя и положение подпитки.
Шаг 6: проектирование с учётом требований производства
На этом этапе важно предусмотреть технологические допуски:
— установить минимально допустимую ширину проводящих дорожек;
— определить радиусы скругления на изгибах;
— предусмотреть защитную маску;
— организовать переходы для связи с нижним слоем (если это необходимо);
— усилить тонкие участки, которые могут быть подвержены механическим повреждениям.
Шаг 7: создание прототипа и проведение измерений
Для проверки коэффициента полезного действия необходимо провести калиброванные измерения в камере с реверберацией.
Практические рекомендации по моделированию
Выбор метода моделирования:
— для тонких проводников рекомендуется использовать метод моментов (MoM), например, в программах FEKO или NEC;
— для широких дорожек на подложке лучше подходят методы конечных элементов (FEM) или конечных разностей во временной области (FDTD), например, в HFSS, CST или openEMS.
Настройка сетки моделирования:
— обеспечить разрешение лучше λ/20 в областях с тонкими деталями;
— для острых углов и мелких ветвей применять локальную адаптивную сетку.
Пограничные условия:
— при моделировании антенн использовать PML или достаточно удалённую воздушную область.
Параметры материалов:
— задать диэлектрическую проницаемость (εr), тангенс угла потерь (tanδ), толщину медного слоя, электрическую проводимость (σ).
Моделирование питания:
— учитывать особенности реального разъёма, поскольку упрощённое представление порта может не выявить потенциальных проблем с согласованием.
Анализ результатов моделирования:
— изучить амплитуду и фазу коэффициента отражения (S11), добротность резонансов, коэффициент полезного действия, реальную излучаемую и затраченную мощность;
— оценить диаграммы направленности в двух- и трёхмерном представлении, поляризацию;
— проанализировать поверхностные плотности тока для выявления зон концентрации тока, что позволит определить места возможных нагревов и потерь.
Методы оптимизации антенн
Можно изменять масштаб отдельных элементов, угол разветвления, ширину полос, расположение точки питания. Для более сложных задач подойдут численные алгоритмы:
— генетические алгоритмы — эффективны для решения многомодальных задач, но требуют значительных вычислительных ресурсов;
— PSO / CPSO — часто обеспечивают более быструю оптимизацию при работе с задачами небольшой размерности;
— CMA-ES — хорошо подходит для работы с непрерывными параметрами.
Практические советы
— сократить пространство переменных: сначала провести глобальную оптимизацию (масштаб, количество итераций), затем перейти к локальной (ширина, радиус скругления);
— ускорить процесс, используя параллельные запуски симулятора на вычислительном кластере или в облачной среде;
— ограничить количество итераций фрактала в оптимизируемой области, поскольку мельчайшие детали могут существенно увеличивать время расчётов, при этом слабо влияя на конечный результат.
Производство, технологии и допуски
При разработке печатных плат следует учитывать:
— минимальную ширину проводящих дорожек (обычно 0,15–0,2 мм, в зависимости от возможностей производителя);
— минимальный зазор (около 0,15 мм);
— точность травления (обычно ±0,05 мм);
— толщину медного слоя (стандарт — 35 мкм (1 oz), для снижения потерь на высоких частотах можно использовать 2 oz);
— уязвимость тонких ветвей (при ширине менее 0,5 мм есть риск механических повреждений, поэтому стоит предусмотреть усиление или увеличить ширину дорожки).
Гибкие подложки и FPC
Они часто применяются в носимых устройствах, при этом важно контролировать, как изгиб влияет на резонансные характеристики.
Допуски на изготовление
Следует иметь в виду, что острые углы и узкие элементы будут скруглены — необходимо моделировать их с радиусом, соответствующим производственным допускам.
Лабораторные измерения и методики
Подготовка к измерениям:
— провести калибровку векторного анализатора цепей (VNA) (методы SOLT, TRL) вблизи точки подключения антенны;
— исключить риск плохого контакта в точке подключения (балун, переход).
Параметры для измерения:
— коэффициент отражения (S11 или S22 для передатчика) — определить частотные пики и провалы, ширину полосы на уровне −10 дБ;
— коэффициент стоячей волны (КСВ), который можно вывести из данных S11;
— коэффициент полезного действия — можно оценить методом сравнения в безэховой и реверберационной камерах или использовать метод Уилера для быстрых оценок;
— диаграмму направленности и усиление — измерить в безэховой камере, сравнивая с эталонной антенной;
— плотность тока — получить косвенно, анализируя распределение поля в камере или с помощью моделирования.
Ограничения и компромиссы
— Миниатюризация может привести к увеличению сопротивления излучающей дорожки и потерям в диэлектрике, что снижает КПД.
— Расширение функционала за счёт аккумуляции резонансов может обеспечить многодиапазонность, но при этом отдельные частотные полосы могут оказаться слишком узкими.
— Мелкие геометрические детали чувствительны к точности травления и расположению соседних модулей.
— Точный аналитический прогноз затруднителен — требуется электромагнитное моделирование и создание прототипа.
— Производственные ограничения могут не позволить точно воспроизвести тонкие ветви и острые углы. В процессе последовательных итераций значение Pn определяется как произведение предыдущего значения Pn−1 на коэффициент 4/3. Таким образом, получаем следующие результаты: P1 составляет 4 сантиметра, P2 — приблизительно 5,333 сантиметра (16/3), а P3 — примерно 7,111 сантиметра (64/9).
Что касается площади, то начальная площадь A0 рассчитывается по формуле и составляет примерно 0,4330127 квадратных сантиметра. Общая формула для вычисления An выглядит следующим образом: An = A0 + A0 ⋅ 1/3 ∑(k=0 до n−1) (4/9)^k = A0 (1 + 1/3 ⋅ (1 − (4/9)^n) / (1 − 4/9)). В предельном случае площадь A∞ стремится к значению примерно 0,692820323 квадратных сантиметра.
С точки зрения практического применения в антенной технике, если размер стороны s сопоставим с габаритами платы, то увеличение периметра может привести к возрастанию эффективного длины ℓeff. Это, в свою очередь, вызовет смещение резонансной частоты в нижнюю область спектра при сохранении исходной площади.
Рассмотрим пример проекта на основе фрактальной структуры Серпинского для работы в диапазонах 2,4–5 ГГц. Цель проекта — создать планарную фрактальную антенну, способную охватить оба указанных диапазона.
Исходя из того, что половина длины волны для диапазона 2,4 ГГц составляет примерно 62,5 мм (λ/2 = c / (2f)), мы принимаем начальную длину стороны s равной 62,5 мм. Это значение соответствует половине длины волны для полуволновой версии антенны.
Проанализируем результаты последовательных итераций:
— после первой итерации длина стороны s1 уменьшается до 31,25 мм (s/2), что приводит к резонансу примерно в диапазоне 4,8–5 ГГц;
— после второй итерации длина стороны s2 сокращается до 15,625 мм (s/4), а резонансная частота смещается в область 9–10 ГГц. Хотя этот диапазон выходит за пределы целевого, он может оказывать влияние на гармоники.
Важно отметить, что фактические резонансные частоты могут варьироваться в зависимости от конструктивных особенностей антенны, используемых материалов и наличия заземляющего слоя. Для точного определения параметров необходимо проводить электромагнитное моделирование и осуществлять экспериментальную подгонку.
Среди литературных источников, которые могут быть полезны при изучении фрактальных структур и их применения в антенной технике, можно выделить следующие работы:
— Мандельброт Б. Б. «Фрактальная геометрия природы»;
— Фальконер К. «Фрактальная геометрия. Математические основы и приложения»;
— Пуэнте С., Ромеу Х. «Фрактальные многодиапазонные антенны на основе треугольника Серпинского» (публикации в IEEE);
— Бест С. Р. «Свойства фрактальных антенн» (IEEE Antennas and Propagation Magazine);
— Ангоэра Ж., Пуэнте С. «Фрактальные антенны: обзор» (Wiley);
— статьи из журнала «Антенны»;
— учебные и обзорные материалы Горячева А. В. и Киселёва В. В. «Компактные и многодиапазонные антенны для беспроводных устройств»;
— сборники трудов конференций РЭНС, МРК, РЭК, содержащие доклады о фрактальных и миниатюрных антеннах.