Радиотехника и электроника, 2022, T. 67, № 1, стр. 84-90

Аналитические условия для выравнивания и уменьшения амплитуд составляющих временного отклика в витке меандровой линии

Р. С. Суровцев *

Томский государственный университет систем управления и радиоэлектроники
634050 Томск, просп. Ленина, 40, Российская Федерация

* E-mail: surovtsevrs@gmail.com

Поступила в редакцию 30.04.2021
После доработки 23.08.2021
Принята к публикации 02.09.2021

Полный текст (PDF)

Аннотация

Отмечена актуальность защиты радиоэлектронной аппаратуры от кондуктивных воздействий, в частности сверхкоротких импульсов (СКИ). Выделен виток меандровой линии (МЛ), позволяющий улучшить защиту от СКИ за счет его разложения на последовательность импульсов меньшей амплитуды. Отмечено, что для быстрой априорной оценки амплитуды СКИ на выходе МЛ в ряде случаев целесообразно применение простых аналитических моделей, а не затратного полноволного анализа. Описан процесс получения аналитических моделей в замкнутой форме для расчета временного отклика в узлах витка МЛ. На основе разработанных моделей получены условия равенства амплитуд импульсов разложения в конце витка МЛ для трех случаев разложения.

ВВЕДЕНИЕ

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

Типичным примером такого элемента является меандровая линия (МЛ), традиционно применяемая для задержки сигналов, поступающих из разных точек печатных плат для их синхронизации в точке приема. Также МЛ нашли применение в приемных и передающих антеннах [1], индукторах [2], многослойных конденсаторах [3], фильтрующих устройствах [4], корректорах групповой задержки и фазовых корректорах [5].

Между тем широко применяемые МЛ, если сделать их в специальной конфигурации, можно использовать также для подавления опасных кондуктивных воздействий, в частности сверхкороткого импульса (СКИ). Так, СКИ, распространяющийся вдоль витка МЛ, может быть разложен на несколько основных импульсов меньшей и одинаковой амплитуды: два импульса в воздушной МЛ [6], три – в микрополосковой [7] и четыре – в несимметричной, с лицевой связью [8]. Действительно, необходимость защиты от такого воздействия отмечалась в указанных работах [6–8], подтвердилась опубликованными фактами различных применений [9], показана, например, для авионики [10], став особенно насущной для беспилотных летательных аппаратов [11]. Специфика результата воздействия именно СКИ состоит в том, что наводки от него могут восприниматься в качестве полезных сигналов, однако они нарушают цифровой обмен, а при более высокой амплитуде проникают через традиционные средства защиты и приводят к выходу из строя [12, 13]. Уменьшение результирующей амплитуды помехи до безопасного уровня с помощью новых (одиночных или последовательно соединенных) устройств, используемых вместо известных или дополнительно с ними, может улучшить защиту. Поэтому активно исследуются новые устройства, позволяющие ослабить влияние СКИ за счет его разложения на импульсы с меньшими амплитудами (минимальными при их выравнивании), например, модальные фильтры [14]. Между тем очень близки к ним защитные МЛ [15, 16], которые могут даже превосходить их: за счет большего числа импульсов разложения, удвоенного пути вдоль МЛ и отсутствия резисторов [68]. Таким образом, исследование и применение МЛ для защиты от СКИ видится перспективным.

Для моделирования полосковых устройств, в том числе МЛ, как правило, используются численные методы. Например, для исследования минимизации искажений полезных сигналов, возникающих в МЛ на печатных платах из-за перекрестных связей, успешно применяется SPICE-моделирование [17], а также широко известные полноволновые методы [18]. Однако использование численных методов не всегда целесообразно, поскольку в ряде частных случаев (с незначительными потерями и дисперсией) вычислительные затраты на моделирование численным методом могут быть крайне высоки. Типичным примером является оптимизация, выполняемая для выравнивания и минимизации амплитуд разложенных импульсов [19]. Тогда для предварительных оценок целесообразно использовать аналитические модели, которые дают быстрые оценки с приемлемой точностью. Кроме того, аналитическая модель, представленная в общем виде, позволяет получить на ее основе более простые и эффективные модели для некоторых частных, но важных структур. Наконец, такие модели позволяют получить аналитические условия для расчета оптимальных параметров без использования вычислительно затратной оптимизации. Поэтому развитие аналитического моделирования не теряет актуальности.

