Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА СЖАТИЯ ОБЪЕМНЫХ УТЕПЛИТЕЛЕЙ ДЛЯ ПРОЕКТИРОВАНИЯ ТЕПЛОЗАЩИТНОЙ ОДЕЖДЫ

Рукавишникова А.С. 1
1 ФБГОУ ВО «Донской государственный технический университет»
При эксплуатации теплозащитной одежды с объемным несвязным утеплителем многие участки подвержены многочисленным деформациям сжатия, которые приводят к снижению исходной величины теплоизоляции. В статье представлена методика исследования процессов сжатия, которая позволяет определять время и скорость деформации при фиксированных нагрузках для утепляющих материалов. В ходе исследований была установлена зависимость между нагрузкой и предельными значениями времени, плотности и объемной деформации. С помощью компьютерных технологий определена зависимость скорости деформации от времени. Разработана квадратичная реологическая 3D модель напряженно-деформированного состояния объемных несвязных утеплителей, описывающая влияние объемной деформации и скорости объемной деформации на давление сжатия. Получена теоретическая зависимость между напряжением и деформацией при одноосном сжатии объемного несвязного утеплителя, учитывающая мгновенную скорость деформации. Разработана математическая модель, которая позволяет определить деформационные характеристики утеплителя и прогнозировать теплозащитные свойства материалов в процессе эксплуатации одежды под действием различных нагрузок. Проведена проверка полученной математической модели на адекватность. Предложенная модель позволяет достоверно описать кинетику сжатия несвязного объемного утеплителя. Она легко адаптируется под различные объемные несвязные утеплители за счет применения коэффициентов уравнения и доступности программного обеспечения.
математическая модель
физико-механические характеристики утеплителей
теплозащитная одежда
деформация сжатия
1. Полимерные материалы в производстве костюма для военнослужащих Арктики (сообщение 3) / Р.Ф. Афанасьева [и др.] // Вестник технологического университета. – 2015. – № 22. – С. 89–93.
2. Исследование релаксационных свойств пуха от температуры / А.Н. Вылкова [и др.] // Инженерный вестник Дона. – 2018. – № 1 [Электронный ресурс]. – URL: http://www.ivdon.ru/ru/magazine/archive/n1y2018/4673 (дата обращения: 22.02.2018).
3. Sovetnikov D.A., Mishakov V.Y., Kirsanova E.A. The study of deformation characteristics of heaters and packages of clothing for military personnel that apply in the arctic zone of the Russian federation // International scientific journal «Theoretical & Applied Science». – 2017. – vol. 48, Issue 04. – P. 145–150.
4. Тунгусова Н.А. Автоматизация проектирования пуховой одежды с использованием математических моделей теплозащитного пакета: автореф. дис. … канд. техн. наук. – Омск, 2009. – 18 с.
5. Иващенко И.Н. Регрессионные модели и оптимизация суммарного теплового сопротивления покровных материалов одежды для нефтяников / И.Н. Иващенко, С.П. Шмалько // Современные проблемы науки и образования. – 2014. – № 3. URL: http://science-education.ru/ru/article/viewid=13517 (дата обращения: 18.02.2018).
6. Лопатченко Т.П. Исследование деформационных свойств объемных утепляющих материалов для проектирования одежды и снаряжения специального назначения / Т.П. Лопатченко, А.С. Рукавишникова. – Ростов н/Д.: Издательский центр ДГТУ, 2014. – 178 с.
7. Рукавишникова А.С., Рукавишников А.А., Лопатченко Т.П. Расчет теплоизоляции комплекта одежды // Свидетельство об официальной регистрации программы для ЭВМ № 2014615366.

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

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

Современные компьютерные технологии позволяют прогнозировать изменение свойств материалов с учетом деформационных характеристик, используя математические модели. Компьютерное моделирование отличается универсальностью, не требует специального оборудования и является достаточно объективным [5]. Таким образом, разработка математической модели является актуальным направлением автоматизации проектирования теплозащитной одежды.

Цель исследования

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

Материалы и методы исследования

