Научная статья на тему 'Проект аналитической системы комплексной оценки экологической безопасности промышленного предприятия'

Проект аналитической системы комплексной оценки экологической безопасности промышленного предприятия Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
67
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭКОЛОГИЧЕСКАЯ БЕЗОПАСНОСТЬ / РАЗРАБОТКА ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ ДЛЯ ПРЕДПРИЯТИЙ / ENVIRONMENTAL SAFETY / SOFTWARE DEVELOPMENT FOR ENTERPRISES

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Лизогуб Анастасия Андреевна, Ляхов Герман Геннадьевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Лизогуб Анастасия Андреевна, Ляхов Герман Геннадьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

On an analytical system for the integrated assessment of the environmental safety in an industrial enterprise

Environmental protection is a global problem many research works have been devoted to finding ways of reducing the anthropogenic load on the environment. This article describes an algorithm of an analytical system for the integrated assessment of the environmental safety of an industrial enterprise.

Текст научной работы на тему «Проект аналитической системы комплексной оценки экологической безопасности промышленного предприятия»

УДК 53.083; 621.317.3

М. А. Никишин, А. В. Румянцев, В. Г. Харюков

БЕСКОНТАКТНЫЙ МЕТОД ИЗМЕРЕНИЯ УДЕЛЬНОГО ЭЛЕКТРОСОПРОТИВЛЕНИЯ ЖИДКИХ МЕТАЛЛОВ В ОБЛАСТИ ВЫСОКИХ ТЕМПЕРАТУР

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

The article offers a new, non-contact method of measuring the specific electrical resistance of liquid metals in the high-temperature range, which is exceptionally simple. The essence of the method is the measurement of the electromotive force induced by a high-frequency electromagnetic field in two circuits of different diameters placed coaxially in the middle of a cylindrical sample, which is a thin-walled ampoule filled with some low-melting metal. The authors substantiate the feasibility of a contactless method for measuring the specific electric conductivity of both solid and liquid metals in the high-temperature range. Result-processing methodology is given as well as the recommendations on the optimal geometry of the contours and the sample for increasing the measurement accuracy.

Ключевые слова: температура, индукционный генератор, электромагнитное поле, скин-слой, электродвижущая сила, удельное электрическое сопротивление, мощность.

Keywords: temperature, induction generator, electromagnetic field, skin-layer, electromotive force, specific electric resistance, power.

Главный недостаток контактного метода измерения удельного электросопротивления металлов, базирующегося на законе Ома, ясен из его названия — контактность. Сложность осуществления контакта образца с проводами — токоподводящими и потенциальными выводами — обусловлена несколькими причинами. Первая из них: если при применении металлов, поддающихся механической обработке, можно вставлять провода в заранее высверленные отверстия с последующей их зачекан-кой, то для металлов типа Ta, W, Re и др., механической обработке не поддающихся, а также в случае исследования жидкого металла, помещенного в тонкостенную ампулу, контакт можно осуществить только

© Никишин М. А., Румянцев А. В., Харюков В. Г., 2018

Вестник Балтийского федерального университета им. И. Канта.

Сер.: Физико-математические и технические науки. 2018. № 3. С. 90—96.

электронно-лучевой сваркой. При отсутствии в лаборатории такой установки возникают трудности, связанные с поисками стороннего исполнителя. Вторая причина: материал образца и проводов — токопод-водящих и потенциальных (диаметрами < 1,0 мм и < 0,1 мм соответственно) — должен быть одним и тем же, что далеко не всегда осуществимо. Поэтому проблема создания бесконтактного (безэлектродного) метода измерения удельного электросопротивления металлов представляется достаточно актуальной, особенно в области высоких температур.

Бесконтактные методы базируются на использовании электромагнитного поля [1 — 3]. В работе [1] предложен метод дифференциального трансформатора; в [2; 3] метод состоит в определении угла поворота твердого образца (или ампулы сферической формы с жидким металлом) под воздействием вращающегося магнитного поля. Эти методы, в отличие от контактного, обременены большим количеством измерительной аппаратуры (особенно в [2; 3]) и довольно трудоемкой процедурой обработки экспериментальных результатов, не говоря уже о сложности создания вращающегося магнитного поля. Поэтому эти методы едва ли можно рассматривать как предпочтительные в сравнении с контактным методом.

