ISSN 0321-3005 ИЗВЕСТИЯ ВУЗОВ. СЕВЕРО-КАВКАЗСКИЙ РЕГИОН. ЕСТЕСТВЕННЫЕ НАУКИ. 2015. № 2
УДК 539.3
КОНТАКТНЫЕ ЗАДАЧИ ДЛЯ СОСТАВНОГО ПОЛУПРОСТРАНСТВА* © 2015 г. Д.А. Пожарский, В.И. Полтинников
Пожарский Дмитрий Александрович - доктор физико-математических наук, профессор, заведующий кафедрой прикладной математики, Донской государственный технический университет, пл. Гагарина, 1, г. Ростов н/Д, 344000, e-mail: pozharda@rambler.ru
Полтинников Виктор Иванович - доцент, кафедра прикладной математики, Донской государственный технический университет, пл. Гагарина, 1, г. Ростов н/Д, 344000.
Pozharskii Dmitrii Aleksandrovich - Doctor of Physical and Mathematical Science, Professor, Head of Applied Mathematics Department, Don State Technical University, Gagarin Sq., 1, Rostov-on-Don, 344000, Russia, e-mail: pozhar-da@rambler.ru
Poltinnikov Viktor Ivanovich - Associate Professor, Applied Mathematics Department, Don State Technical University, Gagarin Sq., 1, Rostov-on-Don, 344000, Russia.
Получены интегральные уравнения (ИУ) трехмерных контактных задач для упругого полупространства, составленного из двух клиновидных слоев, соединенных скользящей заделкой. Клиновидный слой, примыкающий к слою, в который вдавливается штамп, несжимаем (резино-металлическое сочленение). Внешняя грань этого слоя свободна от напряжений либо подчинена условиям скользящей заделки. Для решения вспомогательных краевых задач о действии заданной нормальной нагрузки применен метод комплексных интегральных преобразований Фурье и Кон-торовича - Лебедева, позволивший свести их к системам ИУ Фредгольма второго рода, решения которых затем вошли в ядра ИУ контактных задач. Для решения контактных задач использован метод Галанова.
Ключевые слова: контактная задача, составное полупространство, клин.
The integral equations are derived for three-dimensional contact problems for an elastic half-space consisting of two wedge-shaped layers with sliding support in between. The wedge-shaped layer joined to the layer in contact with a punch is incompressible (rubber-metal conjunction). The outer face of this layer is supposed to be stress-free or subjected to sliding support. To solve auxiliary boundary-value problems for a given normal load the method of Fourier and Kontorovich - Le-bedev complex integral transformations is used. These auxiliary problems are finally reduced to systems of Fredholm integral equations of the second kind, solutions of which go into the kernels of the integral equation of the contact problems. For solving the contact problems the Galanov's method is used.
Keywords: contact problem, composite half-space, wedge.
Ранее изучались аналогичные плоские задачи [1], рассматривался случай, когда клиновидный слой под штампом несжимаем [2].
Рассмотрим, используя цилиндрические координаты, упругое полупространство, составленное из клиновидных слоев Qi={re [0,<ю); фе[а-л,-а]; ze(-<»,<»)} и Q2={re [0,<ю); фе[-а,а]; ze(-<»,<»)}, которые имеют углы раствора л-2а и 2а соответственно и соединены скользящей заделкой по грани ф=-а. Материал клина Q1 несжимаем (коэффициент Пуассона v1=0,5, модуль сдвига G1). Материал клина Q2 имеет коэффициент Пуассона v и модуль сдвига G. В грань ф=а под действием силы P без перекоса вдавливается жесткий штамп. Осадка штампа равна 8, а форма основания описывается функцией frz). Внешняя грань клина Q1 ф=а-тс свободна от напряжений (задача А) или находится в условиях скользящей заделки (задача Б). При известных величинах а, v, G, G1, 8 и заданной функции fy,z) требуется определить контактное давление CT^r^z^-q^z), а также область контакта Q. Затем можно определить связь между величинами P и 8.
Для вывода интегральных уравнений (ИУ) контактных задач А и Б используются метод комплексных интегральных преобразований Фурье и Конторовича-Лебедева и техника сведения краевой задачи теории упругости для трехмерного клина к системе функциональных уравнений со сдвигом аргумента [3]. Функционально-сдвиговые уравнения удается свести к системе двух ИУ Фредгольма второго рода относительно вспомогательных функций Фп(т) (п = 1, 2; Х>0, 0<т«»)
Фп (т) + 6(1-у)0(т)[ВДФ1(т) + (1)
ад
+ W2(т)Ф2(т)] - (1 - 2у)|Ь„ (т, и)Фп №и = 0
лт
= + 6(1 - у)б(т)[ Щ (т) + W2 (т)]еЬ — Кт (Хх) +
+ (1 - 2v) J Ln (т, u)ch
0
nu 2
Kiu (Xx)du,
-г- .. ~чПТ - nu , Ln (т, u) = 2ch — sh — Wn (u) x
3
x|
shnt gn (t )dt
0 (chnt + chnT)(chnt + chnu)
e =
Gl
G
* Работа выполнена по госзаданию Минобрнауки РФ (проект 2941).
зо
ISSN 0321-3005 ИЗВЕСТИЯ ВУЗОВ. СЕВЕРО-КАВКАЗСКИЙ РЕГИОН. ЕСТЕСТВЕННЫЕ НАУКИ. 2015. № 2
W (г) = +
gn (г) =
сЪ2ат + cos 2а sh2ar + г sin 2а
sin 2а
cthar
0(г) = 0(г) =
ch2ar + cos 4а [Шаг
2(sh2(rc - 2а)г - г2 sin2 2а)
sh2(rc - 2а)г -г sin 4а sh2(rc - 2а)г -г sin 4а
(задача А),
(задача Б).
еИ2(л - 2а)х - cos 4а
Здесь и ниже при n = 1 берется верхний знак (функция), а при n = 2 — нижний. Правые части уравнений (1) экспоненциально убывают при х ^ да. Для численного решения системы (1) применим метод механических квадратур.
Контактные задачи затем сводятся к ИУ относительно q(r,z):
Л q(x, y)K(r, z, x, y)dxdy = ^^ [8- f (r, z)], (2)
1 -v
(r, z) еД К (r, z, x, y) =
14 f
= - + — IIj
R Л 00 [
(3)
^(г) - W2(x)
2
- rtfr^
Кг (Xx) +
^СрФ^г) - W^Q^) 2А(лг/2)
Кг (Xr) cos X( z - y)dX^dx.
Здесь функции Ф„(т) (и=1, 2) удовлетворяют системе (1). Главный член R 1 в ядре (3) соответствует случаю однородного упругого полупространства [3]. При е=0 ядро (3) совпадает с известным [3] для однородного клина с одной свободной от напряжений гранью, символом ядра которого служит функция
^(г) - W2(^ 2
sh4aг + г sin 4а А4аг - г2(1 - cos 4а) -1
В частном случае, когда у=0,5, система (1) имеет точное решение
ф (т) = +еб(т)[ад+сЬ ™
" 2 + ее(т)[^1(т) - ВД] 2 гт При у=0,5, е=1, 2а=л/2 (две одинаковые несжимаемые четверти пространства, соединенные скользящей заделкой) символ ядра ИУ (2)
^!(г) - W2(г) ß(г)[Wl(г) + W2(г)]2
2
4 + 20(г)[^1(г) - W2 (г)]
(4)
sh2^-2г2сh 2(лг/2)
2shлг(sh2(лг/2) -г2) sh^
(задача А)
или =
2sh2(лг/2) -г2
(задача Б).
Функции в правой части формулы (4) совпадают с известными символами ядер интегральных уравнений плоской и пространственной контактных задач, рассмотренных ранее. Задача А в случае (4) соответствует задачам А [1] и В [2]. Задача Б в случае (4) соответствует задачам Б [1] и А [2]. При этом в соответствующих формулах [1] следует взять e=1, а=Р=тс/2, а в [2, (1.5)—(1.9)] следует положить v=0,5, s=1, а=Р=тс/2, ц=0.
Для численного решения ИУ (2) применим метод Галанова нелинейных граничных ИУ, позволяющий одновременно определить область контакта и контактные давления [3]. Предположим, что область контакта содержится в прямоугольнике S, который покроем равномерной сеткой. При расчете значений ядра (3) в узлах сетки его особенности сглаживались путем регуляризации [3]. Как показывают расчеты, в случае задачи А значения контактного давления и вдавливающей силы несколько меньше, чем в задаче Б (при тех же значениях задаваемых параметров).
Литература
1. Пожарский Д.А. Смешанные задачи теории упругости для составного плоского клина // Изв. вузов. Сев.-Кавк. регион. Естеств. науки. 2008. № 5. С. 36-38.
2. Александров В.М., Пожарский Д.А. Пространственные контактные задачи с трением для составного упругого клина // Прикладная математика и механика. 2010. Т. 74, вып. 6. С. 969-977.
3. Александров В.М., Пожарский Д.А. Неклассические пространственные задачи механики контактных взаимодействий упругих тел. М., 1998. 288 с.
References
1. Pozharskii D.A. Smeshannye zadachi teorii upru-gosti dlya sostavnogo ploskogo klina [Mixed elasticity problem for a composite plane wedge]. Izvestiya vuzov. Severo-Kavkazskii region. Estestvennye nauki, 2008, no 5, pp. 36-38.
2. Aleksandrov V.M., Pozharskii D.A. Prostranstven-nye kontaktnye zadachi s treniem dlya sostavnogo upru-gogo klina [The spatial contact problem with friction for the composite elastic wedge]. Prikladnaya matematika i mekhanika, 2010, vol. 74, no 6, pp. 969-977.
3. Aleksandrov V.M., Pozharskii D.A. Neklassi-cheskie prostranstvennye zadachi mekhaniki kontaktnykh vzaimodeistvii uprugikh tel [Non-classical spatial problems of mechanics of contact interactions of elastic bodies]. Moscow, 1998, 288 p.
Поступила в редакцию
2 апреля 2015 г.
Д