ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Публикационная активность

(сведения по итогам 2017 г.)
2-летний импакт-фактор РИНЦ: 0,500
2-летний импакт-фактор РИНЦ без самоцитирования: 0,405
Двухлетний импакт-фактор РИНЦ с учетом цитирования из всех
источников: 0,817
5-летний импакт-фактор РИНЦ: 0,319
5-летний импакт-фактор РИНЦ без самоцитирования: 0,264
Суммарное число цитирований журнала в РИНЦ: 6012
Пятилетний индекс Херфиндаля по цитирующим журналам: 404
Индекс Херфиндаля по организациям авторов: 338
Десятилетний индекс Хирша: 17
Место в общем рейтинге SCIENCE INDEX за 2017 год: 527
Место в рейтинге SCIENCE INDEX за 2017 год по тематике "Автоматика. Вычислительная техника": 16

Больше данных по публикационной активности нашего журнале за 2008-2017 гг. на сайте РИНЦ

Вход


Забыли пароль? / Регистрация

Добавить в закладки

Следующий номер на сайте

4
Ожидается:
16 Декабря 2018

Автоматизированная система расчета температурных полей в спаях стекловидного диэлектрика с металлом

Статья опубликована в выпуске журнала № 1 за 2002 год.[ 20.03.2002 ]
Аннотация:
Abstract:
Авторы: Малюков С.П. (malyukov@fep.tsure.ru) - Технологический институт Южного федерального университета в г. Таганроге, г. Таганрог, Россия, доктор технических наук, Бородицкий М.П. () - , ,
Ключевое слово:
Ключевое слово:
Количество просмотров: 6300
Версия для печати
Выпуск в формате PDF (1.30Мб)

Размер шрифта:       Шрифт:

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

Наиболее значимым свойством для теории и практики высокотемпературных покрытий является их термостойкость, для изучения которой требуется расчет температурных полей и напряжений, возникающих в покрытии в различных условиях службы изделий [1]. Зная температурные напряжения, можно определить допустимые для данного покрытия границы применения.

Покрытие имеет значительно меньшую толщину по сравнению с покрываемой деталью, поэтому найдем температурные поля в наиболее простом варианте – покрытие по плоскости, достаточно полно описывающее распределение температур. Рассмотрим задачу охлаждения (нагревания) полосы металла, покрытой с двух сторон стекловидным диэлектриком. Математическая постановка задачи (рис. 1):

,                                       (1)

,                                 (2)

Подпись:  
Рис. 1. Схема полосы металла с покрытием стекловидного ди-электрика
, z=h,                                                (3)

, z=-h1,                                          (4)

, z=0,                                                          (5)

, z=0,                                               (6)

, t=0.                                                    (7)

Здесь q1 и q2 – избыточные температуры в покрытии и полосе; a1 и a2 – температуропроводность покрытия и полосы; t – время нагревания (охлаждения); z – координата; h и h1 – толщина покрытия и половина толщины полосы; , – теплопроводность покрытия, x – коэффициент теплоотдачи на границе раздела покрытие–внешняя среда; b2 – теплопроводность полосы; b3 – отношение теплопроводности полосы и покрытия; q0(z) – начальная избыточная температура полосы и покрытия по отношению к внешней среде.

Для решения применяем метод разделения переменных Фурье [2] при дополнительных предположениях:  .                                              (8)

Разложив полученное решение в ряд по степеням h и ограничившись конечным числом слагаемых, количество которых определяется погрешностью начальных условий, получаем достаточно простые уравнения для нахождения собственных чисел задачи (1)-(7).

Сделаем в уравнении (1) и граничном условии (3) замену переменных .                                    (9)

Тогда уравнения (1) и (3) запишутся следующим образом:

 0 < x < 1;                              (10)

 (x = 1).                                             (11)

Методом Фурье [3] определяем частное решение уравнения (10) в виде: , тогда .     (12)

Для удовлетворения граничным условиям, подставляя (12) в (11), приходим к следующему равенству:

.       (13)

Решение задачи (2):

.  (14)

Из (5):

.    (15)

Подставляем (14) в граничное условие (4), затем (14) и (12) в граничное условие (6). Окончательно получаем следующую систему уравнений:

     (16)

В (16) учтено, что .

Анализируя систему (16), замечаем, что , то есть  и , то есть  .

Аналогично: , где через  обозначим коэффициенты разложения при  выражения  по степеням .

Подставляем в (16) разложения по степеням  функции  и  и положим в полученных равенствах . В результате приходим к системе уравнений для нахождения нулевых членов асимптотики по  функции  и  и .

                    (17)

Первое уравнение системы (17) является трансцендентным и имеет бесчисленное множество корней .

Для того чтобы удовлетворить начальным условиям (7), ищем решение задачи (1)–(7) в виде:

.

Возвращаясь к переменной z (9) и учитывая, что , то есть, что

,   (18)

окончательно получим:

,   (19)

Учтем равенство (18), тогда:

.

Подставляем последнее соотношение в выражение для :

   (20)

Учитываем начальное условие (7):

 ,

получим: ,                                  (21)

где: Un(z)=

   (22)

Во второй формуле равенства (22) сделаем замену:

.                                                              (23)

Введем единую переменную x:

(24)

Легко можно проверить, что функции  ортогональны, то есть:

