Постановка математической задачи и методы решения разработки антистатической одежды для защиты от повышенных температур
Автор: Пулатова С. У., Кодирова С. Х., Ташпулатов С. Ш., Сафаров И. А., Акбаров Р. Д., Черунова И. В.
Журнал: Вестник Алматинского технологического университета @vestnik-atu
Рубрика: Техника и технологии
Статья в выпуске: 2 (132), 2021 года.
Бесплатный доступ
Статья посвящена исследованию и разработке математической модели процесса электризации текстильных материалов в системе «человек - одежда - среда» в жарких климатических условиях. Разработанная математическая модель описывает процессы электризации в поперечных сечениях пакета одежды и на поверхностном слое кожи человека и позволяет рассчитать параметры антистатической одежды для защиты от повышенных температур. Установлено, что в области распределения температуры по длине одежды в четырех различных вариантах было обнаружено, что в I - варианте , во II - варианте , в III - варианте , в IV - варианте , в V - варианте , в VI - варианте длина одежды в этой области увеличилась на 10% по сравнению с длиной одежды во II - варианте.
Математическая модель, электростатическая безопасность, цилиндр, система, деформация одежды, теплообмен
Короткий адрес: https://sciup.org/140256805
IDR: 140256805 | DOI: 10.48184/2304-568X-2021-2-12-23
Текст научной статьи Постановка математической задачи и методы решения разработки антистатической одежды для защиты от повышенных температур
Вопросы обеспечения электростатической безопасности определяются вариантом состава и параметрами материалов. Совокупность этих взаимозависимых характеристик является входной информацией для выполнения инженерных расчетов в области конструирования одежды.
Основы математического моделирования определяются первоначально представлениями процесса электризации текстильных материалов в рамках системы «Человек - Одежда - Среда». Явление статической электризации, возникающее в процессах изготовления текстильных материалов и изделий из них, а также при эксплуатации готовых изделий, ведет к снижению качества продукции, доставляет неудобства в работе, нарушает технологические процессы производства, а при возникновении электрического разряда представляет угрозу человеку и может приводить к несчастным случаям.
Модели, описывающие процесс возникновения электрического заряда, его влияние на производственный процесс и на организм человека, рассматривались в работах исследователей [1-10]. Автор [1] отмечает трудность интерпретации и сопоставления получаемых результатов при исследовании электризации вследствие большого разнообразия процессов взаимодействия между контактирующими материальными средами. Такие взаимодействия как трение, сжатие, растяжение, сдвиг и т.д., увеличивают площадь контакта двух сред. А в зависимости от скорости протекания этих процессов усиливают или ослабляют обменные процессы, что способствует увеличению максимальной величины статического заряда, образующегося на контактирующих поверхностях. Поэтому, механические воздействия играют большую роль в процессах электризации, что отмечают все исследователи.
Модель процесса электризации при контакте двух диэлектриков с переходом электро- тава человека в положениях стоя и сидя. Прогнозирование тепло и влагопотерь осуществляется при помощи моделирования процесса теплопередачи с поверхности сложных стенок коленного сустава, при сгибании ноги, геометрически представленной в виде цилиндра и сферы.
Работа [10] посвящена разработке математической модели воздействия ветра на пакет одежды с объемными несвязными утеплителями. Для решения этой задачи ис- нов, где рассмотрены условия перехода заряда пользована цилиндрическая модель туло- при контакте двух диэлектриков с различными энергиями уровня Ферми, рассмотрена в работе [2]. Модель процесса электризации, основанная преимущественно на диффузном переходе ионов с одной поверхности на другую, предложена в работах [3,4].
Автором [5, 6] разработана математическая модель системы «Человек-Специаль- ная защитная одежда от низких температур и статического электричества-Среда», позво- ляющая рассчитывать параметры пакета одежды с антиэлектростатическим эффектом.
Работа [7] посвящена расчету оптимальных параметров толщины теплозащитных слоев одежды на основе оптимизационной математической модели «Человек-Одежда-Среда» и обосновано, что для оптимизационного моделирования целесообразно использовать модель с «овальным» сечением туловища, так как это сводит последующие решения к осе симметричным.
В работе [8] разработана математи- ческая модель распределения температурного поля в слоях «Человек-Одежда-Среда», учитывающая зависимости деформации теплофизических и геометрических характеристик системы в процессе проникновения сырой нефти в пакет материалов одежды.
В работе [9] разработана математическая модель системы «локальная зона человека теплозащитная одежда окружающая среда», дающая возможность прогнозировать тепловое состояние коленного сус-