Использование электромагнитного поля, создаваемого высокочастотным индукционным генератором (частота / ~ 0,5 мГц), позволяет предложить уникальный по простоте бесконтактный метод измерения удельного электросопротивления металлов. Суть метода — измерение электродвижущей силы (ЭДС), наводимой высокочастотным электромагнитным полем в двух круговых контурах диаметрами йи > йк2 > й (й — диаметр образца), которые размещены коаксиально с цилиндрическим образцом в плоскостях, перпендикулярных к его оси и расположенных посередине образца. Образец может представлять собой либо однослойный, либо двухслойный металлический цилиндр в виде тонкостенной ампулы с помещенным в нее легкоплавким металлом, удельное электросопротивление которого нужно определить. Измеренные значения ЭДС позволяют рассчитать удельное электросопротивление р исследуемого: в первом случае — тугоплавкого твердого металла, во втором — легкоплавкого жидкого металла. Диапазон температуры — от точки Кюри до максимальной, ограниченной мощностью индукционного генератора. Измерительная аппаратура: вольтметр, к которому поочередно подключаются выводы контуров; оптический пирометр. Иллюстративная схема установки представлена на рисунке.

Для реализации предлагаемого метода необходимо найти связь между наводимыми в контурах ЭДС. Для этого воспользуемся очевидным с физической точки зрения положением — вводимая в образец на единицу его длины мощность зависит только от напряженности создаваемого индуктором магнитного поля и электрических свойств материалов (тугоплавкого металла либо ампулы и металла, ее заполняющего). Тогда, приравнивая формулы для мощности, содержащие величину наводимых на одном витке круговых контуров ЭДС е\ и е2, можно установить искомую связь между е1 и е2.

91

92

Рис. Иллюстративная схема установки (индуктор и детали крепления не показаны): 1 — ампула с жидким металлом (диаметр 10 — 12 мм; длина 80—100 мм); 2 — круговой контур диаметром йк2, количество витков — Ш2; 3 — круговой контур диаметром йц, количество витков — Ш1; (йк1 > йк2; Ш2 > Шг); К — двухполюсный сдвоенный переключатель; V — вольтметр; угол между выводами контуров — 90°

Задача по нахождению формулы для мощности на единицу длины в ее общем виде — для двухслойного цилиндра, содержащего среды с разными р — решена в работе [4, с. 38—47]. Там же, как частное, приведено решение и для однослойного цилиндра (достаточно положить рг = р2 и йг= й2).

При решении задачи считались выполненными следующие предпосылки:

1) создаваемое индуктором высокочастотного генератора электромагнитное поле однородно по длине образца и по радиусу вне образца. Это условие эквивалентно предположению о больших в сравнении с диаметрами длинах образца и индуктора;

2) зависимость компонент электромагнитного поля от времени является гармонической;

3) температурные градиенты отсутствуют — образец нагрет равномерно.

Граничные условия задачи: равенство напряженности полей на границах раздела.

В [4] приводятся формулы для мощностей, вводимых в каждую из сред: W1 — в жидкий металл, — в ампулу.

В общем виде результат приравнивая формул для мощности дает следующие выражения в случае использования однослойного и двухслойного — W1 — цилиндров:

42 _ Пл,

¥( х, Ь, у, N1) ¥( х, Ь, у, N 2)'

(1)

(2)

где е = - йФ/йЬ связывает электродвижущую силу, наведенную на одном витке контура с магнитным потоком Ф = В-Б, пронизывающим площадь

