CC BY
127
МАТЕМАТИЧЕСКИЕ И ЕСТЕСТВЕННЫЕ НАУКИ
УДК 539.3 й01 10.12737/4547
Осесимметричный изгиб круглой многослойной пластины на упругом основании
V *
сложной структуры
С. М. Айзикович, С. С. Волков, А. В. Мелконян
Получено в аналитическом виде приближённое решение задачи об изгибе круглой многослойной пластины постоянной толщиныI, лежащей на упругом основании сложной структуры/. Пластина изгибается под действием осесимметричной распределённой нагрузки и реакции со стороны/ основания. Упругое основание представляет собой непрерывно-неоднородный по толщине слой (покрытие), лежащий на однородном полупространстве (подложке). Модуль Юнга в зоне сопряжения покрытия и подложки имеет существенный скачок. Для пластины/ рассмотрены/ два случая граничных условий: условия закреплённого и свободного края. Построенное приближённое аналитическое решение задачи эффективно в широком диапазоне как геометрических параметров (толщина неоднородного слоя и радиус пластины), так и физических параметров (гибкость пластины/ и упругие свойства покрытия и подложки). Методом интегральных преобразований контактная задача сводится к решению системы! интегро-дифференциальных уравнений. Полученные формулы/ могут быть использованы/ для расчёта характеристик контактного взаимодействия многослойной пластины с основанием сложной структуры в зависимости от граничных условий и характера нагрузки на пластину.
Ключевые слова: неоднородные материалы, многослойная пластина, функционально-градиентное покрытие, осесимметричная задача, аналитические методы, приближённое аналитическое решение.
Введение. Применение функционально-градиентных материалов существенно влияет на все характеристики контактного взаимодействия [1-3]. В случае, когда во взаимодействие вовлечены тонкие гибкие элементы (пластины), на перераспределение контактных давлений в зоне контакта влияет как неоднородность свойств взаимодействующих элементов, так и жёсткость пластины. Это нужно учитывать при расчёте взаимодействия гибких элементов (пластин) с неоднородными структурами, что приводит к необходимости рассмотрения контактной задачи о взаимодействии пластины и неоднородного основания.
Простейшая осесимметричная задача контактного взаимодействия как для слоистого, так и непрерывно-неоднородного покрытия упругого полупространства рассмотрена в работе [4]. Задача об изгибе пластины на упругом изотропном и однородном основании рассматривалась в работах [5, 6]. Решение строилось путём представления контактных напряжений в виде степенного ряда, с последующим определением коэффициентов разложения из бесконечной алгебраической системы уравнений.
Методом ортогональных многочленов такая задача решалась в работах [7, 8], а методом коллокации по чебышёвским узлам — в работах [9, 10]. При этом возникала необходимость построить решение некоторых бесконечных систем линейных алгебраических уравнений и ставилась проблема исследования сходимости полученного решения к точному. В работах [11, 12] для решения задачи применялись асимптотические методы типа «больших Л» и специальных ортогональных многочленов, что позволило получить основные характеристики решения в нескольких формах, каждая из которых эффективна в своей области изменения характерных параметров задачи.
* Результаты работы получены при выполнении проекта, поддержанного грантом РФФИ № 13-08-90916-мол_ин_нр.
5
Отметим, что большинство известных решений эффективны только для жёстких пластин. И очень немногие, в частности, представленные в [11, 12], эффективны или для гибких, или для жёстких пластин, каждое в своей области.
Интерес к решению задачи и её актуальность сохраняется и в настоящее время. Так, в работе [13] решение строилось с использованием разложения напряжения в двойной ряд Фурье. Аналогичный подход использовался в работе [14]. Andrea R. D. Silva с соавторами развил численные методы решения задачи [15]. В работе [16] было получено приближённое аналитическое решение осесимметричной задачи об изгибе пластины со свободным краем, лежащей на функционально-градиентном основании. Для этого использован двусторонне асимптотический метод [17].
В настоящей работе метод [17] использован для решения задачи об изгибе многослойной пластины.
Постановка задачи. Круглая многослойная пластина радиуса R и постоянной толщины h лежит на поверхности z = 0 упругого полупространства, состоящего из неоднородного мягкого слоя (покрытия) толщины H (-H < z < 0) и однородного полупространства (подложки), причём h¡ — постоянная толщина каждого слоя пластины. С полупространством связана цилиндрическая система координат r,p,z; координата готнесена к радиусу пластины R, координата z— к толщине покрытия H, ось z нормальна поверхности и проходит через центр пластины. Под действием осесимметричной распределённой нагрузки p* (r) и реакции со стороны слоя пластина изгибается. Прогиб
пластины обозначается функцией w* (r).
Коэффициенты Ламе мягкого слоя меняются с глубиной по закону:
, ч ÍA,(z), -1 < z < 0,
Л (z)4
v ; |Л2, -»< z <-1,
, ч ÍM.(z), -1 < z < 0, M (z) = -j
v ' |M2, -»<z <-1.
При расчётах для задания характера неоднородности покрытия удобно использовать модуль Юнга, а коэффициент Пуассона зафиксировать и считать постоянным. Известно [18], что
e = м3Л±2М, V = . Л
Л + М ' 2 (Л + М)'
В качестве параметра, характеризующего отличие упругих свойств слоя от подложки, введём величину
где Е1 и Е2 — модули Юнга покрытия и подложки соответственно. Уравнение изгиба многослойной пластины имеет вид:
Ц^ (г) = р' (г)-^2 (г),
, (г) = (г)-Я,--1 (г), 0 < г < 1, (2)
(г) = як-1- (г)-а*(г),
где Ц = Di
С 1 С2
с1г2 + гсСг
— дифференциальный оператор изгиба /-й пластины; q* (г) — контактные
напряжения под пластиной; w/ (г) — прогиб каждого слоя; О, — цилиндрическая жёсткость /-й пластины.
Для напряжений, возникающих между слоями пластины, выполнено следующие условие сопряжения:
q/_1 (r)=q/ (г), 0 < г < 1. (3)
Если все слои жёстко сцеплены друг с другом, то w1 (г) = w2 (г) =... = шк (г) = w* (г). Учитывая это обстоятельство, а также равенство (3), система уравнений (2) может быть представлена в виде:
(г) = р (г)-д (г), где О = £ О,, Ц, =
( /-/2
с 1 с
сСг2 + гсСг
, 0 < г < 1.
(4)
Таким образом, вместо многослойной платины можно рассматривать однослойную пластину с цилиндрической жёсткостью, равной сумме цилиндрических жёсткостей каждого слоя в многослойной пластине.
Если пластина неоднородна по толщине, её модуль Юнга Епл (z) при этом изменяется по толщине пластины, а коэффициент Пуассона vпл постоянный, то, согласно [19]:
1 ь
О = ^ ! ^ ^ .
(5)
Введём следующую замену переменных: Л = Н/Я , w*(r) = w (г')Я., р'(г) = р(г')ОЯ3, д* (г) = д (г')ОЛ3. Знак штриха далее опускаем.
Рассматривается два случая граничных условий: первый, когда края пластины закреплены, в этом случае функция прогиба пластины и/(г) удовлетворяет условиям:
дw
~дг
= 0, — (Аw)
дг '
= 0,
и второй случай, когда имеют место условия свободного края:
Г д^ Упл дWЛ
дг2 г дг
= 0, —(Аw) дг
= 0,
(6)
(7)
где vпл — коэффициент Пуассона пластины; А — оператор Лапласа.
Методом интегральных преобразований контактная задача сводится к решению системы уравнений:
Ц^г) = р(г)-д(г), 0 < г < 1, (8)
(9)
|£(а)!(аЛ(аг)Са = sw(г), 0 < г < 1,
0
да
|С(а)Л (аг)Са = 0, г > 1,
, 0
где 1-(и) — трансформанта ядра; 5 = 0Л3О-1 — параметр, характеризующий изгибную жёсткость
пластины.
г=1
г=1
г=1
г=1
Построение решения. Для однородной пластинки со свободными краем, лежащей на неоднородном полупространстве, в работе [16] построено решение задачи в аналитическом виде. Учитывая (4), (5), это решение можно использовать и для неоднородных пластин.
Построим решение задачи в аналитическом виде, в случае, когда края пластины закреплены. Для этого представим функцию прогибов в виде ряда по формам собственных колебаний круглой пластины с закреплённым краем, аналогично работе [20]:
» 1 ™ (0=Х™тФт = (р) Фт (Р) Р^ (10)
т=0 о
где Фт (г) = ^^ (кг); к — корни уравнения Зх (к) = 0, а Ат = 72 /Зй (кт) для т е N .
Учитывая линейность задачи, разложим функцию контактных напряжений с(г) в ряд следующего вида:
да
ЯГ)=^тЯт (Г), 0 < Г < 1. (11)
т=0
Здесь функции цт (г) находятся из решения парного интегрального уравнения (9) с заданной правой частью вида (10).
Для построения решения рассмотрим следующее парное интегральное уравнение:
да ^ да Л
[С(а)г(аЛ) Л (аг)^а = 0611 + (иг) 1, 0 < г < 1,
0 ^ к=: 7 (12)
да
[С(а)^0 (аг)^а = 0, г >:
0
Из результатов работы [21] следует, что при построении решения задачи двухсторонне асимптотическим методом в качестве аппроксимации главной части трансформанты ядра можно использовать аналитическое выражение вида
N и2 + а2
N (и)-пи;^, (13)
/=1 и + О/
где Ьк, мк, 6 — некоторые константы.
В работе [22] получено аналитическое решение парного интегрального уравнения (12) в
виде:
Т (г) = 200[|0)6 {г1 (0)^-= + |с,Ф (г,а/Л-1 ) + £ и;Ф (г,/ц] )Г (Лц; )}. (14)
Здесь введено обозначение:
сЬМ Л sh Atdt
г
Сумма в (14), содержащая функцию Ф (г,/'м), соответствует решению уравнения (12) с правой частью в виде ряда по функциям Бесселя. Функцию фт (г) можно представить в следующем виде:
Фт (г)=и?з0 (мт г),
где ит = Ат; мт = кт.
Таким образом, решение парного интегрального уравнения (9) с правой частью (10) имеет
вид:
Ф(г,А) = "1А -А[ ГАШ' .
1 ' л/Т^Т2 г Л^^-Т2
г- N
qN(0=2V2n-s LN(0)(1 -rp +£c>(ra-1)
i=0 _
, m = 1,2,
д"т (г) = 2П-1 т[£ (Лкт)Ф(г,/кт) +£С;"Ф(г,а/Л1)
/=0 _
Здесь постоянные Ст определяются из системы линейных алгебраических уравнений:
¿С,0с (3;Л-1АЛ-1) + т (0) ЛЬ- = 0, к = 1,2,. ,т,
/=1
¿С,то (а,Л1,ЬкЛ1) + в (ктЛГ1 ) = 0, к = 1,2,. ,т; т = 1,2,.,
/=1
где
а (а, Ь) = [а sh а + Ь сЬ| а_ (Ь2 - а2)-1, в (а, Ь) = # (Ла) а (/а, Ь). Полученные значения контактных напряжений д^ (г) и заданную внешнюю нагрузку, при-
ложенную к пластине, можно представить в виде следующих рядов:
да А
qN (r) = Ту?Ф, (r), y,m = к (p)ф, (p)pdp,
j=0
Р(Г)=Т РтФm (r\ Pm =J P (P) Фт (P) PdP.
(15)
(16)
Подставив разложения (10), (15), (16) в (8), придём к бесконечной системе линейных алгебраических уравнений для определения коэффициентов т т , которую можно записать в следующем виде:
Wm -akmYwjEm = PmC, т = 0,1,2,.; a = -1,
;=о
N
Em = 2nlAjSAm [L (Akj)x(kj,km)+TCmX(anA1,km)
n=0
где j = 1,2,.; m = 0,1,2,.;
x (a, b) = b 1 ^cos a sin b + a (sin (a - b)(a - b)-1 - sin (a + b)(a + b)-1 )j;
X (a,b) = (a2 + b2) 1 (ashacosb + bsins cha). В частности,
(17)
Em = 2n-1V2sAm
E° = n-14s
N
LN (0)km1sinkm +ТСПХ(anA1,km)
n=1
N
L1 (0) + A TCa^sh (anA1)
Ej = 2n-^V2 A,s
N
L (Akj )k,-1 sink, + A ^j^sh (anA1)
Используя метод редукции, сведём решение бесконечной системы алгебраических уравнений (17) к решению системы линейных алгебраических уравнений следующего вида:
Wmk +Tw,Em = Pm, т = 0,1,2,. ,M.
j=0
После определения коэффициентов wm (m = 0,1,...,M) для фиксированного значения M из (18) находим контактные напряжения
м
Q{r)=YwmQm (r), 0 < r < 1, (19)
m=0
и прогибы пластины
м
W (r)=^WmPm (r), 0 < r < 1. (20)
m=0
Полученные формулы могут быть использованы для расчёта характеристик контактного взаимодействия многослойной пластины с основанием сложной структуры в зависимости от граничных условий и характера нагрузки на пластину. Так как неоднородность пластины и многослой-ность в рамках рассмотренной модели влияет только на величину цилиндрической жёсткости D, будем считать, что увеличение количества слоёв в пластине влияет на значение параметра D и, соответственно, на значение параметра s.
Рассмотрим случай влияния на распределение контактных напряжений под пластиной равномерно распределённой нагрузки. В случае, когда края пластинки не закреплены по периметру, при вычислении контактных напряжений использовался результат работы [16]. Для закреплённых краёв использовалась формула (19).
Пусть пластина лежит на мягком покрытии (в = 100), модуль Юнга которого меняется по
закону:
E1 (z)= E0 • ф (z) , ф (z)= Ф0 +(Ф0 -1)-z . Показатель неоднородности ф0 = const > 1 характеризует отношение модуля Юнга на поверхности покрытия (z = 0) к модулю Юнга подложки. Ограничимся рассмотрением случая ф0 = 3,5, что соответствует сочетанию мягкого (алюминий, серебро, медь, свинец и т. д.) и жёсткого (железо, сталь, палладий, молибден и т. д.) металлов.
При численной реализации было выявлено, что граничные условия для жёстких пластин практически не оказывают влияние на распределение контактных напряжений. Причём диапазоны параметров, для которых это справедливо, следующие: для Л = 4 — s < 0,1, для Л = 0,05 — s < 0,07.
На рис. 1, 2 построены графики контактных напряжений в случае гибких пластин для указанного выше закона изменения модуля Юнга с глубиной. Из графиков видно, что условия закрепления краёв пластины оказывают влияние на распределение контактных напряжений под пластиной.
Заключение. Получено приближённое аналитическое решение задачи об изгибе круглой многослойной пластины постоянной толщины, лежащей на упругом неоднородном основании сложной структуры.
q(r)
О 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 г1
| Закрепленныекрая — • — Свободные края |
Рис. 1. Распределение контактных напряжений под гибкой (5 = 4) пластиной при большой зоне контакта относительно толщины покрытия
Закрепленные края — ■ — Свободные края
Рис. 2. Распределение контактных напряжений под гибкой пластиной при малой зоне контакта относительно толщины покрытия
Построенная математическая модель эффективна в широком диапазоне значений как геометрического (Л), так и физических (s, в) параметров задачи и может быть использована для расчётов неоднородных пластин Кирхгофа при различных граничных условиях на краю пластинки. Библиографический список
1. Hager, A. M. Short-Fibre Reinforced, High-Temperature Résistant Polymers for a Wide Field of Tribological Applications / A. M. Hager, M. Davies // Advances in Composite Tribology / K. Friedrich, ed. — Amsterdam : Elsevier, 1993. — Pp. 107-157.
2. Friedrich, K. Wear of polymer composites / K. Friedrich, R. Reinicke, Z. Zhang // Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology. — 2002. — Vol. 216, iss. 6. — Pp. 415-426.
3. Моделирование фрикционного взаимодействия композиционных покрытий триботехнического назначения / И. Г. Горячева [и др.] // Трение и износ. — 2012. — Т. 33, № 6. — С. 557-565.
4. Васильев, А. С. Кручение упругого полупространства с многослойным покрытием периодической структуры / А. С. Васильев, Е. В. Садырин, М. Е. Васильева // Вестник Дон. гос. техн. ун-та. — 2013. — № 5/6 (74). — С. 6-13.
5. Горбунов-Посадов, М. И. Расчёт балок и плит на упругом полупространстве // Прикладная математика и механика. — 1940. — Т. 4, вып. 3. — С. 61-80.
6. Ишкова, А. Г. Об изгибе полосы и круглой пластины, лежащих на упругом полупространстве // Инженерный сборник. — 1960. — Т. 23. — С. 171-181.
7. Гребенщиков, В. Н. Расчёт круглой пластинки на упругом полупространстве // Теория расчёта и надёжность приборов : сб. тр. II Саратовской обл. конф. молодых учёных. — 1969. — С. 48-51.
8. Александров, В. М. Универсальная программа расчёта изгиба балочных плит на линейно-деформируемом основании / В. М. Александров, Л. С. Шацких // Труды 7-й Всесоюзной конференции по теории оболочек и пластин. — Москва, 1970. — С. 46-51.
9. Шацких, Л. С. К расчёту изгиба плиты на упругом слое // Известия Академии наук СССР. Механика твёрдого тела. — 1972. — № 2. — С. 170-176.
10. Александров, В. М. Эффективное решение задачи о цилиндрическом изгибе пластинки конечной ширины на упругом полупространстве / В. М. Александров, И. И. Ворович, М. Д. Солодовник // Известия Академии наук СССР. Механика твёрдого тела. — 1973. — № 4. — С. 129-138.
11. Александров, В. М. Асимптотическое решение задачи о цилиндрическом изгибе пластинки конечной ширины на упругом полупространстве / В. М. Александров, М. Д. Солодовник // Прикладная механика. — 1974. — Т. 10, вып. 7. — С. 77-83.
12. Босаков, С. В. К решению контактной задачи для круглой пластинки / С. В. Босаков // Прикладная математика и механика. — 2008. — Т. 72, № 1. — С. 59-61.
13. Kashtalyan, M. Effect of a functionally graded interlayer on three-dimensional elastic deformation of coated plates subjected to transverse loading / M. Kashtalyan, M. Menshykova // Composite Structures. — 2009. — Vol. 89, № 2. — Pp. 167-176.
14. Kashtalyan, M. Three-dimensional elasticity solution for bending of functionally graded rectangular plates / M. Kashtalyan // European Journal of Mechanics A/Solids. — 2004. — Vol. 23, № 5. — Pp. 853-864.
15. Silva, A. R. D. Numerical methods for analysis of plates on tensionless elastic foundations / A. R. D. Silva, R. A. M. Silveira, P. B. GoncRalves // International Journal of Solids and Structures. — 2001. — Vol. 38, № 10-13. — Pp. 2083-2100.
16. Митрин, Б. И. Распределение контактных напряжений под круглой пластиной, лежащей на мягком слое / Б. И. Митрин, С. С. Волков // Вестник Дон. гос. техн. ун-та. — 2013. — № 5/6 (74). — С. 14-25.
17. Айзикович, С. М. Асимптотическое решение одного класса парных уравнений / С. М. Айзикович // Прикладная математика и механика. — 1990. — Т. 54, вып. 5. — С. 872-877.
18. Лурье, А. И. Теория упругости / А. И. Лурье. — Москва : Наука, 1970. — 824 с.
19. Белубекян, М. В. К вопросу колебаний неоднородной по толщине пластинки / М. В. Бе-лубекян // Известия национальной академии наук Армении. Механика. — 2002. — Т. 55, № 3. — С. 34-41.
20. Цейтлин, А. И. Об изгибе круглой плиты, лежащей на линейно деформируемом основании / А. И. Цейтлин // Известия АН СССР. Механика твёрдого тела. — 1969. — № 1. — С. 99-112.
21. Айзикович, С. М. Осесимметрическая задача о вдавливании круглого штампа в упругое, неоднородное по глубине полупространство / С. М. Айзикович, В. М. Александров // Известия АН СССР. Механика твёрдого тела. — 1984. — Т. 19, № 2. — С. 73-82.
22. Айзикович, С. М. Асимптотическое решение задачи о взаимодействии пластины с неоднородным по глубине основанием / С. М. Айзикович // Прикладная математика и механика. — 1995. — Т. 59, вып. 4. — С. 688-697.
Материал поступил в редакцию 12.03.2014.
References
1. Hager, A. M., Davies, M. Short-Fibre Reinforced, High-Temperature Resistant Polymers for a Wide Field of Tribological Applications. In K. Friedrich, ed.: Advances in Composite Tribology, Amsterdam, 1993, pp. 107-157.
2. Friedrich, K., Reinicke, R., Zhang, Z. Wear of polymer composites. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2002, vol. 216, iss. 6, pp. 415-426.
3. Goryacheva, I. G., Myshkin, N. K., Torskaya, E. V., Kornev, Y. V., Gutsev, D. M., Kudritskiy, V. G., Kovaleva, I. N. Modeling friction of tribological composite coatings. Journal of Friction and Wear, 2012, vol. 33, no. 6, pp. 407-414.
4. Vasilyev, A. S., Sadyrin, E. V., Vasilyeva, M. E. Krucheniye uprugogo poluprostranstva s mno-gosloynym pokrytiyem periodicheskoy struktury. [Torsion of elastic half-space with multilayered coating of periodic structure.] Vestnik of Don State Technical University, 2013, no. 5/6 (74), pp. 6-13 (in Russian).
5. Gorbunov-Posadov, M. I. Raschet balok i plit na uprugom poluprostranstve. [Calculation of beams and plates on an elastic half-space.] Prikladnaya matematika i mekhanika [Journal of Applied Mathematics and Mechanics], 1940, vol. 4, iss. 3, pp. 61-80 (in Russian).
6. Ishkova, A. G. Ob izgibe polosy i krugloy plastiny, lezhashchikh na uprugom poluprostranstve. [A strip and a circular plate bending on an elastic half-space.] Inzhenernyy sbornik, 1960, vol. 23, pp. 171— 181 (in Russian).
7. Grebenschikov, V. N. Raschet krugloy plastinki na uprugom poluprostranstve [Calculation of circular plate on elastic half-space.] Teoriya rascheta i nadejnost' priborov : sb. tr. II Saratovskoy obl. konf. mol. uchenyh. [Theory of calculation and reliability of devices : Proc. II Young Scientists Conf.] Saratov, 1969, pp. 48—51 (in Russian).
8. Aleksandrov, V. M., Shatskikh, L. S. Universalnaya programma rascheta izgiba balochnyh plit na lineino-deformiruemom osnovanii. [Universal program for calculation of plate bending on linear deformable foundation.] Trudy 7 Vsesoyuznoi conferencii po teorii obolochek i plastinok. [Proc. VII All-Soviet Conference on Theory of Plates and Shells.] Moscow, 1970, pp. 46—51 (in Russian).
9. Shatskikh, L. S. K raschetu izgiba plity na uprugom sloe. [About a plate bending on an elastic layer.] Izv. Akad. Nauk SSSR, Mekh. Tverd. Tela [Mechanics of Solids], 1972, vol. 2, pp. 170—176 (in Russian).
10. Aleksandrov, V. M., Vorovich, I. I., Solodovnik, M. D. Effektivnoe reshenie zadachi o ci-lidricheskom izgibe plastinki konechnoy shiriny na uprugom poluprostranstve. [Effective solution of problems on cylindrical bending of plates of finite width on an elastic half-space.] Izv. Akad. Nauk SSSR, Mekh. Tverd. Tela [Mechanics of Solids], 1973, vol. 4, pp. 129—138 (in Russian).
11. Aleksandrov, V. M., Solodovnik, M. D. Asymptotic problem of the cylindrical bending of a plate of finite breadth in an elastic half-space. Soviet Applied Mechanics, 1974, vol. 10, iss. 7, pp. 749—754.
12. Bosakov, S. V. The solution of the contact problem for a circular plate. Journal of Applied Mathematics and Mechanics, 2008, vol. 72, no. 1, pp. 59—61 (in Russian).
13. Kashtalyan, M., Menshykova, M. Effect of a functionally graded interlayer on three-dimensional elastic deformation of coated plates subjected to transverse loading. Composite Structures, 2009, vol. 89, iss. 2, pp. 167—176.
14. Kashtalyan, M. Three-dimensional elasticity solution for bending of functionally graded rectangular plates. European Journal of Mechanics A/Solids, 2004, vol. 23, iss. 5, pp. 853—864.
15. Silva, A. R. D., Silveira, R. A. M., Goncalves, P. B. Numerical methods for analysis of plates on tensionless elastic foundations. International Journal of Solids and Structures, 2001, vol. 38, iss. 10—13, pp. 2083—2100.
16. Mitrin, B. I., Volkov, S. S. Raspredeleniye kontaktnykh napryazheniy pod krugloy plastinoy, lezhashchey na myagkom sloye. [Contact stress distribution under circular plate on a soft layer.] Vestnik of Don State Technical University, 2013, no. 5/6 (74), pp. 14—24 (in Russian).
17. Aizikovich, S. M. An asymptotic solution of a class of coupled equations. Journal of Applied Mathematics and Mechanics, 1990, vol. 54, iss. 5, pp. 719—724.
18. Lourie, A. I. Teoriya uprugosti. [Theory of Elasticity.] Moscow, 1970, 824 p. (in Russian).
19. Belubekyan, M. V. K voprosu kolebaniy neodnorodnoy po tolshchine plastinki. [On the problem of vibration of the plate with nonuniform thickness.] Izvestiya natsionalnoy akademii nauk Armenii. Mekhanika. [Mechanics. Proceedings of National Academy of Sciences of Armenia], 2002, vol. 55, no. 3, pp. 34—41 (in Russian).
20. Tseitlin, A. I. Ob izgibe krugloy plity, lejaschey na lineyno-deformiruemom osnovanii. [Bending of a circular plate lying on a linearly deformable foundation.] Izvestiya AN SSSR. Mekhanika tverdogo tela [Mechanics of Solids], 1969, vol. 1, pp. 99—112 (in Russian).
21. Aizkovich, S. M., Aleksandrov, V. M. Axially symmetric problem of the indentation of a circular die into an elastic nonuniform in depth half-space. Mechanics of Solids, 1984, vol. 19, iss. 2, pp. 73—82 (in Russian).
22. Aizikovich, S. M. Asymptotic solution of the problem of the interaction of a plate with an inhomogeneous in depth foundation. Journal of Applied Mathematics and Mechanics, 1995, vol. 59, iss. 4, pp. 661-669 (in Russian).
AXISYMMETRIC BENDING OF CIRCULAR SANDWICH PLATE ON ELASTIC FOUNDATION WITH COMPLICATED STRUCTURE*
S. M. Aizikovich, S. S. Volkov, A. V. Melkonyan
The approximate closed-form solution to the problem on a circular multilayer (sandwich) constant-thickness plate bending on an elastic foundation with a complicated structure is obtained. The plate is bent under the axisymmetric distributed load, and the foundation reaction. The elastic foundation is a uniformly irregular in thickness layer (coating) based on a homogeneous half-space (substrate). Young's modulus value at the interface of the coating and the substrate has a significant leap. Two different cases of the boundary conditions are considered for the plate: fixed and free edge conditions. The constructed approximate analytical solution to the problem is effective for a wide range of both geometric (the inhomogeneous layer thickness and the plate radius) and physical parameters (plate flexibility, and elastic properties of the coating and the substrate). The contact problem is reduced to the system of the integro-differential equation solution through the integral transformation method. The obtained formulas can be used for calculating the contact interaction characteristics between a multilayer plate and a foundation with a complex structure in various cases of the boundary conditions, and various loads applied to the plate.
Keywords: inhomogeneous materials, sandwich plate, functionally-graded coating, axisymmetric problem, analytical methods, approximate analytical solution.
* The research is supported by RFFI grant no. 13-08-90916-mol_in_nr.