Для достижения поставленной цели необходимо провести испытания утепляющих материалов и провести оценку результатов кинетики сжатия. В качестве объектов исследований выбраны объемные несвязные утеплители – гусиный пух категории «Экстра» (96 % пух, 4 % мелкое перо) и hollowfiber из синтетических волокон. Исследуемые теплоизоляционные материалы рассматриваются как вязкоупругое тело, которое в результате многократного сжатия изменяет свои физико-механические свойства [6]. Изменение объемной плотности утеплителя зависит от его реологических свойств, давления и времени воздействия на утеплитель.

Для проведения испытаний была разработана и изготовлена экспериментальная установка, реализующая метод одноосного сжатия. Получение информации о величине и скорости сжатия утеплителя основано на использовании датчика линейного перемещения. Данные передаются на персональный компьютер и обрабатываются программой [7], которая позволяет рассчитывать реологические характеристики исследуемого материала. В методике проведения исследований были учтены потери на трение исследуемого материала о боковые стороны цилиндра установки. При сжатии утеплителя выполняли расчет его текущей высоты в установке hi по формуле (1), высотной деформации εh по формуле (2) и текущее значение плотности ρi, учитывая насыпную высоту h0, изменение высоты образца Δh за время t:

ruk01.wmf, (1)

ruk02.wmf. (2)

При этом ruk03.wmf, где εv – объемная деформация. Диапазон выбранных для исследования давлений от 0,041 до 2,0 кПа, что соответствует давлению человека на материалы при эксплуатации, а также воздействию ветровой нагрузки. Методика проведения исследования представлена на рис. 1.

rukav1.tif

Рис. 1. Схема проведения исследования

Результаты исследования и их обсуждение

Согласно результатам исследований полученные зависимости εv(t, σсж) и ρ(t, σсж) влияния давления σсж и времени сжатия t на величину объемной деформации εv и плотности ρ показали, что существуют предельные значения времени tпред, плотности ρпред и объемной деформации εv пред. При этом при повышении нагрузки уменьшается предельное значение времени и объемной деформации, а также повышается предельное значение плотности.

Предельные нелинейные зависимости εvпред(σсж), ρпред(σсж), tпред(σсж) могут быть описаны сигмоидальным уравнением

ruk04.wmf (3)

где a, b, c, d – коэффициенты уравнений.

Нелинейные зависимости εv(t) при σсж = const могут быть описаны логистическим уравнением

ruk05.wmf (4)

При этом параметры логистического уравнения a, b, c, d зависят от сжимающего напряжения σсж. Используя программу TableCurve, установили связь от 0..tпред влияния времени сжатия на изменение величин, обеспечивающих переход от динамики к статике объемной деформации при фиксированной σсж, используя логистическую зависимость, учитывающую предельные значения:

ruk06.wmf (5)

где t = 0..tпред – время действия нагрузки, с,

εv = 0..εv пред – объемная деформация при t = 0..tпред,

tпред – предельное значение времени сжатия, характеризующее переход от динамического состояния к статическому, с,

εv пред – предельное значение объемной деформации, достигнутое при tпред,

Δεv – разность значений мгновенной и предельной деформации,

tn – время перехода кривой εv от выпуклой к вогнутой его части, с,

d – определяет степень кривизны, εv.

При этом параметры логистического уравнения зависят от сжимающего напряжения σсж и описываются сигмоидальной функцией:

ruk07.wmf (6)

где f, g, k, i – коэффициенты уравнений.

Логистическую зависимость (5) записали в виде, учитывающем влияние σсж на коэффициенты уравнения:

ruk08.wmf (7)

Для определения скорости деформации ?v продифференцировали объемную деформацию (4) по времени t с использованием программы Mathcad 15:

ruk09.wmf. (8)

Используя полученное уравнение (8), построили зависимости ?v (t) при фиксированных σсж в среде «Statistica» (рис. 2).

rukav2.tif

Рис. 2. Зависимость ?v от t при фиксированных значениях σсж, Па. 1 – 41, 2 – 83, 3 – 124, 4 – 165, 5 – 248, 6 – 330, 7 – 413, 8 – 826, 9 – 1652, 10 – 2064