Б = пй2(Ык2 -1)/4; е — электродвижущая сила, наведенная на одном витке кругового контура, В; п = о/й, где о = (пцоц//р)- 1/2 — эффективная толщина скин-слоя в среде; ц0 = 4п-10- 7 — магнитная постоянная, Гн/м; ц = 1 — магнитная проницаемость сред; / — линейная частота генератора, Гц; х = (р1/р2)1/2 — отношение удельных электросопротивлений жидкого металла и материала ампулы; в = (й2 - й1)/й2 — геометрический параметр; Мы = (йк1 /й2) — отношение диаметра контура к внешнему диаметру образца (1 = 1; 2).

Формула (1) для однослойного цилиндра содержит две неизвестные: е1/е2 и п, содержащую р металла. Если измерить ЭДС е1 и е2, то можно сразу найти величину п, то есть р исследуемого металла. Первый вариант реализован нами экспериментально для тугоплавкого металла — молибдена. Остановимся на рассмотрении второго варианта.

Формула (2) содержит три неизвестных: е\/еъ х = (рг/р2)1/2 и П1, также содержащее неизвестное р1. Первое неизвестное может быть определено, как и в первом варианте, экспериментально, а избавиться от неизвестного п1 можно переходом к п2 согласно связи между ними: П1 = хп2/(1 - в) . Знаменатель в формуле (2) отличается от числителя только тем, что в нем фигурирует Ык2, а не Ми. Представим формулу (2) в виде, удобном для вычислений:

где

[^(х, 3, У)+(К-1)о2(х, 3, п]2 +[^(х, 3,п)+(N1 -Щх, 3, п)]2 ^(х, 3, п)+(N¡2 1)^2 (х, 3, ^)12 +[^(х, 3, п)+(N2 -%(х, 3, /

(х, 3, п) = т(т-х33) / + 2(1 - 3)/4 + 3 (2 + 3), 2(1 3)

(3)

3 (2 + 3 )х2

4(1-3)

а2(х,3,п)= /+1+3+32,

М х, 3, п) = ^(3-3) /2 - 2(1-3) /з,

Ь2(х,3,п) = /2 + ^1Т

2 4(1-3) 2 4(1-3) п

В выражении (4) функции Дх, в, п) имеют вид (1 = 1; 2; 3; 4)

/(х, 3, п) = ^---1,

1 хп 4(1-3)

/,(х, 3, п) = ^ + хп ■ х2п 2

(4)

хп 4(1- 3) 2(1- 3) 2 2 (5)

/з( х, 3, п)=я х, 3, п),

2 2

/4 (х, 3, п) = -Г-3-2 /2 (х, 3, п). (1 3)

93

Для сокращения записи в формулах (3), (4) и (5) индекс у щ опущен и в формуле (4) у функций /1 не приведены аргументы ^(х, в, п) = Л

Процедура вычислений состоит в следующем: при известных величинах ц(й2, р2, А) и в находятся согласно выражения (5) для Дх), затем вставляются в (4), что позволяет найти я^х) и ^(х), которые вносятся в формулу (3). В результате в числителе и знаменателе будут стоять уравнения, содержащие хп с п = 0; 1; 2; 3; 4, то есть уравнения 4-й степени. Как показали конкретные вычисления, членами с х в 3-й и 4-й степенях можно пренебречь вследствие малости величин коэффициентов, на которые они умножаются, в сравнении с коэффициентами при дру-94 гих степенях х. В конечном итоге формула (3) будет иметь вид

- 2 2

_ А0 + Д X + Д * (6)

В0 + В1х + В2 х2'

Формула (6) позволяет по экспериментально найденной величине е1/е2 найти значение х^) путем решения квадратного уравнения и, тем самым, значение р1(Т) = х2(Т)- р2(Т), соответствующее измеренной температуре Т{. Таким образом, второй вариант становится в экспериментальном отношении идентичным первому, но величина х определяется посредством решения уравнения (6) для каждой температуры Т, что усложняет процедуру нахождения функциональной зависимости р1(Т) во всем доступном измерению диапазоне температур.

Для получения формулы (1) для сплошного цилиндра в развернутом виде достаточно в (4) — (5) положить в = 0 и х = 1:

Р(V,Ми) = (N2 -1 + 2/4)2 + 4Ц. (7)

Подстановка (7) в (1) позволяет по экспериментально найденной величине е1/е2 сразу вычислить значение п(Т и, тем самым, искомую величину р(Т).

Диаметры контуров следует задавать, исходя из следующих позиций. Наличие образца в магнитном поле индуктора приводит к радиальной неоднородности магнитного поля вблизи образца [5]. Минимальное расстояние от образца, когда влиянием неоднородности поля можно пренебречь, составляет Я > 1,3 Яобраща. Влияние отдельных витков индуктора начинает сказываться на расстоянии Я < Яиндуктора - 5/1, где к — зазор между витками индуктора. Таким образом, для минимального и максимального радиусов витков круговых контуров имеем выражение:

1,3Яобр азца < Як < Яиндуктора 5/1. (8)

Отметим, что предлагаемый метод не имеет ограничений как на диапазон температур сверху, так и на величину удельного электросопротивления материалов образца. Минимальная температура, при которой может использоваться метод, — температура точки Кюри для сред, то есть когда ц = Ц2 = 1.

Для увеличения ЭДС на керамический каркас катушки следует намотать несколько витков проволоки. Количество витков — т\ и т2, при этом желательно иметь отношение т2/ т1 ~ 2, что обеспечит примерно одинаковые по величине показания вольтметра: И2(йк2, т2) и ¡¿(¿и, т1). Отношение ЭДС, наведенных на одном витке, находится по формуле

ех/ег = (т2/та) • (и/^). (9)

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

1) вместо одной — две среды с разными р1 (жидкого металла) и р2 (материала ампулы), которые должны быть известны;

2) магнитное поле должно проникать в первую среду, что налагает требование на толщину стенки ампулы: АЯ/о2 << 1, что равносильно требованию АЯ ^ 0 или в ^ 0;

