CC BY
15
Пожарский Д.А., Золотов Н.Б. Контактные задачи для трансверсально изотропного слоя // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2022. - № 2. - С. 105-113. DOI: 10Л5593/регт.тесУ2022.2Л0
Pozharskii D.A., Zolotov N.B. Contact problems for a transversely isotropic layer. PNRPUMechanics Bulletin, 2022, no. 2, pp. 105-113. DOI: 10.15593/perm.mech/2022.2.10
пермскии политех
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2,2022 PNRPU MECHANICS BULLETIN
http s://ered.pstu.ru/index .php/mechanics/index
Научная статья
DOI: 10.15593/perm.mech/2022.2.10 УДК 539.3
КОНТАКТНЫЕ ЗАДАЧИ ДЛЯ ТРАНСВЕРСАЛЬНО-ИЗОТРОПНОГО СЛОЯ Д.А. Пожарский1, Н.Б. Золотов1'2
Донской государственный техническим университет, Ростов-на-Дону, Россия 2Южный федеральный университет, Ростов-на-Дону, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 07 февраля 2022 г. Одобрена: 23 июня 2022 г. Принята к публикации: 04 июля 2022 г.
Ключевые слова:
теория упругости, контактные задачи, интегральные уравнения, трансверсальная изотропия, слой, скользящая заделка.
Рассматриваются две пространственные, одна осесимметричная и две плоские контактные задачи для трансверсально-изотропного упругого слоя, одна грань которого находится в условиях скользящей заделки. В пространственных и плоских контактных задачах плоскости изотропии могут быть параллельны или перпендикулярны граням слоя. В случае осевой симметрии плоскости изотропии параллельны граням слоя. При помощи интегрального преобразования Фурье контактные задачи сводятся к интегральным уравнениям относительно контактного давления, из которых предельными переходами можно получить известные уравнения соответствующих задач для изотропного слоя. Для решения пространственных задач с неизвестной областью контакта применяется метод нелинейных граничных интегральных уравнений, который позволяет одновременно определить контактные давления и область контакта. При выделении главной части ядра интегрального уравнения пространственной задачи, когда плоскости изотропии перпендикулярны граням слоя, используется полученное ранее в форме, свободной от квадратур, ядро интегрального уравнения соответствующей контактной задачи для трансверсально-изотропного полупространства. Интегральное уравнение осесимметричной задачи при помощи метода парных уравнений сводится к интегральному уравнению Фредгольма второго рода, для численного решения которого применяется метод механических квадратур. Получены замкнутые решения плоских контактных задач, основанные на специальных аппроксимациях символов ядер. Точность аппроксимаций возрастает с увеличением анизотропии. При этом уровень анизотропии можно охарактеризовать отличием отношения корней характеристического уравнения от единицы, соответствующей изотропному случаю. Расчеты механических характеристик и погрешностей аппроксимаций сделаны для известных трансверсально-изотропных материалов.
© ПНИПУ
© Пожарский Дмитрий Александрович - д.ф.-м.н., проф., зав. каф., e-mail: [email protected], : 0000-0001-6372-1866.
Золотов Никита Борисович - асс., e-mail: [email protected], : 0000-0002-2193-0616.
Dmitrii A. Pozharskii - Doctor of Physical and Mathematical Sciences, Professor, Head of Department, e-mail: [email protected], : 0000-0001-6372-1866.
Nikita B. Zolotov - Assistant, e-mail: [email protected], : 0000-0002-2193-0616.
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
CONTACT PROBLEMS FOR A TRANSVERSELY ISOTROPIC LAYER D.A. Pozharskii, N.B. Zolotov
Don State Technical University, Rostov-on-Don, Russian Federation Southern Federal University, Rostov-on-Don, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 07 February 2022 Approved: 23 June 2022 Accepted for publication: 04 July 2022
Keywords:
contact problems, elasticity theory, integral equations, transversal isotropy, layer, sliding support.
Two spatial, one axisymmetric and two plane contact problems are considered for a transversely isotropic elastic layer with one face subjected to sliding support. In the spatial and plane contact problems, the planes of isotropy may be either parallel or perpendicular to the layer faces. In the case of axial symmetry, the planes of isotropy are parallel to the layer faces. By using Fourier integral transforms, the contact problems are reduced to integral equations with respect to the contact pressure, the limiting cases of which are the well-known equations of the corresponding problems for an isotropic layer. For solving the spatial problems with unknown contact domains, the nonlinear boundary integral equations method is used, which make it possible to determine the contact pressure and the contact domain simultaneously. To extract the kernel principal part of the spatial problem integral equation when the isotropy planes are perpendicular to the layer faces, it is used the kernel of the integral equation of the corresponding contact problem for a transversely isotropic half-space obtained earlier without quadratures. The integral equation of the axially symmetric problem is reduced to a Fredholm integral equation of the second kind with the help of the method of pair equations, and the method of mechanical quadratures is used for numerical solutions. Plane problems are solved in a closed form based on special approximations of the kernel symbols. The approximations accuracy grows as anisotropy increases. Here, the anisotropy level can be characterized by the difference between ratio of a characteristic equation roots and unit because the unit value corresponds to the isotropic case. Mechanical characteristics as well as errors of the approximations are calculated for well-known transversely isotropic materials.
© PNRPU
В современном мире контактные задачи возникают постоянно, например, при касании пальцем экрана смартфона [1], пальпационной томографии [2]. Задачи механики анизотропных материалов, контактные задачи важны для развития современной промышленности и занимают важное место в современной мировой науке [3-9]. Аналитические методы и подходы к решению контактных задач были развиты в исследованиях[10-14]. Контактные задачи для функционально-градиентных материалов исследовались в работах [15-18]. Проблемы множественного, дискретного и периодического контакта рассматривались в работах [19-21]. Учет дополнительных факторов при контактном взаимодействии (температура, адгезия) проводился в статьях [22-25]. Изучались контактные задачи для трансверсально-изотропного полупространства как с параллельными [26] или перпендикулярными [27; 28] границе плоскостями изотропии, так и с произвольно ориентированными плоскостями изотропии [29]. Отметим, что для случая параллельных границе плоскостей ядро интегрального уравнения контактной задачи отличается только множителем от соответствующего уравнения для изотропного случая, тогда как для случая перпендикулярных границе плоскостей ядро усложняется и теряет симметрию, осесимметричная формулировка становится невозможной [27]. В работе [28] ядро интегрального уравнения «нетрадиционной» [27] контактной задачи для трансверсально-изотропного полупространства получено в форме, свободной от квадратур, на основе теории обобщенных функций. Исследовались как взаимодействие жестких штампов на транс-
тропном полупространстве [30], так и контакт двух упругих транстропных тел [31]. Асимптотические решения контактных задач для изотропного слоя приведены в работе [11]. Задача о действии кругового штампа на транс-тропный слой (плоскости изотропии параллельны его граням) на основе оригинального метода последовательных приближений рассматривается в исследовании [32]. Интегральные уравнения контактных задач для транс-тропного тела связаны с интегральными уравнениями задач о трещинах в подобных телах [33]. В настоящей статье фундаментальное решение для транстропного полупространства в форме, свободной от квадратур [28], используется для выделения главной части ядра интегрального уравнения трехмерной контактной задачи для транстропного слоя, что позволяет эффективно рассчитать ядро. В связи с тем что эффективность регулярного асимптотического метода решения контактных задач [11], как правило, падает с ростом анизотропии, для плоских задач предлагается метод специальной аппроксимации символа, приводящий к замкнутому решению, эффективность которого увеличивается при возрастании анизотропии.
1. Пространственные контактные задачи
Рассмотрим трансверсально-изотропный упругий слой {|х| < ж, У|<ж, 0 < г < К}, в верхнюю грань которого без перекоса вдавливается жесткий штамп под действием силы Р. Форма основания штампа описывается функцией / (х, у), осадка штампов равна 5. Нижняя
P
P
Рис. 1. Схема контакта в задачах А (а) и Б (b) Fig. 1. Scheme of contact in problems A (a) and B (b)
грань слоя лежит без трения на недеформируемом основании (скользящая заделка). Плоскости изотропии z = const параллельны граням слоя (задача А, рис. 1, а). Рассмотрим аналогичную задачу Б для слоя {0 < x < h, у\ < да, |z| < да}, которая отличается от задачи А направлением осей координат и тем, что плоскости изотропии z = const перпендикулярны граням слоя (рис. 1, b).
Для обеих задач закон Гука имеет вид [4]
. dux , Л „ л ,дЫу . дм
стx = Ац-^ + (А -2AJ—^ + ,
dx ду dz
„ „ ,диг л dUy duz
ст = (Д. -2А66)—х + Аи—у + А13—^
у V 11 66 / _ 11^ 13
дх ду dz
(1)
. дих 4 диу дм ст z = Аъ-^+А1з^г+Азз^ дх дУ да
л дих , диу Хху = Аб^—; + А^ —
ху 66
ду
дх
x,z = А^ + А ^
yz 44
dz
ду
= дих дuz
Txz = А44 + А44
dz дх
В частном случае изотропного тела в формулах (1) следует положить (О — модуль сдвига, V — коэффициент Пуассона)
А11 = Азз =
2G (1 -v) 1 - 2v '
А13 =
2Gv 1 - 2v'
А44 = А66 = G. (2)
Граничные условия задачи А имеют вид (О - область контакта)
z = 0: Uz = 5- f (х, у), (х, у) еП;
стz = о, (х у) ^ = х^ = о,
z = h : uz = - „ = - „ = 0 ,
z xz yz у
а для задачи Б
х = 0: их =5- f (у, z), (у, z) еП;
Стх = ° (у, z) iQ; ^ = -ху = ^
х = h : их = - ^ = - ху = 0 .
(3)
(4)
При заданных параметрах упругости Ати, известных величинах h, 5, функции f (х,у) (задача А) или f (у,2) (задача Б) требуется определить область контакта П и контактное давление ctz (х, у, 0) = -q (ху) ((ху) еП, задача А) или стх (0,у, z) = -q (y,z) ((у^) еП, задача Б), затем при использовании интегрального условия равновесия штампов может быть найдена сила P.
Для сведения задач с граничными условиями (3) и (4) к интегральным уравнениям следует рассмотреть вспомогательные задачи типа Буссинеска о действии заданной сосредоточенной силы на грани слоя при разных ориентациях плоскостей изотропии. Решения таких задач находятся при помощи двойных преобразований Фурье. Получив вспомогательные решения, интегрируем выражения для нормальных перемещений по области контакта и удовлетворяем первому граничному условию (1.3) или (1.4). Для задачи А придем к интегральному уравнению ((ху) еП)
И q(5, л)К(х -5, у - r[)d%dл = 2^9h[5 - f (х, у)],
K(х,у) = - jj ^^cosf s^lcosfД! dsdt, (5) n 00 w V h) У h)
4( w) =
У1 -Y2
y1cth( w / y1) -y 2cth(w / у 2)
e = А А33 Аз, w = Vs2 +t2.
А11(У1 +y 2 )
Для задачи Б аналогичное уравнение имеет вид ((у^) еП)
jj q(5, л)К (у -5, z - Л )d 5d л = 2пА66 h[5 - f (у, z)],
П
К1(у, z) = - JJm(s, t)cos f s- I cos f t— I dsdt,
^ 00 v h ) v h )
M(s,t) =
(m2 - т1)УзsXC2 D
(6)
b
а
D = m2hC2cth^j - m^C^cthC2 + 4(m1 - m2)t2C^C2C3cthC3
2 - a
A , h = (ml +1)у2s2 + 2t2 (l = 1,2),
m, =
А11У2 - A
A13 + A44
с, = ^2s2 +12 (l = 1,2,3).
Здесь Yz (1=1,2) — корни характеристического уравнения (возможно, комплексно сопряженные, Re Yг>0) [4; 27]
У4A11A44 -у2[AuA33 -A13(A13 + 2A44)] + A33A44 = 0, У 3 ЧA44 / Л6 .
(7)
Ввиду положительной определенности матрицы упругих постоянных 9>0 [5]. При переходе к изотропному материалу в формулах (5) следует учесть соотношения (2) и пределы
lim L (w) = L(w) = ~1 , lim 9 = —, (8) ti ^У2 sh2w + 2w t, ^2 1 -v
что приводит к известному интегральному уравнению контактной задачи для изотропного слоя [11]. Соответствующий предельный переход возможен и из формул (6). При к^ж функции-символы ядер (5) и (6) переходят в известные для транстропного полупространства
[27; 28; 32] (L,(w) переходит в 1, в M (s, t) исчезают котангенсы).
Для символа в задаче А выполняется асимптотика
A(w) -1 = 0(exp(-2 / t,)) (w ^ ж, T1 > T2). (9)
Символ в задаче Б обладает важными асимптотическими свойствами: поведение в бесконечности по одному из аргументов связано с поведением в нуле по другому аргументу. Именно, если один из аргументов фиксирован,
в
M(s,t) L(t) (s ^ 0),
В =
y3(m2 - m1)
2(у3 (m2 -m^ -m2у^ + шху2)'
M(s, t)s ^ B2 (s ^ro),
(m2 - m1) У1У 2
В, =
(10)
2 y32(m1 + m2 + 2)(У2 ^^
M (s, t) ^ B L2(s) (t ^ 0);
s
M(s, t)t ^ B1 (t ^t»),
У1 -У 2
L2(. s) =
y1cth(y 2 s) -у 2cth(y1s)
Для решения задач А и Б применим метод нелинейных граничных интегральных уравнений (ГИУ) [14; 34; 35], в соответствии с которым уравнения (5), (6) по неизвестной области О сводятся к нелинейным уравнениям по заданному прямоугольнику со сторонами 2Ь и
2с (Ь > с), содержащему О. При численном решении нелинейных уравнений используется модифицированный метод Ньютона, область О определяется по точкам разбиения 5, в которых решение положительно. В ходе метода требуется рассчитать ядра (5), (6) в узлах решетки. Для эффективного проведения расчетов ниже выделим из ядер главные члены.
Введем безразмерные обозначения
x y z x*= —, y* =—, z* =—, * b 7 b b
. . q( x, y) q(y, z) c
q* (x*, y*) = , q* (y*, z*) ', s=- , (11)
S*=-, J = h, * b b
f*(x*,y*) = f*(y*,z*) = ^, Q^Q*
и т.д., звездочки далее опускаем. Параметр X характеризует относительную толщину слоя.
Используя известную технику [14], асимптотику (9) и обозначения (11), преобразуем ядро (5) для задачи А к форме
К(х,у) = К+Х}[!,(«)-1У0¿и, К Чх2 + у2, (12)
где ^ (и) — функция Бесселя [36].
Для задачи Б, используя ядро соответствующей контактной задачи для полупространства в форме свободной от квадратур [28], представим ядро в виде
К,( у, г) = М,(у, г) +^ }} М (*, 1) -
-M1(s, t)]cos I s j-j cos I tj1 dsdt,
M 1(s,t) =
(m2 - m1)у2s 2C1C 2 D1
(13)
D1 = m2h12C2 - m1h22C1 + 4(m1 - m2)t2C1C2C3
Важно, что в формулах (12) и (13) функции в квадратных скобках экспоненциально убывают на бесконечности.
Для расчетов выберем штампы в форме круговых параболоидов /(х, у) = А(х2 + у2) (задача А) и /(у, г) = = А (у2 + г2) (задача Б). Возьмем следующие транс-тропные материалы (параметры упругости Ар см. [4, с. 22-23]): древесина (ель Дугласа, Yl2 = 13,79, Y2 = 0,1227, Yll = 0,7101, т1 = 76,22, ш2 = 1/ т1), гнейс влажный (Y2 = 2,627, Y2 = 0,2999, Yз2 = 0,5660, т1 = 4,085, т2 = 1/ т1). На рис. 2 и 3 показаны зависимости безразмерной интегральной характеристики контактных
а b
Рис. 2. Графики зависимостей P0 (7) для древесины (а) и гнейса влажного (b) для задач А (сплошные линии) и Б (пунктир)
при 5=А=В=10-3
Fig. 2. Plots of relation P0 (7) for wood (a) and gneiss rock water-saturated (b) for problems A (solid lines) and B (dashed) for 5=А=В=10-
а b
Рис. 3. Графики зависимостей P0 (5) для древесины (а) и гнейса влажного (b) для задач А (сплошные линии), Б (пунктир)
при 7=1, А=В=10-3
Fig. 3. Plots of relation P0 (5) for wood (a) and gneiss rock water-saturated (b) for problems A (solid lines) and B (dashed) for 5=А=В=10-
давлений Р0 = Р / 0 от относительной толщины слоя X
и осадки штампа 8 соответственно (в = 1). Как видим, вдавливающая сила может быть больше в одной из задач в зависимости от материала слоя. Интегральная характеристика уменьшается с ростом X в связи с усилением влияния скользящей заделки нижней грани слоя (рис. 2) и возрастает с ростом осадки (рис. 3).
2. Осесимметричная контактная задача
Частный случай задачи А, когда выполнены условия осевой симметрии, назовем задачей В. Используя обозначения (11) при в = 1, вводя полярные координаты [11], преобразуем интегральное уравнение задачи В к безразмерной форме
а л
{<70(рЖ0(р, г)рф= —[8-/о(г)], 0 < г < а, (14)
К^ r) = j L1(u) J0 ( и 7 J0 fи 71du,
где а - радиус области контакта, д0 (р) - контактное давление, /0 (г) - форма основания штампа.
Вводя новые функции А(и) и ф(?) по формулам
А(и) = ^(и)| ^0(Р)^0 ^и р^ р^р = А(и)} Ф(')cos ^и -XX) Ж , сведем уравнение (14) к парным уравнениям
j А(и) J0 Vи 7 ) du = 2-[5 - f0 (r)], 0 < r < a, (15)
j А(и)—U—J0 ( и— Idu = 0, a < r <да.
j v \{u) 0V *J
Применяя для решения парных уравнений (15) известный метод [10], приходим к интегральному уравнению Фредгольма второго рода
ф(я)-""^j )j[1 - Li(u )]cos ^ и 7 ) cos f и dudt = F (s),
0 < s < a,
F (s) = -
8 - /o(°) - s J / (s sin y)dу
(16)
для численного решения которого применим метод механических квадратур.
При условии ограниченности давления на краю области контакта неизвестный радиус а можно определить из условия [10]
ф(а) = 0.
(17)
Зная решение уравнения (16), безразмерную интегральную характеристику контактных давлений можно найти по формуле [10]
а
P0 = 2я J ф^ )dt.
(18)
Расчеты интегральной характеристики по формулам (16)-(18) показывают хорошее совпадение с численными результатами для задачи А, полученными при помощи метода нелинейных ГИУ (отличие около 1 %). На рис. 4 показаны зависимости радиуса области контакта а, найденного из условия (17), от осадки 5 для кругового параболоида (/0(г) = Аг2). При X = ж (полупространство) упругие характеристики не влияют на зависимость а (5), поскольку ядро интегрального уравнения (16) тождественно обращается в нуль. В случае слоя (X = 1) область контакта расширяется, а значения а (5) для древесины превышают соответствующие значения для изотропного материала.
0,9 -
0,6 -
0,3 -
0,2
0,6
1,0
1038
Рис. 4. Графики зависимостей a (8), задача В, при A=10-3: Х=ю - сплошная линия; Х=1 - пунктир для изотропного материала, точки для древесины
Fig. 4. Plots of relation a (8) in problem C for A=10-3: solid line for Х=ю; dashed for isotropic material (X=1), dotted for wood (X=1)
3. Плоские контактные задачи
Пусть трансверсально-изотропные слои, изображенные на рис. 1, а и 1, b, взаимодействует с жесткими полосовыми штампами соответственно по заданным областям
{| х |< а, | у |<ж} (задача Г) и {| г |< а, | у |<ж} (задача Д). Формы оснований штампов в задача Г и Д описываются соответственно функциями / (х) и / (г). Интегральные уравнения контактных задач Г и Д относительно соответственно контактных давлений q (х) и q (г) можно вывести из уравнений (5) и (6), используя теорию обобщенных функций, а для задачи Д еще и асимптотику (10) функции М (5, 1) при 1 ^ 0 . После введения безразмерных обозначений Х=К/а,
Г) х, = х / а, (х,) = [5-/ (х)]/а, ф(х,) = q(x)/9,
(19)
Д) х.= г/а, я,(х.) = [5-/(г)]/а, ф(х.) = q(г)/(Б2А66),
(звездочки далее опускаем) интегральные уравнения обеих плоских контактных задач можно записать в форме (п=1 для задачи Г, п=2 для задачи Д)
|ф©к
1-х
d^=ng(х) (| х |< 1),
к (t) = J Ln (u)cos(ut)du.
(20)
Здесь функции-символы ядер определяются по формулам (5) и (10).
Замечание. Если для слоя, показанного на рис. 1, Ь, полосовой штамп направить вдоль оси г, то в каждом плоском сечении, перпендикулярном оси штампа, будет изотропия упругих свойств. В этом случае на основании асимптотики (10) функции М (5,1) при 5 ^ 0 интегральное уравнение плоской контактной задачи также будет иметь вид (19), где следует в качестве символа взять функцию Ь (и) из формул (8), которая соответствует контактной задаче для изотропного слоя [11]. В результате интегральное уравнение будет отличаться от известного [11] только множителем, учитывающим анизотропию (он исключается введением обозначений типа (20)), для решения такой задачи можно рекомендовать асимптотические методы [11].
Регулярный асимптотический метод [11], основанный на разложении решения в ряд по степеням малого геометрического параметра 1/Х, может применяться при Х>1 в контактной задаче для изотропного слоя, при Х> тах(|Y1|,^2|) в задаче Г или Х> тах(|Y-1|,^-'|) в задаче Д, когда сходятся соответствующие разложения ядер (20). С ростом анизотропии (действительные корни Y1 >> 1 >> Y2 > 0) область применимости метода в задачах Г и Д существенно снижается. Заметим, что в этом случае членом с одним из котангенсов в функциях-символах Ьп (и) ядер (20) можно пренебречь (Y2сШ(м / Y2) в задаче Г и Y гС&^м) в задаче Д), что позволяет получить решение уравнения (20) в замкнутой форме. Именно при учете асимптотического поведения символов Ьп (и) в нуле и бесконечности аппроксимируем их при 0 < и < ж выражением
a
L(u) И th(Bu); Г) B =
-, Д) B = _Ü!Z2_ . (21)
Y1 +Т2 Yl +Y2
Интегральное уравнение (20), (21) сводится к сингулярному интегральному уравнению с ядром Коши, имеющему точное решение [14]. Для штампов с плоским основанием (/ (х) =/=сош^) приближенные решения контактных задач Г и Д имеют вид
Ф( х) = КГ) [2 ((2 А.) - сИ(2 А. х) )]Л
K (e)
J ф( x)dx =
fK (e) K (e) '
(22)
А=л/(2£7), d =expA, e = exp(-2 А), К1 (e) = К (VT-ё1),
где К (e) — полный эллиптический интеграл. Относительная погрешность решения (22) не превышает погрешности аппроксимации (21) и уменьшается с ростом анизотропии. Важно, что решение (22) справедливо для любой относительной толщины слоя 7.
Корни характеристического уравнения (7) и погрешность аппроксимации (21)
Roots of the characteristic equation (7) and error of the approximation (21)
Материал Y1 Y2 Y1/Y2 8, %
Изотропный 1,000 1,000 1,00 9,5
Гнейс влажныш 1,621 0,5476 2,96 7,1
Композит (60 % волокон) 4,724 0,6888 6,86 4,1
Ламинат T650/950-1 5,413 0,6458 8,38 3,5
Древесина (ель Дугласа) 3,713 0,3504 10,6 2,8
В таблицу включено несколько известных трансвер-сально-изотропных материалов ([4, табл. 1.4]), а также изотропный материал, указана относительная погрешность в аппроксимации (21), практически одинаковая для задач Г и Д. Как видно из данных таблицы, относительная погрешность в монотонно снижается с ростом параметра у1 / у2 (с увеличением анизотропии).
Заключение
Показана связь интегральных уравнений пространственных контактных задач для трансверсально-изотроп-ного слоя с уравнениями соответствующих плоских задач и задачи при осевой симметрии. Можно сделать вывод о существенном влиянии ориентации плоскостей изотропии на характеристики контакта. Для случая перпендикулярных граням слоя плоскостей изотропии фундаментальное решение соответствующей задачи для транстропного полупространства в форме, свободной от квадратур, позволяет эффективно регуляризовать ядро интегрального уравнения контактной задачи для слоя. При уменьшении X (относительной толщины слоя) возрастает отличие контактных характеристик от соответствующих величин для случая трансверсально-изотропного полупространства. В случае осевой симметрии радиус заранее неизвестной круговой области контакта для слоя оказывается больше, чем в соответствующей задаче для полупространства, а для транстропного слоя область контакта может быть больше, чем для изотропного слоя. В то время как эффективность регулярного асимптотического метода для рассмотренных плоских контактных задач снижается при увеличении анизотропии, эффективность метода, основанного на специальной аппроксимации ядра интегрального уравнения, наоборот, возрастает с ростом анизотропии.
Библиографический список
1. Argatov I.I., Borodich F.M. A macro model for electroad-hesive contact of a soft finger with a touchscreen // IEEE Transactions on Haptics. - 2020. - Vol. 13, no. 3. - P. 504-510. DOI: 10.1109/TOH.2020.2969628
2. Argatov I.I., Jin X.Q., Keer L.M. Collective indentation as a novel strategy for mechanical palpation tomography // Journal of the Mechanics and Physics of Solids. - 2020. - Vol. 143. -P. 104063. DOI: 10.1016/jjmps.2020.104063
3. Лехницкий С.Г. Теория упругости анизотропного тела. -М.: Наука, 1977. - 416 с.
4. Ding H., Chen W., Zhang L. Elasticity of transversely isotropic materials. - Dordrecht: Springer, 2006. - 435 p.
5. Pan E., Chen W. Static Green's functions in anisotropic media. - New York: Cambridge University Press, 2015. - 356 p.
6. Popov V.L., Heß M. Method of dimensionality reduction in contact mechanics and friction. - Berlin: Springer, 2015. - 265 p. DOI: 10.1007/978-3-642-53876-6
7. Argatov I., Heß M., Pohrt R., Popov V.L. The extension of the method of dimensionality reduction to non-compact and non-axisymmetric contacts // ZAMM. - 2016. - Vol. 96, no. 10 -P. 1144-1155. DOI: 10.1002/zamm.201600057
8. Barber J.R. Contact mechanics. - Berlin: Springer, 2018. -585 p. DOI: 10.1007/978-3-319-70939-0
9. Argatov I. From Winkler's foundation to Popov's foundation // Facta Universitatis. Series: Mechanical Engineering. - 2019. Vol. 17, no 2. - P. 181-190. DOI: 10.22190/FUME190330024A
10. Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. - М.-Л.: Изд-во АН СССР, 1963. 368 с.
11. Ворович И.И., Александров В.М., Бабешко В.А. Неклассические смешанные задачи теории упругости. - М.: Наука, 1974. - 456 с.
12. Александров В.М., Коваленко Е.В. Задачи механики сплошных сред со смешанными граничными условиями. - М.: Наука, 1986. - 336 с.
13. Александров В.М., Ромалис Б.Л. Контактные задачи в машиностроении. - М.: Машиностроение, 1986. - 176 с.
14. Alexandrov V.M., Pozharskii D.A. Three-dimensional contact problems. - Dordrecht: Kluwer, 2001. - 406 p.
15. Аналитические решения смешанных осесимметрич-ных задач для функционально-градиентных сред / С.М. Айзи-кович, В.М. Александров, А.С. Васильев, Л.И. Кренев, И.С. Трубчик. - М.: Физматлит, 2011. - 192 с.
16. Напряженно-деформированное состояние упругого мягкого функционально-градиентного покрытия при внедрении сферического индентора / С.С. Волков, А.С. Васильев, С.М. Айзикович, Н.М. Селезнев, А.В. Леонтьева // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 4. - С. 20-34. DOI: 10.15593/peim.mech/2016.4.02
17. Васильев А.С., Волков С.С., Айзикович С.М. Приближенное аналитическое решение задачи о вдавливании проводящего штампа в электроупругое полупространство с неоднородным покрытием // Доклады Академии наук. - 2018. -Т. 478, № 1. - С. 34-39. DOI: 10.7868/S0869565218010073
18. Yastrebov V.A. Anciaux G., Molinari J.-F. From infinitesimal to full contact between rough surfaces: evolution of the contact area // International Journal of Solids and Structures. -2015. - Vol. 52. - P. 83-102. DOI: 10.1016/j.ijsolstr.2014.09.019
19. Goryacheva I.G., Tsukanov I.Y. Modeling of normal contact of elastic bodies with surface relief taken into account // Journal of Physics: Conference Series. - 2018. - Vol. 991, no. 1 -P. 012028. DOI: 10.1088/1742-6596/991/1/012028
20. Пожарский Д.А. Периодические контактные и смешанные задачи теории упругости (обзор) // Известия вузов. Северо-Кавказский регион. Естественные науки. - 2021. -№ 2. - С. 22-33. DOI: 10.18522/1026-2237-2021-2-22-33
21. Yakovenko A., Goryacheva I. The periodic contact problem for spherical indenters and viscoelastic half-space // Tribology International. - 2021. - Vol. 161. - P. 107078. DOI: 10.1016/j .triboint.2021.107078
22. Индикация термоупругой неустойчивости скользящего контакта с помощью заглубленной пьезокерамической прослойки / В.Б. Зеленцов, Б.И. Митрин, А.Г. Сукиязов, С.М. Ай-зикович // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. -№ 1. - С. 63-84. DOI: 10.15593/perm.mech/2017.1.05
23. Goryacheva I.G., Makhovskaya Y. Combined effect of surface microgeometry and adhesion in normal and sliding contacts of elastic bodies // Friction. - 2017. - Vol. 5, no. 3. - P. 339-350. DOI: 10.1007/s40544-017-0179-1
24. He X., Li Q., Popov V.L. Simulation of adhesive contact of soft microfibrils // Lubricants. - 2020. - Vol. 8, no. 10. - P. 94. DOI: 10.3390/lubricants8100094
References
1. Argatov I.I., Borodich F.M. A macro model for elec-troadhesive contact of a soft finger with a touchscreen. IEEE Transactions on Haptics, 2020, vol. 13, no. 3, pp. 504-510. DOI: 10.1109/TOK2020.2969628
2. Argatov I.I., Jin X.Q., Keer L.M. Collective indentation as a novel strategy for mechanical palpation tomography. Journal of the Mechanics and Physics of Solids, 2020, vol. 143, pp. 104063. DOI: 10.1016/j.jmps.2020.104063
3. Lekhnitskii S.G. Teoriia uprogosti anizotropnogo tela [Theory of elasticity of anisotropic solid]. Мoscow, Nauka, 1977, 416 p.
4. Ding H., Chen W., Zhang L. Elasticity of transversely iso-tropic materials. Dordrecht, Springer, 2006, 435 p.
5. Pan E., Chen W. Static Green's functions in anisotropic media. New York, Cambridge University Press, 2015, 356 p.
6. Popov V.L., Heß M. Method of dimensionality reduction in contact mechanics and friction. Berlin, Springer, 2015, 265 p. DOI: 10.1007/978-3-642-53876-6
7. Argatov I., Heß M., Pohrt R., Popov V.L. The extension of the method of dimensionality reduction to non-compact and non-
25. He X., Li Q., Popov V.L. Strength of adhesive contact between a rough fibrillar structure and an elastic body: influence of fibrillar stiffness // Journal of Adhesion. - 2021. - Published online June 2021. - P. 1-14. DOI: 10.1080/00218464.2021.1939017
26. Vasiliev A.S., Volkov S.S., Aizikovich S.M. Indentation of an axisymmetric punch into an elastic transversely-isotropic half-space with functionally graded transversely-isotropic coating // Materials Physics and Mechanics. - 2016. - Vol. 28, no. 1-2. - P. 11-15.
27. Fabrikant V.I. Non-traditional contact problem for transversely isotropic half-space // Quarterly Journal of Mechanics and Applied Mathematics. - 2011. - Vol. 64, no. 2. - P. 151-170. DOI: 10.1093/qjmam/hbq029
28. Давтян Д.Б., Пожарский Д.А. Действие полосового штампа на трансверсально изотропное полупространство // Прикладная математика и механика. - 2012. - Т. 76, вып. 5. -С. 783-794.
29. Fabrikant V.I. Contact problem for an arbitrarily oriented transversely isotropic half-space // Acta Mechanica. - 2017. -Vol. 228, no. 4. - P. 1541-1560. DOI: 10.1007/s00707-016-1788-x
30. Бедоидзе М.В., Пожарский Д.А. Взаимодействие штампов на трансверсально изотропном полупространстве // Прикладная математика и механика. - 2014. - Т. 78, вып. 4. -С. 576-582.
31. Пожарский Д.А. Контакт транстропных тел в теории Герца // Прикладная механика и техническая физика. - 2018. -Т. 59, № 3. - С. 121-128. DOI: 10.15372/PMTF20180313
32. Fabrikant V.I. Contact and crack problems in linear elasticity. - Sharjah: Bentham, 2010. - 1030 p.
33. Fabrikant V.I. New approach to interface crack problems in transversely isotropic materials // ZAMP. - 2021. - Vol. 72. -P. 86. DOI: 10.1007/s00033-020-01445-y
34. Галанов Б.А. Метод граничных уравнений типа Гам-мерштейна для контактных задач теории упругости в случае неизвестных областей контакта // Прикладная математика и механика. - 1985. - Т. 49, вып. 5. - С. 827-835.
35. Пожарский Д.А. Фундаментальные решения статики упругого клина и их приложения. - Ростов-на-Дону: ООО «ДГТУ-Принт», 2019. - 312 с. https://www.rfbr.ru/ rffi/ru/books/o_2089067
36. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. - М.: Наука, 1971. - 1108 с.
axisymmetric contacts. ZAMM, 2016, vol. 96, no. 10, pp. 1144-1155. DOI: 10.1002/zamm.201600057
8. Barber J.R. Contact mechanics. Berlin, Springer, 2018, 585 p. DOI: 10.1007/978-3-319-70939-0
9. Argatov I. From Winkler's foundation to Popov's foundation. Facta Universitatis. Series: Mechanical Engineering, 2019, vol. 17, no 2, pp. 181-190. DOI: 10.22190/FUME190330024A
10. Ufliand Ia.S. Integral'nye preobrazovaniia v zadachakh teorii uprugosti [Integral transforms in problems of the elasticity theory]. Moscow, Leningrad, Academy of Sciences of USSR Publishing, 1963, 368 p.
11. Vorovich I.I., Aleksandrov V.M., Babeshko V.A. Nek-lassicheskie smeshannye zadachi teorii uprugosti [non-classical mixed problems of the elasticity theory]. Moscow, Nauka, 1974, 456 p.
12. Aleksandrov V.M., Kovalenko E.V. Zadachi mekhaniki sploshnyh sred so smeshannymi granichnymi usloviiami [Problems with mixed boundary conditions in continuum mechanics]. Moscow, Nauka, 1986, 336 p.
13. Aleksandrov V.M., Romalis B.L. Kontaktnye zadachi v mashinostroenii [Contact problems in machine-building]. Moscow, Mashinostroenie, 1986, 176 p.
14. Alexandrov V.M., Pozharskii D.A. Three-dimensional contact problems. Dordrecht, Kluwer, 2001, 406 p.
15. Aizikovich S.M., Aleksandrov V.M., Vasil'ev A.S., Kre-nev L.I., Trubchik I.S. Analiticheskie resheniia smeshannykh ose-simmetrichnykh zadach dlia funktsional'no-gradientnykh sred [Analytical solutions of mixed axially symmetric problems for functionally-graded media]. Moscow, Fizmatlit, 2011, 192 p.
16. Volkov S.S., Vasiliev A.S., Aizikovich S.M., Seleznev N.M., Leontieva A.V. Napriazhenno-deformirovannoe sostoianie uprugogo miagkogo funktsional'no-gradientnogo pokrytiia pri vnedrenii sfericheskogo indentora [Stress-strain state of an elastic soft functionally-graded coating subjected to indentation by a spherical punch]. PNRPU Mechanics Bulletin, 2016, no. 4, pp. 20-34. DOI: 10.15593/perm.mech/2016.4.02
17. Vasiliev A.S., volkov S.S., Aizikovich S.M. Approximated analytical solution of a problem on indentation of an electro-elastic half-space with inhomogeneous coating by a conductive punch. Doklady Physics, 2018, vol. 63, no. 1, pp. 18-22. DOI: 10.1134/S1028335818010020
18. Yastrebov V.A. Anciaux G., Molinari J.-F. From infinitesimal to full contact between rough surfaces: evolution of the contact area. International Journal of Solids and Structures, 2015, vol. 52, pp. 83-102. DOI: 10.1016/j.ijsolstr.2014.09.019
19. Goryacheva I.G., Tsukanov I.Y. Modeling of normal contact of elastic bodies with surface relief taken into account. Journal of Physics: Conference Series, 2018, vol. 991, no. 1, pp. 012028. DOI: 10.1088/1742-6596/991/1/012028
20. Pozharskii D.A. Periodicheskie kontaktnye i smeshannye zadachi teorii uprugosti (obzor) [Periodic contact and mixed problems of the elasticity theory (review)]. Izvestiia Vuzov. Severo-Kavkazskii Region. Estestvennye Nauki, 2021, no, 2, pp. 22-33. DOI: 10.18522/1026-2237-2021-2-22-33
21. Yakovenko A., Goryacheva I. The periodic contact problem for spherical indenters and viscoelastic half-space. Tribology International, 2021, vol. 161, pp. 107078. DOI: 10.1016/j .triboint.2021. 107078
22. Zelentsov V.B., Mitrin B.I., Sukiyazov A.G., Aizikovich S.M. Indikatsiia termouprugoi neustoichivosti skol'ziaschego kon-takta s pomosch'yu zaglublennoi p'ezokeramicheskoi prosloiki [Indication of thermoelastic instability of sliding contact using embedded piezoceramic interlayer]. PNRPU Mechanics Bulletin, 2017, no. 1, pp. 63-84. DOI: 10.15593/perm.mech/2017.1.05
23. Goryacheva I.G., Makhovskaya Y. Combined effect of surface microgeometry and adhesion in normal and sliding con-
tacts of elastic bodies. Friction, 2017, vol. 5, no. 3, pp. 339-350. DOI: 10.1007/s40544-017-0179-1
24. He X., Li Q., Popov V.L. Simulation of adhesive contact of soft microfibrils. Lubricants, 2020, vol. 8, no. 10, pp. 94. DOI: 10.3390/lubricants8100094
25. He X., Li Q., Popov V.L. Strength of adhesive contact between a rough fibrillar structure and an elastic body: influence of fibrillar stiffness. Journal of Adhesion, 2021, published online June 2021. pp. 1-14. DOI: 10.1080/00218464.2021.1939017
26. Vasiliev A.S., volkov S.S., Aizikovich S.M. Indentation of an axisymmetric punch into an elastic transversely-isotropic halfspace with functionally graded transversely-isotropic coating. Materials Physics and Mechanics, 2016, vol. 28, no. 1-2, pp. 11-15.
27. Fabrikant V.I. non-traditional contact problem for transversely isotropic half-space. Quarterly Journal of Mechanics and Applied Mathematics, 2011, vol. 64, no. 2, pp. 151-170. DOI: 10.1093/qjmam/hbq029
28. Davtyan D.B., Pozharskii D.A. The action of a strip punch on a transversely isotropic half-space. Journal of Applied Mathematics and Mechanics, 2012, vol. 76, no. 5, pp. 558-566.
29. Fabrikant V.I. Contact problem for an arbitrarily oriented transversely isotropic half-space. Acta Mechanica, 2017, vol. 228, no. 4, pp. 1541-1560. DOI: 10.1007/s00707-016-1788-x
30. BeDOIdze M.V., Pozharskii D.A. The interaction of punches on a transversely isotropic half-space. Journal of Applied Mathematics and Mechanics, 2014, vol. 78, no. 4, pp. 409-414. DOI: 10.1016/j.jappmathmech.2014.12.012
31. Pozharskii D.A. Contact of transversely isotropic bodies in the Hertz theory. Journal of Applied Mechanics and Technical Physics, 2018, vol. 59, no. 3, pp. 491-497. DOI: 10.1134/ S0021894418030136
32. Fabrikant V.I. Contact and crack problems in linear elasticity. Sharjah, Bentham, 2010, 1030 p.
33. Fabrikant V.I. New approach to interface crack problems in transversely isotropic materials. ZAMP, 2021, vol. 72, pp. 86. DOI: 10.1007/s00033-020-01445-y
34. Galanov B.A. The method of boundary equations of the Hammerstein-type for contact problems of the theory of elasticity when the regions of contact are not known. Journal of Applied Mathematics and Mechanics, 1985, vol. 49, no. 5, pp. 634-640.
35. Pozharskii D.A. Fundamental'nye resheniia statiki uprugogo klina i ikh prilozheniia [Fundamental solutions of statics of an elastic wedge and their applications]. Rostov-on-Don, DSTU-Print, 2019, 312 p. https://www.rfbr.ru/rffi/ru/books/o_2089067
36. Gradshtein I.S., Ryzhik I.M. Tablitsy integralov, summ, riadov i proizvedenii [Tables of integrals, sums, series and products]. Moscow, Nauka, 1971, 1108 p.
Финансирование. Работа выполнена при финансовой поддержке Российского научного фонда (грант № 22-21-00013). Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов.
Financing. This study was supported by the Russian Science Foundation (Grant No. 22-21-00013). Conflict of interest. The authors declare no conflict of interest.
