2460
Фундаментальные проблемы теоретической и прикладной механики Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2460-2462
УДК 539.3
НЕКОТОРЫЕ МЕТОДЫ РЕШЕНИЯ ПЛОСКИХ ЗАДАЧ ТЕОРИИ УПРУГОСТИ НЕОДНОРОДНЫХ ТЕЛ
© 2011 г. Н.А. Роганова1, Г.З. Шарафутдинов1
'Московский государственный индустриальный университет 2НИИ механики Московского госуниверситета им. М.В. Ломоносова
Поступила в редакцию 24.08.2011
Рассматривается класс задач, связанный с деформированием трубы, неоднородность материала которой зависит от условий его синтеза или действия различных причин. Также разработана математическая модель для анализа напряженно-деформированного состояния в двумерных задачах теории упругости неоднородных тел. Получены новые формы основных уравнений, предложен метод последовательных приближений при решении этого класса задач. Представлены общие выражения для напряжений в виде рядов, доказана их сходимость. При помощи функций комплексного переменного получены аналитические выражения для последовательных приближений.
Ключевые слова: теория упругости, плоская задача, неоднородные материалы.
Неоднородность упругих характеристик деформируемого изотропного материала зададим в виде переменного модуля сдвига |1, являющегося непрерывной функцией координат, и постоянным или кусочно-постоянным коэффициентом Пуассона V.
Рассматривается класс задач, связанный с деформированием толстостенной цилиндрической трубы с внутренним радиусом а и внешним радиусом Ь, неоднородность материала которой зависит от условий его синтеза или действия различных причин: агрессивной среды, радиации и т.п. Пусть на внутреннюю и внешнюю поверхности трубы действуют равномерные давления р и д соответственно. Помимо этого полагаем, что в продольном направлении действует сила Для определенности будем считать, что внутри трубы находится агрессивная среда, проникновение ко -торой в тело трубы приводит к неоднородному изменению упругих характеристик материала. В силу осевой симметрии характеристики материала трубы зависят лишь от радиальной координаты; при этом радиальное перемещение и = и(г). Окружные перемещения положим равными нулю, осевую деформацию £- считаем постоянной. Ус -тановлено [1], что радиальное перемещение и = = и(г) удовлетворяет уравнению:
ц(г)
- +
|Дг) + йф)
йг2 1 - 2v ц(г) 1 - V г2
йи
йг
и = -
йг ф,(г) 1 - V йг
8,
Уравнение (1) решается методом прогонки при любом типе краевых условий, после чего нетрудно определить компоненты тензоров деформаций и напряжений. Рассматривается также альтернативный подход, в соответствии с которым при действии агрессивной среды внутри трубы вводится, согласно [2], понятие диффузионного фронта с границей раздела г = с. При этом предполагается, что внутренний слой трубы имеет переменный модуль сдвига ц(г) и V = 0.5, в то время как внешний неповрежденный слой ма -териала трубы сохраняет свои однородные свойства. Задача сводится к задаче теории упругости для неоднородного тела в области а < г < с и к задаче Ламе в области с < г < Ь. Таким образом, во внешней области имеем решение задачи Ламе, а во внутренней области [1]
<5г = 2шу + 2ш2[А(г) -В(г)],
<за = 2т + 2ш2 [А(г) + В(г)],
а2 = 2т1 + 2т2 А(г) + 3|о,( г )8 2,
А(г) = |
8 7 т2
и =—- + —,
2 г
В (г) =4»,
(1)
где постоянные т1?т2 определяются путем удовлетворения краевым условиям.
Разработана также математическая модель для исследования напряженно-деформированного состояния в двумерных задачах теории упругости неоднородных тел. Рассматривается плос-
г
г
г
кая задача неоднородного упругого тела в напряжениях. Граничные условия представим в виде агу/у |^ = Pj, где L — граница некоторой односвязной области S; pt — компоненты вектора внешнего напряжения. Неоднородность материала зададим в виде V = const, Ц = Ц0М(х j, x2), 0 < < M(x1, x2) < 1, (Xj, x2) e S, где Ц0 — размерный множитель. Вводятся безразмерные компоненты тензора напряжений Ту, связанные с размерными компонентами Су этого тензора соотношениями Ojj = 2ц0М(х1, х2)т^. Уравнения равновесия, выраженные относительно величин Ту, принимают вид
дтп дт
+ -
12
д ln M
д ln M
дх дх,
= -т
11'
1
дт12
+
2
дт 22
дх
-т
12
д ln M
дх2 д ln M
= -т
т12 д т22
дх1 дх2 дх1 дх2
(2)
(1 -V)
дх1
д
дх2
д
2
кдх1 дх22 j
д ln M т11—------+ т
(т11 + т22) = д ln M
дх1
12
дх2
д ln M д ln M
Ті о---------------------+ Т
Т(к+1) _
11
= £(-!)'
ln "M д2Ф
к -п+1
(k+1)
п =0
k
п!
дх22
22
ln nM г2Ф
к-п +1
т
(k+1) =
п=0
k
п!
дх12
12
= Z (-1)
n+1 ln nM д2ф к - п+1
п=0
п!
дх1дх2
= (1 - v)V2 Z (-1)
(1 - v)V4Фk+1 = ln п+1 M (п +1)!
k-1
V2Ф, +
к-п
+ 2-
д
2
дх1дх д
дх
д
12
2
2
1
2
п=0 ' к-1
,ln п+1 M д2Фк_
Х(-1)п
п=0 (п +1)! дх1дх2
к-1
Z (-1)
п=0 (п +1)!
ln п+1 M д2Фк -п
дх22
дх22
к-1
Z(-1)
п=0 (п + 1)!
lnп+1 M д2Фк-п
дх12
Доказана сходимость выражений для безразмерных компонент тензора напряжений [3]. По-
лученная система упрощена путем применения функций комплексного переменного [4]. Такой подход позволил получить аналитические формы первого, второго и последующих приближений.
Применим указанный подход к решению следующей задачи. Пусть в неоднородном упругом пространстве имеется круговая цилиндрическая полость бесконечной протяженности и пространств о подвергается одноосному растяжению в направлении, перпендикулярном продольной оси цилиндрического выреза, заданными на бесконечности напряжениями. Будем считать, что на бесконечности заданы напряжения оп = р, а поверхность кругового выреза радиуса Я свободна от внешних воздействий. Неоднородность материала зададим при помощи двухпараметрического семейства функций вида
При учете закона Гука и уравнений равновесия условие совместности приводится к виду
ц( х1, х2) = —exp m
R
2 V
х2 + х2
V X1 + x2 j
42 Г + т22 д . (3)
\^1 д%2
Для решения системы уравнений (2), (3) применим метод последовательных приближений. Общее решение системы запишется в виде
х2 + х22 = г2 > Я2, где Ц0 - размерный модуль сдвига, параметры т и д определяют вид неоднородности материала. Коэффициент Пуассона V считаем постоянным. При решении этой задачи в качестве нулевого приближения выбирается тривиальное решение. Заметим, что при определении комплексных потенциалов удобно применить интеграл типа Коши.
В качестве иллюстрации на рис. 1 и 2 приведены зависимости компонент с^Р и с(3> от радиуса для а = п/2 при р = 1, |а0 = 1, Я = 2, V = 0.33.
(3)
О,
0.4
0.3
0.2
0.1
0
/ / , і /
1 1 / '"■ц
і/ / 1 ‘
ІУ
2
4
Рис. 1
д
3
5
r
Сплошная линия — решение для однородного тела; пунктирная — решение для д = 0.5, т = ^[ё, штрихпунктирная — для д = —0.5, т = 1.
Список литературы
1. Шарафутдинов Г.З. Некоторые осесимметричные задачи для упругой неоднородной толстостенной трубы // Вестн. МГУ Сер. 1. Математика, механика. 2008. №2. С. 34—39.
2. Локощенко А.М. Ползучесть и длительная прочность металлов в агрессивных средах. М.: Изд-во Моск. ун-та, 2000. 178 с.
3. Роганова Н.А., Шарафутдинов Г.З. Приближенный метод решения плоских задач теории упругости неоднородных тел // Машиностроение и инженерное образование. 2009. №3 (20). С. 63—71.
4. Роганова Н.А., Шарафутдинов Г.З. Применение функций комплексного переменного в задачах плоской деформации неоднородных тел // Известия МГИУ 2008. №1 (10). С. 75—84.
SOME METHODS OF ANALYZING PLANE ELASTICITY PROBLEMS OF INHOMOGENEOUS BODIES
N.A. Roganova, G.Z. Sharafutdinov
We consider a class of problems related to the deformation of a pipe of an inhomogeneous material, its inhomogeneity being caused by the conditions of its synthesis or by other effects. A mathematical model is developed for analyzing the stressed-strained state in two-dimensional elasticity problems of inhomogeneous bodies. New forms of basic equations are derived, and the method of successive approximations for solving this class of problems is suggested. General expressions for the stresses in the form of series are given, and their convergence is verified. Using complex variable functions, analytical expressions for the successive approximations are obtained.
Keywords: theory of elasticity, plane problem, nonhomogeneous materials.