3) ввиду отсутствия модели черного тела в эксперименте может измеряться только яркостная температура Т8. Это означает, что для нахождения истинной температуры нужно предварительно экспериментально найти монохроматическую степень черноты внешней поверхности материала ампулы £\Т(Т).

Чтобы убедиться в возможности реализации предлагаемого метода в случае использования жидкого металла, рассмотрим конкретный пример.

Зададим следующую геометрию системы: й2 = 10 мм; в = 0,02; йк\ = 22 мм; йк2 = 18 мм; Мк1 = йкг/йч = 2,2; Мк2 = йк2/й2 = 1,8. При Т = 1700 К и / = 400 кГц величина п для ниобиевой ампулы будет равна п = П2 = 0,06498 (рыь = = 66,8-10-6 Ом-см). В качестве жидкого металла возьмем олово с рэп = = 75,6-10-6 Ом-см [6], тогда х = 1,06381. Подставляя эти данные в соответствующие формулы, получим в итоге £г(Т)/£2(Т) = 1,67223.

Возможна более простая оценка. Подставляя в формулу(7) для сплошного цилиндра пэп = 0,0707 или 0,06923 (при й2 = 9,8 или 10 мм), найдем £-[(Т)/е2(Т) = 1,66176 и 1,67061 соответственно. Отличие от точного решения (6) будет 0,6% и 0,1%, что свидетельствует о действительно малом влиянии материала ампулы на результаты расчета. Отметим, что при той же геометрии системы для молибдена в виде сплошного цилиндра экспериментальная величина £г(Т)/е2(Т) = 1,67974, откуда следует правильность полученных оценок для £1(Т1)/е2(Т1) при применении тонкостенной ампулы с жидким металлом. Тем не менее реальную, а не теоретическую работоспособность предлагаемого бесконтактного метода для измерения удельного электросопротивления жидкого металла следует подтвердить экспериментально.

С целью повышения точности метода рекомендуется: брать отношение Ми/Мк2 = 1,4, а не 1,22, как в данной работе; каждую катушку делать не односекционной, а двухсекционной, что даст три значения е1, е2 и е3 (на каждой из секций и их сумме) с последующим их усреднением.

95

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

1. Батыгин Ю. В., Гнатов А. В., Барбашова М. В. и др. Бесконтактный способ измерения удельной электропроводности листовых металлов // Электротех-шка i электромехашка. 2012. № 1. С. 69—73.

2. Способ и устройство для бесконтактного измерения удельного электрического сопротивления металлического сплава методом вращающегося магнитного поля : пат. 2531056 Рос. Федерация. № 201312919l/28 ; заявл. 25.06.2013 ; опубл. 20.10.2014, Бюл. № 29.

3. Рябина А. В., Кононенко В. И., Ражабов А. А. Безэлектродный метод измере-- ния электросопротивления металлов в твердом и жидком состояниях и уста-

96 новка для его реализации // Расплавы. 2009. № 1. С. 36—42.

4. Румянцев А. В. Вопросы сушки и полимеризации (источники нагрева). Калининград, 1975.

5. Бабат Т. М. Индукционный нагрев металлов и его промышленное применение. М. ; Л., 1965.

6. Макаренко И. Н., Труханова Л. Н., Филиппов Л. П. Тепловые свойства молибдена при высоких температурах // Теплофизика высоких температур. 1970. Т. 8, № 2. С. 445 — 447.

Об авторах

Максим Александрович Никишин — асп., Балтийский федеральный университет им. И. Канта, Россия.

E-mail: nikishin.maxim.a@gmail.com

Альберт Владимирович Румянцев — канд. физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Россия.

E-mail: albert37@list.ru

Владимир Геннадьевич Харюков — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Россия.

E-mail: kharyukov@mail.ru

The authors

Maxim A. Nikishin, PhD Student, I. Kant Baltic Federal University, Russia. E-mail: nikishin.maxim.a@gmail.com

Dr Albert V. Rumyantsev, Professor, I. Kant Baltic Federal University, Russia. E-mail: albert37@list.ru

Dr Vladimir G. Kharyukov, Associate Professor, I. Kant Baltic Federal University, Russia.

E-mail: kharyukov@mail.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.