Представим несколько примеров таких исследований. Например, хорошо известен подход к анализу перекрестных помех, задержки распространения сигнала и импульсных искажений в межсоединениях [20]. В [21] авторы приводят выражения в замкнутой форме для определения передаточных функций N связанных межсоединений с произвольными импедансами на концах. Заслуживают внимания также модели, основанные на численном обратном преобразовании Лапласа [22, 23], и аналитические модели для периодических многокаскадных структур из одиночных и симметричных связанных линий [24].

Между тем, насколько известно автору, к настоящему моменту нет окончательных аналитических моделей в замкнутой форме для временного отклика на концах витка МЛ на импульсное воздействие. Этот факт ограничивает получение аналитических оценок для важных приложений МЛ, например, в качестве устройств защиты от кондуктивных СКИ. Однако можно предположить, что подходы и методы, представленные в [2024], могут быть эффективно использованы для заполнения этих пробелов. В частности, модели из [24] успешно модифицированы в работах [25, 26] для сложных структур, но также могут быть использованы для витка МЛ. Это позволит сделать вклад в решение проблемы эффективного автоматизированного синтеза оптимальных МЛ.

Целью данной работы является вывод простых аналитических моделей для быстрой и априорной оценки амплитуд импульсов разложения, а также условий их выравнивания в конце витка МЛ.

1. ИСХОДНЫЕ ДАННЫЕ

На рис. 1 представлена схема соединений исследуемого витка МЛ. Он состоит из одного опорного и двух параллельных ему и друг другу сигнальных проводников длиной l, соединенных на одном конце. Один из сигнальных проводников соединен с источником электродвижущей силы (ЭДС) E с внутренним адмиттансом Y0, а другой – с нагрузкой Y0. Поскольку сигнал в конце витка МЛ представлен последовательностью его составляющих, то сначала нужно получить их амплитуды в аналитическом виде. Для этого удобно использовать аналитические модели для отрезка симметричной связанной линии передачи, полученные на основе моделей для одиночной линии с адмиттансом Y1 и задержкой τ1 (рис. 2) [24].

Рис. 1.

Схема соединений витка МЛ; прямоугольником обозначен отрезок регулярной связанной линии.

Рис. 2.

Эквивалентная схема одиночного отрезка линии передачи с окончаниями [24].