Зависимости ?v(t) (рис. 2) носят экстремальный характер. Минимальные значения ?v(t) наблюдаются при критическом времени tкр = 0,17–0,22. Повышение давления приводит к увеличению абсолютного значения скорости деформации. С увеличением времени t > tкр значения функции |?v(t)| уменьшаются и стремятся к 0. Построили с помощью программы «Statistica» квадратичную 3D модель σсж(?v, εv) при условиях t < tпред и εv < εv пред для натурального и синтетического утеплителя (рис. 3) и определили ее параметры:

ruk10.wmf (9)

где а0, а01, а02, а11, а12, а22 – параметры уравнения сжатия.

Влияние εv на σсж носит экстремальный характер. Максимальные значения σсж наблюдаются для перо-пуховой смеси «Экстра» при εv пред = –1,4 (рис. 3, а). Максимальные значения σсж для hollowfiber наблюдаются при εv пред = –0,7 для σсж = –200 и при εv пред = –0,2 для σсж = –1300 (рис. 3, б). Параметры уравнения уплотнения при εv > εv пред, εv < εv пред и для всего исследуемого диапазона значений представлены в таблице.

rukav3a.tif rukav3b.tif

а б

Рис. 3. 3D модель σсж(?v, εv): а – перо-пуховая смесь, б – hollowfiber

Параметры уравнения сжатия

Вид утеплителя

Исследуемый диапазон

а0

а01

а02

а11

а12

а22

Перо-пуховая «Экстра»

εv = 0..–3,6

–739,99

–1085,35

–3,51

–385,59

–2,62

–7,37

Hollowfiber

εv = 0..–1,6

–471,18

–1231,70

–6,31

–905,55

–73,54

–8,69

Из представленных графиков видно, что характер зависимости синтетического утеплителя значительно отличается от характера зависимости натурального, которые имеют симметричный вид относительно предельной объемной деформации. Характер зависимости σсж(?v, εv) для hollowfiber значительно меняется с увеличением объемной деформации выше критической. Данное явление характерно для исследуемого диапазона значений давления сжатия.

Зависимость влияния времени нагружения и деформации материала утеплителя на напряжение σ могут быть описаны уравнениями изохронных кривых. Для описания зависимости σсж(εv, t) для несвязного объемного утеплителя уравнение изохронных кривых представили в виде

ruk11.wmf (10)

где φ(εv) – кривая мгновенного деформирования,

ψ(t) – зависимость, учитывающая время.

Подставляя данные, рассчитанные по формуле (4) при t = 0 в TableCurve, построили кривую мгновенного деформирования (уплотнения) φ(εv), используя импульсную зависимость, и определили коэффициенты уравнения:

ruk12.wmf (11)

где а1, b1, c1, d1, e1 – коэффициенты уравнения,

зависимость, учитывающую время ruk13.wmf, с учетом (10), записали в виде

ruk14.wmf (12)

Отсюда в TableCurve построили кривую мгновенного деформирования и определили зависимость εv(t) и коэффициенты:

ruk15.wmf, (13)

где a2, b2, c2, d2, e2 – коэффициенты уравнения.

С учетом (11) и (13) уравнение (10) принимает вид

ruk16.wmf (14)

Данное уравнение позволяет адекватно описать результаты экспериментальных исследований кинетики сжатия несвязного объемного утеплителя.

Заключительный этап анализа полученной математической модели предусматривает проверку ее адекватности. Уравнение (14) можно считать адекватным с 95 % доверительной вероятностью, учитывая, что Fp < Fтабл. Данная модель позволяет достоверно описать результаты экспериментальных исследований кинетики сжатия несвязного объемного утеплителя. Модель легко адаптируется под различные объемные несвязные утеплители за счет применения коэффициентов уравнения и доступности программного обеспечения.

Выводы

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

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

Построена квадратичная реологическая 3D модель, описывающая влияние объемной деформации и скорости объемной деформации на давление сжатия.

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


Библиографическая ссылка

Рукавишникова А.С. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА СЖАТИЯ ОБЪЕМНЫХ УТЕПЛИТЕЛЕЙ ДЛЯ ПРОЕКТИРОВАНИЯ ТЕПЛОЗАЩИТНОЙ ОДЕЖДЫ // Современные наукоемкие технологии. – 2018. – № 4. – С. 113-117;
URL: https://top-technologies.ru/ru/article/view?id=36969 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674