CC BY
63
Вестник МГТУ, том 13, №4/2, 2010 г.
стр.931-936
УДК 681.586.621.37:543.275.1
Радиочастотное резонаторное устройство для измерения физических свойств жидкостей
А.А. Маслов , А.С. Совлуков ’ , В.Я. Фатеев , В.В. Яценко , А.Р. Власова
1 Политехнический факультет МГТУ, кафедра автоматики
и вычислительной техники
2
Институт проблем управления им. В.А. Трапезникова РАН, г. Москва
Аннотация. Рассматривается радиочастотное резонаторное измерительное устройство для определения физических свойств жидкостей. Устройство содержит отрезок коаксиальной линии с U-образным внутренним проводником. Электронный блок устройства может быть расположен на значительном расстоянии от контролируемого объекта. Предложена математическая модель чувствительного элемента устройства, позволяющая рассчитывать выходную характеристику устройства и оптимизировать конструктивные параметры чувствительного элемента с точки зрения повышения его чувствительности к измеряемым физическим свойствам жидкости и обеспечения однозначности измерений.
Abstract. Radiofrequency (RF) resonator measuring device for determination of physical properties of liquids has been proposed. The device contains RF sensor that is a section of modified coaxial line with U-shaped inner conductor. Distant location of electronic unit of the device from a monitored object can be provided. Mathematical model of the RF sensor has been proposed. It gives ability to calculate output characteristic of the device and to optimize parameters of the RF sensor in order to increase its sensitivity to measured physical properties of a liquid and to provide unambiguity of measurement results.
Ключевые слова: радиочастотное резонаторное устройство, коаксиальная линия, U-образный внутренний проводник, определение физических свойств жидкости
Key words: radiofrequency resonator, coaxial line, U-shaped inner conductor, determination of physical properties of liquids
1. Введение
Во многих отраслях промышленности необходимо с высокой точностью определять различные физические свойства (концентрацию смеси веществ, влагосодержание, плотность и др.) жидкостей, находящихся в емкостях (технологических резервуарах, измерительных ячейках и т.п.). В частности, в пищевой промышленности требуется измерять концентрацию коптильных препаратов, водно-спиртовых растворов и др.
Известны различные устройства для определения физических свойств жидкостей, основанные на измерении электрофизических параметров (диэлектрической проницаемости или/и тангенса угла диэлектрических потерь) жидкостей с применением радиоволновых ВЧ и СВЧ резонаторов, содержащих контролируемую жидкость (Брандт, 1963; Викторов и др., 1978; 1989). Такие измерительные устройства имеют, однако, ограниченную область применения, обусловленную возможностью контроля жидкостей, являющихся достаточно хорошими диэлектриками: для диэлектрических жидкостей с большими диэлектрическими потерями добротность резонаторов может быть малой величиной и, соответственно, амплитуда информативных сигналов мала для надежной регистрации.
В известном устройстве (Ершов и др., 2007) содержится располагаемый в контролируемой жидкости чувствительный элемент (ЧЭ), представляющий собой отрезок длинной линии в виде совокупности металлического цилиндра и расположенного внутри него и параллельно ему внутреннего проводника U-образной формы, подключенных к электронному блоку для возбуждения в отрезке длинной линии электромагнитных колебаний и измерения его резонансной частоты. Однако данное устройство имеет ограниченную область применения, что обусловлено необходимостью расположения электронного блока непосредственно вблизи ЧЭ и контролируемого объекта: в качестве электронного блока применяется автогенератор, в частотозадающую цепь которого подключен рассматриваемый отрезок длинной линии. Частота автогенератора, зависящая функционально от резонансной частоты отрезка длинной линии, служит информативным параметром. На практике же существует много задач, связанных с необходимостью расположения электронного блока на значительном расстоянии (несколько метров и более) от контролируемого объекта. Схема автогенератора с рассматриваемым ЧЭ не позволяет это сделать.
Рассматриваемое здесь устройство характеризуется расширенной областью применения и обеспечивает возможность проведения измерений при удалении электронного блока от контролируемого объекта и уменьшение влияния на погрешность измерения параметров кабелей, связывающих генератор и блок для измерения резонансной частоты с соответствующими индуктивностями связи.
931
Маслов А. А. и др. Радиочастотное резонаторное устройство...
Рис. 1. Функциональная схема измерительного устройства
1 - контролируемая жидкость,
2 - внутренний проводник отрезка длинной линии,
3 - наружный проводник отрезка длинной линии,
4 - нагрузочные индуктивности,
5 - индуктивности связи,
6 и 7 - высокочастотные трансформаторы,
8 - генератор частотно-модулированных колебаний,
9 - блок измерения резонансной частоты
2. Функциональная схема измерительного устройства
На рис. 1 изображена функциональная схема измерительного устройства. Здесь чувствительный элемент в виде отрезка длинной линии размещается в емкости с контролируемой жидкостью 1; он может быть также заполнен этой жидкостью для проведения измерений. Такой ЧЭ содержит внутренний проводник 2, имеющий U-образную форму, и наружный проводник (металлический экран) 3. Внутренний проводник 2 может быть покрыт диэлектрической оболочкой определенной толщины, изготовленной из фторопласта, полиэтилена и др., что важно при контроле жидкостей с большими диэлектрическими потерями.
Каждый из концов данного отрезка длинной линии подключен к соответствующей нагрузочной индуктивности 4. Эти нагрузочные индуктивности 4 вместе с соответствующими индуктивностями связи 5 образуют высокочастотные трансформаторы 6 и 7, которые работают на частотах ~ 1 МГц и более высоких частотах. К одной из индуктивностей связи 5 подключен генератор частотно-модулированных колебаний 8, с помощью которого в данном отрезке длинной линии - полуволновом резонаторе -возбуждают электромагнитные колебания (обычно в мегагерцовом диапазоне частот электромагнитных волн) на основной резонансной частоте. К другой индуктивности связи 5 подключен блок измерения резонансной частоты 9. Величина индуктивной связи между первичной (индуктивностью связи 5) и вторичной (нагрузочной индуктивностью 4) обмотками каждого из высокочастотных трансформаторов 6 и 7 выбирается из условий обеспечения большого значения нагрузочной (т.е. при подключении внешних элементов к отрезку длинной линии) добротности резонатора (отрезка длинной линии) и достаточно большого значения амплитуды сигнала на выходе данного резонатора, т. е. на входе блока измерения резонансной частоты 9.
В частности, при измерении физических свойств жидкости, являющейся диэлектриком с большими диэлектрическими потерями, собственная (т. е. в отсутствие внешних элементов) добротность резонатора (отрезка длинной линии) имеет малую величину (~ 10 ^ 20). В этом случае основной задачей является обеспечение достаточно большого значения амплитуды выходного сигнала, поступающего в блок измерения резонансной частоты 9. Поэтому в данном случае требуется обеспечить достаточно сильную индуктивную связь между обмотками каждого трансформатора. Это достигается расположением обмотки индуктивности связи 5 каждого из трансформаторов 6 и 7 в непосредственной близости от обмотки соответствующей нагрузочной индуктивности 4 (например, одна из этих обмоток располагается поверх другой обмотки или между витками другой обмотки).
3. Принцип действия измерительного устройства
При совпадении частоты частотно-модулированных колебаний и резонансной частоты полуволнового резонатора амплитуда этих электромагнитных колебаний достигает максимума. Этот максимум фиксируется в блоке измерения резонансной частоты 9, и в нем происходит преобразование измеряемой частоты в какой-либо из электрических сигналов (напряжение, ток) или цифровой код. При расположении полуволнового отрезка длинной линии в контролируемой жидкости 1 его резонансная частота изменяется. Поскольку эта резонансная частота функционально зависит от электрофизических параметров контролируемой жидкости 1, в частности, от диэлектрической проницаемости, которая, в свою очередь, зависит от физических свойств этой жидкости (концентрации смеси веществ, влагосодержания, плотности и др.), то по ее значению определяют соответствующее значение измеряемого физического свойства жидкости.
В общем случае ЧЭ в виде отрезка длинной линии можно представить в виде эквивалентного LC-контура с сосредоточенными параметрами, резонансная частота которого (Ор определяется по формуле
ю
р
1
4lc
(1)
932
Вестник МГТУ, том 13, №4/2, 2010 г. стр.931-936
Измеряемый физический параметр чаще всего воздействует на электрическую емкость С данного эквивалентного LC-контура, в результате чего она изменяется на величину АС, а резонансная частота - на величину Аа>р. Дифференцируя функцию (1) по С, можно получить приближенную формулу, связывающую величины Аа>р и АС
А(°р _ _4(С (2)
Юр 2С
Величина А&р/ор определяет частотное перекрытие датчика, т.е. относительное изменение информативного параметра (резонансной частоты) от относительного изменения измеряемого параметра.
Следует заметить, что емкость С состоит из рабочей емкости Ср и паразитной емкости Сп. Паразитную емкость составляют емкости различных элементов конструкции датчика, не участвующие в преобразовании "измеряемый параметр ^ С ^ op", например, емкости проходных изоляторов, фиксаторов и линий связи. Очевидно, что максимальную чувствительность и минимальную погрешность измерения можно получить только при максимально возможном значении АС/(2С). Поэтому необходимо стремиться к уменьшению влияния паразитной емкости Сп на величину частотного перекрытия (2). В частности, уменьшения влияния емкости линий связи можно добиться при подключении их к точкам резонатора, в которых величина электрического поля минимальная. В полуволновом резонаторе такими точками являются короткозамкнутые концы этого резонатора или концы резонатора, нагруженные на индуктивность LH.
В этом случае входное сопротивление Zex отрезка длинной линии, имеющего длину l и нагруженного на противоположном конце на индуктивность LH, можно определить из следующего выражения (Викторов и др., 1978):
Zex = JW tge(l+x<>), (3)
где W - волновое сопротивление длинной линии; в = со/Уф - коэффициент фазы; со - угловая частота; Уф -фазовая скорость распространения электромагнитных волн в вакууме; х0 - эквивалентное удлинение отрезка длинной линии, нагруженного на индуктивность LH, которое выражается следующей формулой:
Хо = (1/p)-aictg(aLJW). (4)
Из формулы (3) видно, что при А/4 < 1+х0 < А/2 входное сопротивление отрезка длинной линии приобретает емкостной характер, при этом эквивалентная электрическая емкость Ср ^ ^ при 1+х0 ^ А/2. Поэтому всегда можно выбрать такие значения длины l отрезка длинной линии или индуктивности LH, при которых Сп<<Ср, и, тем самым, обеспечить требуемое минимальное влияние Сп на частотное перекрытие (3) и, следовательно, на точность измерения физической величины во всем диапазоне ее изменения.
Отметим, что эксперименты с различными жидкостями подтверждают эффективность устройства, рассмотренного в статье (Ершов и др., 2007), для измерения различных физических свойств жидкостей. Например, в экспериментах применялась конструкция отрезка длинной линии длиной 100 мм, внутренним диаметром металлического цилиндра из нержавеющей стали 20 мм с U-образным внутренним проводником, покрытым диэлектрической (полихлорвиниловой) оболочкой. Диаметр U-образного медного проводника вместе с оболочкой равен 2 мм, а без оболочки - 1,2 мм. Части U-образного проводника расположены симметрично относительно продольной оси металлического цилиндра; расстояние между этими частями внутреннего проводника равно 7 мм. Такая же конструкция ЧЭ применима и с иным, рассмотренным выше, способом возбуждения и съема электромагнитных колебаний.
4. Расчет выходной характеристики измерительного устройства и оптимизация конструктивных параметров чувствительного элемента
Выходной характеристикой рассматриваемого измерительного устройства является зависимость резонансной частоты чувствительного элемента от измеряемого параметра, в данном случае от концентрации смеси веществ. Для получения такой зависимости, выраженной либо функционально, либо графически, прежде всего, необходимо построить теоретическую модель ЧЭ, которую можно представить в виде эквивалентной электрической схемы, изображенной на рис. 2 и состоящей из индуктивностей L1, L2 и отрезка длинной линии, имеющего волновое сопротивление W и длину l. В точке подключения L1 (или L2, в зависимости от выбора точки подключения электронного блока) полное комплексное сопротивление этой схемы Zt будет складываться из параллельного соединения индуктивности L1 и входного сопротивления отрезка длинной линии, нагруженного на оконечную индуктивность L2.
Рис. 2. Эквивалентная электрическая схема чувствительного элемента
933
Маслов А.А. и др. Радиочастотное резонаторное устройство...
Рис. 3. Поперечное сечение коаксиальной линии передачи с изолированным эксцентричным внутренним проводником.
1 - экран;
2 - изолирующая оболочка внутреннего проводника;
3 - внутренний проводник
Рис. 4. Эквивалентная схема элементарного участка длинной линии
Сопротивление ZL индуктивности L1 определяется по формуле
Zl = jaLu (5)
а для сопротивления Z1 отрезка длинной линии в точке подключения L1, нагруженного на оконечную индуктивность L2, справедлива следующая формула (Ефимов, 1964)
Z1 = W(jaL2chyl + Wshyt)l(Wchyl + ja/^shyO, (6)
где у — коэффициент распространения электромагнитной волны в отрезке длинной линии.
Учитывая (5) и (6), полное сопротивление Zt может быть рассчитано по формуле
Zt = ZZlJZ + Zl). (7)
Поскольку объектом исследования являются, как правило, жидкости с диэлектрическими потерями, то у и W в формуле (6) являются комплексными величинами и могут быть представлены в следующей форме:
гЧ(R + jaL)(G + jaC), (В)
W = V(R + jaL)/(G + jaC), (9)
где R - погонное сопротивление длинной линии, обусловленное омическими потерями в проводниках этой линии; L - погонная индуктивность длинной линии; G, C - погонные эквивалентные проводимость жидкости, заполняющей ЧЭ, и емкость длинной линии, соответственно.
Внутренний проводник длинной линии имеет изоляционное покрытие, и поэтому жидкость, заполняющая ЧЭ, является двухслойной. На рис. 3 изображено поперечное сечение длинной линии с двухслойным диэлектриком. Для упрощения теоретического рассмотрения ЧЭ на этом рисунке представлен в виде отрезка длинной линии с одним эксцентричным внутренним проводником. Длина этой линии должна вдвое превышать длину двухпроводной линии, изображенной на рис. 1. Это обусловлено требованием обеспечения эквивалентности основных параметров этих двух линий. Исходя из этого, элементарный участок длинной линии (рис. 3) можно представить в виде эквивалентной схемы (рис. 4), на которой обозначено: G,, C, - погонные проводимость и емкость области длинной линии, занимаемой изоляционным покрытием; G/, Сi - погонные проводимость и емкость области длинной линии, занимаемой контролируемой жидкостью. В этом случае комплексная погонная проводимость линии будет составлять
Gd = (Gi + jaC)(Gt + jaC/)l(Gi + G/ + jaCt + jaC). (10)
В соответствии с (Ефимов, 1964), погонные проводимости G/ и Gi можно определить из соотношений:
G/ = aC/-tg^/ и Gi = oCi-tg^i, (И)
где tg^b tgS/ - тангенс угла диэлектрических потерь в изолирующей оболочке проводника и в контролируемой жидкости, соответственно.
Величиной Gi в (10) можно пренебречь, так как проводимость изоляции в большинстве случаев намного меньше, чем проводимость измеряемой среды, и тогда можно записать
Gd = jaCG + jaC/)l(G/ + jaC, + jaC/), (12)
или, подставляя (11) в (12),
Gd =jaCil{1 + jCi l [C(j + tg£/)]}. (13)
Для составления полной математической модели ЧЭ необходимы еще формулы для Cl и Ci (Иоссель и др., 1969).
Погонная емкость конденсатора с эксцентричным внутренним проводником равна
934
Вестник МГТУ, том 13, №4/2, 2010 г.
стр.931-936
Ci = 2n£)£larch[(rM2 + rs2 - a2)l(lrwrs)\, (14)
где rw - радиус изолирующей оболочки внутреннего проводника; а - расстояние между центрами внутреннего проводника и внешней оболочки.
Погонная емкость коаксиального конденсатора, состоящего из внутреннего проводника с радиусом rc и воображаемого внешнего проводника с радиусом rw равна
Ci = 2n£o£j\n(rJrc). (15)
Погонную индуктивность L можно определить, используя формулу для волнового сопротивления линии с эксцентричным внутренним проводником и воздушным заполнением W0 (Мейнке/Гундлах, 1960)
Wo = 60arch[(r/ + rc2 - a2)l(2rsr/)\, (16)
и формулу для погонной емкости линии с эксцентричным внутренним проводником и воздушным заполнением C0, используя следующее соотношение:
Co = 2rc£/arch[(r/ + rc2 - a2)l(2rsr/)\. (17)
Используя формулы (16), (17), а также формулу (9), записанную для длинной линии без потерь, т.е. без членов Яи G, можно получить формулу для L:
L = W02C0 = 2-10"7arch[(r/ + rc2 - a2)l(2rsr/)\. (18)
Резонансная частота ЧЭ определяется по максимуму резонансной кривой или по экстремальному значению полного входного сопротивления Zt. Это входное сопротивление имеет максимальное значение при равенстве нулю его мнимой части. Также резонансную частоту можно найти, приравняв нулю производную от реальной составляющей входного сопротивления Zt. Следовательно, для определения резонансной частоты необходимо решить относительно т одно из уравнений
Im(Zt) = 0 либо d(Re(Zt))ldo = 0. (19)
Решение каждого из этих уравнений в виде функциональных зависимостей со = ON), где N -измеряемый параметр (например, концентрация бинарного раствора), представляет значительные трудности, поскольку уравнения (19) являются трансцендентными. Поэтому были использованы методы численного анализа, заложенные в программе Matlab 7.
Приведенная выше совокупность формул вполне достаточна для составления на основе Matlab 7 соответствующей программы вычисления резонансной частоты в зависимости от N. В качестве контролируемой жидкости был взят водный раствор хлористого натрия (NaCl), концентрация которого изменялась в диапазоне 0,05^5 моль!л. Значения £l и tg^ этого раствора при температуре 20 °С в зависимости от его концентрации N, выраженной в мольЛитр раствора, приведены в таблице (Ахадов, 1977). При этом значения tgO которые в этом справочнике отсутствуют, были вычислены по известной формуле (Ахадов, 1977)
tgc) = £'l£,
где £ , £ - действительная и мнимая части диэлектрической проницаемости контролируемой жидкости, соответственно.
Таблица
N 0,05 0,2 0,5 1 1,5 2 2,5 3 3,5 4 4,5 5
£ 61,9 60,7 58,4 55,2 52,4 49,8 47,5 45,2 42,9 40,5 38,5 36,6
£ 31,5 32,7 35,6 38,8 42,5 45,8 48,5 51,0 53,0 53,5 54,2 55,4
tgSl 0,51 0,54 0,61 0,70 0,81 0,92 1,02 1,13 1,24 1,32 1,41 1,51
На основе числовых данных таблицы и конкретных конструктивных параметров ЧЭ, а именно: rs = 20 мм, rc = 1,2 мм, l = 200 мм, a = 5 мм, с” помощью программы Matlab 7 было получено семейство соответствующих графиков нормированных зависимостей o(N)lo(Nmin) (рис. 5 и рис. 6), где Nmin=0,05 моль!л - минимальное значение концентрации раствора. На этих графиках величины rw и £ являются варьируемыми и изменяются в диапазоне: rw =1,6^4 мм, £ = 2,2^10.
Из рассмотрения рис. 5 и рис. 6 видно, что максимальная чувствительность достигается при rw = 2,5мм и £ = 6. Следовательно, эти параметры и являются оптимальными для данной конструкции ЧЭ. Однако, следует отметить, что физически проще варьировать значение rw, чем £, так как в последнем случае необходимо подбирать диэлектрические материалы с соответствующим значением £.
Таким образом, используя предложенную математическую модель данного ЧЭ, с помощью программы Matlab можно достаточно легко рассчитать выходные характеристики измерительного устройства для различных значений конструктивных параметров этого ЧЭ.
935
Маслов А.А. и др. Радиочастотное резонаторное устройство...
Рис. 5. Семейство выходных характеристик измерительного устройства для различных значений диэлектрической проницаемости е, изолирующей оболочки при rw=2,0 мм центрального провода
Рис. 6. Семейство выходных характеристик измерительного устройства при различных значениях радиуса изолирующей оболочки rw для е,=6
5. Заключение
Таким образом, рассмотренное радиочастотное устройство обеспечивает расширение области применения, позволяя производить измерения физических свойств жидкости при значительном удалении электронного блока от контролируемого объекта и при минимальном влиянии кабелей связи на точность измерения. Предложенная математическая модель чувствительного элемента позволяет оптимизировать его конструктивные параметры с точки зрения повышения чувствительности к измеряемым физическим свойствам жидкости и обеспечения однозначности измерений. Это устройство может быть применено для контроля концентрации смесей (растворов) жидкостей, в том числе влагосодержания жидкостей, и других физических свойств веществ.
Литература
Ахадов Я.Ю. Диэлектрические свойства бинарных растворов. М., Наука, 400 с., 1977.
Брандт А.А. Исследование диэлектриков на сверхвысоких частотах. М., Физматгиз, 404 с., 1963.
Викторов В.А., Лункин Б.В., Совлуков А.С. Высокочастотный метод измерения неэлектрических величин. М., Наука, 280 с., 1978.
Викторов В.А., Лункин Б.В., Совлуков А.С. Радиоволновые измерения параметров технологических процессов. М., Наука, 208 с., 1989.
Ершов А.М., Маслов А.А., Совлуков А.С., Фатеев В.Я., Яценко В.В. Радиочастотная
концентратометрия технологических водных растворов. Датчики и системы, № 11, c.17-21, 2007.
Ефимов И.Е. Радиочастотные линии передачи. М., Советское Радио, 600 с., 1964.
Иоссель Ю.Я., Кочанов Э.С., Струнский М.Г. Расчет электрической емкости. М., Энергия, Ленингр. отделение, 240 с., 1969.
Мейнке Х., Гундлах Ф. Радиотехнический справочник. М., Госэнергоиздат, 416 с., 1960.
936
