CC BY
52
УДК 678.01+519.7
DOI: http://dx.doi.org/10.20914/2310-1202-2016-1-173-181 Профессор О.В. Карманова
(Воронеж. гос. ун-т. инж. технол.) кафедра химии и химической технологии органических соединений и переработки полимеров
профессор С.Г. Тихомиров, аспирант А.М. Скачков, студент А.Г. Иванов
(Воронеж. гос. ун-т. инж. технол.) кафедра информационных и управляющих систем. тел. (473) 249-92-37 E-mail: [email protected]
Professor O.V. Karmanova
(Voronezh state university of engineering technology) Department of chemistry and chemical technology organic compounds and polymer processing
professor S.G. Tikhomirov, graduate A.M. Skachkov, student A.G. Ivanov
(Voronezh state university of engineering technology) Department of information and control systems. phone (473) 249-92-37 E-mail: [email protected]
Моделирование упруго-прочностных свойств эластомеров
Modeling of elastic-strength properties of elastomers
Реферат. Разработана модель «состав-свойство», учитывающая структурную неоднородность полимерных композиций. В качестве объектов исследования использовали модельные композиции на основе бутадиен-стирольного каучука СКС-30АРК и высоковязкого заструктурированного полимера (высокомолекулярный наполнитель - ВМН) в различных соотношениях, а также мягчителей (индустриальное масло И-12А, масло ПН-6, низкомолекулярный полибутадиен ПБН), что позволяло получать образцы, значительно отличающиеся по вязкости. На основе модельных полимерных композиций изготовлены резиновые смеси и получены их вулканизаты. Определены физико-механические свойства - условная прочность при растяжении, относительное удлинение при разрыве, твердость по Шору А.. При разработке модели принято допущение, что для описания физико-механических свойств полимеров справедливо правило логарифмической аддитивности. Свойства полимерной композиции (ПК) определялись одним доминирующим компонентом (полимерная композиция, состоящая из каучука и высокомолекулярного наполнителя) и дополнительными компонентами (мягчите-ли). Для оценки параметров модели использовали алгоритм идентификации, который включал четыре этапа. Реализация данного алгоритма проводилась с использованием методов планирования эксперимента. Оценку неизвестных параметров в уравнении осуществляли с использованием метода наименьших квадратов. Оценка качества модели, проведенная по критериям Фишера, поворотных точек, Дарбина-Уотсона, R/S-критерию показала, что модель адекватно описывает изменение физико-механических свойств в зависимости от состава полимерных композиций. По расчетным данным построены графические 3d-зависимости физико-механических свойств от состава полимерных композиций, позволяющие оценить вклад доминирующего компонента и дополнительных компонентов (в том числе в комбинации) в изменение показателей. Установлено, что при введении в каучук суммарно более 50 % компонентов (ВМН и мягчителей) снижается условная прочность при растяжении и резко возрастает относительная погрешность расчетов по модели.
Summary. Model "structure-property", which takes into account the structural heterogeneity of polymer compositions has been developed. Experimental compositions based on styrene-butadiene rubber SCS 30ARK and crosslinked, high viscosity polymer (high-molecular filler - VMN) in different proportions, as well as softeners (industrial oil I-12A, oil PN-6, low-molecular polybutadiene PBN) were investigated. Samples that differ significantly in viscosity were obtained. The rubber blends and vulcanizates, based on the experimental of polymer compositions, were prepared. Physico-mechanical properties - tensile strength, elongation at break, Shore A. A hardness were determined. For describe the physical and mechanical properties of polymers logarithmic additivity rule was used. The properties of the polymer composition (PQ were determined by a single dominant component (resin composition consisting of a high rubber and a filler) and further components (softeners). Identification algorithm consists of four steps. The implementation of this algorithm is carried out using experimental design techniques. Estimation of the unknown parameters in the equation was carried out using the method of least squares. Quality evaluation of the model was conducted with the criteria Fisher, turning points, the Durbin-Watson, R / S-criterion. It is found that the model adequately describes the change of physico-mechanical properties depending on the composition of polymer compositions. 3d graphics of the physical-mechanical properties of the polymer compositions were built. This allowed us to estimate the contribution of the dominant component and optional components (including combinations thereof) to change the parameters. It has been established that the introduction of rubber in total more than 50% of the components (BMH and softeners) reduced conventional tensile strength and dramatically increases the relative error of model calculations. (BMH and softeners) reduced conventional tensile strength and dramatically increases the relative error of model calculations.
Ключевые слова: бутадиен-стирольные каучуки, математическая модель, физико-механические показатели.
Keywords: styrene-butadiene rubber, mathematical model, physical and mechanical properties.
Для цитирования Карманова О.В., Тихомиров С.Г., Скачков А.М., Иванов А.Г. Моделирование упруго-прочностных свойств эластомеров // Вестник Воронежского государственного университета инженерных технологий. 2016. №1. С. 173-181. doi:10.20914/2310-1202-2016-1-173-181.
© Карманова О.В., Тихомиров С.Г., Скачков А.М., Иванов А.Г., 2016
For cite
Karmanova O.V., Tikhomirov S.G., Skachkov A.M., Ivanov A.G. Modeling of elastic-strength properties of elastomers. Vestnik voronezhskogo gosudarstvennogo universiteta inzhenernyh tekhnologij [Proceedings of the Voronezh state university of engineering technologies]. 2016, no. 1, pp. 173-181. (In Russ.). doi: 10.20914/ 2310-1202-2016-1-173-181.
Введение. Полимерные материалы в современных условиях являются одними из перспективных конструкционных материалов, применяемых в различных отраслях науки и техники. На протяжении многих лет установление связи структуры полимера с его механическими свойствами остается актуальной задачей, так как это позволяет не только предсказывать направление изменения комплекса механических свойств при изменении структуры, но и получать новую информацию о структуре полимера, если определены его физико-механические свойства (ФМС). Вследствие протекания физических и химических процессов в ходе переработки полимеров в изделия изменяется их структура, а присутствие в полимерной композиции помимо полимера наполнителей и других целевых добавок, обусловливает образование коллоидно-химической структуры [1]. Кроме того, в условиях современного производства при повышении скоростей переработки полимерных композиций увеличивается воздействие внешних механических сил на макромолекулы и структурные образования, что сопровождается механическими превращениями и может привести к ухудшению ФМС. В этой связи механическая прочность полимерных композиций в процессе ее переработки в изделие становится определяющим фактором установления оптимальных параметров технологического процесса.
В настоящее время особую актуальность приобретают исследования, связанные с разработкой и использованием математических моделей для прогнозирования поведения изделий в эксплуатации. В технологии эластомеров основные свойства будущего изделия формируются на стадии изготовления композиций: комплекс технологических свойств и вулканизационных
характеристик эластомерных композиций связан с физико-механическими показателями вулкани-затов и техническими свойствами готовых изделий. Варьируя состав полимерной композиции, можно прогнозировать и регулировать свойства эластомеров на всех стадиях технологического процесса их получения и эксплуатации.
Целью исследования явилась разработка модели «состав-свойство», учитывающей структурную неоднородность ПК.
В качестве объектов исследования использовали модельные композиции на основе бутадиен-стирольного каучука СКС-30АРК и высоковязкого заструктурированного полимера (высокомолекулярный наполнитель - ВМН) в различных соотношениях, а также мягчителей (индустриальное масло И-12А, масло ПН-6, низкомолекулярный полибутадиен ПБН), что позволяло получать образцы, значительно отличающиеся по вязкости. На основе модельных полимерных композиций с использованием микросмесителя при температуре 60±5 оС были изготовлены резиновые смеси и получены их вулканизаты при температуре 155 оС и продолжительности вулканизации 20 минут. Физико-механические свойства - условную прочность при растяжении, относительное удлинение при разрыве - определяли на разрывной машине РМИ-60 в соответствии с ГОСТ 270-75, твердость по Шору А - с ГОСТ 263-75.
Теоретическая часть. Принято допущение, что для вязких систем, таких как полимеры с молекулярной массой больше некоторой критической величины для описания физико-механических свойств (ФМС) справедливо правило логарифмической аддитивности [2], которое можно представить в виде:
п / \ п У = СП )=Срц( *,)П
1=1
1=2
д9„ (0) 1 5 (0)
ЭХ:
+--
2
дХ
х 1 + 0( х3)),
(1)
у П (1 -
1=2
V» (0)
д%1 (°) дх.
1
+ —
V» (0)
54,(0)
дх2
) «70,П(1 + а»х» + Ь»х2),
»=2
х
2
1
1
2
п
где У0, = оя(х1) Пф»(0) - свойство полимер-
»=2
ной композиции при отсутствии в ее составе ВМН и пластифицирующих добавок; а^, Ь - коэффициенты зависимости (1); С - константа, характеризующая структуру полимерных цепей, , - свойство (условная прочность при растяжении,
относительное удлинение, твердость по Шору А), компонент (добавки-мягчители).
Следует отметить, что правило логарифмической аддитивности применимо в определенных диапазонах изменения температуры, напряжения и молекулярной массы [3]. При больших напряжениях, высоких темпера-
турах, значительной доле наполнителя оно нарушается из-за глубокого изменения надмолекулярной структуры бутадиен-стирольного сополимера [1]. В том случае, когда данное правило выполняется, компоненты действуют на ФМС независимо друг от друга.
Свойства полимерной композиции (ПК) определяются одним доминирующим компонентом, в нашем случае это полимерная композиция, состоящая из каучука и высокомолекулярного наполнителя. Под доминирующим компонентом понимается компонент, содержание которого в композиции составляет не менее 50 % масс., остальные компоненты являются дополнительными [4].
Оценку параметров зависимости (1) осуществляли с использованием алгоритма идентификации с выбором доминирующего компонента согласно [2]. Процесс оценки параметров модели состоял из следующих этапов:
1. Определение исследуемого ФМС вул-канизатов исходного каучука (без добавления ВМН и добавок-мягчителей) как среднего значения по нескольким параллельным измерениям:
1 к -
4 = ТЕ1 = 1,М, (2) К j=l
где i, K - номер и число измеренных значений и свойств полимерной композиции при наличии параллельных измерений; i, M - номер и число значений свойство полимерной композиции в общей выборке.
2. Оценка влияния доминирующего компонента в соответствии с уравнением:
Y0i (X1, X2 ) = Y0i (Xl) ' (1 + а1,2 X2 + Ь1,2 X2) , (3)
где xi - массовая доля исходного каучука, Х2 - массовая доля ВМН.
3. Исследование ПК после введения дополнительных компонентов согласно (1) и оценка параметров зависимости.
4. Определение свойств вулканизатов полимерных композиций с учетом влияния всех компонентов с целью оценки адекватности полученной модели.
Реализация данного алгоритма проводилась с использованием методов планирования эксперимента. Оценку неизвестных параметров в уравнении осуществляли с использованием метода наименьших квадратов, исходя из условия минимума среднеквадратичной ошибки [5]:
м / \2 R = Y(y; - Yp (х, х2)) -> min
Относительную погрешность вычисляли по формуле:
Уэ - Yp
5 =
Ys
х100%
Оценку адекватности модели проводили по критериям Фишера, поворотных точек, Дарбина-Уотсона, Я^-критерию.
Экспериментальная часть и обсуждение результатов. Моделирование влияния ВМН и пластифицирующих добавок на условную прочность при растяжении вулка-низатов проведено в соответствии с выбранным алгоритмом.
Влияние ВМН на доминирующий компонент в соответствии с уравнением (1) описывается зависимостью:
fp (x^ X2 ) = fp0 (l
+ a, 2 X2 + 2 X2
), (4)
где fpo = 16,6 МПа - экспериментальное значение прочности при 100 % каучука СКС-30 АРК.
Результаты идентификации параметров зависимости (4) представлены в таблице 1.
Т а б л и ц а 1 Влияние содержания доминирующего компонента ВМН в полимерной композиции на условную прочность при растяжении вулканизатов
Состав полимерной композиции, % Средние значения условной прочности Относительная
при растяжении по погрешность,
результатам 5 измерений, МПа %
Xi X2 fp (X1> X2') эксп fp(X1,X2) расч
100 0 16,6 16,6 0
90 10 13,4 14,8 9,22
80 20 13,2 13,4 1,29
70 30 13,2 12,4 6,19
60 40 12,4 12,0 3,88
50 50 11,4 11,9 4,15
Среднее значение относительной 4,94
погрешности
Проверка качества модели подтвердила ее адекватность: коэффициент детерминированности Л2=0,9624>0,9000, критерий Фишера Fрасч = 38,41 > Fтабл = 9,55, что свидетельствует о высокой степени приближения кривой к экспериментальным данным. Следовательно, гипотеза о случайной природе оцениваемых характеристик отклоняется и признается их статистическая значимость и надежность.
Моделирование влияния дополнительных компонентов индустриального масла И-12А (хз), масла ПН-6 (х*) и ПБН (хз) на условную прочность при растяжении вулкани-затов описывается зависимостями:
fp (х!, х2 , хз ) = fp (х1, х2 )(1 + «1,3хз + Ь1,3х32 ), (5) !р (х1, х2 , х4 ) = Л (х1, х2 X1 + «1,4х4 + Ь1,4х4 ), (6) ^ (х1, х2 , х5 ) = ^ (xl, х2 )(1 + «1,5 х5 + Ь1,5 х52 ). (7)
Результаты параметрической идентификации зависимостей (5-7) приведены в таблице 2 и обобщены на рисунке 1. За базовую линию принято значение условной прочности при растяжении 12 МПа, удовлетворяющее требованиям технических норм для резинотехнических изделий общего назначения. Выявлено, что вид и дозировка мягчителя оказывают влияние на изменение условной прочности при растяжении образцов, которая уменьшалась в ряду «ПБН>ПН-6>И-12А» при одном и том же содержании ВМН. Стабилизирующее действие ароматического масла ПН-6 на проч-
ность по сравнению с минеральным И-12А объясняется его лучшей совместимостью с бу-тадиен-стирольным сополимером. Введение мягчителей более 20 % отрицательно сказывается на прочностных показателях, однако добавлением ВМН удается компенсировать снижение прочности (рисунок 1). Композициони-рованием ВМН и мягчителей удается обеспечить базовый уровень прочности для большинства исследуемых образцов.
В ходе статистического анализа модели (5) установлено, что значение коэффициента детерминированности R2 = 0,8279 < 0,9000 свидетельствует о достаточно высокой степени приближения расчетной кривой к экспериментальным данным. Оценка качества модели по критерию Фишера показала, что Fрасч = 33,66 > Fmабл = 3,81. Следовательно, гипотеза о случайной природе оцениваемых характеристик отклоняется и признается их статистическая значимость и надежность.
Аналогичные выводы были сделаны для моделей (6) и (7): R2 = 0,7503 и 0,5887, Fрасч = 21,03 и 10,02, соответственно.
Таблица 2
Результаты определения условной прочности при растяжении по зависимостям (5) - (7)
оотношения компонентов, %
хг х2 хз Х4 Х5
Средние значения условной прочности при растяжении вулканизатов по результатам 5 измерений, МПа
Относительная погрешность, %
fp
(хг,х2,хз) эксп
fp
(хг ,х2,хз)расч
15,0
15,5
90 5 5 5 5
80 10 10 10 10
70 15 15 15 15
50 25 25 25 25
90 6 4 4 4
80 12 8 8 8
70 18 12 12 12
50 30 20 20 20
90 7 3 3 3
80 14 6 6 6
70 21 9 9 9
50 35 15 15 15
90 9 1 1 1
80 18 2 2 2
70 27 3 3 3 Среднее значение относительной погрешности
13,7 12,2 8,4
15.1
14.6
13.0 8,6 16,0
14.3
13.2 9,7
16.7
15.4
13.1
14.0 12,2 7,7 15,4
14.1 12,6 9,4 15,3 14,0 12,8 10,6 15,0
13.6
12.7
fp
(хг,х2,х4)эксп
fp
(хг,х2,х4)расч
16,7
16,1
15.2
13.3 9,8 17,0
15.0 14,3
10.1 16,0
15.3 13,5
10.4 14,3
14.5 13,3
15.0
13.5 9,1 15,9 14,9
13.6
10.7 15,6 14,6
13.6
11.7
15.1 13,9 13,0
fp
(хг,х2,х5) эксп
fp
(хг ,х2,Х5)расч
17,8
16,6
16,6
14.0
9.6
18.5 16,3 13,7
10.1 17,0 16,2 14,2
9.7 15,0
14.6 13,5
15,7
13,9
8,5
16,3
15,5
14,2
10,7
16,0
15,11
14.1 12,1
15.2
14.1
13.2
fp fp (хг,х2,хз) (хг,х2,х4)
3,43 2,43 0,0 9,09 2,05 3,88 3,39 8,22 4,71 2,22 3,42 8,18 11,69 12,88 3,33 5,30
3,57 1,06 1,21 7,49 6,96 0,84 4,92
5.29
2.30 4,67 0,67 10,76 5,14 4,61 2,49 4,30
fp
(хг,х2,х5) 7,17 6,00 0,57 12,63 13,48 5,48 3,40 5,65 6,50 7,20 0,43 19,50 1,26
3.82 1,94
6.83
Вестник.ВГУИТ, №1, 2016_
а)
Ху, %
б)
Xz, %
в)
X;, % 1
Рисунок 1. Зависимости изменения расчетных значений прочности при растяжении вулканизатов (/) от содержания доминирующего (х2) и дополнительных (х3, х4, х5) компонентов: а - И-12А, б - ПН-6, 3 - ПБН
На заключительном этапе проведено моделирование совместного действия пластифицирующих добавок И-12А, ПН-6 и ПБН на условную прочность при растяжении вулкани-затов полимерных композиции (таблица 3). Оценка параметров искомой зависимости осуществлена в соответствии с уравнением:
fp (x1 , X2 , X3, X 4 , X 5 ) _ = fp (x1 , X2 )fp (x1 , X2 , X3 )x X fp (x1 , X2 , X4 )fp (x1 , X2 , X5 )
(8)
Результаты моделирования обобщены на рисунке 2 и в таблице 3.
- эксп расч
Рисунок 2. Зависимости ранжированных значений условной прочности при растяжении (f) полимерных композиций в порядке ее возрастания от содержания тройной комбинации дополнительный компонентов И-12А, ПН-6 и ПБН
Проведена дополнительная оценка адекватности модели (8) с использованием критерия поворотных точек, критерия Дарби-на-Уотсона и R/S- критерия.
Случайность остаточного компонента по критерию поворотных точек подтверждается, если выполняется неравенство:
М >
2(n - 2) 2 /16и - 29
3
90
(9)
где М - число поворотных точек, п - число измерений в выборке.
Используя формулу (9) получены оценки: «9 > 6,221 - верно», следовательно, случайность остаточного компонента по критерию поворотных точек подтверждается, модель адекватна.
При определении автокорреляции остатков критерием Дарбина-Уотсона:
-£t-i)2
d = ^-
n
(10)
получаем d = 1,24; di = 1,24 < d2 = 2,76. Следовательно, автокорреляция остатков отсутствует, модель адекватна.
Исследование независимости распределения остаточного компонента по R/S- критерию проведено по уравнению
Ц/ _ _Smax Smm
'S
n f n Л2
nZ S-i]T
i=1 V i=1 n(n -1)
(11)
8
6
14
12
с 10
4
2
4
6
8
10
12
14
16
18
n
t=1
Т а б л и ц а 3
Влияние содержания дополнительных компонентов ПБН, ПН-6 и И-12А при их совместном
введении в полимерные композиции на условную прочность при растяжении вулканизатов
№ опыта Состав полимерной композиции, % Средние значения условной прочности при растяжении по результатам 5 измерений, МПа Относительная погрешность, %
х1 х2 хз Х4 х5 /р эксп /р расч
1 90 5 1,7 1,7 1,7 16,2 16,4 1,01
2 80 10 3,3 3,3 3,3 15,0 15,9 5,62
3 70 15 5 5 5 13,8 15,3 9,68
4 50 25 8,3 8,3 8,3 10,0 13,7 26,99
5 90 6 1,3 1,3 1,3 16,0 16,0 0,16
6 80 12 2,7 2,7 2,7 14,6 15,4 5,33
7 70 18 4 4 4 12,6 14,8 14,66
8 50 30 6,7 6,7 6,7 9,8 13,5 27,43
9 90 7 1 1 1 15,5 15,7 1,41
10 80 14 2 2 2 13,3 14,9 10,84
11 70 21 3 3 3 12,3 14,2 13,48
12 50 35 5 5 5 9,7 13,2 26,66
13 90 9 0,3 0,3 0,3 16,6 15,1 10,17
14 80 18 0,7 0,7 0,7 14,3 13,9 2,83
15 70 27 1 1 1 14,0 13,0 7,33
16 50 45 1,7 1,7 1,7 11,0 12,4 16,27
Среднее значение относительной погрешности 11,24
Рассчитанное значение критерия R/S=3,32 лежит в пределах заданного табличного интервала (2,94; 4,49). Следовательно, независимость распределения остаточного компонента подтверждается, а модель адекватна.
В ходе анализа экспериментальных данных (таблица 3) установлено, что при введении в каучук СКС-30АРК суммарно более 50 % компонентов (ВМН (х2) и мягчителей (х3, Х4, Х5)), то есть при соотношении компонентов Х1:(Х2+Х3+Х4+Х5)=50:50, снижается условная прочность при растяжении и резко возрастает относительная погрешность расчетов по модели (номера опытов 4, 8, 12, 16). Это объясняется снижением межмолекулярного взаимодействия в полимерной матрице при высоких дозировках масел и низкомолекулярного полибутадиена, общее содержание которых составило для опытов 4, 8, 12 от 15 до 25 %. С другой стороны высокое содержание ВМН при малых дозировках мягчителей (45%, опыт №16) обусловливает увеличение неоднородности полимерной основы за счет агломерации частиц ВМН и, как следствие, способствует неравномерному нагружению макромолекул.
Влияние доминирующего компонента (ВМН) на относительное удлинение при разрыве и твердость вулканизатов полимерных композиций в соответствии с уравнением (1) описывается зависимостями, соответственно:
х2 ) = е0 (1 + а2,2х2 + Ь2,2Х22 ) (12)
НА (xl, Х2 ) = НА0 (1 + а3,2Х2 + Ь3,2Х22 ) (13)
где £о = 360 % - экспериментальное значение относительного удлинения при разрыве вулканиза-та при 100 % содержания каучука СКС-30 АРК, Нао = 75 усл. ед. для того же этого образца. Результаты идентификации параметров зависимости (12) представлены в таблице 4 и рисунке 3. Статистический анализ моделей (12) - (13) показал, что значение коэффициента детерминированности Е2 = 0,9977 > 0,9000 для относительного удлинения и Е2 = 0,9985 > 0,9000 свидетельствует о достаточно высокой степени приближения расчетных значений к экспериментальным данным. Оценка качества моделей по критерию Фишера показала, что для модели относительного удлинения FрCKH = 647,49 > Fmабл = 9,55, для модели твердости Fрасч = 978,66 > Fmабл = 9,55, то есть признается их статистическая значимость и надежность.
е(х1 , Х2, Х3) — Е(Х1 , Х2 )(1 ^ а2 3 Х3 ^ ¿2 3 Х3 ),(14) е(х1 , х2, х4) — е(х1 , х2 )(1 + а2,4 х4 + Ь2,4 х42) ,(15) г(х1, х2, х5) — а(х1, х2 )(1 + а2 5х5 + Ь2 5х25 ).(16)
НА (х1, х2 , х3 ) — НА (х1, х2 )(1 + а3,3х3 + Ь3,3х3 ) (17) на (х1 , х2, х4 ) — на (х1 , х2 )(1 + а3,4х4 + Ь3,4х4 ) (18) н а (х1, х2, х5 ) — н а (х1, х2 )(1 + а3,5 х5 + Ь3,5 х5 ) (19)
Т а б л и ц а 4
Влияние доминирующего компонента на относительное удлинение при разрыве и
твердость вулканизатов
Состав полимерной композиции, % Относительное удлинение при разрыве Твердость по Шору А
Средние значения относительного удлинения по результатам 5 измерений, % Относительная погрешность, % Средние значения твердости по результатам 5 измерений, усл. ед. Относительная погрешность, %
хг х2 В(х1,х2)зксп б(х1,х2) расч На(Х:,Х2)зксп На(х:,х2) расч
100 0 360 360 0 75 75 0
90 10 305 318 4,10 81,3 81 0,36
80 20 280 277 1,14 86,6 87 0,48
70 30 240 236 1,49 90,9 91 0,96
60 40 200 197 1,55 94,2 94 0,90
50 50 155 158 2,02 96,4 96 0,44
Среднее значение относительной погрешности 2,06 0,63
а)
90
и
> 80
60
20 Х2, % 30
б)
0 10 20 30 40 50 60 Х2, %
Рисунок 3. Зависимости относительного удлинения при разрыве (а) вулканизатов и твердости по Шору А (б) от содержания ВМН (х2) 1 - эксперимент; 2 - расчет
Результаты параметрической идентификации зависимостей (14) - (19) приведены
в таблице 5 и обобщены на рисунках 4-5. За базовые значения приняты: по относительному удлинению - 300 %, по твердости - 70 усл. ед. Установлено, что использование в композиции ВМН более 15 % снижает относительное удлинение и повышает твердость, особенно для образцов с индустриальным маслом.
Таблица 5
Результаты статистического анализа влияния дополнительных компонентов Х3 , х*, Х5 на относительное удлинение при разрыве
Состав полимерной композиции, % Среднее значение относительной погрешности, % Коэффициент детерминированности ^^ (0,9000) Критерий Фишера Fрасч> Fтаб
Относительное удлинение при разрыве
хг х2 Хз 3,96 0,7371 14,72>3,81
хг х2 х4 5,02 0,5011 7,03>3,81
хг х2 Х5 6,62 0,8366 35,83>3,81
Твердость
хг х2 Хз 3,18 0,8747 48,89>3,81
хг х2 х4 2,12 0,8658 45,17>3,81
хг х2 Х5 2,97 0,8309 34,40>3,81
Влияние тройной комбинации мягчи-телей на изменение относительного удлинения при разрыве и твердость по Шору А. описывается зависимостями, соответственно
НА (х1>
х т, х2, х3, х*, х5
х3 1, х2, х4 М^ х2 , х5 ), ) = НА (х1, х2 )НА (х1, х2, х3 )НА (х1, х2, х4 )НА (х1, х2, х5 )
£>(х!, х2, х3 , х*, хх5 ) ^(х!, х2 )^(х1, х2 •
(20) (21)
400
300
200
100
0
0
10
40
50
00
70
Таблица 5 Результаты статистического анализа влияния дополнительных компонентов Х3 , Х4, Х5 на свойства
Состав полимерной композиции, % Среднее значение относительной погрешности, % Коэффициент детерминированности R2 (0,9000) Критерий Фишера Fpac4> Fma6
Относительное удлинение при разрыве
X1 X2 хз 3,96 0,7371 14,72>3,81
X1 X2 Х4 5,02 0,5011 7,03>3,81
X1 X2 Х5 6,62 0,8366 35,83>3,81
Твердость
X1 X2 хз 3,18 0,8747 48,89>3,81
X1 X2 Х4 2,12 0,8658 45,17>3,81
X1 X2 Х5 2,97 0,8309 34,40>3,81
е,%
%
а)
Х2, % б)
400-500 300-400 200-300 100-200 0-100
в)
Рисунок 4. Зависимости изменения расчетных значений относительного удлинения при разрыве вулканизатов (г) от содержания доминирующего (х2) и дополнительных (хз, х4, х5) компонентов: а - И-12А, б - ПН-6, в - ПБН
а)
в)
Рисунок 5. Зависимости изменения расчетных значений твердости по Шору А вулканизатов (На) от содержания доминирующего (x2) и дополнительных (X3, X4, X5) компонентов: а - И-12А, б - ПН-6, в - ПБН
Дополнительная оценка адекватности моделей (20) - (21) с использованием критерия поворотных точек подтвердила их адекватность: для относительного удлинения получено 11 > 6,221; для твердости - 10 > 6,221. Определение автокорреляции остатков критерием Дарбина-Уотсона показало, что автокорреляция остатков отсутствует, модели адекватны: для относительного удлинения получено d1 = 1,68 < d2 = 2,32 ; для твердости -d1 = 1,68 < d2 = 2,32 . Оценка независимости распределения остаточного компонента по R/S- критерию показала, что рассчитанные значения критерия R/S =4,27 для относительного удлинения и R/S =3,14 попадают внутрь табличного интервала (2,94; 4,49), что свидетельствует о независимости распределения остаточного компонента и подтверждает адекватность модели.
Найденные параметры модели «состав-свойство» в матричном виде имеют вид:
' 0 -1,24 10-2 2,54 10-3 1,09 10-2 1,91 10-2
А = (а ) = V J /3x5 0 -1,18 10 2 5,63 10-3 5,85 10-3 2,58 10-
0 V 9,06 10-3 -1,3410-2 - 9,24 10-3 -1,03 10
B=(b,jL = '0 1,35 10 4 - 7,40 10 4 - 9,00 10 4 -1,30 10-3
0 1,12 10-5 - 2,95 10 4 1,04 10-5 - 3,10 10 4
0 - 6,70 10-5 3,56 10 4 2,18 10 4 2,36 10 4
Таким образом, моделирование упруго-прочностных свойств вулканизатов, содержащих ВМН и мягчители разной природы, позволяет получать эластомеры с прогнозируемыми свойствами и корректировать составы композиций в зависимости от предъявляемых требований.
ЛИТЕРАТУРА
REFERENCES
1. Вострокнутов Е. Г. , Новиков М.И., Новиков В.И. Переработка каучуков и резиновых смесей (реологические основы, технология, оборудование). М.: Химия, 2005. 391 с.
2. Тихомиров С.Г., Карманова О.В., Скачков А.М.,. Дьяков А.А Моделирование технологических свойств полимерной композиции с выделением доминирующего компонента. Промышленное производство и использование эластомеров. 2015. №3. С.
3. Бухонов Б.П., Ветохин В.Н., Грановская Г.Л., Тихомиров С.Г. и др. Молекулярное строение и вязкость. дивинильного каучука // Высокомолекулярные соединения.1992. №5. С. 11-14.
4. Николаева С.В. Системный анализ многокомпонентных пищевых объектов и технологий в условиях информационной неопределенности: дисс. ...док. тех. наук; 05.13.01. М., 2013. 394 с.
5. Математическая статистика. Под ред. В.С. Зарубина, А.П. Крищенко / В.Б. Горяинов, И.В. Павлов, Г.М. Цветкова, и др. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. - 424 с.
1. Vostroknutov E.G., Novikov M.I., Novikov V.I. Processing of rubber and rubber compounds (rheolog-ical bases, technology, equipment).). M .: Chemistry, 2005. 391 p. . (In Russ.).
2 Tikhomirov S.G., Karmanova O.V., Skachkov A.M., Dyakov A.A. Modeling of technological properties of the polymer composition with the release of the dominant component // Industrial production and use of elastomers, 2015, №3, p. 16-18. . (In Russ.).
3. Buhonov BP, Vetohin V.N., Granovskaya G.L., Tikhomirov S.G., Shabanov S.V., Lurie V.L., Garshin A.P. Molecular structure and viscosity divinyl rubber // High-molecular compounds. M. T. 34. - №5 - 1992. - p. 11-14. . (In Russ.).
4. Nikolaev S.V. System analysis of multi-component food items and technologies in information neoprene-division: dis. Doctor of Technical Sciences: 05.13.01. M., 2013. 394. . (In Russ.).
5 Zarubin V.S., Kryschenko A.P., Goryai-nov V.B.,. Pavlov I.V. Mathematical statistics. M .: Publishing House of the MSTU named N.E. Bauman, 2008. 424 p. (In Russ.).
