CC BY
38
УДК 539.3
ТРЕХМЕРНАЯ КОНТАКТНАЯ ЗАДАЧА О ВЗАИМОДЕЙСТВИИ УПРУГОГО СЛОЯ С ДВУМЯ ШТАМПАМИ ПРИ УЧЕТЕ ТРЕНИЯ*
Д.А. ПОЖАРСКИЙ, А.А. МОЛЧАНОВ
(Донской государственный технический университет)
Изучены пространственные контактные задачи для упругого слоя конечной толщины, в одну грань которого симметрично вдавливаются два одинаковых жестких эллиптических штампа с учетом трения при разных типах граничных условий на другой грани. Задачи сведены к интегральным уравнениям относительно контактного давления, которые решены методом Галанова.
Ключевые слова: теория упругости, контактная задача, слой, трение, метод Галанова.
Введение. Исследованы трехмерные контактные задачи теории упругости при учете трения о взаимодействии слоя с двумя симметричными штампами, расположенными на одной его грани. Другая грань слоя находится в условиях жесткой или скользящей заделки. Штампы имеют форму эллиптических параболоидов, начинают удаляться друг от друга или сближаться. Области контакта неизвестны. Ранее аналогичные задачи с трением рассматривались для случая одного штампа на слое [1 - 4] и на полосе [5]. В работе [6] рассматривается пространственная задача о контакте с упругим слоем системы двух симметричных эллиптических штампов с плоской подошвой (асимптотический метод решения).
Постановка задачи. В декартовых координатах рассмотрим слой |х, у е(-да, да), г е[0, й]} толщиной h, нижняя грань которого z = 0 находится в жесткой или скользящей заделке (задачи А и Б соответственно). Упругий материал слоя имеет коэффициент Пуассона V и модуль сдвига G. Верхняя грань слоя z = h взаимодействует с двумя одинаковыми штампами. Формы основания штампов имеют вид эллиптических параболоидов и описываются функциями
д±(х,у) = ^О^ЯО+у2/^), Я > Я (1)
Между поверхностью слоя и штампов действуют силы кулоновского трения с коэффициентом трения ц. Штампы начинают достаточно медленно двигаться вдоль оси х так, что задачи симметричны относительно оси у. Силы трения направлены против движения. При ц > 0 штампы удаляются друг от друга, а при ц < 0 — начинают сближаться. К штампам приложены симметричные касательные силы Т, нормальные силы Р. Пусть осадка штампов равна 5, а перекос отсутствует. Симметричные по у области контакта неизвестны (область 0+ при х < 0).
При известных величинах G, V, ц, h, с, Я,, Я2 и 5 и заданной функции 1(г^ требуется определить контактные давления а^ху^) = -^х,у), (х,у)е^ = 0+и0-, а также сами области контакта 0±. Затем можно найти, например, величину Риз условия равновесия штампа
Р = Л q( х, у^хёу. (2)
о+
Из аналогичных интегральных условий можно найти величину Т, а также плечи приложения сил Ри Т.
Решение задачи. Предположим, что области контакта априори содержатся в прямоугольниках 5± = {|х±с| < а, у| < Ь}, Ь > а, с > а.
Для вывода интегральных уравнений (ИУ) контактных задач А, Б используется интегральное преобразование Фурье и закон Кулона. В результате придем к ИУ на двух участках контакта, которое после введения безразмерных обозначений
. х , у . й . с а 5 . Ь п Ь
х =—, у = —, А, = —, с = —, в0 =—, 5 =—, А =------, В =----, (3)
Ь Ь Ь Ь 0 Ь Ь 2Rl 2R2 к }
*
Работа поддержана грантом РФФИ 09-01-00004.
q'(x',y') = *\ , P1 , Є = --, Q' oQ, Q±oQ±, S±o S± (4)
q(xy) P і = P e = —
2лЄ , 2лЄЬ2 ’ 1 -v
можно записать в виде (штрихи далее опускаем):
—+є X-.4+1 т Г x-i, y-n
R R2 Я І Я Я
Ц q(i, П)
did' = f (x,y), (x, y) є Q, (5)
f (x, y) = 5- A(x ± c)2 - By2, (x,y) eQ±, (6)
R±= [(x ±i)2 + (y -n)2]1/2, е = ц(1 - 2v)/(2 - 2v), (7)
да _______ _______
T(t,т) = j[Lj(u) - 1]J0(u\/t2 + x2 )du +—, S I[L2(u) - 1]J^j(Wt2 + x2 )du. (8)
о vt2 + т2 о
Здесь Jn(u) — функции Бесселя. В ядре ИУ (4) выделена главная часть.
Для задачи А (жесткая заделка)
т , ч 2к sh2u - 4u _ „
A(u) =--------------------------------------------------------------------------2-2, к = 3-4v, (9)
2к ch2u + 4u +1 + к
2к ch2u - 4(1 - 2v)-1 u2 - 2к
L2 (u)— ~ ~ , (10)
2к ch2u + 4u +1 + к
а для задачи Б (скользящая заделка)
т ch2u -1 sh2u - 2(1 - 2v)-1 u
L,(u) =---, L,(u) =-----------------------------------------------------------------. (11)
sh2u + 2u sh2u + 2u
Введенный в (3) безразмерный параметр X характеризует относительную толщину упругого слоя с учетом симметрии задач q(-x,y) = q(x,y). Тогда уравнение (5) сводится к ИУ на одном участке контакта
Ц q(i n)
1 1 x-i x + i 1 Гx-i y-n
----+ — + є—^ + є—^ + -ТІ------5- ,--1 | +
R R R2 R2 Я І Я Я
+1ТГx+s y-n
(12)
didn = 8- A(x - c)2 - By2, (x,y) є Q_.
X у X X После замен:
х,= х - с, = с, 4 (х,, у) = д(х, у), О,-о-О-, £, -о- 5- (13)
ИУ (12) можно переписать в форме (звездочки далее опускаем):
Ц 4^, Л)К(І, Л, х, у) =£ (х, у), (х, у) є О,
(14)
1 1 x -i x + i + 2с
K (i, n, x, y) = — + — + є—— + є-------- ----+
R - R1 R- R12
+1ТГx-i y-nl + 1ТГx + ^ + 2c y-n
(15)
Я ^ Я Я ^ Я ^ X Я у
g(x,у) = 5- Ах2 -Ву2, R1 = [(х + ^ + 2с)2 + (у-^)2]Ш. (16)
При численном решении уравнения (14) применим метод нелинейных граничных ИУ типа Гаммерштейна, предложенный Галановым [7, 8]. Уравнение (14) дополним условиями неотрицательности контактного давления в области контакта, отсутствия контакта и обращения в нуль давления в дополнительной области 5\0, записав их все в виде системы
| К(Ы,М^(ЩсЫ = g (М), q(M) > 0, М еП,
5
| К(Ы,М^(> g (М), ^(М) = 0, М е 5 \ П,
(17)
я
где введены обозначения М = (х, у), N = (£,^).
Идея метода состоит в представлении искомого давления в форме
q = q(M) = q+ (М) + q - (М), (18)
где введены нелинейные операторы
q+ (М) = sup{q+ (М),0}, q- (М) = т%(М),0}. (19)
При учете (17) интегральное неравенство (16) будет удовлетворено в результате решения нелинейного операторного уравнения типа Гаммерштейна:
© р = 0 (М е О), © р = р~ + К,р+ - g, (20)
где р, = р, (М), р; = р; (М), g = g (М),
К, р + ={ К (N, М) р+ (N)dN. (21)
я
Можно доказать, что система (17) эквивалентна уравнению (20) [7]. Исследованы вопросы существования и единственности решения уравнения типа (20) [7]. Для численного решения уравнения (20) применим метод М. А. Красносельского, основанный на последовательных приближениях [7].
Прямоугольник 5 покроем сеткой из т узлов с учетом отсутствия симметрии по координате х. Ясно, что с > е0, е0 < 1. Для проверки точности расчета ядра (15) можно использовать интеграл [9]
Dn
Гехр(-Си)J (Du)du = , =-----, —, п = 0;1.
і 4сг+о2 [с+4сг+о2 г
(22)
Исследуя поведение выражений Ln (и) -1 (п = 0;1) при и для функций (9) - (11), положим в (22) С = 2. При выполнении условия
Х>^/4 (1 + с )2 +1/ (4 л/5) (23)
для расчета ядра (15) можно использовать квадратурную формулу Гаусса по 32 узлам. Численные эксперименты. В табл. 1, 2 для задач А и Б соответственно даны значения контактного давления с(х,0) и вдавливающей штамп силы Р. Расчеты, результаты которых приведены в таблицах, сделаны при V = 0,3, 5 = 0,004, А0 = 0,1, £ = 0,005, є0 = с = 0,15 и разных значениях X и ц.
Таблица 1
Значения давления и вдавливающей силы Р в задаче А
X ц С(х,0) х103 Р<103
х=-0,075 -0,0375 0 0,0375 0,075
1 0 1,77 3,00 3,39 3,32 2,81 0,640
1 0,2 1,92 3,06 3,38 3,25 2,66 0,628
1 -0,2 1,64 2,94 3,39 3,39 2,95 0,651
0,5 0 3,12 3,96 4,32 4,25 3,77 0,954
0,5 0,2 3,23 4,00 4,31 4,18 3,64 0,941
0,5 -0,2 3,00 3,91 4,32 4,31 3,88 0,966
0,3 0 4,63 5,42 5,70 5,55 4,99 1,42
0,3 0,2 4,74 5,48 5,71 5,52 4,91 1,41
0,3 -0,2 4,51 5,35 5,67 5,58 5,07 1,42
Для задачи Б значения давления и силы меньше, чем для задачи А, как и должно быть. В обеих задачах контактное давление меньше на той стороне области контакта, которая ближе к участку, расположенному между штампами.
Таблица 2
Значения давления и вдавливающей силы Р в задаче Б
X 0(х,0) х103 Рх103
*=-0,075 -0,0375 0 0,0375 0,075
1 0 1,60 2,89 3,28 3,20 2,68 0,595
1 0,2 1,74 2,95 3,27 3,13 2,53 0,585
1 -0,2 1,46 2,82 3,27 3,27 2,82 0,605
0,5 0 2,79 3,65 4,02 3,95 3,46 0,851
0,5 0,2 2,92 3,70 4,01 3,89 3,34 0,842
0,5 -0,2 2,66 3,59 4,01 4,00 3,57 0,858
0,3 0 4,00 4,81 5,13 5,02 4,50 1,22
0,3 0,2 4,13 4,89 5,17 5,02 4,46 1,23
0,3 -0,2 3,88 4,73 5,07 5,01 4,53 1,21
Заключение. Решены новые пространственные контактные задачи с неизвестной областью контакта для упругого слоя о взаимодействии двух одинаковых эллиптических в плане штампов с учетом трения при различных граничных условиях на другой грани. При использовании метода нелинейных граничных интегральных уравнений с учетом симметрии задач определены области контакта, давления в этих областях, связи между силами и осадками штампов. Сделаны расчеты при разных значениях относительной толщины слоя и коэффициента трения.
Библиографический список
1. Чебаков М.И. Пространственная контактная задача для слоя с учетом трения в неизвестной области контакта / М.И. Чебаков // Доклады РАН. - 2002. - Т. 383. - № 1. - С. 67-70.
2. Чебаков М.И. Трехмерная контактная задача для слоя с учетом трения в неизвестной области контакта / М.И. Чебаков // Известия РАН. Механика твердого тела. - 2002. - № 6.
- С. 59-68.
3. Чебаков М.И. Пространственные контактные задачи для слоя с учетом сил трения в зоне контакта / М.И. Чебаков, Х. Лоренц // Современные проблемы механики сплошной среды: тр. 6-й междунар. науч. конф. 19-23 октября 2000 г. - Ростов н/Д: Изд-во СКНЦ ВШ, 2000.
- С. 232-235.
4. Чебаков М.И. Учет сил трения в пространственной контактной задаче для закрепленного слоя / М.И. Чебаков // Современные проблемы механики сплошной среды: тр. 7-й междунар. науч. конф. памяти акад. РАН И.И. Воровича, 22-25 октября 2001 г. - Ростов н/Д: Изд-во СКНЦ ВШ, 2001. - С. 205-209.
5. Александров В.М. Аналитические методы в контактных задачах теории упругости / В.М. Александров, М.И. Чебаков. - М.: Физматлит, 2004. - 301 с.
6. Соболь Б.В. Пространственная задача о контакте системы штампов с упругим слоем / Б.В. Соболь, И.М. Пешхоев // Экологический вестник научных центров ЧЭС. - 2011. - № 1.
- С. 69-76.
7. Галанов Б.А. Метод граничных уравнений типа Гаммерштейна для контактных задач теории упругости в случае неизвестных областей контакта / Б.А. Галанов // Прикладная математика и механика. - 1985. - Т. 49. - Вып. 5. - С. 827-835.
8. Александров В.М. Неклассические пространственные задачи механики контактных взаимодействий упругих тел / В.М. Александров, Д.А. Пожарский. - М.: Факториал, 1998. - 288 с.
9. Прудников А.П. Интегралы и ряды. Специальные функции / А.П. Прудников, Ю.А. Брыч-ков, О.И. Маричев. - М.: Наука, 1983. - 752 с.
Материал поступил в редакцию 30.05.11.
References
1. Chebakov M.I. Prostranstvennaya kontaktnaya zadacha dlya sloya s uchyotom treniya v neiz-vestnoj oblasti kontakta / M.I. Chebakov // Doklady' RAN. - 2002. - T. 383. - # 1. - S. 67-70.
- In Russian.
2. Chebakov M.I. Tryoxmernaya kontaktnaya zadacha dlya sloya s uchyotom treniya v neizvestnoj oblasti kontakta / M.I. Chebakov // Izvestiya RAN. Mexanika tvyordogo tela. - 2002. - # 6.
- S. 59-68. - In Russian.
3. Chebakov M.I. Prostranstvenny'e kontaktny'e zadachi dlya sloya s uchyotom sil treniya v zone kontakta / M.I. Chebakov, X. Lorencz // Sovremenny'e problemy' mexaniki sploshnoj sredy': tr. 6-j mezhdunar. nauch. konf. 19-23 oktyabrya 2000 g. - Rostov n/D: Izd-vo SKNCZ VSH, 2000.
- S. 232-235. - In Russian.
4. Chebakov M.I. Uchyot sil treniya v prostranstvennoj kontaktnoj zadache dlya zakreplyonnogo sloya / M.I. Chebakov // Sovremenny'e problemy' mexaniki sploshnoj sredy': tr. 7-j mezhdunar. nauch. konf. pamyati akad. RAN I.I. Vorovicha, 22-25 oktyabrya 2001 g. - Rostov n/D: Izd-vo SKNCZ VSH, 2001. - S. 205-209. - In Russian.
5. Aleksandrov V.M. Analiticheskie metody' v kontaktny'x zadachax teorii uprugosti / V.M. Aleksandrov, M.I. Chebakov. - M.: Fizmatlit, 2004. - 301 s. - In Russian.
6. Sobol' B.V. Prostranstvennaya zadacha o kontakte sistemy' shtampov s uprugim sloem / B.V. Sobol', I.M. Peshxoev // E'kologicheskij vestnik nauchny'x centrov CHE'S. - 2011. - # 1.
- S. 69-76. - In Russian.
7. Galanov B.A. Metod granichny'x uravnenij tipa Gammershtejna dlya kontaktny'x zadach teorii uprugosti v sluchae neizvestny'x oblastej kontakta / B.A. Galanov // Prikladnaya matematika i mexanika. - 1985. - T. 49. - Vy' p. 5. - S. 827-835. - In Russian.
8. Aleksandrov V.M. Neklassicheskie prostranstvenny'e zadachi mexaniki kontaktny'x vzaimo-dejstvij uprugix tel / V.M. Aleksandrov, D.A. Pozharskij. - M.: Faktorial, 1998. - 288 s. - In Russian.
9. Prudnikov A.P. Integraly' i ryady'. Special'ny'e funkcii / A.P. Prudnikov, Yu.A. Bry'chkov, O.I. Marichev. - M.: Nauka, 1983. - 752 s. - In Russian.
THREE-DIMENSIONAL CONTACT PROBLEM ON INTERACTION BETWEEN ELASTIC LAYER AND TWO PUNCHES UPON FRICTION
D.A. POZHARSKIY, A.A. MOLCHANOV
(Don State Technical University)
Spatial contact problems are investigated for the finite elastic layer in one face of which two identical rigid elliptical punches are symmetrically indented taking into account friction under various types of boundary conditions on the other face. The problems are reduced to the integral equations relating to the contact pressure. The problems are solved by Galanov's method.
Keywords: theory of elasticity, contact problem, layer, friction, Galanov's method.