Для ясности дальнейшего изложения, представим эти модели. Компоненты отклика для дальнего конца, учитывающие проходящую волну V0(t), и ближнего конца, учитывающие отражение от начала $V_{1}^{'}\left( t \right)$ и конца отрезка линии, имеют следующий вид:

(1)
где Vвх(t) = E(t)/2.

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

(2)
$\begin{gathered} {{V}_{T}}\left( t \right) = {{V}_{0}}\left( t \right) + \sum\limits_{k = 1}^{{{k}_{{отр}}}/2} {{{V}_{{выхK}}}\left( t \right)} , \\ {{V}_{R}}\left( t \right) = V_{1}^{'}\left( t \right) + V_{1}^{{''}}\left( t \right) + \sum\limits_{k = 1}^{\left( {{{k}_{{отр}}} - 1} \right)/2} {{{V}_{{отрK}}}\left( t \right)} , \\ \end{gathered} $
где

$\begin{gathered} {{V}_{{выхK}}}\left( t \right) = \frac{{2{{Y}_{0}}}}{{{{Y}_{0}} + {{Y}_{1}}}}\frac{{2{{Y}_{1}}}}{{{{Y}_{1}} + {{Y}_{{\text{2}}}}}} \times \\ \times \,\,{{V}_{{вх}}}\left( {t - (2k + 1){{{{\tau }}}_{1}}} \right)\prod\limits_{i = 1}^k {\frac{{{{Y}_{1}} - {{Y}_{2}}}}{{{{Y}_{1}} + {{Y}_{2}}}}\frac{{{{Y}_{1}} - {{Y}_{0}}}}{{{{Y}_{0}} + {{Y}_{1}}}}} , \\ {{V}_{{отрK}}}\left( t \right) = \frac{{2{{Y}_{0}}}}{{{{Y}_{0}} + {{Y}_{1}}}}\frac{{2{{Y}_{1}}}}{{{{Y}_{0}} + {{Y}_{1}}}}\frac{{{{Y}_{1}} - {{Y}_{2}}}}{{{{Y}_{1}} + {{Y}_{2}}}} \times \\ \times \,\,{{V}_{{вх}}}\left( {t - 2(k + 1){{{{\tau }}}_{1}}} \right)\prod\limits_{i = 1}^k {\frac{{{{Y}_{1}} - {{Y}_{2}}}}{{{{Y}_{1}} + {{Y}_{2}}}}\frac{{{{Y}_{1}} - {{Y}_{0}}}}{{{{Y}_{0}} + {{Y}_{1}}}}} . \\ \end{gathered} $

Выражения (2) для одиночной линии применимы для симметричной (в поперечном сечении и по нагрузкам) связанной линии, если записать их отдельно, заменив индекс “1” в Y1 и τ1 на индексы “e” и “o” для четной и нечетной мод соответственно. Тогда компоненты отклика для четной и нечетной мод позволяют получить отклик в каждом узле связанной линии как

(3)
$\begin{gathered} {{V}_{1}}(t) = \frac{1}{2}\left[ {V_{R}^{e}(t) + V_{R}^{o}(t)} \right],\,\,\,\,{{V}_{2}}(t) = \frac{1}{2}\left[ {V_{R}^{e}(t) - V_{R}^{o}(t)} \right], \\ {{V}_{3}}(t) = \frac{1}{2}\left[ {V_{T}^{e}(t) + V_{T}^{o}(t)} \right],\,\,\,\,{{V}_{4}}(t) = \frac{1}{2}\left[ {V_{T}^{e}(t) - V_{T}^{o}(t)} \right], \\ \end{gathered} $
где V1(t) и V3(t) – отклики в начале и конце активного проводника, а V2(t) и V4(t) – пассивного.

2. МОДЕЛИ ДЛЯ ВЫЧИСЛЕНИЯ ВРЕМЕННОГО ОТКЛИКА

Структура на рис. 1 – это связанная линия с закороченными на дальнем конце проводниками без окончаний. Поэтому на дальнем конце (рис. 2, выражения (2)) Y2 = ∞ для нечетной моды и Y2 = = 0 – для четной. Тогда выражения (2) для каждой моды существенно упростятся. Используя (3), получим окончательно для отклика в узлах V1, V2, V3:

(4)
$\begin{gathered} {{V}_{1}}\left( t \right) = \frac{{{{V}_{{вх}}}(t)}}{2}\left( {\frac{{{{Y}_{0}} - {{Y}_{e}}}}{{{{Y}_{0}} + {{Y}_{e}}}} + \frac{{{{Y}_{0}} - {{Y}_{o}}}}{{{{Y}_{0}} + {{Y}_{o}}}}} \right) + \\ + \,\,2{{Y}_{0}}\left( {\frac{{{{Y}_{e}}{{V}_{{вх}}}(t - 2l{{\tau }_{e}})}}{{{{{({{Y}_{e}} + {{Y}_{0}})}}^{2}}}}} \right.\left. { - \frac{{{{Y}_{o}}{{V}_{{вх}}}(t - 2l{{\tau }_{o}})}}{{{{{({{Y}_{o}} + {{Y}_{0}})}}^{2}}}}} \right) + 2{{Y}_{e}}{{Y}_{0}} \times \\ \times \,\,\sum\limits_{i = 2}^{{{k}_{{отр}}}} {{{{( - 1)}}^{{i + 1}}}{{V}_{{вх}}}(t - 2l{{\tau }_{e}}i)} {{({{Y}_{e}} + {{Y}_{0}})}^{{ - (1 + i)}}}{{({{Y}_{0}} - {{Y}_{e}})}^{{i - 1}}} - \\ - \,\,2{{Y}_{o}}{{Y}_{0}}\sum\limits_{i = 2}^{{{k}_{{отр}}}} {{{V}_{{вх}}}(t - 2l{{\tau }_{o}}i)} {{({{Y}_{o}} + {{Y}_{0}})}^{{ - (1 + i)}}}{{({{Y}_{0}} - {{Y}_{o}})}^{{i - 1}}}, \\ \end{gathered} $
(5)
$\begin{gathered} {{V}_{2}}\left( t \right) = \frac{{{{V}_{{вх}}}(t)}}{2}\left( {\frac{{{{Y}_{0}} - {{Y}_{e}}}}{{{{Y}_{0}} + {{Y}_{e}}}} - \frac{{{{Y}_{0}} - {{Y}_{o}}}}{{{{Y}_{0}} + {{Y}_{o}}}}} \right) + \\ + \,\,2{{Y}_{0}}\left( {\frac{{{{Y}_{e}}{{V}_{{вх}}}(t - 2l{{\tau }_{e}})}}{{{{{({{Y}_{e}} + {{Y}_{0}})}}^{2}}}}} \right. + \left. {\frac{{{{Y}_{o}}{{V}_{{вх}}}(t - 2l{{\tau }_{o}})}}{{{{{({{Y}_{o}} + {{Y}_{0}})}}^{2}}}}} \right) + 2{{Y}_{e}}{{Y}_{0}} \times \\ \times \,\,\sum\limits_{i = 2}^{{{k}_{{отр}}}} {{{{( - 1)}}^{{i + 1}}}{{V}_{{вх}}}(t - 2l{{\tau }_{e}}i)} {{({{Y}_{e}} + {{Y}_{0}})}^{{ - (1 + i)}}}{{({{Y}_{0}} - {{Y}_{e}})}^{{i - 1}}} + \\ + \,\,2{{Y}_{o}}{{Y}_{0}}\sum\limits_{i = 2}^{{{k}_{{отр}}}} {{{V}_{{вх}}}(t - 2l{{\tau }_{o}}i){{{({{Y}_{o}} + {{Y}_{0}})}}^{{ - (1 + i)}}}{{{({{Y}_{0}} - {{Y}_{o}})}}^{{i - 1}}}} , \\ \end{gathered} $
(6)
$\begin{gathered} {{V}_{3}}\left( t \right) = V4\left( t \right) = \frac{{2{{Y}_{0}}{{V}_{{вх}}}(t - l{{\tau }_{e}})}}{{({{Y}_{e}} + {{Y}_{0}})}} + \\ + \,\,2{{Y}_{0}}\sum\limits_{i = 2}^{{{k}_{{отр}}}} {{{{( - 1)}}^{{i + 1}}}{{V}_{{вх}}}\left( {t - l{{\tau }_{e}}(2i - 1)} \right)} \times \\ \times \,\,{{({{Y}_{e}} + {{Y}_{0}})}^{{ - i}}}{{({{Y}_{0}} - {{Y}_{e}})}^{{i - 1}}}. \\ \end{gathered} $

Для проверки (4)–(6) вычислены временные отклики в узлах V1, V2, V3 (рис. 3) на импульсное воздействие в системе TALGAT численным методом в частотной области [27]. Параметры поперечного сечения витка, источника и нагрузки взяты из [7]. Полученные формы сигналов полностью совпали.

Рис. 3.

Формы напряжения в узлах V1 (а), V2 (б), V3 (в) схемы рис. 1, вычисленные по выражениям (4)–(6) (штриховые кривые) и в TALGAT (сплошные), полностью совпадают.

3. УСЛОВИЯ ДЛЯ ВЫРАВНИВАНИЯ АМПЛИТУД СОСТАВЛЯЮЩИХ ВРЕМЕННОГО ОТКЛИКА

Модель (5) пригодна для воздействия произвольной формы, но для определенности рассмотрим воздействие СКИ. Для выравнивания амплитуд импульсов разложения (исходного СКИ) в конце витка МЛ сначала необходимо определить нормированные амплитуды каждой составляющей временного отклика в конце МЛ. Из (5) легко получить аналитические выражения, определяющие амплитуды первого (перекрестной наводки, Vc), второго (нечетной моды, Vo) и третьего (четной моды, Ve) импульсов:

(7)
$\begin{gathered} {{V}_{c}} = \frac{{{{Y}_{0}}({{Y}_{o}} - {{Y}_{e}})}}{{({{Y}_{o}} + {{Y}_{0}})({{Y}_{e}} + {{Y}_{0}})}},\,\,\,\,{{V}_{o}} = \frac{{2{{Y}_{0}}{{Y}_{o}}}}{{{{{({{Y}_{o}} + {{Y}_{0}})}}^{2}}}}, \\ {{V}_{e}} = \frac{{2{{Y}_{0}}{{Y}_{e}}}}{{{{{({{Y}_{e}} + {{Y}_{0}})}}^{2}}}}. \\ \end{gathered} $

Далее можно сформулировать условия для выравнивания амплитуд импульсов и вывести окончательные выражения для нескольких частных случаев.

3.1. Выравнивание амплитуд двух импульсов

Если погонные задержки четной и нечетной мод линии одинаковы (например, при однородном диэлектрическом заполнении), то в конце линии наблюдаются только два импульса: наводки и основного сигнала. Амплитуда второго импульса определяется как сумма амплитуд четной и нечетной мод линии. Тогда из условия равенства амплитуд первого и второго импульсов Vc = = Vo + Ve с учетом (7) получим

(8)
$\frac{{{{Y}_{o}} + 3{{Y}_{e}}}}{{{{Y}_{0}}}} + \frac{{3{{Y}_{0}} + {{Y}_{e}}}}{{{{Y}_{o}}}} + \frac{{{{Y}_{0}} - {{Y}_{o}}}}{{{{Y}_{e}}}} = - 8.$

Рассмотрим частный случай минимизации отражений сигнала от концов проводников, когда ${{Y}_{0}} = \sqrt {{{Y}_{e}}{{Y}_{o}}} .$ Тогда после подстановки $k = \sqrt {{{{{Y}_{o}}} \mathord{\left/ {\vphantom {{{{Y}_{o}}} {{{Y}_{e}}}}} \right. \kern-0em} {{{Y}_{e}}}}} $ с физическим смыслом коэффициента связи, получим уравнение 4-й степени

${{k}^{4}} - 2{{k}^{3}} - 8{{k}^{2}} - 6k - 1 = 0,$

имеющее один физический корень k = $\sqrt 5 $ + 2 ≈ ≈ 4.236, который определяет нормированную к амплитуде Vвх(t) амплитуду напряжения на выходе витка МЛ как (k – 1)/(k + 1) = ($\sqrt 5 $ – 1)/2 ≈ 0.618.

Для проверки полученных аналитических моделей и оценок вычислены матрицы и временные отклики в узлах V1 и V2 на импульсное воздействие в системе TALGAT. Параметры воздействия, поперечного сечения МЛ в воздухе, генератора и нагрузки взяты из [6]. Вычисленные методом моментов [27] матрицы C и L структуры и вычисленная из них матрица Z имеют вид:

$\begin{gathered} {\mathbf{С}} = \left[ {\begin{array}{*{20}{c}} {153.16}&{ - 136.99} \\ { - 136.99}&{153.16} \end{array}} \right]\,\,{{{\text{пФ}}} \mathord{\left/ {\vphantom {{{\text{пФ}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{L}} = \left[ {\begin{array}{*{20}{c}} {363.35}&{325.00} \\ {325.00}&{363.35} \end{array}} \right]\,\,{{{\text{нГн}}} \mathord{\left/ {\vphantom {{{\text{нГн}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{Z}} = \left[ {\begin{array}{*{20}{c}} {108.93}&{97.43} \\ {97.43}&{108.93} \end{array}} \right]\,\,{\text{Ом}}. \\ \end{gathered} $

Как разность и сумма по строке матрицы Z, вычислены волновые сопротивления нечетной и четной мод линии (Zo = 11.50 Ом, Ze = 206.36 Ом), обеспечивающие условие $k = \sqrt {{{{{Z}_{e}}} \mathord{\left/ {\vphantom {{{{Z}_{e}}} {{{Z}_{o}}}}} \right. \kern-0em} {{{Z}_{o}}}}} \approx 4.236$. На рис. 4 представлены формы напряжений в узлах V1 и V2. Первые два импульса в конце витка имеют одинаковые амплитуды, около 61.8% от амплитуды Vвх(t).

Рис. 4.

Формы напряжения в узлах V1 (штриховая) и V2 (сплошная) витка МЛ в воздухе при k ≈ 4.236.

3.2. Выравнивание амплитуд трех импульсов

Если погонные задержки четной и нечетной мод линии различны (в неоднородном диэлектрическом заполнении), то их импульсы можно разнести по времени. Тогда, приравняв их амплитуды (Vo = Ve из (7)), получим

(9)
${{Y}_{0}} = \sqrt {{{Y}_{e}}{{Y}_{o}}} .$

Из условия равенства амплитуд первого и второго импульсов (Vc = Vo) получим

(10)
$\frac{{{{Y}_{0}} - 3{{Y}_{e}}}}{{{{Y}_{o}}}} + \frac{{{{Y}_{e}}}}{{{{Y}_{o}}}} = 1.$

Учитывая (9), при подстановке $k = \sqrt {{{{{Y}_{o}}} \mathord{\left/ {\vphantom {{{{Y}_{o}}} {{{Y}_{e}}}}} \right. \kern-0em} {{{Y}_{e}}}}} $ получим кубическое уравнение

${{k}^{3}} - {{k}^{2}} - 3k - 1 = 0,$

имеющее один физический корень k = $\sqrt 2 $ + 1 ≈ ≈ 2.414, который определяет нормированную амплитуду напряжения на выходе витка как (k – ‒ 1)/(k + 1) = $\sqrt 2 $ – 1 ≈ 0.414.

Для проверки моделей и оценок, параметры воздействия, поперечного сечения микрополосковой МЛ, источника и нагрузки взяты из [7]. Вычисленные методом моментов матрицы C и L структуры и вычисленная из них матрица Z имеют вид:

$\begin{gathered} {\mathbf{C}} = \left[ {\begin{array}{*{20}{c}} {232.06}&{ - 138.12} \\ { - 138.12}&{232.06} \end{array}} \right]\,\,{{{\text{пФ}}} \mathord{\left/ {\vphantom {{{\text{пФ}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{L}} = \left[ {\begin{array}{*{20}{c}} {390.34}&{309.03} \\ {309.03}&{390.34} \end{array}} \right]\,\,{{{\text{нГн}}} \mathord{\left/ {\vphantom {{{\text{нГн}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{Z}} = \left[ {\begin{array}{*{20}{c}} {50.55}&{35.73} \\ {35.73}&{50.55} \end{array}} \right]\,\,{\text{Ом}}{\text{.}} \\ \end{gathered} $

Из матрицы Z вычислены волновые сопротивления нечетной и четной мод линии (Zo = 14.82 Ом, Ze = 86.28 Ом), которые дают $k = \sqrt {{{{{Z}_{e}}} \mathord{\left/ {\vphantom {{{{Z}_{e}}} {{{Z}_{o}}}}} \right. \kern-0em} {{{Z}_{o}}}}} \approx 2.413$, что очень близко к $\sqrt 2 $ + 1≈ 2.414. На рис. 5 представлены формы напряжения в узлах V1 и V2. Первые три импульса в конце витка имеют одинаковую амплитуду, около 41.4% от амплитуды Vвх(t).

Рис. 5.

Формы напряжения в узлах V1 (штриховая) и V2 (сплошная) витка микрополосковой МЛ при k ≈ ≈ 2.413.

3.3. Дополнительное уменьшение амплитуд трех импульсов

Для дополнительного уменьшения амплитуды СКИ на выходе витка МЛ сначала надо обеспечить равенство моментов прихода импульса поздней моды (третий импульс на рис. 3б) и отраженного от начала линии импульса ранней моды отрицательной полярности (четвертый импульс на рис. 3б). Для этого нужно обеспечить условие [28]

${{\tau }_{{{\text{max}}}}} = 2{{\tau }_{{{\text{min}}}}},$
где τmax и τmin – погонные задержки поздней и ранней мод. Пусть далее для определенности они будут погонными задержками четной (τe) и нечетной (τo) мод, хотя может быть и наоборот.

Теперь дополнительно надо определить амплитуду четвертого импульса разложения (отраженного импульса нечетной моды, VRo) на выходе МЛ (см. рис. 3б), также нормированную относительно амплитуды Vвх(t):

(12)
${{V}_{{Ro}}} = - \frac{{2{{Y}_{o}}{{Y}_{0}}({{Y}_{o}} - {{Y}_{0}})}}{{{{{({{Y}_{o}} + {{Y}_{0}})}}^{3}}}}.$

Приравняв Vc и Vo из (7), получим

(13)
${{Y}_{0}} = \frac{{{{Y}_{o}}({{Y}_{o}} - 3{{Y}_{e}})}}{{{{Y}_{o}} + {{Y}_{e}}}}.$

Приравняв Vo и Ve + V из (7) и (12), получим

(14)
$\frac{{{{Y}_{e}}}}{{{{{({{Y}_{e}} + {{Y}_{0}})}}^{2}}}} = \frac{{2{{Y}_{o}}}}{{{{{({{Y}_{o}} + {{Y}_{0}})}}^{3}}}}.$

Подставив (13) в (14), получим уравнение

(15)
${{Y}_{e}}^{2} - {{Y}_{o}}^{2} + 4{{Y}_{o}}{{Y}_{e}} = 0,$

имеющее два физических корня:

(16)
${{{{Y}_{o}}} \mathord{\left/ {\vphantom {{{{Y}_{o}}} {{{Y}_{e}}}}} \right. \kern-0em} {{{Y}_{e}}}} = \sqrt 5 + 2 \approx 4.236,$
(17)
${{{{Y}_{o}}} \mathord{\left/ {\vphantom {{{{Y}_{o}}} {{{Y}_{e}}}}} \right. \kern-0em} {{{Y}_{e}}}} = \sqrt 5 --2 \approx 0.236.$

Отметим, что второй корень обратен первому, что означает взаимозамену значений Yo и Ye в (16) или (17), поэтому далее для определенности используем первый корень. Выразив Yo из (16) и подставив в (13), получим

(18)
${{Y}_{0}} = {{Y}_{e}}.$

Таким образом, при выполнении (16) достаточно выполнить (18) вместо (13).

Рассмотрим случай, когда (16) и (18) выполняются, так что Vc = Vo = Ve + VRo. Тогда вычислить нормированную амплитуду первых трех импульсов можно, например, по второму выражению из (7) после подстановки в него (13) с учетом (16):

(19)
${{V}_{o}} = {{(\sqrt 5 --1)} \mathord{\left/ {\vphantom {{(\sqrt 5 --1)} 4}} \right. \kern-0em} 4} \approx 0.309.$

Таким образом, при выполнении (11), (16) и (18) амплитуда первых трех импульсов на выходе линии должна составлять 30.9% от амплитуды Vвх(t). Для проверки этого вычислен временной отклик в узле V2 на воздействие СКИ численным методом в системе TALGAT. Параметры поперечного сечения выбраны такими же, как в [28], поскольку они обеспечивают приблизительное выполнение (11) и равенства амплитуд первых трех импульсов. Вычисленные методом моментов матрицы C и L структуры и вычисленная из них матрица Z имеют вид:

$\begin{gathered} {\mathbf{C}} = \left[ {\begin{array}{*{20}{c}} {1100.42}&{ - 396.78} \\ { - 396.78}&{1100.42} \end{array}} \right]\,\,{{{\text{пФ}}} \mathord{\left/ {\vphantom {{{\text{пФ}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{L}} = \left[ {\begin{array}{*{20}{c}} {219.04}&{172.95} \\ {172.95}&{219.04} \end{array}} \right]\,\,{{{\text{нГн}}} \mathord{\left/ {\vphantom {{{\text{нГн}}} {\text{м}}}} \right. \kern-0em} {\text{м}}}, \\ {\mathbf{Z}} = \left[ {\begin{array}{*{20}{c}} {14.575}&{9.027} \\ {9.027}&{14.575} \end{array}} \right]\,\,{\text{Ом}}{\text{.}} \\ \end{gathered} $

Квадратный корень собственных значений произведения матриц C и L дает погонные задержки четной (τe = 16.608 нс/м) и нечетной мод (τo = = 8.307 нс/м). Тогда 2τo = 16.614 нс/м, что отличается от τe = 16.608 нс/м на 0.04%. Из матрицы Z вычислены волновые сопротивления нечетной и четной мод линии: Zo=5.548 Ом, Ze=23.603 Ом. Тогда Ze/Zo = Yo/Ye = 4.254, что отличается от значения в (16) на 0.423%. Полученное выше значение Ze=23.603 Ом отличается от сопротивления 23 Ом, использованного в [28], на 2.555%. Таким образом, условия (11), (16) и (18) выполняются с хорошей точностью.

Далее был вычислен отклик численным методом в системе TALGAT на воздействие СКИ с амплитудой ЭДС 1 В, длительностью плоской вершины 400 пс, а фронта и спада – по 50 пс. Длина линии l = 45 мм. На рис. 6 показана форма напряжения в узле V2 микрополосковой МЛ, из которой видно, что сигнал в конце линии представлен последовательностью из трех импульсов с близкими амплитудами: перекрестной наводки (0.156 В), нечетной моды (0.157 В), а также суммы четной моды и отраженной нечетной моды (0.154 В), составляющими 31.2, 31.4 и 30.8% от амплитуды Vвх(t) соответственно, что отличается от 30.9% из (19) менее чем на 0.5%.

Рис. 6.

Формы напряжения в узлах V1 (штриховая) и V2 (сплошная) витка микрополосковой МЛ при выполнении (11), (16) и (18).

Таким образом, полученное аналитически значение амплитуды в конце линии подтверждается моделированием численным методом, а выполнение (11), (16) и (18) позволяет выровнять и дополнительно уменьшить амплитуды импульсов в конце витка МЛ.

ЗАКЛЮЧЕНИЕ

В работе получены и проверены аналитические модели для расчета отклика на произвольное кондуктивное воздействие во временной области в витке МЛ. Разработанные модели позволили получить условия равенства составляющих отклика в конце витка МЛ на примере воздействия СКИ для трех случаев его разложения: на два импульса (при равных погонных задержках); на три импульса (при разных погонных задержках); на три импульса с дополнительным уменьшением амплитуды (за счет наложения на импульс поздней моды отраженного импульса ранней моды отрицательной полярности). Аналитические модели и условия подтверждены численным моделированием на основе метода моментов и модифицированного узлового метода.

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

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

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

Список литературы

  1. Misman D., Aziz M.Z.A. Abd., Husain M.N., Soh P.J. // Proc. 2009 Europ. Conf. Antennas and Propagation. Berlin. 27–23 Mar. N.Y.: IEEE, 2009. P. 2420.

  2. Acuña J.E., Rodriguez J.L., Obelleiro F. // Int. J. RF and Microwave Computer-Aided Engineering. 2001. V. 11. № 4. P. 219.

  3. Wang X.C., Sun Y.Y., Zhu J.H. // IEEE Trans. 2017 V. EMC-59. № 3. P. 996.

  4. Shlepnev Y., Neves A., Dagostino T., Mcmorrow S. // Proc. 2009 Design Conf. Santa Clara. 2–5 Feb. S.L.C.: UBM Electronics, 2009. P. 2581.

  5. Zhang Q., Bandler J.W., Caloz C. // IEEE Trans. 2013 V. MTT-18. № 4. P. 4040.

  6. Surovtsev R.S., Gazizov T.R., Zabolotsky A.M. // Proc. 2015 Int. Siberian Conf. Control and Communications (SIBCON). Omsk. 21–23 May. N.Y.: IEEE, 2015. P. 7147227.

  7. Surovtsev R.S., Nosov A.V., Zabolotsky A.M. // Proc. 2015 16th Int. Conf. Young Specialists on Micro/Nanotechnologies and Electronic Devices. Erlagol, 29 Jun.–3 Jul. N.Y.: IEEE, 2015. P. 175.

  8. Nosov A.V., Surovtsev R.S. // Proc. 2020 21st Int. Conf. Young Specialists on Micro/Nanotechnologies and Electronic Devices. Chemal. 29 Jun.–3 Jul. N.Y.: IEEE, 2020. P. 149.

  9. Петкау О., Тарабцев А., Дерябин А. и др. // Безопасность объектов топливно-энергетического комплекса. 2014. № 2. С. 74.

  10. Газизов Т.Р. // Успехи совр. радиоэлектрон. № 2. 2004. С. 37.

  11. Макаренко С.И. // Системы управления, связи и безопасности. 2020. № 2. С. 101.

  12. Здухов Л.Н., Парфёнов Л.Н., Тарасов О.А., Чепелев В.М. // Технологии ЭМС. 2018. № 2. С. 22.

  13. Гизатуллин З.М., Гизатуллин Р.М. // РЭ. Т. 61. № 5. 2016. С. 500.

  14. Газизов А.Т., Заболоцкий А.М., Газизов Т.Р. // РЭ. 2018. Т. 63. № 3. С. 292.

  15. Gazizov A.T., Zabolotsky A.M., Gazizov T.R. // IEEE Trans. 2016. V. EMC-58. № 4. P. 1136.

  16. Surovtsev R.S., Nosov A.V., Zabolotsky A.M., Gazizov T.R. // IEEE Trans. 2017. V. EMC-59. № 6. P. 1864.

  17. Wu R.B., Chao F.L. // IEEE Trans. 1995. V. CPMT-18. № 4. P. 644.

  18. Bhobe A.U., Holloway C.L., Piket-May M. // Proc. 2001 Int. Symp. EMC. Montreal. 13–17 Aug. N.Y.: IEEE, 2001. P. 805.

  19. Nosov A.V., Surovtsev R.S., Gazizov T.T., Gazizov T.R. // Proc. 2017 Int. Multi-Conf. Engineering, Computers and Information Sciences (SIBIRCON). Novosibirsk. 18–22 Sept. N.Y.: IEEE, 2017. P. 453.

  20. You H., Soma M. // IEEE Trans. 1990. V. CS-37. № 8. P. 1019.

  21. You H., Soma M. // Proc. IEEE Custom Integrated Circuits Conf. Boston 13–16 May 1990. N.Y.: IEEE, 1990. P. 11.2.

  22. Xiao F., Murano K., Kami Y. // IEICE Trans. Commun. 2009 V. 92B. № 6. P. 1953.

  23. Griffith J.R., Nakhla M.S. // IEEE Trans. 1990 V. MTT-38. № 10. P.1480.

  24. Gu Q., Kong J.A. // IEEE Trans. 1986. V. MTT-34. № 9. P. 952.

  25. Gazizov T.R., Leontiev N.A. // Proc. 4th Int. Symp. Antennas and EM Theory (ISAE 97). Xi’an, 19–22 Aug. Beijing: Int. Academic Publ., 1997. P. 444.

  26. Gazizov T.R., Leontiev N.A. // Proc. 1997 Sino-Japanese Joint Meet. Opt. Fiber Sci. and EM Theory. Wuhan, 14–16 Oct. P. 322.

  27. Kuksenko S.P. // IOP Conf. Ser.: Mater Sci. and Engineering. 2019. V. 560. P. 012110.

  28. Nosov A.V., Surovtsev R.S., Gazizov T.T., Gazizov T.R. // Proc. 2017 18th Int. Conf. Young Specialists Micro/Nanotech. Elec. Devices. Erlagol. 29 Jun.–3 Jul. N.Y.: IEEE, 2017. P. 119.

Дополнительные материалы отсутствуют.