Независимое управление напряжениями во вращающемся диске

Тип работы:
Реферат
Предмет:
Механика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 539. 37
НЕЗАВИСИМОЕ УПРАВЛЕНИЕ НАПРЯЖЕНИЯМИ ВО ВРАЩАЮЩЕМСЯ ДИСКЕ Туктамышев В. С., Лохов В. А.
ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет», Пермь, Россия
(614 990, г. Пермь, Комсомольский проспект, 29), e-mail: helpinvader@list. ru_
В данной работе алгоритм независимого управления напряжениями в системах с собственными (неупругими) деформациями реализован в рамках задачи понижения эксплуатационных нагрузок во вращающемся диске с неравномерным распределением поля температур. При этом под независимым управлением напряжениями подразумевается создание в диске заданных напряжений без изменения его полных деформаций. Соответствующее целевое напряжённое состояние определяется из отдельной оптимизационной задачи. Поиск собственных деформаций, обеспечивающих достижение решения поставленной задачи, осуществляется с помощью ранее разработанной методики независимого управления напряжениями, которая, в свою очередь, основывается на теореме о декомпозиции собственной деформации. В качестве собственных деформаций в представленной работе выбраны деформации фазовых переходов, возникающих в материалах с эффектом памяти формы, для диска, нагруженного поверхностным давлением. В предыдущих исследованиях показано, что распределение таких деформаций вычисляется с помощью простых соотношений. Однако это распределение оказывается свободным от напряжений, что не соответствует условию задачи. Одним из решений данной проблемы является создание составного диска, в каждой из частей которого может быть реализован эффект памяти формы. В работе показано, что увеличение количества составных элементов в диске лучше приближает его напряжённое состояние к требуемому при условии отсутствия полных деформаций.
Ключевые слова: собственные деформации, независимое управление напряжениями, эффект памяти формы, вращающийся диск.
INDEPENDENT STRESS CONTROL IN THE ROTATING DISC Tuktamyshev V.S., Lokhov V.A.
Perm National Research Polytechnic University, Perm, Russia (614 990, Perm, Komsomolsky st, 29), email: helpinvader@list. ru_
In this paper, an algorithm independent stress control in systems with imposed eigenstrain is implemented within the frame of the problem lowering operational loads in a rotating disk with a non-uniform distribution of temperature field. In this case, under the independent stress control means the creation of a disk desired stress field without changing the total strains. Corresponding target stress state is determined from the individual optimization problem. Search the eigenstrain that will promote the solution of the problem, by using the previously developed method of independent stress control, which in turn is based on a theorem of decomposition of the eigenstrain. As the eigenstrain in the present study selected strain of phase transitions that occur in materials with shape memory for the disc loaded with a surface strength. Previous studies have shown that the distribution of such deformations is calculated using the simple relations. However, this distribution is free of stress that does not meet the conditions of the problem. One solution to this problem is to provide a composite disc, in which each of the parts may be implemented shape memory effect. It is shown that an increase in the number of elements in the disk, it is better approximates the desired stress, in the absence of total strains. Keywords: eigenstrain, independent stress control, shape memory effect, rotating disk.
Введение
Возможность эффективной эксплуатации отдельных элементов современных технических систем требует разработки новых алгоритмов управления напряжённо-деформированным состоянием. Одной из актуальных проблем в данном направлении является создание методики независимого управления напряжениями в системах с неупругими деформациями. При этом под независимым управлением напряжениями
подразумевается создание в системе заданного поля напряжений с сохранением её деформированного состояния. Достижение цели в такой задаче осуществляется за счёт возможности создания в рассматриваемом теле неупругих деформаций, в качестве которых могут выступать деформации фазовых переходов, температурные деформации, пьезоэлектрические, ростовые и т. д.
В работе [5] представлен алгоритм независимого управления напряжениями, основанный на использовании теоремы о декомпозиции собственной деформации [7], которая утверждает, что любую собственную деформацию, существующую в теле, можно разложить на две составляющие: свободную от напряжений и свободную от полных деформаций. Теорема о декомпозиции, в свою очередь, является продолжением исследований о разложении гильбертовых пространств в задачах механики [3, 4, 9].
В представленной работе приведена демонстрация методики независимого управления напряжениями, где в качестве решения модельной задачи, приближённой к реальной технической проблеме, рассматривается вопрос понижения уровня механических напряжений во вращающемся диске с неравномерным по радиусу распределением температурного поля. Подобная задача может быть актуальной, например, при проектировании диска авиационной турбины.
Постановка задачи
Однородный диск постоянной толщины с радиусом Л1, выполненный из изотропного материала с эффектом памяти формы, имеет отверстие радиуса Л0. В процессе эксплуатации диск вращается с постоянной угловой скоростью ю, при воздействии стационарного неоднородного температурного поля, вызванного перепадом температур Т0 — Т1 между
внутренним и внешним контурами диска (рис. 1).
СО
В целях повышения рабочих характеристик при помощи собственных деформаций требуется создать такие преднапряжения в диске, которые, во-первых, позволили бы максимально снизить уровень эксплуатационных термоупругих напряжений и, во-вторых, не оказывали бы влияния на первоначальные размеры диска.
Рис. 1. Вращающийся диск под воздействием неоднородного температурного поля.
Рассматриваемую задачу можно разбить на
три части:
1) задача определения «рабочих» напряжений ар в диске-
2) поиск оптимального поля собственных напряжений р0, необходимого для требуемого
снижения напряжений & lt-7 —
3) создание в диске собственных деформаций, свободных от полных деформаций, которые, с учётом технологических ограничений, соответствуют напряжениям р0. 1. Решение первой части производится в рамках теории термоупругости. В полярных координатах при условии осесимметричного плоско-напряжённого состояния, распределения компонент тензора напряжений будет иметь вид [10]:

_ю, _Т 7 + & lt- ,
гг гг & gt-
& lt-Р _ _ю + _Т
(1)
где
а"_ (3 + П& gt-рЮ (г2 — Л2 — г2),
Бг2
& lt-
рю 8г2
(2)
(Л02(Л12 + г2)(3 +п) + г2(Л12(3+п) — г2(1 + 3п))),
вклад напряжений, вызванных центробежными силами инерции.
Г& gt-2/г>-2
& lt-1 _аЕ (То — Т1)
1п (Л0/ г). + (Ло2 — г2)'-
7ее _ (То — Т1)
1п (Л0/ Л) г2 (- Л02) г 1п (г) $ (Ло2 + г2)
(3)
1
1п (Ло/ /?) г2 (Л2 — Но2) 1п (Ло/ Ф
температурные напряжения (рис. 2а).
2. Необходимые собственные напряжения можно найти, исходя из следующих соображений.
С учётом собственных напряжений, суммарные радиальные и окружные напряжения в диске будут складываться из составляющих от центробежных сил при неоднородном нагреве и от собственных напряжений:
& lt- _ & lt- + р ,
гг гг г гг?
(4)
& lt-ее _ & lt-ее + рее ,
где ргг, рее — собственные напряжения. Эти напряжения могут быть представлены с помощью следующей функции напряжений Эйри ф (г):
рг (г) _1 ф (г), рее (г) _ ф
г аг
(5)
Вид функции ф (г) подбирается в соответствии с граничными условиями для ргг, рее:
Ф (г) _ г (г — Л0)(г — ВД (г), С (г) е С1 [ ?0, Л ], (6)
где ?(г) — произвольная гладкая функция с непрерывной производной на отрезке [Ц0, ].
Представим функцию Z® в виде следующего разложения:
Z® = К + + ^ (7)
r r
Тогда соотношения (5) перепишутся:
Prr ® = (r — Ro)(r — R^K + ^ + |,
k k2
2
rr
pee ® = (3r2 — 2(R0+R1)r+Ro R1) f k0++ ^ J + r (r — R0)(r — R1) f — -2
(8)
(9)
Перепишем напряжения (4) в виде следующих зависимостей:
Srr (Л ?0, К k2) = sFr (r) + Prr (Л ^ ^ k2) ,
see (Л ^ ^ k2) = sie ® + Pee (r, ^ ^ k2)•
Таким образом, для получения окончательной формы компонент prr, pee и srr, see требуется определить коэффициенты k0, k1, k2.
Как уже упоминалось выше, собственные напряжения в диске должны максимально понижать уровень эксплуатационных напряжений, т. е. необходимо решить следующую задачу минимизации:
f (^ К k2) = maxP (°rr (r, ^ kl, k2), see (r, ^ kl, k2)) ®inf (kl, k2, кз e D), (10)
R? r? Rj
Решение (10) позволит определить коэффициенты k0, k1, k2, соответствующие оптимальному полю собственных напряжений p0, вычисляемых по формуле (8). 3. После определения собственных напряжений p0 необходимо перейти к алгоритму решения задачи о независимом управлении напряжениями. Данная методика предполагает учёт природы собственной деформации создаваемой в диске, а также технологических ограничений, имеющих место в процессе реализации этих деформаций.
Реально существующая собственная деформация должна быть записана в форме выражения:
e * = ta* X k. (11)
k=1
Определение элементов Xk (k = 1. ^) будет производиться позже.
Согласно теореме о декомпозиции, необходимая для создания напряжений p0 собственная деформация, свободная от полных деформаций, имеет вид [7]
eS=-C& quot-1. p0. (12)
Тогда задача независимого управления напряжениями сводится к поиску минимума следующего функционала:
Il | а 2
Y = s * - s j ® inf (a), при условии e * i Hs.
0. OS 0. 1
0. 12 0. 14 016 0. 18 0.2 0. 22 0. 24
Решение
Как уже отмечалось ранее, диск должен быть выполнен из материала с эффектом памяти формы. Таким эффектом обладает, например, сплав никелида титана (Л№, нитинол), для одной из разновидностей которого в работе [6] (Т — 49. 16%, № - 49. 75%) приведены следующие физико-механические свойства:
• плотность материала у = 4450 кг/м3,
• модуль упругости Е = 75 ГПа (Т = 20 — 200°С),
• коэффициент Пуассона V = 0.3 ,
• коэффициент температурного расширения, а = 11×10−6 1/ °С,
• предел текучести, а = 450 — 800 МПа
(Т = 0 — 200°С),
• температурный интервал прямого мартенситного перехода: Ы5 = -92° С (в начале),
Ы{ =-106° С (в конце),
• температурный интервал обратного мартенситного перехода: Д = -44° С (в начале), А{ =-26° С (в конце).
Для оставшихся параметров, входящих в соотношения (2, 3, 8), приняты следующие значения [8]:
• радиус диска = 240 мм,
• радиус отверстия Л = 80 мм, • температура на внешней поверхности
200 100 о
-100 -200
-300
р0, МПа
_______
_1_ б) ---
*. г / р0г -- ~0

Г м Рее
0. 08 0.1 0. 12 0. 14 0. 16 0. 16 0.2 0. 22 0. 24
200
150
100
ьо


Grr -Оее

в) ---

r, м
0. 08 0.1 0. 12 0. 14 0. 16 0. 18 0.2 0. 22 0. 24
Рис. 2. Распределение эксплуатационных (a), остаточных (б) и суммарных (в) напряжений.
диска Т = 200 °C,
температура на внутренней поверхности диска Т0 = 100 °C, частота вращения диска п = 10 000 об/мин.
Численное решение (10) методом Нелдера-Мида определяет оптимальные значения искомых коэффициентов:
k0 = 1. 06×108 Па/м2, k1 = 3. 74×108 Па/м, k2 = 0. 76×108 Па.
Соответствующее этим значениям распределение оптимальных собственных напряжений, а также поле суммарных напряжений в диске приведены на рис. 2. При этом поле деформаций фазовых переходов, соответствующее собственным напряжениям р0, вычисляется по формуле (12).
В современной литературе существует достаточно большое количество математических моделей, описывающих эффект памяти формы. Среди них одной из наиболее простых и в то же время хорошо согласующихся с экспериментальными данными является модель, описанная Мовчаном в работе [1].
Пусть деформации, вызванные эффектом памяти формы, реализуются за счёт поверхностных давлений рх и p2, приложенных к внутреннему и внешнему контурам диска
соответственно. Такие деформации исследованы в работе [2] и определяются следующими выражениями:
С = - (-1) s", eq0= ^ (-1) 500, (14)
a0 a0
где srr и s00 — компоненты девиатора тензора напряжений
1 (r2 — 3Г22) r2 (3rf — r2)
^ = Р 3г2 (г2 — г!2) + Р23г2 (г2 — г!2)'-
У 7 У 7 (15)
= г2 (г2 + 3Г2) Г22 (3/12 + г2)
500 = Р Зг2 (г2 — Г2) + Р23г2 (г2 — г2) •
Однако деформации (14) оказываются свободными от напряжений, что, вообще говоря, нежелательно, согласно условию задачи.
Одним из решений данной проблемы является создание составного диска, в каждой из частей которого за счёт поверхностного давления реализуется эффект памяти формы. Под воздействием этого нагружения элементы охлаждаются до температуры более низкой, чем конечная температура прямого мартенситного перехода [1]. При этом в составных частях накапливаются собственные (фазовые) деформации. После снятия нагрузки упругие деформации становятся равными нулю, и элементы вставляются друг в друга без натяга. При нагревании составного диска до температуры более высокой, чем конечная температура обратного мартенситного перехода, существующие собственные деформации исчезают, следовательно, составные части будут стремиться к восстановлению первоначальных размеров, что в свою очередь приведёт к появлению собственных напряжений.
С учётом формул (14) и (15) деформации фазовых переходов для составных частей диска можно выразить следующим образом:
* г2 (г2 — 3г22) * г22 (3г2 — г2) _
г1гг р11−2/2 2 + Р22/2 2, 1 0,
3 г (гп — г1) 3 г (г12 — г 1)
'- '- (16)
* _. г2 (г2+3г22). г2 (3г2+г2) _ 1 е'-ее~ Р113г/ (г2 —) + Р123г/ (г2 — г2), ^ 1& quot-0,
где 0 — количество составных элементов, г — радиальная координата составного элемента с номером 1, г 1, г2 — внутренний и внешний радиусы этого элемента (г 1? г? г2),
Р*1 _ Р 1 — (^ - 1), Р*2 _ Р, 2 ^ (^ - 1), 1 _ 1. 0,
ао ао
Р11 — давление на внутреннюю поверхность элемента с номером 1, Р12 — давление на
внешнюю поверхность того же элемента (1 _ 1. 0).
В системе существует 0 подобластей, в которых могут быть созданы независимые друг от друга собственные деформации, тогда согласно формуле (11)
е* _? е*, (17)
11
где е* - собственные деформации в соответствующей подобласти. Тензор е* (1 _ 1. 0)
можно представить в виде вектора-столбца, поскольку в данной задаче ненулевыми являются лишь две его компоненты (индекс «Т» означает транспонирование):
е* _{е*гг, е*ее}Г, 1 _ 1. 0. (18)
Из соотношений (16) видно, что величинами собственных деформаций можно управлять посредством изменения давлений Р11 и Р12 1 _ 1. 0, которые входят в параметры
Рл и Р*2 соответственно. Таким образом, в каждой подобласти существует два независимых
параметра управления. Введём следующее обозначение:
Р*1 _ Р*, Р*2 _ Р*, Р*1 _ Рз*, Р*2 _ Р^ к
Элементы Р* (к _ 1. 20) с нечётным номером к соответствуют внутреннему давлению,
элементы с чётным номером — внешнему давлению.
В связи с этим выражение (17) можно переписать следующим образом
20
е* _? р*Х к, (4. 4о)
к _1
где согласно (16) и (18)

(Г -Зг22) (Г+3гД)
3г2 (? — г2!)'- 3г2 (г22 — г2:)
и Хк
о, МПа
ЖТГ… __А'- Г*
** м*
А ж… У- & quot-"-

/ ?
* i _____ /
? /
0. 08 0.1 0. 12 0. 14 0. 16 0. 1 В 0.2 0. 22 0. 24 Г, м
г22 (зг2 -г2) Г22 (3^+г2 зг2 (г22 -г21), зг/ (г2: -г22)
нечётных номеров
Рис. 3. Распределение собственных напряжений.
для чётных и соответственно.
Распределение требуемых собственных напряжений и получаемых в результате управления зависит от количества 0 составных элементов диска. При этом для любого числа 0 полные деформации остаются равными нулю, так как поиск оптимального решения производится среди собственных деформаций, свободных от полных деформаций.
На рис. 3. представлены распределения собственных напряжений для диска, состоящего из двадцати равных по толщине составных элементов. Собственные напряжения, получаемые в результате вычислений (о — радиальные напряжения, А — окружные), практически полностью совпадают с распределением р0 (рис. 2б). Однако использование большого количества составных частей несёт в себе существенный недостаток, заключающийся в сложности изготовления таких дисков.
Заключение
Теорема о декомпозиции собственной деформации открывает возможность независимого управления напряжениями, что является актуальным для некоторых прикладных задач. В данной работе приведён пример использования алгоритма независимого управления напряжениями в задаче понижения эксплуатационных нагрузок вращающегося диска, находящегося под действием неоднородного поля температур. В рамках представленной проблемы показаны возможности достижения требуемого поля напряжений при условии отсутствия полных деформаций.
Работа выполнена при финансовой поддержке РФФИ (проект № 12−01−31 404_мол_а).
Список литературы
к
1. Мовчан А. А. Микромеханические определяющие уравнения для сплавов с памятью формы // Проблемы машиностроения и надежности машин. — 1994. — № 6. — С. 47−53.
2. Лохов В. А., Туктамышев В. С. Исследование условий отсутствия механических напряжений в системах с собственными деформациями // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. — 2013. — № 2. — С. 198−207.
3. Поздеев А. А., Няшин Ю. И., Трусов П. В. Остаточные напряжения: теория и приложения. — М.: Наука, 1982. — 112 с.
4. Соколов А. Г, Стружанов В. В. Об одной задаче оптимизации напряженного состояния в упругом теле // Прикл. мат. и механика. — 2001. — Т. 65. — № 2. — С. 317−322.
5. Туктамышев В. С., Лохов В. А., Няшин Ю. И. Независимое управление напряжениями в неоднородных системах с собственными деформациями // Математические методы и физико-механические поля. — 2012. — № 4. — С. 131−142.
6. Cubela D. The research of technological parameters influence on the process of nitinol fabrication and plastic deformation // Metalurgija. — 2006. — Vol. 45. — № 1. — P. 61−66.
7. Nyashin Y, Lokhov V, Ziegler F. Decomposition method in linear elastic problems with eigenstrain // ZAMM — Z. Angew. Math. Mech. — 2005. — Vol. 85. — № 8. — P. 557−570.
8. Nyashin Y., Shishlyaev V. Analytic creep durability of rotating uniform disks // International Journal of Rotating Machinery. — 1998. — Vol. 4. — № 4. — P. 249−256.
9. Rafalski P. Orthogonal projection method. II. Thermoelastic problem // Bull. Acad. Polon. Sci. Ser. Sci. Tech. — 1969. — Vol. 17. — № 2. — P. 69−74.
10. Ziegler F. Mechanics of Solids and Fluids, 2nd ed. — New York, Springer, 1995. — 911 p. Рецензенты:
Акулич Ю. В., д.ф. -м.н., профессор кафедры теоретической механики, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Пермский национальный исследовательский политехнический университет», г. Пермь.
Кротов Л. Н., д.ф. -м.н., доцент, зав. кафедрой прикладной физики, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Пермский национальный исследовательский политехнический университет», г. Пермь.

ПоказатьСвернуть
Заполнить форму текущей работой