, где .                 (25)

Сделаем в выражении (21) замену (23), получим:

,                          (26)

где  

Умножим ряд (26) на  и проинтегрируем в пределах . Тогда, учитывая (25), получим:

.            (27)

Из последнего равенства находим коэффициент , тем самым удовлетворено начальное условие (7) и задача (1)–(7) полностью решена.

По формулам (19), (20) можно установить температурные поля в покрытии (стекле) и металле (полосе) при различных условиях теплообмена на границе покрытия (окружающая среда).

Подпись:  
Рис. 2. Распределение температуры в покрытии и металле при охлаждении:
а, б – система А2 – Ст.3; Ci=100(а) и 10(б); в, г – система А1 – ВТ1; Ci=2(в) и 0,7(г); цифры у кривых – значения Fo, где  

Расчеты приведены для покрытий, составы которых представлены в таблице 1. Составы стекловидных диэлектриков и теплофизические параметры взяты из работы [4].

Покрытия наносили на полосы технического титана ВТ–1 и стали СТ–3 методом центрифугирования суспензий стекловидного диэлектрика. Испытания проводили сбрасыванием образцов, нагретых до определенной температуры, в резервуар с холодной водой.

При проведении пробных расчетов по формулам (19) и (20) определили, что изменение температуры в покрытии со временем определяется в основном величиной критерия, который характеризует скорость теплообмена на границе и скорость притока тепла к границе из объема покрытия:

,                                                            (28)

где h – толщина покрытия; l1 – теплопроводность покрытия; .

Таким образом, используя полученные формулы (19), (20), возможно построение автоматизированной системы расчета температурных полей в спаях стекловидного диэлектрика с металлом.

Назначение разрабатываемого алгоритма – определить толщину покрытия и материал, из которого оно изготовлено. При этом металл, на который наносится покрытие, выбирается изначально. Целевой функцией в данном алгоритме является коэффициент Ci, который задается пользователем в начале работы программы. Он зависит от толщины покрытия и от теплофизических данных материалов (таблица 2): l1 – теплопроводность покрытия; æ – коэффициент теплоотдачи на границе покрытия (внешняя среда). Коэффициент Ci также зависит от констант, которые определяются составом материала (таблица 3).

Подпись: Таблица 1
Состав исследуемых покрытий (вес. %)
Покрываемый материал	Шифр покрытия	SiO2	Al2O3	B2O3	TiO2	Li2O	Na2O	K2O	ZrO2	BaO	Прочие
Титан ВT-1	А1	65	10	5	5	8	5	-	2	-	-
	А2	60	10	-	10	8	10	-	2	-	-
	А3	60	10	10	10	5	-	-	5	-	-
	А4	38	-	6,5	-	-	-	-	-	44	11,5
	А5	15	20	30	5	-	-	-	10	-	20
	А7	60	-	15	10	-	5	10	-	-	-
Сталь Ст. 3	А9	44	5	15	10	-	15	8	-	-	3


Для поиска состава материала и толщины покрытия, которые описывает целевая функция в виде коэффициента Ci, применяем алгоритм эвристического перебора [5].

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

Подпись: Таблица 2				  Таблица 3
Теплофизические данные 				Константы материалов
Ма-териал	l1, oC
Вт/(м• oC)	æ
Вт/(м•К)		Сист. по-крытие – металл	 

ВТ1	10,88	21,9		А5 – ВТ1	0,275
Ст. 3	46,05	46		А4 – ВТ1	0,245
А1	0,71	0,058		А7 – ВТ1	0,326
А2	0,79	0,052		А9 – Ст. 3	0,20
А3	0,79	0,049	H = 5,74 • 103 Вт/(м2 • oC)
А4	0,80	0,055
А5	0,82	0,056
А7	1,01	0,047
А9	1,01	0,051

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

Для реальных покрытий (при небольших толщинах) значение критерия Ci колеблется в пределах от 0,1 до 100. При Ci < 0,1 температуру по всему сечению покрытия можно принять постоянной, то есть градиентом температуры можно пренебречь, а при Ci > 100 температура поверхности и окружающей среды выравнивается, и формулы (19) и (20) значительно упрощаются.

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

1.   Журавлев Г.И. Химия и технология термостойких неорганических покрытий. - М.: Химия. – 1975. - 198 с.

2.   Смирнов В.И. Курс высшей математики. - М.: Наука. - 1979. - Т.4. - 652 с.

3.   Олвер Ф. Введение в асимптотические методы и специальные функции. - М.: Наука - 1978. - 375 с.

4.   Малюков С.П. Стекловидные диэлектрики в производстве магнитных головок. (Монография). - Изд–во ТРТУ, 1998.

5.   Капустин Н.М., Васильев Г.Н. Автоматизация конструкторского и технологического проектирования. - М.: Высшая школа, 1986. - Кн. 6.

6.   Малюков С.П., Обжелянский С.А. Алгоритм формирования математической модели синтеза стекловидных диэлектриков для магнитных головок. // Изв. ТРТУ.- Таганрог. - №4. - 2001.


Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=721
Версия для печати
Выпуск в формате PDF (1.30Мб)
Статья опубликована в выпуске журнала № 1 за 2002 год.

Возможно, Вас заинтересуют следующие статьи схожих тематик: