CC BY
10
Тарлаковский Д.В., Лам Нгуен Ван Распространение нестационарных антисимметричных кинематических возмущений от сферической полости в среде Коссера // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2020. - № 4. - С. 201-210. DOI: 10Л5593/perm.mech/2020.4Л7
Tarlakovskii D.V., Lam Nguyen Van Propagation of non-stationary antisymmetric kinematic perturbations from a spherical cavity in Cosserat medium. PNRPU Mechanics Bulletin, 2020, no. 4, pp. 201-210. DOI: 10.15593/perm.mech/2020.4.17
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 4,2020 PNRPU MECHANICS BULLETIN
https://ered.pstu.ru/index.php/mechanics/index
DOI: 10Л5593/perm.mech/2020.4Л7 УДК 539.3
РАСПРОСТРАНЕНИЕ НЕСТАЦИОНАРНЫХ АНТИСИММЕТРИЧНЫХ КИНЕМАТИЧЕСКИХ ВОЗМУЩЕНИЙ ОТ СФЕРИЧЕСКОЙ ПОЛОСТИ В СРЕДЕ КОССЕРА
Д.В. Тарлаковский1' 2, Нгуен Ван Лам2
Научно-исследовательский институт механики Московского государственного университета, Москва, Россия
Московский авиационный институт, Москва, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 20 октября 2020 г. Принята: 20 ноября 2020 г. Опубликована: 30 декабря 2020 г.
Ключевые слова: среда Коссера, пространство со сферической полостью, антисимметрия, функции влияния, сферические функции, преобразование Лапласа, малый параметр.
Рассматривается заполненное линейно-упругой средой Коссера пространство со сферической полостью при заданных нестационарных антисимметричных поверхностных возмущениях, под которыми понимается соответствующий аналог классических антиплоских деформаций. Движение среды описывается системой из трех уравнений относительно ненулевых компонент вектора перемещения и потенциалов поля вращения, записанной в сферической системе координат с началом в ее центре полости. Начальные условия полагаются нулевыми.
Для решения задачи используются разложения искомых функций в ряды по полиномам Лежандра и Гегенбауэра, а также преобразование Лапласа по времени. В результате задача сводится к независимым системам обыкновенных дифференциальных уравнений с оператором Лапласа относительно коэффициентов рядов. Формулируется утверждение о структуре общего решения этой системы. Изображения коэффициентов рядов представлены в виде линейных комбинаций граничных условий с коэффициентами - трансформантами поверхностных функций влияния, в явные формулы для которых входят функции Бесселя полуцелого индекса. В силу сложности этих выражений для определения оригиналов в линейном приближении используется метод малого параметра, в качестве которого принимается коэффициент, характеризующий связь полей перемещений и поворотов. Затем с учетом связи функций Бесселя с элементарными функциями изображения записываются в виде линейных комбинаций экспонент с коэффициентами - рациональными функциями параметра преобразования. Дальнейшая процедура обращения преобразования Лапласа проводится с использованием вычетов. Показано, что имеются три волновых фронта, соответствующие модифицированной с учетом свободного вращения волне сдвига и двум волнам вращения.
Представлены примеры расчетов для зернистого композита из алюминиевой дроби в эпоксидной матрице.
©ПНИПУ
© Тарлаковский Дмитрий Валентинович - д.ф.-м.н., проф., зав. лаб., e-mail: tdvhome@mail.ru
: 0000-0002-5694-9253. Нгуен Ван Лам-асп., e-mail: nvlammai2019@gmail.com, : 0000-0001-5695-7176.
Dmitrii V. Tarlakovskii - Doctor of Physical and Mathematical Sciences, Professor, Head of the Laboratory, e-mail: tdvhome@mail.ru, : 0000-0002-5694-9253.
Nguyen Van Lam - PhD Student, e-mail: nvlammai2019@gmail.com, : 0000-0001 -5695-7176.
Эта статья доступна в соответствии с условиями лицензии 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)
PROPAGATION OF NON-STATIONARY ANTISYMMETRIC KINEMATIC PERTURBATIONS FROM A SPHERICAL CAVITY IN COSSERAT MEDIUM
D.V. Tarlakovskii1' 2, Nguyen Van Lam2
1Scientific Research Institute of Mechanics of the Moscow State University, Moscow, Russian Federation Moscow Aviation Institute, Moscow, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 20 October 2020 Accepted: 20 November 2020 Published: 30 December 2020
Keywords:
Cosserat medium, Space with a spherical cavity, Antisymmetry, Influence functions, Spherical functions, Laplace transform, Small parameter
We consider a space filled with a linearly elastic Cosserat medium with a spherical cavity under given nonstationary antisymmetric surface perturbations, which are understood as the corresponding analogue of classical antiplane deformations. The motion of a medium is described by a system of three equations with respect to nonzero components of the displacement vector and potentials of the rotation field, written in a spherical coordinate system with the origin at its center of the cavity. The initial conditions are assumed to be zero.
To solve the problem, we use decomposition of functions to Legendre and Gegenbauer polynomials, as well as the Laplace transform in time. As a result, the problem is reduced to independent systems of ordinary differential equations with the Laplace operator for the coefficients of the series. A statement about the structure of the general solution of this system is formulated. Images of the series coefficients are presented in the form of linear combinations of boundary conditions with coefficients - transformants of surface influence functions, the explicit formulas for which include the Bessel functions of a half-integer index. Due to the complexity of these expressions, to determine the originals in the linear approximation, the method of a small parameter is used, which is taken as a coefficient characterizing the relationship between the displacement and rotation fields. Then, taking into account the connection between the Bessel functions and elementary functions, the images are written in the form of linear combinations of exponentials with coefficients - rational functions of the transformation parameter. The further procedure for inverting the Laplace transform is carried out using residues. It is shown that there are three wave fronts corresponding to a shear wave modified with allowance for free rotation and two rotation waves.
Examples of calculations for a granular composite of aluminum shot in an epoxy matrix are presented.
©PNRPU
Введение
С развитием современной науки и техники требуется точное знание процессов деформирования материалов с усложненной структурой. Этому требованию отвечают модели упругих моментных сред, к которым, в том числе, относится модель Коссера. Общая теория такой несимметричной теории упругости впервые была разработана братьями Коссера (Э. и Ф. Коссера) [1].
В настоящий момент моментная теория упругости привлекает внимание многих исследователей. В [2] рассмотрена динамическая связанная осесимметричная задача микрополярной теории термоупругости для изотропных слоя, полупространства или пространства. А в работах [3-5] получены решения аналогичных нестационарных осесимметричных задач для однородной изотропной среды псевдокоссера. В статьях [6,7] построены решения двумерных нестационарных задач для упругих моментных полупространств и полуплоскости. В работах [8, 9] исследованы осесимметричные задачи для упругих тел со сферическими границами. В [10, 11] исследованы особенности распространения плоских периодических и уединенных волн. B работе [12] рассматривается задача о распространении поверхностных волн в среде Коссера (случай полупространства), в [13] исследовано распро-странене нестационарных поверхностных возмущений для
полуплоскости, заполненной псевдоконтинуумом Коссера. Континуальная модель слоистой среды изучалась Н.В. Зволинским и К.Н. Шхинек в работе [14].
В статье [15] рассматривается начально-краевая задача для тонкой пластинки с позиций общей трехмерной несимметричной теории упругости. В статье [16] исследуется динамическая задача для микрополярных упругих тел при помощи использования метода собственных значений. В [17] рассматривается динамическая связанная осесимметричная задача микрополярной теории упругости для бесконечной в радиальном направлении изотропной среды. В статье [18] получен ряд новых аналитических решений статических и динамических волновых задач линейной упругой среды Коссера. В цикле работ М.А. Кулеша рассматриваются задачи о распространении объемных продольных и поперечных волн [19]. В статье [20] изучаются одномерные динамические уравнения микрополярных упругих тонких балок со свободным вращением, со стесненным вращением и «малой сдвиговой жесткостью». В том числе рассматриваются свободные колебания балок при шарнирном опирании на концах. В статье [21] динамическая задача моментной теории упругости о трещине конечной длины при нормальной нагрузке на берегах методом интегральных преобразований сводится к системе сингулярных интегральных уравнений относительно пе-
ремещений и поворотов, которая решается численно. Нелинейные моментные теории упругости рассматриваются в работах [22-26]. В статье [27] рассматриваются поверхностные волны в упругом слоистом полупространстве с периодическим изменением жесткости по глубине.
Целью работы является постановка задач о распространении нестационарных антисимметричных кинематических возмущений от сферической полости в среде Коссера.
1. Постановка задачи
Рассмотрим однородное изотропное упругое мо-ментное пространство со сферической полостью радиусом R, занятое средой Коссера [28]. Используем сферическую систему координат r, 9, где r > 0,0 <9<п, 0 <$< 2п, с началом в центре полости и оронормиро-ванным базисом er, e9, ed .
Полагаем, что движение является аналогом рассматриваемой в теории упругости антиплоской деформации - антисимметричным относительно полупрямой 9 = 0, что соответствует следующим полям перемещений u и вращения ш:
u = u(r,9,t), ш = ю, (r,9,t)er + Ю9 (r,9,t)e9, ¥ = Y(r,9,t), Ф = Ф(г,9,t).
При этом координаты вектора ш связаны со скалярным Ф и векторным ¥ = ¥(r, 9, t) потенциалами следующим образом:
дФ 1 д (ш . 9) 1 дФ 1 d(r¥) т. = —+-——(sin9), Ю9 =——----^—-. (1)
dr r sin 8 Э8
r д8 r dr
Такое движение среды описывается тремя уравнениями [28]:
ди=4 д (и)-^ Л (ф) дt р
д 2Ф
4а д2 Ф
2а,
(2)
dt 2 = с4 ЛФ- -у Ф -д-Г = clD (Ф)+Т (и - 2ф)
где
D(и) = Ли--U—, Лф = -2Afr2 ^
r sin 8 r dr l dr
A f sin8*
r2 d81 d8
P + 2Y c2 = , + а J P
Здесь г - время; д, р и 3 - упругая постоянная Ламе, плотность и мера инерции среды при вращении (плотность момента инерции); а, в, у, £ - физические параметры континуума Коссера, характеризующие наличие моментных эффектов; с2, с3, с4 и с2а - скорости распространения волн; А - оператор Лапласа для скалярной функции.
Полагаем, что в начальный момент времени г = 0 возмущения отсутствуют:
I 0 ди
и = 0,— Ь=° dt
= 0, Ф = 0,
' it=0
дФ
"д7
=0,
1=0 = 0
"дТ
(3)
= 0,
а на границе полости заданы кинематические возмущения
и
=1 = Ü0 (8, Т), fflr|r=1 ="r 0 (8, т), ®8 L =^80 (8, т). (4)
Решение начально-краевой задачи (1)-(4) разыскивается в классе ограниченных функций.
Выражения для компонент тензоров напряжений и моментных напряжений в силу громоздкости здесь не приводятся, и далее ограничимся нахождением только кинематических параметров среды.
Отметим, что частный случай граничных условий (4) соответствует жестко сцепленному с полостью единичного радиуса абсолютно твердому шару, поворачивающемуся вокруг оси 8 = 0 по закону Q(t) . Действительно, при этом вектор поворота & и правые части последних двух равенств в (4) имеют вид:
& = Q(cos8er -sin8e8),
Or0 = Q cos 8, QJ
, т) = -Q sin 8.
Вектор перемещения поверхности шара и и правая часть первого равенства в (4) определяются так:
и = ОбЩ8ед, и0 = ОБЩ8 .
Далее везде будем использовать безразмерные величины (при одинаковом начертании они обозначены штрихами, которые в последующем изложении опущены):
' г с1г , и , и0 , Ф
г = —, т = —, и = —, и0=—, Ф =—,
R
Ф' =
R
R
R
2 2 2
у2 = О- у2 = О- у2 = £l 10 2 ' 11 2 ' 12 2
^ + 2u , а pL + -2
а = —, U = V, П1а=Г1 +а,
Pci
J
где X - упругая постоянная Ламе; с1 - скорость распространения волн растяжения-сжатия.
В этих величинах кинематические соотношения (1) и граничные условия (4) сохраняют свой вид, а уравнения (2) и начальные условия (3) записываются следующим образом (точками здесь и далее обозначены производные по безразмерному времени т):
Ф = у02 ЛФ-4аиФ, и = п+аD(и)-2oD(Ф), Ф = у-2D(Ф) + 2аи(и-2Ф).
(5)
= Ф| = Ф| = ф| = ф| = 0. (6)
1т=0 т=0 1т=0 т=0 v '
t=0
L
С1 =
2. Представление решения в виде рядов
Потенциалы, перемещения и координаты вектора вращения представляем в виде рядов по полиномам Лежандра Рп (х) и Гегенбауэра СП- (х) [30]:
Ф I = £ (Т" I P (cos е),
«r ) n=0 V «
¥
V «е )
¥ n
= - sin е£
n=1
Подставляя их в (1), получаем
V «n )
cTi (cos е).
дФ ¥
«„ =-г"-n(n +1)—- (n >0),
«е„ =
dr ' r
Фп -1
r r dr
(7)
(n >1).
При этом уравнения (5) и условия (4), (6) переходят в следующие равенства:
фn =Yo2Anфn -4аифn (n > 0),
un =П+аАnun - 2oAn¥n, ¥n = Y-2An¥n + 2au( -2¥n) ( > 1).
ЮД=1 = « on (T)(" > 0) un|r=1 = u0 n (T) «enlr=1 = «H>" (T) (n > 1).
ф| = ф| = 0 ( > 0), ul =
n 1т=0 n |x=0 V )' n |x=0
= 0 (n > 1).
= u = ¥ = ¥
n 1т=0 n 1т=0 n |t=0
(8)
(9)
(10)
Здесь использованы разложения в ряды правых частей условий (4):
U (е, т) = - sin е£ u0" (t)c"-21 (cos e),
n=1
Or 0 (e, T) = £ Юг 0n (t)P" (cos e),
n=0
Oe0 (e, T) = - sin e£ Юе0n (т) c"- (cos e).
n=1
Соотношения (7)-(10) при каждом n вместе с требованием ограниченности образуют независимые начально-краевые задачи.
3. Изображения общего решения
Для решения этих задач применяем к ним преобразование Лапласа по времени (S - параметр; индекс «L» обозначает трансформанту). В результате получаем краевые задачи относительно изображений:
s^L =Y02AnФL - 4a^L ( > 0).
s2uL = n+AuL -2aAn¥„L, s2¥L = 2auuL + (2a" - 4au) ¥L (n > 1).
(11) (12)
«d 1 =®l0" (s) (n >0 ) uL\ 1=uL (s)
lr=1 ir=1
юЦ=1 =«l0" (s)(n > 1),
где
ЭФ!
¥L
«L = ^-n( +1) — ( > 0),
dr r
ф 1 d(r ¥ L)
«L. =ф^ - (n > 1).
(14)
r r dr
Общее решение уравнения (11) имеет вид [31-33]
ФL = £Ci0'(s)Zn (А,0r), ^ =у^s2 + 4au (Re^ > 0),
1=1
Z1n (z) = А=К„+У2 (z) Z2n (z) = "Ц+1/2 (z).
Vz vz
(15)
Здесь С(0) - постоянные интегрирования; (2) и К (2) - модифицированные функции Бесселя первого и второго рода порядка V.
Для построения общего решения системы из второго и третьего уравнений (12) используем достаточно просто доказываемое утверждение.
Пусть матрица системы обыкновенных дифференциальных уравнений
L (y) = Ay, L (y) = £bm-k (x)y(k)(x)
k=0 T
(16)
У = (((x),У2(x),...,у" (x)) , A = (ai,-l (x,aü )
где а^ не зависит от х, имеет простой спектр с собственными значениями \ и собственными векторами ¥к (к = 1,2,...,п) . Тогда ее общее решение имеет вид
у = Tz, Т = (¥1,...,¥п),
Т т
= (21 (х), 22 (х),..., 2п (х)) , = x сыеы (х)
I=1
где {^и (х)} ((= 1,2,...,т) - фудаментальные системы решений уравнений Ь (к )=^к2к, а Сы - произвольные постоянные.
В свете этого утверждения систему уравнений записываем в матричном виде
СА,.
Л,L I ( uL I
u" = B
¥L
B =
( s2
¥ L
V ' n ) V n ) 0
Л
-2au (s2 + 4au)
С =
(n+a -2a1 0 Y-2
v
и приводим ее к канонической форме (16):
( 1 \ и„
¥
V п
= А
( 1\ ип
¥
V п У
А=С-1В=
П+а
(у2252 -4а2и) 2а(2 + 4аи) -2аи< <(2 + 4аи)
Далее строим характеристическое уравнение матрицы А относительно собственных значений X2. Как оказывается, оно сопадает с аналогичным уравнением для осесимметричной задачи [30], и его корни определяются так:
х2,2 =
В0 + аВ1 Б0 -2аЦ +а2Б2
2 (1+ау?)
В0 = (У? + У2) *2, В, = у22 (у?/ + 4и), Б0 = (у? - у2 )2 ^4 а =У2 ((-У2) (у,252 + 4и),
а =у2
у2 ((52 + 4и)2 - 16иу452
Соответствующие Х^ собственные векторы есть решения систем линейных алгебраических уравнений
(А-Х22Е) = 0, У, = (,у2к) .
Их выбираем так, чтобы имела место линейная независимость и при а = 0:
Уп = У2 ( + 4аи) - X2, У21 = 2ау2и,
уи = -2ау;! у2 (2 + 4аи), (17)
У22 =У2 ( - 4а2у2и)-Х2 (1 + ау2).
Следовательно, в соответствии с (16) общее решение системы уравнений (12) имеет вид
(иь ) 2
¥
V п У
£ г I)(5 1 (V).
к ,1=1 V У2к У
Далее, подставляя (15) в (14), с учетом свойств функций Бесселя получаем следующие результаты:
< =Х01с(0)(5)Хп, (Х0г) +
1=1
+ п (п + 1)£ у2кс(:, )(5 )Х А (ХкГ),
к ,1=1
«п =-Х0 £ сп0)(5 К (Х г)+
1=1
+ £ У2к Х С, )(5 )У + 2,п (ХкГ ),
к ,1=1
иЬ = £ Уи^)(51 (ХкГ).
(18)
Здесь использованы обозначения
X п (2) = 2'п, (2) = 1 \п2ы (2) + (-1) 21^ (2)] ,
(19)
У1п (2 )= — 1* (2 ), У +2,п (2 )= 4 (2 )-Хп (2 ).
4. Изображения решения задачи
Поскольку функция Бесселя /п+1/2 (2) неограничена
в окрестности бесконечно удаленной точки [30], то в силу (15) и (19) таковыми же являются и функции 11п (2), Х2п (2), У2п (2). Следовательно, в (18) необходимо положить сШ($) = 0 (у = 0,1,2). Подставляя теперь с учетом этого равенства (18) в (13), получаем алгебраические уравнения относительно постоянных интегрирования:
) Х1(
(5 ) = (5 )
А. =
АпСп = Вп (п > 1), (ХсХ, (Х0) п(п + 1)у21Х1У1п (Х) п(п + 1)у22Х2У1п (Х2)) (20) -ХсУ.п (Х0) у21Х1У3п (Х) у22Х2У3п (Х2)
0 У1111п (Х1 ) У1211п (Х2 )
Г сп0)(5 'Ю10п (Х0 ))
Сп = С(1)(5 ) , Вп = ®1оп (Х0 )
Г ч и0п (Х0) ^
Подставляя их решения в (14), получаем изображения искомых функций:
(г, 5) = 01т (г, 5Ксп (5) + + (Г, 5 К„ (5 ) + ОЬил (Г, 5 У0п ( ),
< (Г, 5) = С^п (Г, 5)<«! (5) + + Ою88п (г, 5Кс (5) + С«8и„ (г, 5) и01„ (5) ,
и1 (Г, 5) = С,1тт (Г, 5Не (5) + + Сю8И (Г, 5) Ю^ (5) + О^ (Г, 5К (5) .
Здесь °Юг™ (Г, 5 Ж--- О1юит (Г, 5) - изображения поверхностных функций влияния, которые определяются так:
ЪпОЮгт =Х0Х1и (Х0Г)Ая11 +п(п + 1)£У2кХкУ1п (ХкГ)А„и+1,
к=1
ДОСг8„ =Х0 Х1и (Х0 Г )Ап +п (п + 1)£ у2 к ХкУ1п (ХкГ )А„ 2, к+1,
к=1
АОЮ„„ = XХ1и (V) Ап31 +п(п +1) £у2кХкУ1п (Х,г) А„3,к+1,
к=1
ДОС8г„ = АУ, (ХГ) Ап11 + £У2кХкУ3„ (Х,г) Ап1,к+1,
к=1
Я„ОЮ88„ = АУ, (Х0Г2! + £ У2кХ^ (ХкГ )А„2,к + 1,
2
к=1
DnGLun = -Vi„ (ог)Аз1 + ЕУ2кV3, (V)An,Ml, (21) к=1
2
DnGLrn = Е У1кАШ,к+1z1n ( V ) к=1 2
DnGurníjn = Еy1kAn2,k+1Z1n (^кГ ), к=1
2
Dn (S )GLn (r, s ) = Е У1кАп3,к+1z1n ( V ), к=1
Dn (s) = |A„| = \ [Xln (^) A.11 (s)- Yln (^) An21 (s)].
В этих равенствах и далее Anjj (s) - алгебраическое
дополнение расположенного в i-й строке и j-м столбце элемента матрицы An.
5. Линейное приближение решения
Получить аналитически оригиналы функций влияния при n > 1 не представляется возможным. Поэтому используем разложения в степенные ряды по малому параметру а, ограничиваясь только линейными слагаемыми. При этом приближенные равенства заменяем точными.
Соответствующие равенства для ^012 и координат (17) собственных векторов имеют вид:
=Y0 I S + ^ j' ^ =у1аS'
^2 =Y2 I s +
'^T ), Y1a = Y1 v1 -aY?.
У11 = ((-Y2)я2 + а(я2 + 4иу2), У21 = 2ау2и, У22 = (Yi! - Y2) я2 - аY2 ((я2 + 4и), уп = -2аYi!Y2я2.
Отсюда с учетом (20) следует, что имеют место следующие асимптотические соотношения:
У12 , У21, An12, An22, Ап31, Ап33 = 0 (а) 0 .
Их использование с учетом связи модифицированных функций Бесселя полуцелого индекса с элементарными функциями [30] позволяет с принятой точностью привести равенства (21) к следующему виду:
GL» (^s )=HL» s)e
-Y1a(r-1)s
Hl«U" (Г, s) =
rn0 (y10
(22)
umu» V 5 У n+1 Г) ( \
r Rn0 (Y1a s)
GL. (r, s ) = a£ H «uL (r, s )e-Yk (r-1)s k=0,2
GLun (r,s) = a£HOk) (r,s)e-Yk(r-1)s,
k=0 2
GL. (r,s) = a£HÜTL (r,s)e
(1k)LÍ„ o^„-Yk(r-1)s
(23)
k=1 2
GLen (r,s) = a£H^ (r,s)e-Yk(r-1)s
где
H^ (r,s) = - "(" + s) (k = 02)s
H (1k)L (r s ) = -■'■'me«» V ) .. n+2 2
Y1rn+2s2P» (Y0s, Y2s)R»0 (Y1s) Y2k« Fletar, s)
Y^2s2p. (Y0s, Y2s) (Y1s Y1Y2 km F,'
(k = 01,2),
7(1k) ' июГИ
H(1k)L (r s ) =__
j r»+1P» (Y0s,Y2s)R»0 (Y1s)
HH (r, s ) = -
Y1Y2km К
(1k)
где
r»+1P» (Y0s,Y2s)R,0 (Y1s)
(k = 1,2),
(k = 1,2), k« =
2 UY1Y2 Y?-Y2 .
KH(r, s ) = [ R»0 (%s )R»3 (Y2 s )--R»3 (Y1s)R»0 (Y2s)]R»1 (Y0rs), F«m„ (r, s ) = P» (Y0 s, Y1s )R"0 (Y2 rs), F^ (r, s ) = »(» + 1)R»0 (Y0 rs )x
x[R»3 (Y1s)R»0 (Y2s)-R»0 (Y1s)R»3 (Y2s)], F^ (r,s) = R»3 (Y1rs)P» (y0s,Y2s), Fe2» (r,s) = -R»3 (Y2rs)P» (Y0s, Ys), FüÜ (r, s) = R» 0 (Y2 s )R» 0 (Y1rs )R» 0 (Y0 s), К12» (r, s) = - R» 0 (v )R»0 (Y2 rs )R» 0 (Y0 s ) К«е» (r,s) = R»0 (Y2s)R»0 (Yr)R»1 (Y0s), Kl (r, s) = -R»0 (Y1s )R»0 (Y2 rs )R»1 (Y0 s).
GL» (r, s) = G<0)l (r, s) + aG<f (r, s),
-Yk (r-1)s
С (r, s)= £ HZ"k)L (r, s)e
(Z = «rr, «re, «er, «ее; l = 0,1).
H»)L (r, s) = -..»+2
С )(r, s)
HZ»k)L (r, s) = -»+2^2
r"+2P» (Y0 s, Y2s )
Z»
2uKK(»lk )(r, s)
r"+2s2P (Y0s, Y2s)
)(r, s ) = Fz(1k )(r, s )P» (y0 s, Y2 s) + + F(ok )(r, s )P«(y0 s, Y2 s).
F«00,)(r,s) = -R»1 (Y0sr)R»3 (Y2s), K0™ (r, s) = n(n +1)R,0 (Y2sr)R,0 (Y0s), Fi00"(r,s) = -,(, + 1)R»1 (Y0sr)R»0 (Y2s), F«0e2»(r,s) = n(n + 1)R»0 (Y2sr)R»1 (Y0s),
F«00»(r,s) = R»0 (Y0sr)R»3 (Y2s),
F«02)(r,s) = -R»0 (Y0s)R»3 (y2sr),
F«00»(r,s) = n(n + 1)R»0 (Y2s)R»0 (Y0sr), F«(¡e»(r,s) = -R»1 (Y0s)R»3 (Y2sr),
(24)
k=0,2
k=1
F^V,s) = R, (YоsrR (s) + Rnl (y„srR (s), FfnV,s) = -n(n + 1)[R,„ (Y2sr)R„5 (y„s)+ R„„ (y„sR (2sr)], F^(r,s) = n(n +1)Rn„ (Y2s)R,2 (Y„sr) + R,i (Y„srR (Y2s)], F(12n(r,s) = -n(n +1)R„„ (Y2rs)R2 (y„s) + R„1 (y„sR (Y2rs)], Fe^r,s) = -R„„ (y„rs)Rn6 (Y2s)-R„3 (Y2s)R»5 (Y„rs), Fiem(r,s) = R„3 (Y2rsR (yоs) + R„o (Y„s)Rn6 (Y2rs), F2(r,s) = -n(n + 1)[R„ (Y2s)R5 (Y„rs) + R„„ (Y„rs)R5 (Y2s)], (r,s) = R„1 (Y„s)Rn6 (Y2rs) + Rn2 (Yos)Rn3 (Y2sr).
Здесь использованы следующие многочлены [33]: n (n + k)!
R„„ (* ) = Z -k, A* =
2k (n-k)!k!'
Яп1 (2) = Я„+.,0 (2)-пЛ„0 (2), Я2 (2) = ЯП4 (2)+ Яи1 (2) , Я^ (2) = Яп1 (2)-К„0 (2),
Я4 (2) = Яи1 (2) + Я„0 (2), К5 (2) = Я (2) + Я„0 (2),Яп6 (2) = Я„4 (2)-Я„0 (2), Рп (х,У) = Яп1 (х)ЯП3 (у) - п(п + 1)Ки0 (х)ЯП0 (у),
Ри(1)(х,у) = Я2 (х)Кп3 (у) +
+Яи1 (х)Яя6 (у)-п(п +1)\Я„5 (х)ЯП0 (у) + Яп0 (х)ЯП5 (у)].
6. Оригиналы решения
Из (22)-(24) следует, что коэффициенты перед экспонентами в изображениях функций влияния являются рациональными функциями параметра 5. Анализ степеней числителей и знаменателей показывает, что среди этих функций только три неправильные дроби Н^,
н!™1 и нЮее«1. Выделение у них целых частей приводит к таким результатам:
HLn (r ,s) =r-1+HL„r (r ,s) H00) (r, s) = r-1 + H(r) (r, s) H(B) (r, s) = r-1 + H(в» (r, s)
(25)
Здесь функции с дополнительным индесом «г» -правильные дроби. Их оригиналы так же, как и оригиналы остальных коэффициентов перед экспонентами, достаточно просто находятся с помощью вычетов. При этом обращение преобразования Лапласа для функций в (25) приводит к следующим равенствам [29, 34]:
ниюп (г т) = г-15(х) + Ниттг (г,х)
НН(г, х) = г-15(х) + Н (1 (Г , х),
нС^ (г,т) = г-15(х) + НЮ81 (г,х),
где 5(т) - дельта-функция Дирака.
Окончательно оригиналы функций влияния в (23)-(25) в соответствии со свойствами преобразования Лапласа имеют следующий вид (Н (т) - единичная функция Хевисайда):
Guraun (r, т) = r-1S [ Х-Y1a (r -^ + Guaunr (r, x),
Guraunr (r, T) = Huraunr [r, T-Y1a(r - 0] H [T-Y1a(r - 0] .
gz» (r,Т) = а z H%) [r,T-Yk (r-1)]H [т-Yk (r-1)]
k=0,2
(z = COru, raBu, afflr, uraB),
Gzn (r,t) = G(n (r,T^aG^ (r,t) (Z = rarr,rarB,raBr,raBB), где
GL (r,t) = r-1S[t-y„ (r-1)] + Go (r,t), GH (r, t) = r-1S [T-Y2 (r -1)] + Gor, t),
Gl- (r, T) = H („) [r, T-Y„ (r -1)] H [r, T-Y„ (r -1) + + H ra*> [r, T-Y2 (r-1) H [r, T-Y2 (r-1)],
g<oL (r, t)=h raee) [r, t-y„ (r -1) H [r, t-y„ (r -1)+
+ Hra„l [r,T- Y2 (r-1)]H [r,T- Y2 (r -1).
Вычисление оригиналов правильных дробей проводится методами компьютерной математики в среде Maple 18.
Для регулярной части одной из функций влияния без громоздких выкладок может быть получен явный вид:
Grarr0r (r, T) = Hrarr0r [r, T-Y„ (r - 1)] H [t-Y„ (r - 1)] ,
HrarrOr (r, T) = - (r -1) r-2 [ Y-1 + 2a и (t + Yor)] e-T/Yo.
7. Примеры
В качестве примера рассмотрим среду в виде композита из алюминиевой дроби в эпоксидной матрице со следующими физическими характеристиками [2, 32]:
ц = 1,89 ГШ; c1 = 2,28 103 ^с ; c2 = 0,929 1 03 ^с;
c3 = 2,48 103 мс; J = 0,429 10-3 к^м ; a = 7,45 МПа .
Скорость c4 находим с помощью формулы, указанной в [35, 36]:
+ а/(рс22) = 2,485 • 103 м/с.
Принимая характерный линейный размер, получаем необходимые безразмерные параметры:
у0 = 0,918; у1 = 2,45; у2 = 0,919; а' = 0,654 10-3; и = 5,10 106.
На рис. 1, 2 представлены результаты расчетов для регулярных составляющих функций влияния ОЮгг1г
и Ою881г . Два фронта волн - т = у, (г -1) и т = У2 (г -1) -
здесь сливаются в один разрыв первого рода в силу того, что для используемой среды имеет место приближенное равенство у0 ~ у2.
k=0
Рис. 1. Распределение регулярных составляющих функций ВЛИЯНИЯ Gmrr1 r и Gmeeir по радиусу в различные моменты времени
Fig. 1. Distribution of the regular components of the influence functions Gmrrir and Gmeeir along the radius at different points in time
Заключение
В линейном приближении по малому параметру аналитически построены поверхностные функции влияния для перемещения и углов поворота в задаче о распространении нестационарных антисимметричных кинематических возмущений от сферической полости в упругом моментном пространстве. Показано, что имеются три волновых фронта, соответствующие модифицированной с учетом свободного вращения волне сдвига и двум волнам вращения. Предложенный алгоритм решения может быть использован как для построения остальных функций влияния, так и для задач с силовыми возмущениями.
GccfSh.lO \--i-1-,-;-г-
"б т - 4 ! »
-5 J..... . 1—2__Г ----I t
Рис. 2. Зависимость регулярной составляющей функции влияния Gmrr1r и Gmee1r от времени при различных значениях радиуса
Fig. 2. Dependence of the regular component of the influence function Gmrrir and Gmeeir on time at different values of the radius
С помощью полученных результатов решения при произвольных возмущениях могут быть представлены в квадратурах (в свертках по времени).
Благодарности
Работа выполнена при финансовой поддержке РФФИ (проект № 18-08-00471).
Acknowledgements
The work was carried out with financial support RFBR (project No. 18-08-00471).
Библиографический список
1. Cosserat E., Cosserat F. Theorie des corps deformables. -Paris: A. Hermann et fils, 1909. - 226 p.
2. Ерофеев В.И. Волновые процессы в твердых телах с микроструктурой. - М.: Изд-во МГУ, 1999. - 328 с.
3. Лай Тхань Туан, Тарлаковский Д.В. Распространение нестационарных кинематических возмущений от сферической полости в псевдоконтинууме Коссера // Механика композиционных материалов и конструкций. - 2011. - Т. 17, № 2. -С. 184-195.
4. Лай Тхань Туан, Тарлаковский Д.В. Распространение нестационарных осе-симметричных возмущений от поверхности шара, заполненного псевдоупругой средой Коссера [Электронный ресурс] // Труды МАИ - 2012. - № 53. - URL: www.trudymai.ru/published.php?ID=29267/.
5. Лай Тхань Туан, Тарлаковский Д.В. Дифракция нестационарных волн на сферической полости в псевдоконтинууме Коссера // Радиоэлектроника. Наносистемы. Информационные технологии. - 2013. - T. 5, № 1. - С. 119-125.
6. Белоносов С.М. Моментная теория упругости. - Владивосток: Дальнаука, 1993. - 148 с.
7. Большаков В.И., Андрианов И.В., Данишевский В.В. Асимптотические методы расчета композитных материалов с учетом внутренней структуры. - Днепропетровск: Пороги, 2008. - 196 с.
8. Аэро Э.Л., Кувшинский Е.В. Континуальная теория асимметрической упругости. Равновесие изотропного тела. -Физика твердого тела. - 1964. - Т. 6, вып. 9. - С. 2689-2699.
9. Аэро Э.Л., Кувшинский Е.В. Основные уравнения теории упругости сред с вращательным взаимодействием частиц. - Физика твердого тела. - 1960. - Т. 2, вып. 7. - С. 1399-1409.
10. Ерофеев В.И., Потапов А.И. Нелинейные продольные волны в упругих средах с моментными напряжениями // Акустический журнал. - 1991. - Т. 37, № 3. - С. 477-483.
11. Ерофеев В. И. Распространение нелинейных сдвиговых волн в твердом теле с микроструктурой // Прикл. механика (Киев). - 1993. - Т. 29, № 4. - С. 18-22.
12. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. О свойствах поверхностных волн в упругой среде Коссера // Математическое моделирование систем и процессов: сб. науч. тр. - Пермь: ПГТУ, 2006. - Вып. 14. - С. 109-113.
13. Суворов Е.М., Тарлаковский Д.В., Федотенков Г.В. Плоская задача об ударе твердого тела по полупространству, моделируемому средой Коссера // ПММ. - 2012. - Т. 76, Вып. 5. - С. 850-859.
14. Зволинский Н.В., Шхинек К.Н. Континуальная модель слоистой среды // Изв. АН СССР. МТТ. - 1984. - № 1. -С. 5-14.
15. Атоян А.А., Саркисян С.О. Задача динамики тонкой пластинки на основе несимметричной теории упругости. -Изв. АН Армении. Мех. - 2004. - Т. 57, № 2. - С. 18-33.
16. Kumar Rajneesh, Singh Ranjit, Chadha Т. К. Eigenvalue approach to second dynamic problem of micropolar elastic solid // Indian J. Pure and Appl. Math. 2003. - Vol. 34, № 5. -P. 743-754.
17. Saxena Hirdeshwar S., Dhaliwal Ranjit S. Application of the eigen-number method to an axisymmetric coupled micropolar thermoelasticity // Bull. Pol. Acad. Sci. Techn. Sci. - 1990. -T. 38, № l. - P. 7-18.
18. Аналитические и численные решения статических и динамических задач несимметричной теории упругости / В.В. Корепанов, М.А. Кулеш, В.П. Матвеенко, И.Н. Шардаков // Физ. мезомех. - 2007. - Т. 10, № 5. - С. 77-90.
19. Анализ волнового решения уравнений эластокине-тики среды коссера в случае плоских объемных волн / М.А. Кулеш, В.П. Матвеенко, М.В. Улитин, И.Н. Шардаков // Прикл. мех. и техн. физ. - 2008. - Т. 49, № 2. - С. 196-203.
20. Birsan Mircea. Thermal stresses in cylindrical Cosserat elastic shells // Eur. J. Mech. A. - 2009. - Vol. 28, № 1. -P. 94-101.
21. Han S.Y., Narasimhan M.N.L., Kennedy T.C. Dynamic propagation of a finite crack in a micropolar elastic solid // Аста месн. - 1990. - Vol. 85, № 3-4. - P. 179-191.
22. Бровко Г.Л., Иванова О.А. Моделирование свойств и движений неоднородного одномерного континуума сложной микроструктуры типа Коссера // Известия РАН. Мех. тверд. тела. - 2008. - № 1. - С. 22-36.
23. Гарагаш И.А., Николаевский В.Н. Механика Коссера для наук о земле // Вычислительная механика сплошных сред. - 2009. - Т. 2, № 4. - С. 44-66.
24. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. Построение аналитического решения волны Лэмба в рамках континуума Коссера // Прикладная механика и техническая физика. - 2007. - Т. 48, № 1. - С. 143-150.
25. Кулеш М.А., Матвеенко В.П., Шардаков И.Н. Построение и анализ аналитического решения для поверхностной волны Рэлея в рамках континуума Коссера // Прикл. мех. и техн. физ. - 2005. - Т. 46, № 4. - С. 116-124.
26. Кунин И.А. Теория упругих сред с микроструктурой. Нелокальная теория упругости. - М.: Наука, 1975. - 416 с.
References
1. Cosserat E., Cosserat F. Theorie des corps deformables. Paris, Librairie Scientifique A. Hermann et Fils, 1909. - 226 p.
2. Erofeev V.I. Volnovye protsessy v tverdykh telakh s mikrostrukturoi [Wave Processes in Solids with Microstructure] (Izd-vo MGU, Moscow, 1999) - 328 p.
3. Lai Tkhan' Tuan, Tarlakovskii D.V. Rasprostranenie nestatsionarnykh kinematicheskikh vozmushchenii ot sfericheskoi polosti v psevdokontinuume Kossera [Propagation of unsteady kinematic perturbations from a spherical cavity in the Cosserat pseudo-continuum] // Mekhanika kompozitsionnykh materialov i konstruktsii. - 2011. - V. 17. - № 2. - pp. 184 - 195.
4. Lai Tkhan' Tuan, Tarlakovskii D. V. Rasprostranenie nestatsionarnykh ose-simmetrichnykh vozmushchenii ot poverkhnosti shara, zapolnennogo psevdouprugoi sredoi Kossera [Elektronnyi zhurn] // Trudy MAI - 2012. № 53. (www.trudymai.ru/published.php?ID=29267/).
5. Lai Tkhan' Tuan, Tarlakovskii D. V. Difraktsiia nestatsionarnykh voln na sfericheskoi polosti v psevdokontinuume Kossera [Diffraction of unsteady waves by a spherical cavity in the Cosserat pseudo-continuum] // Radioelektronika. Nanosistemy. Informatsionnye tekhnologii. - 2013. V. 5, № 1. - pp. 119-125.
6. Belonosov S.M. Momentnaia teoriia uprugosti [Moment theory of elasticity]. - Vladivostok: Izd-vo Dal'nauka, 1993. - 148 p.
7. Bol'shakov V.I., Andrianov I.V., Danishevskii V.V. Asimptoticheskie metody rascheta kompozitnykh materialov s
27. Muhlhaus H.-B., Triantafyllidis Th. Surface waves in a layered half-space with beding stiffness // Ground Motion and Eng. Seismol. Amsterdam e. a. - 1987. - P. 277-290.
28. Новацкий В. Теория упругости. - М.: Мир, 1975 -872 с.
29. Волны в сплошных средах / А.Г. Горшков, А.Л. Мед-ведский, Л.Н. Рабинский, Д.В. Тарлаковский. - М.: ФИЗМАТЛИТ, 2004. - 472 с.
30. Абрамовиц М., Стиган И. Справочник по специальным функциям. - М.: Наука, 1979. - 832 c.
31. Lam V. Nguyen, Tarlakovskii D.V. Propagation of Non-stationary Axisymmetric Perturbations from a Spherical Cavity in Cosserat Medium // Advanced Structured Materials. - Vol. 122. Nonlinear Wave Dynamics. - Springer Nature Switzerland AG, 2020. - P. 273-292.
32. Gauthier R.D., Jahsman W.E. A quest for micropolar elastic constants. Part II // Arch. Mech. - 1981. - Vol. 33, № 5. -P. 717-737.
33. Горшков А.Г., Тарлаковский Д.В. Нестационарная аэрогидроупругость тел сферической формы. - М.: Наука. Гл. ред. физ.-мат. лит., 1990. - 264 с. Перевод: Gorshkov A.G., Tarlakovskiy D.V. Transient Aerohydroelasticity of Spherical Bodies. - Berlin - Heidelberg - New York: Springer-Verlag, 2001. - 289 p.
34. Ван Дер Поль, Бреммер Х. Операционное исчислние на основе двустороннего преобразования Лапласа. - М.: ИЛ, 1952. - 506 с.
35. Деев В.М. Системный анализ уравнений пространственной задачи несимметричной теории упругости в перемещениях // Математическое моделирование в естественных науках: тез. докл. 10-й Всерос. конф. молодых ученых. -Пермь, 2001. - С.14.
36. Герасимов С.И., Ерофеев В.И., Солдатов И.Н. Волновые процессы в сплошных средах. - Саров: Изд-во РФЯЦ-ВНИИЭФ. - 2012. - 260 с.
uchetom vnutrennei struktury [Asymptotic methods for calculating the composite materials taking into account the internal structure]. Dnepropetrovsk, Porogi Publ., 2008. - 196 p.
8. Aero E.L., Kuvshinskii E.V. Kontinual'naia teoriia asimmetricheskoi uprugosti. Ravnovesie izotropnogo tela [Continual theory of asymmetric elasticity]. - Fizika tverdogo tela. -1964. - V.6. No 9. - pp. 2689 - 2699.
9. Aero E.L., Kuvshinskii E.V. Osnovnye uravneniia teorii uprugosti sred s vrashchatel'nym vzaimodeistviem chastits [Basic equations of the theory of elasticity of media with rotational interaction of particles]. - Fizika tverdogo tela. - 1960. - V. 2. Vyp. 7. - pp. 1399 - 1409.
10. Erofeev V.I., Potapov A.I. Nelineinye prodol'nye volny v uprugikh sredakh s momentnymi napriazheniiami [Nonlinear longitudinal waves in elastic media with moment stresses] // Akusticheskii zhurnal. - 1991. - V. 37. № 3. - pp. 477-483.
11. Erofeev V.I. Rasprostranenie nelineinykh sdvigovykh voln v tverdom tele s mikrostrukturoi [Propagation of nonlinear shear waves in a solid with a microstructure] // Prikl. mekhanika (Kiev). - 1993. - V. 29. № 4. - pp.18 - 22.
12. Kulesh M.A., Matveenko V.P., Shardakov I.N. O svoistvakh poverkhnostnykh voln v uprugoi srede Kossera [Properties of surface waves in an elastic Cosserat medium] // Matematicheskoe modelirovanie sistem i protsessov: Sbornik nauchnykh trudov. - Perm': PGTU, 2006. - Vyp. 14. - pp. 109-113.
13. Suvorov E.M., Tarlakovskii D.V., Fedotenkov G.V. Ploskaia zadacha ob udare tverdogo tela po poluprostranstvu, modeliruemomu sredoi Kossera [The plane problem of the impact of a rigid body on a half-space modeled by the Cosserat medium] // PMM. - 2012. - V. 76. - Vyp. 5. - pp. 850-859.
14. Zvolinskij N.V., Shhinek K.N. Kontinual'naja model' sloistoj sredy [Continual model of a layered medium] // Izv. AN SSSR. Mekhanika tverdogo tela. 1984, no 1, pp. 5-14
15. Atoian A.A., Sarkisian S.O. Zadacha dinamiki tonkoi plastinki na osnove nesimmetrichnoi teorii uprugosti [The problem of the dynamics of a thin plate based on the asymmetric theory of elasticity]. - Izv. AN Armenii. Mekh. - 2004. - V. 57. № 2. - pp. 18 - 33.
16. Kumar Rajneesh, Singh Ranjit, Chadha T. K. Eigenvalue approach to second dynamic problem of micropolar elastic solid // Indian J. Pure and Appl. Math. 2003. - V. 34. № 5. - PP. 743-754.
17. Saxena Hirdeshwar S., Dhaliwal Ranjit S. Application of the eigen-number method to an axisymmetric coupled micropolar thermoelasticity // Bull. Pol. Acad. Sci. Techn. Sci. - 1990. -V. 38. № l. - pp. 7 - 18.
18. Korepanov V.V., Kulesh M.A., Matveenko V.P., Shardakov I.N. Analiticheskie i chislennye resheniia staticheskikh i dinamicheskikh zadach nesimmetrichnoi teorii uprugosti [Analytical and numerical solutions of static and dynamic problems of asymmetric theory of elasticity]// Fiz. mezomekh. - 2007. - V. 10. - № 5. -pp. 77-90.
19. Kulesh M.A., Matveenko V.P., Ulitin M.V., Shardakov I.N. Analiz volnovogo resheniia uravnenii elastokinetiki sredy kossera v sluchae ploskikh ob"emnykh voln [Analysis of the wave solution of the equations of elastokinetics of the cosserat medium in the case of plane bulk waves] // Prikl. mekh. i tekhn. fiz. - 2008. -V. 49. - № 2. - pp. 196-203.
20. Birsan Mircea. Thermal stresses in cylindrical Cosserat elastic shells // Eur. J. Mech. A. - 2009. - V. 28. № 1. - P. 94-101.
21. Han S. Y., Narasimhan M. N. L., Kennedy T. C. Dynamic propagation of a finite crack in a micropolar elastic solid // AcTa MecH. - 1990. - V. 85. №3 - 4. - P. 179-191.
22. Brovko G. L., Ivanova O. A. Modelirovanie svoistv i dvizhenii neodnorodnogo odnomernogo kontinuuma slozhnoi mikrostruktury tipa Kossera [Modeling the properties and motions of an inhomogeneous one-dimensional continuum of a complex Cosserat-type microstructure] // Izvestiia RAN. Mekh. tverd. tela. - 2008. - № 1. - pp. 22-36.
23. Garagash I.A., Nikolaevskii V.N. Mekhanika Kossera dlia nauk o zemle [Cosserat Mechanics for Geosciences]// Vychislitel'naia mekhanika sploshnykh sred. - 2009. - V.2. -№ 4. - pp. 44-66.
24. Kulesh M. A., Matveenko V. P., Shardakov I. N. Postroenie analiticheskogo resheniia volny Lemba v ramkakh kontinuuma Kossera [Construction of the Analytical Solution of
the Lamb Wave in the Framework of the Cosserat Continuum] // Prikl. mekh. i tekhn. fiz. - 2007. - V. 48. - № 1. - pp. 143-150.
25. Kulesh M.A., Matveenko V.P., Shardakov I.N. Postroenie i analiz analiticheskogo resheniia dlia poverkhnostnoi volny Releia v ramkakh kontinuuma Kossera [Construction and Analysis of an Analytical Solution for the Surface Rayleigh Wave within the Framework of the Cosserat Continuum] // Prikl. mekh. i tekhn. fiz. - 2005. - V. 46. - № 4. - pp. 116-124.
26. Kunin I.A. Teoriia uprugikh sred s mikrostrukturoi. Nelokal'naia teoriia uprugosti [Theory of Elastic Media with Microstructure. Nonlocal Theory of Elasticity]. - M.: Nauka, 1975. - 416 p.
27. Muhlhaus H.-B., Triantafyllidis Th. Surface waves in a layered half-space with beding stiffness // Ground Motion and Eng. Seismol. Amsterdam e. a. - 1987. - pp. 277 - 290.
28. Novatsky V. Teoriya uprugosti [Theory of Elasticity]. Moscow, Mir, 1975. 872 p.
29. Gorshkov A.G., Medvedskii A.L., Rabinskii L.N., Tarlakovskii D.V. Volny v sploshnykh sredakh [Waves in Continuous Media]. Moscow, FIZMATLIT publ., 2004. 472 p.
30. Abramowitz M, Stegun I. (Eds.) Spravochnik pospezialnym funkziyam [Handbook of mathematicalfunctions]. Ttransl. from engl. Moscow, Nauka Publ.,1979, 832 p.
31. Lam V. Nguyen and Tarlakovskii D.V. Propagation of Non-stationary Axisymmetric Perturbations from a Spherical Cavity in Cosserat Medium // Advanced Structured Materials, V. 122. Nonlinear Wave Dynamics. - Springer Nature Switzerland AG, 2020. - pp. 273 - 292.
32. Gauthier R.D., Jahsman W.E. A quest for micropolar elastic constants. Part II // Arch. Mech. - 1981. - V. 33. № 5. -pp. 717-737.
33. Gorshkov A.G., Tarlakovskii D.V. Nestatsionarnaia aerogidrouprugost' tel sfericheskoi formy [Transient aerohydroelasticity of spherical - shaped bodies]. - M.: Nauka. Gl. red. fiz.-mat. lit., 1990. - 264 p.
34. Van Der Pol'., Bremmer Kh. Operatsionnoe ischislnie na osnove dvustoronnego preobrazovaniia Laplasa [Transient aerohydroelasticity of spherical - shaped bodies]. - M.: IL. 1952. - 506 p.
35. Deev V.M. Sistemnyi analiz uravnenii prostranstvennoi zadachi nesimmetrichnoi teorii uprugosti v peremeshcheniiakh [System analysis of the equations of the spatial problem of the asymmetric theory of elasticity in displacements] // Matematicheskoe modelirovanie v estestvennykh naukakh. Tez.dokl. 10-i Vseros. konf. molodykh uchenykh. - Perm', 2001. - 14p.
36. Gerasimov S.I., Erofeev V.L., Soldatov I.N. Volnovye protsessy v sploshnykh sredakh [Wave Processes in Continuous Media]. Sarov. RFNC - VNIIEF Publ. 2012. 260 p.
