CC BY
34
Пожарский Д.А., Пожарская Е.Д. Контактные задачи для упругого неоднородного тела с цилиндрической шахтой // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 4. - С. 200-208. DOI: 10.15593/perm.mech/2018.4.18
Pozharskii D.A., Pozharskaya E.D. Contact problems for an elastic inhomogeneous body with a cylindrical cavity. PNRPU Mechanics Bulletin, 2018, no. 4, pp. 200-208. DOI: 10.15593/perm.mech/2018.4.18
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 4, 2018 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01: 10.15593/регш.шееЬ/2018.4.18 УДК 539.3
КОНТАКТНЫЕ ЗАДАЧИ ДЛЯ УПРУГОГО НЕОДНОРОДНОГО ТЕЛА С ЦИЛИНДРИЧЕСКОЙ ШАХТОЙ
Д.А. Пожарский, Е.Д. Пожарская
Донской государственный технический университет, Ростов-на-Дону, Россия
О СТАТЬЕ АННОТАЦИЯ
Изучается осесимметричная задача упругого равновесия непрерывно неоднородного пространства с цилиндрической полостью, когда коэффициент Пуассона является произвольной достаточно гладкой функцией радиальной координаты, а модуль сдвига постоянный. При этом модуль упругости Юнга также является переменным по радиальной координате. Предложено общее представление решения, которое приводит к векторному уравнению Лапласа и скалярному уравнению Пуассона, правая часть которого зависит от коэффициента Пуассона. При помощи интегрального преобразования Фурье построены в квадратурах точные общие решения уравнений Лапласа и Пуассона. Получены интегральные уравнения двух осесимметричных контактных задач о взаимодействии поверхности полости (шахты) с жестким цилиндрическим вкладышем, вставленным в нее с натягом. В первой задаче контакт считается абсолютно гладким, для решения интегрального уравнения первого рода относительно контактного давления используется сингулярный асимптотический метод, эффективный для относительно длинных вкладышей. Во второй задаче учитывается шероховатость поверхности шахты, которая моделируется дополнительной прослойкой винклеровского типа, для решения интегрального уравнения второго рода применяется метод коллокации, эффективный для относительно коротких подкрепляющих вкладышей. Контактное давление на границе области контакта имеет характерную корневую особенность в первой задаче и принимает конечное значение во второй задаче. Для однородного материала отмечается близость интегральных характеристик контактных давлений, получаемых в обеих задачах, при малых показателях шероховатости (коэффициентах постели) в определенном диапазоне относительных длин вкладышей. Показано, что учет шероховатости снижает влияние неоднородности на распределение контактных давлений. Расчеты сделаны для случаев, когда коэффициент Пуассона и модуль упругости возрастают или убывают при удалении от поверхности полости.
©ПНИПУ
Получена: 27 июля 2018 г. Принята: 7 декабря 2018 г. Опубликована: 28 декабря 2018 г.
Ключевые слова:
контактные задачи, теория упругости, неоднородное тело, цилиндрическая полость, асимптотический метод, шероховатость.
© Пожарский Дмитрий Александрович - д.ф.-м.н., проф., зав. каф., e-mail: [email protected], 0000-0001-6372-1866
Пожарская Елизавета Дмитриевна - студентка, e-mail: [email protected], 0000-0002-5745-6135
Dmitrii A. Pozharskii - Doctor of Physical and Mathematical Sciences, Professor, Head of Department, e-mail: [email protected], 0000-0001-6372-1866
Elizaveta D. Pozharskaya - Undergraduate Student, e-mail: [email protected], 0000-0002-5745-6135
Эта статья доступна в соответствии с условиями лицензии 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 AN ELASTIC INHOMOGENEOUS BODY WITH A CYLINDRICAL CAVITY
D.A. Pozharskii, E.D. Pozharskaya
Don State Technical University, Rostov-on-Don, Russian Federation
ABSTRACT
An axially symmetric elastic equilibrium problem is investigated for a continuously inhomo-geneous space with a cylindrical cavity when Poisson's ratio is being an arbitrary fairly smooth function with respect to radial coordinate while shear modulus is constant. For this case Young's modulus is also variable with respect to the radial coordinate. A general solution is suggested which leads us to a vector Laplace equation and a scalar Poisson equation whose right-hand side depends on Poisson's ratio. As a result, exact general solutions of the Laplace and Poisson equations are constructed in integral forms with the help of Fourier transformations. Then integral equations of two axially symmetric contact problems are derived on the interaction between the cavity surface and a rigid cylindrical insert fitted with interference. In the first problem the contact is supposed to be absolutely ideal, a singular asymptotical method is used here to solve the integral equation of the first kind with respect to the contact pressure, which is effective for fairly long inserts. In the second problem the mine surface is supposed to be rough simulated by an extra Winkler type, a collocation method is used for solving the integral equation of the second kind, which is effective for fairly short inserts. The contact pressure has typical square root singularities at end-points in the first problem while it takes finite values at those points in the second problem. For a homogeneous material, the integral characteristics of the contact pressures in both problems are close for small coefficients of roughness for some values of the insert relative length. It is shown that the rough surface distributes contact pressures more uniformly removing effect of nonhomogeneity. The calculations are made for the cases when Poisson's ratio and Young's modulus increase or decrease from the surface of the cavity.
©PNRPU
ARTICLE INFO
Received: 27 July 2018 Accepted: 7 December 2018 Published: 28 December 2018
Keywords:
contact problems, elasticity theory, inhomogeneous body, cylindrical cavity, asymptotical method, roughness.
Введение
Контактные задачи теории упругости для неоднородных (функционально градиентных) материалов [1-10] важны для оценки контактной прочности различных конструкций и деталей. Интегральные уравнения контактных задач взаимосвязаны с интегральными уравнениями задач о трещинах [11]. При рассмотрении однородных тел цилиндрической формы, контактирующих с бандажами и вкладышами, для решения возникающих интегральных уравнений развивались асимптотические методы [12-17]. При решении статических и динамических контактных задач для неоднородных цилиндрических тел превалировали численные методы, обобщенный метод фиктивного поглощения [18-22]. Изучались волновые процессы и задачи идентификации упругих параметров для цилиндров [23, 24]. Учитывались сложные явления на поверхности контакта, включая микрогеометрию и адгезию [25, 26].
При переменном коэффициенте Пуассона точные фундаментальные решения были получены для полупространства и слоя (коэффициент Пуассона непрерывно меняется по глубине) [27, 28] и для плоского клина (коэффициент Пуассона зависит от угловой координаты) [29]. Изучался частный случай изменения коэффициента Пуассона с глубиной полупространства [30]. Для неоднородного цилиндра, когда модуль упругости и коэффициент Пуассона зависят от радиальной коорди-
наты, предлагалось искать решение по методу возмущений [31]. В настоящей работе для пространства с цилиндрической шахтой, когда коэффициент Пуассона непрерывно зависит от радиальной координаты, на основе подхода А.Н. Бородачева [28] получено точное фундаментальное решение, символ ядра интегрального уравнения контактной задачи построен в квадратурах. Для решения контактной задачи предлагается сингулярный асимптотический метод (гладкий контакт) и метод коллокации (шероховатая поверхность контакта).
1. Фундаментальное решение
В цилиндрических координатах г, 2 рассмотрим осесимметричную деформацию упругого пространства с цилиндрической полостью {р<г<да, 21<да}, когда модуль сдвига О постоянный, а коэффициент Пуассона у(г) — произвольная достаточно гладкая функция, удовлетворяющая условию -1<у(г)<0,5 [28]. В этом случае модуль продольной упругости материала Е(г)=2О(1+у(г)) — положительная функция угловой координаты. При отсутствии массовых сил векторное уравнение равновесия в перемещениях для данной модели принимает форму
Ди + Ега<1(у (г)<у и) = 0, (1)
и = (и, w), и = и(г, г), w = w(r, г), у(г) = [1 - 2у(г )]-1 .
Разыскивая общее решение уравнения (1) в форме
u = B + gradé, B = (B, B2 ), Bl = Bl (r, z),
B2 = B2(r, z), b = b(r, z)
(2)
и внося представление (2) в уравнение (1), получаем векторное уравнение Лапласа и скалярное уравнение Пуассона
АВ = 0, АЬ = ^(г^ В, ^(г) = -[2(1 - у(г))]-1 . (3)
Для однородного материала, у(г)=у=сош1, представление общего решения (2), (3) в точности совпадает с известным решением Фрайбергера [32], которое для этого случая при
В1 = 4(1 — у)Вг, Ь = -гВг —В0, В2 = Вг = 0 (4)
сводится к классическому представлению Папкови-ча—Нейбера в цилиндрических координатах через три гармонические функции Bo, Br, Bz [33].
Проекции векторного уравнения Лапласа (3) приводят к уравнениям
АЯ, — В = 0, АВ2 = 0, А = 4 +1 - + 4 . (5)
г дг г дг дх
При помощи интегрального преобразования Фурье, предполагая, что задача симметрична по z, учитывая убывание перемещений на бесконечности, общее решение уравнений (5) получим в виде
да да
В1 = |с1 ^(Рг)СО5(РГ)^Р, В2 = |с2 К0($г)ътф2^ р, (6)
0 0
где Kn(г) - модифицированные цилиндрические функции Бесселя (п=0,1), также называемые функциями Макдональда [34, 35].
Из соотношений (6) найдем дивергенцию в правой части уравнения Пуассона (3), а затем, используя метод вариации произвольных постоянных, построим частное решение этого уравнения в форме (Ь(г) - другая модифицированная цилиндрическая функция Бесселя)
b = J (C - C2)[£^0(Pr) - E2/0(Pr)]pccs(pz)jp,
0
r r
E = j s^(s)I0 (Ps)K (Ps)ds, E = j s^(s)K02 (Ps)ds.
(7)
Величины С1 и ^ в формулах (6) и (7) определяются из двух граничных условий на поверхности цилиндрической полости. На основании закона Гука для напряжений имеем выражения
du dw
т„ = G|--1--
dz dr
а„ = 2Gy(r ) |(1 -v(r)) dU + v(r ) ( - + dW II.
dr l r dz
(8)
Пусть к поверхности полости (г=р) при z=0 приложена осесимметричная нормальная сосредоточенная нагрузка q. Тогда граничные условия имеют вид (5(г) -5-функция Дирака)
r = p : хГ2 = 0, стг = -qô(z).
(9)
На основании формул (2)—(8) и условий (9) приходим к системе двух уравнений относительно ^ и решив которую, найдем используемое в контактных задачах нормальное перемещение точек поверхности полости:
u(p, z) =
qp r K12(P p)cos(Pz)dP
2nG 0 Kf(Pp) - 2F
(10)
F = P21 s^(s)K2(Ps)ds ,
которое при у(г)=у=сош1 совпадает с известным для однородного материала [12—15]
qp(1 -v) '
u (p, z) = Kf(Pp)cos(Pz)dP
(11)
nG j 2(1 - v)Kf (Pp) - P2p2 (К2 (Pp) - К2 (Pp))
При выводе формулы (11) использовано значение интеграла [35]
да 2
JsKl (Ps)ds = p-[Ki2 (Pp) - К2 (Pp)].
p
(12)
2. Контактная задача
Рассмотрим осесимметричную контактную задачу о взаимодействии без учета сил трения и шероховатости упругого неоднородного (с переменным коэффициентом Пуассона) пространства с цилиндрической полостью радиуса р с жестким вкладышем длины 2a, цилиндрическая поверхность которого описывается уравнением г=р+5 (5-натяг). Требуется определить контактное давление стг(р^) = — q(z) в области контакта |z|<a. Учитывая выражение (10) и условие контакта и(рХ)=5 (И^), после введения безразмерных величин
x = —, u = Pp, X = p, a a
Ф( x ) = ,
G
(13)
v„ = v(œX f = -
придем относительно функции ф(х) к интегральному уравнению
j ф©£ ^^ d%=nf (|x |< 1)
(14)
с ядром
p
к (г) = | К (и)соч(и£)йи, К (и) -
Р0 (и) = — | "лООК02 I — I
(1 -Уо)-1 К2(и) 2(К2(и) - 2^о (и)),
(15)
Допустим, что функция л (г) (3) разлагается в ряд по нечетным отрицательным степеням радиальной координаты вида
3 5
л(г) = — + ^2 — + ••• (р< г <да)
г г
1
(16)
2(1 -Уя )
Тогда, подставляя разложение (16) в формулу (15) для Я>(и), используя формулу (12), интеграл [35]
г л* = 2 [Р2р2 (К!2 (Рр) - к02 (рр>)+ р 5 рь
1
2
+ 1 К2(рр) -ррКо(рр)К,(рр)
(17)
и рекуррентное значение интеграла [35]
^ = _4тР2ГК^Ё£) л* +
гКр(Р*)±_ .2 |К^^+^е-!!
J е"1+2 (тА- П3 J (т +1)3
(.т +1) ; "
р2р2 К1(Рр) - .+1 Ко(рр)
Р2р2 -
(т +1)2
К2(Рр)^ (т = 2, 4,--.),
(18)
можно последовательно вычислять возникающие интегралы по переменной 5. В результате можно показать, что при удержании любого конечного числа слагаемых в разложении (16) асимптотическое поведение в бесконечности и нуле символа ядра К(и) вида (15) принципиально не меняется и имеет вид
К (и) = с0и 1 + си 2 + О (и 3) (и ^ да), К (и) = й0 + О(и) (и ^ 0).
(19)
Для закона изменения (20) значение коэффициента Пуассона в бесконечности (вдали от границы полости) не зависит от значения ль При л 1=0 имеем однородный упругий материал. При приближении к границе полости из бесконечности (от значения у») коэффициент Пуассона и модуль Юнга убывают при ^1>0 и возрастают при "Л1<0. При этом на границе полости коэффициент Пуассона достигает значения
у(р) = 1+-
0,5
(21)
С целью приведения коэффициента с0 в разложении (19) к единице домножим обе части интегрального уравнения (14) на
1-0,5^, л = 4(1- у„ И и будем рассматривать интегральное уравнение
(22)
Г Ф©к
х-X
й^ = я(1-0,5л^./" (| х |< 1), (23)
для которого символ ядра, используя интегралы (12) и (17), запишем в форме
К (и) =
1 - 0,5л
2(1 -у„ )-и2(1-2ли2)(а2-1)-ли 2(а2-2иО0)'
Й0 =
Кр(и) К(и) •
(24)
При л=0 символ (24) совпадает с известным для однородного пространства с цилиндрической шахтой [12 15].
Для символа ядра (24), используя известные асимптотические формулы для функций Макдональда [35], получим асимптотику
^ 1 1 3
= 1--+ -
3
2и 8и2 8и
+ 63 , + О| (и ^+да) 128и4 I и5 ]
и найдем коэффициенты в разложениях (19):
(25)
Поведение (19) обеспечивает возможность применения сингулярного асимптотического метода, эффективного при малых значениях параметра X [14-16] для решения интегрального уравнения (14). Безразмерный параметр X, введенный в формулах (13), характеризует относительную длину вкладыша. Малые значения X соответствуют относительно длинным вкладышам.
Далее будем удерживать в разложении (16) только два первых члена. Тогда
1 - 2у + 1,25л 1 - 0,5л _ с = 1, с =--2-,-1, й„ =-. (26)
1 - 0,5л
2(1 -У„ )
На основании асимптотического поведения символа (24) ядро (15) интегрального уравнения (23) имеет логарифмическую особенность, которую можно выделить в форме
к (г) = - 1п| г | +|
иК (и) -1 + ехр(-и)
С0$(иг)йи, (27)
у(г) = 1+
0,5
л„ +л1р3 г-3
л„ =■
1
(20)
2(1 -У„ )
ую = У (да).
где использовано значение интеграла [13,15] г1 ехр( и) С05(и/)йи = - 1п | г |.
(28)
0
2 да
р
-1
2
X
и
0
и
0
Зная главную часть ядра (27), можно утверждать, что решение интегрального уравнения первого рода (23) должно иметь корневую особенность на границе области контакта x = +1.
Ниже используем сингулярный асимптотический метод для решения интегрального уравнения (23) с символом ядра (24). Метод предполагает построение решения в виде суперпозиции вырожденного (проникающего) решения и решений типа пограничного слоя, несущих корневые особенности на краях области контакта. При малых значениях X заменим интегральное уравнение (23) на эквивалентную систему трех интегральных уравнений [13—15]:
1±_5Ы *—5 х
d 5 =
— 1 J
«I1-5!—p(5)
к| 5 |d5 (— 1 <*<x),
Н^ к (d 5=
= л:Х + J
«(1-5p(5)
к I | d5 (— x< * < 1),
(29)
P'x) к (^x5 1 d 5=«X (— ю< * <ю)
при условии
ф( *) = (1 — 0,5 л/ X X
1 + *) ( 1 — *) ( * ю| — М ~)—p U
(I* 1< 1). (30)
Из последнего интегрального уравнения (29) при помощи интегрального преобразования Фурье сразу найдем вырожденное решение, не учитывающее корневые особенности на краях области контакта (см. формулы (19), (26)),
p(*X ') = d0
(31)
Первые два интегральных уравнения (29), служащие для построения решений типа погранслоя, заменами переменных сводятся к одному интегральному уравнению вида
да да Г 2 ^
|ю(х)к(г — т)Л = л+ | [ю(х) — d01]kI - — г — х Idх . (32)
0 2/X \Х /
Можно показать, что интеграл в правой части интегрального уравнения (32) экспоненциально убывает при X—^0, для решения этого интегрального уравнения можно использовать метод последовательных приближений, на каждом шаге которого возникает интегральное уравнение Винера—Хопфа. При малых значениях X ограничимся нулевым приближением решения интегрального
уравнения (32), пренебрегая интегралом в правой его части. Для эффективной факторизации символа ядра в соответствии с его свойствами (19), (26) используем аппроксимацию
K (и) —^exP
и2 + C2 B
D
%/и2 + E2
D = ~rexP| — ! = do, E > 1.
(33)
Значения постоянных, входящих в формулу (33), и погрешность этой аппроксимации на действительной оси 6 приведены в табл. 1, рассчитанные при Уда=0,3 и разных Известно [13], что погрешность решения интегрального уравнения (23), связанная с заменой точного выражения символа ядра на приближенное (33), не превосходит погрешности аппроксимации 6.
Применяя метод Винера—Хопфа, решение типа по-гранслоя интегрального уравнения (32), (33) в нулевом приближении найдем в виде
ю(s) = ■
1 (exp(—Bs) +C erf <JBS) + / (s)
Tdo ^ V^s Vb v > )
1 (s) = —D J
TT J
с
+C erf(^ B(s — x))
exp(—B(s — x)) + s — x)
K (Ex)d x.
(34)
Таблица 1
Значения параметров аппроксимации (32)
Table 1
Parameters of the approximation (32)
л B C D E 0 (%)
—0,08 1,110 1,129 —0,288 1,810 5,5
—0,04 1,030 1,069 —0,343 1,610 5,0
0 1,070 1,081 —0,400 1,610 4,5
0,04 1,000 1,052 —0,459 1,810 4,5
0,08 1,020 1,071 —0,521 2,010 4,5
Решение (34) имеет корневые особенности на краях области контакта.
Формулы (30), (31) и (34) дают сингулярное асимптотическое решение интегрального уравнения (23). На основании этих формул получим интегральную характеристику контактных давлений
1
N = J ф( *)d*
в форме
No = (1 — 0,5л)/ ^
2C
d0B
с-2B+C) ) -
X
ю
—1
—1
X
0
X
—1
exp(- ßQ
nß
С , 2J(С) К = 2
- —+
d0 J
J (s) = — 1
n 0
С (s-т ß
exp(-B(s - т)) +
^ßp» (1-С + C(s-т)
K (ET)dт . (36)
Погрешность решения (30), (31), (35) и (36) при X<1 не превосходит (5+6)%. Это решение можно рекомендовать, когда длина вкладыша существенно больше диаметра цилиндрической полости.
В табл. 2 приведены значения интегральной характеристики No/f, рассчитанные по формуле (36) при vM=0,3 и разных значениях X и При уменьшении параметра X контактные давления и интегральная характеристика существенно возрастают.
Таблица 2
Значения интегральной характеристики (35)
Table 2
Integral characteristics (35)
2 1 0,5 0,25
л
-0,08 2,95 4,67 8,00 14,6
-0,04 2,71 4,29 7,34 13,4
0 2,47 3,90 6,67 12,2
0,04 2,20 3,49 5,97 10,9
0,08 1,94 3,07 5,25 9,61
3. Учет шероховатости
Допустим, что поверхность полости является не идеально гладкой, как в предыдущем разделе, а шероховатой. Известно [36], что шероховатость поверхности в контактных задачах может моделироваться покрытием винклеровского типа. Пусть поверхность шахты армирована таким покрытием, которое дает дополнительный вклад в нормальное перемещение в области контакта (снова используем размерные обозначения):
U = Aq( z), A = const,
(37)
где д(2) - искомое контактное давление. Условие контакта имеет вид
и,+ uq = 8
(|z|< a),
(38)
где и - упругое перемещение точек поверхности под действием контактного давления, которое определяется с использованием фундаментального решения (10); 5 -натяг. Для изменяющегося коэффициента Пуассона снова принимаем закон (20).
Используя условие (38), безразмерные величины (13), а также обозначение
AG
A = ■
(39)
(1- Уда )а
сведем контактную задачу для шахты с покрытием к интегральному уравнению относительно безразмерного контактного давления ф(х):
A.c.<p(X) +11 <(%)k (■j d% = cf (| X |< 1),
да
k(t) = 1K (u)cos(ut)du, c, = 1-0,5r|.
(40)
Здесь символ ядра К(и) определяется по формуле (24).
В отличие от уравнения (23) интегральное уравнение (40) является уравнением второго рода. Несмотря на логарифмическую особенность ядра, его решение уже не обладает корневыми особенностями и принимает конечные значения на краях области контакта. Для численного решения интегрального уравнения (40) используем метод коллокации на основе квадратурной формулы Гаусса по 32 узлам [34]. Метод эффективен при не слишком малых значениях параметра X, когда для расчета ядра можно применять формулу (27). При расчетах брали X>l.
В табл. 3 приведены значения интегральной характеристики Ы0//, рассчитанные по формуле (35) при Ую=0,3 и разных значениях Д, X и л.
Таблица 3
Значения характеристики (34) по методу коллокации
Table 3
Characteristics (34) in collocation method
4 2 1 4 2 1
л A = 0,05 A = 0,01
-0,08 1,68 2,31 3,34 1,76 2,46 3,64
-0,04 1,67 2,30 3,33 1,75 2,45 3,63
0 1,66 2,30 3,32 1,74 2,44 3,62
0,04 1,65 2,29 3,31 1,73 2,43 3,60
0,08 1,64 2,28 3,30 1,72 2,42 3,59
Заключение
При необходимости аппроксимация (33) может быть уточнена путем добавления множителей по типу аппроксимации Паде [16].
Увеличение контактных давлений и интегральной характеристики при уменьшении параметра X для гладкого и шероховатого контакта (см. табл. 2 и 3) связано с возрастанием площади области контакта (при фиксированном радиусе шахты р):
0
_ 4л
2 = X '
(41)
P
Как показывает численный анализ (см. табл. 2 и 3), при л<0, когда коэффициент Пуассона возрастает при приближении к поверхности шахты от значения в бесконечности у» (при этом модуль Юнга также возрастает) контактные давления и их интегральная характеристика, отнесенные к натягу вкладыша /, больше, чем при убывающем в направлении из бесконечности коэффициенте Пуассона и модуле Юнга (л>0).
При малых значениях параметра Д , характеризующего коэффициент постели покрытия, решение интегрального уравнения (40) вдали от краев области контакта становится близко к вырожденному (проникающему) решению (31). Решение уравнения (40) возрастает вблизи краев области контакта тем больше, чем меньше значение Д , что согласуется с эффектом
регуляризации интегрального уравнения (23) по А.Н. Тихонову. При увеличении параметра A (возрастании степени шероховатости) контактные давления и их интегральная характеристика снижаются (см. табл. 3). Сравнение табл. 2 и 3 показывает, что покрытие существенно снижает влияние неоднородности на интегральную характеристику и распределение контактных давлений.
Благодарности
Работа выполнена при финансовой поддержке РФФИ (грант №18-01-00017).
Acknowledgements
The work was carried out with the financial support from the Russian Foundation for Basic Research (Project Nr. 18-01-00017).
Библиографический список
1. Axisymmetric contact problem of the theory of elasticity for inhomogeneous layers / A.S. Vasiliev, S.S. Volkov, S.M. Aizikovich, Y.-R. Jeng // ZAMM. - 2014. - Vol. 94. - No. 9 - P. 705712. DOI: 10.1002/zamm.201300067
2. Айзикович С.М., Васильев А.С., Волков С.С. Осесим-метричная контактная задача о вдавливании конического штампа в полупространство с неоднородным по глубине покрытием // Прикладная математика и механика. - 2015. -Т. 79, вып. 5. - С. 710-716.
3. Axisymmetric problem on the indentation of a hot circular punch into an arbitrarily non-homogeneous half-space / L.I. Kre-nev, S.M. Aizikovich, Y.V. Tokovyy, Y.-C. Wang // International Journal of Solids and Structures. - 2015. - Vol. 59 - P. 18-28. DOI: 10.1016/j.ijsolstr.2014.12.017
4. Напряженно-деформированное состояние упругого мягкого функционально-градиентного покрытия при внедрении сферического индентора / С.С. Волков, А.С. Васильев, С.М. Айзикович, Н.М. Селезнев, А.В. Леонтьева // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 4. - С. 20-34. DOI: 10.15593/perm.mech/2016.4.02
5. Torsion of a circular punch attached to an elastic half-space with a coating with periodically depth-varying elastic properties / A.S. Vasiliev, M.V. Swain, S.M. Aizikovich, E.V. Sadyrin // Archive of Applied Mechanics. - 2016. - Vol. 86. - No. 7 - P. 12471254. DOI: 10.1007/s00419-015-1089-1
6. Influence of a soft FGM interlayer on a contact stresses under a beam on an elastic foundation / S.M. Aizikovich, B.I. Mitrin, N.M. Seleznev, Y.-C. Wang, S.S. Volkov // Structural Engineering and Mechanics. - 2016. - Vol. 58. - No. 4 - P. 613-625. DOI: 10.12989/sem.2016.58.4.613
7. 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 -P. 11-15.
8. Индикация термоупругой неустойчивости скользящего контакта с помощью заглубленной пьезокерамической прослойки / В.Б. Зеленцов, Б.И. Митрин, А.Г. Сукиязов,
С.М. Айзикович // Вестник Пермского национального исследовательского политехнического университета. Механика. -2017 - № 1. - С. 63-84. DOI: 10.15593/perm.mech/2017.1.05
9. Indentation of functionally graded coating on an elastic substrate by a sphero-conical indenter / L.I. Krenev, E.V. Sadyrin, S.M. Aizikovich, T.I. Zubar // Springer Proceedings in Phycics. -2017. - Vol. 193 - P. 397-405. DOI: 10.1007/978-3-319-56062-5_33
10. Vasiliev A.S., Volkov S.S., Aizikovich S.M. Approximated analytical solution of contact problem on indentation of elastic half-space with coating reinforced with inhomogeneous interlayer // Materials Physics and Mechanics. - 2018. - Vol. 35. - No. 1 -P. 175-180. DOI: 10.18720/MPM.3512018_20
11. Aizikovich S.M., Galybin A.N., Krenev L.I. Semi-analytical solution for mode I penny-shaped crack in a soft inhomoge-neous layer // International Journal of Solids and Structures. -2015. - Vol. 53 - P. 129-137. DOI: 10.1016/j.ijsolstr.2014.10.010
12. Александров В.М., Белоконь А.В. Асимптотическое решение одного класса интегральных уравнений и его применение к контактным задачам для цилиндрических упругих тел // Прикладная математика и механика. - 1967. - Т. 31, вып. 4. -С. 704-710.
13. Александров В.М., Ромалис Б.Л. Контактные задачи в машиностроении. - М.: Машиностроение, 1986. - 176 с.
14. Александров В.М., Пожарский Д.А. Об одном асимптотическом методе в контактных задачах // Прикладная математика и механика. - 1999. - Т. 63. - Вып. 2. - С. 295-302.
15. Alexandrov V.M., Pozharskii D.A. Three-dimensional contact problems. - Dordrecht: Kluwer, 2001. - 406 p.
16. Пожарский Д.А. Контактная задача для полого цилиндра // Прикладная математика и механика. - 2017. - Т. 81. -Вып. 6. - С. 727-733.
17. Пожарский Д.А., Золотов Н.Б. К одной задаче Белоко-ня А.В. // Вестн. Дон. гос. техн. ун-та. - 2017. - Т. 17, №. 2. -С. 7-11. DOI: 10.23947/1992-5980-2017-17-2-7-11
18. Калинчук В.В., Белянкова Т.И. Динамика поверхности неоднородных сред. - М.: Физматлит, 2009. - 316 с.
19. Калинчук В.В., Белянкова Т.И. Динамическая контактная задача для заполненной жидкостью преднапряженной цилиндрической трубы // Прикладная математика и механика. - 2009. - Т. 73. - Вып. 2. - С. 289-302.
20. Абрамович М.В., Колосова Е.М., Чебаков М.И. Контактная задача при наличии сил трения в зоне контакта для трехкомпонентного цилиндрического основания // Прикладная математика и механика. - 2014. - Т. 78. - Вып. 2. -С. 262-269.
21. Белянкова Т.И., Калинчук В.В., Лыжов В.А. Особенности динамики трехслойного полого цилиндра // Экологический вестник научных центров ЧЭС. - 2015. - № 4. - С. 19-32.
22. Finite-element modeling of a damaged pipeline repaired using the wrap of a composite material / A.A. Lyapin, M.I. Che-bakov, A. Dumitrescu, G. Zecheru // Mechanics of Composite Materials. - 2015. - Vol. 51. - No. 3 - P. 333-340. DOI: 10.1007/s11029-015-9504-9
23. Ватульян А.О., Юров В.О. Волновые процессы в полом цилиндре в поле неоднородных предварительных напряжений // Прикладная механика и техническая физика. - 2016. -Т. 57, № 4. - С. 182-191. DOI: 10.15372/PMTF20160418
24. Identification of inhomogeneous elastic properties of isotropic cylinder / I.V. Bogachev, R.D. Nedin, A.O. Vatulyan, O.V. Yavrunyan // ZAMM. - 2017. - Vol. 97. - No. 3 - P. 358364. DOI: 10.1002/zamm.201600179
25. 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
26. Goryacheva I.G., Tsukanov I.Y. Modeling of normal contact of elastic bodies with surface relief taken into account // Jour-
References
1. Vasiliev A.S., volkov S.S., Aizikovich S.M., Jeng Y.-R. Axisymmetric contact problem of the theory of elasticity for in-homogeneous layers. ZAMM, 2014, vol. 94, no. 9, pp. 705-712. DOI: 10.1002/zamm.201300067
2. Aizikovich S.M., Vasil'ev A.S., volkov S.S. The axisym-metric contact problem of the indentation of a conical punch into a half-space with coating inhomogeneous in depth. Journal of Applied Mathematics and Mechanics, 2015, vol. 79, no. 5, pp. 500505. DOI: 10.1016/j.jappmathmech.2016.03.011
3. Krenev L.I., Aizikovich S.M., Tokovyy Y.V., Wang Y.-C. Axisymmetric problem on the indentation of a hot circular punch into an arbitrarily non-homogeneous half-space. International Journal of Solids and Structures, 2015, vol. 59, pp. 18-28. DOI: 10.1016/j.ij solstr.2014.12.017
4. 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]. PNRPUMechanics Bulletin, 2016, no. 4, pp. 2034. DOI: 10.15593/perm.mech/2016.4.02
5. Vasiliev A.S., Swain M.V., Aizikovich S.M., Sadyrin E.V. Torsion of a circular punch attached to an elastic half-space with a coating with periodically depth-varying elastic properties. Archive of Applied Mechanics, 2016, vol. 86, no. 7, pp. 1247-1254. DOI: 10.1007/s00419-015-1089-1
6. Aizikovich S.M., Mitrin B.I., Seleznev N.M., Wang Y.-C., volkov S.S. Influence of a soft FGM interlayer on a contact stresses under a beam on an elastic foundation. Structural Engineering and Mechanics, 2016, vol. 58, no. 4, pp. 613-625. DOI: 10.12989/sem.2016.58.4.613
7. Vasiliev A.S., volkov S.S., Aizikovich S.M. Indentation of an axisymmetric punch into an elastic transversely-isotropic half-
nal of Physics: Conference Series. - 2018. - Vol. 991. - No. 1 -P. 012028. DOI: 10.1088/1742-6596/991/1/012028
27. Бородачев А.Н., Дудинский В.И. Жесткий штамп на упругом полупространстве с изменяющимся по глубине коэффициентом Пуассона // Прикладная механика. 1985. - Т. 21, № 8. - С. 34-39.
28. Бородачев А.Н. Упругое равновесие неоднородного по толщине слоя // Прикладная механика. 1988. - Т. 24, № 8. -С. 30-35.
29. Пожарский Д.А. Упругое равновесие неоднородного клина с переменным коэффициентом Пуассона // Прикладная математика и механика. - 2016. - Т. 80. - Вып. 5. - С. 614-621.
30. Кузнецов Е.А. Давление круглого цилиндра на полупространство с переменным по глубине коэффициентом Пуассона // Изв. АН СССР. МТТ. - 1985. - № 1. - С. 73-86.
31. Ломакин В.А. Теория упругости неоднородных тел. -М.: Ленанд, 2014. - 376 с.
32. Gurtin M.E. The linear theory of elasticity. Handbuch der Physik. vol. VIa/2. - Berlin: Springer, 1972 - P. 1-295.
33. Лурье А.И. Пространственные задачи теории упругости. - М.: ГИТТЛ, 1955. - 491 с.
34. Справочник по специальным функциям / под ред. М. Абрамовица и И. Стиган. - М.: Наука, 1979. - 832 с.
35. Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. Специальные функции. - М.: Наука, 1983. - 752 с.
36. Александров В.М., Мхитарян С.М. Контактные задачи для тел с тонкими покрытиями и прослойками. - М.: Наука, 1983. - 487 с.
space with functionally graded transversely-isotropic coating. Materials Physics and Mechanics, 2016, vol. 28, no. 1-2, pp. 11-15.
8. Zelentsov V.B., Mitrin B.I., Sukiyazov A.G., Aizikovich S.M. 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
9. Krenev L.I., Sadyrin E.V., Aizikovich S.M., Zubar T.I. Indentation of functionally graded coating on an elastic substrate by a sphero-conical indenter. Springer Proceedings in Phycics, 2017, vol. 193, pp. 397-405. DOI: 10.1007/978-3-319-56062-5_33
10. Vasiliev A.S., volkov S.S., Aizikovich S.M. Approximated analytical solution of contact problem on indentation of elastic half-space with coating reinforced with inhomogeneous interlayer. Materials Physics and Mechanics, 2018, vol. 35, no. 1, pp. 175180. DOI: 10.18720/MPM. 3512018_20
11. Aizikovich S.M., Galybin A.N., Krenev L.I. Semi-analytical solution for mode I penny-shaped crack in a soft inho-mogeneous layer. International Journal of Solids and Structures, 2015, vol. 53, pp. 129-137. DOI: 10.1016/j.ijsolstr.2014.10.010
12. Aleksandrov V.M., Belokon' A.V. Asymptotic solution of a class of integral equations and its application to contact problems for cylindrical elastic bodies. Journal of Applied Mathematics and Mechanics, 1967, vol. 31, no. 4, pp. 718-724.
13. Aleksandrov V.M., Romalis B.L. Kontaktnye zadachi v mashinostroenii [Contact problems in machine-building]. Moscow, Mashinostroenie, 1986, 176 p.
14. Aleksandrov V.M., Pozharskii D.A. An asymptotic method in contact problems. Journal of Applied Mathematics and Mechanics, 1999, vol. 63, no. 2, pp. 283-290.
15. Alexandrov V.M., Pozharskii D.A. Three-dimensional contact problems. Dordrecht, Kluwer, 2001, 406 p.
16. Pozharskii D.A. Contact problem for a hollow cylinder. Journal of Applied Mathematics and Mechanics, 2017, vol. 81, no. 6, pp. 499-503. DOI: 10.1016/j.jappmathmech. 2018.03.020
17. Pozharskii D.A., Zolotov N.B. K odnoi zadache Beloko-nia A.V. [To one Belokon's problem]. Vestnik Donskogo gosudar-stvennogo tekhnicheskogo universiteta - DSTUBulletin, 2017, vol. 17, no. 2, pp. 7-11. DOI: 10.23947/1992-5980-2017-17-2-7-11
18. Kalinchuk V.V., Belyankova T.I. Dinamika poverkhnosti neodnorodnykh sred [Surface dynamics of inhomogeneous media]. Moscow, Fizmatlit, 2009, 316 p.
19. Belyankova T.I., Kalinchuk V.V. The dynamic contact problem for a prestressed cylindrical tube filled with a fluid. Journal of Applied Mathematics and Mechanics, 2009, vol. 73, no. 2, pp. 209-219. DOI: 10.1016/j.jappmathmech.2009.04.011
20. Abramovich M.V., Kolosova Ye.M., Chebakov M.I. The contact problem when there are friction forces in the contact area for a three-component cylindrical base. Journal of Applied Mathematics and Mechanics, 2014, vol. 78, no. 2, pp. 181-186. DOI: 10.1016/j.jappmathmech.2014.07.011
21. Belyankova T.I., Kalinchuk V.V., Lyzhov V.A. Osoben-nosti dinamiki trekhsloinogo pologo tsilindra [Peculiarities of dynamics of a three-layered hollow cylinder]. Ecological Bulletin of BEC Scientific Centers, 2015, no. 4, pp. 19-32.
22. Lyapin A.A., Chebakov M.I., Dumitrescu A., Zecheru G. Finite-element modeling of a damaged pipeline repaired using the wrap of a composite material. Mechanics of Composite Materials, 2015, vol. 51, no. 3, pp. 333-340. DOI: 10.1007/s11029-015-9504-9
23. Vatulyan A.O., Yurov V.O. Wave processes in a hollow cylinder in an inhomogeneous prestress field. Journal of Applied Mechanics and Technical Physics, 2016, vol. 57, no. 4, pp. 731 -739. DOI: 10.1134/S0021894416040180
24. Bogachev I.V., Nedin R.D., Vatulyan A.O., Yavrunyan O.V. Identification of inhomogeneous elastic properties of iso-tropic cylinder. ZAMM, 2017, vol. 97, no. 3, pp. 358-364. DOI: 10.1002/zamm.201600179
25. 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
26. 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
27. Borodachev A.N., Dudinskii V.I. Rigid punch on an elastic semispace with a depth-varying Poisson ratio. Soviet Applied Mechanics, 1985, vol. 21, no. 8, pp. 753-757.
28. Borodachev A.N. Elastic equilibrium in a layer inhomogeneous with depth. Soviet Applied Mechanics, 1988, vol. 24, no. 8, pp. 753-758.
29. Pozharskii D.A. The elastic equilibrium of an inhomogeneous wedge with varying Poisson's ratio. Journal of Applied Mathematics and Mechanics, 2016, vol. 80, no. 5, pp. 433-438. DOI: 10.1016/j.jappmathmech.2017.02.010
30. Kuznetsov E.A. Pressure of a curcular cylinder against a half-space with Poisson's ratio that varies with depth. Mechanics of Solids, 1985, vol. 20, no. 1, pp. 68-80.
31. Lomakin V.A. Teoriia uprugosti neodnorodnykh tel [Theory of elasticity of inhomogeneous bodies]. Moscow, Lenand, 2014, 376 p.
32. Gurtin M.E. The linear theory of elasticity. Handbuch der Physik. vol. VIa/2. Berlin, Springer, 1972, pp. 1-295.
33. Lur'e A.I. Prostranstvennye zadachi teorii uprugosti [Spatial problems in the elasticity theory]. Moscow, GITTL, 1955, 491 p.
34. Handbook of Mathematical Functions. Edited by M. Abra-mowitz and I. Stigan. New York, Dover Publications, 1964, 832 p.
35. Prudnikov A.P., Brychkov Yu.A., Marichev O.I. Integraly i riady. Spetsial'nye funktsii [Integrals and series. Special functions]. Moscow, Nauka, 1983, 752 p.
36. Aleksandrov V.M., Mkhitaryan S.M. Contact problems for bodies with thin coatings and interlayers [Kontaktnye zadachi dlia tel s tonkimi pokrytiiama i prosloikami]. Moscow, Nauka, 1983, 487 p.
