CC BY
10
80
Евразийский Союз Ученых (ЕСУ) # 9 (18), 2015 | ТЕХНИЧЕСКИЕ
уголь
выработанное пространство обрушенная порода направление движения горной массы
направление движения людей
направление движения
материалов
струя свежего воздуха
струя отработанного
воздуха
шлюз
ленточный конвейер гезенк (горный бункер) скребковый конвейер кроссинг
Рис. 3. Схема обрудования для создания и поддержания искусственной атмосферы в стендах Условные обозначения
1 - нагнетательный трубопровод (подающий)
5 - станция регенерации искусственной атмосферы 9 - автоматические клапаны, управляемые
датчиками содержания в искусственной атмосфере кислорода
Рис. 4.
Схема стенда для экспериментов с животными
Рис. 5
Схема стенда для исследования жизнедеятельности человека и производственных процессов в искусственной атмосфере
Рис. 6. Схема создания искусственной атмосферы в замкнутом пространстве.
Условные обозначения
1 - нзгаггатгльныйтрубопровсд (подающий)
2 - всасывающий трубопровод для исходящих
газов
3- -шлюзовая перемычка
4 - клапан для подачи в замкнутое пространство смеси газов, составляющих искусственную атмосферу
5 - станция регенерации искусственной атмосферы
6 - убирающаяся перемычка
7 - клапан для отсоса из замкнутого пространства смеси газов
8 - вентилятор
9 - автоматические клапаны, управляемые датчиками содержания в искусственной атмосфере кислорода
10 - автоматические клапаны для стравливания атмосферного давления из искусственной атмосферы
11 - автоматические клапаны, присоединённые к адсорберу, для выделения из искусственной атмосферы кислорода и для стравливания во всасывающий трубопровод остальных газов искусственной атмосферы
12 - замкнутое пространство
14 - задвижка
15 - вакуумный насос
16 - обратные клапаны
ОПЫТ УЧЕТА НЕЛИНЕЙНОСТИ СВОЙСТВ МАТЕРИАЛОВ ПРИ РЕШЕНИИ ЗАДАЧ ИНДУКЦИОННОГО НАГРЕВА СРЕДСТВАМИ ПРОГРАММЫ FLUX
Клочкова Н. Н., Обухова А. В., Проценко А. Н.
Канд. тех наук, доценты кафедра «Электроснабжение промышленных предприятий», СамГТУ
Евразийский Союз Ученых (ЕСУ) # 9 (18), 2015 I ТЕХНИЧЕСКИЕ НАУКИ
81
АННОТАЦИЯ
В статье обсуждаются возможности, предоставляемые программным комплексом FLUX, для учета нелинейности свойств материалов исследуемых объектов. Зависимость тепловых, магнитных и электрических свойств от температуры и напряженности магнитного поля рассматриваются на примере модели индукционной установки специального назначения - азимутального одновиткового индуктора для нагрева заготовки сферической формы. Проведен ряд вычислительных экспериментов для определения влияния на результат решения задач различного вида функций, описывающих нелинейную зависимость. В качестве результатов приведено распределение температуры по радиусу и по экватору сферы при учете нелинейности различных теплофизических свойств материалов.
ABSTRACT
The article discusses the possibilities the software complexes, catfish FLUX, to account for the nonlinearity of material properties of the objects. The dependence of the thermal, magnetic and electrical properties of temperature and magnetic field are considered on the model of the induction installation of special purpose - azimuthal single-turn inductor for heating the workpiece spherical shape. A number of computational experiments to determine the effect on the result of solving problems of various kinds offunctions that describe the nonlinear relationship. As the results of the temperature distribution given by the radius at the equator of the sphere and taking into account the nonlinearity of various thermophysical properties of materials.
Ключевые слова: индукционный нагрев, многопараметрический анализа, электромагнитные и тепловые поля, метод конечных элементов
Keywords: induction heating, multivariable analysis, electro-magnetic and thermal fields, finite element method
Как известно, физические свойства материалов меняются в зависимости от состояния самого материала - в первую очередь агрегатное состояние и температура, а так же от окружающих условий -внешние магнитные и электрические поля.
Различные свойства материалов:
электропроводность, магнитные свойства,
диэлектрические или тепловые (теплопроводность, теплоемкость) характеризуются уравнениями, которые представлены в таблице 1.
Таблица 1.
Уравнения, характеризующие различные свойства материалов
Материал Уравнение Свойство материала
Магнитный В = pH р: магнитная проницаемость Г н/м
Диэлектрик D = гЁ е: диэлектрическая проницаемость Ф/м
Проводник / = аЁ ст: проводимость Ом"1м"1
Т еплопроводность ф = — [k^gradT [k]: матрица теплопроводности Вт/м/градус
Объемная теплоемкость ДС = СРАТ pCP(T): удельная объемная теплоемкость
Вследствие взаимосвязи и взаимовлияния свойств материала и условий окружающей среды, точное представление свойств и законы их изменения могут оказаться очень сложным. Как правило, невозможно выразить всю сложность законов поведения в одной модели, которая бы одновременно учитывала различные свойства. Материал можно
моделировать несколькими способами, и выбор зависит от изучаемых явлений и условий функционирования.
Flux предоставляет пользователю несколько различных моделей для описания магнитных свойств материалов (табл. 2).
Таблица
2
Модели описания магнитных свойств материалов
Модель B(H) B(H,T)
Линейная В(Ю = Рар„н (1) В(Н, Т) = РоН + Ро (рго — 1)НК(Т) (2)
С учетом насыщения НЮ = р0 Н + 2Harctg( п ^ * °н )(3) В(Н, Т) = р0Н + ^ arctg (tQlr0~1^°Н) К(Т) (4)
C учетом насыщения + регулирование перегиба НЮ = П | Г На +1—J(На +1)2—4На (1—а) "0 *s 2 (1—а) На = Ро (5) JS В(Н,Т) = До И +1„ “•+1—-l(«-2f—4«-( 1—а> к (Т) (6) на = Ро Ло
Несколько математических моделей электрических свойств представлены в таблице 3.
82
Евразийский Союз Ученых (ЕСУ) # 9 (18), 2015 | ТЕХНИЧЕСКИЕ
Таблица 3.
Математические модели описания электрических свойств материалов
Постоянное сопротивление (константа) J(E) = (1/р)Е
Линейная функция температуры p(T) = ра(l + a(T - To))
Экспоненциальная функция температуры T—To P(T) = Pa + Pbe T
Модель теплопроводности в программе Flux представлена в двух видах
Представление модели теплопроводности в программе Flux
Линейная *(T) = Яа (l + a(T — To))
Экспоненциальная T—To A.(T) = Яа + ЯЬе т
Таблица 4
Теплоемкость представлена объемной теплоемкости, отличающихся характером
теплоемкостью, зависящей от температуры Ср = температурной зависимости (табл.5).
/(Г). Возможно использование нескольких моделей
Таблица 5
Модели теплоемкости
Линейная рСрСО = рСра(1 + а(Г — То))
Экспоненциальная г-То рСДГ) = рС„а + рС„йе *
Константа + Г ауссиан 1 СТ-г^ рС„(Т) = рС„а + FGauss(T), Gauss(T) = —-=.е 2! о- ) р р а-у2л
Экспонента + Г ауссиан Т рСрСО = pCpi + (рСро - pCpj)e * + FGaussCT), 1 1гг-гс\2 Gauss(T) =—= e 2! a )
Как видно из сказанного, арсенал имеющихся в распоряжении пользователя моделей предоставляет возможность моделировать материалы с практически любым набором свойств.
Очевидно, повышенный интерес с точки зрения взаимодействия и взаимовлияния свойств различных материалов и сред представляет технология индукционного нагрева, основанная именно на воздействии магнитного поля на металл, в результате которого изменяются физические свойства металла -закалка или нагрев перед деформацией.
Использование различных моделей свойств материалов рассмотрено на примере физического описания геометрической модели индукционной установки специального назначения рассмотренной в [1, 2].
Физическое описание модели начинается с образования физического приложения [3]. Рассматриваемая задача индукционной закалки требует образования физического приложения типа 3D Steady State AC Magnetic coupled with Transient Thermal application которое позволяет совместно решать стационарную электромагнитную и нестационарную тепловую задачи.
Расчетная область представляется в виде двух вложенных параллелепипедов, внутренний объем представляет основную расчетную область, а внешний объем зону перехода в бесконечность (рис. 1).
Благодаря этому удовлетворяются
противоречивые требования - моделирование бесконечности с помощью ограниченного объема.
Построение расчетной области начинается с помощью команды Geometry/Infinite box/New и завершается процедурой построения поверхностей и объемов расчетной области.
Рис. 1. Внешний вид «бесконечного» параллелепипеда.
Следующий этап физического описания геометрической модели - создание объемных областей (Volume region), представляющих определенные части индукционной установки. При создании областей им назначаются тепловая и магнитная характеристики (табл.6), а так же указывается материал, из которого сделана часть установки, представленная областью. В рассматриваемой задаче созданы следующие области: Air - расчетная область, Copper - медь трубки индуктора, Core - нагреваемая заготовка, Slip -охватывающий сфероид [1], Water - внутренний объем трубки индуктора.
Созданные области необходимо связать с
Евразийский Союз Ученых (ЕСУ) # 9 (18), 2015 I ТЕХНИЧЕСКИЕ НАУКИ
83
соответствующими объемами геометрической модели с помощью команды Physics/Assign region to geometric entities/Assign region to volumes (completion mode).
Для построения моделей содержащих проводники, присоединенные к источнику тока, FLUX предоставляет возможность объединения уравнений внешней электрической цепи и уравнений электромагнитного поля. В этом случае решается совместная система уравнений [2].
Внешняя электрическая цепь строится с помощью редактора схем, запускаемого по команде Physics/Circuit/Circuit editor context.
Параметры элементов электрической цепи назначаются в редакторе цепи или в дереве данных в основном окне проекта. В данном случае назначены следующие значения: V=100B, C=25*10-<^;
Следующий шаг - описание механического аспекта модели по команде Physics/Mechanical set/New. Используя эту команду, образуем две части механической системы: подвижная часть (Mobile) -сфера с охватывающим сфероидом и неподвижная (Fixed) - остальные элементы системы. Подвижная часть имеет параметры, показанные на рисунке 2.
а) б)
Рис. 2. Параметры подвижной части механической системы: а-характер движения, вращение вокруг
оси OZ; б - скорость вращения 100об/мин
После образования механической системы необходимо связать ее части с соответствующими областями (regions) геометрической модели.
Принадлежность геометрической области определенной части механической системы
указывается в окне Edit volume region в закладке Mechanical Set (рис. 3).___________
^ Edit Volume Region[COREJ
Name of the region*_____________________________
I
Comment_________________________________________
Magneto thermal \ Appearance ' Mechanical Set \ Evaluated information \ I
Mechanical set of the region___________________
[mobile Ч1Л
Рис. 3. Диалоговое окно: назначения механических свойств области.
Для окончательного формирования
приложения (application) необходимо описать
сценарий решения. Сценарий имеет следующие параметры:
- частота тока 20кГц;
- контроль переходного процесса по положению механической системы;
- пределы: 0°, 90°;
- шаг 30°.
В ходе решения ряда задач исследовалось влияние учета нелинейности свойств материалов с использованием различных моделей на результат решения. Оценивалось также время решения.
Решение проведено на компьютере с процессором Intel® Core™ i5 CPU и оперативной памятью 16 Гб.
В качестве результатов рассмотрено распределение температуры по радиусу сферы, а также по экватору.Для сравнения, был проведен расчет нескольких задач с учетом нелинейности различных свойств стали.
Результаты сравнения приведены на рис. 4 и в
табл. 7
Таблица 7
Сравнительные результаты решения задач с различной нелинейностью.
Свойство материала Модель нелинейности Время решения сек
Электросопротивление Линейная 6824
Экспонента 7352
Т еплопроводность Линейная 7778
Экспонента 9787
Теплоемкость Линейная 778
Экспонента 18774
Все свойства Линейная 1784
Экспонента 10793
Зависимость В(Н) (5) табл. 2. 2884359 (33,4 суток)
84
Евразийский Союз Ученых (ЕСУ) # 9 (18), 2015 | ТЕХНИЧЕСКИЕ
Приведенные результаты показывают, что учет различных нелинейностей влияет очень по разному и может потребовать очень больших затрат времени, поэтому полный учет нелинейностей целесообразно проводить только на заключительных этапах вычислительных экспериментов.
Из рисунка также видно, что наибольшее влияние на результат оказывает нелинейность магнитных свойств, но решение таких задач требует
чрезвычайно больших затрат времени. В задачах проектирования индукционных установок и режимов их работы большой интерес представляют именно электромагнитные явления и их нелинейность. Поэтому, для решения таких задач с помощью FLUX за приемлемое время необходимы дополнительные меры. Самый простой путь - увеличение вычислительных мощностей.
Рис. 7. Распределение температуры по экватору и радиусу сферы в зависимости от учитываемой нелинейности: а) - нелинейная зависимость В(Н); b) - нелинейность других свойств (табл. 7).
Список литературы
1. Клочкова Н.Н, Обухова А.В., Проценко А.Н. Проектирование индукционной нагревательной установки для тел сферической формы средствами программного пакета Flux. // IX Международная научно-практическая конференция. Научная
публикация: Отечественная наука в эпоху изменений: постулаты прошлого и теории нового времени. Россия
Екатеринбург 16-17 мая 2015., часть 2, стр 151-154
2. Клочкова Н.Н, Обухова А.В., Проценко А.Н. Моделирование индукционной установки специального назначения средствами программного пакета FLUX. // Вестник Саратовского ТГУ Вестник СГТУ. 2015. №2 (79).
3. http://www.cedrat.com
ГИДРОАКУСТИЧЕСКИЙ МОДЕМ ДЛЯ АВТОНОМНОГО НЕОБИТАЕМОГО
ПОДВОДНОГО АППАРАТА
Попко Артем Олегович
ведущий инженер лаборатории автоматизации АНПА ОАО «Концерн «НПО «Аврора», кандидат военных наук.
Ключевые слова: автономный необитаемый подводный аппарат, гидроакустический модем,
гидроакустическая связь Аннотация.
В статье автор освещает вопросы создания гидроакустического модема для нужд телеметрического контроля и управления АНПА.
Keywords: autonomous underwater vehicle, hydroacoustic link modem, sonar communication.
Abstracts.
In the topic author wrote about realization of hydroacoustic modem for remotely controlled AUV.
Аппаратная часть гидроакустического модема состоит из двух модулей: корабельного и
автономного. Автономный модуль собирает информацию с АНПА и передает ее корабельному модулю. Между двумя модулями реализуется полудуплексная связь на несущей частоте 25 кГц в
полосе 5 кГц, причем нисходящий поток (в сторону корабельного модуля) является основным и высокоскоростным (для передачи данных), а восходящий поток - низкоскоростным (для управления). Вариант исполнения аппаратной части модема представлен на рисунке 1 .