^Ө-Ц\/^—77iwi--O/ х
вища человека, покрытого равномерным слоем одежды с объемным утеплителем.
Из анализа рассмотренных моделей можно сделать вывод, что большинство моделей явления электризации, предложенных исследователями, касаются теплозащитной одежды для пониженных температур. При математическом моделировании изменения электростатического заряда во времени рассматривался только процесс стекания заряда посредством электропроводимости, а эмиссия электрического заряда с поверхности материала в воздушную прослойку и зависи -мость электрической проводимости материала от его влажности не учитывались.
Материалы и методы исследований
Исходя из вышеизложенного, основной целью данной работы является разра -ботка математической модели антистати -ческой одежды для работы в жарких климатических условиях. В этой связи тело человека рассматривается как деформируемое цилиндрическое тело конечного размера, а одежда как мягкая оболочка (рис 1) Исследуемая задача сводится к моделированию деформаций цилиндра с мягкой оболочкой. Данная постановка задачи позволяет исследовать напряженно-деформированное состояние цилиндра при тепловом воздействии с учетом влияния электростатического поля. Уравнения динамики цилиндра с тонкой скрепленной оболочкой принимает следующий вид
, „ _ л -7=№2/уц где: к-объемный модуль; G-модуль сдвига; Ө -температура;


Рисунок. 1. Расчетная схема для разработки математической модели антистатической одежды:
λ и µ - коэффициенты Ламе; α -коэффициент линейного термического расширения; χ= λ/ c ; c - удельная теплоемкость при постоянной деформации; λ -постоянная величина; Q- величина, которая зависит от количества тепла; γ -плотность на единицу объема материала цилиндрического тела; w(wr, wϕ, wz) - вектор перемещений цилиндра.
Если одежда рассматривается как многослойная цилиндрическая оболочка, тогда дифференциальные уравнения оболочки записываются в виде
cf v 1—VcTv 1+v і?м 5w mH2^' ----- -I----I----1------— I
082 2 do? 2 ОаОб Од S' OP
^M 1—V^U 1+V С? V ^^^ mb'2 ^u "Oa?^ 2 cV+ 2 ОаЗӨ+ Tl?
, ди д\' 1V7№
V=^-a=--,nl где
масса
многослойной оболочки, отнесенная к единице площади нейтральной поверхности (А); b′ -радиус нейтральной поверхности. Система уравнений (2) отличается от соответствующих уравнений оболочки только жесткостями на растяжение и изгиб. Полные уравнения движения оболочки цилиндра в цилиндрической системе координат имеет следующий вид:


Вынужденные уравнения электростатики имеют следующий вид (рис. 2)

Рисунок.2. Цилиндр в электростатическом поле
Электростатическое поле учитивается через граничные условия при z= ± h , которое имеет следующий вид :


Таким образом, математическая постановка модели системы «человек-одежда-среда» дает возможность оптимального выбора определения состояния одежды при внешних воздействиях.
Методика и алгоритм решения задачи деформации одежды ограниченной длины при наличии теплового потока.
Дана одежда ограниченной длины . Площадь поперечного сечения предполагается постоянной по длине. Модуль упругости материала одежды - E , коэффициент теплового расширения материала стержня - а. Одежда с внутренней стороны свободно держится, или оболочечный цилиндр и внутренняя полость тела имеют скользящий контакт. Вместе с этим на поверхности поперечного сечения приложен тепловой поток q(Вт / см2) . Через площадь поперечного сечения происходит теплообмен с окружающей средой. (рис. 1). Требуется по длине тела с заполнителем от поля температуры и электромагнитное поле.
Естественно, возникающие при этом поля температур, перемещений, деформаций и напряжений будут распределены по длине конструкции нелинейным образом. В связи с этим, в малой части длины конструкции, поле распределения физических величин примем в виде кривой второго порядка. Например, рассмотрим поле распределения температуры в интервале 0 < х < £ , где £ << L - конечная длина рассматриваемого частично теплоизолированного стержня; £ -
£ < (ОД см) длина части одежда,
Предположим, что поперечное сечение постоянно по длине стержня, тогда в интервале 0 < х < £ поле распределения температуры Т = Т(х) представим в виде кривой второго порядка, определить поле распределения температуры
, при 0 < х < £ ,
где а,а,а - некоторые константы, значения которых пока неизвестны. Для нахождения значения этих констант, рассматриваемую часть стержня делим пополам. И в этой части фиксируем три узла (рис.
3). Глобальные координаты этих узлов соот- ветственно , при этом

