CC BY
8
№ 3 (108)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2023 г.
ТЕХНОЛОГИЯ МАТЕРИАЛОВ И ИЗДЕЛИЙ ТЕКСТИЛЬНОЙ И ЛЕГКОЙ ПРОМЫШЛЕННОСТИ
DOI - 10.32743/UniTech.2023.108.3.15122
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КРАСКОПЕРЕХОДА ГОФРИРОВАННОГО КАРТОНА МЕТОДОМ ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА
Ешбаева Улбосин Жамаловна
д-р техн.х наук, проф. Намангаского инженерно- технологического института, Республика Узбекистан, г. Наманган E-mail: Guli- 6 7@mail.ru
Нишанов Акбаржон Мухаматджанович
д-р филос. по техн. наукам, PhD Намангаского инженерно- технологического института Республика Узбекистан, г. Наманган
Джалилов Анвар Абдугафарович
д-р филос. по техн. наукам, доц. Ташкентского института текстильной и легкой промышленности,
Республика Узбекистан, г. Ташкент, E-mail: anvar-matbaa@mail. ru
MATHEMATICAL SIMULATION OF THE INK TRANSITION OF THE CORPORATED CARDBOARD BY THE EXPERIMENT PLANNING METHOD
Ulbosin Eshbaeva
Doctor of Technical Sciences, Associate Professor of the Namangan engineering-technological institute, Republic of Uzbekistan, Namangan
Akbarjon Nishanov
Assistent Professor of Namangan engineering-technological institute Uzbekistan, Republic of Uzbekistan, Namangan
Anvar Djalilov
Assistent Professor of Tashkent Institute of Textile and Light Industry Republic of Uzbekistan, Tashkent
АННОТАЦИЯ
В статье создана математическая модель в виде уравнения, связывающего параметр оптимизации с факторами для планирования эксперимента. Было проанализировано влияние скорости печати и гладкости материала на краскопереход. Доказано, что максимальный краскоперенос достигается при скорости флексографского печатного оборудования 160 отт./час и гладкости бумаги 70 с.
ABSTRACT
The article created a mathematical model in the form of an equation that relates the optimization parameter to the factors for planning the experiment. The effect of print speed and material smoothness on the ink transition was analyzed. It has been proven that maximum ink transfer is achieved with a flexographic press speed of 160 prints per hour and paper smoothness of 70 seconds.
Библиографическое описание: Ешбаева У.Ж., Джалилов А.А., Нишанов А.М. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КРАСКОПЕРЕХОДА ГОФРИРОВАННОГО КАРТОНА МЕТОДОМ ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА // Universum: технические науки : электрон. научн. журн. 2023. 3(108). URL: https:// 7un iversum. com/ru/tech/arch ive/item/15122
№ 3 (108)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2023 г.
Ключевые слова: математическая модель, параметр оптимизации, полином, гофрированный картон, флексографская печать.
Keywords: mathematical model, optimization parameter, polynomial, corrugated cardboard, flexographic printing.
Введение. Целью планирования эксперимента является создание математической модели в виде уравнения, связывающего параметр оптимизации с факторами. Такое уравнение называют также функцией отклика [1]. В общем виде функция отклика, являющаяся и параметром оптимизации (у), может быть представлена выражением
y = fxi, Х2, Xk),
(1)
где XI, Х2, ..., Хк- независимые переменные факторы.
Наиболее простой моделью является полином. Полином линеен относительно неизвестных коэффициентов, что упрощает обработку наблюдений. Полином может быть первой, второй и более высокой степени. Коэффициенты полинома вычисляют по результатам опытов.
Экспериментальное исследование. На первом этапе планирования - определении направлени движения к оптимуму и крутого восхождения по поверхности отклика - налиболее целесообразно неизвестную функцию отклика аппроксимировать полиномом первой степени. Полином первой степени в общем виде выражается уравнением
У = Ьо + ЪХ1 + Ь2Х2 + ... + ЬкХк + Ь]2Х]Х2 + Ъ13Х1Х3 + ...
... + Ь12...к Х1Х2...Хк (2)
Для двух факторов это уравнение имеет вид:
У = Ьо + Ь1Х1 + Ь2Х2 + Ь12 Х1Х2 (3)
Для удобства записи условий эксперимента и обработки экспериментальных данных уровни факторов кодируют. Кодированное значение фактора XI определяют по выражению:
* = ^ (4)
где хь - натуральное значение /-фактора;
х0 - натуральное значение основного уровня /-фактора;
— интервал варьирования /-фактора. В кодированном виде верхний уровень обозначают + 1, нижный -1, а основной 0. Число N всех сочетаний уровней факторов, и следовательно, и число опытов в полном факторном эксперименте, определяется выражением
N = mk
(5)
где т - число уровней каждого фактора; к - число факторов.
Полный факторный эксперимент позволяет количественно оценить линейные эффекты и все эффекты взаимодействия. Для полного факторного эксперимента типа 22 уравнение регрессии с учетом эффектом взаимодействия представляется выражением (3).
Проведем математическое планирование и оптимизацию краскоперехода гофрированного картона [2-4]. В качестве параметра оптимизации (у) примем краскопереход, а входными факторами: Х1 - скорость печати (отт/час), Х2 - гладкость бумаги (с) (табл. 1).
Таблица 1.
Уровни и интервалы варьирования факторов
с
Фaктoры Кодовое обозначение Интервалы варьирования Уровни факторов
верхний +1 основной 0 нижний -1
Скорость, отт./час Х1 40 160 120 80
Гладкость бумаги, с Х2 25 70 45 20
В табл. 2 даны матрицы планирования и результаты опытов.
Таблица 2.
Матрица планирования и результаты опытов
Номер опыта Х0 Х1 Х2 Х1Х2 у
1 + - - + 29
2 + + - - 43
3 + - + - 46
4 + + + + 51
№ 3 (108)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2023 г.
Рассмотрим методику обработки эксперимента при отсутствии дублирования [5-6]. Обработку результатов эксперимента в этом случае произведем по следующему алгоритму.
1. Вычисление дисперсии Бу воспроизводимости эксперимента. Для этого необходимо выполнить несколько параллельных опытов в нулевой точке (в центре плана) и вычислить дисперсию Б у воспроизводимости эксперимента (табл. 3).
S2 = ¿у
ЪИ^и-уП
п0-1
(6)
где по - число параллельных опытов в нулевой точке; Уи - значение параметра оптимизации в и-м опыте; у - среднее арифметическое значение параметра оптимизации в по параллельных опытах.
Таблица 3.
Вспомогательная таблица для расчета дисперсии Б2 параметра оптимизации у
1
Номер опыта в центре плана yu У yu - У (yu - У)2 S2 Sy
1 45 0,7 0,49
2 46 Yl=iyu 1,7 2,89 И=1(Уи-^2 8,67 ¿*и=ыи yj = > = 4,335 п0 — 1 3-1
3 42 3 133 =—=443 -2,3 5,29
3 ^уи = 133 и=1 3 ^(У«-У)2 = 8,67 и=1
Примечание: где по - число опытов в центре плана; уи - значение параметра оптимизации в цетре плана.
Таблица 4.
Вспомогательная таблица для расчета дис^рсии адекватности Б2
ad
Номер опыта yj yj yj - у (yj - yj)2
1 29 31,25 -2,25 5,0625
2 43 40,75 2,25 5,0625
3 46 43,75 2,25 5,0625
4 51 53,25 -2,25 5,0625
у - среднее арифметическое значение параметра оптимизации в по параллельных опытах.
В табл. 3 даны результаты расчета Б2 =9,335. 2. Вычисление коэффициентов модели. Свободный член Ьо определяют по формуле:
Ьо=$ ъи*
(7)
Коэффициенты регрессии, характеризующие линейные эффекты, вычисляют по выражению
bi=1YIi=iXijyj
(8)
Расчеты по формулам (7), (8) и (9) дали следующие значения коэффициентов:
Ьо = 42,25; Ь = 4,75; Ь2 = 6,25; Ьп = -2,25.
3. Проверка статистической значимости коэффициентов уравнении регрессии. Проверку значимости коэффициентов произведем способом сравнения абсолютной величины коэффициента с доверительным интервалом.
Вычислим дисперсию коэффициентов регрессии
Коэффициенты регрессии, характеризующие эффекты взаимодействия, определяют по формуле
ЬЦ = iJY'j=iXijxljyj
(9)
где /, I - номера факторов;
j - номера строки или опыта в матрице планирования;
Уj - значение параметра оптимизации в у-м опыте; Xij, ху - кодированные значения (±1) факторов I и I в 7-м опыте.
s2{bt} = ^s2
(10)
где Б2{Ь1} - регрессия ^коэффицентининг диспер-сияси;
N - режалаштириш матрицасида ^аторлар ёки синовлар сони.
Б2[Ь{} = 1 • 4.335 = 1.084; Б{Ь1} = 1.041 4
Из формулы (10) следует, что дисперсии всех коэффициентов равны.
№ 3 (108)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2023 г.
Доверительный интервал ДЬ/ находят по формуле:
ДЬ = ±гт S{Ь/} (11)
где ^ - табличное значение / - критерия для принятого уровня значимости (5% ном уровне) и числе степеней свободы / которое определяют по выражению j=n0-1=3-1=2.
С учетом приведенного ^ = 4,3, тогда
ДЬ/ = ±4,3 х 1,041 = ±4,48
Таким образом, |Ь1 > |ДЬ/|; |Ь2| > |ДЬ/|; |Ь12| < |ДЬ±|. Поэтому с учетом статистической чзанимости коэффициентов Ь1 и Ь2 получим модель в виде полинома первой степени:
тогда Ет = 18,5 [1], а расчетное значение F-критерия Фишера равно
20,25
Fp =-= 4,67
р 4,335
Так как Ер > Ет, то модель, представленная уравнением (12) адекватна.
Таким образом, получена математическая модель краскоперехода гофрированного картона в виде уравнения (12).
Перейдем от кодированных Х1 и Х2 значений факторов к натуральным, если обозначим через V -скорость печати оттиска, g - гладкость бумаги. Для этого запишем кодированные значения факторов через натуральные в соответствии с зависимостями:
y = 42,25 + 4,75xi + 6,25x2
(12)
У-Уо E1
*2
з-з о
E2
(15)
4. Определение дисперсии Б2 адекватности по
формуле (табл. 4):
2 =?j=i(yj-9j)2 _Zj=i(yj-9j)2
f
N-(k+1)
(13)
где yj - наблюденные значение параметра оптимизации в j-опыте;
у, - значение параметра оптимизации, вычисленное по модели для условий j-ного опыта;
j - числов степеней свободы, которое для линейной модели определяется по выражению j = N -(к+1), где к - число факторов.
5. Проверка гипотезы адекватности модели по Б-критерию Фишера:
F = -ад
(14)
Если расчетное значение Ер меьше табличного Ет при принятом уровне значимости и соответствующих чисел степеней свободы, то модель считают адекватной. При Ер > Ет гипотеза адекватности отвергается.
Степени свободы для и соответственно равны:
j = N - (к+1) = 4 - (2 + 1) = 1; f= по - 1 = 3 - 1) = 2
где £1 и е2 - интервалы варьирования соответственно факторов х1 и х2;
vo, go - основные уровни факторов в натуральных выражениях.
_у — 120 _д — 45
Хл — ; —
40
25
Тогда получим зависимость краскоперехода (К) от скорости оттиска (V) и гладкости бумаги в виде:
К = 42,25 + 4,75 —+ 6,25 5-45
25
и после преобразований
К = 16,75 + 0,119v + 0,25g
(16)
Уравнение (16) адекватно, поэтому его можно использовать как интерполяционную формулу для вычисления краскоперехода К. Максимальный краско-переход осуществляется при скорости оттиска 160 (отт/час) и гладкости бумаги 70 (с) [7-10].
Наиболее наглядно результаты проведенного полного факторного эксперимента можно изобразить в трехмерном пространстве (рис. 1), где по двум осям отложена значения факторов х1 и х2, а по третьей -значения параметра оптимизации у.
-
№ 3 (108)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
март, 2023 г.
Рисунок 1. Графическое изображение полного факторного эксперимента (ПЭФ) 22
в трехмерном пространстве
Заключение
Условия опытов задаются комбинацией уровней Х1 и Х2 (точки 1, 2, 4 и 3), а результаты (у1, у2, у4 ва у3), отложены параллельно оси у. Точка 0 соответствует центру эксперимента и является началом координат кодированной системы. Контур у1,у2,у4 и у3
отсекает часть поверхности отклика, которая была исследована в пределах интервалов варьирования факторов.
Доказано, что максимальный краскопереход достигается при скорости флексографского печатного оборудования 160 отт./час и гладкости бумаги 70 с.
Список литературы:
1. Спиридонов А.А. Планирование эксперимента при исследовании технологических процессов. -М.: Машиностроение. -1981. -184 с.
2. Ешбаева У.Ж., Джалилов А.А., Шин И.Г. Оценка толщины красочного слоя печатной продукции из опытной образцова бумаги // Universum: технические науки. - 2019. - №. 2 (59). - С. 22-26.
3. Ешбаева У.Ж., Джалилов А.А. Анализ поверхностной обработки композиционной упаковочной бумаги и картона // Universum: технические науки. - 2020. - №. 8-2 (77). - С. 5-9.
4. Ешбаева У.Ж. Печатно-технические свойства новых видов бумаг, содержащих химические волокна. Дисс. на соис. уч. степ. канд. тех. наук // Ташкент. ТИТЛП. - 2008. - С. 230.
5. Ешбаева У .Ж., Джалилов А. А., Рафиков А.С. - Бумага из текстильных отходов. Монография. LAP LAMBERT Academic Publishing, 2018
6. Eshbaeva U.J., Shin I.G., Djalilov A.A. OPTIMIZATION OF COLOR PERCEPTION PROCESS IN THE PRINT PRODUCT BY THE STEEP CLIMBING METHOD BY BOX-WILSON //Technical science and innovation. -2019. - Т. 2018. - №. 4. - С. 37-44.
7. Ешбаева У.Ж., Джалилов А.А. Анализ поверхностной обработки композиционной упаковочной бумаги и картона // Universum: технические науки. - 2020. - №. 8-2 (77). - С. 5-9.
8. Джалилов А.А. Спектроскопическое исследование свойств многослойных целлюлозных композиционных материалов для упаковки //Universum: технические науки. - 2020. - №. 5-2 (74). - С. 5-9.
9. Eshbaeva U. et al. Mathematical Modeling of Ink Transition of Print Product // European Journal of Molecular & Clinical Medicine. - 2021. - Т. 8. - №. 1. - С. 709-717.
10. Eshbaevа U., Nishonov A., Saodatov A. Development of mathematical models of print quality by Box-Wilson // E3S Web of Conferences. - EDP Sciences, 2021. - Т. 304. - С. 03023.
