УДК 539.3 Вестник СПбГУ. Математика. Механика. Астрономия. 2023. Т. 10 (68). Вып. 1
МЯС 74М20, 26А33
Удар жесткого шара по бесконечной пластинке Кирхгофа — Лява с учетом объемной и сдвиговой релаксации*
М. В. Шитикова
Московский государственный строительный университет, Российская Федерация, 129337, Москва, Ярославское шоссе, 26 Воронежский государственный технический университет, Российская Федерация, 394006, Воронеж, ул. 20-летия Октября, 84
Для цитирования: Шитикова М. В. Удар жесткого шара по бесконечной пластинке Кирхгофа — Лява с учетом объемной и сдвиговой релаксации // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2023. Т. 10(68). Вып. 1. С. 139154. https://doi.org/10.21638/spbu01.2023.112
Рассматривается задача о нормальном низкоскоростном ударе жесткого шара по бесконечной вязкоупругой пластинке Кирхгофа — Лява. Динамическое поведениие вяз-коупругой пластинки описывается моделью стандартного линейного тела с дробными производными. Параметр дробности, определяющий порядок дробной производной, учитывает изменение вязкости материала пластинки в зоне контакта в процессе удара. Местное смятие материала пластинки и контактная сила определяются по обобщенной теории Герца. Используя алгебру операторов Ю. Н. Работнова, а также учитывая действие объемной и сдвиговой релаксации, удается получить интегральное уравнение относительно местного смятия контактирующих тел. Приближенное решение этого уравнения позволяет найти временные зависимости не только для контактного смятия, но и для контактной силы.
Ключевые слова: низкоскоростной удар, вязкоупругая пластинка Кирхгофа — Лява, модель стандартного линейного тела с дробными производными, алгебра дробных операторов Работнова.
1. Введение. Задачи ударного взаимодействия относятся к наиболее сложным динамическим задачам механики деформируемого твердого тела [1], для решения которых используются различные аналитические, численные и экспериментальные методы.
Одним из наиболее эффективных подходов, позволяющих аналитически исследовать динамические свойства соударяемых тел, является волновой подход. В обзорной статье [2], посвященной анализу работ, использующих этот подход, отмечалось, что одной из первых в этом направлении была работа К. Зинера (^ Zener) [3], где рассмотрена задача об ударе упругого шара по упругой бесконечной пластинке Кирхгофа —Лява. Позже подход Зинера был обобщен в работе [4], где изучалась задача об ударе жесткого шара по вязкоупругой бесконечной пластинке Кирхгоф-фа — Лява. Однако в этой работе не дано четкого опеределения вязкоупругой мо-
* Работа выполнена при финансовой поддержке Министерства науки и высшего образования РФ (грант №Р20Ы-2020-0007).
© Санкт-Петербургский государственный университет, 2023
дели, которая использовалась для описания динамического поведения бесконечной пластинки.
В обзоре [5] отмечается, что для описания вязкоупругих свойств соударяемых тел в последнее время с успехом используются различные модели вязкоупругости, содержащие производные дробного порядка [6, 7]. Некоторые из этих задач включены в «Справочник по дробным производным и их приложениям» [8] и «Энциклопедию по механике сплошных сред» [9—11].
Во всех предыдущих исследованиях ударного взаимодействия вязкоупругих тел с использованием моделей с дробными операторами (см. работы [9-12], а также и другие статьи, процитированные в них), предполагалось, что объемные деформации остаются упругими, что позволяет пренебречь объемной релаксацией по сравнению с релаксацией при сдвиге. Однако это предположение не всегда является обоснованным, поскольку для композитных материалов объемная релаксация может являться весьма существенной [13-15].
В данной работе этот недостаток устранен: в качестве вязкоупругой модели, описывающей поведение пластинки при ударном воздействии, выбрана трехмерная модель стандартного линейного тела с дробными производными, которая позволяет описать и сдвиговую, и объемную релаксации. Хотя трехмерные модели вязкоупругости с дробными операторами были известны и ранее [14, 16-18], однако в задачах ударного взаимодействия используются, насколько известно автору, впервые.
2. Основные соотношения трехмерной модели стандартного линейного тела с дробными производными. Для девиаторных частей а'^ и е'^ тензоров напряжений а^ и деформаций е^, а также для их шаровых частей акк и екк для модели стандартного линейного тела с производными дробного порядка можно записать следующие соотношения:
4 + т^'а'ц =2^0 {е'ц + т^В"'е^) , акк + т]2В7акк = 3Ко (екк + т^В1 екк)
(1) (2)
где тб1 и те2 — времена сдвиговой и объемной релаксации; та\ и та2 — времена сдвиговой и объемной ретардации; мо и Ко — релаксированные модуль сдвига и объемный модуль,
— дробная производная Римана — Лиувилля [19]; 7 (0 < 7 < 1) — параметр дробности; у(Ь) — произвольная функция; Г(1 — 7) — гамма-функция.
Если ввести в рассмотрение безразмерный дробный оператор Ю. Н. Работно-ва [20]
= л , * (г = 1,2) (4)
1 + т7В
и учесть формулу [17]
т?гВ~'-Э1 (т£ = 1—э! (т?г) (г = 1, 2),
(5)
справедливость которой устанавливается непосредственной проверкой, то соотношения (1) и (2) можно переписать в виде
2^с
1
Мс
Мо
К"!)
(6)
а
И
оо
&кк = 3К0
1 -
К^-Ка
К
э; (т?2)
ъ3>
(7)
где введены обозначения для нерелаксированных модуля сдвига ¡¡ж и объемного модуля Кж:
^7 „Л
(8)
/Лоо МО -у ?
Те1
К со = Ко .
е2
Используя выражения (6) и (7), запишем операторы сдвига ¡1 и объемного расширения — сжатия К
¡ = ¡ж [1 -^1 э*7 (Тб1)] , (9)
К = Кж [1 - э; (т?2)\
где
¡ж — ¡¡о
^е2
Кур-Ко
к0
(10)
(11)
¡ж Кж
В дальнейшем нам понадобятся операторы ¡-1 и К-1, обратные к операторам ¡1 и К. Чтобы их найти, нужно воспользоваться определением обратных операторов
■ 1=1, К-1 ■ К= 1,
а также учесть формулу [17]
э* (т7 V Э* (т'-) =
эу ^ еъ; эу ^ аъ'
.7 \ _ ег 7
э; (т7) — т7 э; (т^ъ)
(г = 1,2),
(12)
(13)
еъ аъ
справедливость которой устанавливается непосредственной проверкой, поскольку
1 А В
э* (т1.). э* (т1.) =- _
7 ^ ег> 7 ^ "> (1 + Г7Д7) (1 + Т^В 7) 1 + Т7£7 1 + Т^В 7 '
+
где
А
еъ - аъ
В
Т 7
т ■
еъ - аъ
Выбирая операторы ¡1 1 и К 1 в виде
= ¡Ж [1 + Vа 1 э*1 (т^)] , (14)
К-1 = КЖ [1+ э; (Та72 )] , (15)
где уаъ и таъ — пока неизвестные константы, подставляя выражения (9) и (14), (10) и (15) в формулу (12) и учитывая соотношение (13), в результате получим
^ъ э*1 (т1ъ)
откуда следует, что
1- V,
ег -у -у
Т ■ — Т ■
аъ - еъ
— Vеъ э*; (т7)
1 - Vа
аъ - еъ
0,
1 - V,
ег -у -у
Т ■ — Т ■
аъ - еъ
(16)
(17)
1 - Vа
7 7 Т ■ — Т ■
аъ - еъ
(18) 141
7
7
7
Т
7
7
Т
7
7
Т
Из уравнений (17) и (18) находим
7
Т ■ =
7 Т ■
7 - U0
1 - Vi 1 - Ve
,7 / \ 7
Tei \ ' Т 4
(Í = 1,2) (19)
= 1-^, ^ =1 + ^ (г = 1,2). (20)
таг } \ тег /
Первая формула в (19) совпадает с формулами (8), а вторая формула в (19) имеет вид (сравни с формулами (11))
Мс — Мо Кс — Ко , ,
=-, гЛт2 =-^-• (21)
Мо Ко
3. Алгебра обобщенных операторов Ю. Н. Работнова. Зная операторы (14) и (15), можно найти оператор податливости J по формуле
-1_ 13К + и _ 1 (~-1 , 1^-1
где E — оператор Юнга.
Подставляя выражения (14) и (15) в формулу (22), находим
J = J.х [1 + U! э* (т^) + П2 э; (т;2)], (23)
где
1 2 1 ^<т2 0 s
"-1 = о-Т~ = О (! +^00)^1, "-2 = —= - (1 -21У00)1Уа2-
3 3 9 3
Чтобы определить E, нужно использовать свойство взаимообратных операторов
J ■ E=1, (24)
а также учесть формулу (13). В результате получим
E = Еж [1 - mi э* (t^) — m2 Э* (t¡2)] , (25)
где величины th и t^ определяются из квадратного уравнения
П1Х ^7—— = 1, x = tl (¿ = 1,2), (26)
7 7
Тя1 — x тя2 — x
или
(1+ ni + n2) x2 — [tYi(1+ n2) + r¿2(1+ ni)] x + tYit¿2 =0, (27)
а величины mi и m-2 находятся из линейной системы уравнений
TY TY
+ Т1 W2 = Ta1 — el Ta1 — te2
(28)
TY TY
ct2 mi + T a2„ m2 = 1
TY tY "1-1 1 TY tY Ta2 — el Ta2 — te2
при условии, что величины и уже известны.
V.
или
Теперь подсчитаем оператор, пропорциональный оператору цилиндрической жесткости, который будет использоваться в дальнейшем при решении задачи на удар
Е Е Е .
+ (29)
1 — V2 2(1 + V) 2(1 — V)'
где
1
2 3Л + Д 9 К '
1
2(3К + ¡¡)
Е =
9 Л1
1 +1 9 Л 1 — V 3К + 41 3Л + ¡1
С учетом выражений (30) формулу (29) перепишем в виде
(30)
Е1
1- V
= ¡V +
9К}1
= 1 + 1
1
причем
где
3 Кро + 4/Хо
9КЖ1<Ж
п 1
ЗЛоо + 4/х0
п 2
ЗЛоо + 4/хс
Оператор V 1 находится в виде
V-1 =
9К{1 Ж+ 4(1
= «ж1 [1 — т э; (¿70 — ^2 э; (^2)]
(31)
5=^щг= I 1 + = [1+(тЛ)+' (32)
(33)
где величины реъ (г = 1, 2) определяются из квадратного уравнения (27), в котором величины П1 и П2 заменены соответственно на П1 и П2, а величины т1 и т2 определяются из линейной системы уравнений (28), в которой величины и ^2 соответственно на р1 и р2.
Наконец, подставляя выражение (33) в формулу (31), окончательно получим
Е1
1- V2
1 — £ пч э; (Ч
3 = 1
Р = т1
еЗ = Те1,
где
¿п
¡ж«ж + 1
П1
т 1
¡ж «ж + 1
П2
т 2
¡ж «ж + 1
Обратный оператор к оператору (34) имеет вид
3
1-Г Ё
= й-
1+£ Сз э; (Та
3=1
Пз
¡¡ж«жVе1 ¡¡ж «ж + 1
где величины Т^, Т^ и Т03 находятся из кубического уравнения
— X р — X — X
П1х п2 х пзх + ^ Ь
ие1
ье2
г/
— 1, х = Т2ъ (г = 1, 2, 3)
е3
(34)
(35)
(36) 143
«
й
эо
«
оо
1
или
(1 - П1 - П2 - Пз) х3+
+ [(% + Сз) П1 + (¿23 + & ) П2 + (& + &) П3 - & - % - -2 +
+ [¿72^з(1 - П1)+ *7з%(1 - П2)+ ^1^2(1 - пз)] х - ¿71^2^3 = 0, (37) а величины £1, (2 и (3 находятся из системы трех линейных уравнений:
С1П1 , С271 , Сз^1
+
+
ГУ та ■ 1-у гр-у 1-у гр-у Ьс1 1а1 Ьс1 1а2 Ьс1 1 а3
(1^2 , С2^2 , С3^2
+
+
ГУ Т>7 -1-у гр-у 1-у гр-у ье2 1а1 ье2 1 а2 ье2 1 а3
С1П3 , С2*7з , Сз7з
+
+
ГУ т>7 ' -П Т"1 ' -П Т~< Ъе3 -1а1 Ъе3 - 1 а2 Ъе3 - 1 а3
-1,
-1,
-1.
(38)
Применяя теорему Виета к уравнению (37), получаем полезную формулу, которая обобщает первую формулу в (20):
/ 3 А \Т
' П * з
3=1
. П
\з=1 /
1 -£ пз.
з=1
(39)
Из формулы (39) видно, что при 7 = 0, т.е. при отсутствии вязкости,
3
Т,Пз =°.
з=1
(40)
Если считать, наоборот, что оператор (35) известен, а оператор (34) необходимо найти, то вместо формулы (39) получаем
/ 3 \ 1
3=1
. П из .
з=1
1 + £ Сз
з=1
(41)
и при 7 = 0
Е<з =0. (42)
з=1
Формула (41) обобщает вторую формулу в (20).
В заключение данного раздела покажем, как представить оператор (4) в виде, удобном для приложений. С этой целью умножим числитель и знаменатель дроби, стоящей справа в формуле (4), на Г'т-7, где
г* (+ _ +>)-У-1
гт=к Чхтг^'
(43)
— дробный интеграл [19].
Учитывая, что D7I7 = I7D7 = 1, перепишем выражение (4) в виде
IY Т-7
Э* (т7) = -—^--(44)
7 ^ ег) 1_(_/7тТ7) 1 ^
Интерпретируя формулу (35) как сумму бесконечно убывающей геометрической
* (т7)
7 V ег)
прогрессии со знаменателем —I7т-7, представим оператор Э* (те7) в виде
э7
* (тЛ) = Е- 1)nT-i7(n+1)I7(n+1) (45)
= 0
э; (rZ) y(t) = fo э7 y(t')àt', (46)
где
=—£ г[7(„+1)] (47)
— дробноэкспоненциальная функция Ю. Н. Работнова [21], которая при y = 1 превращается в обычную экспоненту
(-t/rei) = — exp(-t/rei).
Тег
4. Постановка задачи и определяющие уравнения. В работе [3] К. Зинером (C. Zener) была рассмотрена задача о нормальном низкоскоростном ударе шара по достаточно протяженной упругой пластинке, используя теорию тонких пластин. При этом предполагалось, что процесс ударного взаимодействия мишени (пластинки) и ударника (шара) заканчивается раньше, чем в область контакта вернутся волны, отраженные от границ пластинки (в [22] отмечалось, что для этого требуется, чтобы масса ударника не была слишком большой по сравнению с массой мишени). Было показано, что прогиб такой пластинки в месте удара шара определяется в виде
^ = lF{t,)i>- (48)
где u(t) — перемещение срединной плоскости; h — толщина пластинки; g — ее плотность; E и v — ее модуль упругости и коэффициент Пуассона; F(t) — контактная сила, которая, согласно теории Герца, связана с местным смятием материала пластинки x(t) при помощи формулы
m = A\^ Е (49)
3 1 — v2
где R — радиус жесткого шара.
Уравнение движения жесткого шара массы m имеет вид
m = -~F(t), (50)
m
где точки над величинами означают производные по времени t.
или
Перемещение центра шара г (г) связано с перемещением центра пластинки и(г) и местным смятием материала пластинки в зоне контакта х(г) следующим соотношением:
х(г) = г(г) - и(г).
(51)
Согласно принципу Вольтерра, для вязкоупругой пластинки упругие постоянные в формулах (48)-(50) нужно заменить операторами (34) и (35). Оператор, входящий в уравнение (48), можно разложить в ряд Неймана и учесть только два первых члена в силу кратковременности ударного процесса, т. е.
1 ~2 \ 1/2 1 — V2 \
Е
сГ1/2
1 + Е Сз э; 1
з=1
1/2
сГ1/2
1 3 ( 1 + 2ЕС3Э; {т:0
2 з=1
Контактную силу с учетом формулы (34) можно представить в виде
Р (г)
3/2
(г) -Е пз (
з=1 -]о
г-г'
х3/2 (t")dt"
Тогда с учетом выражений (52) и (53) уравнение (48) примет вид
. (52)
(53)
х001и(г)
з=1 ^о V *
Ф
3 г"
- пз э -
з=1 о
з
/\ г г'
¿г'+
/ ¿о
г'' - г'''
где азоо — ш2 у Зд .
Интегрируя дважды уравнение (50), находим
г
с3/2(г'' )-
х3/2 (г''')сг'''
&'') ¿г', (54)
1 / Ейс
Ф) = -— /
т ) о
Подставляя в уравнение (55) выражение (53), имеем
. Ал/Шоо Г*
3 , г'
х3/2(г') -Т, пз э7
г' г'
ге
х3/2 (г'')сг''
(55)
(г -г')сг'. (56)
з=1 -0 х
Подставляя теперь соотношения (54) и (56) в (51) и учитывая начальные усло-
с(0) = 0, х(0) = У0,
(57)
где Уо — скорость шара в момент удара, получим интегральное уравнение для определения функции х(г):
г
оо
3
о
х(г) =
3т
г3/2
з г'
Пз
х3/2
(г') -Е
з=1 з , г ^з
з=1 3 г'
Пз
'о
г' - г'
х3/2(г'' )А''
(г-г' )3£-
о
з=1
о
гг
г' - г"
г ез
/\ г г'
'аз / ./0
з з
з=1
Пз
г'' г'''
х3/2(г'' )&'' х3/2 (г'')-х3/2 (гт)д£''
¿г'-
&'') ¿г'. (58)
5. Приближенные вычисления. Поскольку ударный процесс является кратковременным процессом, то дробноэкспоненциальную функцию (46) в уравнении (58) можно заменить первым членом ряда (47) [23]
г
г7
-1
Эт1
(59)
Если при решении уравнения (58) использовать метод последовательных приближений, а в качестве первого приближения выбрать выражение
х(г) = Уог,
то, подставляя его в правую часть уравнения (58) и учитывая, что
1'\г - гу-1^')3'2^ = ±1*^(1- (¿')1/2^' =
(60)
7 \3 5
находим
= у^ь -
16 л/Д4х> ^3/2^7/2 4УД^00^3/Ч(1/3-1/57)^/2+-у
105 т
5
(5/2 + 7)
т 7 (5/2 + 7 )(7/2 + 7)
зес
7
+ § ^ ^ - - (61)
7
5/2 + 7 7/2+ 7;
где
<5 - Д
~~ л! ' 7 Т-Т •
Г(^) з=1 гз
ГМ^ та
\ I ! з=1 аз
При 7 = 0 из соотношения (61) с учетом формул (40) и (42) получаем зависимость местного смятия упругой пластинки от времени:
хЩг) = ъг - - ¡^У^.
105 т
(62)
оо
-
о
— т
о
-
--ае
2
-
о
Выражение (62) можно получить и из соотношения (58), если учесть формулы (40) и (42), а также, что [см. формулу (59)]
lim Э7 = lim Э7 f-^^) = S(t -1'), (63)
0 7 V tj ) Y^0 7 V Taj J
где ö(t — t') — дельта-функция Дирака.
Полагая в выражении (62) функцию x(0) равной нулю, находим время контакта жесткого шара и упругой бесконечной пластинки:
где
1 V2v/^£ßoo Дифференцируя соотношение (62) по времени t
внивая по
ция x(0)(t) достигает максимального значения
« ,65,
и приравнивая полученное выражение нулю, находим время г^ах, при котором функ-
4°L«40) 1 _ ^ vi^.co) (66)
max 2 I 45 2 I v '
.(0) ( 1 \2/3 (0) ,.(0) ч
где to = —Lj= , а само максимальное значение Xmaxlimax) имеет вид
2 Vo )
м«(I - • (er)
2 у 5 189 тавоо 2 у
Используя выражение (61) и пренебрегая членами более высокого порядка малости по времени г, получим соотношение х(1\г) для произвольных значений величины 7 (0 < 7 < 1)
« - - - аеоо^2А7(1/5"1/77)^.
105 т 5 7
(68)
Приравнивая нулю выражение (68), находим
Дифференцируя соотношение (68) по времени г, имеем
« у0 _ А 3/2^/2_
¿г 15 т
_ ae^Vb3/2^/2 - аеоо^2А7(1/5~1/7>у)(5/2 + >у)^. (70)
Y
Приравнивая выражение (70) нулю, получим
«а. - е - (4«.) -. (71)
Подставляя выражение (71) в соотношение (68), находим
-м (t™ ) - rtC°) ) - УпД (1/5~1/77) Г- (- +-Л
"maxl max/ _ хтахГтах^ ^О^у 3 \2 /
• (72)
Теперь получим те же самые величины при значении 7 =1. Из выражения (68) при 7 =1 находим
*<!>(*) « ^ - - + 1-ооД^ (73)
где Д1 = Е F2"-
j=i
Соотношение (73) позволяет определить все основные характеристики ударного взаимодействия:
= + (74)
= (75)
= - ^VoA^f. (76)
Из полученных выражений видно, что с увеличением вязкости время контакта шара и пластинки увеличивается, а максимумы у контактной силы уменьшаются и смещаются влево по временной оси, что находится в полном соответствии с результатами, полученными в работе [4] (см. рис. 2 в [4]).
6. Дискуссия и выводы. В работе [4] при решении аналогичной задачи на удар не конкретизирована вязкоупругая модель материала, из которого состоит пластинка. В данной же работе для описания вязкоупругих свойств материала пластинки была использована модель стандартного линейного тела с дробными производными, что позволило подтвердить справедливость основных результатов, полученных Дж. В. Филлипс и Х. Х. Кальвит (J. W. Phillips и H. H. Calvit) [4]. Только в качестве основных параметров в данной работе используется не тангенс угла механических потерь, а порядок дробной производной. Однако можно показать, что эти величины связаны друг с другом.
С этой целью подействуем оператором J, заданным формулой (23), на функцию егшt, предварительно заменив в формуле (23) дробные производные (3) на
Используя формулы (3) и (45), можно получить
„гшЬ
э;+ (т7 ) в*- = 6 п7 = ]Г(-1 )пт-7(п+1)/у(п+1)е<ш4 {г = г 2) (78)
1 + Tai n
+ n=0
или
то
(tYÍ) e™ = e^Y.(-!)" (i^i)-(n+1), (79)
n=0
где
rt (t-t')f-1 w
W) = / v J ^ y(f)dt'
— дробный интеграл [19].
Сумму, стоящую в формуле (79), можно интерпретировать как сумму бесконечно убывающей геометрической прогрессии со знаменателем d = — {iwTai) 1, т. е.
£(-1)" (^гН("+1) = = ТТ7~~-v7' (80)
Таким образом,
(~< ^icot eiwí í еер + совф . sin^ \
=4+ УТаг) е
aiJ \ + {гшта1)1 + в>1.+2сжф ее;7 + ее7, + 2 cos V у
(81)
где eeai = iOTfji, ф =
Учитывая формулу (81) в соотношении
Jeiut = J (iu)eiut, (82)
получим выражение для комплексной податливости
JM = J'H - = boo Ul + v00)^+^77+[1tKlMl/)7]COS^ +
3 [ + + 2 cos ф
, , ае72 + ае~27 + [1 + («Waее2)7] cos )
+ <2 +2 cos ф J
- Д Joo (2(1 + Уоо) K^i/^i)7 - + (1 _ 2Уоо) [(^2/^Г - 1] sin 1 (83)
3 l Ж-7 + + 2 cos ф Ж-27 + ж72 + 2 cos ф J
где ж« = wrei.
Зная реальную и мнимую части комплексной податливости, входящие в выражение (83), можно найти тангенс угла механических потерь, поскольку
J 'И
Чтобы показать зависимость максимальной величины tg¿ от параметра дробности y, в работе [14] был рассмотрен численный пример при следующих значениях параметров вязкоупругой модели: = 0.3, цо/^ж = Ко/Кто = 0.8,
150
те\/ те2 = та1 /та2 = 103, ш = 1. На рис. 1 в [14] приведен тангенс угла механических потерь tgJ как функция от 1пже2, где в качестве параметра взята величина 7. Как и утверждалось, с увеличением параметра дробности 7 от 0 до 1 (это соответствует увеличению вязкости) максимальное значение величины tgJ также увеличивается. Это видно как на пиках, связанных со сдвиговой релаксацией (левые пики на рис. 1 в [14] ), так и на пиках, связанных с объемной релаксацией (правые пики на рис. 1 в [14]), при этом сдвиговый пик был в 5.5 раза больше объемного.
Таким образом, в данной работе была рассмотрена задача о нормальном ударе жесткого шара по бесконечной вязкоупругой пластинке Кирхгофа — Лява, демпфирующие свойства которой описываются трехмерной моделью стандартного линейного тела с дробными производными, что позволило учесть объемную и сдвиговую релаксацию материала в процессе ударного взаимодействия.
Литература
1. Гольдсмит В. Удар. Теория и физические свойства соударяемых тел, пер. с анг. Москва, Изд-во лит. по строительству (1965).
2. Rossikhin Yu.A., Shitikova M.V. Transient response of thin bodies subjected to impact: Wave approach. The Shock and Vibration Digest 39 (4), 273-309 (2007). https://doi.org/10.1177 /0583102407080410
3. Zener C. The intrinsic inelasticity of large plates. Physical Review 59 (8), 669-673 (1941). https://doi.org/10.1103/PHYSREV.59.669
4. Phillips J.W., Calvit H.H. Impact of a rigid sphere on a viscoelastic plate. Journal of Applied Mechanics 34 (4), 873-878 (1967). https://doi.org/10.1115/L3607850
5. Shitikova M. V. Wave theory of impact and Professor Yury Rossikhin contribution in the field (A memorial survey). Journal of Materials Engineering and Performance 28 (6), 3161-3173 (2019). https://doi.org/10.1007/s11665-018-3824-6
6. Шитикова М. В. Обзор вязкоупругих моделей с операторами дробного порядка, используемых в динамических задачах механики твердого тела. Известия РАН. Механика твердого тела 1, 3-40 (2022). https://doi.org/10.31857/S0572329921060118
7. Shitikova M. V., Krusser A. I. Models of viscoelastic materials: a review on historical development and formulation. Advanced Structured Materials 175, 285-326. (2022). https://doi.org/10.1007/978-3-031-04548-6-14
8. Rossikhin Yu.A., Shitikova M.V. Fractional calculus in structural mechanics. In: Handbook of Fractional Calculus with Applications. Vol. 7: Applications in Engineering, Life and Social Sciences. Part A, 159-192, D. Baleanu, A.M.Lopes (eds). Berlin, De Gruyter (2019). https://doi.org/10.1515/9783110571905-009
9. Rossikhin Yu.A., Shitikova M.V. Classical beams and plates in a fractional derivative medium, Impact response. In: Encyclopedia of Continuum Mechanics, 294-305. H. Altenbach, A. (Ochsner (eds). Germany, Springer-Verlag GmbH (2020). https://doi/org/10.1007/978-3-662-53605-6-86-1
10. Rossikhin Yu. A., Shitikova M. V. Fractional derivative Timoshenko beams and Uflyand-Mindlin plates, Impact response. In: Encyclopedia of Continuum Mechanics, 962-971, H. Altenbach, A. Ochsner (eds). Germany, Springer-Verlag GmbH (2020). https://doi/org/10.1007/978-3-662-53605-6-87-1
11. Rossikhin Yu.A., Shitikova M.V. Collision of two spherical shells. Fractional operator models. Encyclopedia of Continuum Mechanics 324-332, H. Altenbach, A. (Ochsner (eds). Germany, SpringerVerlag GmbH (2020). https://doi/org/10.1007/978-3-662-53605-6-88-1
12. Rossikhin Yu. A., Shitikova M. V. Two approaches for studying the impact response of viscoelastic engineering systems: An overview. Computers and Mathematics with Applications 66 (5), 755-773 (2013). https://doi.org/10.1016/jxamwa.2013.01.006
13. Мешков С. И., Пачевская Г. Н., Постников В. С. К выявлению объемной релаксации при динамических испытаниях. Ученые записки Ереванского государственного университета 2, 10-14 (1968).
14. Мешков С. И., Пачевская Г. Н. К учету объемной релаксации методом внутреннего трения. ПМТФ 2, 80-82 (1967).
15. Shitikova M.V., Popov 1.1., Rossikhin Yu.A. Theoretical and experimental evidence of the bulk relaxation peak on the loss tangent versus frequency diagrams for concrete. Mechanics of Advanced, Materials and Structures (2022). https://doi.org/10.1080/15376494.2021.2012857
16. Makris N. Three-dimensional constitutive viscoelastic laws with fractional order time derivatives. Journal of Rheology 41, 1007-1020 (1997). https://doi.org/10.1122/L550823
17. Rossikhin Yu. A., Shitikova M. V. Features of fractional operators involving fractional derivatives and their applications to the problems of mechanics of solids. In: Fractional Calculus: History, Theory and Applications, ed. by R. Daou, M.Xavier. Chapter 8, 165-226 (2015). New York, Nova Science Publ.
18. Alotta G., Barrera O., Cocks A., Di Paola M. On the behavior of a three-dimensional fractional viscoelastic constitutive model. Mecca,nica, 52 (9), 2127-2142 (2017). https://doi.org/10.1007/s11012-016-0550-8
19. Самко С. Г., Килбас А. А., Маричев О. И. Интегралы и производные дробного порядка и некоторые их приложения. Минск, Наука и техника (1987).
20. Rossikhin Yu.A., Shitikova M.V. Centennial jubilee of Academician Rabotnov and contemporary handling of his fractional operator. Fractional Calculus and Applied Analysis 17 (3), 674-683 (2014). https://doi.org/10.2478/s13540-014-0192-2
21. Работнов Ю.Н. Равновесие упругой среды с последействием. ПММ 12 (1), 53-62 (1948).
22. Джонсон К. Механика контакного взаимодействия, пер. с англ. Москва, Мир (1989).
23. Работнов Ю. Н. Элементы наследственной механики твердых тел. Москва, Наука (1977).
Статья поступила в редакцию 26 мая 2022 г.;
доработана 24 июля 2022 г.; рекомендована к печати 8 сентября 2022 г.
Контактная информация:
Шитикова Марина Вячеславовна — д-р физ.-мат. наук, проф.; mvs@vgasu.vrn.ru
Impact of a rigid sphere on an infinite viscoelastic
Kirchhoff—Love plate considering volume and shear relaxations*
M. V. Shitikova
Moscow State University of Civil Engineering,
26, Yaroslavskoe shosse, Moscow, 129337, Russian Federation
Voronezh State Technical University,
84, ul. 20-letiya Oktyabrya, Voronezh, 394006, Russian Federation
For citation: Shitikova M. V. Impact of a rigid sphere on an infinite viscoelastic Kirchhoff— Love plate considering volume and shear relaxations. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2023, vol. 10(68), issue 1, pp. 139-154. https://doi.org/10.21638/spbu01.2023.112 (In Russian)
The problem of a low-velocity normal impact of a rigid sphere upon an infinite viscoelastic Kirhhoff—Love plate is considered. The dynamic behaviour of the viscoelastic plate is described by the fractional derivative standard linear solid model. The fractional parameter defining the order of the fractional derivative governs the variation in the viscosity of plate's material within the contact domain during the impact process. The local bearing of the plate material under sphere's indentation, as well as the contact force are defined via the generalized Hertzian contact theory. Using the algebra of Rabotnov's fractional-order operators and taking the volume and shear relaxations into account, the integral equation for the local bearing of the contacting bodies has been obtained. Its approximate solution allows one to find the time dependence of the local indentation and the contact force.
Keywords: low-velocity impact, viscoelastic Kirhhoff—Love plate, fractional derivative standard linear solid model, algebra of Rabotnov's fractional operators.
*The work is supported by Ministry of Science and High Education of the Russian Federation (project no. FZGM-2020-0007).
References
1. Goldsmith W. Impact. The theory and physical behaviour of colliding bodies. London, Arnold (1960). [Rus. ed.: Goldsmith W. Udar. Teoriy i fizicheskie svoistva soudaryemych tel. Mir Publ. (1975)].
2. Rossikhin Yu.A., Shitikova M.V. Transient response of thin bodies subjected to impact: Wave approach. The Shock and Vibration Digest 39 (4), 273-309 (2007). https://doi.org/10.1177 /0583102407080410
3. Zener C. The intrinsic inelasticity of large plates. Physical Review 59 (8), 669-673 (1941). https://doi.org/10.1103/PHYSREV.59.669
4. Phillips J.W., Calvit H.H. Impact of a rigid sphere on a viscoelastic plate. Journal of Applied Mechanics 34 (4), 873-878 (1967). https://doi.org/10.1115/L3607850
5. Shitikova M. V. Wave theory of impact and Professor Yury Rossikhin contribution in the field (A memorial survey). Journal of Materials Engineering and Performance 28 (6), 3161-3173 (2019). https://doi.org/10.1007/s11665-018-3824-6
6. Shitikova M.V. Fractional operator viscoelastic models in dynamic problems of mechanics of solids: A review. Isvestiy Rossiiskoi akademii nauk. Mekhanika tverdogo tela 1, 3-40 (2022). https://doi.org/10.31857/S0572329921060118 (In Russian). [Eng. transl.: Mechanics of Solids 57 (1), 1-33 (2022). https://doi.org/10.3103/S0025654422010022].
7. Shitikova M. V., Krusser A. I. Models of viscoelastic materials: a review on historical development and formulation. Advanced Structured Materials 175, 285-326 (2022). https://doi.org/10.1007/978-3-031-04548-6-14
8. Rossikhin Yu.A., Shitikova M.V. Fractional calculus in structural mechanics. In: Handbook of Fractional Calculus with Applications. Vol. 7: Applications in Engineering, Life and Social Sciences. Part A, 159-192, D. Baleanu, A.M.Lopes (eds). Berlin, De Gruyter (2019). https://doi.org/10.1515/9783110571905-009
9. Rossikhin Yu.A., Shitikova M.V. Classical beams and plates in a fractional derivative medium, Impact response. In: Encyclopedia of Continuum Mechanics, 294-305. H. Altenbach, A. (Ochsner (eds). Germany, Springer-Verlag GmbH (2020). https://doi.org/10.1007/978-3-662-53605-6-86-1
10. Rossikhin Yu. A., Shitikova M. V. Fractional derivative Timoshenko beams and Uflyand-Mindlin plates, Impact response. In: Encyclopedia of Continuum Mechanics, 962-971, H. Altenbach, A. (Ochsner (eds). Germany, Springer-Verlag GmbH (2020). https://doi.org/10.1007/978-3-662-53605-6-87-1
11. Rossikhin Yu.A., Shitikova M.V. Collision of two spherical shells. Fractional operator models. Encyclopedia of Continuum Mechanics 324-332, H. Altenbach, A. (Ochsner (eds). Germany, SpringerVerlag GmbH (2020). https://doi.org/10.1007/978-3-662-53605-6-88-1
12. Rossikhin Yu. A., Shitikova M. V. Two approaches for studying the impact response of viscoelastic engineering systems: An overview. Computers and Mathematics with Applications 66 (5), 755-773
(2013). https://doi.org/10.1016/jxamwa.2013.01.006
13. Meshkov S. I., Pachevskaya G. N., Postnikov V. S. To the detection of volume relaxation during dynamic tests. Proceedings of Erevan State University 2, 10-14 (1968). (In Russian)
14. Meshkov S. I., Pachevskaya G.N., Allowance for bulk relaxation by the method of internal friction. Journal of Applied Mechanics and Technical Physics 2, 80-82 (1967). (In Russian) [Engl. transl.: Journal of Applied Mechanics and Technical Physics 8 (2), 47-48 (1967). https://doi.org/10.1007/BF00918032].
15. Shitikova M.V., Popov 1.1., Rossikhin Yu.A. Theoretical and experimental evidence of the bulk relaxation peak on the loss tangent versus frequency diagrams for concrete. Mechanics of Advanced, Materials and Structures (2022). https://doi.org/10.1080/15376494.2021.2012857
16. Makris N. Three-dimensional constitutive viscoelastic laws with fractional order time derivatives. Journal of Rheology 41, 1007-1020 (1997). https://doi.org/10.1122/L550823
17. Rossikhin Yu. A., Shitikova M. V. Features of fractional operators involving fractional derivatives and their applications to the problems of mechanics of solids. In: Fractional Calculus: History, Theory and Applications, ed. by R. Daou, M. Xavier. Chapter 8, 165-226. New York, Nova Science Publ. (2015).
18. Alotta G., Barrera O., Cocks A., Di Paola M. On the behavior of a three-dimensional fractional viscoelastic constitutive model. Meccanica, 52 (9), 2127-2142 (2017). https://doi.org/10.1007/s11012-016-0550-8
19. Samko S. G., Kilbas A.A., Marichev O.I. Fractional Integrals and Derivatives and Some of Their Applications, New York, Gordon & Breach Sci. Publishers (1993).
20. Rossikhin Yu.A., Shitikova M.V. Centennial jubilee of Academician Rabotnov and contemporary handling of his fractional operator. Fractional Calculus and Applied Analysis 17 (3), 674-683
(2014). https://doi.org/10.2478/s13540-014-0192-2
21. Rabotnov Yu. N. Equilibrium of an elastic medium with after-effect. Fractional Fractional Calculus and Applied Analysis 12 (1), 53—62 (1948). (In Russian) [Engl. transl.: Fractional Calculus and Applied Analysis 17 (3), 684-696 (2014). https://doi.org/10.2478/s13540-014-0193-1].
22. Johnson K.L. Contact Mechanics. Cambridge, Cambridge University Press (1985). [Rus. ed.: Johnson K.L. Mechanika kontactnogo vsaimodeistvia. Moscow, Mir Publ. (1989)].
23. Rabotnov Yu. N. Elements of Hereditary Solid Mechanics. Moscow, Mir Publ. (1980).
Received: May 26, 2022 Revised: July 24, 2022 Accepted: September 8, 2022
Author's information:
Marina V. Shitikova — mvs@vgasu.vrn.ru