.
i
— |
||
x |
||
j k |
||
l/2 |
Рисунок 3. Квадратичный конечный элемент с тремя узлами
При этом в местной системе координат координаты трёх узлов определяются следующим образом: х = 0;

xk = £ .
Далее введём следующие обозначения, характеризующие значение температуры в узлах i , j , k .
Т1=Т(х = х1У Т Т(х х
У Тк=Т(х = хк)
Таким образом, для первой и второй половины стержня находим значения истинного термоупругого напряжения
45о£ 55ой ІООой
=-25ой;
45а£ 145о£ ЮОай ------------=-------=- 2 5 ой.
Полученные численные результаты (на 100% точно) совпадают с аналитическими решениями, приведенными в работах [1,2]. Эти сравнения показывают, что точность результата метода приведенного универ-
сального вычислительного алгоритма является самой высокой.
В пределах этого элемента апрокси-^тм мируя температуру квадратным полиномом, получим следующее выражение
Г=Г(х) = ^ (х)Г, + ^ (х)^ + % (х)Тк
,
где
- функции формы
для одномерного квадратичного конечного элемента с тремя узлами, и они выражаются
через формулу (8). Тогда в пределах этого элемента градиент температуры определяется следующим образом
^=^^г+^^г+^^г
.
Как видно, здесь надо учитывать удлинение и расширение одежды за счет
температуры.
Дана одежда
ограниченной
кк поперечного сечения правого конца температура окружающей среды
,
.
длины . Площадь поперечного сечения постоянна по длине. Коэффициент теплообмена с окружающей средой через площадь
Требуется определить поле распределения температуры по одежде. Расчетная схема задачи приведена ниже (рис. 4).
hk,T h,T i oci

q
k
i
j
Рисунок. 4. Одежда под влиянием теплового потока и растягивающей силы
Рассмотрим данный стержень как квадратичный конечный элемент с тремя узлами .
Дифференциальное уравнение распределения температуры данного случая будет таким

Т
Отсюда определим
Тк =
^хх^^і^осі ^кТоскУ ^^ ^Кхх*^^*^^Ч^ ^оск
(AA^+AA+W2)^
Из последней формулы видно, что на значение температуры прямо пропорционально значению коэффициента тепло проводности материала одежды, тепловому потоку подведенному к боковой поверхности
одежды. Теперь, подставляя найденное зна чение в выражение (12), найдем значе- т.
ние температуры -^ в точке x = x^ стержня
oti +^ock V^k
^оск
, \ік-ғ<ЛтгМ*к.-ггМ*ҺіМ1')
ядгУ (ДбК^1 ^Х^К^и^К^^+ЗҺ^У)
nmK^uK^^Vb
Далее построен вычислительный алгоритм определения деформационного состояния одежды ограниченной длины при наличии теплового потока и теплообмена. Выявлен соответствующий функционал, выражающий полную тепловую энергию рассматриваемой одежды с учетом заданного теплового потока, теплообмена и растягивающей силы. Минимизируя последние по узловым значениям температур, построена разрешающая система линейных алгебраических уравнений.
Результаты и их обсуждение
Дана одежда ограниченной длины. Площадь поперечного сечения
постоянна нец одеждs жестко конце приложена
по длине. Верхний ко-защемлен, а на нижнем осевая растягивающая
сила
стержня проводности
. Модуль упругости материала
, коэффициент тепло-
к (Вт/\.см
Кхх ,(
расширения-стержня
-
материала стержня
), коэффициент теплового
.
В
задана
верхнем конце температу-
т = ти = ^ = т, ра .
боковой поверхности и
Через
участки
1 — л2 — л3
площади попереч-
ного сечения нижнего конца стержня проис-
ходит теплообмен с окружающими их средами. Здесь температуры окружающей среды т т т обозначим через 1 ociVoc^ocs , а коэффициенты теплообмена с окружающей сре- дой соответственно

Боковая по-
верхность
0 < х < хр х2 < х < х3, х4 < х < Z теплоизолирована участков одежды
(рис.5)

Рисунок. 5. Расчетная схема задачи
Требуется определить величину удлинения одежды от воздействия осевой растя-ГЛЛПЛтіПСЧЛ ! Р Т = Т (Х = 0) = Тд гивающей силы P , и температуры окружающей сре-
I T т т
. Здесь надо учитывать, что
некоторые участки одежды теплоизолированы по боковой поверхности. Для этого рассматриваемый стержень длиной L будем
средой, дискретизируем квадратичными
1 ^ _ х2 - х1 _ х4 —х3
■^4 ^5
элементами с длиной .
Теперь будем исследовать следующие положения, где длина рассматриваемой
L = одежды равна , площадь поперечного сечения стержня ,
-
дискретизировать квадратичными конечны-
ми элементами с тремя узлами.
Теплоизолированную по боковой по-
модуль упругости материала одежды
2 "10 )
, коэффициент тепло-
верхности
0 < х < Хр х2< х < х3, х4< х < L
участки одежды
К — 72 (хВтн / СЛІ С ^) проводности - , коэффи-
дискретизируем на равные элементы. Пусть
элемента
будет
длина каждого

А происходит теплообмен с
циент теплового расширения материала одежда - , х^іб^лі) х2 = 32 (елі) х3 = 48(см)
и, ,
,
х . = 64 (елі )
, x 5 = L .
участки, где окружающей
I-вариант К — ^2 —^з — 10 (Sm / см С), Тос11 — Т0С2 — Тосз — 20 ( С)
Число конечных элементов 800. Число узловых точек 1601. Длина конечных элементов одинакова: Л/ = 80 (см) /800 = 0,1 (см) одинакова:
В таком положении значения температуры в 1601 узловых точках одежды приведены в табл. 1. А поле распределения темпе- сунке 6. В этом рисунке по оси абцисс расположены номера узловых точек одежды, а по оси ординат - значения температуры.
Из рис. 6 и табл. 1 видно, что в этом варианте поле распределения температуры по длине одежды выражается через равномерную непрерывную кривую.
ратуры по длине одежды приведено на ри-
Таблица 1 - Закон распределения поля температур в фиксированных точках стержня при
^oci - ^ос2 - ^осз - 20 ( С)
T1= |
T350= |
T700= |
T1050= |
T1400= |
60,0000000 |
23,9191529 |
20,0457544 |
20,0019617 |
20,0000559 |
T50= |
T400= |
T750= |
T1100= |
T1450= |
54,8458907 |
21,7101294 |
20,0386337 |
20,0008571 |
20,0000478 |
T100= |
T450= |
T800= |
T1150= |
T1500= |
49,5865956 |
20,7468463 |
20,0315131 |
20,0003769 |
20,0000397 |
T150= |
T500= |
T850= |
T1200= |
T1550= |
44,3273004 |
20,3276018 |
20,0243924 |
20,0001714 |
20,0000316 |
T200= |
T550= |
T900= |
T1250= |
T1600= |
39,0680053 |
20,1469982 |
20,0172718 |
20,0000906 |
20,0000235 |
T250= |
T600= |
T950= |
T1300= |
T1601= |
33,8087101 |
20,0734752 |
20,0101511 |
20,0000721 |
20,0000233 |
T300= |
T650= |
T1000= |
T1350= |
|
28,5494150 |
20,0528751 |
20,0044952 |
20,0000640 |
Здесь удлинение одежды от поля распреления температуры по длине одежды
/Vг - ct jT(x )^ - 0,0248791929(см) соответственно [2] будет равно ° .


Длина стержня
Рисунок 8. Закон распределения поля температуры по длине одежды при
Tocv-Toc^Tom=A Длина стержня Рисунок 9. Закон распределения поля температуры по длине одежды при ^оа = Тоса = Тхз = 50( О п ^i — ^2 — ^з — ІО^ли/сл; С),7 — 2ОСЗ — 30( С) ІІ-вариант. Пусть Для этого варианта полученные числовые результаты показаны в таблице 2. Поле распределения температуры по длине одежды выражается через равномерную непрерывную кривую. Оно приведено на рисунках 6,7. А также в этом варианте удлинение стержня от поля температуры будет . равно Д/у —0,0336593947(см) . Сравнивая результаты, полученные в І-ІІ вариантах видим, окружающей среды когда температура увеличивается на 50% , то удлинение одежды увеличивается на 132,84%. Таблица 2 - Закон распределения поля температур в фиксированных точках стержня при T1= 60,0000000 T350= 32,9393647 T700= 30,0343158 T1050= 30,0014712 T1400= T50= 56,1344181 T400= 31,2825971 T750= 30,0289753 T1100= 30,0006428 30,0000419 T100= 52,1899467 T450= 30,5601347 T800= 30,0236348 T1150= 30,0002827 T1450= T150= 48,2454753 T500= 30,2457014 T850= 30,0182943 T1200= 30,0001285 30,0000358 T200= 44,3010040 T550= 30,1102487 T900= 30,0129538 T1250= 30,0000679 T1500= T250= 40,3565326 T600= 30,0551064 T950= 30,0076133 T1300= 30,0000541 30,0000298 T300= 36,4120612 T650= 30,0396563 T1000= 30,0033714 T1350= 30,0000480 T1550= 30,0000237 T1600= 30,0000176 T1601= 30,0000175 IІI-вариант. Пусть ^і - ^2 - ^з - 10(5wj !см CYTq^ - Tom - Тот - 40 ( С) Тогда, сравнивая с I вариантом, видим, что при увеличении температуры окружаю- щей среды в два (pmTQC = гОСОдоГдс =40(‘,С')) раза , удли- Численные результаты, полученные в этом варианте, приведены в табл. 3. Здесь поле распределения температуры по длине стержня показано в виде рис. 8. нение одежды увеличивается на 177,99%. IV-вариант. ПустьҚ. —һ^—һз —10(Sm/см СУуТ^^ —Т^ —^осз — 50( С) Здесь по сравнению с I-вариантом температура окружающей среды ется (сот Тос - 20 (°C) до Тос - 50 (°C)) зультате удлинение одежды Д1= 0,0512197982 (см ) дет . меня В ребу- По сравнению с I-вариантом удлинение одежды увеличивается на 198,51%. Закон распределения поля температуры по длине одежды приводится в виде табл. 4 и рис. 9. Таблица 3 - Закон распределения поля температур в фиксированных точках одежды при Тос,=Тоа=Тоа=Д^С) T1= 60,0000000 T350= 41,9595765 T700= 40,0228772 T1050= 40,0009808 T1400= T50= 57,4229454 T400= 40,8550647 T750= 40,0193169 T1100= 40,0004286 40,0000279 T100= 54,7932978 T450= 40,3734231 T800= 40,0157565 T1150= 40,0001885 T1460= T150= 52,1636502 T500= 40,1638009 T850= 40,0121962 T1200= 40,0000857 40,0000231 T200= 49,5340026 T550= 40,0734991 T900= 40,0086359 T1250= 40,0000453 T1500= T250= 46,9043551 T600= 40,0367376 T950= 40,0050756 T1300= 40,0000361 40,0000198 T300= 44,2747075 T650= 40,0264375 T1000= 40,0022476 T1350= 40,0000320 T1550= 40,0000158 T1600= 40,0000117 T1601= 40,0000117 Таблица 4 - Закон распределения поля температур в фиксированных точках одежды при Tqci =^oa = T0C3 = 50 ( C) T1= T350= 50,9797882 T700= 50,0114386 T1050= 50,0004904 T1400= 60,0000000 T400= 50,4275324 T750= 50,0096584 T1100= 50,0002143 50,0000140 T50= T450= 50,1867116 T800= 50,0078783 T1150= 50,0000942 T1450= 58,7114727 T500= 50,0819005 T850= 50,0060981 T1200= 50,0000428 50,0000119 T100= T550= 50,0367496 T900= 50,0043179 T1250= 50,0000226 T1500= 57,3966489 T600= 50,0183688 T950= 50,0025378 T1300= 50,0000180 50,0000099 T150= 56,0818251 T200= 54,7670013 T250= 53,4521775 T300= 52,1373537 T650= 50,0132188 T1000= 50,0011238 T1350= 50,0000160 T1550= 50,0000079 T1600= 50,0000059 T1601= 50,0000058 Заключение, выводы Таким образом, на основе анализа существующих математических моделей электризации защитной одежды определены критерии для разработки новой математической модели системы «Человек - Одежда – Среда», позволяющие описать процессы электризации в одежде в жарких климатических условиях. На основе полученных дифференциальных уравнений были найдены значения температуры для специальной одежды ограниченной длины в узловых точках. При этом, количество узловых точек равно 1601, количество связанных элементов равно 800. В области распределения температуры по длине одежды в четырех различных вариантах было обнаружено, что в I варианте AlT0,0249(см), во ii варианте A/r = 0,0337(см), в in варианте AlT = 0,0413( см), в IV варианте Мт =0,0512 (си) л г T , см , в V варианте AlT =0,0600(см), в ух варианте длина одежды в этой области увеличилась на 10% по сравнению с длиной одежды во II – варианте. Разработанная математическая модель описывает процессы электризации в поперечных сечениях пакета одежды и на поверхностном слое кожи человека и позволяет рассчитать параметры антистатической одежды для защиты от повышенных температур. Предложенная математическая модель позволяет определить условия проведения экспериментов, оптимизировать основные параметры, а также оценить эффективность подбора материалов и конструктивных параметров одежды на этапе проектирования. На основе разработанной математической модели была изготовлена экспериментально-опытная партия специальной антистатической одежды на предприятиях ООО " Al-Hakim Plus" и ООО "Okean-Lyuхе"(Бухарская область). Разработанная новая антистатическая спецодежда была положительно оценена потребителями данной спецодежды - сотрудниками ООО "Региональные электрические сети» Бухарской области и специалистами швейных предприятий. Годовой экономический эффект от внедрения в производство новой спецодежды составил 309 060 000 рублей.
Список литературы Постановка математической задачи и методы решения разработки антистатической одежды для защиты от повышенных температур
- Лёб, Л. Статическая электризация - М.-Л.: Госэнергоиздат, 1963.-408 с.
- Ландау Л.Д., Лифшиц Е.М Механика. Электростатика. Краткий курс теоретической физики - М.: Наука, 1969.-272 с.
- Меркулова А.В., Исследование и разработка специальной антиэлектростатической одежды для защиты от пониженных температур. Дисс...канд. техн. наук по спец. 05.19.04. Шахты, Юж.-Рос. гос. ун-т экон. и сервиса 2007. -184 с.
- Незгода, В.Ю. Изменение влажности, температуры и образование электростатических зарядов на хлопчатобумажной ленте при ее вытягивании // Текстильная промышленность, №2.- 1988. - С.42-43
- Черунова, И.В. Совершенствование методов проектирования специальной одежды для горноспасателей: дис....канд. техн. наук: 05.19.04 / И.В. Черунова. - М., 2001 - 193 с. библиогр.: с. 183-193
- Черунова, И.В. Теоретические основы комплексного проектирования специальной теплозащитной одежды: дис.... докт. техн. наук: 05.19.04. - Шахты: Юж.-Рос. гос. ун-т экон. и сервиса., 2008. - 394 с.
- Кудрявцев, В. И. Усовершенствованная технология проектирования теплозащитной одежды на основе уточненных моделей теплообмена. Автореферат дис....канд. техн. наук: 05.19.04 / Новочеркасск -2004.- 24с.
- Куренова, И. В. Разработка и исследование специальной нефтезащитной одежды с модифицированным пакетом материалов. автореферат дис.... канд. техн. наук: 05.19.04 / Юж.-Рос. гос. ун-т экон. и сервиса. - Шахты, 2013 год,- 18 с.
- Герасименко, М.С Разработка методики оптимизации параметров специальной теплозащитной одежды для астрономов, работающих в условиях гиподинамии: автореферат дис.... канд. техн. наук: 05.19.04 / Юж.-Рос. гос. ун-т экон. и сервиса. - Шахты, 2005. - 16 с.
- Лебедева, Е. О. Исследование и разработка пакета специальной теплозащитной одежды с повышенной устойчивостью к ветру: автореферат дис. … канд. техн. наук: 05.19.04 / Юж.-Рос. гос. ун-т экон. и сервиса. - Шахты, 2006. - 24 с.