CC BY
63
УДК 539.3
DOI: 10.17586/0021-3454-2016-59-8-657-663
МОДЕЛИРОВАНИЕ ПРОЦЕССА НАМАТЫВАНИЯ НИТЕЙ НА ЖЕСТКИЙ ЦИЛИНДР
В. М. Мусалимов1, Ю. С. Монахов1, А. Ю. Кутьин2, Г. А. Соловьева1
1 Университет ИТМО, 197101, Санкт-Петербург, Россия E-mail: musvm@ya.ru
2Производственная компания „Индустрия цвета", 196084, Санкт-Петербург, Россия
Изложены положения теории деформирования спирально-анизотропного тела, которая применяется при определении параметров, формирующих напряженно-деформированное состояние таких тел. Приведены примеры применения данной теории для кабельных конструкций и других тел, обладающих спиральной анизотропией. Особое внимание уделено моделированию процесса наматывания длинномерных материалов. Рассматриваемая модель учитывает взаимосвязь модуля упругости и коэффициента Пуассона спирально-анизотропных материалов, физико-механические и кинематические параметры нити и системы намотки. Проведенный таким образом анализ состояния тела намотки позволяет обеспечить выбор оптимального режима наматывания нитей на жесткое основание.
Ключевые слова: спирально-анизотропный стержень, композиционные материалы, окружная намотка, натяжение нити, упругие постоянные
Все существующие в природе твердые тела можно разделить на две большие группы: изотропные и анизотропные. Изотропные тела имеют одинаковые свойства во всех направлениях, а анизотропные — различные свойства в разных направлениях. Задачи теории упругости для случаев прямолинейной анизотропии рассматривались в работах [1—5] и др. Исследование механики тел, обладающих спиральной анизотропией, началось относительно недавно. Сам термин „спирально-анизотропный стержень" (САС) впервые был введен в работе [6], посвященной задаче о растяжении и кручении сильно закрученной пряжи, которая моделировалась упругим спирально-анизотропным цилиндром.
Различные аспекты механики спирально-анизотропного тела для решения практических задач рассматривались многими исследователями [7—9].
Одним из приложений теории спирально-анизотропного тела к решению практических задач является расчет кабельных конструкций. В работах [10—12] представлена модель спирально-анизотропной среды, которая может быть использована при описании характера деформирования гибких кабелей под действием растягивающей нагрузки, в частности для определения модуля упругости второго рода; выведены уравнения, связывающие внешние силы и моменты с линейными и угловыми деформациями САС; приведены результаты разработки экспериментально-аналитического метода вероятностной оценки модулей спирально-анизотропных стержней. Впоследствии эта модель была применена для построения динамической модели магнитостатического вариометра, чувствительный элемент которого может быть представлен в виде спирально-анизотропного цилиндра [13]. Также эта модель была использована при проектировании устройства для изготовления чувствительных элементов [14], в частности, были разработаны оптимизационные методы оценки физико-механических свойств и упругих характеристик готовых образцов [15, 16].
Теория деформации спирально-анизотропного тела [10, 11] может применяться, в том числе, и для моделирования процессов наматывания длинномерных материалов. Для
корректной регулировки параметров намотки нити на жесткое цилиндрическое основание необходимо учитывать давление, возникающее на границе слоев нитей. Схема формирования слоя нитей представлена на рис. 1, где Ы\ и Ы2 — давление на соответствующие слои намотки; Ль Л2 — толщина 1-го и 2-го слоев; Т — натяжение нити; ю — скорость вращения оправки; с1 — диаметр оправки.
ю
Рис. 1
При давлении на границе первого слоя и основания следует учитывать физико-механические характеристики нити, а также кинематические параметры процесса наматывания, а именно: натяжение нити, скорость вращения оправки, угол намотки (а), модуль упругости нити (E), модуль сдвига нити (G) и коэффициент Пуассона (v).
Определение модуля упругости нитей производится экспериментально с использованием специальных разрывных машин и по определенным методикам (например, ISO 6939:1988, ГОСТ 6611.2-73). Для минимизации напряжений, возникающих в захватах машины при проведении испытаний, используются специальные кнехтовые захваты [17].
Для определения модуля сдвига нитей и коэффициента Пуассона, что является достаточно сложной задачей, используют экспериментальные зависимости. Один из подходов к определению данных величин — представление нити в виде спирально-анизотропного стержня [15]. Для изотропных тел справедливо соотношение
G =
E
2 (1 + v) '
(1)
Для анизотропных тел это соотношение не выполняется, так как Е, 0, V — независимые параметры. Однако, принимая для анизотропных тел 20/Е = т , можно вывести зависимость т = / (V) [16] и таким образом определить неизвестные характеристики, отражающие упругие свойства тела.
Давление на границе первого слоя нитей и жесткого цилиндрического основания можно представить как функцию, зависящую от указанных выше величин:
N1 = / (Т, ш, а, Е, 0, V). (2)
Выделим в первом сформированном слое намотки бесконечно малый элемент АВ (рис. 2). Окружные (ое) и радиальные (ог) напряжения считаем положительными, поскольку они вызывают растяжение материала. Кроме того, полагаем, что при деформации отсутствует осевое смещение витков нитей (актуально для фланцевой катушки, ткацкого
навоя, сновального вала и т.п.), распределение напряжений симметрично относительно оси г и зависит только от переменной г.
<Jr+d<Jr
> -А йг {- ■ ^
. иг ,
г \ . I Ме
и-/
V.
Рис. 2
Для выделенного элемента получено уравнение равновесия в цилиндрических коорди-
натах:
иг -ие
й иг йг
= 0.
(3)
Решение данного уравнения, известное как решение задачи Ламе для плоской деформации [17], имеет следующий вид:
а г = а + ь\, (4)
г
здесь р1 > г > р0, где р0 — радиус цилиндрического основания, р1 — радиус первого слоя намотки; а и Ь — постоянные, полученные при интегрировании уравнения (3), их численные значения определяются граничными условиями на внутренней и наружной поверхностях слоя намотки.
При г = р0 напряжение ог слоя намотки достигнет максимального значения, равного N1, но с противоположным знаком. Появление знака „минус" объясняется выбором в качестве положительного растягивающего напряжения (см. рис. 2). В противном случае при г = р1 напряжение ог = 0. С учетом этого получим формулы для определения радиального ог и окружного ое напряжений в первом слое намотки в зависимости от его давления:
иг =
(
Р1
р2
1 -
2 ^ г2 J
ад =
(
Р1
р2
1+
2 ^ г2 J
(5)
В то же время следует заметить, что сформированный слой намотки рассматривается как спирально-анизотропное тело. В случае криволинейной анизотропии, как отмечено в работе [2], целесообразно закон Гука записывать в системе криволинейных ортогональных координат, выбранной так, чтобы координатные направления в каждой точке совпадали с эквивалентными (в отношении упругих свойств). Рассматриваемый случай спиральной анизотропии позволяет связать со слоем намотки геликоидальную ортогональную систему координат (рис. 3), обладающую отмеченным выше качеством, а именно — совпадением координатных линий с упругоэквивалентными направлениями [11].
Компоненты тензора напряжений в двух системах координат связаны следующими соотношениями (рис. 4) [11]:
r=const
ог =ог,
• 2 2 а^=ае sin а + 2t6z sinа + аz cos а,
2 • 2 ал = Oq cos а - 2i6z sina+oz sin а,
xr5 = = xrQ sin а + xrz cos а,
тгл = тцг = -Tr0 cos а + xrz slna,
= = (а z - Oq ) sin а cos а + i6z (sin2 а - cos2 a);
ar =ar,
2 2 Oq=o5 sin а- 2х5л sin а cos а + ал cos а,
л
2 • 2 а z =а^ cos а + 2т^л sin а cos а + ал sin а,
xrQ = TQr = sin а - Tr^ cos а, xrz = xzr = cos а - sin а,
^ez = ^ze = К " ал)sina cosa + Чц (sin2 a - cos" a
I
z ■
z=const Л 65 Л eQ
N. Л е С.^ Vr
Q У
Рис. 3
\ r
4 '
4
—--v— x I
la
n
Рис. 4
/
■o5
(6)
(7)
Аналогичные системе (6) зависимости существуют между компонентами тензора деформаций в цилиндрической и геликоидальных системах координат:
sr sr,
• 2 2 s^=8q sin a-yQz sin а cos а + sz cos а,
2 -2 s^=8q cos a + yQz sin а cos а + sz sin а,
■ 2~ —4 )■
(8)
У = y^r = УrQ slna - уrz cos а,
Угл = Улг = -УrQ cosа - Уrz sinа
У5л = Ул5 = 2(sz -sQ)sinа cosa + yQz (sin" a-cos" a) По условиям задачи TrQ = Trz = Tqz = 0 , yrQ = jrz = Yqz = 0 и s z = 0 , поэтому при совместном решении 2-го и 3-го уравнений системы (6) получим взаимосвязь между отдельными компонентами тензора напряжений в двух рассматриваемых системах координат:
1
O =^ctg a-aQ
sin2 а
- 2 ■
(9)
Из последнего уравнения системы (7) получим соотношение
, .sin а cos а
^л = (ал )_
sin2 а- cos2 а
(10)
а из последнего уравнения системы (8) — выражение
= -2se sin a cos а = -2(s^ + 8Л )sin а cos а . (11)
Выражения (9)—(11) подставим в уравнение обобщенного закона Гука для спирально-анизотропной среды
= 1
У^Л = G ,
(12)
связывающее деформацию сдвига с напряжением, заменяя при этом в уравнении (11) деформации на соответствующие им напряжения. В результате с учетом выражения (9) найдем один из компонентов тензора напряжений в геликоидальной системе координат:
С
С1
— о„
С
Ст
(14)
где
С1 = ctg а + hp
V—1 1 11 22
---, hp =---, G„ = G(sin a —cos а),
' Ev Gn p Gn E n V '
21 С2 = (cosec а — 2)hp, С3 = —.
Ev
На основе 3-го уравнения системы (6) и выражения (14) определим нормальный компонент тензора напряжений в цилиндрической системе координат, выраженный через параметры спирально-анизотропной среды намотки:
z =
-ае cos2 а sin2 а
- C4G0 - С5аг,
(15)
где
^ С2 2 . 2 ^ С3 2
С4 = —cosec а-ctg а , С5 = cosec а .
С
С,
-1 м
Для того чтобы получить взаимосвязь напряжения и осевой нагрузки Тф, прикладываемой к слою намотки со стороны фланцев, необходимо воспользоваться известными формулами [17]:
Р1
Тф = = 2тйт ,
£ Ро
где £ — площадь поперечного сечения слоя намотки.
С учетом (15) и (5) выражение (16) примет следующий вид:
тф = ф(Ро, Р1, а>Е, О) N1,
(16)
где
Ф = л(С4 -С5)р2 + 2л(С4 + С5)р2
р2
2 2 Р1 -ро
ln
'рО
Ро
Соответственно, зная величину осевой нагрузки Тф, можно определить физико-механические характеристики нити. При наматывании каждого последующего слоя решается аналогичная задача о деформации спирально-анизотропного слоя намотки при наличии спирально-анизотропного заполнителя с учетом новых граничных условий [11].
На основе представленного анализа можно осуществить выбор оптимального режима наматывания нитей на жесткое основание.
СПИСОК ЛИТЕРАТУРЫ
1. Лехницкий С. Г. Кручение анизотропных и неоднородных стержней. М.: Наука, 1971. 310 с.
2. Лехницкий С. Г. Теория упругости анизотропного тела. М.: Наука, 1977. 415 с.
3. Житков П. Н. Плоская задача теории упругости неоднородного ортотропного тела в полярных координатах // Тр. Воронежского гос. ун-та. 1954. Т. 27. С. 20—29.
4. Амбарцумян С. А. Общая теория анизотропных оболочек. М.: Наука, 1974. 324 с.
5. Ашкенази Е. К. Анизотропия машиностроительных материалов. Л.: Машиностроение, 1969. 240 с.
6. Thwaits J. J. The elastic deformation of a rod with helical anisotropy // Intern. J. Mech. 1977. Vol. 19, N 3. P. 161—169.
7. Hearle J. W. S., Grosberg P., Backer S. Structural Mechanics of Fibres, Yarns and Fabrics. N. Y.: Wiley-Interscience, 1969.
8. PlattM. M. Mechanics of elastic performance of textile materials: iii. some aspects of stress analysis of textile structures — continuous-filament yarns // Textile Research Journal. 1950. N 20. P. 1—15.
9. Treloar L. R. G. The stress-strain properties of multi-ply cords. Part I: Theory // J. Textile Inst. 1965. N 56. T477—T488.
10. Мусалимов В. М., Мокряк С. Я., Соханев Б. В., Шиянов В. Д. Определение упругих характеристик гибких кабелей на основе модели спирально-анизотропного тела // Механика композитных материалов. 1984. № 1. С. 136—141.
11. Мусалимов В. М. Механика деформируемого кабеля. СПб: СПбГУ ИТМО, 2005. 203 с.
12. Perechesova A. D., Soloveva G. A., Kalapyshina I. I. Hough transform for the calculation of twist angle of aramid torsion // Proc. WSCG 2015 — Poster Papers Proc. 2015. P. 72—75.
13. Сергушин П. А. Моделирование динамики магнитостатического вариометра // Научно-технический вестник СПбГУ ИТМО. 2007. Вып. 37. С. 333—338.
14. Sergushin P., Perechesova A., Petrishchev M. The torsion magnetic variometer with kevlar-hanger-based sensor // Proc. of ISMTII-2009, 29 June - 2 July 2009, St. Petersburg. Vol. 4. P. 411—414.
15. Мусалимов В. М., Заморуев Г. Б., Перечесова А. Д. Расчет физико-механических характеристик винтовых элементов спирально-анизотропных стержней // Изв. вузов. Приборостроение. 2012. Т. 55, № 6. С. 24—30.
16. Perechesova A. Calculation of elastic constants of the torsion bars with helical anisotropy using the methods of optimization theory // Proc. of the 23rd Intern. Congress of Theoretical and Applied Mechanics. Beijing, China, 2012. SM04-050.
17. Тимошенко С. П., Гудьер Д. Теория упругости. М.: Наука, 1975. 576 с.
Сведения об авторах
Виктор Михайлович Мусалимов — д-р техн. наук, профессор; Университет ИТМО; кафедра мехатрони-
ки; E-mail: musvm@ya.ru
Юрий Сергеевич Монахов — канд. техн. наук; Университет ИТМО; кафедра мехатроники;
E-mail: monahovu@gmail.com Алексей Юрьевич Кутьин — канд. техн. наук, доцент; Производственная компания „Индустрия
цвета"; программист; E-mail: kutin_alex@mail.ru Галина Александровна Соловьева — аспирант; Университет ИТМО; кафедра мехатроники;
E-mail: ymatem@rambler.ru
Рекомендована кафедрой Поступила в редакцию
мехатроники Университета ИТМО 05.04.16 г.
Ссылка для цитирования: Мусалимов В. М., Монахов Ю. С., Кутьин А. Ю., Соловьева Г. А. Моделирование процесса наматывания нитей на жесткий цилиндр // Изв. вузов. Приборостроение. 2016. Т. 59, № 8. С. 657—663.
MODELING OF THE PROCESS OF WINDING THREADS ON A RIGID CYLINDER
V. M. Musalimov1, Yu. S. Monakhov1, A. Yu. Kutin2, G. A. Solovieva1
1ITMO University, 197101, St. Petersburg, Russia E-mail: musvm@ya.ru 2Printing House "Industry of Color", 196084, St. Petersburg, Russia
Principles of the theory of deformation of helically-anisotropic body as applied to determination of parameters characterizing the stress strain state of the body are presented. Examples of the theory application to cable constructions and other helically-anisotropic bodies are presented. A special attention is paid to modeling the process of layered structure formation during winding. The model accounts for relationship between the elastic modulus and Poisson's ratio of helically-anisotropic materials, physical, mechanical, and kinetic parameters of the yarn and winding system. The presented analysis of winding body state allows to optimize the process of winding yarn on a rigid base.
Keywords: helically-anisotropic rod, composite material, circumferential winding, yarn tension, elastic constants
Data on authors
Victor M. Musalimov — Dr. Sci., Professor; ITMO University, Department of Mechatronics;
E-mail: musvm@ya.ru Yury S. Monakhov — PhD; ITMO University, Department of Mechatronics;
E-mail: monahovu@gmail.com Alexey Yu. Kutin — PhD, Associate Professor; Printing House "Industry of Color"; Programmer;
E-mail: kutin_alex@mail.ru Galina A. Solovieva — Post-Graduate Student; ITMO University; Department of Mechatronics;
E-mail: ymatem@rambler.ru
For citation: Musalimov V. M., Monakhov Yu. S., Kutin A. Yu., Solovieva G. A. Modeling of the process of winding threads on a rigid cylinder // Izv. vuzov. Priborostroenie. 2016. Vol. 59, N 8. P. 657—663 (in Russian).
DOI: 10.17586/0021-3454-2016-59-8-657-663
