CC BY
69
Труды БГТУ, 2016, № 6, с. 127-131
127
УДК 621.373.52:539.1
М. А. Анкуда, И. О. Оробей
Белорусский государственный технологический университет
МОДЕЛЬ АВТОДИНА С ПАРАМЕТРИЧЕСКИМ ВОЗБУЖДЕНИЕМ
В данной статье рассмотрен анализатор для оценки качественных параметров технологических процессов, который основан на эффекте ядерного магнитного резонанса. Приведено математическое описание протекающих процессов в электронных компонентах первичного преобразователя, которое позволяет анализировать динамику процессов и оценивать влияние внешних воздействий. Описана работа параметрического усилителя мощности, в котором энергия от изменения реактивного параметра преобразуется в энергию сигнала. Усиление в предложенной параметрической схеме происходит за счет внесения генератором накачки в сигнальный контур условного отрицательного сопротивления, которое зависит от амплитуды накачки. В данном генераторе в качестве переменного реактивного элемента колебательного контура используется емкость p-n перехода полупроводникового варикапа. В качестве генератора накачки применяется схема возведения в квадрат, а петля обратной связи всего параметрического генератора реализуется через подачу сигнала на управляющий электрод варикапа после определенного усиления. Разработанный автодинный генератор обеспечивает повышенную точность измерений за счет увеличения помехоустойчивости, которая достигается благодаря параметрической обратной связи с частотой накачки, равной удвоенной частоте сигнала. Решение предложенной модели генератора реализовано с помощью метода медленно меняющихся амплитуд.
Ключевые слова: автодин, колебательный контур, параметрический резонанс, ядерный магнитный резонанс, генератор накачки, шумовой сигнал.
M. A. Ankuda, I. O. Orobei
Belarusian State Technological University
THE MODEL OF AUTODYNE WITH PARAMETRIC EXCITATION
In this article we have studied the analyzer to evaluate the quality parameters of the technological processes, which are based on the effect of nuclear magnetic resonance. The mathematical description of the processes in the electronic components of the primary device, which allows you to analyze the dynamics of the process and assess the impact of external influences is presented. The functioning of the parametric amplifier, in which the energy from the changes reactive parameter is converted into the energy of the signal is described. Parametric amplification in the proposed scheme is due to the introduction of the generator pumping conditional negative resistance in signal circuit that is independent of the pump amplitude. In this generator as a variable reactance element is used oscillation circuit with capacitance p-n transition semiconductor varicap. The pump generator is used by the squaring circuit. The feedback loop of all parametric oscillator is realized through a signal to the control electrode varicap after a certain gain. The designed autodyne provides improved measurement accuracy by increasing the noise immunity, which is achieved due to the parametric feedback from the pump frequency equal to twice the signal frequency. The solution of proposed model is implemented by the method of slowly varying amplitudes.
Key words: autodyne, oscillation circuit, parametric resonance, nuclear magnetic resonance, generator of pumping, noise signal.
Введение. Для оценки качественных параметров технологических процессов широкое распространение получили анализаторы, которые основаны на изучении интенсивности сигналов ядерного магнитного резонанса (ЯМР). Математическое описание процессов в электронных компонентах первичных преобразователей ЯМР-измерителей позволяет более точно изучить динамику происходящих процессов, провести анализ сигналов и оценить реакцию при различном внешнем воздействии.
Основная часть. Метод автодинного детектирования спиновых резонансов достаточ-
но актуален в технике стационарной спектроскопии, в частности ЯМР [1]. Современные методы обработки сигналов позволяют существенно расширить возможности автодинных радиоспектрометров. Автодин, по сути, представляет собой генератор слабых высокочастотных колебаний, амплитуда которых зависит от внесенных в колебательный контур потерь за счет ядерной динамической восприимчивости образца [2]. Поглощение энергии радиочастотного поля на резонансной частоте ядерной спиновой системы исследуемого образца приводит к уменьшению добротности контура
и, как следствие, вызывает изменение амплитуды высокочастотных (ВЧ) колебаний. Для эффективного поглощения спиновой системой энергии радиочастотного поля необходимо, чтобы добротность катушки с исследуемым веществом была максимальной.
Обычно параллельный колебательный контур автодинного детектора характеризуется индуктивностью Ь, емкостью С и сопротивлением Я, которые при своей независимости от сигналов в контуре образуют линейную систему, описываемую следующим дифференциальным уравнением:
1 г 7 и ^ёи , — I иЖ + — + С— = ?вых, (1)
Ь* Я ж
где и - напряжение колебательного контура; ¿вых - ток с выхода положительной обратной связи усилителя.
Если какой-либо из параметров контура зависит от величины сигнала в нем, то уравнение (1) становится нелинейным, а для нелинейной системы нет принципиального различия в получении резонанса либо воздействием на систему внешним ЭДС, либо параметрическим возбуждением. Периодическое изменение реактивного параметра при достаточной интенсивности и соответствующей частоте вызывает параметрическое возбуждение колебаний в системе [3]. Подобная система представляет собой усилитель мощности, в котором энергия от изменения реактивного параметра преобразуется в энергию сигнала. Усиление в параметри-
ческом усилителе происходит за счет внесения генератором накачки в сигнальный контур условного отрицательного сопротивления, которое зависит от амплитуды накачки.
Однако, помимо создания соответствующей амплитуды колебаний генератора накачки, необходимых для изменения емкости контура, требуется также обеспечить совпадение фаз сигналов накачки и контура. Причем частоту накачки следует выбирать в 2 раза больше частоты сигнала, чтобы обеспечить подвод мощности к контуру [3].
Таким образом, схема автодина с параметрическим возбуждением может быть условно разделена на несколько функциональных узлов: колебательный контур, генератор накачки и цепь обратной связи.
Схема автодина приведена на рисунке.
В рассматриваемом генераторе в качестве переменного реактивного элемента колебательного контура используется емкость р-п перехода полупроводникового варикапа УБ1; в качестве генератора накачки применяется схема возведения в квадрат, которая построена на БЛ1. Поскольку в качестве опорного сигнала для параметрического изменения емкости используется непосредственно сигнал колебательного контура на удвоенной частоте, которая получается с выхода схемы умножения, то тем самым обеспечивается совпадение фаз. Обратная связь реализуется через подачу сигнала после усиления на схеме БЛ2 на управляющий электрод варикапа УБ1.
С2
I
С1
Х1 Х2 XX • п и У1 У2
УР УЫ
W г
КЗ
БЛ1 ЛБ835
Ь1
Е 1
УБ1 ~|
"X
+5
УБ1 ББ212
Я8 С8
Ъ.
С3
С4
"X
Я4 Я.5
К8
УР
г оит
-м +Ш С УЫ
+5
X
БЛ2 ЛБ8099ЛСР
£
С6
С7
С5
Электрическая принципиальная схема автодина
Для составления математической модели рассмотрим колебательный контур, емкость которого представляет собой некоторую изменяемую величину. При этом зависимость емкости варикапа УБ1 от приложенного обратного напряжения является обратно пропорциональной величиной С = 1 / (1 + иобр). Поскольку обратное напряжение представлено гармоникой на частоте накачки и изменение емкости будет однозначно пропорционально сопоставляться изменению напряжения, то значение емкости варикапа можно записать в виде
C = -
Cn
m =
1 + m cos2®t
AC C - C
^^ = VD1max ^VDlmin
2C " C + C
^ VDlmax ^VDlmm
(2)
где С0 - значение емкости варикапа при фиксированном обратном напряжении; т - коэффициент модуляции емкости; СУ01тах - значение емкости варикапа при максимальном отклонении от фиксированного обратного напряжения; Суо1тт - значение емкости варикапа при минимальном отклонении от фиксированного обратного напряжения.
При определенных условиях подобный колебательный контур может стать параметрическим генератором.
Для получения математической модели в динамике продифференцируем уравнение (1):
1 1 du d 2u
—u +--+ C—— = 0.
L1 R1 dt dt2
(3)
Тогда подставляя уравнение (2) в (3) для колебательного контура можно будет получить:
1
—г L1
R1
CoC1
d 2u
C1(1 + m cos 2 ®t) + C0
= 0.
(4)
Уравнение (4) будет иметь место при условии, что т < 1, что является справедливым, так как АС - небольшое отклонение от С0. Если пренебречь потерями на Ш и принять значение С1 равным значению С0, то уравнение (4) можно записать в следующем виде:
ю2 (2 + mcos2rat)u -
d u
~dF
= 0,
(5)
где
ю0 =
л/LíCO
Решим (5) методом медленно меняющихся амплитуд. При малой нелинейности емкости и малой глубине модуляции параметра можно ограничиться рассмотрением гармонического приближения и искать решение п(?) как:
u = a(t) sin ®t + b (t) cos ®t,
(6)
где а(0 и Ъ(?) - медленно меняющиеся во времени амплитуды искомого колебания.
Медленность изменения а и Ъ во времени определяется тем, что их производные по времени являются величинами первого порядка малости по сравнению с произведениями юа и юЪ:
—«юа, — «юЬ. da da
С учетом (7) первая производная:
du ,
— « аю cos rat - bra sin rat + dt
da db +sin rat--+ cos rat—
в первом приближении принимает вид
du
(7)
dt
- ara cos rat - bra sin ra t.
(8)
Для второй производной:
í -ra2 a sin rat -ra2b cos rat +
d 2u
dt2
db
+ racos rat--ra sin rat--+
db dt
2
.."a . da
+--—sin rat +--—cos rat +
dt dt
da db + racos rat--ra sin rat—
пренебрежем слагаемыми второго порядка малости и оставим слагаемые первого порядка малости. Получим:
d2u ( 2 ~ db , . —— ra a + 2га— |sinrat-dt2
i 2, - da , + | -ra b + 2ra — I cos ra t. dt
(9)
Слагаемые первого порядка малости оставлены в (9) и ими пренебрегли в (8). Это объясняется тем, что схема обладает малыми потерями, поэтому амплитуда второго слагаемого левой части исходного уравнения относительно мала по сравнению с амплитудами первого и третьего слагаемых.
1
Подставив (6), (9) в уравнение (5), получим: (2 + m cos2at )(a sin at + b cos at) -
i 2 ~ db | . -I a a + 2ю— I sin at +
da
+ I-ю2Ь + 2ю— I cos at = 0. (10)
V dt J
Поскольку на систему не действует внешняя периодическая сила и частота автоколебаний равна ю0, а не ю, то (10) примет вид
ю2 (2 + m cos2a0t )a sin a0t -da
-ю02а sin a0t + 2ю0—cos a0t = 0. (11)
dt
Введем безразмерный параметр т = a0t, а производную будем брать по времени. Тогда решение уравнения (11) можно будет записать в виде
K1 ,
a (t) = Be
(12)
где
K1 = (1 + mcos2x)sin т, K 2 = 2ю0 cos т.
В момент времени t = 0 амплитуда определяется шумом Найквиста:
a(t = 0) = yj 4kTRAf,
(13)
где a(t = 0) - амплитуда напряжения при t = 0; к = 1,38 • 10-23 - постоянная Больцмана; T -температура, К; R - общее сопротивление, Ом; А/ - ширина полосы.
В итоге общее решение примет следующий вид:
®0 (1+mcos2x) tg т
u (t) = y¡ 4kTRA/e 2 t sin a0t. (14)
Для получения временной зависимости сигналов ЯМР на выходе автодина нужно выявить влияние на его работу движения вектора намагниченности в катушке регистрации. При помещении образца в постоянное поле B0, направленное по 0z, и наложении на него переменного вращающегося с частотой ю в плоскости 0xy поля Br уравнениями движения вектора намагниченности являются уравнения Блоха [2]:
dMx y(Bz -MzBy) + Mjr = 0, (15)
dt
dM. dt
T2
-j(MzBx -MXBZ)+My = 0, (16)
T o
-у(МхВу -МуВх) = , (17)
где Мх, Му, М2 - мгновенные значения проекций векторов намагниченности результирующего поля в образце на оси 0х, 0у, 0г, А/м; Вх, Ву, В2 -мгновенные значения проекций векторов индукции результирующего поля в образце на оси 0х, 0у, 02, Вб/м2; М0 - равновесная намагниченность (М0 = Х0В0), А/м; Вг = В0; Ть Т2 - соответственно спин-решеточная и спин-спиновая релаксации вещества образца, с.
Общее решение (15)-(17) затруднено в связи с необходимостью прохождения резонансных условий, что требует изменения во времени В2 = В0 или частоты ю. Поэтому переходят в ВСК с частотой ю вокруг 0г, полагая
Mx = ucosa t - vsina t, My = + (cosa t + usina t),
(18)
где и - намагниченность, вращающаяся синфаз-но с Вг, А/м; V - намагниченность, опережающая Вг на п/2, А/м.
С учетом (18) решение в комплексной форме даст вид
V + ju = -|у| BXM0 X
t (t'-1
J exp —--jjAa(t")dt"
v T t
(19)
Спектр ЯМР снимают прохождением резонансных условий, изменением поля В0 при постоянной частоте ю или изменением ю при постоянном В0. Уравнение (19) анализируется при разных вариантах изменения В0.
Заключение. Таким образом, была рассмотрена возможность создания генератора автоколебаний с параметрическим возбуждением; приведена его математическая модель и способ определения ее решения.
Разработанный автодинный генератор обеспечивает повышенную точность измерений в сравнении с классическими схемами автоди-нов, которые работают на усилителях с положительной обратной связью, за счет увеличения помехоустойчивости. Предложенная модель генератора имеет параметрическую обратную связь с частотой накачки, равной удвоенной частоте сигнала. Поскольку шум является обратно пропорциональной величиной частоте сигнала, то подобное удвоение необратимо приводит к повышению соотношения полезного сигнала к шумовой составляющей. Последнее является неоспоримым преимуществом при создании схем регистрации сигналов ЯМР.
Литература
1. Яценко А. В., Евдокимов С. В. Датчик сигналов ЯМР для исследования веществ с большим значением // ПТЭ. 2003. № 1. С. 64-66.
2. Леше А. Ядерная индукция. М.: Изд-во иностранной лит-ры, 1963. 684 с.
3. Мандельштам Л. И., Папалекси Н. Д. О параметрическом возбуждении электрических колебаний // Журн. техн. физики. 1934. № 3. С. 5-29.
References
1. Yatsenko A. V., Yevdokimov S. V. The sensor signals of NMR for research of materials with high. Pribory i tekhnika eksperimenta [Instruments and Experimental Technique], 2003, no. 1, pp. 64-66 (In Russian).
2. Leshe A. Yadernaya induktsiya [Nuclear induction]. Moscow, Izdatel'stvo inostrannoy literatury Publ., 1963. 684 p.
3. Mandel'shtam L. I., Papaleksi N. D. Parametric excitation of electrical oscillations. Zhurnal tekhnicheskoy fiziki [Journal of Technical Physics], 1934, no. 3, pp. 5-29 (In Russian).
Информация об авторах
Анкуда Максим Анатольевич - ассистент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: M.Ankuda@belstu.by
Оробей Игорь Олегович - кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: I.Orobey@belstu.by
Information about the authors
Ankuda Maksim Anatol'yevich - assistant lecturer, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: M.Ankuda@belstu.by
Orobei Igor Olegovich - PhD (Engineering), Assistant Professor, Assistant Professor, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: I.Orobey@belstu.by
Поступила 01.03.2016
