Моделирование теплового процесса при сварке полиэтиленовых труб при естественно низких температурах Текст научной статьи по специальности «Механика и машиностроение»

УДК 539.3

О. А. Аммосова, Р. А. Каюмов, И. З. Мухамедова МОДЕЛИРОВАНИЕ ТЕПЛОВОГО ПРОЦЕССА ПРИ СВАРКЕ ПОЛИЭТИЛЕНОВЫХ ТРУБ ПРИ ЕСТЕСТВЕННО НИЗКИХ ТЕМПЕРАТУРАХ

Ключевые слова: сварка полиэтиленовых труб, тепловой процесс, предварительный подогрев, грат, сварное соединение,

скорость охлаждения, теплоизоляционная камера.

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

Keywords: welding ofpolyethylene pipes, heat process, preheating, beads, weld, the rate of cooling, thermo camera.

A mathematical model of heat welding of polyethylene pipes taking into account the phase transition in the temperature range is developed. The method of the calculation of Length ofpreheating of the heating at welding of the polymeric pipes at low temperature of the air is suggested. The mathematical model of the heat process of the welding of the polymeric pipes, taking into account forming and heat influence of a burr on the temperature field dynamics is suggested. Possibility of the regulation to velocities of the cooling by variation size heat-insulating cameras is shown.

Введение

В строительстве газораспределительных сетей полиэтиленовые трубы (ПЭ) заняли уже довольно прочные позиции. В тоже время их применение в северных регионах представляется наиболее проблематичным и сдерживается, главным образом, действующими нормативными документами. Прочность сварного соединения во многом определяется температурным режимом в процессе сварки. Температурное поле в конце нагрева и скорость охлаждения свариваемых труб на этапе осадки оказывают существенное влияние на формирование надмолекулярной структуры материала сварного шва и, соответственно, прочность соединения. При сварке полиэтиленовых труб в условиях низких температур окружающего воздуха (ОВ) динамика температурного поля не обеспечивает необходимую прочность сварного шва. Для расширения диапазона температур ОВ, при котором допускается проведение сварки полиэтиленовых труб, практически не используются методы управления динамикой температурного поля. Регулирование температурным режимом и обеспечение такого же температурного поля при оплавлении и такого же темпа охлаждения при осадке, как и при допустимых температурах ОВ позволит провести сварку полимерных труб при температурах ОВ ниже нормативных без изменения основных параметров сварки (давления на торцы, технологической паузы, величины и скорости осадки).

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

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

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

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

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

« (т » >1 = 1 f h <T >f К (>•' T >f ) • <■)

0 < t < tm; 0 < r < rr; 0 < z < l, где T(r,z,t) - температура в момент времени t; с(Т)-коэффициент теплоемкости; р(Т) - удельная плотность; Х(Т) - коэффициент теплопроводности; tm -время расчета; r, z - цилиндрические координаты; Г - свободная поверхность грата и части трубы. Индекс i=l для твердой фазы материала трубы; i=2 для жидкой фазы материала трубы; i=3 для воздуха.

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

Т (r,z,0> = То . (2)

На одном торце задается температура нагретого инструмента:

T(r,0,t) = тн . (3)

Учитывая низкую теплопроводность полиэтилена, примем, что на некотором удалении ъ = 1 от зоны сварки температура трубы не изменяется на протяжении всего времени протекания процесса. Таким образом, на другом торце трубы задается условие:

Т (г.Н ) = То (4)

На внешней поверхности трубы происходит конвективный теплообмен с окружающей сре-

„ , дГ дои: X,—

дг

= -а

Г

(Tl Г- T0 )

(5)

При r = 0 условие ограниченности решения:

lim Xr — = 0 (6)

r ^0 dr

На границе фазового перехода записывается условие Стефана:

0(t)+0 ~ Tl<!

= T

\<b(t)-0 Ф-

дФ

(XigradT - X2gradT,gradФ) - Lp-^ = 0 •

(7)

где Ф(г,ъД) - уравнение положения границы раздела фаз в момент времени 1; ТФ - температура фазового перехода; Ь - удельная теплота фазового перехода.

В полиэтилене не существует такой четко выраженной границы, фазовый переход происходит в интервале температур. В случае фазового перехода в интервале температур в математической модели необходимо учитывать промежуточную фазу между твердым и жидким веществом [2]. Границы промежуточной фазы определяются температурами соли-дуса Т и ликвидуса Ть . В работе [2, 3] температурное поле в этом случае предлагается определять из уравнения теплопроводности с эффективным коэффициентом теплоемкости:

C (T )p(T )дт =iA (rA(T )дТ) + —( X(T )дТ V ; v ' dt r дг дг dz I V dz

(8)

где

0 < t < tm

С (T ) =

0 < r < rr

Cr, T < TL

0 < z < l,

q - l

d У

dT '

tl < T < tS

(9)

C2, T > Ts

-2> 1 - ^

Т(Т)- объемная доля твердой фазы. При предположении прямой пропорциональности объема расплавленного полиэтилена количеству выделившейся теплоты, *¥(Т) определяется по формуле:

y(T ) = 1 --

J q(u)du

Ts_

Tl

J q (udu )

(10)

где q (Т) - зависимость теплового потока от температуры, регистрируемая дифференциальным сканирующим калориметром, отнесенная к единице массы вещества.

Задача (1)-(10) решалась численно методом

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

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

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

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

Во-вторых, согласно расчетам, приведенным в работе [8], напряжения, обусловленные сваркой, повышаются с понижением температуры окру-

S

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

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

Расчетная схема стенки трубы с нагревателем представлена на рис. 1. Предполагая однородность температурного поля по окружности полиэтиленовой трубы, рассмотрим отрезок его стенки длиной N с нагревателем и эластичным элементом толщиной N1, обеспечивающей идеальный тепловой контакт соприкасающихся поверхностей. Нагреватель обеспечивает постоянную температуру Тн на внешней поверхности эластичного элемента, поэтому изменение температуры в нагревателе не рассчитывается. Значение Тн меньше температуры плавления полиэтилена.

Рис. 1 - Схема стенки трубы 1 с нагревателем 3 и эластичным элементом 2

В области в, ограниченной свободными поверхностями стенки трубы и эластичного элемента запишем уравнение теплопроводности:

дТ 1 5

(11)

dt

„ dT Л d L dT rX,— |+—I X,— r dr l dr I dz l dz

т,х е О, 0 < г < tm На внешней границе эластичного элемента Ъ задана постоянная температура:

Т()\г = Т2, 0 < г < т (12)

На свободных поверхностях стенки и эластичного элемента С происходит конвективный теплообмен с окружающей средой с температурой Т0:

X —

, dn

= -а(Т (г,2,г)с - То), о < г < гт (13)

с

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

T (r,N + N1 ,t ) = T0 (14)

В начальный момент времени в стенке трубы и эластичном элементе распределение температуры однородно и равно температуре ОВ:

T (r,z,0)= To (15)

Индекс i=1 для полиэтилена, i=4 для эластичного материала. Задача решалась методом конечных разностей сведением задачи к цепочке одномерных и использованием однородной разностной схемы со вторым порядком аппроксимации по пространственным переменным. Граничные условия вида (13) аппроксимировали со вторым порядком.

Численное моделирование проводилось для трубы ПЭ 100 ГАЗ SDR 9 110x12,3. Расчеты проводились при следующих данных: r1 = 0,0416; r2=0,055 м; Х1=0,38; Х2=0,29 Вт/(кг-К); р1=954; р2=700 кг/м3; с1=1900; с2=1900 Дж/(кг-К). Длина подогреваемой области нагревательного инструмента равнялась 6,5 см. Температура окружающего воздуха минус 40 °С.

При подогреве инструментом через определенное время температура в стенке трубы становится равной температуре инструмента по толщине и по длине части трубы. Для удаления инструмента предварительного подогрева и приведения труб в контакт с другим нагревательным инструментом необходимо некоторое время, за которое прогретая часть остывает. Поэтому температура инструмента для подогрева специально выбирается 60 °С, что выше допустимой для сварки температуры ОВ, но ниже температуры размягчения. Продолжительность подогрева выбирается таким образом, чтобы распределение температуры стенки трубы стало неоднородным по радиальной переменной. Такой прием позволяет сделать достаточно длительную паузу для обеспечения более однородного распределения температуры (выравнивания) в допустимом для сварки диапазоне.

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

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

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

Как видно из расчетов за температуру окружающего воздуха можно принять температуру равную 20 °С. Тогда дальнейший технологический процесс сварки ПЭ труб происходит как при допустимых температурах. Согласно [1] продолжительность нагрева торцов труб при температуре ОВ равной 20 °С составляет 100 с. т.°с

50

40

го

у

\\

ч /

- ч - - _ Г т . 1

■ — ■

41.6 43,3

45

45,6 48,3

50

51,7 53.3 55 I, 1,11,1

Рис. 2 - Распределение температуры в стенке трубы по радиусу:

1-3 -соответственно расстояние от торца 0,01; 0,03; 0,06 м;

1-3 - в момент завершения подогрева (90 с); 1'—3' - в момент завершения выравнивания (190 с)

Анализ распределения температуры по длине трубы в некотором радиусе (рис. 3) показал, что при проведение предварительного подогрева температурное поле (кривая 1) в зоне сварки становится идентичным температурному полю при температуре ОВ 20 °С (кривая 2). Как видно из рис. 3 при проведении сварки по стандартной технологии при температуре ОВ -40 °С без предварительного подогрева (кривая 3) температурное поле в трубе не достигает необходимого уровня.

Т,°С

250

200

150

50

-50

1

\

.... ж

\

\

0,32

0.04

0.05

0,0В г, мм

Рис. 3 - Распределение температуры по длине трубы в момент завершения этапа нагрева (100 с), радиус трубы 0,0483 м: 1- Т0 = -40 °С с подогревом, 2 - Т0 = 20 °С, 3 - Т0 = - 40 °С без подогрева

Охлаждение сварного соединения в теплоизоляционной камере. Определение размера теплоизоляционной камеры

После оплавления (нагрева) кромок свариваемых труб следует этап осадки. Этап осадки наи-

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

Осадка при температурах воздуха ниже регламентируемых приводит к возрастанию скорости охлаждения сварного соединения. Слишком быстрое охлаждение вызывает в зоне сварки значительные внутренние напряжения, которые не успевают сглаживаться до того, как к швам прикладывается дополнительные монтажные нагрузки [9]. Также при слишком быстром охлаждении процесс кристаллизации может оборваться на промежуточной стадии, тогда как при формировании структуры в условиях медленного охлаждения разрушившиеся центры кристаллизации успевают восстанавливаться [10]. Поэтому регулирование скорости охлаждения при температурах воздуха ниже регламентируемых является актуальной проблемой.

Пренебрегая продолжительностью технологической паузы и временем, затрачиваемым на течение расплава при осадке, труба считается укороченной на величину осадки. Для упрощения расчетов принимается допущение, что поперечное сечение образующегося грата имеет прямоугольную форму. Для величины осадки 2,6 мм поперечное сечение бралось в форме прямоугольника со сторонами 3,3 и 5,3 мм. Начало координат по длине трубы ъ смещалось на величину осадки и ставилось на месте стыка труб.

В начале процесса охлаждения распределение температуры в грате предполагается равномерным. Значение температуры в грате определяется из условия теплового баланса - равенства количеств теплоты в области уходящей в грат и в образовавшемся грате. Таким образом, распределение температуры в трубе с гратом в начале процесса охлаждения получается в результате определения температурного поля в конце нагрева. Несмотря на то, что в расчетах пренебрегается продолжительность технологической паузы, предполагается, что при отделении нагревательного инструмента от свариваемых труб и смыкании оплавленных деталей температура воздуха внутри трубы становится равной температуре окружающего воздуха и вновь повышается за счет теплоотдачи от нагретой части трубы на стадии охлаждения. При этом постановка задачи аналогична этапу нагрева, но условие (3) на одном торце изменяется на условие отсутствия теплового потока (условие симметрии температурного поля) в режиме охлаждения:

дГ_

аГ

= 0.

(16)

я=0

На свободной поверхности грата задается условие конвективного теплообмена с окружающим воздухом.

Для расчета использовались те же обозначения и исходные данные, что и для нагрева.

Так как скорость охлаждения увеличивается с понижением температуры ОВ, предлагается использование цилиндрической камеры. Размеры камеры определяются на основе теоретического моделирования процесса охлаждения сварного соединения с теплоизоляционной камерой из условия обеспечения допустимой скорости охлаждения. Примем следующие допущения: предполагается, что температура воздуха внутри камеры достаточно быстро (в пределах расчетного шага по времени) становится однородной, камера может быть изготовлена из тонкой пленки с очень низкой теплопроводностью (X=0,0001), тогда толщиной стенки камеры можно пренебречь.

Учитывая эти допущения температуру в камере находим из соотношений:

Q = T|к-Tkam (tК, (17)

Q = c3P3Vkam ^Tkam - Tkam | , (18)

v

где Vkam - объем камеры, T - температура в камере на нижнем временном шаге, Т - температура в камере, ak - коэффициент теплообмена в камере, К -внешняя поверхность грата и части трубы, находящейся внутри теплоизоляционной камеры.

Осредненная температура в камере, полученная в результате расчета, зависит от геометрических размеров теплоизоляционной камеры. Расчеты показывают, что увеличение длины камеры при фиксированной высоте способствует снижению температуры в камере, что соответствует физическому представлению. При этом увеличивается не только объем камеры, но и увеличивается поверхность трубы с более низкой температурой, что приводит к снижению температуры в теплоизоляционной камере. Однако, увеличение высоты (радиуса) камеры при фиксированной длине также способствует снижению температуры воздуха внутри камеры, но не столь значительно. Поэтому, принимая во внимание возможности изготовления и практическое использование теплоизоляционной камеры, один из геометрических размеров, например, высоту можно зафиксировать и, варьируя длиной камеры, определить размер камеры, обеспечивающий допустимую скорость охлаждения сварного соединения. При высоте камеры равной 0,03 м расчетами найдена полудлина камеры L/2 равная 0,03 м, обеспечивающая допустимую скорость охлаждения в интервале низких температур ОВ (-45, -15) °С.

Характерные зависимости температуры от времени при охлаждении сварного соединения в теплоизоляционной камере при различных температурах ОВ приведены на рис. 4. Для сравнения показаны полученные расчетами предельные кривые изменения температур по времени при темпера-

турах ОВ -15 и 45 °С (кривые 2, 1 соответственно) при охлаждении соединения в естественных условиях, а также изменение температуры по времени при температуре ОВ -40 °С при естественном охлаждении (кривая 5). Зависимость температуры по времени при охлаждении соединения в теплоизоляционной камере лежит между предельными кривыми.

т,°с

' \ \ 1

г\; \

\

X /2

■

■ ,

О 120 240 360 4В0 600 720 840 Э50

X. С

Рис. 4 - Изменение температуры по времени в точке (0,0553; 0,001)при различных температурах ОВ, (1-3, 5) - без камеры, 4 - с камерой: 1- 45 °С; 2 - -15 °С; 3 - 20 °С; 4 - -40 °С с подогревом; 5 - -40°С без подогрева

Заключение

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

Работа выполнена при поддержке Министерства образования и науки РФ в рамках ФЦП «Научные и научно-педагогические кадры инновационной России» на 20092013 годы (соглашение 14.B37.21.2018).

Литература

1. СП 42-103-2003. Проектирование и строительство газопроводов из полиэтиленовых труб и реконструкция изношенных газопроводов, М.: Полимергаз, ФГУП ЦПП, 2004. - 86 с.

2. Авдонин Н.А. Математическое описание процессов кристаллизации. - Рига: Зинатне, 1980. - 180 с.

3. Вабищевич П.Н. Численные методы решения задач со свободной границей. - М.: Изд-во МГУ, 1987.- 164 с.

4. Самарский А.А. Теория разностных схем. - М.: Наука, 1977. - 656 с.

5. Старостин Н.П. Контактная сварка полимерных труб оплавлением при низких температурах окружающей среды. 1. Математическое моделирование теплового процесса / Н.П. Старостин, О.А. Аммосова // Сварочное

производство. 2007. - № 4. - С.17-20.

6. Старостин Н.П. Контактная сварка полиэтиленовых труб оплавлением при низких температурах окружающей среды. Ч.2. Исследование процесса охлаждения / Н.П. Старостин, О.А. Аммосова // Сварочное производство. 2008. - № 9. - С. 31-34.

7. Гринюк В.Д. Микроструктура и качество стыковых сварных соединений полиэтилена / В.Д.Гринюк, А.А.Шадрин, А.В.Золотарь // Автоматическая сварка. -1990. - № 11. - (452). С. 23-26.

8. Нестеренко Н. П. Моделирование температурных полей и напряжений в полиэтиленовых трубах при сварке нагретым инструментом / Н.П. Нестеренко, И.К. Сен-ченков, О.П.Червинко, М.Г. Менжерес // Автоматическая сварка. 2009. - № 2. - С. 11-15.

9. Каргин В.Ю., Шурайц А.Л. Сварка и контроль газопроводов из полимерных материалов. В помощь сварщикам и специалистам сварочного производства. - М.: ОАО «Приволжск. кн. изд-во», 2003. - 330 с.

10. Крюкова И.М. Влияние температуры расплава на свойства полиэтилена в крупногабаритных изделиях / И.М.Крюкова, И.И.Сквирская, В.Я.Ушаков, Б.В.Шмаков

// Пластические массы. 1998. - № 6. - С. 38-39.

11. Каюмов Р.А. Моделирование поведения пленочно-тканевого материала при воздействии эксплуатационных факторов/ Р. А. Каюмов, А.М. Сулейманов, И.З.Мухамедова // Механика композиционных материалов и конструкций. 2005. - т.11. - №4. - С.519-530.

12. Каюмов, Р.А. Физически нелинейная модель тканевых материалов, обработанных в потоке высокочастотной плазма/ Р.А.Каюмов, И.Ш.Абдуллин, М.Н.Серазутдинов, В.В.Хамматова, И.З.Мухамедова, Е.В.Кумпан, А.М.Шакирова // Вестник Казанского технологического университета. №6; Федеральное агентство по образованию, Казан. гос. технол. ун-т. - Казань: КГТУ, 2010. - С. 42 - 53.

13. Каюмов, Р.А. Математическое моделирование деформирования текстильных материалов с содержанием полимерных волокон /Р.А.Каюмов, В.В.Хамматова// Вестник Казанского технологического университета. № 14; Федеральное агентство по образованию, Казан. гос. тех-нол. ун-т. - Казань: КГТУ, 2012. - С. 154 - 156.

© О. А. Аммосова- к.т.н, научный сотрудник ФГБУН Институт проблем нефти и газа Сибирского отделения Российской академии наук, ammosova_o@mail.ru, И. З. Мухамедова - к.ф.-м.н., доц. каф. сопротивления материалов и основ теории упругости КГАСУ, muhamedova-inzilija@mail.ru; Р. А. Каюмов - д.ф.-м.н., проф., зав. каф. сопротивления материалов и основ теории упругости КГАСУ, проф. каф. дизайна КНИТУ, kayumov@rambler.ru.